Search

Methodological advancements have been made in in situ observations of water stable isotopes that have provided valuable insights into ecohydrological processes. The continuous measurement capabilities of laser-based analyzers allow for high temporal resolutions and non-destructive minimally invasive study designs of such in situ approaches. However, isotope analyzers are expensive, heavy, and require shelter and access to electrical power, which impedes many in situ assays. Therefore, we developed a new inexpensive technique to collect discrete water vapor samples in the field via diffusion-tight inflatable bags that can later be analyzed in the lab. In a series of structured experiments, we tested different procedural settings, bag materials, and closure types for diffusion tightness during storage as well as for practical handling during filling and extraction. To facilitate reuse of sampling bags, we present a conditioning procedure using ambient air as primer. In order to validate our method, direct measurements through hydrophobic in situ probes were compared to repeated measurements of vapor sampled with our bags from the same source. All steps are summarized in a detailed standard operating procedure (SOP). This procedure represents the preparation and measurement of calibration and validation vapor standards necessary for processing of unknown field-collected vapor samples in the foreseen application. By performing pertinent calibration procedures, accuracy was better than 0.4 ‰ for δ18O and 1.9 ‰ for δ2H after 1 d of storage. Our technique is particularly suitable when used in combination with minimally invasive water vapor sampling in situ probes that have already been employed for soils and tree xylem. It is an important step towards minimally invasive monitoring of stable isotope distributions and also time series in virtually undisturbed soils and trees without the need to have an analyzer in the field. It is therefore a promising tool for many applications in ecohydrology and meteorology.

The vast majority of studies investigating the source depths in the soil of root water uptake with the help of stable water isotopes implicitly assumes that the isotopic signatures of root water uptake and xylem water are identical. In this study we show that this basic assumption is not necessarily valid, since water transport within a plant's xylem is not instantaneous. However, to our knowledge, no study has yet tried to explicitly assess the distribution of water transport velocities within the xylem. With a dye tracer experiment, we were able to visualize how the transport of water through the xylem happens at a wide range of velocities which are distributed unequally throughout the xylem. In an additional virtual experiment we could show that, due to the unequal distribution of transport velocities throughout the xylem, different sampling approaches of stable water isotopes might effectively lead to xylem water samples with different underlying age distributions.

Soil hydrologic processes play an important role in the hydro-pedo-geomorphological feedback cycle of landscape evolution. Soil properties and subsurface flow paths both change over time, but due to a lack of observations, subsurface water flow paths are often not properly represented in soil and landscape evolution models. We investigated the evolution of subsurface flow paths across a soil chronosequence in the calcareous glacier forefield at the Griessfirn glacier in the Swiss Alps. Young soils developed from calcareous parent material usually have a high pH value, which likely affects vegetation development and pedogenesis and thus the evolution of subsurface flow paths. We chose four glacial moraines of different ages (110, 160, 4 900, and 13 500 years) and conducted sprinkling experiments with the dye tracer Brilliant Blue on three plots at each moraine. Each plot was divided into three equal subplots, and dyed water was applied with three different irrigation intensities (20, 40, and 60 mm h−1) and an irrigation amount of 40 mm. Subsequent excavation of soil profiles enabled the tracing of subsurface flow paths. A change in flow types with increasing moraine age was observed from a rather homogeneous matrix flow at 110 and 160 years to heterogeneous matrix and finger-shaped flow at 4 900 and 13 500 years. However, the proportion of preferential flow paths is not necessarily directly related to the moraine age but rather to soil properties such as texture, soil layering, organic matter content, and vegetation characteristics such as root length density and biomass. Irrigation intensity had an effect on the number of finger-shaped flow paths at the two old moraines. We also found that flow paths in this calcareous material evolved differently compared to a previous study in siliceous material, which emphasizes the importance of parent material for flow path evolution. Our study provides a rare systematic dataset and observations on the evolution of vertical subsurface flow paths in calcareous soils, which is useful to improve their representation in the context of landscape evolution modeling.

Spatially explicit quantification on design storms is essential for flood risk assessment and planning. Due to the limited temporal data availability from weather radar data, design storms are usually estimated on the basis of rainfall records of a few precipitation stations only that have a substantially long time coverage. To achieve a regional picture, these station-based estimates are spatially interpolated, incorporating a large source of uncertainty due to the typical low station density, in particular for short event durations. In this study we present a method to estimate spatially explicit design storms with a return period of up to 100 years on the basis of statistically extended weather radar precipitation estimates, based on the ideas of regional frequency analyses and subsequent bias correction. Associated uncertainties are quantified using an ensemble-sampling approach and event-based bootstrapping. With the resulting dataset, we compile spatially explicit design storms for various return periods and event durations for the federal state of Baden Württemberg, Germany. We compare our findings with two reference datasets based on interpolated station estimates. We find that the transition in the spatial patterns of the design storms from a rather random (short-duration events, 15 min) to a more structured, orographically influenced pattern (long-duration events, 24 h) seems to be much more realistic in the weather-radar-based product. However, the absolute magnitude of the design storms, although bias-corrected, is still generally lower in the weather radar product, which should be addressed in future studies in more detail.

Intermittent streams represent a substantial part of the total stream network, and their occurrence is expected to increase due to climate change. Thus, it is of high relevance to provide detailed information on the temporal and spatial controls of streamflow intermittency to support management decisions. This study presents an event-based analysis of streamflow responses in intermittent streams in a mesoscale catchment with a temperate climate. Based on the streamflow responses, precipitation events were classified into flow or no-flow classes. Response controls like precipitation, soil moisture, and temperature were used as predictors in a random forest model to identify the temporally changing factors that explain streamflow intermittency at the event scale. Soil moisture was the most important predictor, but the predictor importance varied with the geology in the catchment. Streamflow responses in the slate geology were controlled by soil moisture in the shallow and deep soil layers, while streamflow in the marl geology was primarily controlled by soil moisture in the upper soil layer. Streamflow responses in catchments underlain by both marl and sandstone were dependent on soil moisture, whereas streamflow in the only catchment with a pure sandstone geology depended on precipitation characteristics. In all slate and marl catchments, streamflow intermittency also varied with soil temperature, which is probably a proxy for seasonal changes in evapotranspiration and an indicator of freezing conditions. Our findings underline the importance of using high temporal resolution data and tailored event definitions that account for the fast changes between flow/no flow in intermittent streams to identify streamflow controls at the event scale.

The direct vapor equilibration laser spectrometry (DVE-LS) method has been developed for obtaining matrix-bound water stable isotope data in soils, the critical zone, and bedrock, deriving therefrom subsurface water flow and transport processes and, ultimately, characterizing, for example, groundwater recharge and vulnerability. Recently, DVE-LS has been increasingly adopted due to its possible high sample throughput, relative simplicity, and cost-efficiency. However, this has come at the cost of a non-unified standard operation protocol (SOP), and several contradictory suggestions regarding protocol details do exist which have not been resolved to date. Particularly, sample container material and equilibration times have not yet been agreed upon. Beside practical constraints, this often limits DVE-LS applicability to interpreting relative isotope dynamics instead of absolute values. It also prevents data comparability among studies or laboratories, and several previous comparisons of DVE-LS with other, more traditional approaches of water extraction and subsequent stable isotope analysis yielded significant discrepancies for various sample matrices and physical states. In a series of empirical tests, we scrutinized the controversial DVE-LS protocol details. Specifically, we tested 10 different easily available and cost-efficient inflatable bags previously employed or potentially suitable for DVE-LS sample collection and equilibration. In storage tests similar to the DVE-LS equilibration process but lasting several weeks, we quickly found heat-sealed bags made of laminated aluminum (Al) sheets to be superior by several orders of magnitude over more frequently used freezer bags in terms of evaporation safety and accompanying adverse isotope effects. For the first time, Al-laminated bags allow the applied equilibration time to be adapted exclusively to sample requirements instead of accepting reduced data quality in a trade-off with material shortcomings. Based on detailed physical considerations, we further describe how to calculate the minimum available container headspace and sample-contained liquid water volume and how their ratio affects analytical precision and accuracy. We are confident that these guidelines will expand DVE-LS applicability and improve data quality and comparability among studies and laboratories by contributing to a more unified, physically well-founded SOP based on more appropriate components.

A better understanding of the reasons why hydrological model performance is unsatisfying represents a crucial part of meaningful model evaluation. However, current evaluation efforts are mostly based on aggregated efficiency measures such as Kling–Gupta efficiency (KGE) or Nash–Sutcliffe efficiency (NSE). These aggregated measures provide a relative gradation of model performance. Especially in the case of a weak model performance it is important to identify the different errors which may have caused such unsatisfactory predictions. These errors may originate from the model parameters, the model structure, and/or the input data. In order to provide more insight, we define three types of errors which may be related to their source: constant error (e.g. caused by consistent input data error such as precipitation), dynamic error (e.g. structural model errors such as a deficient storage routine) and timing error (e.g. caused by input data errors or deficient model routines/parameters). Based on these types of errors, we propose the novel diagnostic efficiency (DE) measure, which accounts for these three error types. The disaggregation of DE into its three metric terms can be visualized in a plain radial space using diagnostic polar plots. A major advantage of this visualization technique is that error contributions can be clearly differentiated. In order to provide a proof of concept, we first generated time series artificially with the three different error types (i.e. simulations are surrogated by manipulating observations). By computing DE and the related diagnostic polar plots for the reproduced errors, we could then supply evidence for the concept. Finally, we tested the applicability of our approach for a modelling example. For a particular catchment, we compared streamflow simulations realized with different parameter sets to the observed streamflow. For this modelling example, the diagnostic polar plot suggests that dynamic errors explain the overall error to a large extent. The proposed evaluation approach provides a diagnostic tool for model developers and model users and the diagnostic polar plot facilitates interpretation of the proposed performance measure as well as a relative gradation of model performance similar to the well-established efficiency measures in hydrology.

Hillslope–stream connectivity controls runoff generation, during events and during baseflow conditions. However, assessing subsurface connectivity is a challenging task, as it occurs in the hidden subsurface domain where water flow can not be easily observed. We therefore investigated if the results of a joint analysis of rainfall event responses of near-stream groundwater levels and stream water levels could serve as a viable proxy for hillslope–stream connectivity. The analysis focuses on the extent of response, correlations, lag times and synchronicity. As a first step, a new data analysis scheme was developed, separating the aspects of (a) response timing and (b) extent of water level change. This provides new perspectives on the relationship between groundwater and stream responses. In a second step we investigated if this analysis can give an indication of hillslope–stream connectivity at the catchment scale. Stream water levels and groundwater levels were measured at five different hillslopes over 5 to 6 years. Using a new detection algorithm, we extracted 706 rainfall response events for subsequent analysis. Carrying out this analysis in two different geological regions (schist and marls) allowed us to test the usefulness of the proxy under different hydrological settings while also providing insight into the geologically driven differences in response behaviour. For rainfall events with low initial groundwater level, groundwater level responses often lag behind the stream with respect to the start of rise and the time of peak. This lag disappears at high antecedent groundwater levels. At low groundwater levels the relationship between groundwater and stream water level responses to rainfall are highly variable, while at high groundwater levels, above a certain threshold, this relationship tends to become more uniform. The same threshold was able to predict increased likelihood for high runoff coefficients, indicating a strong increase in connectivity once the groundwater level threshold was surpassed. The joint analysis of shallow near-stream groundwater and stream water levels provided information on the presence or absence and to a certain extent also on the degree of subsurface hillslope–stream connectivity. The underlying threshold processes were interpreted as transmissivity feedback in the marls and fill-and-spill in the schist. The value of these measurements is high; however, time series of several years and a large number of events are necessary to produce representative results. We also find that locally measured thresholds in groundwater levels can provide insight into the connectivity and event response of the corresponding headwater catchments. If the location of the well is chosen wisely, a single time series of shallow groundwater can indicate if the catchment is in a state of high or low connectivity.

The fields of eco-hydrological modelling and extreme flow prediction and management demand detailed information of streamflow intermittency and its corresponding landscape controls. Innovative sensing technology for monitoring of streamflow intermittency in perennial rivers and intermittent reaches improves data availability, but reliable maps of streamflow intermittency are still rare. We used a large dataset of streamflow intermittency observations and a set of spatial predictors to create logistic regression models to predict the probability of streamflow intermittency for a full year as well as wet and dry periods for the entire 247 km2 Attert catchment in Luxembourg. Similar climatic conditions across the catchment permit a direct comparison of the streamflow intermittency among different geological and pedological regions. We used 15 spatial predictors describing land cover, track (road) density, terrain metrics, soil and geological properties. Predictors were included as local-scale information, represented by the local value at the catchment outlet and as integral catchment information calculated as the mean catchment value over all pixels upslope of the catchment outlet. The terrain metrics catchment area and profile curvature were identified in all models as the most important predictors, and the model for the wet period was based solely on these two predictors. However, the model for the dry period additionally comprises soil hydraulic conductivity and bedrock permeability. The annual model with the most complex predictor set contains the predictors of the dry-period model plus the presence of tracks. Classifying the spatially distributed streamflow intermittency probabilities into ephemeral, intermittent and perennial reaches allows the estimation of stream network extent under various conditions. This approach, based on extensive monitoring and statistical modelling, is a first step to provide detailed spatial information for hydrological modelling as well as management practice.

Preferential flow strongly controls water flow and transport in soils. It is ubiquitous but difficult to characterize and predict. This study addresses the occurrence and the evolution of preferential flow during the evolution of landscapes and here specifically during the evolution of hillslopes. We targeted a chronosequence of glacial moraines in the Swiss Alps to investigate how water flow paths evolve along with the soil-forming processes. Dye tracer irrigation experiments with a Brilliant Blue FCF solution (4 g L−1) were conducted on four moraines of different ages (30, 160, 3000, and 10 000 years). At each moraine, three dye tracer experiments were conducted on plots of 1.5 m ×1.0 m. The three plots at each moraine were characterized by different vegetation complexities (low, medium, and high). Each plot was further divided into three equal subplots for the application of three different irrigation amounts (20, 40, and 60 mm) with an average irrigation intensity of 20 mm h−1. The day after the experiment five vertical soil sections were excavated, and the stained flow paths were photographed. Digital image analysis was used to derive average infiltration depths and flow path characteristics such as the volume and surface density of the dye patterns. Based on the volume density, the observed dye patterns were assigned to specific flow type categories. The results show a significant change in the type of preferential flow paths along the chronosequence. The flow types change from a rather homogeneous matrix flow in coarse material with high conductivities and a sparse vegetation cover at the youngest moraine to a heterogeneous infiltration pattern at the medium-age moraines. Heterogeneous matrix and finger flow are dominant at these intermediate age classes. At the oldest moraine only macropore flow via root channels was observed in deeper parts of the soil, in combination with a very high water storage capacity of the organic top layer and low hydraulic conductivity of the deeper soil. In general, we found an increase in water storage with increasing age of the moraines, based on our observations of the reduction in infiltration depth as well as laboratory measurements of porosity. Preferential flow is, however, not only caused by macropores, but especially for the medium-age moraine, it seems to be mainly initiated by soil surface characteristics (vegetation patches and microtopography).

Soil moisture at the catchment scale exhibits a huge spatial variability. This suggests that even a large amount of observation points would not be able to capture soil moisture variability. We present a measure to capture the spatial dissimilarity and its change over time. Statistical dispersion among observation points is related to their distance to describe spatial patterns. We analyzed the temporal evolution and emergence of these patterns and used the mean shift clustering algorithm to identify and analyze clusters. We found that soil moisture observations from the 19.4 km2 Colpach catchment in Luxembourg cluster in two fundamentally different states. On the one hand, we found rainfall-driven data clusters, usually characterized by strong relationships between dispersion and distance. Their spatial extent roughly matches the average hillslope length in the study area of about 500 m. On the other hand, we found clusters covering the vegetation period. In drying and then dry soil conditions there is no particular spatial dependence in soil moisture patterns, and the values are highly similar beyond hillslope scale. By combining uncertainty propagation with information theory, we were able to calculate the information content of spatial similarity with respect to measurement uncertainty (when are patterns different outside of uncertainty margins?). We were able to prove that the spatial information contained in soil moisture observations is highly redundant (differences in spatial patterns over time are within the error margins). Thus, they can be compressed (all cluster members can be substituted by one representative member) to only a fragment of the original data volume without significant information loss. Our most interesting finding is that even a few soil moisture time series bear a considerable amount of information about dynamic changes in soil moisture. We argue that distributed soil moisture sampling reflects an organized catchment state, where soil moisture variability is not random. Thus, only a small amount of observation points is necessary to capture soil moisture dynamics.

Understanding components of the total streamflow is important to assess the ecological functioning of rivers. Binary or two-component separation of streamflow into a quick and a slow (often referred to as baseflow) component are often based on arbitrary choices of separation parameters and also merge different delayed components into one baseflow component and one baseflow index (BFI). As streamflow generation during dry weather often results from drainage of multiple sources, we propose to extend the BFI by a delayed-flow index (DFI) considering the dynamics of multiple delayed contributions to streamflow. The DFI is based on characteristic delay curves (CDCs) where the identification of breakpoint (BP) estimates helps to avoid rather subjective separation parameters and allows for distinguishing four types of delayed streamflow contributions. The methodology is demonstrated using streamflow records from a set of 60 mesoscale catchments in Germany and Switzerland covering a pronounced elevation gradient of roughly 3000 m. We found that the quickflow signal often diminishes earlier than assumed by two-component BFI analyses and distinguished a variety of additional flow contributions with delays shorter than 60 d. For streamflow contributions with delays longer than 60 d, we show that the method can be used to assess catchments' water sustainability during dry spells. Colwell's predictability (PT), a measure of streamflow periodicity and sustainability, was applied to attribute the identified delay patterns to dynamic catchment storage. The smallest dynamic storages were consistently found for catchments between approx. 800 and 1800 m a.s.l. Above an elevation of 1800 m the DFI suggests that seasonal snowpack provides the primary contribution, whereas below 800 m groundwater resources are most likely the major streamflow contributions. Our analysis also indicates that dynamic storage in high alpine catchments might be large and is overall not smaller than in lowland catchments. We conclude that the DFI can be used to assess the range of sources forming catchments' storages and to judge the long-term sustainability of streamflow.

The spatial and temporal controls of preferential flow (PF) during infiltration are still not fully understood. As soil moisture sensor networks allow us to capture infiltration responses in high temporal and spatial resolution, our study is based on a large-scale sensor network with 135 soil moisture profiles distributed across a complex catchment. The experimental design covers three major geological regions (slate, marl, sandstone) and two land covers (forest, grassland) in Luxembourg. We analyzed the responses of up to 353 rainfall events for each of the 135 soil moisture profiles. Non-sequential responses (NSRs) within the soil moisture depth profiles were taken as one indication of bypass flow. For sequential responses maximum porewater velocities (vmax⁡) were determined from the observations and compared with velocity estimates of capillary flow. A measured vmax⁡ higher than the capillary prediction was taken as a further indication of PF. While PF was identified as a common process during infiltration, it was also temporally and spatially highly variable. We found a strong dependence of PF on the initial soil water content and the maximum rainfall intensity. Whereas a high rainfall intensity increased PF (NSR, vmax⁡) as expected, most geologies and land covers showed the highest PF under dry initial conditions. Hence, we identified a strong seasonality of both NSR and vmax⁡ dependent on land cover, revealing a lower occurrence of PF during spring and increased occurrence during summer and early autumn, probably due to water repellency. We observed the highest fraction of NSR in forests on clay-rich soils (slate, marl). vmax⁡ ranged from 6 to 80 640 cm d−1 with a median of 120 cm d−1 across all events and soil moisture profiles. The soils in the marl geology had the highest flow velocities, independent of land cover, especially between 30 and 50 cm depth, where the clay content increased. This demonstrates the danger of treating especially clay soils in the vadose zone as a low-conductive substrate, as the development of soil structure can dominate over the matrix property of the texture alone. This confirms that clay content and land cover strongly influence infiltration and reinforce PF, but seasonal dynamics and flow initiation also have an important impact on PF.

The water isotopic composition of throughfall is affected by complex diffusive exchange with ambient water vapour, evaporative enrichment of heavy isotopes, and mixing processes in the tree canopy. All interception processes occur simultaneously in space and time, generating a complex pattern of throughfall depth and water isotopic composition. This pattern ultimately cascades through the entire hydrologic system and is therefore crucial for isotope studies in catchment hydrology, where recharge areas are often forested, while reference meteorological stations are generally in the open. For the quasi real-time observation of the water isotopic composition (δ18O and δ2H) of both gross precipitation and throughfall, we developed an approach combining a membrane contactor (Membrana) with a laser-based Cavity Ring-Down Spectrometer (CRDS, Picarro), obtaining isotope readings every 2 s. A setup with two CRDS instruments in parallel analysing gross precipitation and throughfall simultaneously was used for the continuous observation of the temporal effect of interception processes on the stable isotopes of water. All devices were kept small to minimize dead volume with time lags of only 4 min for water from the rainfall collectors to the isotope analysers to increase the temporal resolution of isotope observations. Complementarily, meteorological variables were recorded at high temporal resolution at the same location. The achieved evolution from discrete liquid or event-based bulk samples to continuous measurements allows for direct comparison of water stable isotope data with common meteorological measurements. Future improvements of the spatial representativeness will make our approach an even more powerful tool towards detailed insight into the dynamic processes contributing to interception during rainfall events.

The young water fraction Fyw, defined as the proportion of catchment outflow younger than approximately 2–3 months, can be estimated directly from the amplitudes of seasonal cycles of stable water isotopes in precipitation and streamflow. Thus, Fyw may be a useful metric in catchment inter-comparison studies that investigate landscape and hydro-climatic controls on streamflow generation. Here, we explore how Fyw varies with catchment characteristics and climatic forcing, using an extensive isotope data set from 22 small- to medium-sized (0.7–351 km2) Swiss catchments. We find that flow-weighting the tracer concentrations in streamwater resulted in roughly 26 % larger young water fractions compared to the corresponding unweighted values, reflecting the fact that young water fractions tend to be larger when catchments are wet and discharge is correspondingly higher. However, flow-weighted and unweighted young water fractions are strongly correlated with each other among the catchments. They also correlate with terrain, soil, and land-use indices, as well as with mean precipitation and measures of hydrologic response. Within individual catchments, young water fractions increase with discharge, indicating an increase in the proportional contribution of faster flow paths at higher flows. We present a new method to quantify the discharge sensitivity of Fyw, which we estimate as the linear slope of the relationship between the young water fraction and flow. Among the 22 catchments, discharge sensitivities of Fyw are highly variable and only weakly correlated with Fyw itself, implying that these two measures reflect catchment behaviour differently. Based on strong correlations between the discharge sensitivity of Fyw and several catchment characteristics, we suggest that low discharge sensitivities imply greater persistence in the proportions of fast and slow runoff flow paths as catchment wetness changes. In contrast, high discharge sensitivities imply the activation of different dominant flow paths during precipitation events, such as when subsurface water tables rise into more permeable layers and/or the river network expands further into the landscape.

For more than two decades, research groups in hydrology, ecology, soil science, and biogeochemistry have performed cryogenic water extractions (CWEs) for the analysis of δ2H and δ18O of soil water. Recent studies have shown that extraction conditions (time, temperature, and vacuum) along with physicochemical soil properties may affect extracted soil water isotope composition. Here we present results from the first worldwide round robin laboratory intercomparison. We test the null hypothesis that, with identical soils, standards, extraction protocols, and isotope analyses, cryogenic extractions across all laboratories are identical. Two standard soils with different physicochemical characteristics along with deionized (DI) reference water of known isotopic composition were shipped to 16 participating laboratories. Participants oven-dried and rewetted the soils to 8 and 20 % gravimetric water content (WC), using the deionized reference water. One batch of soil samples was extracted via predefined extraction conditions (time, temperature, and vacuum) identical to all laboratories; the second batch was extracted via conditions considered routine in the respective laboratory. All extracted water samples were analyzed for δ18O and δ2H by the lead laboratory (Global Institute for Water Security, GIWS, Saskatoon, Canada) using both a laser and an isotope ratio mass spectrometer (OA-ICOS and IRMS, respectively). We rejected the null hypothesis. Our results showed large differences in retrieved isotopic signatures among participating laboratories linked to soil type and soil water content with mean differences compared to the reference water ranging from +18.1 to −108.4 ‰ for δ2H and +11.8 to −14.9 ‰ for δ18O across all laboratories. In addition, differences were observed between OA-ICOS and IRMS isotope data. These were related to spectral interferences during OA-ICOS analysis that are especially problematic for the clayey loam soils used. While the types of cryogenic extraction lab construction varied from manifold systems to single chambers, no clear trends between system construction, applied extraction conditions, and extraction results were found. Rather, observed differences in the isotope data were influenced by interactions between multiple factors (soil type and properties, soil water content, system setup, extraction efficiency, extraction system leaks, and each lab's internal accuracy). Our results question the usefulness of cryogenic extraction as a standard for water extraction since results are not comparable across laboratories. This suggests that defining any sort of standard extraction procedure applicable across laboratories is challenging. Laboratories might have to establish calibration functions for their specific extraction system for each natural soil type, individually.

Glaciers play an important role in high-mountain hydrology. While changing glacier areas are considered of highest importance for the understanding of future changes in runoff, glaciers are often only poorly represented in hydrological models. Most importantly, the direct coupling between the simulated glacier mass balances and changing glacier areas needs feasible solutions. The use of a complex glacier model is often not possible due to data and computational limitations. The Δh parameterization is a simple approach to consider the spatial variation of glacier thickness and area changes. Here, we describe a conceptual implementation of the Δh parameterization in the semi-distributed hydrological model HBV-light, which also allows for the representation of glacier advance phases and for comparison between the different versions of the implementation. The coupled glacio-hydrological simulation approach, which could also be implemented in many other semi-distributed hydrological models, is illustrated based on an example application.

Transpiration is a key process in the hydrological cycle, and a sound understanding and quantification of transpiration and its spatial variability is essential for management decisions as well as for improving the parameterisation and evaluation of hydrological and soil–vegetation–atmosphere transfer models. For individual trees, transpiration is commonly estimated by measuring sap flow. Besides evaporative demand and water availability, tree-specific characteristics such as species, size or social status control sap flow amounts of individual trees. Within forest stands, properties such as species composition, basal area or stand density additionally affect sap flow, for example via competition mechanisms. Finally, sap flow patterns might also be influenced by landscape-scale characteristics such as geology and soils, slope position or aspect because they affect water and energy availability; however, little is known about the dynamic interplay of these controls.We studied the relative importance of various tree-, stand- and site-specific characteristics with multiple linear regression models to explain the variability of sap velocity measurements in 61 beech and oak trees, located at 24 sites across a 290 km2 catchment in Luxembourg. For each of 132 consecutive days of the growing season of 2014 we modelled the daily sap velocity and derived sap flow patterns of these 61 trees, and we determined the importance of the different controls.Results indicate that a combination of mainly tree- and site-specific factors controls sap velocity patterns in the landscape, namely tree species, tree diameter, geology and aspect. For sap flow we included only the stand- and site-specific predictors in the models to ensure variable independence. Of those, geology and aspect were most important. Compared to these predictors, spatial variability of atmospheric demand and soil moisture explains only a small fraction of the variability in the daily datasets. However, the temporal dynamics of the explanatory power of the tree-specific characteristics, especially species, are correlated to the temporal dynamics of potential evaporation. We conclude that transpiration estimates on the landscape scale would benefit from not only consideration of hydro-meteorological drivers, but also tree, stand and site characteristics in order to improve the spatial and temporal representation of transpiration for hydrological and soil–vegetation–atmosphere transfer models.

Few studies have quantified the differences between celerity and velocity of hillslope water flow and explained the processes that control these differences. Here, we asses these differences by combining a 24-day hillslope sprinkling experiment with a spatially explicit hydrologic model analysis. We focused our work on Watershed 10 at the H. J. Andrews Experimental Forest in western Oregon. Celerities estimated from wetting front arrival times were generally much faster than average vertical velocities of δ2H. In the model analysis, this was consistent with an identifiable effective porosity (fraction of total porosity available for mass transfer) parameter, indicating that subsurface mixing was controlled by an immobile soil fraction, resulting in the attenuation of the δ2H input signal in lateral subsurface flow. In addition to the immobile soil fraction, exfiltrating deep groundwater that mixed with lateral subsurface flow captured at the experimental hillslope trench caused further reduction in the δ2H input signal. Finally, our results suggest that soil depth variability played a significant role in the celerity–velocity responses. Deeper upslope soils damped the δ2H input signal, while a shallow soil near the trench controlled the δ2H peak in lateral subsurface flow response. Simulated exit time and residence time distributions with our hillslope hydrologic model showed that water captured at the trench did not represent the entire modeled hillslope domain; the exit time distribution for lateral subsurface flow captured at the trench showed more early time weighting.

The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al.(2017).

We combine ecohydrological observations of sap flow and soil moisture with thermodynamically constrained estimates of atmospheric evaporative demand to infer the dominant controls of forest transpiration in complex terrain. We hypothesize that daily variations in transpiration are dominated by variations in atmospheric demand, while site-specific controls, including limiting soil moisture, act on longer timescales. We test these hypotheses with data of a measurement setup consisting of five sites along a valley cross section in Luxembourg. Both hillslopes are covered by forest dominated by European beech (Fagus sylvatica L.). Two independent measurements are used to estimate stand transpiration: (i) sap flow and (ii) diurnal variations in soil moisture, which were used to estimate the daily root water uptake. Atmospheric evaporative demand is estimated through thermodynamically constrained evaporation, which only requires absorbed solar radiation and temperature as input data without any empirical parameters. Both transpiration estimates are strongly correlated to atmospheric demand at the daily timescale. We find that neither vapor pressure deficit nor wind speed add to the explained variance, supporting the idea that they are dependent variables on land–atmosphere exchange and the surface energy budget. Estimated stand transpiration was in a similar range at the north-facing and the south-facing hillslopes despite the different aspect and the largely different stand composition. We identified an inverse relationship between sap flux density and the site-average sapwood area per tree as estimated by the site forest inventories. This suggests that tree hydraulic adaptation can compensate for heterogeneous conditions. However, during dry summer periods differences in topographic factors and stand structure can cause spatially variable transpiration rates. We conclude that absorption of solar radiation at the surface forms a dominant control for turbulent heat and mass exchange and that vegetation across the hillslope adjusts to this constraint at the tree and stand level. These findings should help to improve the description of land-surface–atmosphere exchange at regional scales.

Climate change is expected to impact the water cycle and severely affect precipitation patterns across central Europe and in other parts of the world, leading to more frequent and severe droughts. Usually when projecting drought impacts on hydrological systems, it is assumed that system properties, like soil properties, remain stable and will not be affected by drought events. To study if this assumption is appropriate, we address the effects of drought on the infiltration behavior of forest soils using dye tracer experiments on six sites in three regions across Germany, which were forced into drought conditions. The sites cover clayey-, loamy- and sandy-textured soils. In each region, we compared a deciduous and a coniferous forest stand to address differences between the main tree species. The results of the dye tracer experiments show clear evidence for changes in infiltration behavior at the sites. The infiltration changed at the clayey plots from regular and homogeneous flow to fast preferential flow. Similar behavior was observed at the loamy plots, where large areas in the upper layers remained dry, displaying signs of strong water repellency. This was confirmed by water drop penetration time (WDPT) tests, which revealed, in all except one plot, moderate to severe water repellency. Water repellency was also accountable for the change of regular infiltration to fingered flow in the sandy soils. The results of this study suggest that the drought history or, more generally, the climatic conditions of a soil in the past are more important than the actual antecedent soil moisture status regarding hydrophobicity and infiltration behavior; furthermore, drought effects on infiltration need to be considered in hydrological models to obtain realistic predictions concerning water quality and quantity in runoff and groundwater recharge.

Several controls are known to affect water quality of stream networks during flow recession periods, such as solute leaching processes, surface water–groundwater interactions as well as biogeochemical in-stream turnover processes. Throughout the stream network, combinations of specific water and solute export rates and local in-stream conditions overlay the biogeochemical signals from upstream sections. Therefore, upstream sections can be considered functional units which could be distinguished and ordered regarding their relative contribution to nutrient dynamics at the catchment outlet. Based on snapshot sampling of flow and nitrate concentrations along the stream in an agricultural headwater during the summer flow recession period, we determined spatial and temporal patterns of water quality for the whole stream. A data-driven, in-stream-mixing-and-removal model was developed and applied for analysing the spatio-temporal in-stream retention processes and their effect on the spatio-temporal fluxes of nitrate from subcatchments. Thereby, we have been able to distinguish quantitatively between nitrate sinks, sources per stream reaches, and subcatchments, and thus we could disentangle the overlay of nitrate sink and source signals. For nitrate sources, we determined their permanent and temporal impact on stream water quality and for nitrate sinks, we found increasing nitrate removal efficiencies from upstream to downstream. Our results highlight the importance of distinct nitrate source locations within the watershed for in-stream concentrations and in-stream removal processes, respectively. Thus, our findings contribute to the development of a more dynamic perception of water quality in streams and rivers concerning ecological and sustainable water resource management.

Determining the soil hydraulic properties is a prerequisite to physically model transient water flow and solute transport in the vadose zone. Estimating these properties by inverse modelling techniques has become more common within the last 2 decades. While these inverse approaches usually fit simulations to hydrometric data, we expanded the methodology by using independent information about the stable isotope composition of the soil pore water depth profile as a single or additional optimization target. To demonstrate the potential and limits of this approach, we compared the results of three inverse modelling strategies where the fitting targets were (a) pore water isotope concentrations, (b) a combination of pore water isotope concentrations and soil moisture time series, and (c) a two-step approach using first soil moisture data to determine water flow parameters and then the pore water stable isotope concentrations to estimate the solute transport parameters. The analyses were conducted at three study sites with different soil properties and vegetation. The transient unsaturated water flow was simulated by solving the Richards equation numerically with the finite-element code of HYDRUS-1D. The transport of deuterium was simulated with the advection-dispersion equation, and a modified version of HYDRUS was used, allowing deuterium loss during evaporation. The Mualem–van Genuchten and the longitudinal dispersivity parameters were determined for two major soil horizons at each site. The results show that approach (a), using only the pore water isotope content, cannot substitute hydrometric information to derive parameter sets that reflect the observed soil moisture dynamics but gives comparable results when the parameter space is constrained by pedotransfer functions. Approaches (b) and (c), using both the isotope profiles and the soil moisture time series, resulted in good simulation results with regard to the Kling–Gupta efficiency and good parameter identifiability. However, approach (b) has the advantage that it considers the isotope data not only for the solute transport parameters but also for water flow and root water uptake, and thus increases parameter realism. Approaches (b) and (c) both outcompeted simulations run with parameters derived from pedotransfer functions, which did not result in an acceptable representation of the soil moisture dynamics and pore water stable isotope composition. Overall, parameters based on this new approach that includes isotope data lead to similar model performances regarding the water balance and soil moisture dynamics and better parameter identifiability than the conventional inverse model approaches limited to hydrometric fitting targets. If only data from isotope profiles in combination with textural information is available, the results are still satisfactory. This method has the additional advantage that it will not only allow us to estimate water balance and response times but also site-specific time variant transit times or solute breakthrough within the soil profile.

Meteorological droughts like those in summer 2003 or spring 2011 in Europe are expected to become more frequent in the future. Although the spatial extent of these drought events was large, not all regions were affected in the same way. Many catchments reacted strongly to the meteorological droughts showing low levels of streamflow and groundwater, while others hardly reacted. Also, the extent of the hydrological drought for specific catchments was different between these two historical events due to different initial conditions and drought propagation processes. This leads to the important question of how to detect and quantify the sensitivity of a catchment to meteorological droughts. To assess this question we designed hydrological model experiments using a conceptual rainfall-runoff model. Two drought scenarios were constructed by selecting precipitation and temperature observations based on certain criteria: one scenario was a modest but constant progression of drying based on sorting the years of observations according to annual precipitation amounts. The other scenario was a more extreme progression of drying based on selecting months from different years, forming a year with the wettest months through to a year with the driest months. Both scenarios retained the observed intra-annual seasonality for the region. We evaluated the sensitivity of 24 Swiss catchments to these scenarios by analyzing the simulated discharge time series and modeled storage. Mean catchment elevation, slope and area were the main controls on the sensitivity of catchment discharge to precipitation. Generally, catchments at higher elevation and with steeper slopes appeared less sensitive to meteorological droughts than catchments at lower elevations with less steep slopes.

The transit time of water is a fundamental property of catchments, revealing information about the flow pathways, source of water and storage in a single integrated measure. While several studies have investigated the relationship between catchment topography and transit times, few studies expanded the analysis to a wide range of catchment properties and assessed the influence of the selected transfer function (TF) model. We used stable water isotopes from mostly baseflow samples with lumped convolution models of time-invariant TFs to estimate the transit time distributions of 24 meso-scale catchments covering different geomorphic and geologic regions in Switzerland. The sparse network of 13 precipitation isotope sampling sites required the development of a new spatial interpolation method for the monthly isotopic composition of precipitation. A point-energy-balance based snow model was adapted to account for the seasonal water isotope storage in snow dominated catchments. Transit time distributions were estimated with three established TFs (exponential, gamma distribution and two parallel linear reservoirs). While the exponential TF proved to be less suitable to simulate the isotopic signal in most of the catchments, the gamma distribution and the two parallel linear reservoirs transfer function reached similarly good model fits to the fortnightly observed isotopic compositions in discharge, although in many catchments the transit time distributions implied by equally well fitted models differed markedly from each other and in extreme cases, the resulting mean transit time (MTT) differed by orders of magnitude. A more thorough comparison showed that equally suited models corresponded to agreeing values of cumulated transit time distributions only between 3 and 6 months. The short-term (< 30 days) component of the transit time distributions did not play a role because of the limited temporal resolution of the available input data. The long-term component (> 3 years) could hardly be assessed by means of stable water isotopes, resulting in ambiguous MTT and hence questioning the relevance of an MTT determined with stable isotopes. Finally we investigated the relation between MTT estimates based on the three different TF types as well as other transit time properties and a range of topographical catchment characteristics. Depending on the selected transfer model, we found a weak correlation between transit time properties and the ratio between flow path length over the flow gradient, drainage density and the mean discharge. The catchment storage derived from MTTs and mean discharge did not show a clear relation to any catchment properties, indicating that in many studies the mean annual discharge may bias the MTT estimates.

According to Dooge (1986) intermediate-scale catchments are systems of organized complexity, being too organized and yet too small to be characterized on a statistical/conceptual basis, but too large and too heterogeneous to be characterized in a deterministic manner. A key requirement for building structurally adequate models precisely for this intermediate scale is a better understanding of how different forms of spatial organization affect storage and release of water and energy. Here, we propose that a combination of the concept of hydrological response units (HRUs) and thermodynamics offers several helpful and partly novel perspectives for gaining this improved understanding. Our key idea is to define functional similarity based on similarity of the terrestrial controls of gradients and resistance terms controlling the land surface energy balance, rainfall runoff transformation, and groundwater storage and release. This might imply that functional similarity with respect to these specific forms of water release emerges at different scales, namely the small field scale, the hillslope, and the catchment scale. We thus propose three different types of "functional units" – specialized HRUs, so to speak – which behave similarly with respect to one specific form of water release and with a characteristic extent equal to one of those three scale levels. We furthermore discuss an experimental strategy based on exemplary learning and replicate experiments to identify and delineate these functional units, and as a promising strategy for characterizing the interplay and organization of water and energy fluxes across scales. We believe the thermodynamic perspective to be well suited to unmask equifinality as inherent in the equations governing water, momentum, and energy fluxes: this is because several combinations of gradients and resistance terms yield the same mass or energy flux and the terrestrial controls of gradients and resistance terms are largely independent. We propose that structurally adequate models at this scale should consequently disentangle driving gradients and resistance terms, because this optionally allows equifinality to be partly reduced by including available observations, e.g., on driving gradients. Most importantly, the thermodynamic perspective yields an energy-centered perspective on rainfall-runoff transformation and evapotranspiration, including fundamental limits for energy fluxes associated with these processes. This might additionally reduce equifinality and opens up opportunities for testing thermodynamic optimality principles within independent predictions of rainfall-runoff or land surface energy exchange. This is pivotal to finding out whether or not spatial organization in catchments is in accordance with a fundamental organizing principle.

Characterizing hydrological processes within tropical wetlands is challenging due to their remoteness, complexity and heterogeneity. In particular, estimates of evaporative water loss are inherently uncertain. In view of the large influence on the local and regional climate, the quantification of evaporation is essential for the determination of the water balance of permanent and intermittent water bodies. Data for tropical wetlands are scarce where their remoteness impedes direct evaporation measurements. Seasonal inundation dynamics affect evaporation processes in tropical wetlands, which can be analysed in two stages: the first stage during the wet season and the second stage during the dry season. As yet no adequate method exists for determining second-stage evaporation in a data-scarce environment that additionally allows for a transfer of simulated actual evaporation (AET) to other locations. Our study aimed at developing a process-based model to simulate first- and second-stage evaporation in tropical wetlands. We selected a set of empirical potential evaporation (PET) models of varying complexity, each based on different assumptions and available data sets, and evaluated the models with pan evaporation observations in the Pantanal of South America, one of the largest tropical wetlands in the world. We used high-resolution measurements of surface and groundwater levels at different locations to determine the water available for evaporation. AET was derived by constraining simulated PET based on available water. The model of best fit was applied to different types of water bodies with varying hydroperiods to capture first- and second-stage evaporation across a range of wetland types. With our new model we could quantify evaporative water loss in the dry and the wet season for different locations in the Pantanal. This new spatially explicit approach represents an improvement in our understanding of the role of evaporation in the water balance of the Pantanal. We recommend the application of this model in other remote tropical wetlands, since only a minimum of input data is necessary.

In January 2011 a rain-on-snow (RoS) event caused floods in the major river basins in central Europe, i.e. the Rhine, Danube, Weser, Elbe, Oder, and Ems. This event prompted the questions of how to define a RoS event and whether those events have become more frequent. Based on the flood of January 2011 and on other known events of the past, threshold values for potentially flood-generating RoS events were determined. Consequently events with rainfall of at least 3 mm on a snowpack of at least 10 mm snow water equivalent (SWE) and for which the sum of rainfall and snowmelt contains a minimum of 20% snowmelt were analysed. RoS events were estimated for the time period 1950–2011 and for the entire study area based on a temperature index snow model driven with a European-scale gridded data set of daily climate (E-OBS data). Frequencies and magnitudes of the modelled events differ depending on the elevation range. When distinguishing alpine, upland, and lowland basins, we found that upland basins are most influenced by RoS events. Overall, the frequency of rainfall increased during winter, while the frequency of snowfall decreased during spring. A decrease in the frequency of RoS events from April to May has been observed in all upland basins since 1990. In contrast, the results suggest an increasing trend in the magnitude and frequency of RoS days in January and February for most of the lowland and upland basins. These results suggest that the flood hazard from RoS events in the early winter season has increased in the medium-elevation mountain ranges of central Europe, especially in the Rhine, Weser, and Elbe river basins.

Stable isotope signatures provide an integral fingerprint of origin, flow paths, transport processes, and residence times of water in the environment. However, the full potential of stable isotopes to quantitatively characterize subsurface water dynamics is yet unfolded due to the difficulty in obtaining extensive, detailed, and repeated measurements of pore water in the unsaturated and saturated zone. This paper presents a functional and cost-efficient system for non-destructive continual in situ monitoring of pore water stable isotope signatures with high resolution. Automatic controllable valve arrays are used to continuously extract diluted water vapor in soil air via a branching network of small microporous probes into a commercial laser-based isotope analyzer. Normalized liquid-phase isotope signatures are then obtained based on a specific on-site calibration approach along with basic corrections for instrument bias and temperature dependent isotopic fractionation. The system was applied to sample depth profiles on three experimental plots with varied vegetation cover in southwest Germany. Two methods (i.e., based on advective versus diffusive vapor extraction) and two modes of sampling (i.e., using multiple permanently installed probes versus a single repeatedly inserted probe) were tested and compared. The results show that the isotope distribution along natural profiles could be resolved with sufficiently high accuracy and precision at sampling intervals of less than four minutes. The presented in situ approaches may thereby be used interchangeably with each other and with concurrent laboratory-based direct equilibration measurements of destructively collected samples. It is thus found that the introduced sampling techniques provide powerful tools towards a detailed quantitative understanding of dynamic and heterogeneous shallow subsurface and vadose zone processes.

The role of spatial variability in water inputs on runoff dynamics has generally not received as much research attention as topography and soils; however, the influence of topography and forest cover on snow surface energy exchanges can result in asynchronous snowmelt throughout a catchment, complicating the space–time patterns of runoff generation. This study investigates temporal variation in the relative importance of spatial controls on the occurrence, duration, and timing of shallow groundwater response, utilizing a highly distributed monitoring network in a snowmelt-dominated montane catchment in western Canada. The study findings indicate that deep-soil hydraulic conductivity is a first-order control on the spatial distribution of sites that generate shallow groundwater response versus sites that experience only deep percolation. Upslope contributing area and slope gradient are first-order controls on the duration of groundwater response during peak-flow, recession-flow, and low-flow periods. Shallow runoff response areas expand and contract throughout these periods and follow the general spatial patterns of topographic convergence. However, spatial controls on the timing, intensity, and quantity of snowmelt and controls on vertical versus lateral flux partitioning in the soil overwhelm the influence of topographic convergence on runoff patterns during early spring freshet periods. The study findings suggest that various topographic indices and topography-based rainfall runoff models would not likely be good predictors of runoff patterns in snowmelt-dominated montane catchments during early phases of the spring freshet, but would increase in importance as the freshet and post-freshet periods proceed.

Pesticides applied onto agricultural fields are frequently found in adjacent rivers. To what extent and along which pathways they are transported is influenced by intrinsic pesticide properties such as sorption and degradation. In the environment, incomplete degradation of pesticides leads to the formation of transformation products (TPs), which may differ from the parent compounds regarding their intrinsic fate characteristics. Thus, the export processes of TPs in catchments and streams may also be different. In order to test this hypothesis, we extended a distributed hydrological model by the fate and behaviour of pesticides and transformation products and applied it to a small, well-monitored headwater catchment in Switzerland. The successful model evaluation of three pesticides and their TPs at three sampling locations in the catchment enabled us to estimate the quantity of contributing processes for pollutant export. Since all TPs were more mobile than their parent compounds (PCs), they exhibited larger fractions of export via subsurface pathways. However, besides freshly applied pesticides, subsurface export was found to be influenced by residues of former applications. Export along preferential flow pathways was less dependent on substance fate characteristics than soil matrix export, but total soil water flow to tile drains increased more due to preferential flow for stronger sorbing substances. Our results indicate that runoff generation by matrix flow to tile drains gained importance towards the end of the modelling period whereas the contributions from fast surface runoff and preferential flow decreased. Accordingly, TPs were to a large extent exported under different hydrological conditions than their PCs, due to their delayed formation and longer half-lives. Thus, not only their different intrinsic characteristics but also their delayed formation could be responsible for the fact that TPs generally took different pathways than their PCs. We suggest that these results should be considered in risk assessment for the export of agricultural chemicals to adjacent rivers and that models should be extended to include both PCs and TPs.

More than 30% of Europe's land surface is made up of karst exposures. In some countries, water from karst aquifers constitutes almost half of the drinking water supply. Hydrological simulation models can predict the large-scale impact of future environmental change on hydrological variables. However, the information needed to obtain model parameters is not available everywhere and regionalisation methods have to be applied. The responsive behaviour of hydrological systems can be quantified by individual metrics, so-called system signatures. This study explores their value for distinguishing the dominant processes and properties of five different karst systems in Europe and the Middle East. By defining ten system signatures derived from hydrodynamic and hydrochemical observations, a process-based karst model is applied to the five karst systems. In a stepwise model evaluation strategy, optimum parameters and their sensitivity are identified using automatic calibration and global variance-based sensitivity analysis. System signatures and sensitive parameters serve as proxies for dominant processes, and optimised parameters are used to determine system properties. By sensitivity analysis, the set of system signatures was able to distinguish the karst systems from one another by providing separate information about dominant soil, epikarst, and fast and slow groundwater flow processes. Comparing sensitive parameters to the system signatures revealed that annual discharge can serve as a proxy for the recharge area, that the slopes of the high flow parts of the flow duration curves correlate with the fast flow storage constant, and that the dampening of the isotopic signal of the rain as well as the medium flow parts of the flow duration curves have a non-linear relation to the distribution of groundwater storage constants that represent the variability of groundwater flow dynamics. Our approach enabled us to identify dominant processes of the different systems and provided directions for future large-scale simulation of karst areas to predict the impact of future change on karst water resources.

The aim of this study was to analyse and predict the mean stable water isotopic composition of the snow cover at specific geographic locations and altitudes. In addition, the dependence of the isotopic composition of the entire snow cover on altitude was analysed. Snow in four Swiss catchments was sampled at the end of the accumulation period in April 2010 and a second time during snowmelt in May 2010 and analysed for stable isotope composition of 2H and 18O. The sampling was conducted at both south-facing and north-facing slopes at elevation differences of 100 m, for a total altitude difference of approximately 1000 m. The observed variability of isotopic composition of the snow cover was analysed with stepwise multiple linear regression models. The analysis indicated that there is only a limited altitude effect on the isotopic composition when considering all samples. This is due to the high variability of the isotopic composition of the precipitation during the winter months and, in particular in the case of south-facing slopes, an enrichment of heavy isotopes due to intermittent melting processes. This enrichment effect could clearly be observed in the samples which were taken later in the year. A small altitudinal gradient of the isotopic composition could only be observed at some north-facing slopes. However, the dependence of snow depth and the day of the year were significant predictor variables in all models. This study indicates the necessity to further study the variability of water isotopes in the snow cover to increase prediction for isotopic composition of snowmelt and hence increase model performance of residence time models for alpine areas in order to better understand the accumulation processes and the sources of water in the snow cover of high mountains.

A network of 45 spatially distributed time-lapse cameras was used to carry out a continuous observation of snow processes and snow cover properties throughout three mid-latitude medium elevation mountain catchments in hourly intervals during the winter seasons of 2010/2011 and 2011/2012. A simple technical modification was conducted to enable the deployment of the standard digital cameras in any location. Image analysis software was applied to extract information about snow depth, surface albedo and canopy interception from the digital images. Furthermore, the distributed design of the camera network made it possible to identify the elevation of the snow rain interface for any precipitation event which is very helpful for the interpretation of winter flooding events resulting from snow melt. Exemplary data for all these analyses is presented to show the potential fields of application of this innovative approach. Study results prove that the application of digital time-lapse photography is an appropriate technique to observe the spatial distribution and temporal evolution of seasonal snow covers in a mountainous environment.

Streamflow recession has been investigated by a variety of methods, often involving the fit of a model to empirical recession plots to parameterize a non-linear storage–outflow relationship based on the dQ/dt−Q method. Such recession analysis methods (RAMs) are used to estimate hydraulic conductivity, storage capacity, or aquifer thickness and to model streamflow recession curves for regionalization and prediction at the catchment scale. Numerous RAMs have been published, but little is known about how comparably the resulting recession models distinguish characteristic catchment behavior. In this study we combined three established recession extraction methods with three different parameter-fitting methods to the power-law storage–outflow model to compare the range of recession characteristics that result from the application of these different RAMs. Resulting recession characteristics including recession time and corresponding storage depletion were evaluated for 20 meso-scale catchments in Germany. We found plausible ranges for model parameterization; however, calculated recession characteristics varied over two orders of magnitude. While recession characteristics of the 20 catchments derived with the different methods correlate strongly, particularly for the RAMs that use the same extraction method, not all rank the catchments consistently, and the differences among some of the methods are larger than among the catchments. To elucidate this variability we discuss the ambiguous roles of recession extraction procedures and the parameterization of the storage–outflow model and the limitations of the presented recession plots. The results suggest strong limitations to the comparability of recession characteristics derived with different methods, not only in the model parameters but also in the relative characterization of different catchments. A multiple-methods approach to investigating streamflow recession characteristics should be considered for applications whenever possible.

Hillslope hydrological dynamics, particularly subsurface flow (SSF), are highly variable and complex. A profound understanding of factors controlling this variability is needed. Therefore we investigated the relationship between variability of shallow water table dynamics and various hillslope characteristics. We ask whether measurable hillslope properties explain patterns of subsurface flow variability. To approach this question, shallow water table dynamics of three adjacent large-scale hillslopes were monitored with high spatial and temporal resolution over 18 months. The hillslopes are similar in terms of topography and parent material, but different in vegetation cover (grassland, coniferous forest, and mixed forest). We expect vegetation to be an important driver of water table dynamics at our study site, especially given the minor differences in topography. Various hillslope properties were determined in the field and via GIS analysis: common topography descriptors, well depth, soil properties via slug tests, and several vegetation parameters. Response variables characterizing the water table response per well were calculated for different temporal scales (entire time series, seasonal scale, event scale). Partial correlation analysis and a Random Forest machine learning approach were carried out to assess the explainability of SSF variability by measurable hillslope characteristics. We found a complex interplay of predictors, yet soil properties and topography showed the highest single explanatory power. Surprisingly, vegetation characteristics played a minor role. Solely throughfall and canopy cover exerted a slightly stronger control, especially in summer. Most importantly, the examined hillslope characteristics explained only a small proportion of the observed SSF variability. Consequently there must be additional important drivers not represented by current measurement techniques of the hillslope configuration (e.g. bedrock properties, preferential pathways). We also found interesting differences in explainability of SSF variability among temporal scales and between both forested hillslopes and the grassland hillslope.

Protection from hydrological extremes and the sustainable supply of hydrological services in the presence of changing climate and lifestyles as well as rocketing population pressure in many parts of the world are the defining societal challenges for hydrology in the 21st century. A review of the existing literature shows that these challenges and their educational consequences for hydrology were foreseeable and were even predicted by some. However, surveys of the current educational basis for hydrology also clearly demonstrate that hydrology education is not yet ready to prepare students to deal with these challenges. We present our own vision of the necessary evolution of hydrology education, which we implemented in the Modular Curriculum for Hydrologic Advancement (MOCHA). The MOCHA project is directly aimed at developing a community-driven basis for hydrology education. In this paper we combine literature review, community survey, discussion and assessment to provide a holistic baseline for the future of hydrology education. The ultimate objective of our educational initiative is to enable educators to train a new generation of "renaissance hydrologists," who can master the holistic nature of our field and of the problems we encounter.

In karst systems, near-surface dissolution of carbonate rock results in a high spatial and temporal variability of groundwater recharge. To adequately represent the dominating recharge processes in hydrological models is still a challenge, especially in data scarce regions. In this study, we developed a recharge model that is based on a conceptual model of the epikarst. It represents epikarst heterogeneity as a set of system property distributions to produce not only a single recharge time series, but a variety of time series representing the spatial recharge variability. We tested the new model with a unique set of spatially distributed flow and tracer observations in a karstic cave at Mt. Carmel, Israel. We transformed the spatial variability into statistical variables and apply an iterative calibration strategy in which more and more data was added to the calibration. Thereby, we could show that the model is only able to produce realistic results when the information about the spatial variability of the observations was included into the model calibration. We could also show that tracer information improves the model performance if data about the spatial variability is not included.

Predicting and understanding subsurface flowpaths is still a crucial issue in hydrological research. We present an experimental approach to reveal present and past subsurface flowpaths of water in the unsaturated and saturated zone. Two hillslopes in a humid mountainous catchment have been investigated. The H2O(liquid) – H2O(vapor) equilibration laser spectroscopy method was used to obtain high resolution δ2H vertical depth profiles of pore water at various points along two fall lines of a pasture hillslope in the southern Black Forest, Germany. The Porewater-based Stable Isotope Profile (PSIP) approach was developed to use the integrated information of several vertical depth profiles of deuterium along transects at the hillslope. Different shapes of depth profiles were observed in relation to hillslope position. The statistical variability (inter-quartile range and standard deviation) of each profile was used to characterize different types of depth profiles. The profiles upslope or with a weak affinity for saturation as indicated by a low topographic wetness index preserve the isotopic input signal by precipitation with a distinct seasonal variability. These observations indicate mainly vertical movement of soil water in the upper part of the hillslope before sampling. The profiles downslope or at locations with a strong affinity for saturation do not show a similar seasonal isotopic signal. The input signal is erased in the foothills and a large proportion of pore water samples are close to the isotopic values of δ2H in streamwater during base flow conditions indicating the importance of the groundwater component in the catchment. Near the stream indications for efficient mixing of water from lateral subsurface flow paths with vertical percolation are found.

The hydrology of ecosystem succession gives rise to new challenges for the analysis and modelling of water balance components. Recent large-scale alterations of forest cover across the globe suggest that a significant portion of new biophysical environments will influence the long-term dynamics and limits of water fluxes compared to pre-succession conditions. This study assesses the estimation of summer evapotranspiration along three FLUXNET sites at Campbell River, British Columbia, Canada using a data-driven soil water balance model validated by Eddy Covariance measurements. It explores the sensitivity of the model to different forest succession states, a wide range of computational time steps, rooting depths, and canopy interception capacity values. Uncertainty in the measured EC fluxes resulting in an energy imbalance was consistent with previous studies and does not affect the validation of the model. The agreement between observations and model estimates proves that the usefulness of the method to predict summer AET over mid- and long-term periods is independent of stand age. However, an optimal combination of the parameters rooting depth, time step and interception capacity threshold is needed to avoid an underestimation of AET as seen in past studies. The study suggests that summer AET could be estimated and monitored in many more places than those equipped with Eddy Covariance or sap-flow measurements to advance the understanding of water balance changes in different successional ecosystems.

We present a new modeling approach analyzing and predicting the Transit Time Distribution (TTD) and the Response Time Distribution (RTD) from hourly to annual time scales as two distinct hydrological processes. The model integrates Isotope Hydrograph Separation (IHS) and the Instantaneous Unit Hydrograph (IUH) approach as a tool to provide a more realistic description of transit and response time of water in catchments. Individual event simulations and parameterizations were combined with long-term baseflow simulation and parameterizations; this provides a comprehensive picture of the catchment response for a long time span for the hydraulic and isotopic processes. The proposed method was tested in three Andean headwater catchments to compare the effects of land use on hydrological response and solute transport. Results show that the characteristics of events and antecedent conditions have a significant influence on TTD and RTD, but in general the RTD of the grassland dominated catchment is concentrated in the shorter time spans and has a higher cumulative TTD, while the forest dominated catchment has a relatively higher response distribution and lower cumulative TTD. The catchment where wetlands concentrate shows a flashier response, but wetlands also appear to prolong transit time.

Rainfall-runoff modelling in ungauged basins is still one of the greatest challenges in hydrological research. The lack of discharge data necessitates the establishment of new innovative approaches to guide hydrological modelling in ungauged basins. Besides the transfer of calibrated parameters from similar gauged catchments, the application of distributed data as a hydrological response in addition to discharge seems to be promising. A new approach to guide hydrological modelling based on explicit simulation of the spatial stream network was tested in four different catchments in Germany. In a first step we used a simplified version of the process-based model Hill-Vi together with regional climate normals to simulate stream networks. The calculation of gravity driven lateral subsurface and groundwater flow is used to identify patterns of stream cells, which were compared to reference stream networks and their degree of spatial agreement was evaluated. Significant differences between good and poor simulations could be distinguished and the corresponding parameter sets relate well with the hydrogeological properties of the catchments. The optimized parameters were subsequently used to simulate daily discharge using an observed time series of precipitation and air temperature. The performance was evaluated against observed discharge and water balance. This approach shows some promising results but also some limitations. Although the model's parsimonious model structure could be further improved regarding discharge recession and evapotranspiration, the performance was similar to regionalisation methods. Stream network modelling, which has minimal data requirements, seems to be a reasonable alternative for model development and parameter evaluation in ungauged basins.

Water table depth and its dynamics on hillslopes are often poorly predicted despite they control both water transit time within the catchment and solute fluxes at the catchment outlet. This paper analyses how relaxing the assumption of lateral homogeneity of physical properties can improve simulations of water table depth and dynamics. Four different spatial models relating hydraulic conductivity to topography have been tested: a simple linear relationship, a linear relationship with two different topographic indexes, two Ks domains with a transitional area. The Hill-Vi model has been modified to test these hypotheses. The studied catchment (Kervidy-Naizin, Western France) is underlain by schist crystalline bedrock. A shallow and perennial groundwater highly reactive to rainfall events mainly develops in the weathered saprolite layer. The results indicate that (1) discharge and the water table in the riparian zone are similarly predicted by the four models, (2) distinguishing two Ks domains constitutes the best model and slightly improves prediction of the water table upslope, and (3) including spatial variations in the other parameters such as porosity or rate of hydraulic conductivity decrease with depth does not improve the results. These results underline the necessity of better investigations of upslope areas in hillslope hydrology.

Preferential flow paths have been found to be important for runoff generation, solute transport, and slope stability in many areas around the world. Although many studies have identified the particular characteristics of individual features and measured the runoff generation and solute transport within hillslopes, very few studies have determined how individual features are hydraulically connected at a hillslope scale. In this study, we used dye staining and excavation to determine the morphology and spatial pattern of a preferential flow network over a large scale (30 m). We explore the feasibility of extending small-scale dye staining techniques to the hillslope scale. We determine the lateral preferential flow paths that are active during the steady-state flow conditions and their interaction with the surrounding soil matrix. We also calculate the velocities of the flow through each cross-section of the hillslope and compare them to hillslope scale applied tracer measurements. Finally, we investigate the relationship between the contributing area and the characteristics of the preferential flow paths. The experiment revealed that larger contributing areas coincided with highly developed and hydraulically connected preferential flow paths that had flow with little interaction with the surrounding soil matrix. We found evidence of subsurface erosion and deposition of soil and organic material laterally and vertically within the soil. These results are important because they add to the understanding of the runoff generation, solute transport, and slope stability of preferential flow-dominated hillslopes.

No abstract available.