Routine stream monitoring data enables the unravelling of hydrological pathways and transfers of agricultural contaminants through catchments

Catchment-scale understanding of water and contaminant fluxes through all pathways is essential to address land use and climate change impacts on freshwater. However, few options exist to obtain this understanding for the many catchments worldwide for which streamflow and low-frequency water chemistry, but little other data exists. We applied the Bayesian chemistry-assisted hydrograph separation and load partitioning model (BACH) to 47 catchments with widely differing characteristics. As BACH relies on concentration differences between pathways, chemodynamic behaviour of a water constituent indicates its likely suitability as tracer. Typical tracers (e.g. silica, chloride) were unavailable, but Electrical Conductivity and a few monitored nutrients proved chemodynamic in most catchments. Using one of two tracer combinations (Total Nitrogen + Electrical Conductivity, Total Nitrogen + Total Phosphorus) allowed in 85 % of the catchments to estimate streamflow contributions by near-surface (NS), shallow groundwater (SGW), and deep groundwater (DGW) pathways and pathway-specific tracer concentrations and yields with acceptable confidence. In 46 catchments, at least two pathways contributed ≥20 % of the streamflow, and all three ≥20 % in 12 catchments, cautioning against the notion of a single 'dominant' pathway. In contrast to hydrometric hydrograph separation, BACH allows differentiation between 'young' (NS + SGW) and 'old' (DGW) water, which is crucial for the understanding of pollution in catchments with strong temporal gradients in land use intensity. Consistent with generally increasing land use intensity, and groundwater denitrification occurring in some catchments, Total Nitrogen (TN) concentrations were in most catchments higher in NS and SGW compared to DGW. In most catchments, the greatest fraction of the TN yield was conveyed by SGW (≈ 40-90 %). Exceptions were wet and hilly catchments under bush, where the NS transferred most of the very low yields, and three young volcanic catchments where the DGW transferred the majority of the yield due to particularly high DGW flow contributions.

Spatial characteristics of hydrochemistry and stable isotopes in river and groundwater, and runoff components in the Shule River Basin, Northeastern of Tibet Plateau

Water resources play a crucial role in constraining the high-quality development of the arid, necessitating an in-depth investigation and understanding of hydrological processes, hydrochemical characteristics, and their influencing factors amidst climate change. This study meticulously examined and analyzed the hydrochemistry and stable isotope composition (δ18O and δD) of river and groundwater within the Shule River Basin (SRB). Results showed that both river (mean: 8.01) and groundwater (mean: 7.92) had alkaline pH values, while average total dissolved solids were measured at 709.25 mg/L in river and 861.88 mg/L in groundwater, indicating predominance of fresh water sources. HCO3-, SO42-, Na+ and Ca2+ were the most abundant ions, influenced by evaporation-crystallization processes and rock weathering. The dominated hydrochemical facies in both river and groundwater were Ca-HCO3 type in the upper (UR) and the middle reaches (MR), while Ca-Mg-Cl type in the lower reaches (LR). The local meteoric water line (LMWL) was defined as δD = 8.01δ18O + 18.48 (R2 = 0.98, n = 163; P < 0 0.001). The more negative δ18O and δD values in river and groundwater were plotted nearby and lower right of the LMWL, implying that the important recharge source of those waters is from precipitation. The relationship between river δ18O and elevation showed an increase of 0.14‰/100 m in the UR, but a negative correlation with a rate of -0.47‰/100 m in the MR and LR. Precipitation, groundwater, baseflow and meltwater accounted for 62.5%, 19.8%, 11.9% and 5.8% of the UR river, respectively, during the observed period, according to the end-member mixing analysis. These runoff components displayed distinct seasonal variations, primarily driven by precipitation during the early and groundwater/baseflow during the rapid and end-stage ablation periods. The observed alterations in hydrological elements present both opportunities and challenges for water resource management across the SRB, and adaptive measures have been proposed based on our study. These findings provide valuable insights into efficient utilization of water resources from water chemistry and environmental isotopes.

Improving the representation of stream water sources in surrogate nutrient models with water isotope data

An important part of meeting nutrient reduction goals in the lower Great Lakes basin and assessing the success of different land management strategies is modeling nutrient losses from agricultural land. This study aimed to improve the representation of water source contributions to streamflow in generalized additive models for predicting nutrient fluxes from three headwater agricultural streams in southern Ontario monitored during the Multi-Watershed Nutrient Study (MWNS). The previous development of these models represented baseflow contributions to streamflow using the baseflow proportion derived using an uncalibrated recursive digital filter. Recursive digital filters are commonly used to partition stream discharge into separate components from slower and faster pathways. In this study, we calibrated the recursive digital filter using stream water source information from stable isotopes of oxygen in water. Across sites, optimization of the filter parameters reduced bias in baseflow estimates by as much as 68 %. In most cases, calibrating the filter also improved agreement between filter-derived baseflow and baseflow calculated from isotope and streamflow data: the average Kling-Gupta Efficiencies using default and calibrated parameters were 0.44 and 0.82, respectively. When incorporated into the generalized additive models, the revised baseflow proportion predictor was more often statistically significant, improved model parsimony, and reduced prediction uncertainty. Moreover, this information allowed for a more rigorous interpretation of how different stream water sources influence nutrient losses from the agricultural MWNS watersheds.

Evaluating the contribution of subsurface drainage to watershed water yield using SWAT+ with groundwater modeling

Drainage outflow from artificial subsurface drains can be a significant contributor to watershed water yield in many humid regions of the world. Although many studies have undertaken to simulate hydrologic processes in drained watersheds, there is a need for a study that first, uses physically based spatially distributed modeling for both surface and subsurface processes; and second, quantifies the effect of surface and subsurface parameters on watershed drainage outflow. This study presents a modified version of the SWAT+ watershed model to address these objectives. The SWAT+ model includes the gwflow module, a new spatially distributed groundwater routine for calculating groundwater storage, groundwater head, and groundwater fluxes throughout the watershed using a grid cell approach, modified in this study to simulate the removal of groundwater by subsurface drains. The modeling approach is applied to the South Fork Watershed (583 km²), located in Iowa, USA, where most fields are drained artificially. The model is tested against measured streamflow, groundwater head at monitoring wells, and drainage outflow from a monitored subbasin. Sensitivity analysis is then applied to determine the land surface, subsurface, and drainage parameters that control subsurface drainage. Simulated drainage flow fractions (fraction of streamflow that originates from subsurface drainage) range from 0.37 to 0.54 during 2001-2012, with lower fractions occurring during years of high rainfall due to the increased volumes of surface runoff. Subsurface drainage comprises the vast majority of baseflow. Results indicate surface runoff and soil percolation parameters have the strongest effect on watershed-wide subsurface drainage rather than aquifer and drain properties, pointing to a holistic watershed approach to manage subsurface drainage. The modeling code presented herein can be used to simulate significant hydrologic fluxes in artificially drained watersheds worldwide.

Hydrograph apportionment of the Chandra River draining from a semi-arid region of the Upper Indus Basin, western Himalaya

Melting of snow and glaciers from the high-altitude Himalayan region is a significant water source to the major Himalayan rivers, especially in the upper Indus Basin (UIB), which contributes up to 70% of river discharge. Considering Indus Basin as a largest irrigation system dependent on snow and glacier melt runoff, it is imperative to study the rivers' current status and water budget. In this study we have performed a tracer-based hydrograph separation to quantify the contribution of seasonal snow, glacier melt, and groundwater to the Chandra River draining from a semi-arid region of the upper Indus basin, western Himalaya. Our study revealed a negligible control of summer (May-September 2017) precipitation and significant control of summer air temperature (May-September 2017) and winter precipitation over the Chandra River discharge, with 1 °C rise in air temperature leading to 22 m³s^-1 (15% of mean) increase in the river discharge (R² = 0.85; n = 541; p < 0.001). The hydrograph separation of the Chandra River suggests groundwater (38.3 ± 5.6%; 96.8 m³s^-1) as a significant source to the river runoff, followed by a direct contribution from glacier melt (30.9 ± 9%; 88.2 m³s^-1) and seasonal snowmelt (30.6 ± 5.7%; 84.2 m³s^-1), respectively, with negligible contribution from rainfall. Although groundwater is a significant contributor to the river runoff, the infiltration of seasonal snowmelt (54%) and glacier melt (46%) mostly contributed to the groundwater recharge. Present study establishes a linkage between seasonal snowmelt, glacier melt, groundwater, and the river runoff and would be useful to better model and predicts the future changes in the water resources of the upper Indus Basin.

Convenient use of electrical conductivity measurements to investigate hydrological processes in Alpine headwaters

Developing effective hydrological models for streamflow generation in Alpine catchments is challenging due to the inherent complexity of the intertwined processes controlling water transfer from hillslopes to streams and along the river network. Over the past decades, studies have proposed complementing traditional hydrological information with environmental tracer data, e.g. stable isotopes or electrical conductivity (EC), for different purposes such as the separation of streamflow components or the estimation of catchment mean residence time. In particular EC has been applied in Alpine environments mainly for hydrograph separation but also, more recently, considered as a possible proxy for streamflow (Q) prediction. The reason is simple: EC data loggers are convenient because of their relative low cost, easiness of installation and low maintenance, unlike traditional water stage gauges. However, EC time series require careful interpretation since electrical conductivity is influenced by a number of geochemical processes not always introduced in the analysis since these can be difficult to parametrize. Likewise, the relationship between EC and Q is very complex because it is characterized by hysteresis loops and often site specific. This study shows how the continuous monitoring of EC in Alpine catchments can be useful specifically for: hydrograph separation, including a proper quantification of uncertainty; process understanding of catchment functioning through the interpretation of hysteresis loops and time lags between EC and Q signals; and finally, water discharge estimation through calibrated functional EC-Q relationships. We discuss advantages and limitations of the use of EC in hydrology and provide information to encourage its use in studies dealing with streamflow generation dynamics in snow-dominated catchments.

Hydrograph separation of subsurface tile discharge

Baseflow is an important component of streamflow and watershed hydrologic budgets, yet quantifying the baseflow fraction of tile drainage has rarely been reported. In this study, we used two common hydrograph separation methods (local minimum method, recursive digital filter) to separate the discharge hydrographs from three drainage district tiles located in Iowa. Based on data collected from 2009 to 2013, annual baseflow ranged from 116 to 162 mm and comprised approximately 60% of the annual discharge. Baseflow was greatest during June (average of 34% of annual baseflow) and the March through August period produced 86% of the total annual baseflow. We found that the two methods of hydrograph separation produced similar results but the digital filter method was less erratic in estimating baseflow fraction. Study results can be used to better quantify hydrologic pathways in tiled landscapes and improve the design, implementation, and evaluation of nutrient reduction strategies.

Hydrochemistry of waters in snowpacks, lakes and streams of Mt. Dagu, eastern of Tibet Plateau

There is little available information on hydrochemistry of waters from glacial source to downstream of glacierized catchments. Here we examine the water chemistry of the snowpacks, lakes and streams at eight sampling sites within glacial basin in Mt. Dagu, east Tibetan Plateau. An air mass trajectory model, correlation analysis, Gibbs model, Piper diagram and hydrograph separation analysis were utilized to investigate the characteristics and solutes sources of these waters. Generally, the TDS (Total dissolved solids; 7.54, 13.95 and 18.70mg/L for snowpacks, lakes and streams respectively) and concentrations of main chemicals in all samples exhibited downstream trend from snowpacks to streams. Of the cations, Ca²⁺ appeared with the highest concentration followed by K⁺ and Na⁺. Of the anions, HCO₃^- was most abundant, followed by Cl^-, SO₄^2- and NO₃^-. For snowpack samples, the air masses arriving at the sampling sites were typically prevailing from the western Tibet Plateau, central Asia and the northern Mongolia plateau. The fine particulate matter in the Mt. Dagu snowpacks was most likely transported long range from three arid regions above-mentioned. High concentrations of SO₄^2- and NH₄⁺ in snowpacks, with twice as much NH₄⁺ as SO₄^2-, implying that the soluble part of the finer particles was transported as a form of ammonium sulfate. Rock weathering determined the ion components of the meltwater and the water could be classified as calcium and bicarbonate type based on the Piper diagram. The chemical contributions from glacier-snow meltwater were 20%-131% for lake and 5%-79% for stream, while the runoff contribution to lake varied from 65.4% to 84.9%, and 66.1% to 81.6% for stream. This study suggested that glacier-snow meltwater was the mainly runoff contributor to lake and stream water and that snowpack solutes derived from eolian additions exert a significant influence on lake and stream chemistry.