Increased urbanization, coupled with the projected impacts of climatic change, mandates further evaluation of the impact of urban development on water flow paths to guide sustainable land use planning. Though the general urbanization impacts of increased storm runoff peaks and reduced baseflows are well known; how the complex, non-stationary interaction of the dominant water fluxes within dynamic urban water stores sustain streamflow regimes over longer periods of time are less well quantified. In particular, there is a challenge in how hydrological modelling should integrate the juxtaposition of rapid and slower flow pathways of the urban ‘karst’ landscape and different approaches need evaluation. In this context, we utilized hydrological and water stable isotope datasets within a modelling framework that combined the commonly used HEC urban runoff model along with a simple hydrological tracer module and transit time modelling to evaluate the spatial and temporal variation of water flow paths and ages within a heavily urbanized 217km 2 catchment in Berlin, Germany. Deeper groundwater was the primary flow component within less urbanized regions of the catchments, with increased direct runoff and shallow subsurface contributions in more urbanized areas near the catchment outlet. The addition of wastewater effluent in the mid-reaches of the catchment was the dominant water supply to the lower stream, and sustained baseflows during the summer months. Water ages from each modelling approach imitated flow contributions and opportunity for mixing with subsurface storage; with older water and lower young water contributions in less urbanized sub-catchments and younger water and higher young water contributions in more urbanized regions. The results form a first step towards more integrated modelling tools for similar peri-urban catchments, given the potential limitations of more simple model frameworks. The results have broader implications for assessing the uncertainty in evaluating urban impacts on hydrological function under environmental change.

Restoring riparian shade helps maintain healthy stream ecosystems and computer models are helpful in guiding restoration. The effects of meandering and orientation on shade to direct and diffuse lighting was investigated using a simplified model of riparian vegetation. Previous studies have shown that straight streams oriented east-west (EW) experience higher lighting than those oriented north-south (NS). Meandering decreased reach-averaged direct lighting in EW valley streams but increased lighting in NS valley streams. When meander amplitude equalled 50% of meander wavelength, lighting of streams in EW and NS-aligned valleys was similar. Meandering and valley orientation had little effect on diffuse lighting. In EW valleys lighting was highest at stream bends and a strategy to reduce lighting quickly would be to prioritise planting of tall vegetation on bends. In NS valleys lighting was lowest at bends, suggesting the opposite planting strategy. Shade exceeding 70%, a criterion for preventing nuisance aquatic plant growths and stressfully high water temperatures, occurs in ‘typical’ meandering streams on cloudless, mid-summer days at mid-latitudes once H/ W exceeds about 2. The model over-simplifies riparian vegetation and the effects of canopy shape and overhang merit further investigation.

Water sources carry chemicals that can have a significant impact on the water environment of a river network, and understanding the contribution of different water sources to the river network can help to manage the pollution of the river network at its source. Hydrological connectivity of a river network affects the self-purification capacity and flood prevention capacity of the river. Thus an isotope tracer approach was applied to figure out the contribution rate of different water bodies to a river network and hydrological connectivity was quantified by introducing retention rate. Changzhou city was selected as the study area because it is an urbanized city with the characteristics of plain river network and it is faced with poor hydrological connectivity due to artificial constructions (dams and pumps) and human activity (urbanization). River water, well water (shallow groundwater), lake water and rainfall were collected during the flood season and nonflood season, and hydrogen and oxygen isotopes were determined. The temporal and spatial variations in hydrogen and oxygen isotopes in different water bodies and the state of the water cycle in different water bodies were analyzed. IsoSource and MixSIAR models were established to analyze the contribution rate of river network water sources in the study area, and their effectiveness was compared. Results of MixSIAR model were selected to evaluate the hydrological connectivity of the river network in the study area, providing a method to quantify the hydrological connectivity of specific river of the river network in Changzhou. This method could also be applied to other urban plain river network area to study its river connectivity.

We discuss a study that aimed to understand the genesis and inflow conditions of abstracted water by intake with induced infiltration located at a mountain river. A simple approach based on a combination of two research techniques was used: two-component water mixing modelling and studies of the variability of concentrations of environmental tracers in a dynamic test. This approach is versatile, easy to apply and modify, and can be good method for controlling surface/groundwater interactions. We used the gas tracers because gas exchange and dissolution of gases in infiltrating rainwater gives it a unique gas signature that is largely retained in groundwater. We focus on understanding river/aquifer interactions at the scale of reach of an intake. To understand these issues, a two-day field hydrogeological experiment based on a pumping test of increasing intensity was conducted. At each pumping stage, groundwater and river samples were collected to determine the concentration of noble gases, CFCs, SF6, stable isotope content, and chemical composition of the water. The study results indicate a short pressure propagation time between the intake and the river, which results in inflow of water already at low water abstract by the intake. Despite limited river water inflow, there was a continuous hydraulic contact river/aquifer, largely independent of changes in pumping intensity. As pumping intensity increased, the river/aquifer hydraulic system tended towards stabilised conditions of mixing and inflow to the intake. In general, a small proportion of river water (up to 26%) was found in the exchange flux between the river and groundwater.

Water bodies such as lakes and reservoirs affect the regional climate by acting as heat sinks and sources through the evaporation of substantial quantities of water over several months of the year. Unfortunately, energy exchange observations between inland water bodies and the atmosphere remain rare in northeastern North America, which has one of the highest densities of lakes in the world. This study helps to fill this gap by analyzing field observations collected from a subarctic hydropower reservoir (50.69°N, 63.24°W) characterized by a mean depth of 44 m and a surface area of 85 km 2. Two eddy covariance (EC) systems, one on a raft and one onshore, were deployed from 27 June 2018 to 12 June 2022. The thermal regime of the reservoir was monitored using vertical chains of thermistors. Results indicate a mean annual evaporation rate of 590 ± 66 mm (~70% of the annual precipitation), with 84% of the evaporation occurring at a high rate from August to freeze-up in late December through episodic pulses. It was difficult to close the energy balance because of the magnitude and the large time lag of the heat storage term. In order to circumvent this problem, we opted to perform calculations for a year that started from the first of March, as the heat storage in the water column was at its lowest at that point and could thus be ignored. From June to December, monthly Bowen ratios increased from near-zero negative values to about 1.5. After September, due to smaller vapor pressure deficits, latent heat fluxes steadily declined until an ice cover sealed up the reservoir. Two opposite diurnal cycles of sensible and latent heat fluxes were revealed during the open water period, with sensible heat fluxes peaking at night and latent heat fluxes peaking in the afternoon.

Since this is a HP today contribution I understand that it should not have an abstract

It has been proposed that conservation laws might not be beneficial for accurate hydrological modeling due to errors in input (precipitation) and target (streamflow) data (particularly at the event time scale), and this might explain why deep learning models (which are not based on enforcing closure) can out-perform catchment-scale conceptual and process-based models at predicting streamflow. We test this hypothesis at the event and multi-year time scale using physics-informed (mass conserving) machine learning and find that: (1) enforcing closure in the rainfall-runoff mass balance does appear to harm the overall skill of hydrological models, (2) deep learning models learn to account for spatiotemporally variable biases in data (3) however this “closure” effect accounts for only a small fraction of the difference in predictive skill between deep learning and conceptual models.

Particle tracers are sometimes used to track sources and sinks of riverine particulate and contaminant transport. A potentially new particle tracer is ~200 nm sized superparamagnetic silica encapsulated DNA (SiDNAFe). The main objective of this research was to understand and quantify the settling and aggregation behaviour of SiDNAFe in river waters. Our results indicated, that in quiescent conditions, more than 60% of SiDNAFe settled within 30 hours, starting with a rapid settling phase followed by an exponential-like slow settling phase in the three river waters we used (Meuse, Merkske, and Strijbeek) plus MilliQ water. From this, we inferred that the rapid SiDNAFe settling was mainly due to homo-aggregation and not due to hetero-aggregation (e.g., with particulate matter present in river water). Incorporating a first-order mass loss term which mimics the exponential phase of the settling in quiescent conditions seems to be an adequate step forward when modelling the transport of SiDNAFe in river injection experiments. Furthermore, we validated the applicability of magnetic separation and up-concentration of SiDNAFe in real river waters, which is an important advantage for carrying out field-scale SiDNAFe tracing experiments.

Sustainable water management in semi-arid agriculture practices requires quantitative knowledge of water fluxes within the soil-vegetation-atmosphere system. Therefore, we used stable-isotope approaches to evaluate evaporation (Ea), transpiration (Ta), and groundwater recharge (R) at sites in Senegal's Groundnut basin and Ferlo Valley pasture region during the pre-monsoon, monsoon, and post-monsoon seasons of 2021. The approaches were based upon (i) the isothermal evaporation model (for quantifying Ea); (ii) water and isotope mass balances (to partition Ea and Ta for groundnut and pasture); and (iii) the piston displacement method (for estimating R). Ea losses derived from the isothermal evaporation model corresponded primarily to Stage II evaporation, and ranged from 0.02–0.09 mm d-1 in the Groundnut basin, versus 0.02–0.11 mm d-1 in Ferlo. At the groundnut site, Ea rates ranged from 0.01 to 0.69 mm d-1; Ta was in the range 0.55–2.29 mm d-1; and the Ta/ETa ratio was 74–90%. At the pasture site, the ranges were 0.02–0.39 mm d-1 for Ea; 0.9–1.69 mm d-1 for Ta; and 62–90 % for Ta/ETa.  The ETa value derived for the groundnut site via the isotope approach was similar to those from eddy covariance measurements, and also to the results from a previous validated HYDRUS-1D model. However, the HYDRUS-1D model gave a lower Ta/ETa ratio (23.2%). The computed groundwater recharge for the groundnut site amounted to less than 2% of the local annual precipitation. Recommendations are made regarding protocols for preventing changes to isotopic compositions of water in samples that are collected in remote arid regions, but must be analysed days later. The article ends with suggestions for studies to follow up on evidence that local aquifers are being recharged via preferential pathways

We examined conditions that form or prevent thermal stratification in river pools using field measurements and statistical and three-dimensional (3D) computational fluid dynamics (CFD) modeling. Our motivation is to identify variables that control stratification for exploitation to enhance or prevent thermal gradients as needed to benefit species in rivers. One study pool (UT) is above water storage reservoirs and receives natural flows, and the other pool (PT) is regulated and receives unnaturally high, cold water in summer on the Trinity River, California. Thermal stratification formed in UT pool in spring at a critical flow of 1.01 m 3/s, peaked at 8.1 oC in summer, and exhibited diurnal formation and destruction under sub-critical flows until fall. At PT pool, the 14.2 m 3/s baseflow caused mixing that prevented stratification and formed a spatially homogenous thermal environment. Statistical modeling indicated the daily range in air and inlet water temperature at UT pool best correlated with the occurrence and strength of stratification but were progressively irrelevant as flows increased above the critical value. The 3D CFD model was verified by predicting the observed critical flow and dynamics of stratification at UT pool and isohytes observed at PT pool. The 3D model was then used to explore the thermal stratification process. Results confirmed low flows are the main variable for stratification to form, and the daily range in inlet water temperature drives the strength of the thermal gradient. The model estimated a critical discharge at PT pool of 2.0 m 3/s, twice that for UT pool owing to its 2.6-times larger area, suggesting critical flows scale with pool size. Results show that releasing critical and lower flows in summer on regulated streams may conserve water and provide thermal gradients that benefit poikilothermic species; alternately, higher than critical flows can prevent stratification where needed to improve water quality.

Flooding is a frequent disaster that has a wide-spread footprint globally with significant financial and societal impacts. With availability of Earth observation data from private and public entities at varying spatial, temporal, and spectral resolution as well as data from crowdsourcing, there is no shortage of models. In fact, models and algorithms are abundant and proliferating. However, the question remains where is a global flood model when we need one? Just because models are available does not mean they are usable or accessible and adequate for emergency managers, first responders and other stakeholders who use the model outputs for preparedness, response and resource planning. Often the issue of usability stems from the fact that the models are not always reproducible or replicable. The accuracy and uncertainty associated with the models and how they change based on the scale of analysis and the resolution of input and output datasets are often not communicated properly to stakeholders so they can be part of their decision-making process. The proliferation of machine learning and data driven models that rely on historical data also adds to this problem. This paper discusses several important issues associated with global flood models and provides recommendations that could be used to increase the usability of these models.

In arid area, the liquid water and water vapor states in soil profiles and fluxes at the upper and bottom interfaces are extremely complex due to heterogeneity of soil textures and the driving forces of heat and matrix potential. In this study, we used Hydrus-1D to simultaneously simulate liquid water, water vapor, and heat transports based on the observed data of atmosphere, soil and groundwater at three soil profiles in an arid area of northwest China. Comparison and contrast of the observed and simulated results at the three soil profiles show that there are diurnal vapor entry and outlet fluxes at the dry surface layer (DSL) of 30 cm in the summer season. The vapor entry and re-evaporation account for about 14% of annual precipitation for the heterogeneity soil profile with a mean groundwater depth of 210 cm. Because of limited soil moisture in this arid area, vapor induced re-evaporation occurs shortly in the early daytime. Moreover, the extent of vapor entry, condensation and re-evaporation strong depends on soil properties and water table depths. The lower water table produces the drier soil surface, allowing more vapor entry, condensation and re-evaporation. Whereas the finer grained soil layers benefits the vapor fixation to produce zero fluxes that substantially inhibit the upward liquid water and vapor fluxes, and thereby reduces soil evaporation. The reduced soil evaporation correspondingly decreases the capillary effect on phreatic evaporation, proven by that soil evaporation decreases slowly with decline of water table and the large extinct depth of phreatic evaporation for the finer grained soil profiles. The estimated extinct depth is 180 cm and 200 cm for the soil profiles consisting of silt loam and loamy sand, respectively, much larger than 100 cm of the sandy soil profile. Additionally, as water table is higher and lower than the extinct depth, the models neglecting the vapor - heat function could respectively overestimate and underestimate soil evaporation.

Floods claim a high toll in fatalities and economic impacts. Despite their societal relevance, there is much more to learn about the projected changes in discharge and flooding. Here we force an operational hydrologic model over the state of Iowa with high-resolution convection-permitting climate-model precipitation to evaluate the response of 140 watersheds to climate change. At the end of the century, under the most aggressive scenario in terms of fossil fuel use, we show that the transition from snow to rainfall and a ~30% increase in extreme precipitation rates lead to a doubling of maximum discharge during the spring and extended the flood season into the fall. Total discharge volumes are also expected to increase. Our results suggest that flood projections based on extreme precipitation increases alone substantially underestimate future risk due to the nonlinearity of the hydrologic response explained by long-term soil moisture memory and its feedbacks with precipitation.

Forested watersheds provide many ecosystem services that have become increasingly threatened by wildfire. Stream nitrate (NO 3 -) concentrations often increase following wildfire and can remain elevated for decades. We investigated the drivers of persistent elevated stream NO 3 - in nine watersheds that were burned to varying degrees 16 years prior by the Hayman fire, Colorado, USA. We evaluated the ability of multiple linear regression and spatial stream network modeling approaches to predict observed concentrations of the biologically active solute NO 3 - and the conservative solute sodium (Na +). Specifically, we quantified the degree to which landscape and stream network characteristics predict stream solute concentrations. No landscape variables were strong predictors of stream Na +. Rather, stream Na + variability was largely attributed to flow-connected spatial autocorrelation, indicating that downstream hydrologic transport was the primary driver of spatially distributed Na + concentrations. In contrast, vegetation cover, measured as mean normalized differenced moisture index (NDMI), was the strongest predictor of spatially distributed stream NO 3 - concentrations. Furthermore, stream NO 3 - concentrations had weak flow-connected spatial autocorrelation and high spatial variability. This pattern is likely the result of spatially heterogeneous wildfire behavior that leaves intact forest patches interspersed with high burn severity patches that are dominated by shrubs and grasses. Post-fire vegetation also interacts with watershed structure to influence stream NO 3 - patterns. For example, severely burned convergent hillslopes in headwaters positions were associated with the highest stream NO 3 ‑ concentrations due to the high proportional influence of hillslope water in these locations. Our findings suggest that reforestation is critical for the recovery of stream NO 3 - concentrations to pre-fire levels and targeted planting in severely burned convergent hillslopes in headwater positions will likely have a large impact on stream NO 3 - concentrations.

Ecohydrological investigations commonly use the stable isotopes of water (hydrogen and oxygen) as a conservative ecosystem tracer. This approach requires accessing and analyzing water constrained within plant and soil matrices. Generally, there are six steps that researchers must pass through to retrieve hydrogen and oxygen isotope values from these plant and soil matrices: (i) sampling, (ii) sample storage and transport, (iii) extraction, (iv) pre-analysis processing, (v) isotopic analysis, and (vi) post-processing and correction. At each of these steps cumulative errors can be introduced which sum to non-trivial magnitudes. These errors can impact subsequent interpretations about water cycling through the soil-plant-atmosphere continuum. But these steps in the research ‘process chain’ are just the tip of the iceberg when it comes to uncertainly in published findings. At each of these discreet steps, there are multiple possible options to select from resulting in, as we will show, tens of thousands of possible combinations used by researchers to go from plant and soil samples to isotopic data. In a newly emerging science, so many options can create interpretive confusion and major issues with data comparability. This points to the need for the development of shared standardized approaches. Here we critically examine the state of the process chain, reflecting on the issues associated with each step, and end with suggestions to move our community towards standardization. We hope that critically assessing this common approach will help us see the current problem in its entirety and facilitate community action toward agreed upon standardized approaches.

Using annual water balance analyses may mask intra-annual variability in runoff generation, which could limit our understanding of the similarities and differences between water- and energy-limited catchments. This may be especially limiting in comparisons between catchments close to the threshold between water- and energy-limitation. For this study, we examined runoff generation as a function of catchment storage in four watersheds, with focus on two that exist close to these thresholds to identify how year-to-year variability in storage resulted in intra-annual variations of runoff generation efficiency. Specifically, we focused on one energy-limited catchment in the humid subtropics and one water-limited in a Mediterranean climate. We used measured and calculated daily water balance components to calculate variations in the relative magnitude of daily storage. We isolated precipitation events to draw connections between storage and runoff generation at intra-annual scales and compared our findings to the same metrics in two intensely energy-limited landscapes. We observed distinct stages in daily storage across water years in watersheds at the threshold, where systems experienced wet-up, plateau, and dry-down stages. During the wet-up, precipitation was partitioned to storage, and runoff ratios ( RR) were low. In the plateau, storage was filled, precipitation was partitioned to runoff, causing high RRs. During the dry-down, storage decreased as precipitation was partitioned to evapotranspiration and runoff, causing low RRs. The critical role of evapotranspiration during the growing season resulted in relatively higher RRs during the wet-up than during the dry-down for a given storage value. Thus the same storage amount was partitioned to evapotranspiration or runoff differently throughout the year, depending on the storage stage. Despite their different positions on opposite sides of the threshold, the similarity between the two focus catchments suggests a potential characteristic behavior of systems at the threshold common to both humid and semi-arid landscapes.

Land use affects evapotranspiration rates and is a primary driver of the catchment water balance. The water balance of two catchments in southeastern Australia dominated by either grazed pasture or blue gum ( Eucalyptus globulus) plantation was studied, focusing on the patterns of evapotranspiration (ET) throughout the year. Rainfall, streamflow, and groundwater levels measured between 2015-2019 were combined to estimate annual ET using a water balance equation. In the pasture, eddy covariance was used to measure ET from the catchment. Sap flow measurements were used to estimate tree transpiration in May 2017 – May 2018 and Feb 2019 – Feb 2021 in two different plots within the plantation. The tree transpiration rates were added to direct evaporation, estimated as a percentage of annual rainfall, to calculate ET from the plantation catchment. ET in the pasture showed strong seasonal cycles with very low ET rates in summer and ET rates in spring that were larger than the transpiration rates in the plantation, where trees transpired consistently throughout the year. The estimated annual ET from the water balance equation were comparable to ET estimated from other measurements. In the pasture, ET on average accounted for 88% of annual rainfall, while ET in the plantation was on average 93% of rainfall, exceeding it in the years with annual rainfall lower than 500 mm. The difference between the ET rates in the plantation and the pasture were approximately 30 to 50 mm y­ -1. The larger ET rates in the plantation are reflected in a gradual decrease in the groundwater storage. The differences in ET rates were thus enough to cause a decrease in groundwater storage in the plantation, while the groundwater levels in the pasture remained stable.

Insight into the rainfall-soil moisture (SM) response to land cover is critical for soil hydrological process modeling and management. In this study, five typical land-cover types (forest, shrub, grass, crop, and bare land) and four rainfall patterns (heavy, intermediate, light, and continuous rains) were selected to assess the effects of SM response characteristics on the Loess Plateau of China. We monitored SM at five depths on each land-cover type at 1-h intervals over the growing season of 2019. The results showed that rainfall patterns and land-cover typestogether determined the SM response process and infiltration efficiency. A minimum accumulated rainfall amount of 5 mm was the threshold to trigger a 10-cm SM response. Rain events with higher intensity and smaller sum triggered a quick surface SM response, while larger amounts could percolate deeper and faster. Land-cover change significantly altered the rainfall-SM response dynamics and rainwater utilization efficiency after 20 years of ecological construction. Revegetation sites (mean values of forest, shrub, and grass) increased the soil wetting depth by 14.7%, shortened the SM response time by 27.3%, and accelerated the SM wetting front velocity by 67.2%, which promoted a 35.2% rainfall transformation rate (RTR) across the 1-m profile over all rainfall events (R 1-13). Moreover, planted forest showed the highest RTR of R 1-13 and the maximal increase in soil water storage, which did not aggravate the soil water deficit across the 1-m profile over the growing season. Therefore, we present evidence that planted forests, instead of shrubs, may be beneficial for water conservation if precipitation is greater than 550 mm. The findings of this study prove the role of revegetation on rainwater infiltration capacity and efficiency and can help improve the management of afforestation in arid and semiarid regions.