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

as


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.

a) Separation of n Time Components
Consider a control volume, for instance, a catchment in a river basin, with the following bulk water balance: 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: 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 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 Furthermore, an iterative balance of the catchment event and pre-event mass and volume flows can be formulated along the time axis; that is, the pre-event water fraction of each rainfall-runoff event is entirely composed of the event water and pre-event water of the previous event.Therefore, for a first backward iteration, we have the following equations: Here, the event components V e , V e−1 . . .V e−τ and pre-event components V p , V p−1 . . .V p−τ are usually unknowns, whereas the volumetric flow rate V and its tracer concentration C V are usually measured physical quantities.The tracer concentrations of the event components 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): 1.The isotopic composition of the event component is significantly different from that of the pre-event component.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 Based on the measured knowns V and C V and the estimated tracer concentrations C e V , C e−1 V , . . ., C e−τ V and C p V , C p−1 V , . . ., C p−τ V , we can iteratively derive for τ ∈ N backward iterations for the following solutions: In addition, from linear equation system (6), we obtain the following total mass balance for each volumetric flow rate V : 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 The above equations show that the quality of the estimators 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 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

Global Journal of Science Frontier Research ( H ) XXIII Issue VI Version I
Year 2023

6
© 2023 Global Journals variables x i (i = 1..m) (assumed to be independent from each other) in the following way: The first-order partial derivatives ∂y j ∂x i can be collected in the following Jacobian n × m matrix: 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: 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).

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) 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: 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 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: 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: 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 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 [: , respectively 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   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.

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.
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, , and continuously added for each backward iteration:  magnified section of event no. 4, event water from the previous events 1-3 still occur at the catchment outlet.Von Freyberg et al. (2018) reported that pre-event water is more efficiently mobilized under wetter conditions, showing 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)).
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  ).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.the event water component (up to 5% in August) and pre-recent water component (up to 20% in August) when daily evapotranspiration rates are considered up to 50% when calculating daily end member concentrations C e Q , C e−1 Q , . . ., C e−3 Q .However, for the Alp catchment, the event water component and recent water component still exceed the pre-recent component in August and September.
The latter result is noteworthy because, on the one hand, it confirms that considerable fractions of summer precipitation become streamflow, even though evapotranspiration fluxes are much larger in the summer (Allen et al. ( 2019)); on the other hand, it shows that a certain limit appears to have been reached in the Alp catchment for the months of August and September.A further increase of hydrological drought in central Europe (Balting et  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(  2016)).

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.Moving forward, there are opportunities to further refine the method, apply it to new contexts, or integrate it with other techniques.For example, the method can be used to examine the event, recent, and pre-recent components in the evapotranspiration (volume flow ET ) or in the groundwater recharge (volume flow dS/dt) of a catchment.In addition, there is the

Acknowledgements
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.

7 ©Figure 1 :Figure 2 :
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)).

10 © 11 ©
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 Site and Data According to Von Freyberg et al. (2022) the 46.4 km 2 Alp catchment is located near the city of Einsiedeln in central Switzerland (Figure2).The catchment spans an elevation range of 1058 m, with the outlet at 840 m above sea level and highest summit (Grosser Mythen) at 1898 m above sea level.The average slope of the Alp catchment is 16 • with a flat valley bottom and very steep slopes of up to 75 • at the south-western catchment boundary.The 0.7 km 2 Erlenbach catchment is located on the eastern side of the Alp valley (elevation range from 1080 to 1520 m above sea level), whereas the 1.6 km 2 Vogelbach catchment is located on the western side (elevation range from 940 to 1480 m above sea level).The bedrock geology of the Alp catchment consists of tertiary flysch (sandstone, limestone, clays, and marls) and subalpine molasse (conglomerates, sandstone, and marls); the valley bottoms are overlain by gravel and 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 2023 Global Journals

Figure 3 :Figure 4 :
Figure 3: Alp catchment: Four years time series (June 2015-May 2019) of daily stream discharge [mm/d] and precipitation [mm/d], temperature [ • C], and snow depth [cm], as well as δ 18 O composition in discharge and precipitation [‰] (Von Freyberg et al. (2022)).Tracer concentrations in precipitation are area weighted according to the hypsometric curve.The snow depths are retrieved from ERA5 as part of the Copernicus Climate Change Services (Hersbach et al. (2019)

Figure 5 :
Figure 5: Vogelbach catchment: Four years time series (June 2015-May 2019) of daily stream discharge [mm/d] and precipitation [mm/d], temperature [ • C], and snow depth [cm], as well as δ 18 O composition in discharge and precipitation [‰] (Von Freyberg et al. (2022)).Tracer concentrations in precipitation are area weighted according to the hypso-metric curve

Figures 3 , 4 ,
and 5 show excerpts of the data set published via Von Freyberg et al. (2022) regarding the Alp main stream and Erlenbach and Vogelbach tributaries.In addition, Figure 6 shows an excerpt of a data set with high frequency measurements (10-minute time interval) published by Von Freyberg et al. (2018) for the Erlenbach tributary.

Figure 8 :
Figure 8: Erlenbach catchment: Precipitation [mm/10 minutes], stream discharge [mm/10 minutes], and separated event water response [mm/10 minutes] (inclusive Gaussian standard error bounds) for rainfall runoff event number 3 between September 19, 2017 -October 21, 2017 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 ApproachGlobal Journal of Science Frontier Research ( H ) XXIII Issue VI Version IYear 2023

Figure 9 :
Figure 9: Monthly Pard é coefficients (1: January -12: December) of stream discharge Q, event water Q e , recent water Q r−3 , pre-recent water Q p−3 , precipitation (J), and evapotranspiration (ET) for the Erlenbach, Vogelbach, and Alp catchments, based on the four-year investigation between June 2015 and May 2019, and an iterative extension of the standard two-component hydrograph separation method.Evapotranspiration data are taken from the reanalysis product ERA5 b) Seasonal Variations of Event, Recent and Pre-Recent Runoff Components

3 Q
. ., C p−are taken directly from the concentrations in the stream 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

Figure 10 :
Figure 10: Monthly Pard é coefficients (1: January -12: December) of event water Q e , recent water Q r−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) in %, based on the four-year investigation between June 2015 and May 2019, and an iterative extension of the standard two-component hydrograph separation method 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 ApproachGlobal Journal of Science Frontier Research ( H ) XXIII Issue VI Version I

Figure 11 :
Figure 11: 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 ´e 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.Evapotranspiration data are taken from the reanalysis product ERA5

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
2. The event component maintains a constant isotopic signature in space and time, and if not, any variations can be accounted for.

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
Freyberg et al. (2022), who provided detailed information on the precipitation and

2017 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
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

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
al. (2021), Satoh et al. (2022)) is unlikely to lead to a greater mobilization and availability of pre-recent runoff components in this area.Consequently, the

Table 1 :
Senf et al. (2020)t water and recent water for the catchment water balance will continue to increase, and with lower rainfall and/or higher temperatures in the summer months, the situation for vegetation, which relies on the remaining summer precipitation, might become even worse than observed in the past few years and decades, as noted bySenf et al. (2018), andSenf et al. (2020).If climate projections prove to be accurate, there will be a significant increase in drought conditions across extensive areas of the Swiss forest by the conclusion of the 21st century.Nevertheless, the duration of drought will increase the most by the year 2100 in those forest site types that are already relatively dry today.The least vulnerable are forest sites in areas with very high precipitation and cool temperatures, on deep or hydromorphic soils.

1) 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 IYear 2023