Winter snow accumulation obviously has an effect on the following catchment runoff. The question is, however, how long this effect lasts and how important it is compared to rainfall inputs. Here we investigate the relative importance of snow accumulation on one critical aspect of runoff, namely the summer low flow. This is especially relevant as the expected increase of air temperature might result in decreased snow storage. A decrease of snow will affect soil and groundwater storages during spring and might cause low streamflow values in the subsequent warm season. To understand these potential climate change impacts, a better evaluation of the effects of inter-annual variations in snow accumulation on summer low flow under current conditions is central. The objective in this study was (1) to quantify how long snowmelt affects runoff after melt-out and (2) to estimate the sensitivity of catchments with different elevation ranges to changes in snowpack. To find suitable predictors of summer low flow we used long time series from 14 Alpine and pre-Alpine catchments in Switzerland and computed different variables quantifying winter and spring snow conditions. In general, the results indicated that maximum winter snow water equivalent (SWE) influenced summer low flow, but could expectedly only partly explain the observed inter-annual variations. On average, a decrease of maximum SWE by 10 % caused a decrease of minimum discharge in July by 6–9 % in catchments higher than 2000 m a.s.l. This effect was smaller in middle- and lower-elevation catchments with a decrease of minimum discharge by 2–5 % per 10 % decrease of maximum SWE. For higher- and middle-elevation catchments and years with below-average SWE maximum, the minimum discharge in July decreased to 70–90 % of its normal level. Additionally, a reduction in SWE resulted in earlier low-flow occurrence in some cases. One other important factor was the precipitation between maximum SWE and summer low flow. When only dry preceding conditions in this period were considered, the importance of maximum SWE as a predictor of low flows increased. We assessed the sensitivity of individual catchments to the change of maximum SWE using the non-parametric Theil–Sen approach as well as an elasticity index. Both sensitivity indicators increased with increasing mean catchment elevation, indicating a higher sensitivity of summer low flow to snow accumulation in Alpine catchments compared to lower-elevation pre-Alpine catchments.
The shift from snowfall to rain is one of the most important effects of
predicted climate change on the hydrological cycle (Laghari et al., 2012;
Berghuijs et al., 2014; Zhang et al., 2015). This shift results in a
decrease of the fraction of solid precipitation (snow/total precipitation,
known as
For the western United States the decrease of
Berghuijs et al. (2014) showed that a higher fraction of precipitation fallen as snow is associated with higher long-term mean streamflow in comparison to catchments with lower snowfall fraction. Higher air temperatures during spring affect the onset of snowmelt in streamflow shifting it towards earlier spring (Barnett et al., 2005; Godsey et al., 2014; Langhammer et al., 2015; Ledvinka, 2015). These changes lead to a higher fraction of annual flow occurring earlier in the water year as evident from many studies across the western United States (Cayan et al., 2001; Stewart et al., 2005; Day, 2009). However, snowmelt and consequent spring streamflow are affected by a wide range of factors, such as topography, vegetation and connected radiation as well as shading effects which might overlay the effect of increasing air temperature (Kliment et al., 2011; Kutlakova and Jenicek, 2012; Pomeroy et al., 2012; Kucerova and Jenicek, 2014). Earlier onset of snowmelt could, for instance, be slowed down by less shortwave radiation due to lower sun inclination in early spring (Lundquist and Flint, 2006).
Speich et al. (2015) demonstrated the sensitivity of catchments in the Swiss Alps to a reduction of snow contribution to total runoff by applying bivariate-mapping techniques. The combination of total runoff and snowmelt appeared to be more sensitive to predicted future changes of air temperature and precipitation than the combination of precipitation and potential evapotranspiration (Speich et al., 2015). Additionally, the elevation band between 1000 and 2500 m a.s.l. was found to be relatively more sensitive to future temperature and precipitation scenarios than lower-elevation bands. Further, Zappa and Kan (2007) demonstrated that the presence of above-average snow resources contributed to mitigating the effects of the 2003 summer drought in some high-elevation areas within the Swiss Alps.
Snow conditions in winter can affect low flows during the subsequent summer, especially in areas with large differences in winter and summer precipitation. The total amount of snow precipitation in winter affects groundwater recharge and hence also runoff during dry summer periods (Earman et al., 2006; Beaulieu et al., 2012; Van Loon et al., 2015). While meteorological drivers and overall catchment storage both affect the drought duration during summer, seasonal storage in snow and glaciers affect the drought deficit (Van Loon and Laaha, 2015). However, snow cannot solely explain the sensitivity to drought, although higher-elevation catchments in the Swiss Alps were found to be less sensitive to drought origin (Staudinger et al., 2015). Additionally, some modelling experiments suggested larger groundwater storages in higher-elevation Swiss catchments which may additionally explain the lower sensitivity of higher-elevation catchments to low flows (Staudinger and Seibert, 2014).
Based on historical records from selected Sierra Nevada catchments in the western United States, every 10 % decrease in snow water equivalent maximum in spring leads to a decrease of 9–22 % in minimum runoff during summer months and the runoff minimum occurs about 3–7 days earlier (Godsey et al., 2014). Higher-elevation catchments showed a longer memory to the previous season's climate variability than lower-elevation catchments (Cayan et al., 1993), and some catchments in the Sierra Nevada mountains were affected by the snowpack of the preceding year during the subsequent summer runoff (Godsey et al., 2014).
The above-mentioned studies show that the influence of snow amount on early spring discharge is widely studied and known. However, we still lack a quantitative assessment of the sensitivity of summer low flows on snow conditions from the preceding winter. In this study we want (1) to quantify how long snowmelt affects runoff after melt-out and (2) to estimate the sensitivity of the catchments to changes in snowpack. We benefit from a recently generated snow water equivalent (SWE) data set which allowed for an in-depth analysis of snowpack changes and detection of melt-out dates. Our study adds to earlier studies, by focusing on the combined effect of snow and liquid precipitation during the warm period and its varying importance for individual catchments. Exploring this combined effect is particularly important in humid regions where annual precipitation is approximately equally distributed over the year, while most studies were performed in climates with more seasonal precipitation and/or smaller precipitation amounts overall (such as in the western United States). Furthermore, we set the sensitivities of low flows to varying snow conditions in the context of simple catchment properties which offers a way to indicate regions that might become more vulnerable to droughts in the future.
Study catchments and selected characteristics (
Location of the study catchments in Switzerland.
We selected 14 Alpine and pre-Alpine catchments in Switzerland with a
catchment area ranging from 0.93 to 1577 km
Predictor and response variables used in analyses.
Daily gridded precipitation and air temperature data (2 by 2 km
Daily values of actual evapotranspiration (ET) were calculated using a radiation-based approach included in the PREVAH model (Gurtz et al., 1999; Viviroli et al., 2009). The PREVAH model uses observed global radiation at a daily temporal resolution as input to the Penman–Monteith equation to compute potential evapotranspiration (PET). Actual evapotranspiration is then computed by reducing PET as a function of soil water deficit. This data set was previously used and evaluated e.g. in Speich et al. (2015). The ET data were available for the period 1980–2009.
Daily discharge data were obtained from the Swiss Federal Office for the Environment (BAFU). Data from 1971 to 2012 were used in all analyses except a few shorter time series, as specified in Table 1.
We selected different predictors related to winter and spring meteorological conditions and water storage conditions in the catchments (Table 2).
These predictors were tested to explain the variability of three variables describing low-flow conditions: (i) minimum 7-day moving average of daily discharge was calculated based on BAFU data; different sizes of the moving window (3, 7 and 15 days) were tested without significant influence on the results; (ii) the day of year (DOY) of 7-day minimum of discharge was calculated from June to September to exclude low flows before snowmelt or after the onset of new winter snow accumulation; (iii) number of days below a specified discharge threshold (25 % quantile of discharge from May to October).
We used nine variables as predictors of future summer low flows (Table 2). The advantage of this choice of predictors is that only SWE, precipitation, air temperature and runoff data are needed for their calculation. These data are available for many regions which allows to test our methods also elsewhere with possible transfer to ungauged catchments.
Late winter conditions were represented by the maximum SWE before snowmelt onset calculated using the SWE data from February to May. We used both the maximum SWE calculated as a catchment mean and the maximum SWE calculated from the highest 50 % of the catchment area, assuming that snowpack at higher elevations melts later and could be more important for summer discharges. The melt-out date was calculated from SWE data for each catchment and year. The melt-out date was defined as the first occurrence of snow free conditions (snow cover fraction less than 10 %) after the day of maximum SWE. The sum of positive SWE changes (sum of new snow) and the sum of positive air temperatures were used as well. Both variables were calculated as a sum from 1 November to 30 April.
While the variables related to snow describe the state of the individual
catchment before snowmelt, total winter precipitation calculated from
1 November to 30 April describes the available water amount from winter
precipitation. Additionally, we calculated the fraction of snowfall to total
winter precipitation (
The DOY with maximum SWE was used to show the dependence of low flows on this variable. In this way, we could investigate whether low flows occur later in the year and whether low flows are higher with later occurrence of maximum SWE.
A current precipitation index
The sums of actual evapotranspiration and precipitation from 1 June to 30 September and the maximum SWE were used to assess the inter-annual variability of these variables.
All parameters were calculated assuming a complete data series and additionally considering only years with below-average spring and summer precipitation. By doing this the effect of spring and summer liquid precipitation on low flows could be separated and thus the effect of snow could be highlighted.
To assess the relations between predictors and response variables we used the Spearman rank correlation coefficient and the bivariate linear regression. Most of the predictors and response variables were expressed as a percentage difference from the mean value, which enabled a comparison between individual catchments. The linear regression was computed from log-transformed variables. Prediction intervals of linear regression were used, which allowed the future observation of the response variable to be estimated. The R software was used for all calculations in this study (R Core Team, 2015).
The slope of regression calculated using the nonparametric Theil–Sen method was used to evaluate our statistical models. The Theil–Sen slope is a median of slopes calculated for each pair of observations (Birsan et al., 2005; Pellicciotti et al., 2010). The higher the value, the steeper the slope of regression and thus the more sensitive is the dependent variable (e.g. minimum discharge) to the change of the independent variable (e.g. maximum SWE). The Theil–Sen linear regression model is suitable for non-normally distributed data with outliers.
Similar to the slope of regression, the elasticity index (Eq. 2) was used to
describe how sensitive the minimum discharge is to the change of SWE. The
climate elasticity is often used to describe sensitivity of streamflow to
the change of climate variables (Andréassian et al., 2015). A similar
concept was used in this study to describe what percentage change of minimum
discharge is caused by a defined percentage change of maximum SWE by the
elasticity,
The relative influence of snow and liquid precipitation during the warm
season on low flows was analysed calculating scores for both maximum SWE and
All analyses were done separately for each catchment and almost all for the
period May to September to highlight the changing importance of snow
contribution to low flows in different catchments and time. Analyses of the
combined effect of snow and liquid precipitation were made only for the
period from June to September, because liquid precipitation (expressed as
Spearman rank correlation coefficients between predictors and response variables were calculated separately for three elevation classes (highest-elevation catchments: above 2000 m a.s.l.; middle-elevation catchments: 1300–2000 m a.s.l.; low-elevation catchments: 850–1300 m a.s.l.) (Fig. 2). Using these three elevation classes showed changing correlations for catchments in different elevation and thus different influence of snow storage on runoff. Spearman rank correlation coefficients were displayed as heat maps together with dendrograms showing clusters (based on hierarchical cluster analysis) of similar predictors and response variables (Fig. 2).
Heat maps showing Spearman rank correlation coefficients for all predictors (rows) and response variables (columns) separately for three elevation groups. Left panel: catchments with mean elevation higher than 2000 m a.s.l.; middle panel: catchments with mean elevation between 1300 and 2000 m a.s.l.; right panel: catchments with mean elevation between 850 and 1300 m a.s.l. Hierarchical cluster analysis and Euclidean distance were used to show similarity of individual predictors and response variables. Grey colour used for NA values.
The maximum SWE (both averaged per catchment and in upper 50 % of the
catchment area) was in most cases the best predictor for higher-elevation
catchments during summer (July and later). Additionally, maximum SWE and the
sum of new SWE were better predictors than winter precipitation in snow-dominated catchments to predict the number of days with low discharge
(number of days with runoff below
Minimum discharge and
The role of spring and summer liquid precipitation (expressed as
Dependence of 7-day minimum discharge on maximum SWE for individual months. Top panels: Ova da Cluozza River representing a high-elevation catchment with a mean catchment elevation of 2361 m a.s.l., correlations from May to September are statistically significant (0.05 level). Middle panels: Simme River, representing a middle-elevation catchment with a mean catchment elevation of 1632 m a.s.l., correlations from May to June are significant. Bottom panels: Sitter River as a representative of a low-elevation catchment with a mean catchment elevation of 1247 m a.s.l., only the correlation in May is significant. Solid lines represent the low flow occurring with a 50 % probability, dotted lines represent the 95 % prediction interval. The boxes represent the 25 and 75 % quantiles, with the thick line showing the median, and the whiskers represent minimum and maximum values.
Despite the significance of the correlations, their values are not high, which indicates that low flows are influenced by more than a single variable (maximum explained inter-annual variability in the group was 60 %). Additionally, some of the predictors are not mutually independent (see dendrograms in Fig. 2). Since our focus was primarily on middle- and high-elevation catchments as well as on summer months, maximum SWE seems to be the best predictor, although differences are not large.
We used maximum SWE as variable to predict 7-day minimum discharge (Fig. 3). Snow influence decreased with monthly progression as shown for selected catchments representing high-, middle- and low-elevation ranges. For the three elevation ranges maximum SWE differed not only in the overall amount but also in its inter-annual variability (Fig. 3).
The relationships between the 7-day minimum discharge and maximum SWE
(Fig. 3) are characterized by a large variability. This indicates that only a
certain portion of low-flow variability can be explained using maximum SWE.
Coefficients of determination (
The relationship between predictor and response variable were in more detail described by the Theil–Sen slope of the regression and by the elasticity index (Figs. 4–6). The Theil–Sen slope and elasticity describe the sensitivity of low flows to both decrease and increase of maximum SWE.
Elasticity index for all catchments classified according to
elevation describing the sensitivity of 7-day minimum discharge on maximum
SWE for individual months. Elevation classes on
The elasticity index for high- (mean elevation higher than 2000 m a.s.l.), middle- (mean elevation between 1300 and 2000 m a.s.l.) and low- (mean elevation between 850 and 1300 m a.s.l.) elevation catchments (Fig. 4) decreased progressing from June to September. The elasticity index for high-elevation catchments was for each month higher than the elasticity index for middle- and low-elevation catchments and the elasticity index of the middle-elevation catchments was higher than for the low-elevation catchments. While the spread of the elasticity indices per elevation class was about equal for high- and middle-elevation catchments the spread for the low-elevation increased progressing from month to month. This means that the general sensitivity to SWE is lower for lower- and middle-elevation catchments than for high-elevation catchments and decreases for each class progressing from June to September.
Theil–Sen slopes for each catchment and for every week from the beginning of May to the end of September, allowed an analysis of the sensitivity in terms of the memory effect of each catchment (Fig. 5). These weekly slopes describe how long water from snowmelt contributes to runoff formation and thus how long snowmelt affects low flows. With this approach, a significant effect of snow on low flows became visible during the whole summer and until September for catchments higher than 2000 m a.s.l. (Fig. 5). Snow affected low flows until July in catchments with mean catchment elevation in the range of 1500 to 2000 m a.s.l. However, snow did not affect summer low flow (July to September) in catchments lower than 1500 m a.s.l. Here, snow affected low flows during May and June only, which is probably caused mostly by lower SWE (maximum less than 250 mm).
There was a clearly longer snowmelt contribution to minimum discharges in higher-elevation catchments even when different melt-out dates (black points in the Fig. 5) were considered. This could be related to more available water released from snow in higher-elevation catchments despite their steeper slopes and shallow soils. The negative correlations in some of the lower-elevation catchments (usually not statistically significant) indicate a mixed effect of snow and liquid precipitation in the warm season.
Additionally, Theil–Sen slopes for each catchment and for every week from the beginning of May to the end of September were calculated and set in context to the melt-out days only for situations with dry preceding conditions, i.e. when liquid precipitation prior to minimum discharge in a specific week was below average (results not shown in the paper). Considering only these situations, the sensitivity of minimum discharges to maximum SWE increased. This was due to the reduced influence of liquid precipitation in the warm season. As a consequence, snow became more important and thus any decrease/increase of snow storage in individual year resulted in a more sensitive response of minimum discharge.
The sensitivity as described by Theil–Sen slopes of individual catchments
strongly depends on catchment properties, such as mean catchment elevation,
maximum SWE and
Spearman rank correlation coefficients for the relation between catchment properties and Theil–Sen slopes (TS), which were computed for assessing the low-flow sensitivity to peak SWE. Statistically significant correlations (at the 0.05 level) are shown in bold.
Dependence of 7-day minimum discharge on maximum SWE for all studied catchments (sorted by elevation from highest to lowest) for individual weeks from the beginning of May (week 19) to the end of September (week 39). April (weeks 15–18) was not included in calculation. Colour key provides Theil–Sen slope values. Red indicates positive effect of SWE on minimum discharge (positive slopes), blue indicates negative effect of SWE on minimum discharge (negative slopes). Black points indicate average week of melt-out, whiskers represent 10 and 90 % quantiles.
The elasticity calculated for the 50 % probability of prediction enables us to describe the impact of future changes of snowpack (Fig. 6). The elasticity index in this study is usually lower than 1, which means that a particular percentage change in maximum SWE causes a lower percentage change of minimum discharge. For catchments higher than 2000 m a.s.l., every decrease of the maximum SWE by 10 % will cause a decrease of minimum discharge in July by 6–9 % (Fig. 6, top right panel). This means that the decrease of minimum discharge is almost proportional to the decrease of SWE in some cases (Ova Da Cluozza and Ova dal Fuorn). For catchments with a mean elevation between 1500 and 2000 m a.s.l., the decrease of minimum discharge ranges from 2 % (Grande Eau) to 5 % (Simme). The lowest catchments are characterized with even lower values indicating that any decrease of maximum SWE will not significantly affect low flows at least from July to September. However, there is some small effect during June (Fig. 6, top left panel). Generally, the sensitivity of low flows to the change of SWE increases with elevation and decreases from June to September. However, the elasticity is not linear and the decrease of low flows accelerates with decreasing SWE.
Elasticity index showing the sensitivity of minimum discharge to
changes in SWE. The index was calculated from the 50 % probability of
prediction. Line colours indicate the catchment group according to mean
elevation (dark brown:
Day of year (DOY) with 7-day minimum discharge against long-term
mean annual maximum SWE. Colour circles represent catchment mean and whiskers
represent 10 and 90 % quantiles. The DOY “1” represents the first
day of calendar year (1.1) and the DOY “365” represents 31.12. The colour
of the circle indicates the catchment group according to mean elevation
(dark brown:
Top panels: 7-day minimum discharge in July against maximum SWE for
years grouped according to the current precipitation index
Score plots indicating the combined effect of snow and liquid
precipitation on low flows in the different months (four plots for June to
September). Points below the one-to-one line indicate catchments with a
stronger effect of maximum SWE on low flows compared to spring and summer
precipitation (expressed as
Maximum SWE for each catchment impacted the DOY with 7-day minimum discharge (Fig. 7). Our hypothesis was that minimum summer discharge would occur later in the year for higher maximum SWE. However, later low-flow occurrence may be additionally influenced by a later melt-out. Low flows occurred in September and October for higher-elevation catchments with a higher SWE maximum (Fig. 7, brown points). In contrast, July and August are typical months for low-flow occurrences for lower-elevation catchments with lower SWE maximum (Fig. 7, green points). On average, every decrease in maximum SWE by 100 mm resulted in discharge minima occurring about 12 days earlier. However, inter-annual variability markedly increases in lower-elevation catchments indicating an increasing role of summer precipitation. Thus, the relation between maximum SWE and the DOY of 7-day minimum discharge cannot be used predictively. Additionally, the Spearman rank correlations between maximum SWE and the DOY with 7-day minimum discharge were significant only in some higher- and middle-elevation catchments (above 1500 m a.s.l.) and their values were rather low (detailed results not shown).
The relation between snow and minimum discharge during the summer period is not often clear and may be overlaid by several other factors, mostly precipitation after melt-out. To demonstrate the combined effect of snow and precipitation on summer low flows, three snow-dominated catchments at high and middle elevations (Ova da Cluozza, Vorderrhein and Lümpenenbach) were selected as typical representatives and then further analysed.
For low
The minimum discharge decreased significantly in years with lower than average maximum SWE and average preceding precipitation compared to years with higher than average SWE maximum and same amount of preceding precipitation (Fig. 8, bottom panels). For the Ova da Cluozza catchment, as an example, and considering only years with above-average SWE maximum, there is a 50 % probability that given an average preceding precipitation there will be a 7-day minimum discharge equal or higher than 107 % of its normal in July. In contrast, considering years with below-average SWE maximum, the 7-day minimum discharge will decrease to 75 % of its normal level. Similar changes were detected both in higher-elevation catchments and lower-elevation catchments, although in the latter this decrease is somewhat smaller.
The combined effect of snow and liquid precipitation on low flows was
analysed using “score plots”. In these plots the position of each
catchment is shown according to its average influence of snow and
precipitation on the 7-day minimum discharge separately for the period from
June to September (Fig. 9). Points located below the
Based on our results it seems that dependencies between predictors and
response variables may be connected to catchment properties and climate
drivers to some degree, such as elevation and thus maximum SWE and
The correlation of the dependencies of summer low flows on catchment elevation can be explained by lower air temperature at higher elevation and thus more snow accumulation and may be supported by results of Birsan et al. (2005) and Staudinger et al. (2015) in Swiss catchments. Staudinger et al. (2015) showed that higher-elevation and steeper catchments were less sensitive to droughts mainly because of an increasing snow influence but also because of potentially larger storages for the higher-elevation catchments of the selection. Our results showed that this sensitivity might increase with decreasing SWE especially in the highest-elevation catchments.
The elevation was also related to the memory effect of individual catchments which was generally longer for the highest-elevation catchments than for middle- or low-elevation catchments. However, even with the highest-elevation catchments, we did not find any significant correlations of snow and minimum discharges in October and later. In contrast, Godsey et al. (2014) found significant correlations even with the previous year's snowpack for some catchments in the western United States.
A longer memory effect in catchments with higher elevation is not only connected to higher snowpack accumulations but also to the simple fact that snowmelt occurs later in spring and persists longer compared to catchments at lower elevations (often until late spring or even early summer). The dependence of the DOY with minimum 7-day discharge on the DOY of maximum SWE was confirmed in our study. Similar dependences were found also in Whitaker et al. (2008), using the timing of the first significant snowmelt event instead of the DOY of maximum SWE. A negative trend in the number of days with discharge below a specified threshold in the case of increasing maximum SWE was proved. A 25 % quantile of discharge from May to October was used in this study. A 10 % quantile was also tested and found to have only minor impact on the results.
As documented by Beaulieu et al. (2012) in British Columbia, snow from headwater parts of catchments contributes significantly to base flow in lower parts of the catchments during summer. Earlier snowmelt onset and thus decrease of minimum streamflow has been observed (Jefferson, 2011) and a further shift of snowmelt towards earlier spring is predicted (Barnett et al., 2005; Bavay et al., 2009; Hanel et al., 2012; Godsey et al., 2014; Blahusiakova and Matouskova, 2015).
Inter-annual variability of seasonal actual evapotranspiration (ET),
seasonal precipitation (
Of course there are many other factors together with snow in winter that influence and can explain low-flow conditions in summer. For instance, evapotranspiration may change from year to year. However, water balance component estimates for the entire Switzerland during the last 100 years show that annual precipitation and runoff vary far more than evaporation (Hubacher and Schädler, 2010). Additionally, our results using data from PREVAH model proved that the inter-annual variability of ET in the warm season is much a lower than inter-annual variability of precipitation and maximum SWE (Fig. 10). Thus, the ET had a smaller impact on the inter-annual variability of 7-day minimum discharges compared to the impact of precipitation and maximum SWE. Although this is a relatively simple approach, it shows that ET likely is a less useful predictor to explain inter-annual variability of low flows compared to the other predictors. However, the within-year variation of ET is of course important to explain the occurrence of low flows (Teuling et al., 2013).
We tested also drainage density to account for landscape draining properties (e.g. Tague and Grant, 2004) and found significant correlation of low-flow sensitivity to the change of maximum SWE (Table 3). Draining properties together with catchment storage properties may help understanding the process causality leading to summer low flow. Combining this kind of catchment properties with the snow information might be useful for prediction. However, in this study we explored and quantified the general dependency and sensitivity to winter snowpack in humid regions.
However, the results presented in this study do not explain the process causality in detail. It means that we quantified the relations based on data we used, but process-based understanding at the catchment scale is limited and has to be further investigated.
The influence of snow conditions on summer low flow will likely decrease due
to predicted air temperature increase during winter and thus the decrease of
We did not explore possible impact of climate change on SWE and minimum
discharges, such as relations between possible warming in the cold season
and minimum discharges in the warm season, or explore whether any SWE decrease
will occur related to it. However, if a SWE decrease occurs, then the
same percentage decrease of SWE in higher-elevation catchments will result
in a stronger percentage decrease of minimum discharges than the same
percentage decrease of SWE in lower-elevation catchments (see Fig. 6 showing
the elasticity). This means that any interpretation going towards the possible
effect of increasing air temperature might be misleading since we did not
explore the relation between temperature increase and SWE decrease (e.g. due
to decrease in snowfall fraction). This interpretation would be difficult
especially in the highest-elevation catchments where air temperature in the
cold period is well below 0
In this study we looked at catchment mean elevations and for some analyses
we also classified catchments as high-, mid- or low-elevation catchments. In
practice this might be oversimplified as there is generally a large gradient
of precipitation and
The correlation between minimum discharge and maximum SWE considering years with little rain was higher than in years with a lot of rain. Low flows are usually higher during years with above-average snow conditions. Even in the case of low antecedent precipitation, low flow was higher than in years with below-average snow conditions. Therefore, snow plays an important role, although below-average snow conditions do not necessarily indicate below-average low flows. Preceding precipitation seems to be more important in this case. Because of the combined effect of snow and summer precipitation on summer low flows, snow-related parameters cannot fully explain the inter-annual variability of low flows in humid regions as documented by Godsey et al. (2014) for strongly seasonal regions even for the highest-elevation catchments. Nevertheless, most of detected trends in our study were significant at less than the 0.05 level, showing the significant influence of snow on summer low flows.
The decrease of maximum SWE and snowfall fraction increased the relative importance of rain during summer. Our results showed that the continuous decrease of maximum SWE and snowfall fraction in the future might increase the sensitivity of catchments at mid- and high elevations to hydrological droughts. This conclusion is in accordance with results of Birsan et al. (2005).
We chose the period from May to September to show the changing importance of snow contribution to low flows in different catchments, both at lower and higher elevations. We also tested the effect of maximum SWE on summer minimum discharge (June–August, not shown in this paper). The results for most of catchments are very similar to existing relations calculated for August because most of summer minimum discharges occur in August. Clearly we see the lowest summer flow as a compelling response variable, given the water management interest and possible issues connected to it. However, for the development of the role of snow compared to liquid precipitation this one response variable is not sufficient.
Our results do not provide a general answer to the question of whether snow storage is more important than precipitation. Due to moderate humid climate in Switzerland with precipitation almost equally distributed in a year (opposite to western United States), the aim was to show the combined effect of snow and liquid precipitation and their changing role in time (in different months) and in catchments with different elevation. Summer precipitation in Switzerland is relatively higher than summer precipitation in the western United States and, as shown in our study, summer precipitation dominates over the effect of snow, especially with an increasing time from the snowmelt period and with decreasing elevation. This combined effect explains the contrary results for the western United States and Switzerland.
We used winter precipitation as a predictor and we expected similar results
as using maximum SWE as predictor. Winter precipitation (from November to
April) is highly correlated with SWE and we expect increasing mutual
correlation for higher-elevation catchments with higher
Our results quantified the effect of snow on minimum discharges when liquid precipitation is below average (or opposite, when SWE is below average) as documented in Fig. 8. This could increase the reliability of predictions of minimum discharge during summer. Additionally, we provided information about sensitivity of low flow in individual catchments to changes in maximum SWE using prediction intervals showing the 50 % probability as well as prediction bands enabling the prediction of future observation. With this approach, it was possible to quantify not only the effect of snow storage on minimum discharge, but also on other low-flow parameters, such as length of the period with minimum discharge, DOY of minimum discharge occurrence and number of days below a specified runoff threshold.
We used new SWE data covering the entire Switzerland. From our study we see a big potential to use these data for instance to regionalize the catchment sensitivity and the length of snowmelt contribution to runoff in poorly gauged areas.
In this study we described and quantified the influence of winter and spring snow conditions on summer low flows in 14 Swiss Alpine and pre-Alpine catchments over the last 42 years. Specifically, we investigated the memory effect related to snow influence in runoff and the sensitivity of the catchments to low-flow reduction due to any change of snowpack. The main conclusions were the following.
Snowmelt significantly affected minimum discharge in May to September (with decreasing importance) in study catchments higher than 2000 m a.s.l., up to, in July and August in mid-elevation catchments and only in June and July in the lowest-elevation catchments. The sensitivity of minimum discharges to maximum annual SWE was higher for catchments at higher elevation when looking on relative changes in minimum discharge and maximum SWE.
Low flows occurred later in the year for years with above-average snow accumulations. However, the Spearman rank correlation between maximum SWE and the DOY with 7-day minimum discharge was significant only in some higher- and middle-elevation catchments (mean catchment elevation above 1500 m a.s.l.). The differences between the catchments were determined by both higher maximum SWE and later melt-out in the higher-elevation catchments.
Maximum SWE showed the best prediction ability from all winter-related predictors used in this study especially for higher-elevation catchments of our selection. Applicable results were achieved also with winter precipitation (November–April). However, winter precipitation is not suitable to describe the catchment sensitivity and they are not corrected for undercatch. Thus, using maximum SWE is recommendable for sensitivity studies.
Snow and summer precipitation had a combined effect on summer low flows, and snow accumulation alone cannot explain the inter-annual variability of low flows even in high-elevation catchments. Snow was a better predictor for the variability of low flows when only years with lower than average preceding precipitation were considered. However, even if both snow and liquid precipitation are considered, there is still some portion of inter-annual variability which cannot be explained by these two predictors.
Summer low flows in the study catchments were significantly sensitive to any SWE changes. Although our study did not address climate change impacts explicitly, a reduced effect of snow on late summer low flows in the highest- and middle-elevation catchments can be expected due to the predicted decrease of snowfall fraction in the future. As a consequence, the sensitivity of catchments at mid- and high elevations to meteorological droughts might increase.
M. Zappa and J. Seibert initiated the project. M. Jenicek developed the methodology (with contributions of Jan Seibert) and performed all analyses. M. Zappa, M. Staudinger and T. Jonas prepared input meteorological data used for analyses. M. Jenicek prepared the paper with substantial contributions from all co-authors.
Support from the Swiss National Research Programme Sustainable Water Management (NRP 61, project DROUGHT-CH) and the Czech Science Foundation (GACR 13-32133S, project Headwaters) is gratefully acknowledged. The authors also thank SCIEX – Scientific Exchange Programme NMS.CH for the support of the first author during his postdoc stay at the University of Zurich. The contribution of M. Zappa was supported by the Swiss National Science Foundation SNF through the Joint Research Projects (SCOPES) Action (Grant IZ73Z0_152506). Many thanks to Tracy Ewen for improving the English of the paper. Edited by: A. Van Loon