# I. Introduction In recent years, there has been a noticeable increase in the severity of droughts in central Europe and various other regions around the world (Balting et al. (2021), Satoh et al. (2022)). These prolonged periods of water scarcity have had significant impacts on the natural water cycle and the ecosystems that rely on it (Trumbore et al. (2015), Neumann et al. (2017), Senf et al. (2020)). As a result, addressing this issue requires a comprehensive reevaluation of our approach to forest and landscape design (Liu et al. (2022), Gvein et al. (2023)), along with the implementation of appropriate hydrological methods. Isotope-based hydrograph separations can be particularly helpful in this regard, as they can provide insights into the travel times of water within the catchment. Over the past few decades, the separation of storm hydrographs using stable isotope tracers has become a standard method for investigating runoff generation processes in catchment hydrology. The early pioneering work was accomplished in the late 1960s and 1970s (Pinder and Jones (1969), Dincer et al. (1970), Martinec et al. (1974), and Fritz et al. (1976), Sklash and Farvolden (1979)), and over the years, the methodology has been progressively expanded and adapted to address the challenges and tasks found in the feld (Klaus and McDonnell (2013), Jasechko (2019)). Hoeg (2019) recently proposed a method that iteratively extends the standard two-component separation, such that n time components are separated by using a single stable isotope tracer. This approach can be used to trace the event water over a much longer period after the initial event, hence expanding the space of addressable use cases for catchment hydrologists. Hoeg (2019) applied the new method to an experimental data set of the mountainous Zastler catchment (18.4 km 2 , Southern Black Forest, Germany) and compared the outcome with previous investigations in that area, showing the influence of antecedent moisture conditions on the event water contributions of subsequent events. This iterative extension uncovered the temporal structure of the pre-event component, and enabled a closer look at the temporal composition of the pre-event water, hence determining the extent to which recent events were involved. A catchment response pattern related to antecedent moisture conditions similar to that of the Zastler catchment was found by Iorgulescu et al. (2007); they used a hydrochemical model based on a parameterization of three runoff components (direct precipitation, acid soil water, and deep groundwater) to predict conservative tracer data in the Haute-Mentue catchment (12.5 km 2 , Swiss Plateau). The authors concluded that the soil water component that corresponds to recent water stored in the upper soil horizons dominates catchment outflow in wet conditions but is virtually absent in dry conditions. James and Roulet (2009) formulated antecedent moisture conditions and catchment morphology as controls on the spatial patterns of runoff generation. Based on stable isotopes, they examined the spatial patterns of storm runoff generation from eight small nested forest catchments ranging in size from 0.07 to 1.5 km 2 (formerly the glaciated terrain of Mont Saint-Hilaire, Quebec), here as a function of antecedent moisture conditions and catchment morphology. For the storms observed under dry conditions, larger magnitudes of new water were generated from the three largest catchments attributable to basin morphology, while the storms observed under wet conditions exhibited no consistent pattern, with larger variability among the smaller catchments. The results illustrated the complexity of the influences of antecedent moisture conditions. For the pre-Alpine Erlenbach tributary (0.7 km 2 ) Von Freyberg et al. (2018) showed that pre-event water as a fraction of precipitation was strongly correlated with all measures of antecedent wetness but not with storm characteristics, implying that wet conditions primarily facilitate the mobilization of old (pre-event) water rather than the fast transmission of new (event) water to streamflow, even at a catchment where runoff coefficients can be large. Time series of the natural isotopic composition (2H, 18O) of precipitation and streamwater can provide important insights into ecohydrological phenomena at the catchment scale. However, multi-year, high-frequency isotope datasets are generally scarce, limiting our ability to study highly dynamic short-term ecohydrological processes. Von Freyberg et al. ( 2022) recently presented a four years of daily isotope measurements in streamwater and precipitation at the Alp catchment in Switzerland and two of its tributaries. Therefore, the current study contributes here in three ways: 1. The classic hydrograph separation is embedded in a discretization of the catchment water and tracer mass balance along the event and preevent time axis. 2. A characteristic seasonal variation of event, recent, and pre-recent runoff components is shown for the pre-Alpine Alp catchment (46.4 km 2 ) and two smaller tributaries (Erlenbach, 0.7 km 2 , and Vogelbach, 1.6 km 2 ). 3. Single rain-runoff events of the Erlenbach catchment are analyzed in more detail to visualize and quantify the rapid mobilization of recent water. # II. Methods # a) Separation of n Time Components Consider a control volume, for instance, a catchment in a river basin, with the following bulk water balance: dS(t) dt = J(t) -ET (t) -Q(t)(1) where S is the time evolution of the water storage, J is the precipitation, ET is the evapotranspiration, and Q is the total stream discharge. Let C be a conservative isotope tracer with the following bulk mass balance: d(C S (t)S(t)) dt = C J (t)J(t) -C ET (t)ET (t) -C Q (t)Q(t)(2) where C S , C J , C ET , and C Q are tracer concentrations of the water storage S, and the volumetric flow rates J, ET , and Q. In addition, there are the time points t 0 , t 1 , . . . , t n and time intervals [t 0 , t 1 [, [t 1 , t 2 [, . . . , [t n-1 , t n [ that describe the start and end of n rainfall-runoff events along the time axis. Furthermore, there is a semantic time measure with the intervals e (event) and p (pre-event) that can be moved across the rainfall-runoff events, whereas the interval p is the range of all intervals just before interval e. For instance, the stream discharge Q during event e is composed of water from the current rainfall-runoff event and the prior rainfall-runoff events, such that Q(t) = Q e (t) + Q p (t) C Q (t)Q(t) = C e Q (t)Q e (t) + C p Q (t)Q p (t)(3) where C e Q and C p Q are bulk tracer concentrations in the event and pre-event components, respectively. The same relations can be applied to the physical variables dS(t)/dt, J(t), and ET (t). Therefore, we can write generally for each volumetric flow rate V that # Global Journal of Science Frontier Research ( H ) XXIII Issue VI Version I Year 2023 are usually the estimated physical quantities. When being applied to all physical variables dS(t)/dt, J(t), ET (t), and Q(t), the linear equation system (6) can be regarded as a discretization of the ordinary differential equation system (1) and (2) along ? +1 rainfall events. In the literature, equation (3) is the standard two-component separation model and has been used in many hydrological investigations (Klaus and McDonnell (2013)). Sklash and Farvolden (1979) and Buttle (1994) mentioned the following criteria, which also apply to the iterative separation model (6): 3. The pre-event component maintains a constant isotopic signature in space and time, and if not, any variations can be accounted for. I would like to add another criterion (Criterion 4) that is usually implicitly considered and demands that both the event water V e and pre-event water V p cannot be less than zero or larger than the total runoff V (Liu et al. (2004)). Given the equations above, this is the case if the tracer concentration in the volumetric flow C V is always between that of the event water C e V and preevent water C p V . In the context of separation model (6), it is required for all the backward iterations ? that C e V < C V < C p V ? C p V < C V < C e V(7) Based on the measured knowns V and C V and the estimated tracer concentrations C e V , C e-1 V , . . . , C e-? V (t) C e-? V (t) -C p-? V (t) , ? ? 1 (11) In addition, from linear equation system (6), we obtain the following total mass balance for each volumetric flow rate V : V (t) = ? i=0 V e-i (t) + V p-? (t)(12) Resolving the system of balance equations being introduced in section 2.1 for the mentioned unknowns ( V e , V e-1 . . . V e-? and V p , V p-1 . . . V p-? ), also called end members, requires appropriate estimators for the end member concentrations (C e V , C e-1 V , . . . , C e-? V and C p V , C p-1 V , . . . , C p-? V ). For instance, when referring to the use case of hydrograph separations, which is addressed by the volumetric flow rate variable Q, the isotope signature C e J of precipitation J can be taken as an estimator for the tracer concentration C e Q regarding end member Q e , whereby changes in the isotope composition of precipitation as a result of evapotranspiration ET or changes in the water storage S (e.g. snow pack with sublimation and re-sublimation processes) must be considered. This context becomes visible if balance equations (1) and ( 2) are restricted to the event water fraction, such that # b) Determination of End Member Concentrations and Error Estimation In a similar way, the isotope composition C Q in the discharge Q, which occurs right before a rainfall-runoff event, can be taken as an estimator for the tracer concentration C p Q of the end member Q p , whereby changes in the isotope composition of the bulk pre-event isotope composition because of interception J, evapotranspiration ET , or changes in the water storage S (rapid mobilization of pre-event water) must be considered, which becomes obvious, if balance equations ( 1) and ( 2) are restricted to the pre-event water Q e (t) = J e (t)Q p (t) = J p (t) -ET p (t) -dS p (t)/dt C p Q (t)Q p (t) = C p J (t)J p (t) -C p ET (t)ET p (t) -d(C p S (t)S p (t))/dt(14) The above equations show that the quality of the estimators C e Q , C e-1 Q , . . . , C e-? Q and C p Q , C p-1 Q , . . . , C p-? Q can be improved by increasing the observation rate in all volume flow rates dS(t)/dt, J(t), ET (t), and Q(t), and in fact, there are various approaches in the literature that have directly or indirectly addressed the time variant transformation of tracer mass precipitation via evapotranspiration losses or water storage changes. For instance, by solving the above set of ordinary differential equations (Laudon et In a natural hydrological system, the accurate determination of tracer concentrations C e Q , C e-1 Q , . . . , C e-? Q and C p Q , C p-1 Q , . . . , C p-? Q remains a difficult task, and this is where the error estimators come into play. For the analysis of the field data, the sensitivity of the model plus the input errors of the known variables are usually included in one measure. For instance, Genereux (1998) and Uhlenbrook and Hoeg (2003) have demonstrated this based on analytic expressions for the case of uncorrelated known variables and assumed uncertainties, that is, a classical Gaussian error propagation. Others, for example, Kuczera Gaussian error propagation means that the uncertainty u j of the unknown variable y j (j = 1..n) is related to the uncertainties u i (i = 1..m) of the known # Investigating the Seasonal Variations of Event, Recent, and Pre-Recent Runoff Components in a Pre-Alpine Catchment using Stable Isotopes and an Iterative Hydrograph Separation Approach # Global Journal of Science Frontier Research ( H ) XXIII Issue VI Version I Year 2023 variables x i (i = 1..m) (assumed to be independent from each other) in the following way: u j = ?y j ?x 1 u 1 2 + ?y j ?x 2 u 2 2 + .. + ?y j ?x m u m 2(15) The first-order partial derivatives ?y j ?x i can be collected in the following Jacobian n × m matrix: J(y 1 ..y n ) = ? ? ? ?y 1 ?x 1 ? ? ? ?y 1 ?xm . . . . . . . . . ?yn ?x 1 ? ? ? ?yn ?xm ? ? ?(16) For instance, in case of ? = 0 for the linear equation system ( 6) with the unknown variables V e and V p and known (respectively estimated) variables V , C V , C e V , and C p V , we get the following: J( V e , V p , t) = ? ? ? ?(17) When looking at the denominators of the single Jacobian entries, we can consider the uncertainties u e V (t) and u p V (t) as functions of 1/(C e V (t) -C p V (t)) 2 . Beyond the uncertainties of the known variables, a high model-driven uncertainty is expected for event water isotope concentrations that are close to the corresponding pre-event water isotope concentration. An iterative separation model ( 6) can be used to trace event water over a much longer period after the initial event being studied. For the obtained result, I use the term separated event water response to emphasize the fact that the traced response relates to the water of exactly one rainfall event. Conceptually, it can be related to the more commonly used term travel time distribution. Basically, it can be a time-varying approximation for this. The travel time distribution is the response or breakthrough of an instantaneous, conservative tracer addition over the entire catchment area. It is the probability distribution that can be derived analytically based on the physical assumptions of the system under investigation. By applying a convolution integral, it can balance the tracer inputs and outputs of equation ( 2), as follows (Niemi (1977)): For non-stationary systems, it makes sense to distinguish between two types of probability density functions: the forward travel time distribution and backward travel time distribution (Niemi (1977)). The forward travel time distribution, -? h (?, t in ), is the probability distribution of the travel times ? conditional on the injection time t in of a volume flow (e.g., precipitation). C Q (t) = t -? C J (t in ) h(t -t in )dt in ,(18) # The backward travel time distribution, ?h (?, t), is the probability distribution of the travel times ? conditional on the exit time t of a volume flow (e.g., discharge). When the system is in a steady state (constant input/output fluxes), then the forward and backward travel time distributions collapse into a single probability density function (Niemi (1977) ? - h (t -t in , t)Q(t) = J(t in )?(t in ) - ? h (t -t in , t in ) (19) where ?(t in ) is a partition function describing the fraction of rainfall J(t in ) that ends up as runoff Q. In addition, an age function ? Q can be defined that describes the ratio between the number of water particles with an age in the interval [?, ? + d?] sampled by Q at time t and the amount of particles with the same age stored in the control volume at that time: ? Q (?, t) = ? - h (?, t) h(?, t)(20) where h(?, t) is the probability distribution of the residence times ? of the water particles stored within the control volume at time t. The age function ? Q (?, t) is an interesting quantity in the sense that the tracer concentrations C e Q , C e-1 Q , . . . , C e-? Q and C p Q , C p-1 Q , . . . , C p-? Q from equation system (6) can be represented as a function of the same. For instance, the pre-event concentration C p Q for the event e can be calculated as follows: C p Q (e) = e-1 -? C J (t in )? Q (t -t in , t) h(t -t in , t)dt in(21) The basic procedure to reconstruct the event water response has already been demonstrated in Hoeg (2019), where the single components Q e , Q e-1 . . . Q e-? were arranged in the following way: Assume we are interested in the event water contribution of event 1 during events 2, 3, and 4. In this case, I can arrange one after the other, in which we have the event water Q e of event 1, the last event water Q e-1 of event 2, the second-to-last event water Q e-2 of event 3, and the third-to-last event water Q e-3 of event 4, as illustrated in Figure 1. When referring to the rainfall events J 1 , J 2 , . . . , J ? +1 and the related catchment responses Q 1 , Q 2 , . . . , Q ? +1 , I can define the separated event water response as follows: H [t 1 ,t 2 [ := ? +1 k=1 Q e-k+1 (Q k )(22) This is the volume flow of the water that entered the catchment at interval [t 1 , t 2 [ with rainfall J 1 , which appears in the stream discharge Q during interval [t 1 , t 1+? +1 [, here as a result of the rainfall events J 1 , J 2 , . . . , J ? +1 and the related catchment responses Q 1 , Q 2 , . . . , Q ? +1 . Furthermore, I postulate that the volume-weighted function of the timevarying separated event water response H [t i ,t i+1 [ can be considered an approximation of the time-varying backward travel time distribution, ?h (?, t), on the time interval t ? [t i , t i+? +1 [ when it comes to all water molecules that entered the catchment (the control volume of system (1)) at t in ? [t i , t i+1 [: ? - h (?, t) ? H [t i ,t i+1 [ t i+? t i H [t i ,t i+1 [ (t)dt(23) , respectively ? - h (?, t) ? H e e+? e H e (t)dt (24) regarding the semantic time intervals e, e + 1, . . . , e + ? on ? backward iterations. 2018) mention that the Alp tributaries are wet throughout most of the year. This is due to the high clay content, the low drainable porosity and shallow soils. The water table is generally close to the soil surface, especially in hollows and flatter areas, where the hydraulic gradient is low or at the bottom of hillslopes because of the large amount of water coming from upslope areas. Surface soil moisture measurements show that soil moisture is lowest in the forested ridge sites and highest in the flatter meadow and wetland sites. Total annual precipitation in the Alp valley is strongly controlled by elevation, averaging 1791 mm/year in the flat northern part near the outlet, and roughly 30% more in the mountainous headwaters of the catchment (2300 mm/year). Snowfall comprises up to one-third of the total precipitation in the headwaters of the Alp, although snowfall is frequently interrupted by rainfall during mild periods in winter, with a corresponding occurrence of rain-on-snow events (Rücker et al. ( 2019)). # Investigating the Seasonal The ratio of current to potential evapotranspiration (ETa/ETp ratio) serves as a valuable indicator of water availability for plants. When this ratio falls below 0.8, it suggests an increased likelihood of drought-related impairments (Allgaier Leuch et al. ( 2017)). In the Alp catchment, we can expect ETa/ETp values to range from 0.61 to 0.9 in the valley bottoms, and from 0.81 to 1.00 in the upslope areas. These figures are based on the long-term average spanning from 1981 to 2010. In the present study, a four-years time series (June 2015-May 2019) of daily stable water isotope data in stream water and precipitation is used for the analysis of seasonal variations. The data set is described by Von Freyberg et al. ( 2022), who provided detailed information on the precipitation and streamwater sampling, sample handling and isotope analysis. Streamwater isotopes were measured in the Alp main stream and in two of its tributaries (Erlenbach and Vogelbach). Precipitation isotopes were measured at two grassland locations in the Alp catchment: in the headwaters at 1228 m a.s.l. and near the outlet at 910 m a.s.l. The data set also includes the daily time series of key hydrologic and meteorologic variables, such as daily streamwater and precipitation fluxes, air temperature, relative humidity, and snow depth. To better classify single rainfall runoff events and seasonal relationships, the hydrological data were supplemented with meteorological data (e.g., amount of snowfall, evapotranspiration, soil moisture) from ERA5. The reanalysis product ERA5 has recently been released by European Centre for Medium-Range Weather Forecasts (ECMWF) as part of Copernicus Climate Change Services (Hersbach et al. (2019)). This product covers the period from 1979 to present. # 3. 4. # 5. Q The calculated event, recent, and pre-recent runoff components can be used to estimate time-varying backward travel time distributions. In the following, I demonstrate this on the basis of a simple example. Five rainfall runoff events, taken from a field study in the Erlenbach catchment between September 2016 and October 2017 (Von Freyberg et al. ( 2018)), serve as the examples. The analysis is based on the separation model ( 6) applying three backward iterations (? = 3). To determine the end member concentrations 13) and ( 14) are not explicitly solved; instead, the pre-event water concentrations are taken directly from the measured isotope compositions in the discharge right before the hydrograph rises. The event water concentrations refer to the volume-weighted isotope composition in the respective precipitation event. Q e noC e Q , C e-1 Q , . . . , C e-3 Q and C p Q , C p-1 Q , . . . , C p-3 Q , differential equations ( For calculating the Gaussian standard errors, an uncertainty at the scale of the discharge measurements in the field, u(Q) = 0.001 [m 3 /s], and of the isotope analysis in the laboratory is adopted, that is, u(C Q ) = 0.09 [?], u(C e Q ) = 0.09 [?], u(C p Q ) = 0.09 [?] , and continuously added for each backward iteration: 2018) reported that pre-event water is more efficiently mobilized under wetter conditions, showing u(C e-1 Q ) = 0.18 [?], u(C p-1 Q ) = 0.18 [?], u(C e-2 Q ) = 0.27 [?], u(C p-2 Q ) = 0.27 [?], u(C e-3 Q ) = 0.36 [?], and u(C p-3 Q ) = 0.36 [?]. Figure 7 shows the separated event water response according to equation (22) and Figure 1 for five rainfall runoff events between September 19, 2017 and October 21, 2017. The volume-weighted version of this time-varying response can be considered an approximation of the time-varying backward travel time distribution, as shown in equation (23). In addition, the so-called rapid mobilization of pre-event water can be detected. For instance, in the that the rapid activation of the pre-event water at Erlenbach (even during small storms) can be explained by generally shallow perched groundwater tables in the aquifer overlying the low permeability bedrock. The median Gaussian standard error for the calculated event, recent, and pre-recent runoff components are u(Q e ) = 81.9% / 0.01 mm 10min , u(Q e-1 ) = 123.4% / 0.10 mm 10min , u(Q e-2 ) = 205.2% / 0.18 mm 10min , u(Q e-3 ) = 231.6% / 0.07 mm 10min , and u(Q p-3 ) = 32.3% / 0.18 mm 10min . For better illustration, Figure 8 shows the separated event water response for a single rainfall runoff event (number 3) together with the Gaussian standard error bounds according to equation (15). # Alp In the following, an analysis of the runoff generation processes in terms of event Q e , recent Q r-3 , and pre-recent Q p-3 runoff components is demonstrated for the pre-Alpine Alp catchment (46.4 km 2 ) and two smaller tributaries (Erlenbach, 0.7 km 2 , and Vogelbach, 1.6 km 2 ), whereas recent water is understood as the contribution of the respective three prior rainfall-runoff events, that is To compensate for altitude effects on the catchment scale, the daily isotope concentrations in precipitation considers the available isotopic samples of the two measuring stations and are area weighted based on to the hypsometric curve of the Alp catchment. The analysis is based on the separation model (6) that applies three backward iterations (? = 3), whereas the entire data set is being processed; that is, 147 rainfall runoff events are being investigated for the Alp catchment, 113 rainfall runoff events for the Erlenbach catchment, and 120 rainfall runoff events for the Vogelbach catchment. Again, the end member concentrations Alp discharge and precipitation. In addition, where necessary, the end member concentrations are partially adjusted to ensure that Criterion 4 (equation ( 7)) is continuously fulfilled. This approach does not reduce or increase the size of the relative error, but it does stabilize the solution overall (Hoeg (2021)). Q r-3 = Q e-1 + Q e-2 + Q e-3(25)C e Q , C e-1 Q , . . . , C e-3 Q and C p Q , C p-1 Q , . . Although the absolute Gaussian standard error bounds are moderately low with values between 0.09 mm/d and 0.61 mm/d (see Table 1), the relative errors associated with this hydrograph separation are relatively high, with values between 38.5% and 622.4%. Therefore, the results may not lead to accurate quantitative conclusions. Nevertheless, qualitative statements about the seasonal runoff formation in the Alpine catchment area and its tributaries Erlenbach and Vogelbach are possible, as I will show below. To visualize the seasonal variations, I calculate the Pardé coefficients for each runoff component, which is the quotient of the long-term (four years) average monthly discharge and long-term (four years) average annual discharge. The precipitation and evapotranspiration regimes are added to better evaluate and classify the monthly and seasonal changes of the runoff components, which are: stream discharge Q, event water Q e , recent water Q r-3 , and pre-recent water Q p-3 . Evapotranspiration data are taken from the reanalysis product ERA5 of the Copernicus Climate Change Service. For the Alp catchment and its tributaries, a nivale runoff and precipitation regime can [-] Alp be shown in Figure 9. In relation to the stream discharge Q, the event water component Q e has its maximum (>32%) during August in general, as shown in Figure 10 and Figure 11. The highest recent water fractions Q r-3 (50-59%) can be expected at the beginning of winter (December, January), whereas the lowest fractions (7-14%) are found at the end of autumn (November). In return, higher pre-recent water fractions Q p-3 can be found in the summer (July, 60-82%) and in the middle of autumn (October, 64-74%). The two tributaries expose an additional peak of pre-recent water in the snowmelt season (April, 47-62%), during which the lowest event water fractions (6-8%) can be expected. For the Alp catchment, the event water component (29-32%) and recent water component (56%) clearly exceed the pre-recent component (17-21%) in August and September after the month July, which is characterized by relatively low precipitation, high evapotranspiration, and low soil moisture. A sensitivity analysis regarding the evapotranspiration rates ET from ERA5 confirms the calculated seasonal variations of the Pardé coefficients both for the Alp catchment and its tributaries, Erlenbach and Vogelbach. The analysis, shown in Figure 12, indicates a shift in proportions towards Added monthly Pardé coefficients (1: January -12: December) of event water Q e , recent water Q r-3 = Q e-1 +Q e-2 +Q e-3 , and pre-recent water Q p-3 for the Erlenbach, Vogelbach, and Alp catchments in relation to the monthly Pardé coefficients of stream discharge (Q). Monthly Pardé coefficients (1: January -12: December) of precipitation (J), and evapotranspiration (ET), based on the four-year investigation between June 2015 and May 2019, and an iterative extension of the standard two-component hydrograph separation method. Sensitivity analysis to assess the impact of ET on the daily end member concentrations (C e Q , C e-1 Q , . . . , C e-3 Q ). This analysis involves considering a rate factor f ET , which can take values of 0.25 or 0.5. To first-order approximate the daily end member concentrations, the daily precipitation rates (J) are directly reduced by 25% or 50% of the daily evapotranspiration rates, similar to equation (13). The evapotranspiration data used in this analysis are obtained from the reanalysis product ERA5. 2022)) is unlikely to lead to a greater mobilization and availability of pre-recent runoff components in this area. Consequently, the Median of the relative and absolute Gaussian standard errors (in % and mm d ) for event water Q e , recent water Q r-3 , and pre-recent water Q p-3 in the Erlenbach, Vogelbach and Alp catchment between June 2015 and May 2019, here based on an iterative extension of the standard two-component hydrograph separation method with the uncertainties of stream discharge, u(Q) = 0.001 [m 3 /s], and of the isotope analysis u( C Q ) = 0.09 [?], u(C e Q ) = 0.09 [?], u(C p Q ) = 0.09 [?], u(C e-1 Q ) = 0.18 [?], u(C p-1 Q ) = 0.18 [?], u(C e-2 Q ) = 0.27 [?], u(C p-2 Q ) = 0.27 [?], u(C e-3 Q ) = 0. # V. Conclusions In the present study, I have shown that the classic two-component hydrograph separation can be iteratively embedded in a discretization of the catchment water and tracer mass balance along the event and pre-event time axis. With this method, it is possible, for instance, to analyze and quantify the rapid mobilization of recent water for single (high-frequency measured) rainfall-runoff events, and to estimate time-varying backward travel time distributions. When applied to longer time series of daily stable water isotope data in stream water and precipitation, the method can be used to analyze seasonal variations of event, recent, and pre-recent runoff components. In relation to the Alp catchment, the relatively high fractions of event water and recent water in stream discharge during August and September represents a rather unexpected result that certainly requires further investigation, but shows that fundamental assumptions, used for instance in runoff recession analysis, need to be questioned and that the sensitivity of the catchment water balance in response to drought situations could be greater than expected. The author would like to thank Jana von Freyberg (ETH -Swiss Federal Institute of Technology Zurich, Department of Environmental Systems Science) for providing the additional excerpt of the high-frequency measurement series from the Erlenbach catchment. My special thanks goes to chief editor Marie V. Carlsen, assistant editor Marian C. Miller, and two anonymous reviewers, who contributed to this work with numerous improvement proposals. potential to incorporate this approach into hydrological models to improve their accuracy and predictive capabilities. Moreover, there are opportunities to apply the method in more complex and detailed ways, potentially uncovering new insights in catchment hydrology. This could include, for instance, more detailed analyses of the impacts of extreme weather events or the effects of different land use practices. ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 0 0 0 0 0 0 C e V C p V 0 0 0 0 0 0 0 -1 1 1 0 0 0 0 0 -C p V C e-1 V C p-1 V 0 0 0 0 0 0 0 -1 1 1 0 0 0 0 0 -C p-1 V C e-2 V C p-2 V 0 0 0 0 0 0 0 -1 1 1 0 0 0 0 0 -C p-2 V C e-3 V C p-3 V ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? V e V p V e-1 V p-1 V e-2 V p-2 V e-3 V p-3 ? ? ? ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? ? ? V C V V 0 0 0 0 0 0 ? ? ? ? ? ? ? ? ? ? ? ? (B.![al. (2002)), by using transfer functions (Weiler et al. (2003), Iorgulescu et al. (2007), and Segura et al. (2012)), age functions (Botter et al. (2011)) or storage selection functions (Harman (2015), Benettin et al. (2017)), or by applying the correlations between tracer fluctuations in precipitation, evapotranspiration and discharge regarding longer observation periods (Kirchner (2019), Kirchner and Allen (2020))).](image-2.png "") ![and Parent (1998), Joerin et al. (2002), Weiler et al. (2003), Iorgulescu et al. (2007), Segura et al. (2012), and Borriero et al. (2023), approximated the expected values based on designed field scenarios and the law of large numbers, which is better known as the Monte Carlo method.](image-3.png "") 712![Figure1: Basic procedure to reconstruct the event water response, here showing the event water contribution of event 1 during events 2, 3, and 4. It is arranged one after the other: the contribution of event water Q e (t) during event 1, the contribution of the last event water Q e-1 (t) during event 2, the contribution of the second-to-last event water Q e-2 (t) during event 3, and the contribution of the third-to-last event water Q e-3 (t) during event4where h(?) is the probability distribution of the travel time ? and t in is the injection time of the tracer. The residence time, travel time, and life expectancy of water particles, along with the associated constituents flowing through watersheds, are three related quantities whose meaning has been well defined (see Rigon et al. (2016)). For instance, Maloszewski and Zuber (1982) introduced a similar convolution integral, as shown in equation (18), but they used the quantity exit age instead of injection time. On the catchment scale, travel time distributions can serve as fundamental catchment descriptors, revealing information about storage, flow pathways, and sources of water in a single characteristic(McGuire and McDonnell (2006)). Assuming steady-state conditions, travel time distributions are often interpreted and applied as time invariant, for example, as a mean over the period under investigation. In addition, travel time distributions are usually inferred by using lumped parameter models that simplify the description of a spatially distributed catchment behavior. Of course, there is evidence and knowledge to the contrary(Hrachowitz et al. (2010), McDonnell et al. (2010), Botter (2012), Heidbüchel et al. (2013)).](image-4.png "7 ©Figure 1 :Figure 2 :") C e Q (t)Q e (t) = C e J (t)J e (t) -C e ET (t)ET e (t) -d(C e S (t)S e (t))/dt(13)fraction, such that Year 202312Global Journal of Science Frontier Research ( H ) XXIII Issue VI Version I© 2023 Global Journals IV. Results and Discussion0.8Erlenbach00.7Q [mm/10min]0.3 0.4 0.5 0.65 10J [mm/10min]0.2150.109/2410/0110/0810/15-6-7C Q-8C J-918-10-11-1209/2410/0110/0810/15Year 202313.1 Q e no.2 Q e no.3 Q e no.4 Q e no.5Global Journal of Science Frontier Research ( H ) XXIII Issue VI Version I0.20.150.10.050© 2023 Global Journals Figure 7: Erlenbach Catchment: Precipitation [mm/10 minutes], Stream Discharge [mm/10 minutes], and Separated Event Water Response [mm/10 minutes] for Five Rainfall Runoff Events Between September 19, 2017 -October 21, a) Estimating the Time-Varying Backward Travel Time Distribution for Single Events Year 202317Global Journal of Science Frontier Research ( H ) XXIII Issue VI Version I© 2023 Global Journals 136 [?], and u(C p-3 Q ) = 0.36 [?]. (Scherler et al. (2016) ). ? ? ? 10. Appendix A.2. Greek Symbols partition function backward iteration travel time ? age functionGlobal Journal of Science Frontier Research ( H ) XXIII Issue VI Version I Year 2023 21© 2023 Global Journals © 2023 Global Journals © 2023 Global Journals ## Acknowledgements ## Appendix A.1. Latin Symbols The data used for this study can be retrieved via Von Freyberg et al. (2018) ## References Références Referencias The iterative separation model ( 6) can be represented in the form of a linear equation system for ? = 3, as follows: * The seasonal origins of streamwater in switzerland STAllen JVon Freyberg MWeiler GRGoldsmith JWKirchner 10.1029/2019GL084552 Geophysical Research Letters 46 2019 * Der Schweizer Wald im Klimawandel: Welche Entwicklungen kommen auf uns zu? BAllgaier-Leuch KStreit PBrang Merkblatt fu?r die Praxis -Swiss Federal Institute for Forest, Snow and Landscape Research 59 2017 * Northern hemisphere drought risk in a warming climate DFBalting AAghakouchak GLohmann MIonita 10.1038/s41612-021-00218-2 Climate and Atmospheric Science 4 2021 * Using sas functions and highresolution isotope data to unravel travel time distributions in headwater catchments PBenettin CSoulsby CBirkel DTetzlaff GBotter ARinaldo 10.1002/2016WR020117 Water Resources Research 53 2017 * Uncertainty in water transit time estimation with storage selection functions and tracer data interpolation ABorriero RKumar TVNguyen JHFleckenstein SRLutz 10.5194/hess-27-2989-2023 Hydrology and Earth System Sciences 27 2023 * Catchment mixing processes and travel time distributions GBotter 10.1029/2011WR011160 Water Resources Research 48 2012 * Catchment residence and travel time distributions: The master equation GBotter EBertuzzo ARinaldo 10.1029/2011GL047666 Geophysical Research Letters 38 2011 * Isotope hydrograph separations and rapid delivery of pre-event water from drainage basins JButtle 10.1177/030913339401800102 Progress in Physical Geography: Earth and Environment 18 1994 * Snowmelt runoff from measurements of tritium and oxygen-18 TDincer BRPayne TFlorkowski JMartinec ETongiorgi 10.1029/WR006i001p00110 Water Resources Research 6 1970 * Global reanalysis: goodbye era-interim, hello era5. European Centre for Medium-Range Weather Forecasts (ECMWF) HHersbach BBell PBerrisford AHoranyi JMSabater JNicolas RRadu DSchepers ASimmons CSoci DDee 10.21957/vf291hehd7 2019 159 * On the balance equations and error estimators for separating n time components of runoff with one stable isotope tracer SHoeg 10.1029/2019WR025555 Water Resources Research 55 2019 * Benchmark tests for separating n time components of runoff with one stable isotope tracer SHoeg 10.5194/hess-2021-213 Hydrology and Earth System Sciences Discussions 2021 2021 * Gamma distribution models for transit time estimation in catchments: Physical interpretation of parameters and implications for time-variant transit time assessment MHrachowitz CSoulsby DTetzlaff IAMalcolm GSchoups 10.1029/2010WR009148 Water Resources Research 46 2010 * Flow, mixing, and displacement in using a data-based hydrochemical model to predict conservative tracer data IIorgulescu KJBeven AMusy 10.1029/2005WR004019 Water Resources Research 43 2007 * Antecedent moisture conditions and catchment morphology as controls on spatial patterns of runoff generation in small forest catchments AJames NRoulet 10.1016/j.jhydrol.2009.08.039 Journal of Hydrology 377 2009 * Global isotope hydrogeology-review SJasechko Reviews of Geo-physics 57 2019 * Uncertainty in hydrograph separations based on geochemical mixing models CJoerin KBeven IIorgulescu AMusy 10.1016/S0022-1694(01)00509-1 Journal of Hydrology 255 2002 * Quantifying new water fractions and transit time distributions using ensemble hydrograph separation: theory and benchmark tests JWKirchner 10.5194/hess-23-303-2019 Hydrology and Earth System Sciences 23 2019 * Seasonal partitioning of precipitation between streamflow and evapotranspiration, inferred from end-member splitting analysis JWKirchner STAllen 10.5194/hess-24-17-2020 Hydrology and Earth System Sciences 24 2020 * Hydrograph separation using stable isotopes: Review and evaluation JKlaus JMcdonnell 10.1016/j.jhydrol.2013.09.006 Journal of Hydrology 505 2013 * Monte carlo assessment of parameter uncertainty in conceptual catchment models: the metropolis algorithm GKuczera EParent 10.1016/S0022-1694(98)00198-X Journal of Hydrology 211 1998 * Oxy-gen 18 fractionation during snowmelt: Implications for spring flood hydrograph separation HLaudon HFHemond RKrouse KHBishop 10.1029/2002WR001510 Water Resources Research 38 2002 * Drought resistance enhanced by tree species diversity in global forests DLiu TWang JPenuelas SPiao 10.1038/s41561-022-01026-w Nature Geoscience 15 2022 * Source waters and flow paths in an alpine catchment, colorado front range, united states FLiu MWWilliams NCaine 10.1029/2004WR003076 Water Resources Research 40 2004 * Determining the turnover time of groundwater systems with the aid of environmental tracers 1. models and their applicability PMaloszewski AZuber 10.1016/0022-1694 Journal of Hydrology 57 1982 * New insights into the run-off mechanism by environmental isotopes JMartinec USiegenthaler HOeschger ETongiorgi Symposium Isotope Techniques in Groundwater Hydrology International Atomic Energy Agency (IAEA) 1974 * How old is streamwater? open questions in catchment transit time conceptualization, modelling and analysis JJMcdonnell KMcguire PAggarwal KJBeven DBiondi GDestouni SDunn AJames JKirchner PKraft SLyon PMal-Oszewski BNewman LPfister ARinaldo ARodhe TSayama JSeibert KSolomon CSoulsby MStewart DTetzlaff CTobin PTroch MWeiler AWestern AWo?rman SWrede 10.1002/hyp.7796 Hydrological Processes 24 2010 * A review and evaluation of catchment transit time modeling KJMcguire JJMcdonnell Journal of Hydrology 330 2006 * Potential of land-based climate change mitigation strategies on abandoned cropland MHGvein XHu JSNaess MD BWatanabe OCavalett MMalbranque GKindermann FCherubini 10.1038/s43247-023-00696-7 Communications Earth and Environment 4 2023 * Time-variable transit time distributions and transport: Theory and application to storagedependent transport of chloride in a watershed CJHarman 10.1002/2014WR015707 Water Resources Research 51 2015 * Separating physical and me-teorological controls of variable transit times in zero-order catchments IHeidbu?chel PATroch SWLyon 10.1002/2012WR013149 Water Resources Research 49 2013 * Age-ranked hydrological budgets and a travel time description of catchment hydrology RRigon MBancheri TRGreen 10.5194/hess-20-4929-2016 Hydrology and Earth System Sciences 20 2016 * Catchment travel time distributions and water flow in soils ARinaldo KJBeven EBertuzzo LNicotina JDavies AFiori DRusso GBotter 10.1029/2011WR010478 Water Resources Research 47 2011 * Monitoring snowpack outflow volumes and their isotopic composition to better understand streamflow generation during rain-on-snow events ARu?cker SBoss JWKirchner JVon Freyberg 10.5194/hess-23-2983-2019 Hydrology and Earth System Sciences 23 2019 * The tim-ing of unprecedented hydrological drought under climate change YSatoh KYoshimura YPokhrel HKim HShiogama TYokohata NHanasaki YWada PBurek EByers HMSchmied DGerten SOstberg SNGosling JE SBoulange TOki 10.1038/s41467-022-30729-2 Nature Communications 13 2022 * Wald im Klimawandel. Grundlagen fu?r Adaptations strategien: Wasserhaushalt von Waldern bei zunehmender Trockenheit MScherler JRemund LWalthert 2016 Swiss Federal Institute for Forest, Snow and Landscape Research * Scaling relationships for event water contributions and transit times in small-forested catchments in eastern quebec CSegura ALJames DLazzati NTRoulet 10.1029/2012WR011890 Water Resources Research 48 2012 * Excess forest mortality is consistently linked to drought across europe CSenf ABuras CSZang ARammig RSeidl 10.1038/s41467-020-19924-1 Nature Communications 11 2020 * Canopy mortality has doubled in europe's temperate forests over the last three decades CSenf DPflugmacher YZhiqiang JSebald JKnorn MNeumann PHostert RSeidl 10.1038/s41467-018-07539-6 Nature Communications 9 2018 * The role of groundwater in storm runoff MGSklash RNFarvolden 10.1016/0022-1694(79)90164-1 Journal of Hydrology 43 1979 * Forest health and global change STrumbore PBrando HHartmann 10.1126/science.aac6759 Science 349 2015 * Quantifying uncertainties in tracer-based hydrograph separations: a case study for two-, three-and five-component hydrograph separations in a mountainous catchment SUhlenbrook SHoeg 10.1002/hyp.1134 Hydrological Processes 17 2003 * Runoff generation in a prealpine catchment: A discussion between a tracer and a shallow groundwater hydrologist HVan Meerveld BFischer MRinderer MSt Ähli JSeibert 10.18172/cig.3349 Cuadernos de Investigaci´on Ge-ogr ´afica 44 2018 * Four years of daily stable water isotope data in stream water and precipitation from three swiss catchments JVon Freyberg ARu¨cker MZappa ASchlumpf BStuder JWKirchner 10.1038/s41597-022-01148-1 Scientific Data 9 2022 * Studying catchment storm response using eventand pre-event-water volumes as fractions of precipitation rather than discharge JVon Freyberg BStuder MRinderer JWKirchner 10.5194/hess-22-5847-2018 Hydrology and Earth System Sciences 22 2018 * How does rainfall become runoff? a combined tracer and runoff transfer function approach MWeiler BLMcglynn KJMcguire JJMcdonnell 10.1029/2003WR002331 Water Resources Research 39 2003 * Climate variability drives recent tree mortality in europe MNeumann VMues AMoreno HHasenauer RSeidl 10.1111/gcb.13724 Global Change Biology 23 2017 * Residence time distributions of variable flow processes AJNiemi 10.1016/0020-708X(77)90026-6 The International Journal of Applied Radiation and Isotopes 28 1977 * Determination of the ground-water component of peak discharge from the chemistry of total runoff GFPinder JFJones 10.1029/WR005i002p00438 Water Resources Research 5 1969