Streamflow Generation Processes Research Articles

JAMES BUTTLE REVIEW: Interflow, subsurface stormflow and throughflow: A synthesis of field work and modelling

AbstractInterflow, throughflow and subsurface stormflow are interchangeable terms that refer to the lateral subsurface flow above a restricting layer of lower hydraulic conductivity that occurs during and following storm events. Interflow (used here) is a more dominant process in steeper catchments with high infiltration capacity soils overlying a more impermeable soil or geologic layer. Interflow as a runoff process was first recognised in the early 1900s, yet hydrologists still struggle to predict its occurrence, persistence, importance, interaction with other streamflow generation processes, and potential to connect to valleys and streams during and following storms. We review the history of interflow research and address some of the challenges in understanding its role in runoff production. We argue that characterising the controls on interflow initiation and occurrence relies on detailed field observations of subsurface properties, which exist only in limited experimental settings. This data shortcoming contributes to our inability to predict interflow or determine its contribution to streamflow more broadly. There remain many opportunities to advance our understanding of interflow that include both modelling and experimental or observational approaches in hydrology.

Temporal Fusion Transformers for streamflow Prediction: Value of combining attention with recurrence

Over the past few decades, the hydrology community has witnessed notable advancements in streamflow prediction, particularly with the introduction of cutting-edge machine-learning algorithms. Recurrent neural networks, especially Long Short-Term Memory (LSTM) networks, have become popular due to their capacity to create precise forecasts and realistically mimic the system dynamics. Attention-based models, such as Transformers, can learn from the entire data sequence concurrently, a feature that LSTM does not have. This work tests the hypothesis that combining recurrence with attention can improve streamflow prediction. We set up the Temporal Fusion Transformer (TFT) architecture, a model that combines both of these aspects and has never been applied in hydrology before. We compare the performance of LSTM, Transformers, and TFT over 2,610 globally distributed catchments from the recently available Caravan dataset. Our results demonstrate that TFT indeed exceeds the performance benchmark set by the LSTM and Transformers for streamflow prediction, when models are trained individually for each basin. Additionally, being an explainable AI method, TFT helps in gaining insights into the streamflow generation processes.

Shifting groundwater fluxes in bedrock fractures: Evidence from stream water radon and water isotopes

Geologic features (e.g., fractures and alluvial fans) can play an important role in the locations and volumes of groundwater discharge and degree of groundwater-surface water (GW-SW) interactions. However, the role of these features in controlling GW-SW dynamics and streamflow generation processes are not well constrained. GW-SW interactions and streamflow generation processes are further complicated by variability in precipitation inputs from summer and fall monsoon rains, as well as declines in snowpack and changing melt dynamics driven by warming temperatures. Using high spatial and temporal resolution radon and water stable isotope sampling and a 1D groundwater flux model, we evaluated how groundwater contributions and GW-SW interactions varied along a stream reach impacted by fractures (fractured-zone) and downstream of the fractured hillslope (non-fractured zone) in Coal Creek, a Colorado River headwater stream affected by summer monsoons. During early summer, groundwater contributions from the fractured zone were high, but declined throughout the summer. Groundwater contributions from the non-fractured zone were constant throughout the summer and became proportionally more important later in the summer. We hypothesize that groundwater in the non-fractured zone is dominantly sourced from a high-storage alluvial fan at the base of a tributary that is connected to Coal Creek throughout the summer and provides consistent groundwater influx. Water isotope data revealed that Coal Creek responds quickly to incoming precipitation early in the summer, and summer precipitation becomes more important for streamflow generation later in the summer. We quantified the change in catchment dynamic storage and found it negatively related to stream water isotope values, and positively related to modeled groundwater discharge and the ratio of fractured zone to non-fractured zone groundwater. We interpret these relationships as declining hydrologic connectivity throughout the summer leading to late summer streamflow supported predominantly by shallow flow paths, with variable response to drying from geologic features based on their storage. As groundwater becomes more important for sustaining summer flows, quantifying local geologic controls on groundwater inputs and their response to variable moisture conditions may become critical for accurate predictions of streamflow.

Snow sampling strategy can bias estimation of meltwater fractions in isotope hydrograph separation

The Isotope based Hydrograph Separation (IHS) has been instrumental in understanding the partitioning of streamflow sources and processes. However, uncertainties persist in the accuracy of IHS estimations and the appropriate definition and sampling of endmembers. To address these uncertainties, we used field data of snowpack, snowfall, and snow meltwater isotopes (δ18O) from Pallas, Northern Finland to estimate the total meltwater contribution during the snowmelt period. We investigated the biases resulting from the application of different sampling strategies for event water endmember. The total meltwater contribution to streamflow was 59.6 % (±2% uncertainty) using the time-variant rolling runoff-corrected melt flux-weighted meltwater 18O isotope value. However, replacing it with either snowfall or winter snowpack 18O isotope weighted average values underestimated the meltwater contribution by 17.8 % or 22.6 %, respectively. Conversely, using time-variant instantaneous meltwater 18O isotope values overestimated the meltwater contribution by only 1.5 %. These discrepancies highlight the importance of choosing the appropriate endmember isotopes in IHS. The large differences in meltwater contribution for a 2-week peak discharge period based on different endmembers can lead to different interpretations of hydrological, ecohydrological, and biogeochemical processes. Thus, to better understand streamflow generation processes, we suggest using rolling runoff-corrected meltwater 18O or 2H isotope values in the IHS. In the absence of meltwater samples, the 18O or 2H isotope values of snowpack samples during the peak melt season may provide reasonable estimates of the meltwater contribution, with some minor underestimations. Our study highlights the importance of appropriate event meltwater endmember selection and sampling methodology for the IHS.

Influence of Critical Zone Architecture and Snowpack on Streamflow Generation Processes: A Mountain‐Meadow Headwater System in a Mediterranean Climate

AbstractObservations from a granitic watershed within a Mediterranean climate reveal the hydrologic and critical zone functioning of a perennial stream headwater and its upslope contributing area within a meadow system in the Sierra Nevada, California. Chemical analysis (diagnostic tools of mixing models, end member mixing analysis, tritium, etc.) and physical data (stream stage, piezometers, soil water, snowpack, etc.) indicate there are two primary pathways of water input into a headwater stream sourced from a mountain meadow. One input is a shallower and younger subsurface pathway with water that resembles snowpack chemistry, and the other a deeper and older subsurface pathway with water that reflects the chemistry of the groundwater derived from the contributing hillslopes. Multi‐year observations reveal that regardless of snowpack amount, during the period of peak hillslope infiltration, shallow and deep pathways in the hillslope behave similarly to initiate headwater streams. However, during summer dry periods, similarities in active pathways within the meadow center are not maintained between high and low snowpack years. With less snow, perennial groundwater discharge within the meadow center is eliminated, becoming only a seasonal source at the meadow's outlet. At the meadow's edge, geophysically observed downslope thinning in saprolite thickness creates reduced lateral transmissivity and initiation points for headwater streams via enhanced groundwater discharge of upslope water. Combined, these findings suggest how loss of snowpack and critical zone structure can together mediate hydrologic function in a wet meadow system in a Mediterranean climate. Creating new understanding about the stability of hydraulic functioning in headwater wet‐meadow systems under a changing climate.

Simulation of the impact of future changes in climate on the hydrology of Bow River headwater basins in the Canadian Rockies

This study diagnoses the impact of projected changes in climate and glacier cover on the hydrology of several natural flowing Bow River headwater basins in the Canadian Rockies: the Bow River at Lake Louise (∼420.7 km2), the Pipestone River near Lake Louise (∼304.2 km2), the Bow River at Banff (∼2192.2 km2) all of which drain the high elevation, snowy, partially glaciated Central Range, and the Elbow River at Calgary (∼1191.9 km2), which drains the drier Front Ranges and foothills, using models created using the modular, flexible, physically based Cold Regions Hydrological Modelling platform (CRHM). Hydrological models were constructed and parameterised in CRHM from local research results to include relevant streamflow generation processes for Canadian Rockies headwater basins, such as blowing snow, avalanching, snow interception and sublimation, energy budget snow and glacier melt, infiltration to frozen and unfrozen soils, hillslope sub-surface water redistribution, wetlands, lakes, evapotranspiration, groundwater flow, surface runoff and open channel flow. Surface layer outputs from Weather Research and Forecasting (WRF) model simulations for the current climate and for the late 21st century climate under a “business-as-usual” scenario, Representative Concentration Pathway 8.5 (RCP8.5) at 4-km resolution, were used to force model simulations to examine the climate change impact. A projected glacier cover under a “business-as-usual” scenario (RCP8.5) was incorporated to assess the impact of concomitant glacier cover decline. Uncalibrated model simulations for the current climate and glacier coverage showed useful predictions of snow accumulation, snowmelt, and streamflow when compared to surface observations from 2000 to 2015. Under the RCP8.5 climate change scenario, the basins of the Bow River at Banff and Elbow River at Calgary will warm up by 4.7 and 4.5 °C respectively and receive 12% to 15% more precipitation annually, with both basins experiencing a greater proportion of precipitation as rainfall. Peak snow accumulation in Bow River Basin will slightly rise by 3 mm, whilst it will drop by 20 mm in Elbow River Basin, and annual snowmelt volume will increase by 43 mm in Bow River Basin but decrease by 55 mm in Elbow River Basin. Snowcovered periods will decline by 37 and 46 days in Bow and Elbow river basins respectively due to suppressed snow redistribution by wind and gravity and earlier melt. The shorter snowcovered period and warmer, wetter climate will increase evapotranspiration and glacier melt, if the glaciers were held constant, and decrease sublimation, lake levels, soil moisture and groundwater levels. The hydrological responses of the basins will differ despite similar climate changes because of differing biophysical characteristics, climates and hydrological processes generating runoff. Climate change with concomitant glacier decline is predicted to increase the peak discharge and mean annual water yield by 12.23 m3 s−1 (+11%) and 11% in the higher elevation basins of the Bow River but will decrease the mean annual peak discharge by 3.58 m3 s−1 (-9%) and increase the mean annual water yield by 18% in the lower elevation basin of the Elbow River. This shows complex and compensatory hydrological process responses to climate change with the reduced glacier contribution reducing the impact of higher precipitation in high elevation headwaters and drier soil conditions and lower spring snowpacks reducing peak discharges despite increased precipitation during spring runoff in the Front Range and foothills headwaters under a warmer climate.

Parameter estimation for models of major rainfall induced floods in ungaged mountain basins of Colorado

Study regionRocky Mountains in Colorado including Front and San Juan ranges Study focusHydrologic modeling of ungaged mountain basins plays an important role in ensuring the safety of Colorado’s dams. Recent research has shown that infiltration-excess runoff, saturation-excess runoff, and subsurface stormflow can all contribute to streamflow generated by extreme rainfall in Colorado’s mountains, and the soil moisture accounting (SMA) method in HEC-HMS has been suggested as an appropriate approach to model these mechanisms. However, SMA requires estimation of parameters that have not been previously considered for large floods. The objective of this work is to develop methods to estimate SMA parameters for modeling the response to extreme rainfall in ungaged mountain basins of Colorado where calibration to observed discharges is not possible. New hydrological insights for the regionThe proposed parameter estimation methods are tested for five basins in the Front Range and three basins in the San Juan Range that have streamflow data available for major flood events. For these historical events, the proposed uncalibrated models produce the same streamflow generation processes as calibrated models, and the predicted peak discharges from the uncalibrated models usually have 25% errors or smaller. For design storms, the peak discharges from the uncalibrated models can differ substantially from the peak discharges from calibrated models but are conservative relative to the envelope of observed peak discharges.

Hydrochemistry of surface waters in a permafrost headwater catchment in the Northeastern Tibetan Plateau

The permafrost headwater catchments in the Tibetan Plateau (TP) have experienced extensive permafrost degradation, which may cause major changes in riverine solutes. However, surface water hydrochemistry and its influencing factors in such catchments are poorly understood. Hydrochemistry data for different surface waters were obtained for the Yakou catchment in the Northeastern TP. The results indicate that the ionic and organic concentrations of frozen soil seeps (FSS) were higher on the north-facing slope compared to the south-facing slope, and that FSS may be involved in streamflow generation processes and in determining the spatial pattern of riverine solutes. The north-facing slope of the catchment has a thin active layer and wet moisture conditions compared to the south-facing slope; hence supra-permafrost water, with high ionic concentrations, can drain to the ground surface as FSS in the riparian zone and then recharge the surface ponding water and the main stream water. The high ionic concentrations of the supra-permafrost water and FSS can be attributed to intense rock weathering and evaporative effect, together with the high mobility of elements and the transport of organic matter. The tributaries, with low ionic concentrations, comprise a mixture of infiltrating precipitation and diluted supra-permafrost water. Carbonate weathering is the dominant weathering type within the catchment, but the weathering of evaporite and silicate is more important on the north-facing and south-facing slopes, respectively, and chemical weathering on the north-facing slope may be enhanced by strong physical erosion during repeated freeze–thaw cycles due to the wet conditions. The results indicate that the surface water hydrochemistry is heterogeneous on the different hillslope units, and that a thicker active layer under climate change may lead to a shift of hydrological and hydrochemical pathways, and thus a decrease in water and solute flux from the hillslopes, with underlying permafrost, to the river channel.

Soil frost controls streamflow generation processes in headwater catchments

The relationship between snowmelt and spring streamflow is changing under warming temperatures and diminishing snowpack. At the same time, the hydrologic connectivity across catchment landscape elements, such as snowpack and surface wetlands, can play a critical role in controlling the routing of snowmelt to streams. The role of hydrologic connectivity is important in headwater regions of the continental northern latitudes, where catchments have low topographic relief and seasonally frozen ground. Nevertheless, the effects of soil frost on the sequence, timing, and magnitudes of hydrologic events that drive the movement of water from a snowpack to a stream are not fully understood. Therefore, we examine two questions: First, what is the flowpath that snow melt and precipitation from spring rain events takes to generate spring streamflow, and second, what hydrologic, climatic, or landscape variables exert the most control on the magnitude of streamflow? Here, we use long-term hydrological records from the two reference basins at the Marcell Experimental Forest in northern Minnesota to analyze the cascading effects across precipitation, snow, water table elevation, soil frost, and streamflow in peatland-dominated headwater catchments. We identify a sequence of fill-and-spill effects across the landscape that control the timing of spring streamflow generation. Then, we use stepwise regression to show that soil frost is a key supporting predictor for both the magnitude of streamflow in the spring as it adds significantly to the predictive power of precipitation and water table elevation. Our results highlight the importance of recognizing the role of soil frost, when present, on the partitioning of snowmelt between overland runoff and water table recharge during the critical snowmelt period, as well as the later partitioning between evapotranspiration and subsurface flows.

Similarity of catchment dynamics based on the interaction between streamflow and forcing time series: Use of a transfer entropy signature

The transfer of the hydrological information between catchments is founded on the definition of hydrological similarity, which is in turn strictly connected to the features to be regionalised. In order to characterise the catchment behaviour in the streamflow generation processes, the similarity should reflect also the interaction between meteorological forcings and river streamflow time series. While previous hydrological research has identified basins with similar meteorological forcings (i.e. similarity of climate) or with similar streamflow time-series (i.e. similarity of runoff response), the present work proposes, for the first time, to quantify the interaction between the entire time-series of different forcing data and streamflow observations, to be considered as a novel hydrological signature and used as a catchment similarity metric.In particular, the present study proposes the use of a multi-variate entropy-based measure, the so-called transfer entropy, a time-asymmetric quantity which analyses the interaction between different signals. The concept of transfer entropy is applied for identifying the dominant hydrological processes occurring in a catchment, measuring the transfer of information from different meteorological forcings over the catchment to the corresponding observed time series of daily streamflow at the basin outlet. The resulting transfer entropy values are then used as signatures to characterise the main catchment dynamics, and a classification of the basins region is obtained assuming that similar values of transfer entropy correspond to hydrologically similar basins.The methodology is tested on a densely-gauged set of more than 200 catchments across Austria and the outcomes of the approach are evaluated against a set of morpho-climatic catchment attributes and typical streamflow signatures.Despite the limitations, the method is able to distinguish the predominant or partial role of snow melt and evapotranspiration across the dataset, helping to assess differences in catchment response time and to highlight the role of very high orographic precipitation in catchments with a dominant snow regime. The study demonstrates the potential of transfer entropy as complementary to consolidated streamflow signatures for assessing hydrological similarity and for quantifying the connection between different catchment processes.

Hydrology and biogeochemistry datasets from Sleepers River Research Watershed, Danville, Vermont, USA

AbstractThe Sleepers River Research Watershed (SRRW) in Danville, Vermont, USA, was established in 1957 by the U.S Department of Agriculture, Agricultural Research Service (ARS). The 111‐km2 area was selected as a representative mixed‐land‐use upland landscape in the northeastern USA, with the objective to study fundamental hydrological processes. The U.S. Geological Survey (USGS) currently operates SRRW and is the fourth federal agency to lead it. SRRW has hosted historic research on snowmelt modelling and streamflow generation processes, and USGS has augmented this hydrological foundation with biogeochemical cycling research. This data note describes five freely available data releases in stable repositories, extending to the late 2010s, from the indicated start date: aqueous chemistry and isotopes (1991); precipitation, water temperature, and streamflow (1991); groundwater levels (1991); ground frost depth (1983); and snow depth and water equivalent (1960). USGS research at SRRW is expected to continue, and we plan to update these releases every 1–3 years.

Hunting for Information in Streamflow Signatures to Improve Modelled Drainage

About half of the Danish agricultural land is drained artificially. Those drains, mostly in the form of tile drains, have a significant effect on the hydrological cycle. Consequently, the drainage system must also be represented in hydrological models that are used to simulate, for example, the transport and retention of chemicals. However, representation of drainage in large-scale hydrological models is challenging due to scale issues, lacking data on the distribution of drain infrastructure, and lacking drain flow observations. This calls for more indirect methods to inform such models. Here, we investigate the hypothesis that drain flow leaves a signal in streamflow signatures, as it represents a distinct streamflow generation process. Streamflow signatures are indices characterizing hydrological behaviour based on the hydrograph. Using machine learning regressors, we show that there is a correlation between signatures of simulated streamflow and simulated drain fraction. Based on these insights, signatures relevant to drain flow are incorporated in hydrological model calibration. A distributed coupled groundwater–surface water model of the Norsminde catchment, Denmark (145 km2) is set up. Calibration scenarios are defined with different objective functions; either using conventional stream flow metrics only, or a combination with hydrological signatures. We then evaluate the results from the different scenarios in terms of how well the models reproduce observed drain flow and spatial drainage patterns. Overall, the simulation of drain in the models is satisfactory. However, it remains challenging to find a direct link between signatures and an improvement in representation of drainage. This is likely attributable to model structural issues and lacking flexibility in model parameterization.

Taking theory to the field: streamflow generation mechanisms in an intermittent Mediterranean catchment

Abstract. Streamflow dynamics for non-perennial networks remain poorly understood. The highly nonlinear unsaturated dynamics associated with the transitions between wetting and drying in non-perennial systems make modelling cumbersome. This has stifled previous modelling attempts and alludes to why there is still a knowledge gap. In this study, we first construct a conceptual model of the physical processes of streamflow generation in an intermittent river system in South Australia, based on the hypothesis that the vertical and longitudinal soil heterogeneity and topography in a basin control short term (fast flows), seasonal (slow flow), and a mixture of these two. We then construct and parameterise a fully integrated surface–subsurface hydrologic model to examine patterns and mechanisms of streamflow generation within the catchment. A set of scenarios are explored to understand the influences of topography and soil heterogeneity across the catchment. The results showed that distinct flow generation mechanisms develop in the three conceptualised areas with marked soil and topographic characteristics and suggested that capturing the order of magnitude for the average hydraulic conductivity of each soil type across the catchment was more important than pinpointing exact soil hydraulic properties. This study augments our understanding of catchment-scale streamflow generation processes, while also providing insight on the challenges of implementing physically based integrated surface–subsurface hydrological models in non-perennial stream catchments.

An overview of the hydrology of non‐perennial rivers and streams

AbstractNon‐perennial rivers and streams are ubiquitous on our planet. Although several metrics have been used to statistically group or compare streamflow characteristics, there is currently no widely used definition of how many days or over what reach length surface flow must cease in order to classify a river as non‐perennial. At the same time, the breadth of climate and geographic settings for non‐perennial rivers leads to diversity in their flow regimes, such as how often or how quickly they go dry. These rivers have a rich and expanding body of literature addressing their ecologic and geomorphic features, but are often said to be ignored by hydrologists. Yet there is much we do know about their hydrology in terms of streamflow generation processes, water losses, and variability in flow. We also know that while they are prevalent in arid regions, they occur across all climate types and experience a diverse set of natural and anthropogenic controls on streamflow. Furthermore, measuring and modeling the hydrology of these rivers presents a distinct set of challenges, and there are many research directions, which still require further attention. Therefore, we present an overview of the current understanding, methodologic challenges, knowledge gaps, and research directions for hydrologic understanding of non‐perennial rivers; critical topics in light of both growing global water scarcity and ever‐changing laws and policies that dictate whether and how much environmental protection these rivers receive.This article is categorized under: Science of Water > Science of Water

LINKING DISSOLVED ORGANIC CARBON SPATIAL HETEROGENEITY TO GROUNDWATER DYNAMICS AND SOIL ORGANIC CARBON CONTENT IN A SUBTROPICAL HEADWATER CATCHMENT

Storage and fluxes of dissolved organic carbon (DOC) are triggered by interactions between hydrological and biogeochemical processes in the landscape. However, the relationship between DOC sources and catchment water storage is poorly understood due to uncertainties about the interactions between shallow groundwater (GW) dynamics and soil organic carbon (SOC). To investigate the relationships between catchment functioning and the spatiotemporal variability of DOC concentrations in shallow GW, we monitored a network of 12 shallow GW wells in a 0.24 km² catchment in Southern Brazil using DOC analysis and SOC measurements together with GW and stream gauging and statistical analyses. The greatest DOC concentration was observed in wells in the zero-order basin with highly dynamic saturation of the soil profile (50% with water table level below 0.5 m, mean DOC = 4.83 ± 2.53 mg L-1), where GW was characterized by a fast rise and slow decay during storm events. In these areas, the relationship between DOC concentration and GW level followed a decay with soil depth and higher DOC concentrations were found when GW level was within 0.1 m of the soil surface (ranging from 2.00 to 14.8 mg L-1). Wells in mid-slope locations were characterized by the lowest SOC content and DOC concentration in GW (1.96 ± 0.91 mg L-1). During the driest period, lower slopes showed higher DOC concentrations than riparian zones due to the persistent near saturation of GW and the occurrence of DOC-enriched water due to the high SOC content in superficial layers. When the riparian zone and hillslopes showed a low water storage deficit and a quasi-permanent hydrological connectivity, DOC produced between rainfall events was transported continuously from saturated areas in hillslopes to riparian zones, with a higher DOC concentration in wells near streams. This study demonstrates that spatial sources of DOC in a subtropical headwater system are linked to streamflow generation processes and GW dynamics.

Summer runoff generation in foothill catchments of the Colorado Front Range

Climatic shifts, disturbances, and land-use change can alter hydrologic flowpaths, water quality, and water supply to downstream communities. Prior research investigating streamflow generation processes in mountainous areas has largely focused on high-elevation alpine and subalpine catchments; less is known about these processes in lower-elevation foothills and montane catchments. In these lower-elevation ecoregions, precipitation shifts seasonally from snow to rain, which can result in differing seasonal flowpaths. We analyzed stream water for electrical conductivity, SiO2, Ca, Mg, Na, Cl, SO4, K, and dissolved organic carbon on both a weekly and storm event basis from April to August 2018 in three small (<10 km2) foothill catchments, and one larger (63.2 km2) catchment extending from the foothills to the subalpine ecoregions, in the Colorado Front Range. Using two end-member hydrograph separations and concentration-runoff relationships, we inferred the dominant catchment-scale flowpaths of precipitation to the streams. We selected catchments with varying land use to investigate the relationship between these characteristics and hydrologic flowpaths. We observed that concentrations of lithogenic constituents generally increased and dissolved organic carbon decreased as seasonal runoff decreased in the three foothill catchments, reflecting a transition from shallow subsurface flowpaths to deeper subsurface flowpaths. Elevated SO4 and Cl concentrations during low-flow periods in two of our catchments suggest that historical or current anthropogenic activities, such as mining, application of road salt, and/or near-stream septic systems, affect local stream and groundwater chemistry. In a foothill catchment with anthropogenic and geologic impervious surfaces, streamflow during storm responses was sourced from faster, surficial flowpaths compared to a less disturbed neighboring catchment, highlighting the influence of anthropogenic land-use on runoff generation. This study provides insight into the fundamental hydrology of foothill catchments and how they may function in the future with human development, precipitation shifts and disturbances.

Combined use of stable isotopes and hydrochemical characteristics to determine streamflow sources in the Jonkershoek catchment, South Africa

ABSTRACT The stable isotopes of water (18O, 2H) and other hydrochemical properties were used to assess seasonal changes between sources of stream water and flow pathways in baseflow conditions for two headwater sub-catchments (∼3 km2) of the Jonkershoek, a mountainous catchment in the Western Cape of South Africa. The sub-catchments differ in land cover, one is dominated by indigenous fynbos vegetation and the other by pine plantation. Stream water, higher elevation springs, and lower elevation groundwater samples were collected monthly from January 2018 to January 2019, and were analysed for stable isotopes, electrical conductivity (EC) and pH. The stream water isotopic values resembled those of groundwater during the dry and wet seasons. Results indicated a steady contribution of spring discharge to streams during the dry season, with rainfall contribution less evident. Wet season flows were attributed to interflow including perennial and ephemeral springs. Spatial variations in EC between the sub-catchments were indicative of a greater proportional input from subsurface water, more evaporation and/or more mineral dissolution in the pine-dominated sub-catchment. The spatiotemporal variability in EC was significantly larger than for the stable isotopes. These findings enabled conceptualization of streamflow generation processes that can support strategic water resource management practices in this region.

Technical note: An improved discharge sensitivity metric for young water fractions

Abstract. Recent virtual and experimental investigations have shown that the young water fraction Fyw (i.e. the proportion of catchment outflow younger than circa 2–3 months) increases with discharge in most catchments. The discharge sensitivity of Fyw has been defined as the rate of increase in Fyw with increasing discharge (Q) and has been estimated by the linear regression slope between Fyw and Q, hereafter called DS(Q). The combined use of both metrics, Fyw and DS(Q), provides a promising method for catchment inter-comparison studies that seek to understand streamflow generation processes. Here we explore the discharge sensitivity of Fyw in the intensively sampled small Mediterranean research catchment Can Vila. Intensive sampling of high flows at Can Vila allows young water fractions to be estimated for the far upper tail of the flow frequency distribution. These young water fractions converge toward 1 at the highest flows, illustrating a conceptual limitation in the linear regression method for estimating DS(Q) as a metric of discharge sensitivity: Fyw cannot grow with discharge indefinitely, since the fraction of young water in discharge can never be larger than 1. Here we propose to quantify discharge sensitivity by the parameter of an exponential-type equation that expresses how Fyw varies with discharge. The exponential parameter (Sd) approximates DS(Q) at moderate discharges where Fyw is well below 1; however, the exponential equation and its discharge sensitivity metric better capture the non-linear relationship between Fyw and Q and are robust with respect to changes in the range of sampled discharges, allowing comparisons between catchments with strongly contrasting flow regimes.

Saturated areas through the lens: 1. Spatio‐temporal variability of surface saturation documented through thermal infrared imagery

AbstractSurface saturated areas are key features in generating run‐off. A detailed characterization of the expansion and contraction of surface saturation in riparian zones and its connectivity to the stream is fundamental to improve our understanding of the spatial and temporal variability of streamflow generation processes. In this first contribution of a series of two papers, we used ground‐based thermal infrared imagery for characterizing riparian surface saturation seasonal dynamics of seven different riparian areas in the Weierbach catchment (0.42 km2), a small forested catchment in Luxembourg. We collected biweekly panoramic images of the seven areas over a period of 2 years. We identified the extension of saturation in each collected panoramic image (i.e., percentage of pixels corresponding to saturated surfaces in each riparian area) to generate time series of surface saturation. Riparian surface saturation in all areas was seasonally variable, and its dynamics were in accordance with lower hillslope groundwater level fluctuations. Surface saturation in the different areas related to catchment outlet discharge through power law relationships. Differences in these relationships for different areas could be associated with the location of the areas along the stream network and to a possible influence of local riparian morphology on the development of surface saturation, suggesting a certain degree of intracatchment heterogeneity. This study paves the way for a subsequent investigation of the spatio‐temporal variability of streamflow generation in the catchment, presented in our second contribution.

Catchment Travel Times From Composite StorAge Selection Functions Representing the Superposition of Streamflow Generation Processes

AbstractCatchment travel times integrate the multitude of hydrological flow processes and provide insights into catchment functioning. StorAge Selection (SAS) functions describe how residence times of water in storage are related to travel times of water in catchment outflows. As such, SAS functions are useful to summarize transport processes in catchments and are ideal to simulate catchment outflows and the concentrations of various solutes and tracers. Recent studies suggested that using one probability distribution function (pdf) for SAS functions may not account for all transport processes in various hydrological systems. In this study we introduced a composite streamflow SAS function defined as a weighted sum of several pdfs (called components), using a uniform pdf for the youngest water and two gamma pdfs for the oldest water. The novel parameterization of this SAS function is nonlinear with respect to catchment storage and to a proxy of storage variations. This composite streamflow SAS function was needed to obtain realistic simulations of the complex high‐resolution (subdaily) stream deuterium (δ2H) dynamics measured for 2 years in the Weierbach, a forested headwater catchment in Luxembourg, whereas various SAS functions using only one component (i.e., pdf) failed. The three components of the composite streamflow SAS function confirmed the superposition of streamflow generation mechanisms with contrasting travel times suggested by previous experimental and modeling studies. Our work suggests that in some catchments, composite SAS functions may offer a more realistic perspective on transport processes derived from high‐resolution tracer data than simpler SAS functions.