The young water fraction (Fyw*), defined as the fraction of catchment outflow with transit times of less than 2–3 months, is increasingly used in hydrological studies that exploit the potential of
isotope tracers. The use of this new metric in catchment intercomparison
studies is helpful to understand and conceptualize the relevant processes
controlling catchment functioning. Previous studies have shown surprising
evidence that mountainous catchments worldwide yield low Fyw*. These low values have been partially explained by isolated hydrological processes, including deep vertical infiltration and long groundwater flow paths. However, a thorough framework illustrating the relevant mechanisms leading to a low Fyw* in mountainous catchments is missing.
The main aim of this paper is to give an overview of what drives Fyw* variations according to elevation, thus clarifying why it generally decreases at high elevation. For this purpose, we assembled a data set of 27 study catchments, located in both Switzerland and Italy, for which we calculateFyw*. We assume that this decrease can be explained by the groundwater storage potential, quantified by the areal extent of Quaternary deposits over a catchment (Fqd), and the low-flow duration (LFD) throughout the period of isotope sampling (PoS). In snow-dominated systems, LFD is strictly related to the snowpack persistence, quantified through the mean fractional snow cover area (FSCA). The drivers are related to the catchment storage contribution to the stream that we quantify by applying a cutting-edge baseflow separation method to the discharge time series of the study sites and by estimating the mean baseflow fraction (Fbf) over the PoS.
Our results suggest that Quaternary deposits could play a role in modulating
Fyw* elevation gradients via their capacity to store groundwater, but subsequent confirmation with further, more detailed geological information is necessary. LFD indicates the proportion of PoS in which the stream is sustained and dominated by stored water coming from the catchment storage. Accordingly, our results reveal that the increase of LFD at high elevations, to a large extent driven by the persistence of winter snowpacks and the simultaneous lack of a liquid water input to the catchments, results in lower Fyw*. In our data set, Fbf reveals a strong complementarity with Fyw*, suggesting that the latter could be estimated as Fyw*≃1-Fbf for catchments without stable water isotope measurements.
As a conclusion, we develop a perceptual model that integrates all the results of our analysis into a framework for how hydrological processes
control Fyw* according to elevation. This lays the foundations for an improvement of the theory-driven models.
Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen ForschungPP00P2_157611PRIN MIUR 2017SL7ABC_005Introduction
Mountainous alpine catchments are often assumed to generate high shares of
rapid surface or subsurface runoff due to the presence of exposed bedrock
and steep landscapes. Consequently, the role of groundwater storage in
high-elevation catchments has been often neglected (Hayashi, 2020). On the contrary, multiple worldwide studies quantified a considerable groundwater input to streamflow in high mountain catchments using tracer or water balance methods (Somers and McKenzie, 2020). Several studies from the Rocky Mountains and Andes show that, on average, about 47 % of groundwater annually sustains the streamflow (Saberi et al., 2019; Somers et al., 2019; Carroll et al., 2018; Harrington et al., 2018; Cowie et al., 2017; Baraer et al., 2009, 2015; Gordon et al., 2015; Frisbee et al., 2011; Liu et al., 2004; Clow et al., 2003). Similar percentages, 49 % and 48 %, are also found in the Himalayas and the Alps, respectively (Chen et al., 2018; Engel et al., 2016; Käser and Hunkeler, 2016; Williams et al., 2016; Wilson et al., 2016; Andermann et al., 2012). It is well known that the water is stored longer than a year or a few years and that stored water plays a key role in streamflow generation processes (McDonnell, 2017; Jasechko, 2019). The study of water age has implications for predicting the timing of nutrient cycles and pollutant transport, since water age and solute dynamics are closely coupled (Li et al., 2021). Nevertheless, water age quantification is not straightforward.
Kirchner (2016a, b) proposed a new metric to quantify the share of catchment outflow with transit times lower than roughly 0.2 years or 2–3 months: the young water fraction. This metric can be conveniently inferred from the dampening effect that a catchment has on the seasonal cycle of stable water isotopes in precipitation, i.e., by estimating the ratio of the amplitudes of the seasonal cycles of stable water isotopes in streamflow and in precipitation (Kirchner, 2016a). In this method, the seasonal cycle of stream water isotope measurements is modeled using a sine wave that can be flow weighted, using the discharge measured at the moment of sampling as a weight, or not (von Freyberg et al. 2018). Isotopes measured in precipitation can be modeled with a sine function weighted according to the volume of precipitation to reduce the influence of low-precipitation periods and to account for temporally aggregated rainfall samples (von Freyberg et al., 2018). Flow-weighted fits to the seasonal tracer cycles predict the flow-weighted average young water fraction (Fyw*) in streamflow, while unweighted fits to the seasonal tracer cycles predict the unweighted one (Fyw) (Kirchner, 2016b). Gallart et al. (2020a) recently highlighted the advantages of the flow-weighted analysis to compensate for subsampled high-flow periods, thus reducing the underestimation of the young water fraction. Hereafter, we will use the symbol “*” for referring to a flow-weighted variable, in order to be consistent with previous studies (von Freyberg et al., 2018; Gallart et al., 2020a).
Fyw* is increasingly used in hydrological studies because it has the advantage of being free from the aggregation errors inherent to mean transit time (MTT) estimates obtained through the classical convolution approach (Kirchner, 2016a). Even more so, Fyw* is an
informative descriptor of catchment hydrological functions, of nutrients
cycles and of pollutant transport (Stockinger et al., 2019; Benettin et al., 2017; Jasechko et al., 2016; Xia et al., 2023). For these reasons, this new metric is useful for catchment intercomparison studies to find what are the main hydroclimatic and landscape characteristics that drive the transit times of water lower than a threshold age, which varies from about 2 to 3 months. Indeed, previous work has tried to study the relationship between Fyw* and catchment characteristics. von Freyberg et al. (2018)
found that young water fractions of 22 Swiss catchments are significantly
positively correlated with selected hydroclimatic indices and with the
fraction of saturated area, suggesting that Fyw* depends on catchment wetness, which promotes rapid flow paths. Interestingly, von Freyberg et al. (2018) found a statistically significant positive correlation with elevation after removing the five snow-dominated catchments, which expressed the smallest Fyw* (von Freyberg et al., 2018). Likewise, Lutz et al. (2018) estimated Fyw* for 24 catchments in Germany and found the smallest values for higher-elevation sites. These results are partially consistent with those of Jasechko et al. (2016), who based on the analysis of 254 watersheds worldwide discovered a reduction of Fyw* in mountainous, steeper terrains. This could be related to deep vertical infiltration caused by fractures generated by high rock stress in complex terrain morphologies or by freely draining soils (i.e., cambisols and luvisols), both associated with high-elevation environments
(Lutz et al., 2018; Jasechko et al., 2016; Gleeson et al., 2014). In addition, topographic roughness increases flow path and, correspondingly, transit time (Gleeson and Manning, 2008; Frisbee et al., 2011; Jasechko et al., 2016). Despite these studies, there is still a lack of a unified framework of how the variation among mountainous catchments results in less young water at high elevation.
An early example from the Swiss Alps showed that high celerity originates
from massive meltwater infiltration that pushes out groundwater reserves:
streamflow following snowmelt is older than meltwater infiltrated in the
current year (Martinec, 1975). The resulting effect on water partitioning between the surface and the subsurface should be analyzed considering the temporal concentration of water input during the snowmelt period, but this remains largely unexplored (Rey et al., 2021). Despite this lack of studies on water partitioning during snowmelt, several studies have
demonstrated the pivotal role of snowmelt in recharging groundwater during
summer in high-elevation environments (Hayashi, 2020; Cochand et al., 2019; Du et al., 2019; Flerchinger et al., 1992).
From a water modeling perspective and thus from a water age perspective,
snowpack storage and groundwater storage can be considered a single entity:
they both constitute catchment storage. Therefore, the analytical estimation
of Fyw* must reflect this “conceptual” decision of whether to consider the snowpack storage to be part of the catchment storage. This point has been previously addressed by von Freyberg et al. (2018). If total precipitation is considered as catchment input (direct input case), the snowpack is implicitly considered to be part of the catchment storage, and Fyw* results from the combination of snowpack and subsurface storage. In this direct input case, Fyw* is computed from the amplitudes of the seasonal cycles of stable isotopes of water in precipitation (AP) and streamflow (As*). If total liquid water input (composed of rainfall and snowmelt before sometimes called equivalent precipitation) is considered to be catchment input, Fyw* is computed based on the amplitudes of the cycles in equivalent precipitation (APeq) and in streamflow (As*). This Fyw* value then results from subsurface storage alone, since snowpack storage is excluded from the catchment storage (von Freyberg et al., 2018). If Fyw* is estimated using a direct input setting (i.e., total precipitation directly), Fyw* is expected to be smaller, since the catchment storage is larger (von Freyberg et al., 2018). Also, Ceperley et al. (2020) investigated the role of water input from snow in Fyw* estimation, concluding that the low values in high alpine snow-dominated catchments result from a combination of snow cover effects and the storage in the subsurface. In the present work, the main aim is not to address how the snowpack affects Fyw* estimation in a single catchment, as this was treated previously (von Freyberg et al., 2018; Ceperley et al., 2020), but to investigate the hydrological processes (also related to the snowpack storage) that lead to variations in Fyw* between catchments located at different elevations with a focus on high-elevation alpine catchments.
Some authors have revealed the possibility of Quaternary deposits (e.g.,
talus, moraine and alluvium) to store groundwater in high-elevation alpine
catchments (Arnoux et al., 2021; Hayashi, 2020; Christensen et al., 2020). The stored water in these deposits can in fact sustain streamflow during the low-flow period (Hayashi, 2020; Arnoux et al., 2021), as supported by the strong positive correlation found by Arnoux et al. (2021) between the fraction of Quaternary deposits and the winter flow index (a low-flow indicator reflecting the groundwater contribution to the stream) for 13 alpine catchments. During winter, the period without liquid water input can last 6 months or more in high-elevation catchments. The occurrence of such long periods of low flows hints towards important amounts of stored water (or old water) that are well connected to the stream network, and it thereby remains accessible throughout the low-flow period (Somers et al., 2019).
To further discuss the role of low flow in Fyw* estimation, let us first consider that Fyw* can be theoretically estimated based on the flow-weighted average of young water fractions (Kirchner, 2016b):
Fyw*=AS*AP≃∑i=1nQtiFywti∑i=1nQti,
where n is the number of time steps (e.g., days) in the period of isotope
sampling, PoS, Q(ti) is the discharge at the time ti (e.g., daily
discharge) and Fyw(ti) is the young water fraction at the time ti (e.g., daily young water fraction). As is clear from Eq. (1), Fyw* becomes low if either Fyw(ti) is low for high flows or if Fyw(ti) is very low for many time steps or both. The low-flow periods correspond to the recession periods in which there is no new rainfall or meltwater input in the catchments. Thus, during these periods, the catchment storage releases stored water (or old water) to the stream sustaining the streamflow (Hayashi, 2020). Thus, we can anticipate that low Q(ti) values imply low Fyw(ti) values. Accordingly, the proportion of the low-flow period during a specified time window should reduce the amount of young water reaching the stream during that time window. Nevertheless, the Fyw(ti) is higher during high-flow (wet) periods (von Freyberg et al., 2018; Wilusz et al., 2017; Gallart et al., 2020b). Thus, the overall effect of the proportion of low-flow and high-flow periods upon Fyw* remains a priori unclear. It is however tempting to think that the duration
of low-flow period or the share of baseflow could explain Fyw* variations at different elevations (since both low-flow duration and the share of baseflow change with elevation). In addition, in high-elevation,
snow-dominated catchments, the persistence of the snowpack is the main
driver of the long low-flow duration, since the low-flow period at high
elevation corresponds to the presence of the seasonal snowpack (corresponding to an absence of liquid water input), while the high-flow period is
generally snowmelt driven. Such snowmelt generally occurs in late spring or
summer, and it is likely to be older than 2–3 months (because the peak snow
fall occurred 3 months earlier). As a result, summer discharge mainly
consists of old water, either of current snowmelt that reaches the stream
via faster surface or subsurface flow paths, or old snowmelt (main component
of groundwater storage) pushed out in the stream by infiltrated rainwater or
meltwater. Part of the snowpack can release young water, but this is a minor
component in catchments with a seasonal snowpack. In contrast, in catchments
with an ephemeral snowpack, it is common to observe intermittent winter
snowmelt that is likely younger than 2–3 months: snowmelt is temporally
close to snowfall. In this case, streamflow receives relatively more young
water from short-lived snowpacks. However, it is still unclear if seasonal
or ephemeral snow cover dynamics can affect the Fyw* (Ceperley et al., 2020).
An innovative focus of our work is on variables that were not previously
considered for explaining elevation gradients of young water fractions. We
specifically exclude catchment size, annual precipitation, bedrock porosity,
pasture cover and open water cover that have been discussed and shown to have
little correlation in the work of Jasechko et al. (2016).
A special case in terms of explanatory variables is mean annual precipitation: Jasechko et al. (2016) did not observe any significant correlation between the Fyw* and annual precipitation in
their worldwide study. Lutz et al. (2018) found that Fyw* decreases with increasing mean annual precipitation, based on 24 catchments in Germany. In contrast, in the relatively wet rainfall-dominated and hybrid catchments studied by von Freyberg et al. (2018), Fyw* was shown to increase with mean monthly precipitation and correspondingly also with elevation. In their study, discharge (unsurprisingly correlated with precipitation) was considered as a proxy of catchment wetness, which favors rapid flow paths and thereby increases Fyw* (von Freyberg et al., 2018). In snow-dominated systems, the use of mean annual precipitation as a proxy for catchment wetness could be misleading because the seasonal snowpack leads to a very dry period of the year despite the high solid water input. In other words, the temporal concentration of the liquid water input is the relevant variable. Indeed, the saturation of the system (i.e., high wetness conditions) can be observed also when the annual
precipitation is low if a large volume of water (stored in the snowpack) is
released in a relatively concentrated time interval. Indeed, despite the
fact that precipitation and, correspondingly, discharge are higher in
snow-dominated than in rainfall-dominated catchments, Fyw* is
generally lower in snow-dominated systems that are potentially wetter than
rainfall-dominated ones. This suggests that the precipitation can only
partially explain the variations of Fyw* and that other variables should be put under observation.
Accordingly, we omit here total annual precipitation as an explanatory variable
of low Fyw* in snow-dominated catchments (but we consider
precipitation for rainfall-dominated and hybrid catchments) and study a new
set of hydrological variables to gain new insights into Fyw* along elevation gradients: the fraction of Quaternary deposits (Fqd), the mean fraction of baseflow (Fbf), the low-flow duration (LFD) and the mean fractional snow cover area (FSCA), defined in detail in Sect. 3.2 and 3.3. We first describe the data set (Sect. 2). Then, we present the Fyw* estimation method (Sect. 3.1) followed by the correlation analysis of the selected hydrological variables with the estimated Fyw*, and we bring these results back into the ongoing scientific discussion of Fyw* (Sect. 4.2–4.6).
Study sites
We analyze 27 study catchments located both in Switzerland and Italy
integrating observations from multiple published data sets (25 catchments)
with new additional observations (2 catchments) (Fig. 1). Geomorphological and hydroclimatic characteristics of the study sites are reported in
Table 1.
Catchment geomorphological and hydroclimatic characteristics.
The catchment area and mean slope are directly calculated in Google Earth
Engine. For the slope calculation, we use the Shuttle Radar Topography
Mission (SRTM) DEM (Farr et al., 2007). We obtained mean elevation and precipitation information of the existing data set directly from published papers (von Freyberg et al., 2018; Ceperley et al., 2020). Discharge (Q), precipitation (P) and isotopic composition (δ18O) data are all referred to the period of sampling (PoS) indicated in this table. The letter in brackets in the first column indicates the hydroclimatic regime: (R) is rainfall dominated, (H) is hybrid and (S) is snow dominated.
Location of the 27 study catchments with indication of the hydroclimatic regime.
Specifically, we assembled the 22 Swiss catchments studied by von Freyberg
et al. (2018) with the three alpine catchments investigated by Ceperley et al. (2020) (Vallon de Nant, Noce Bianco at Pian Venezia and Bridge Creek catchment) into a single data set. Hereafter, we refer to these catchments with the ID reported in the above-mentioned published papers (Table 1). We also consider two additional high-elevation catchments located near the Nivolet Pass (Valsavaranche, Aosta Valley, Italy) (Gisolo et al., 2022). In this alpine environment, we monitor the mainstream, called “Dora del Nivolet”, and a secondary river called “Source”. Hereafter we refer to these catchments with the IDs DOR and SOU, respectively. A detailed description of the DOR and SOU catchments is reported in the Supplement.
The von Freyberg et al. (2018) data set includes catchments with areas between 0.7 and 351 km2 and mean elevations between 472 and 2369 m a.s.l. With the five catchments added here, the complete data set includes catchment areas between 0.14 and 359 km2 and spans mean elevation between 472 and 3049 m a.s.l. The mean monthly precipitation ranges between 61.3 and 168.7 mm per month, while mean discharge ranges between 28.6 and 138.9 mm per month. The mean slope ranges from 4 to 34∘, and our study sites reveal an increase of steepness with elevation (Fig. 2a). Precipitation increases with elevation until 1500 m a.s.l., above which it decreases (Fig. 2c), highlighting a change of precipitation regime as described by previous studies (Santos et al., 2018). The five catchments added to the initial data set of von Freyberg et al. (2018) allow the analysis to explore the high-elevation regions (mean elevation > 1500 m a.s.l.) that were previously poorly represented. Most of the study catchments reveal a sedimentary bedrock but dolomitic and metamorphic bedrocks, characteristic of high-elevation sites, are also included in our data set. Moreover, the presence of unconsolidated Quaternary deposits is widespread among our study catchments: only two catchments (BCC and SOU) do not reveal this type of geology. The complete data set now explores case studies from the Swiss plateau and pre-alpine area; from the Jura; and from five different alpine regions, including the northern part of the Swiss Alps, the southern Swiss Alps (Alpi Ticinesi), the western Italian Alps (Alpi Graie), the Rätische Alps and the Dolomites. Overall, this represents a good range of geologies as well as of climatic conditions.
(a) Mean slope against mean elevation. Vertical bars represent the mean slope standard deviation, horizontal bars represent min–max elevation range, (b) boxplot of the mean slope values, (c) mean precipitation and discharge against elevation, and (d) boxplots of the mean precipitation and discharge values. Here and later: the boxplots show the median and the interquartile range (IQR), the whiskers are defined as the IQR multiplied by 1.5 and outliers are plotted with red “+” markers.
Boxplots of mean monthly flow for all the study catchments grouped
according to their flow regime (rainfall dominated, hybrid, snow dominated).
Colored areas represent the monthly flow of each study catchment belonging
to the relative regime.
In order to be consistent with previous studies (von Freyberg et al., 2018; Staudinger et al., 2017), we classify the 23 Swiss catchments according to the hydroclimatic regimes proposed by Staudinger et al. (2017) which group the regimes defined by Weingartner and Aschwanden (1992) in three categories:
rainfall dominated (R), hybrid (H) and snow dominated (S). For the four
Italian catchments, where the aforementioned classification schemes cannot
be rigorously applied, we use that proposed by Stoelzle et al. (2020). This classification scheme is based on mean and maximum catchment elevation, periods of typical low flow, snow onset and beginning of snowmelt and was already used by Stoelzle et al. (2020) to classify catchments outside the Swiss borders (e.g., German catchments). According to this classification scheme, the four Italian catchments (DOR, SOU, BCC and NBPV) are all categorized as snow dominated (S). The classification of BCC is also consistent with the one given in a previous study without considering the application of a formal classification scheme (Penna et al., 2016). Across the three considered streamflow regimes, a shift of the monthly hydrograph peak (computed using discharge data in the PoS) from winter to summer months is observed (Fig. 3): this flow peak shifting is a clear sign of the increasing predominance of snowmelt in the streamflow generation processes. Our data set includes NBPV, whose area is 42 % glacier covered and consequently exhibits a characteristic glacier-dominated streamflow regime with a monthly peak in late summer (Zuecco et al., 2019; Carturan, 2016). NBPV has been classified as snow dominated following the Stoelzle et al. (2020) classification scheme.
Nevertheless, its characteristics suggest it may belong to a fourth category
of glacier-dominated catchments. Unfortunately, this category has not been
considered by the aforementioned classification scheme, and the definition of
the classifiers for a new category is outside the scope of this work. In
this catchment, the effect of glacier melt on Fyw* cannot be neglected, and this was partially discussed by Ceperley et al. (2020). In our data set, also the Dischmabach (DIS) and the Vallon de Nant (VdN) catchments are 2 % and 3 % glacier covered, but we assume that the effect on Fyw* is negligible when compared with that of the seasonal snowpack.
Material and methodsYoung water fraction estimation from seasonal cycles of stable water isotopes in precipitation and stream water: the “direct” input
Kirchner (2016a) designed the young water fraction as the proportion of the transit time distribution younger than a threshold age (τyw). By assuming that the transit time distribution mathematical form is the regularized lower incomplete gamma function for all the study catchments, the theoretical young water fraction (FywT) can be expressed as
FywT=Pτ<τyw=∫0τywτα-1βαΓ(α)e-τβdτ,
where α and β are the shape and scale factor, respectively.
By using thought experiments, Kirchner (2016a) has demonstrated that for a given shape factor α (ranging from 0.2 to 2) and across a wide range of scale factors β, the theoretical young water fraction can be accurately predicted by the amplitude ratios of seasonal sine curves fitted to stream water and precipitation isotope values by considering a τyw of 2–3 months. Operatively, we model seasonal isotope (e.g., δ18O) cycles in stream water and precipitation as reported in Eqs. (3a) and (3b):
3aδ18OS(t)=ASsin2πft-φS+kS,3bδ18OP(t)=APsin2πft-φP+kP,
where δ18O (‰) is the isotopic composition of water sampled at the time t (expressed in decimal years), A (‰) is the amplitude of the seasonal isotope cycle, φ (in radians, with 2π rad = 1 year) is the phase, f (yr-1) is the frequency and k (‰) is the vertical offset of the seasonal isotope cycle. The subscript S refers to stream water, while the subscript P refers to precipitation. The sine wave is fitted to the isotopes measured in precipitation weighted according to the volume of precipitation, reducing the influence of low-precipitation periods and accounting for temporally aggregated rainfall samples (von Freyberg et al., 2018); the sine fit of stream water isotope measurements can be discharge-weighted, using the discharge measured at the moment of sampling as weights, or not (von Freyberg et al., 2018). The sine curves of Eqs. (3a) and (3b) are fitted on the isotope measurements using the iteratively re-weighted least squares (IRLS) method (for reducing the influence of outliers), which leads to estimates of A, φ and k parameters. A function for performing a sine fit using IRLS, based on the IRLS function made available by Kirchner and Knapp (2020), is available in the Supplement.
Accordingly, depending on the unweighted or the flow-weighted fit, an
unweighted amplitude (AS) or a flow-weighted amplitude (AS*) can be obtained, respectively. Such amplitudes can be used to calculate the time-weighted (Eq. 4a) or the flow-weighted (Eq. 4b) young water fractions (Fyw or Fyw*, respectively) via the “amplitude ratio approach”:
4aFyw=ASAP,4bFyw*=AS*AP.
Gallart et al. (2020a) highlighted the advantages of the flow-weighted analysis (generally yielding AS* greater than AS) to compensate for subsampled high-flow periods, which would otherwise lead to a young water fraction underestimate. Accordingly, in this work, we calculate the flow-weighted young water fractions for all the study catchments by applying Eq. (4b). We obtain the standard errors (SEFyw*) of the estimated Fyw* using Eq. (5):
SEFyw*=1APSEAS*2+AS*AP2⋅SEAP2,
where SEAS* and SEAP are the standard errors of the regression coefficients AS* and AP. The analytical choice of using the amplitude (AP) fitted to precipitation isotopes, instead of the amplitude (APeq) fitted to the equivalent precipitation (i.e., rain plus snowmelt) isotopes, for estimating Fyw* implies that the snowpack (and/or the glacier) is considered as part of the catchment storage. Thus, the damped seasonal cycle observed in the stream is given by the mixing of precipitation with snowpack (and/or the glacier) and subsurface storage (the last two considered as a single entity), as illustrated in Fig. 4.
Schematic representation of the “direct input” approach for
estimating Fyw*. Light blue arrows indicate that meltwater coming from the snowpack preferentially infiltrates. The term C refers to the isotopic composition. The subscript P refers to precipitation, S refers to stream, SP refers to snowpack, Sub refers to subsurface storage and St refers to catchment storage.
Since we have assumed that the transit time distribution belongs to the family of gamma distributions, we can determine the parameter α of such a distribution by solving the following implicit expression for α
(Eq. 6):
φS*-φP=αarctanAS*/AP-2α-1.
We optimized to find the best solution of Eq. (6) using the best-fitting
parameters (φS*, φP, AS*, AP) with α spanning a wide interval between 0.01 and 20. For the relevant math, the reader is referred to Kirchner (2016a). We estimate the uncertainty of α assuming all the fitted parameters having a Gaussian distribution with a standard deviation equal to the regression error (Eqs. 3a and 3b). We generate 1000 random samples of such parameters, and we estimate the optimal α for each parameter set and then compute their standard deviation.
As in past studies (Gallart et al., 2020a; Lutz et al., 2018), we use the parameter α to estimate τyw with the following second-order polynomial fit (Kirchner, 2016a):
τywT≈0.0949+0.1065α-0.0126α2,
where T is the period of the tracer cycle (for a seasonal cycle, T=1 year). We estimate the uncertainty of τyw by calculating the standard deviation of the threshold ages computed using the 1000 optimal
α values previously obtained. A code for estimating α and
τyw with their uncertainties is made available in the Supplement.
The comparison of Fyw* among different catchments is potentially subject to a bias given by possibly different threshold ages ranging between 2 and 3 months. Accordingly, we couple each Fyw* with the corresponding τyw to illustrate what the term “young” means for each study catchment.
Snow cover persistence quantified through the mean fractional snow cover area (<inline-formula><mml:math id="M164" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi mathvariant="normal">SCA</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>)
In this paper, we quantify the snowpack persistence by calculating the mean
fractional snow cover area (FSCA). It is calculated for each catchment over the period 1 October 2017–30 September 2021 (hereafter defined as PoC, i.e., period of calculation) by using the collection of Sentinel-2 L2A satellite images available in Google Earth Engine (Gorelick et al., 2017). Temporally, this relatively recent satellite has increased the visitation frequency to a subweekly temporal resolution and increased the spatial resolution to 20 m for snow cover (Gascoin et al., 2019). High temporal resolution makes Sentinel-2 images preferable to Landsat images, which are available only once every 16 d and whose total number is often further reduced because of cloudiness (Hofmeister et al., 2022). The PoC generally differs from the PoS for the 27 study catchments. This is because Sentinel-2 L2A satellite images are not available before March 2017. For each image available in the PoC, we calculate the normalized difference snow index (NDSI) as suggested in the work of Dozier (1989):
NDSI=rgreen-rSWIRrgreen+rSWIR,
where rgreen is the reflectance in the green band (Sentinel-2 band 3) and rSWIR is the shortwave infrared reflectance band (Sentinel-2 band 11). We classify as snowy pixels those with an NDSI value > 0.4 (Dozier, 1989). Based on the pixel-by-pixel snow classification, we compute the snapshot fractional snow cover area (fSCA) according to Di Marco et al. (2020) and Hofmeister et al. (2022):
fSCA=NsnowNtot-Nclouds,
where Nsnow is the number of snow cover pixels according to the applied NDSI threshold method, Ntot is the total number of pixels within the catchment area and Ncloudsis the number of pixels classified as clouds and water bodies (Hofmeister et al., 2022). We identify the cloudy pixels directly using the Sentinel-2 band “Scene Classification Map”. We operatively calculate Nsnow, Ntot and
Nclouds using a Google Earth Engine code.
By using this procedure for calculating fSCA, we sometimes obtain
fSCA>1. The NDSI threshold method is generally able to distinguish between snow and no-snow pixels (Aalstad et al., 2020). Accordingly, clouds and snow have similar reflectance in the green band, but clouds highly reflect in the shortwave infrared band, while snow reflectance is low in this band. Thus, the Nsnow estimation is generally accurate. On the other hand, it is necessary to disregard the false positive pixels deriving from cloud detection (i.e., snow classified as clouds). If fSCA>1, we calculate fSCA as Nsnow/Ntot, since this is the only heuristic solution that guarantees no overestimation. Moreover, by looking at sample Sentinel-2 images during the summer periods for all the catchments, we impose fSCA=0 during July and August, since when fSCA≠0, this usually results from clouds falsely identified as snow: imposing fSCA=0 clearly leads to fewer errors (only missing occasional summer snowfall events of very shallow depth) than falsely accounting for (far more) frequent clouds. The NBPV catchment is an exception: we do not impose fSCA=0 during July and August, since it generally has snow over the glacier also during summer. Finally, we compute the mean fractional snow cover area (FSCA) for each catchment by averaging all fSCA values available for all snow images in the PoC, without interpolation between the time steps.
Fraction of Quaternary deposits, low-flow duration and the groundwater contribution to the stream
Similarly to Arnoux et al. (2021), we calculate the portion of the catchment area occupied by Quaternary deposits (Aqd) (available from government geological data sets) with respect to the total catchment area (A). Thus, we calculate the fraction of Quaternary deposits (Fqd) as reported by Eq. (10):
Fqd=AqdA.
Additionally, we use the same winter flow index (WFI) as Arnoux et al. (2021), as indicated by Eq. (11):
WFI=QNM7Qmean,
where QNM7 is the minimum discharge over 7 consecutive days during the winter period (from November to June) and Qmean is the mean annual discharge. We calculate it for the 27 study catchments during the PoS. To relate WFI to low flow, we apply the recent baseflow separation technique described by Duncan (2019) to the discharge time series of the 27 study catchments (within the PoS indicated in Table 1). In short, this method comprises a single backward pass through the data to fit an exponential master baseflow recession curve (Eq. 12a), followed by the single forward pass (Eqs. 12b and 12c) of the Lyne and Hollick (1979) algorithm. This allows the smoothing of the connection between segments of the master recession by simulating a gradual groundwater recharge during the runoff event (Duncan,
2019):
12aMti-1=Mti-ck+c12bQqti=kQqti-1+Mti-Mti-11+k212cQbfti=Mti-Qqti,
where M(ti), Qq(ti) and Qbf(ti) are the master recession value, the quick recession flow and the baseflow at time ti, respectively. In this study, we consider daily time steps (i.e.,
ti-ti-1=1 d). This method has two parameters: k is the recession
constant, c is a constant flow added to the exponential decay component. We
set the recession constant k=0.925 (Nathan and McMahon, 1990): we add no constant flow to the exponential decay (i.e., in terms of the method by Duncan (2019), c=0). A code with the implementation of the Duncan (2019) baseflow filter has been made available in the Supplement.
To express the catchment storage contribution to streamflow in a form that
is directly comparable to the Fyw*, we define the baseflow fraction (Fbf) as reported in Eq. (13):
Fbf=1n∑i=1nQbftiQti,
where Qbf(ti) is the baseflow (mm d-1) at the time ti (obtained as indicated by Eq. 12c) and Q(ti) is the discharge (mm d-1) at the time ti. We tested the uncertainty of k by drawing random samples (10 000) from a normal distribution spanning Nathan and McMahon's (1990) recommended range for k from 0.90 to 0.95, with a mean of 0.925 and a standard deviation equal to 25 % of the range. Thereby, we obtain 10 000 values of Fbf for each catchment of which we compute the standard deviation.
As introduced in Sect. 1, Fyw* can be low if the snapshot young water fraction Fyw(ti) is very low for many time steps. If we consider the discharge (Q) as a proxy for the catchment wetness, we can reliably assert that Fyw(ti) is low for low Q(ti). Thus, another important variable is the duration of the low-flow period. In this study, we define a low-flow period (TLow) as follows:
TLow=∀ti:QbftiQti≥0.85.
Thus, a low-flow period is defined here as a period when 85 % of the total
flow is composed of baseflow (i.e., baseflow dominated). Accordingly, we
define the low-flow duration (LFD) as the proportion of the time steps (e.g., days) in the PoS that can be considered as a low-flow period according to Eq. (14).
Sinusoidal cycles of both precipitation and streamflow fitted to the δ18O data (using the IRLS method) for six representative study catchments. Amplitudes (‰) and phases (year) are indicated in the figure. Please note that both δ18O data and sinusoidal cycles of precipitation and stream water are vertically shifted of kP and kS, respectively.
Results and discussion – towards a harmonious and exhaustive framework of the hydrological processes that drive the young water fraction variations with elevation
We present and discuss hereafter the Fyw* and τyw estimates (Sect. 4.1) and the identified relations between Fyw* and the studied explanatory variables (Sect. 4.2–4.5), followed by the perceptual model that describes the main processes driving the Fyw* variations with mean catchment elevation and that harmonizes our results with previous studies (Sect. 4.6).
Young water fractions (<inline-formula><mml:math id="M234" display="inline"><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mi mathvariant="normal">yw</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>) and corresponding threshold ages
Assembling Fyw* values determined by different authors who very likely used different source codes could possibly result in a bias. Indeed, differences in Fyw* among catchments could be driven by the different methods rather than the physical factors. Therefore, the same approach must be applied to all the study catchments to remove the bias introduced by the estimation method of Fyw*. For all the study catchments, sinusoidal cycles were fitted to both precipitation and stream water δ18O data by using the IRLS regression (results for six representative study catchments in Fig. 5; complete results in Fig. S2). We estimate Fyw* via Eq. (4b) by using the best-fitting amplitudes of seasonal cycles. The best-fit amplitudes (AP, AS*), phases (φP, φS*) and corresponding standard errors are reported in Table 2.
Summary table with all the relevant quantities estimated for the 27 study catchments.
The Fyw* values achieved in this study are consistent with published Fyw* values of Ceperley et al. (2020) and with direct-input Fyw* of von Freyberg et al. (2018). Accordingly, a two-sample Kolmogorov–Smirnov test accepts the null hypothesis that the new and past Fyw* estimates are from the same continuous distribution at the 5 % significance level. The Fyw* estimates are reported in Fig. 6a against the mean catchment elevation and are also listed in Table 2. Fyw* increases with mean catchment elevation until 1500 m a.s.l., which corresponds to the elevation above which all catchments are snow dominated (with NBPV detected as an outlier as will be discussed in Sect. 4.5). This pattern is also reflected by the median Fyw* within each flow regime: the median Fyw* is 0.13 for rainfall-dominated catchments, 0.29 for hybrid catchments and 0.12 for snow-dominated catchments. Such results are consistent with previous studies that have shown the tendency toward low Fyw* in mountainous catchments (Ceperley et al., 2020; von Freyberg et al., 2018; Lutz et al., 2018; Jasechko et al., 2016).
(a)Fyw* as function of mean catchment elevation. Points dimension is proportional to τyw(b), obtained with Eq. (14) of Kirchner (2016a), as function of the shape factor α.
Even though we remove the bias introduced by the estimation method of
Fyw*, the application of Eq. (4b) implicitly introduces another bias if we want to use Fyw* for intercomparison purposes (which is the goal of this work). By computing the amplitude ratio, we estimate FywT (Eq. 2), without defining a corresponding τyw (which can be estimated via Eq. 7). Therefore, part of the scatter of Fyw* between catchments
might be because of different τyw rather than physical factors, also if the τyw is expected to vary modestly between 2 and 3 months (Kirchner, 2016a). Accordingly, past studies estimated Fyw* using the amplitude ratio approach without information about the corresponding τyw (Stockinger et al., 2019; von Freyberg et al., 2018; Jasechko et al., 2016). Nevertheless, we estimate τyw and α for each study catchment: the resulting estimates are reported in Fig. 6b and are listed in Table 2. Our estimates of the α parameter (0.19≤α≤2.1) are consistent with the shape factor range (0.2≤α≤2) investigated by Kirchner (2016a).
Consequently, the τyw obtained by applying Eq. (7) falls between 1.38 and 3.16 months. As expected, τywvaries in a narrow range consistent with the explanation provided by Kirchner (2016a). However, in order to have a fully coherent metric for all the catchments, the optimal procedure would be to set τyw and to calculate the young water fraction corresponding to this τyw.
Nevertheless, establishing a constant τyw for all the catchments could be a tricky choice. Indeed, by setting a τyw higher than that obtained via Eq. (7), we are improperly using the TTD to estimate the young water fraction. From this reasoning, the only solution is to set τyw equal to the overall minimum τyw. This choice ensures that the TTD is used properly to estimate the young water fraction in all the sites. What we could expect is that all the young water fractions would be lower with respect to those obtained with the amplitude ratio approach. However, changes in young water fraction depending on changes in τyw are unintuitive, since they vary with the TTD shape. Accordingly, a constant reduction of τyw would change the
area under the transit time pdf differently based on the α value.
Thus, the overall effect of the reduced τyw upon the young
water fraction remains a priori unknown.
Since in this study we are using the amplitude ratio approach, our
Fyw* estimates refer to the proportion of runoff younger than an inconsistent threshold age. This variation by catchment (albeit limited) is the main limitation of this work.
The role of Quaternary deposits
In line with the results of Arnoux et al. (2021), we find a negative statistically significant correlation between Fyw* and WFI (ρSpearman=-0.5, p value < 0.01; see Fig. S6), suggesting (unsurprisingly) that more groundwater contribution to streamflow increases the water age. WFI and Fqd values for all the study catchments are reported in Table 2. To analyze the relationship of
Fyw* with Quaternary deposits, we exclude the SOU and BCC catchments, since they show Fqd=0 (see Table 2). The inclusion of catchments with Fqd=0 would bias the analysis, since an absent parameter cannot modulate the share of groundwater and thus the young water fraction in the stream.
By focusing first only on the snow-dominated catchments, a linear fit on the
data returns a negative slope of -0.36 (R2=0.52), indicating a reduction of Fyw* with increasing Fqd (Fig. 7a). Moreover, we find a Spearman rank correlation coefficient of -0.6 with a p value of 0.13, meaning a negative but not statistically significant correlation between Fyw* and Fqd. This result can be explained by considering several factors. First, water storage in Quaternary deposits is not the only groundwater storage contribution to the stream in such environments: additional storage is provided by the bedrock fractures (Gleeson et al., 2014; Jasechko et al., 2016; Martin et al., 2021), possibly caused by rock stress and high erosion rates and by the bedrock geology, which has influence on groundwater retention capacity (Hayashi, 2020). Second, the area covered by Quaternary deposits could be an insufficient proxy of the groundwater storage potential: the knowledge of the thickness of these deposits (i.e., their volume) and the bedrock topography are crucial factors for controlling groundwater storage (Arnoux et al., 2021; Hayashi, 2020), but corresponding data are not available to date.
(a) Young water fraction against fraction of Quaternary deposits. Points dimension is proportional to τyw. (b) Fraction of Quaternary deposits against mean catchments elevation.
Fyw* values of the hybrid catchments reveal a weak positive correlation with Quaternary deposits (ρSpearman=0.13, p value = 0.74), while for rainfall-dominated catchments they show a negative correlation (ρSpearman=-0.52, p value = 0.2); however, both correlations are not statistically significant. These weak correlations suggest that Fqd represents only a limited part of the catchment geology responsible for groundwater flow and that it can only be considered as a first-order measure of geological groundwater storage.
We furthermore observe that Fqd decreases with mean catchment elevation in our data set (Fig. 7b), revealing a negative statistically significant correlation (ρSpearman=-0.5, p value < 0.01). This negative correlation reflects the fact that Fqd decreases when the mean slope increases (Arnoux et al., 2021).
To conclude, we stress that more catchments and more geological information
would be required to statistically validate these observations about the role of the groundwater storage potential for explaining young water fraction variations.
The role of groundwater flow (baseflow) in <inline-formula><mml:math id="M520" display="inline"><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mi mathvariant="normal">yw</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>Baseflow under different hydroclimatic regimes
The baseflow time series resulting from the baseflow separation of Duncan (2019) for six representative study catchments (two of each regime) are reported in Fig. 8 (complete results in Fig. S3). This figure shows the effect of groundwater recharge from rain and snowmelt through the “smoothed” baseflow proposed by Duncan (2019). This “smoothing” simulates a delayed storage contribution to the stream following the recharge phase during an input event. This recharge phase promotes the system wetness, thus favoring increasing quick flow. The increasing quick flow during events also leads to an increase of Fyw(ti), as found previously (von Freyberg et al., 2018). However, the relative amount of baseflow remains high during events: the mean baseflow fraction during the high-flow period is 0.49 and 0.52 for hybrid and rainfall-dominated catchments, while it is 0.63 for snow-dominated catchments. In agreement with worldwide stable-isotope-based hydrograph separation results (Jasechko, 2019), this
outcome underlines the mobilization of stored water (i.e., old water) during
rainfall and snowmelt events, and this process seems to be particularly
relevant in high-elevation catchments.
Baseflow separation for six representative study catchments using
the Duncan (2019) filter. The black area represents the daily discharge, while the colored area represents the estimated daily baseflow. The darker color represents a time step in which at least 85 % of the daily discharge is composed by baseflow.
Accordingly, in snow-dominated systems, the snowmelt largely transits
through the groundwater store (Hayashi, 2020; Cochand et al., 2019; Du et al., 2019; Flerchinger et al., 1992; Martinec, 1975), as schematized in Fig. 4, and the very high baseflow in high mountain catchments during summer is a direct sign of meltwater infiltration and percolation to groundwater that pushes old snowmelt (the main groundwater storage component) out to the stream network, as also found by Martinec (1975). This is also supported by the fact that groundwater, in such catchments, often has the isotopic signature of snowmelt (Michelon et al., 2023; Pavlovskii et al., 2018).
When examining the overall flow (and not only at the high-flow periods),
Fbf is generally lower for hybrid catchments (mean of Fbf=0.67) than for rainfall-dominated (mean of Fbf=0.74) and snow-dominated catchments (mean of Fbf=0.83). The values of Fbf for all the study sites are reported in Table 2. In the BCC
catchment, the Fbf (0.87) is consistent with the previous findings of Penna et al. (2016) who found, using stable water isotopes, that on average between 80 % and 98 % of the discharge in BCC is composed of pre-event water (assumed to represent groundwater). On average, the Fbf computed over the entire PoS is higher than that computed during the high-flow periods. This result suggests, unsurprisingly, that the largest percentage of base flow is released during low-flow periods. Accordingly, the variations of Fbf with elevation among different catchments (Fig. 9b) can be explained considering the changes in low-flow duration (LFD) with elevation, as will be discussed in Sect. 4.4.
(a) Young water fraction plotted against fraction of baseflow: vertical and horizontal bars represent ± standard deviation. Gray area represents the 95 % prediction bounds of a linear regression of
Fyw* on Fbf. Points dimension is proportional to τyw. (b) Fraction of baseflow and young water fraction against mean elevation. Bars with black edge indicate Fbf (left axis), while bars with gray edge indicate Fyw*(right axis). Vertical bars represent ± standard deviation.
Baseflow filters were already applied in previous studies, and their results
were correlated with Fyw*. For example, von Freyberg et al. (2018) found a strong positive correlation (ρSpearman=0.73, p value < 0.001) between Fyw* and the quick-flow index (QFI), calculated as the mean ratio between (Q-Qbf) and Q, where Q is the daily discharge and Qbf is the daily baseflow calculated in their paper with the Lyne and Hollick (1979) baseflow filter. By relating the Fbf to Fyw*, we have found a strong negative correlation (ρSpearman=-0.73, p value < 0.001), as shown in Fig. 9a,
consistent with the results of von Freyberg et al. (2018).
In snow-free systems, Fyw* is by definition related to Fbf: baseflow is composed of groundwater and groundwater is the dominant source of old water in such systems (in absence of large lakes). In snow-influenced systems, through the “direct input” approach for estimating Fyw*, we consider the snowpack (i.e., a temporarily old water storage) as part of the catchment storage. However, the share of snowmelt (with age > 3 months) that flows off quickly as surface or fast subsurface runoff will not show up in Fbf. In other words, Fbf is not able to take into account all the snowmelt but only the part of meltwater that infiltrates and recharges the groundwater storage, which is a large portion of the overall snowmelt.
The complementarity between the fraction of baseflow (<inline-formula><mml:math id="M557" display="inline"><mml:mrow><mml:msub><mml:mi>F</mml:mi><mml:mi mathvariant="normal">bf</mml:mi></mml:msub></mml:mrow></mml:math></inline-formula>) and the young water fraction (<inline-formula><mml:math id="M558" display="inline"><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mi mathvariant="normal">yw</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>)
A by-product of this work is that the Fbf, estimated with the Duncan (2019) baseflow filter, is roughly the complementary term of Fyw* (Fig. 9a and b), which is an important result for catchments where isotope measurements are missing. In such catchments, Fyw* could potentially be estimated without the application of the amplitude ratio approach as follows:
Fyw*≃1-Fbf.
Some of our case studies show considerable “residuals” of 1-(Fbf+Fyw*) (Fig. 9b). This is partially due to the uncertainty of the parameters used for estimating Fbf. In this regard, Duncan (2019) suggests some calibration guidelines to obtain optimal parameters set for baseflow estimation per catchment. In this work, we did not use the calibration guidelines, but we simply used the recession parameter proposed by Nathan and McMahon (1990) in order to achieve factual and reproducible results. In addition, the estimation of baseflow during an event is generally less rigorous than during the recession phase (Duncan, 2019), affecting the Fbf estimation. Moreover, Fyw* values are influenced by the sampling rate: the higher the frequency of sampling is, the higher the young water fraction is (Gallart
et al., 2020a; Stockinger et al., 2016). Thus, the young water fraction
calculated with the amplitude ratio approach generally underestimates the
“theoretical” young water fraction, and we simply compensate by computing
the flow-weighted young water fraction (Fyw*). In hybrid and
snow-dominated catchments, these “residuals” can also be explained by
considering that the Fbf does not include surface runoff or fast
lateral subsurface flow of meltwater, likely older than the estimated threshold ages, following a snowmelt event. On the other hand, these residuals might also be related to the non-linear recession behavior of
catchments, which was shown by Santos et al. (2018) to dominate Swiss low-elevation (i.e., rain dominated) catchments, when the exponential recession assumption of the baseflow filter necessarily leads to less reliable results (Duncan, 2019).
Low-flow duration (LFD) and <inline-formula><mml:math id="M569" display="inline"><mml:mrow><mml:msubsup><mml:mi>F</mml:mi><mml:mi mathvariant="normal">yw</mml:mi><mml:mo>*</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula>
The values of LFD for all the study sites are reported in Table 2. Specifically, LFD is lower for hybrid catchments (median of LFD = 0.39), and it is increasingly higher for rainfall (median of LFD = 0.50) and snow-dominated catchments (median of LFD = 0.62). In hybrid catchments, the presence of rain and snowmelt events spanning large parts of the year and the relatively low evapotranspiration (compared to rainfall-dominated catchments) due to reduced temperatures (Goulden et al., 2012) dramatically reduces the duration of low-flow periods, and this is also visible from the recurring discharge peaks (Fig. 8). In low-lying, rain-dominated catchments, evapotranspiration and precipitation are respectively higher and lower than in hybrid catchments, leading to longer low-flow periods (usually during summer and autumn). Under current climate and according to our data set, in snow-dominated catchments, we observe longer winter low-flow periods (streamflow decreasing below 0.5 to 1 mm d-1 for the highest locations; see Fig. S7) on an annual scale than in hybrid catchments. To gain additional insights into the high LFD in snow-dominated catchments and the low LFD in hybrid catchments, it is necessary to further consider the role of snowpack persistence, discussed in the following section. The variations of LFD with elevation are shown in Fig. 11b.
Low-flow periods are typically baseflow dominated (or old water dominated).
Accordingly, as anticipated in Sect. 4.3, the variation of Fbf between catchments reflects the proportion of the low-flow duration during the PoS. We observe that the higher the LFD is, the higher the Fbf is: in fact, they are strongly positively correlated
(ρSpearman=0.97, p value < 0.001) as shown in Fig. 10. The negative correlation between LFD and Fyw* is lower (ρSpearman=-0.74, p value < 0.001; Fig. 11a) but nevertheless suggests that LFD is an important predictor for Fyw*.
Fbf against the low-flow duration (LFD).
(a)Fyw* against the low-flow duration, LFD. Boxplots of Fyw* for catchments belonging to the same regime are plotted in correspondence to the median LFD. Points dimension is proportional to τyw. (b) LFD against mean elevation.
Time series of fSCA for six representative study catchments (two for each hydroclimatic regime), illustrating the gradual increase of the FSCA passing from rainfall-dominated to snow-dominated catchments.
The role of snowpack persistence
We explore next the presence of an ephemeral or seasonal snowpack as a
relevant factor for the time concentration of liquid water input and for
LFD. We consider the FSCA, calculated as reported in Sect. 3.2, as a proxy of the snowpack persistence. The FSCA values for all the study sites are reported in Table 2. The underlying fSCA time series for six representative study catchments (two for each hydroclimatic regime) are reported in Fig. 12 (for complete results, see Fig. S4). All the catchments characterized by a seasonal snow cover (i.e.,
snow dominated) reveal a high FSCA (>0.40, median of
FSCA=0.51). Gradually smaller FSCA values correspond to increasingly more ephemeral snowpacks with intermittent snowmelt events
during the winter season (Petersky and Harpold, 2018), as reflected by the spiky fSCA time series of hybrid and rainfall-dominated catchments (Fig. 12).
(a) Young water fraction against FSCA. The gray area represents the perceptual bell-shaped behavior of Fyw* with increasing FSCA. Points dimension is proportional to τyw. (b)FSCA against mean elevation.
Our results exhibit a bell-shaped behavior of Fyw* with varying FSCA (Fig. 13a). Specifically, we observe a general increase of Fyw* for FSCA values roughly below 0.3. This result can be explained considering that especially in hybrid catchments (median of FSCA=0.28), but partially also in rain-dominated catchments (median of FSCA=0.13), streamflow receives relatively more young water from ephemeral snowpacks. These short-lived snowpacks melt during the winter season resulting in only a short delay between precipitation input and melt (i.e., no water aging in the snowpack), and correspondingly meltwater flows off quickly into the stream (reducing LFD, Fig. 14), e.g., in presence of a frozen surface soil layer. In fact, ephemeral and slightly thick snowpacks do not protect the underlying soil from freezing (Harrison et al., 2021; Rey et al., 2021). Even for low-elevation locations (<≈1500 m a.s.l.), freezing conditions are regularly observed during winter (Keller et al., 2017).
Low-flow duration (LFD) against FSCA.
For FSCA values roughly higher than 0.3, we observe a decrease of
Fyw* with FSCA; here all the catchments of our data set are snow dominated. The mechanisms at play here are as follows: (i) in catchments with seasonal snowpacks, streamflow receives snowmelt in spring and summer that is at least partly older than 2–3 months (because part of the snow fell more than 3 months before the melt occurs); and (ii) the building up of a persistent, deep snowpack can promote deep vertical infiltration during the main melt period, either by insulating the soil and thereby preventing/reducing freezing (Harrison et al., 2021; Rey et al., 2021; Jasechko et al., 2016) or by gradual soil thawing during the melt period (Rey et al., 2021; Scherler et al., 2010). The temporal dynamic of snow accumulation and melt supports the pivotal role of snowmelt in recharging groundwater during summer in high-elevation environments (Cochand et al., 2019; Du et al., 2019; Flerchinger et al., 1992). A similar result was also
found for dolomitic catchments (such as BCC and OVA) by Lucianetti et al. (2020), who discovered that different proportions of rain and snow contribute to the recharge of springs in the Dolomites, with a gradually higher meltwater contribution in springs with increasing elevation. This role of snowmelt supports our analytical choice of computing Fyw* through the “direct input” approach, thus considering the snowpack as part of the catchment storage. In addition, the potentially large shares of meltwater that recharge groundwater via deep vertical infiltration also result in old water sustaining winter baseflow (Fig. S5): the persistent snowpack and the absence of a liquid water input favor a groundwater storage release that creates a longer winter low-flow period that increases LFD (Fig. 14), thus reducing Fyw*, as discussed in Sect. 4.4.
FSCA is strongly correlated with the mean catchment elevation in our data set (ρSpearman=0.97, p value < 0.01; Fig. 13b). A posteriori, we could have considered mean elevation instead of FSCA as a proxy for snowpack persistence. However, a priori, it could be approximative to describe the snow cover persistence only with the increasing elevation: the persistence of snow in a catchment also depends on overall topographic and climatic characteristics, specifically relating to snow and aspect (Painter et al., 2023). In fact, catchments with very different characteristics (e.g., different elevation ranges and different areas) can reveal a similar mean elevation, but the snowpack persistence could considerably change. This is the reason why we focused on FSCA that integrates these physical factors.
Perceptual model of the hydrological processes that drive the
young water fraction variations with elevation. This model emerges from our
analysis and harmonizes these results with those of previous studies. For
snow-dominated and hybrid catchments, we indicate the dominant processes,
occurring during summer and during winter, that lead to low and high
Fyw*, respectively.
The above mechanisms are unable to explain the hydrological function of the
glacier-dominated NBPV catchment, which has a very high Fyw* and is an outlier among the snow-dominated catchments (Fig. 13a). The high Fyw* of the high-elevation glacier-covered (42 %) catchment can be explained considering that the glacier melt produces high amounts of streamflow that transit the glacier system very quickly during the summer, given generally fast englacial and subglacial flow paths and the often limited water storage capacity in the glacier forefield (Müller et al., 2022; Saberi et al., 2019; Jansson et al., 2003). Schmieder et al. (2019) also found a high young water fraction in an Austrian glacier-covered (35 %) catchment, leading them to the conclusion that the basin behaves locally like a “Teflon basin” with quickly transmitted ice melt.
Process interplay along elevation: perceptual model
The identified key drivers of young water fractions for rainfall-dominated,
hybrid and snow-dominated catchments can conveniently be summarized into a
perceptual model of the involved hydrological processes and their seasonal
interplay (Fig. 15).
High-elevation catchments are characterized by long winter low-flow periods,
resulting from the build-up of a seasonal snowpack, and are sustained by the
emptying of groundwater (or old water) storage. Accordingly, such storage
releases stored water, mainly old meltwater, for prolonged periods where the
snowpack can last for several months (typically from December to early April) before releasing water during the melting season. Such seasonal snowpack can protect the underlying soils from freezing, thus promoting meltwater infiltration and groundwater recharge. From this viewpoint, snowpack is considered as part of the catchment storage, and there is a thin line between groundwater and meltwater in snow-dominated catchments. Snowmelt or rain events push out old meltwater to the stream during summer, suggested by the relatively high amount of daily baseflow during the melting season. During this period, the high catchment wetness might even lead to saturation and thereby favor fast flow paths of meltwater or rainwater, which in turn can temporarily increase the young water fraction. Despite this increase during high-flow periods, the prevailing winter low-flow periods in such systems lead to a reduction of the average annual young water fraction.
In catchments with an ephemeral snowpack, at lower elevations, snowmelt events occur regularly during winter such that water released from the
corresponding short-lived snowpack is likely younger than 3 months. Moreover, ephemeral snowpacks do not protect the underlying soils from freezing, and rapid flow paths can emerge during episodic or long-term soil surface freezing, by increasing the young water fraction. The high Fyw* of such systems is also explained by the simultaneity of
snowmelt and rain events during extended parts of the year (leading to large
volume of annual precipitation) and the relatively low (compared to
rainfall-dominated catchments) evapotranspiration. Both processes increase
the catchment's wetness and reduce the low-flow period's length.
Finally, at the lowest elevations, lower amounts of precipitation and higher
evapotranspiration favor longer low-flow periods, mainly sustained by old
groundwater from alluvial aquifers, which lead to both a Fyw* and a catchment wetness reduction. Further, the relatively flat topographies at the lowest elevations favor slow flow paths increasing the transit times of water.
How well current hydrological models can represent the interplay of these
processes along elevation gradients is left for future research, but our
perceptual model builds a solid basis for an improvement of theory-driven
models (Clark et al., 2016).
Conclusion
This study proposes a conceptualization of the processes behind changes in
young water fraction (Fyw*) with elevation, defined here following Kirchner (2016a) as stream water that is younger than a threshold age of about 2–3 months. The analysis is focused on amplitude-ratio-based young water fractions for a set of 27 study catchments located in Switzerland and Italy, which span a wide range of geological and hydroclimatic conditions. The young water fraction estimates (Fyw*) obtained from the phase and amplitude information of seasonal isotope cycles correspond to different young water threshold ages for the different catchments, which is a limitation of this work. However, the threshold ages vary only modestly from about 1.5 to 3 months.
Our analysis focuses on mountainous catchments to fill the knowledge gap,
referring to the surprisingly low young water fractions at high elevations
(>1500 m a.s.l.), but we have also considered catchments at lower elevations to obtain a complete picture of the dominant hydrological processes at different elevations.
We have focused on variables and processes that were not previously considered for explaining elevation gradients of young water fraction. We
have investigated the role of (i) groundwater storage potential, (ii) catchment storage contribution to the stream, (iii) low-flow duration and
(iv) snowpack persistence. Our results suggest that (ii), (iii) and (iv) are interconnected: low-flow periods are generally sustained by old water deriving from the catchment storage, and the length of such periods is driven by the snowpack persistence at high elevations. The proportion of low-flow periods during the period of isotope sampling strongly influences the amount of old water contributing to the stream, thus reducing the estimated Fyw*. Consequently, the low-flow duration, which varies with elevation, can be retained as a driver of the Fyw* changes with elevation. Given the importance of low-flow periods, we have also investigated the role of groundwater storage potential, represented here by the portion of catchment area covered by Quaternary deposits. Our results suggest that an exhaustive description of the groundwater storage potential should be completed with more detailed geological information, e.g., the geology and topography of bedrock, the fraction of fractured bedrock and the deposits' thickness, which is challenging to retrieve from a geological data set. We have brought together the results of this analysis in a perceptual model that describes a framework for how hydrological processes control the Fyw* according to elevation, laying the foundations for an improvement of theory-driven models.
The strong complementarity between Fyw* and the mean fraction of baseflow obtained for our data set suggests that Fyw* could be estimated starting from automated baseflow separation techniques for catchments in which stable water isotope measurements are not available. This complementarity should however be validated in future work, by considering, e.g., alternative baseflow separation techniques and different
hydroclimatic conditions.
Finally, the conceptualization of the hydrological processes described in
this paper do not fit the high young water fraction of the single glacier-dominated catchment of our data set. In conclusion, we encourage
future studies to compare and to collect isotopic data from glacier-dominated catchments to better understand the processes in such systems that, under glacier retreat due to climate change, will see a gradual transition to purely snow-dominated systems.
Data availability
Time series of both δ2H and δ18O in streamflow and precipitation, complemented with MeteoSwiss daily precipitation data, for the 22 Swiss catchments investigated by von Freyberg et al. (2018), are available in the data repository Zenodo at 10.5281/zenodo.4057967 (Staudinger et al., 2020). Meteorological, hydrological and isotope data of VdN, BCC and NBPV catchments are available at
https://onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1002/hyp.13937&file=hyp13937-sup-0009-Supinfo2.zip (Ceperley et al., 2020).
Daily discharge data for the ERL, LÜM and VOG catchments are provided by
the Swiss Federal Institute for Forest, Snow and Landscape research (WSL), Birmensdorf, Switzerland. Streamflow data for the AAB and GUE
catchments are provided by the Office for Waste, Water, Energy and Air (WWEA) of the Canton of Zurich and by the Office for Water and Waste of the Canton of Bern, respectively. Daily discharge data of the remaining 17 Swiss catchments studied by von Freyberg et al. (2018) are provided by the Swiss Federal Office for the Environment (FOEN).
The .shp of the AAB, GUE, ERL, LÜM and VOG catchment boundaries are
available from the data repository Zenodo at
10.5281/zenodo.4057967 (Staudinger et al., 2020). The .shp of NBPV, BCC and VdN catchments are provided by Giulia Zuecco and Anthony Michelon (University of Lausanne, Switzerland) as personal communication. The DOR and SOU catchment boundaries are delineated in a GIS environment using the 10 m resolution digital elevation model (DEM) available from the Aosta Valley Regional Geoportal. Finally, the catchment boundaries of the remaining 17 Swiss catchments investigated by von Freyberg et al. (2018) are directly obtained from the Swiss Federal Office for the Environment (FOEN).
Quaternary cover for all Swiss catchments has been calculated using the
Geological Atlas of Switzerland (GeoCover data set; 1:25000 scale) available from the Federal Office of Topography swisstopo. For the DOR and SOU catchments, the vectorized Valsavaranche geological map (1:100000 scale) is provided by the Cartography Office of SCT Geoportal. For the NBPV and BCC catchments, the .shp of unconsolidated sediments is provided by Giulia Zuecco.
DOR and SOU data are available from Alessio Gentile upon reasonable request.
The supplement related to this article is available online at: https://doi.org/10.5194/hess-27-2301-2023-supplement.
Author contributions
AG, NC, BS and SF identified the research gap, defined the methodology, developed the perceptual model and prepared the paper. AG analyzed the data set. DG, DC, MP and SF collected the water samples for the DOR and SOU catchments. GZ analyzed spatial data related to NBPV and BCC catchments. All authors revised the paper and gave final approval to the submitted version.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
This publication is part of the project NODES, which has received funding from the MUR-M4C2 1.5 of PNRR with grant agreement no. ECS00000036. We warmly thank the COST Action CA19120 – WATSON (WATer isotopeS in the critical zONe) for the acceptance of the application procedure for one virtual mobility (VM) and one short-term scientific mission (STSM). Both activities allowed to speed up the planning and conceptualization of this work as well as to stimulate the collaboration, the sharing of data and ideas. We acknowledge the support of the Valsavarenche Municipality and the Gran Paradiso National Park. The authors thank Chiara Marchina (University of Padova, Italy) for the isotopic analyses of DOR and SOU samples. Finally, we thank Jana von Freyberg and the anonymous referee for their comments that helped to improve the paper significantly.
Financial support
This research has been supported by the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (grant no. PP00P2_157611), by the PRIN MIUR 2017SL7ABC_005 WATZON Project and by the MIUR – Dipartimento di Eccellenza DIST department funds.
Review statement
This paper was edited by Rohini Kumar and reviewed by Jana von Freyberg and one anonymous referee.
ReferencesAalstad, K., Westermann, S., and Bertino, L.: Evaluating satellite retrieved
fractional snow-covered area at a high-Arctic site using terrestrial
photography, Remote Sens. Environ., 239, 111618, 10.1016/j.rse.2019.111618, 2020.Andermann, C., Longuevergne, L., Bonnet, S., Crave, A., Davy, P., and Gloaguen, R.: Impact of transient groundwater storage on the discharge of
Himalayan rivers, Nat. Geosci., 5, 127–132, 10.1038/ngeo1356, 2012.Arnoux, M., Brunner, P., Schaefli, B., Mott, R., Cochand, F., and Hunkeler,
D.: Low-flow behavior of alpine catchments with varying quaternary cover
under current and future climatic conditions, J. Hydrol., 592, 125591, 10.1016/j.jhydrol.2020.125591, 2021.Baraer, M., McKenzie, J. M., Mark, B. G., Bury, J., and Knox, S.:
Characterizing contributions of glacier melt and groundwater during the dry
season in a poorly gauged catchment of the Cordillera Blanca (Peru), in:
Advances in Geosciences, 4th EGU Alexander von Humboldt Conference “The
Andes: Challenge for Geosciences” – 4th Alexander von Humboldt
International Conference on The Andes: Challenge for Geosciences, 28 November 2008, Santiago de Chile, Chile, 41–49, 10.5194/adgeo-22-41-2009, 2009.Baraer, M., McKenzie, J., Mark, B. G., Gordon, R., Bury, J., Condom, T., Gomez, J., Knox, S., and Fortner, S. K.: Contribution of groundwater to the
outflow from ungauged glacierized catchments: a multi-site study in the tropical Cordillera Blanca, Peru, Hydrol. Process., 29, 2561–2581,
10.1002/hyp.10386, 2015.Benettin, P., Bailey, S. W., Rinaldo, A., Likens, G. E., McGuire, K. J., and
Botter, G.: Young runoff fractions control streamwater age and solute
concentration dynamics, Hydrol. Process., 31, 2982–2986, 10.1002/hyp.11243, 2017.Carroll, R. W. H., Bearup, L. A., Brown, W., Dong, W., Bill, M., and Willlams, K. H.: Factors controlling seasonal groundwater and solute flux from snow-dominated basins, Hydrol. Process., 32, 2187–2202,
10.1002/hyp.13151, 2018.Carturan, L.: Replacing monitored glaciers undergoing extinction: a new
measurement series on La Mare Glacier (Ortles-Cevedale, Italy), J. Glaciol., 62, 1093–1103, 10.1017/jog.2016.107, 2016.Ceperley, N., Zuecco, G., Beria, H., Carturan, L., Michelon, A., Penna, D.,
Larsen, J., and Schaefli, B.: Seasonal snow cover decreases young water
fractions in high Alpine catchments, Hydrol. Process., 34, 4794–4813,
10.1002/hyp.13937, https://onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1002/hyp.13937&file=hyp13937-sup-0009-Supinfo2.zip
(last acces: 8 June 2023), 2020.Chen, Z., Hartmann, A., Wagener, T., and Goldscheider, N.: Dynamics of water
fluxes and storages in an Alpine karst catchment under current and potential
future climate conditions, Hydrol. Earth Syst. Sci., 22, 3807–3823, 10.5194/hess-22-3807-2018, 2018.Christensen, C. W., Hayashi, M., and Bentley, L. R.: Hydrogeological
characterization of an alpine aquifer system in the Canadian Rocky Mountains, Hydrogeol. J., 28, 1871–1890, 10.1007/s10040-020-02153-7, 2020.Clark, M. P., Schaefli, B., Schymanski, S. J., Samaniego, L., Luce, C. H.,
Jackson, B. M., Freer, J. E., Arnold, J. R., Moore, R. D., Istanbulluoglu,
E., and Ceola, S.: Improving the theoretical underpinnings of process-based
hydrologic models, Water Resour. Res., 52, 2350–2365, 10.1002/2015WR017910, 2016.Clow, D. W., Schrott, L., Webb, R., Campbell, D. H., Torizzo, A., and
Dornblaser, M.: Ground Water Occurrence and Contributions to Streamflow in
an Alpine Catchment, Colorado Front Range, Groundwater, 41, 937–950,
10.1111/j.1745-6584.2003.tb02436.x, 2003.Cochand, M., Christe, P., Ornstein, P., and Hunkeler, D.: Groundwater Storage in High Alpine Catchments and Its Contribution to Streamflow, Water Resour. Res., 55, 2613–2630, 10.1029/2018WR022989, 2019.Cowie, R. M., Knowles, J. F., Dailey, K. R., Williams, M. W., Mills, T. J.,
and Molotch, N. P.: Sources of streamflow along a headwater catchment elevational gradient, J. Hydrol., 549, 163–178,
10.1016/j.jhydrol.2017.03.044, 2017.Di Marco, N., Righetti, M., Avesani, D., Zaramella, M., Notarnicola, C., and
Borga, M.: Comparison of MODIS and Model-Derived Snow-Covered Areas: Impact
of Land Use and Solar Illumination Conditions, Geosciences, 10, 134,
10.3390/geosciences10040134, 2020.Dozier, J.: Spectral signature of alpine snow cover from the landsat thematic mapper, Remote Sens. Environ., 28, 9–22, 10.1016/0034-4257(89)90101-6, 1989.Du, X., Fang, M., Lv, H., Cheng, T., Hong, P., and Liu, C.: Effect of snowmelt infiltration on groundwater recharge in a seasonal soil frost area:
a case study in Northeast China, Environ. Monit. Assess., 191, 151, 10.1007/s10661-019-7285-7, 2019.Duncan, H. P.: Baseflow separation – A practical approach, J. Hydrol., 575, 308–313, 10.1016/j.jhydrol.2019.05.040, 2019.Engel, M., Penna, D., Bertoldi, G., Dell'Agnese, A., Soulsby, C., and Comiti, F.: Identifying run-off contributions during melt-induced run-off events in a glacierized alpine catchment, Hydrol. Process., 30, 343–364, 10.1002/hyp.10577, 2016.Farr, T. G., Rosen, P. A., Caro, E., Crippen, R., Duren, R., Hensley, S.,
Kobrick, M., Paller, M., Rodriguez, E., Roth, L., Seal, D., Shaffer, S.,
Shimada, J., Umland, J., Werner, M., Oskin, M., Burbank, D., and Alsdorf, D.: The Shuttle Radar Topography Mission, Rev. Geophys., 45, RG2004, 10.1029/2005RG000183, 2007.Flerchinger, G. N., Cooley, K. R., and Ralston, D. R.: Groundwater response
to snowmelt in a mountainous watershed, J. Hydrol., 133, 293–311,
10.1016/0022-1694(92)90260-3, 1992.Frisbee, M. D., Phillips, F. M., Campbell, A. R., Liu, F., and Sanchez, S.
A.: Streamflow generation in a large, alpine watershed in the southern Rocky
Mountains of Colorado: Is streamflow generation simply the aggregation of
hillslope runoff responses?, Water Resour. Res., 47, W06512, 10.1029/2010WR009391, 2011.Gallart, F., Valiente, M., Llorens, P., Cayuela, C., Sprenger, M., and Latron, J.: Investigating young water fractions in a small Mediterranean mountain catchment: Both precipitation forcing and sampling frequency matter, Hydrol. Process., 34, 3618–3634, 10.1002/hyp.13806, 2020a.Gallart, F., von Freyberg, J., Valiente, M., Kirchner, J. W., Llorens, P.,
and Latron, J.: Technical note: An improved discharge sensitivity metric for
young water fractions, Hydrol. Earth Syst. Sci., 24, 1101–1107,
10.5194/hess-24-1101-2020, 2020b.Gascoin, S., Grizonnet, M., Bouchet, M., Salgues, G., and Hagolle, O.: Theia
Snow collection: high-resolution operational snow cover maps from Sentinel-2
and Landsat-8 data, Earth Syst. Sci. Data, 11, 493–514,
10.5194/essd-11-493-2019, 2019.Gentile, A.: GEE code that accompanies the paper: “Towards a conceptualization of the hydrological processes behind changesof young water fraction with elevation: a focus on mountainous alpine catchments”, Google Earth Engine, https://code.earthengine.google.com/8239cfe7aab498180e5c42475023cb80?noload=true
(last access: 8 June 2023), 2023.Gisolo, D., Previati, M., Bevilacqua, I., Canone, D., Boetti, M., Dematteis,
N., Balocco, J., Ferrari, S., Gentile, A., Nsassila, M., Heery, B., Vereecken, H., and Ferraris, S.: A Calibration Free Radiation Driven Model
for Estimating Actual Evapotranspiration of Mountain Grasslands (CLIME-MG),
J. Hydrol., 610, 127948, 10.1016/j.jhydrol.2022.127948, 2022.Gleeson, T. and Manning, A. H.: Regional groundwater flow in mountainous
terrain: Three-dimensional simulations of topographic and hydrogeologic
controls, Water Resour. Res., 44, W10403, 10.1029/2008WR006848, 2008.Gleeson, T., Moosdorf, N., Hartmann, J., and van Beek, L. P. H.: A glimpse
beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability
and porosity, Geophys. Res. Lett., 41, 3891–3898, 10.1002/2014GL059856, 2014.Gordon, R. P., Lautz, L. K., McKenzie, J. M., Mark, B. G., Chavez, D., and
Baraer, M.: Sources and pathways of stream generation in tropical proglacial
valleys of the Cordillera Blanca, Peru, J. Hydrol., 522, 628–644,
10.1016/j.jhydrol.2015.01.013, 2015.Gorelick, N., Hancher, M., Dixon, M., Ilyushchenko, S., Thau, D., and Moore,
R.: Google Earth Engine: Planetary-scale geospatial analysis for everyone,
Remote Sens. Environ., 202, 18–27, 10.1016/j.rse.2017.06.031, 2017.Goulden, M. L., Anderson, R. G., Bales, R. C., Kelly, A. E., Meadows, M.,
and Winston, G. C.: Evapotranspiration along an elevation gradient in California's Sierra Nevada, J. Geophys. Res., 117, G03028, 10.1029/2012JG002027, 2012.Harrington, J. S., Mozil, A., Hayashi, M., and Bentley, L. R.: Groundwater
flow and storage processes in an inactive rock glacier, Hydrol. Process., 32, 3070–3088, 10.1002/hyp.13248, 2018.Harrison, H. N., Hammond, J. C., Kampf, S., and Kiewiet, L.: On the
hydrological difference between catchments above and below the intermittent-persistent snow transition, Hydrol. Process., 35, e14411,
10.1002/hyp.14411, 2021.Hayashi, M.: Alpine Hydrogeology: The Critical Role of Groundwater in Sourcing the Headwaters of the World, Ground Water, 58, 498–510,
10.1111/gwat.12965, 2020.Hofmeister, F., Arias-Rodriguez, L. F., Premier, V., Marin, C., Notarnicola,
C., Disse, M., and Chiogna, G.: Intercomparison of Sentinel-2 and modelled
snow cover maps in a high-elevation Alpine catchment, J. Hydrol. X, 15, 100123, 10.1016/j.hydroa.2022.100123, 2022.Jansson, P., Hock, R., and Schneider, T.: The concept of glacier storage: a
review, J. Hydrol., 282, 116–129, 10.1016/S0022-1694(03)00258-0, 2003.Jasechko, S.: Global Isotope Hydrogeology – Review, Rev. Geophys., 57, 835–965, 10.1029/2018RG000627, 2019.Jasechko, S., Kirchner, J. W., Welker, J. M., and McDonnell, J. J.:
Substantial proportion of global streamflow less than three months old, Nat. Geosci., 9, 126–129, 10.1038/ngeo2636, 2016.Käser, D. and Hunkeler, D.: Contribution of alluvial groundwater to the
outflow of mountainous catchments, Water Resour. Res., 52, 680–697,
10.1002/2014WR016730, 2016.Keller, D. E., Fischer, A. M., Liniger, M. A., Appenzeller, C., and Knutti,
R.: Testing a weather generator for downscaling climate change projections
over Switzerland, Int. J. Climatol., 37, 928–942, 10.1002/joc.4750, 2017.Kirchner, J. W.: Aggregation in environmental systems – Part 1: Seasonal
tracer cycles quantify young water fractions, but not mean transit times, in
spatially heterogeneous catchments, Hydrol. Earth Syst. Sci., 20, 279–297, 10.5194/hess-20-279-2016, 2016a.Kirchner, J. W.: Aggregation in environmental systems – Part 2: Catchment mean transit times and young water fractions under hydrologic nonstationarity, Hydrol. Earth Syst. Sci., 20, 299–328,
10.5194/hess-20-299-2016, 2016b.Kirchner, J. W. and Knapp, J. L. A.: Technical note: Calculation scripts for
ensemble hydrograph separation, Hydrol. Earth Syst. Sci., 24, 5539–5558, 10.5194/hess-24-5539-2020, 2020.Li, L., Sullivan, P. L., Benettin, P., Cirpka, O. A., Bishop, K., Brantley, S. L., Knapp, J. L. A., van Meerveld, I., Rinaldo, A., Seibert, J., Wen, H.,
and Kirchner, J. W.: Toward catchment hydro-biogeochemical theories, WIREs
Water, 8, e1495, 10.1002/wat2.1495, 2021.Liu, F., Williams, M. W., and Caine, N.: Source waters and flow paths in an
alpine catchment, Colorado Front Range, United States, Water Resour. Res., 40, W09401, 10.1029/2004WR003076, 2004.Lucianetti, G., Penna, D., Mastrorillo, L., and Mazza, R.: The Role of Snowmelt on the Spatio-Temporal Variability of Spring Recharge in a Dolomitic Mountain Group, Italian Alps, Water, 12, 2256, 10.3390/w12082256, 2020.Lutz, S. R., Krieg, R., Müller, C., Zink, M., Knöller, K., Samaniego, L., and Merz, R.: Spatial Patterns of Water Age: Using Young Water Fractions to Improve the Characterization of Transit Times in Contrasting Catchments, Water Resour. Res., 54, 4767–4784, 10.1029/2017WR022216, 2018.Lyne, V. and Hollick, M.: Stochastic Time-Variable Rainfall-Runoff Modeling,
in: Institution of Engineers Australia National Conference, September 1979, Barton, Australia, Institute of Engineers Australia, 89–92,
https://www.researchgate.net/profile/Vincent-Lyne/publication/272491803_Stochastic_Time-Variable_Rainfall-Runoff_Modeling/links/54f45fb40cf299c8d9e6e6c1/Stochastic-Time-Variable-Rainfall-Runoff-Modeling.pdf
(last access: 8 June 2023), 1979.Martin, C., Kampf, S. K., Hammond, J. C., Wilson, C., and Anderson, S. P.:
Controls on Streamflow Densities in Semiarid Rocky Mountain Catchments, Water, 13, 521, 10.3390/w13040521, 2021.Martinec, J.: Subsurface flow from snowmelt traced by tritium, Water Resour. Res., 11, 496–498, 10.1029/WR011i003p00496, 1975.McDonnell, J. J.: Beyond the water balance, Nat. Geosci., 10, 396, 10.1038/ngeo2964, 2017.Michelon, A., Ceperley, N., Beria, H., Larsen, J., Vennemann, T., and Schaefli, B.: Hydrodynamics of a high Alpine catchment characterized by four
natural tracers, Hydrol. Earth Syst. Sci., 27, 1403–1430,
10.5194/hess-27-1403-2023, 2023.Müller, T., Lane, S. N., and Schaefli, B.: Towards a hydrogeomorphological understanding of proglacial catchments: an assessment of groundwater storage and release in an Alpine catchment, Hydrol. Earth Syst. Sci., 26, 6029–6054, 10.5194/hess-26-6029-2022, 2022.Nathan, R. J. and McMahon, T. A.: Evaluation of automated techniques for base flow and recession analyses, Water Resour. Res., 26, 1465–1473,
10.1029/WR026i007p01465, 1990.Painter, K. J., Gentile, A., and Ferraris, S.: A stochastic cellular automaton model to describe the evolution of the snow-covered area across a
high-elevation mountain catchment, Sci.e Total Environ., 857, 159195, 10.1016/j.scitotenv.2022.159195, 2023.Pavlovskii, I., Hayashi, M., and Lennon, M. R.: Transformation of snow
isotopic signature along groundwater recharge pathways in the Canadian Prairies, J. Hydrol., 563, 1147–1160, 10.1016/j.jhydrol.2017.09.053, 2018.Penna, D., van Meerveld, H. J., Zuecco, G., Dalla Fontana, G., and Borga, M.: Hydrological response of an Alpine catchment to rainfall and snowmelt events, J. Hydrol., 537, 382–397, 10.1016/j.jhydrol.2016.03.040, 2016.Petersky, R. and Harpold, A.: Now you see it, now you don't: a case study of
ephemeral snowpacks and soil moisture response in the Great Basin, USA,
Hydrol. Earth Syst. Sci., 22, 4891–4906, 10.5194/hess-22-4891-2018, 2018.Rey, D. M., Hinckley, E.-L. S., Walvoord, M. A., and Singha, K.: Integrating
observations and models to determine the effect of seasonally frozen ground
on hydrologic partitioning in alpine hillslopes in the Colorado Rocky Mountains, USA, Hydrol. Process., 35, e14374, 10.1002/hyp.14374, 2021.Saberi, L., McLaughlin, R. T., Ng, G.-H. C., La Frenierre, J., Wickert, A.
D., Baraer, M., Zhi, W., Li, L., and Mark, B. G.: Multi-scale temporal
variability in meltwater contributions in a tropical glacierized watershed,
Hydrol. Earth Syst.Sci., 23, 405–425, 10.5194/hess-23-405-2019, 2019.Santos, A. C., Portela, M. M., Rinaldo, A., and Schaefli, B.: Analytical
flow duration curves for summer streamflow in Switzerland, Hydrol. Earth Syst. Sci., 22, 2377–2389, 10.5194/hess-22-2377-2018, 2018.Scherler, M., Hauck, C., Hoelzle, M., Stähli, M., and Völksch, I.:
Meltwater infiltration into the frozen active layer at an alpine permafrost
site, Permafrost Periglac. Process., 21, 325–334, 10.1002/ppp.694, 2010.Schmieder, J., Seeger, S., Weiler, M., and Strasser, U.: `Teflon Basin' or
Not? A High-Elevation Catchment Transit Time Modeling Approach, Hydrology,
6, 92, 10.3390/hydrology6040092, 2019.Somers, L. D. and McKenzie, J. M.: A review of groundwater in high mountain
environments, WIREs Water, 7, e1475, 10.1002/wat2.1475, 2020.Somers, L. D., McKenzie, J. M., Mark, B. G., Lagos, P., Ng, G.-H. C., Wickert, A. D., Yarleque, C., Baraër, M., and Silva, Y.: Groundwater
Buffers Decreasing Glacier Melt in an Andean Watershed – But Not Forever,
Geophys. Res. Lett., 46, 13016–13026, 10.1029/2019GL084730, 2019.Staudinger, M., Stoelzle, M., Seeger, S., Seibert, J., Weiler, M., and Stahl, K.: Catchment water storage variation with elevation, Hydrol. Process., 31, 2000–2015, 10.1002/hyp.11158, 2017.Staudinger, M., Seeger, S., Herbstritt, B., Stoelzle, M., Seibert, J., Stahl, K., and Weiler, M.: The CH-IRP data set: fortnightly data of δ2H and δ18O in streamflow and precipitation in Switzerland (Version 2), Zenodo [data set], 10.5281/zenodo.4057967, 2020.Stockinger, M., Reemt Bogena, H., Lücke, A., Stumpp, C., and Vereecken,
H.: Time variability and uncertainty in the fraction of young water in a
small headwater catchment, Hydrol. Earth Syst. Sci., 23, 4333—4347, 10.5194/hess-23-4333-2019, 2019.Stockinger, M. P., Bogena, H. R., Lücke, A., Diekkrüger, B.,
Cornelissen, T., and Vereecken, H.: Tracer sampling frequency influences
estimates of young water fraction and streamwater transit time distribution,
J. Hydrol., 541, 952–964, 10.1016/j.jhydrol.2016.08.007, 2016.Stoelzle, M., Schuetz, T., Weiler, M., Stahl, K., and Tallaksen, L. M.: Beyond binary baseflow separation: a delayed-flow index for multiple streamflow contributions, Hydrol. Earth Syst. Sci., 24, 849–867,
10.5194/hess-24-849-2020, 2020.von Freyberg, J., Allen, S. T., Seeger, S., Weiler, M., and Kirchner, J. W.:
Sensitivity of young water fractions to hydro-climatic forcing and landscape
properties across 22 Swiss catchments, Hydrol. Earth Syst. Sci., 22, 3841–3861, 10.5194/hess-22-3841-2018, 2018.
Weingartner, R. and Aschwanden, H.: Abflussregimes als Grundlage zur
Abschätzung von Mittelwerten des Abflusses, Hydrologischer atlas der
Schweiz, Tafel 5.2, https://hydrologischeratlas.ch/produkte/druckausgabe/fliessgewasser-und-seen/tafel-5-2
(last access: 8 June 2023), 1992.Williams, M. W., Wilson, A., Tshering, D., Thapa, P., and Kayastha, R. B.:
Using geochemical and isotopic chemistry to evaluate glacier melt contributions to the Chamkar Chhu (river), Bhutan, Ann. Glaciol., 57, 339–348, 10.3189/2016AoG71A068, 2016.Wilson, A. M., Williams, M. W., Kayastha, R. B., and Racoviteanu, A.: Use of
a hydrologic mixing model to examine the roles of meltwater, precipitation
and groundwater in the Langtang River basin, Nepal, Ann. Glaciol., 57, 155–168, 10.3189/2016AoG71A067, 2016.Wilusz, D. C., Harman, C. J., and Ball, W. P.: Sensitivity of Catchment Transit Times to Rainfall Variability Under Present and Future Climates,
Water Resour. Res., 53, 10231–10256, 10.1002/2017WR020894, 2017.Xia, C., Zuecco, G., Chen, K., Liu, L., Zhang, Z., and Luo, J.: The
estimation of young water fraction based on isotopic signals: challenges and
recommendations, Front. Ecol. Evol., 11, 1114259, 10.3389/fevo.2023.1114259, 2023.Zuecco, G., Carturan, L., De Blasi, F., Seppi, R., Zanoner, T., Penna, D.,
Borga, M., Carton, A., and Dalla Fontana, G.: Understanding hydrological
processes in glacierized catchments: Evidence and implications of highly
variable isotopic and electrical conductivity data, Hydrol. Process., 33, 816–832, 10.1002/hyp.13366, 2019.