Hyporheic Research Articles

Toward a Universal Model of Hyporheic Exchange and Nutrient Cycling in Streams

AbstractIn this paper we demonstrate that several ubiquitous hyporheic exchange mechanisms can be represented simply as a one‐dimensional diffusion process, where the diffusivity decays exponentially with depth into the streambed. Based on a meta‐analysis of 106 previously published laboratory measurements of hyporheic exchange (capturing a range of bed morphologies, hydraulic conditions, streambed properties, and experimental approaches) we find that the reference diffusivity and mixing length‐scale are functions of the permeability Reynolds Number and Schmidt Number. These dimensionless numbers, in turn, can be estimated for a particular stream from the median grain size of the streambed and the stream's depth, slope, and temperature. Application of these results to a seminal study of nitrate removal in 72 headwater streams across the United States, reveals: (a) streams draining urban and agricultural landscapes have a diminished capacity for in‐stream and in‐bed mixing along with smaller subsurface storage zones compared to streams draining reference landscapes; (b) under steady‐state conditions nitrate uptake in the streambed is primarily biologically controlled; and (c) median reaction timescales for nitrate removal in the hyporheic zone are 0.5 and 20 hr for uptake by assimilation and denitrification, respectively. While further research is needed, the simplicity and extensibility of the framework described here should facilitate cross‐disciplinary discussions and inform reach‐scale studies of pollutant fate and transport and their scale‐up to watersheds and beyond.

A Review on Storage Process Models for Improving Water Quality Modeling in Rivers

Water quality is intricately linked to the global water crisis since the availability of safe, clean water is essential for sustaining life and ensuring the well-being of communities worldwide. Pollutants such as industrial chemicals, agricultural runoff, and untreated sewage frequently enter rivers via surface runoff or direct discharges. This study provides an overview of the key mechanisms governing contaminant transport in rivers, with special attention to storage and hyporheic processes. The storage process conceptualizes a ubiquitous reactive boundary between the main channel (mobile zone) and its surrounding slower-flow areas (immobile zone). Research from the last five decades demonstrates the crucial role of storage and hyporheic zones in influencing solute residence time, nutrient cycling, and pollutant degradation. A review of solute transport models highlights significant advancements, including models like the transient storage model (TSM) and multirate mass transport (MRMT) model, which effectively capture complex storage zone dynamics and residence time distributions. However, more widely used models like the classical advection–dispersion equation (ADE) cannot hyporheic exchange, limiting their application in environments with significant storage contributions. Despite these advancements, challenges remain in accurately quantifying the relative contributions of storage zones to solute transport and degradation, especially in smaller streams dominated by hyporheic exchange. Future research should integrate detailed field observations with advanced numerical models to address these gaps and improve water quality predictions across diverse river systems.

Special Seepage Paths Among Nested Groundwater Flow Systems Linking Surface Water Bodies

ABSTRACTA surface water body fed by groundwater is normally known as a terminal place of groundwater flow systems originating from precipitation recharge on highlands. The theory of Tóth predicted that these flow systems form a hierarchically nested structure of groundwater circulation in a composite basin. In this study, we will report new flow paths among groundwater flow systems that were unknown in Tóth's theory, identified as special seepage paths linking different surface water bodies. These seepage paths do not start from the groundwater table but can transmit water between lakes or streams that already serve as discharge zones of traditional local flow systems. As indicated in theoretical models and two real‐world cases, special seepage paths are developed if some parametric conditions are satisfied, especially when surface water bodies cut deeply below the water table or are large enough. Different surface water bodies or different river reaches can directly exchange water, chemicals and heat through deep seepage paths even when both surface and subsurface water divides exist between them. Special seepage paths may play a role in the regional scale hyporheic flow or contribute to inter‐basin groundwater flow. The knowledge of special seepage paths could greatly improve our conventional perception of surface water‐groundwater interaction, groundwater age and geochemical and heat transport at the river basin scale.

Groundwater Geochemistry in the Karst-Fissure Aquifer System of the Qinglian River Basin, China

The Qinglian River plays a significant role in China’s national water conservation security patterns. To clarify the relationship between hydrogeochemical properties and groundwater quality in this karst-fissure aquifer system, drilling data, hydrochemical parameters, and δ2H and δ18O values of groundwater were analyzed. Multiple indications (Piper diagram, Gibbs diagram, Na+-normalized molar ratio diagram, chloro-alkaline index 1, mineral saturation index, and principal component analysis) were used to identify the primary sources of chemicals in the groundwater. Silicate weathering, oxidation of pyrite and chlorite, cation exchange reactions, and precipitation are the primary sources of dissolved chemicals in the igneous-fissure water. The most relevant parameters in the karst water are possibly from anthropogenic activities, and other chemicals are mostly derived from the dissolution of calcite and dolomite and cation exchange reactions. Notably, the chemical composition of the deep karst water from the karst basin is mainly influenced by the weathering of carbonate and cation exchange reactions and is less affected by human activities. The hydrogeochemical properties of groundwater in the karst hyporheic zone are influenced by the dissolution of carbonates and silicates, evaporation, and the promotion effect of dissolution of anorthite or Ca-containing minerals. Moreover, the smallest slope of the groundwater line from the karst hyporheic zone among all groundwater groups revealed that the mixing effects of evaporation, isotope exchange in water–rock interaction or deep groundwater recharge in the karst hyporheic zone are the strongest. The methods used in this study contribute to an improved understanding of the hydrogeochemical processes that occur in karst-fissure water systems and can be useful in zoning management and decision-making for groundwater resources.

Hydrological controls of a riparian wetland based on stable isotope data and model simulations.

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ2H and δ18O) measurement, as well as pH and oxidation-reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2-4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ2H and δ18O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. A monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

Investigation of the hyporheic zone of two gravel-bed rivers after reservoir draining

Investigation of the hyporheic zone of two gravel-bed rivers after reservoir draining

Spatio-temporal and depth-oriented evaluation of hyporheic zone processes in headwater catchments

The hyporheic zone (HZ) is vital for nutrient cycling and overall stream ecosystem health across groundwater-surface water interfaces, with biogeochemical and physical variations influencing exchange processes. To gain a depth-oriented insight into the various hyporheic functioning, isotopic (18O and 2H) and chemical analyses, along with physical parameters, were performed in two streams named Ahna and Losse in North Hesse, Germany. Multi-level interstitial probes were installed to extract subsurface pore water samples up to 0.45 m depth. We identified three distinct HZ in Ahna upstream, downstream and Losse. The Ahna upstream HZ was a less permeable barrier with low water fluxes but reactive, while the downstream HZ was less reactive with higher permeability and vertical fluxes. The Losse HZ, remained non-reactive but well-mixed with good hydraulic circulation due to fast downward flows. Isotopic composition revealed distinctive features, indicating increased evaporation and intricate mixing with other water sources and intensified subsurface flow in Ahna downstream compared to upstream. Hydro-chemical dynamics in Ahna exposed increased ion concentrations downstream, driven by anthropogenic factors along the course. Upstream denitrification was dominant under low oxygen conditions due to the slow infiltration and low mixing rates. The primary water source for both streams and HZs was groundwater, with the potential for rainwater to serve as a secondary source. This study shows the spatial heterogeneity of HZ processes and highlights that the specific HZ conditions need to be examined locally.

The impact of organic matter and iron/calcium coupling on phosphorus retention in the hyporheic zone of the Danjiangkou area tributary: Evidence from bonding recognition

The coupling between organic matter (OM) and minerals considerably influences the phosphorus (P) cycle within the hyporheic zone, but the role of different geological mineral-organic complexes (MOCs) on P burial during hyporheic exchange remains under-explored. This study investigates the effects of OM and iron (Fe)/calcium (Ca) coupling on P migration within the hyporheic zone of an agricultural tributary to the Danjiangkou Reservoir. These relationships were explored by measuring hyporheic flow (q), organic and inorganic P forms, and sediment PO4-P adsorption capacity [following treatment with fulvic acid (FA), Fe-OM, or Ca-OM]. Multivariate statistical analysis, X-Ray Diffraction, Fourier-transform Infrared Spectroscopy, and X-ray Photoelectron Spectroscopy were employed to elucidate the underlying mechanisms. Results indicate that upward hyporheic flow transports dissolved porewater P into surface water, contributing 11.27–12.13 % of the total P flux. MOCs associated with Fe(III)/Ca silicate minerals, along with FA and labile OM, were identified as key OM fractions influencing P migration, contributing 5–24 %, 10–11.7 %, and 6–14.9 % to the overall flux, respectively. FA and labile OM facilitate P release, whereas MOCs enhance P retention. Ca-OM is the most efficient PO4-P adsorption [adsorption capacity (AC): 0.8980–0.9524 mg/g], followed by Fe-OM (AC: 0.5120–0.7020 mg/g), original sediment (AC: 0.4368–0.5596 mg/g), and FA (AC: 0.2657–0.2769 mg/g). Cation bridges, primarily formed by -OH and -NH2 groups within Ca-OM (outer-sphere complexes), promote greater P adsorption than Fe-OM (inner-sphere complexes, mainly associated with -COOH). However, Fe-OM-P exhibits a more stable structure. In high P environments, P adsorption onto Ca-OM may induce the release of labile OM, temporarily retaining P through resorption onto labile OM. Hyporheic flow with higher pH and Eh values promotes MOC formation, underscoring their significant P retention capacity. Therefore, strategic MOC use within the hyporheic zone is crucial for mitigating surface water eutrophication.

Ecologically friendly remediation of groundwater sulfamethoxazole contamination: Biologically synergistic degradation by thermally modified activated carbon-activated peracetic acid in porous media

This study examines the efficacy and environmental impact of peracetic acid (PAA) activated by thermally modified activated carbon (AC600) for degrading antibiotics in actual groundwater. Laboratory-scale experiments evaluated the system's effects on contaminant degradation, ecological balance, and substance cycling in the hyporheic zone. Our findings demonstrated the effectiveness of the AC600/PAA system in removing sulfamethoxazole (SMX) from groundwater porous media. Additionally, the AC600/PAA system synergistically interacted with the in-situ microbiota of the hyporheic zone, producing more fragmented degradation products without increasing mixed toxicity. Bacterial abundance increased post-reaction, with notable alterations in the bacterial community and enhanced bacterial metabolism. Key genera such as Lysobacter thrived in the treated environment, playing critical roles in microbiota modification and SMX degradation. The pH remained stable before and after the reaction, while dissolved organic carbon content increased. Overall, our results highlight the promising potential of PAA activation by carbonaceous materials as a low-impact, ecologically friendly technology for in-situ remediation of organic pollutants in groundwater, characterized by high compatibility and biosynthesis.

Interaction and bacterial effects of microplastics pollution on heavy metals in hyporheic sediments of different land-use types in the Beiluo River Basin

Microplastics (MPs) pollution is ubiquitous, causing serious ecological damage by threatening the growth and health of living organisms. This study investigated the vertical and horizontal distribution of MPs, MPs-heavy metals (MPs-HMs) accumulation, contamination assessment and microbial biodiversity in hyporheic sediments of different land-use types. MPs abundance in shallow sediments (0-30 cm) was significantly higher than that in deep sediments (30-60 cm), with fewer large MPs in the deep sediments. Blue, fiber, and <500 μm were the dominant MPs types, and polystyrene, polylactic acid, and polyvinyl chloride were the dominant polymers in the Beiluo River Basin. The average concentrations of HMs detected in MPs were all much higher than the same metals in the sediments. The pollution loading index of MPs was higher in areas with a greater proportion of anthropogenic land use, and MP-HM were present to varying degrees in the vertical distribution (PN>1). Critically, bacterial diversity of anthropogenic land use was smaller than that of natural land use. High MP-HM concentrations reduced the abundance of cyanobacteria, nitrospirota, acidobacteriota, and planctomycetota, whereas desulfobacterota, chloroflexi, myxococcota, actinobacteriota, and proteobacteria have developed tolerance to MP-HM. Overall, our findings contribute to the understanding of the relationship between different land-use types and the spatial distribution of MPs and MP-HM, which is critical to manage and mitigate the hyporheic zone pollution.

Numerical Simulation and Verification of Water and Heat Transport in the Lateral Hyporheic Zone of River-Groundwater under Water Exchange Conditions

Numerical Simulation and Verification of Water and Heat Transport in the Lateral Hyporheic Zone of River-Groundwater under Water Exchange Conditions

Large eddy simulations of zinc ions transfer to turbulent flows from hyporheic zone

Large eddy simulations of zinc ions transfer to turbulent flows from hyporheic zone

Longitudinal Dispersion and Hyporheic Exchange of Neutrally Buoyant Microplastics in the Presence of Waves and Currents

An experimental study was conducted to identify the behaviour of neutrally buoyant microplastics (specific density, 0.94) in different hydrodynamic conditions while focusing on combined wave–current conditions and the mixing across the hyporheic zone. For in-water-column microplastics, it was observed that the streamwise dispersion of neutrally buoyant microplastics is comparable to solute dye in both slow open-channel flow conditions and combined wave–current conditions. However, for in-bed microplastics, when compared to soluble tracers, the longer timespans associated with the hyporheic exchange process allowed the density effects to enhance the vertical exchange when compared to solutes.

Mechanistic definition and prediction of the mass exchange coefficient between rivers and hyporheic zones: The [formula omitted] of two [formula omitted]s

Solute transport in interconnected rivers and hyporheic zones is typically modeled through dual-domain models where first-order solute mass transfer between the two domains, ΩR and ΩHZ, is represented by a coefficient α. The transient storage model (TSM) is an example of such an approach. In practice, α is determined by fitting the tails of solute tracer breakthrough curves using a TSM. This approach has led to ambiguity regarding α’s physical meaning and transferability. Here, we investigated the physical basis for α and tested it with virtual experiments through the fully coupled multiphysics model hyporheicFoam for the ΩR−ΩHZ system. hyporheicFoam explicitly simulated coupled flow and solute transport over a kilometer with centimeter-scale resolution. Model results were analyzed to calculate α following its theoretical definition directly. Using the determined α within a TSM enables accurate reproduction of solute transport, underscoring α’s physical relevance and precision.

Uncovering the hidden world of riverbed sediments: The role of sediment heterogeneity and cross-bar channel fills in the hydrogeochemical dynamics of the hyporheic zone

Groundwater-surface water interaction (hyporheic exchange) is critical in numerous hydrogeochemical processes; however, hyporheic exchange is difficult to characterize due to the various spatial (e.g., sedimentary architecture) and temporal (e.g., stage fluctuations) variables that influence it. This interdisciplinary study brings forth novel insights by integrating various methodologies including geophysical surveys, physical and chemical sediment characterization, and water chemistry analysis to explore the interplay of the numerous facets governing hyporheic zone processes within a compound bar deposit. The findings reveal distinct sedimentary facies and geochemical zones within the compound bar, driven by the sedimentary architecture. Cross-bar channel fills are identified as critical structures influencing hydrogeochemical dynamics, acting as baffles to groundwater flow and modulating nutrient transformations. Geophysical imaging and hydrogeochemical analyses highlight the complex interplay between sediment characteristics and subsurface hydraulic connectivity, emphasizing the role of sediment heterogeneity in controlling hyporheic exchange and solute mixing. The study concludes that sediment heterogeneity, particularly the presence of cross-bar channel fills, plays a pivotal role in the hydrogeochemical dynamics of the hyporheic zone. These structures significantly influence hyporheic flow paths, solute residence times, and nutrient cycling, underscoring the necessity to consider the fine-scale sedimentary architecture in models of hyporheic exchange. The findings contribute to a deeper understanding of riverine ecosystem processes, offering insights that can inform management strategies for water quality and ecological integrity.

On the role of the Froude number on flow, turbulence, and hyporheic exchange in open-channel flow through boulder arrays

In this paper, the results of numerical simulations of open-channel flow through boulder arrays at varying Froude numbers are reported. The simulations aim at clarifying the role of the Froude number on flow, turbulence, and hyporheic exchange. At low and intermediate Fr, the boulder top is above the water surface, and time-averaged streamwise flow velocity, Reynolds shear stresses, and the turbulent kinetic energy (TKE) are relatively low in the wake of boulders. Conversely, at high Fr values, the boulders are submerged, hence the flow separates at the boulder crest, creates vertical recirculation, and reattaches on the bed downstream, resulting in an area of elevated Reynolds shear stresses and TKE downstream of the boulders. Two dominant turbulence structures are observed: (i) flapping of boulder wakes with a characteristic length of 2.1 times the boulder diameter (D) at low and intermediate Fr and (ii) an upstream oriented hairpin vortex with a length scale of 1.0D at high Fr. These turbulence structures influence hyporheic exchange downstream of boulders within a limited region of x/D&amp;lt;2.0. In other locations, hyporheic flow is driven by downwelling flow immediately upstream of boulders with a wavelength larger than 2.9D. Finally, the normalized time-averaged hyporheic flux increases with increasing Fr, but it decreases at higher Fr values once the overtopping flow disrupts the formation of the boulder wake.

Hotspots of Dissolved Arsenic Generated from Buried Silt Layers along Fluctuating Rivers.

Previous studies along the banks of the tidal Meghna River of the Ganges-Brahmaputra-Meghna Delta demonstrated the active sequestration of dissolved arsenic (As) on newly formed iron oxide minerals (Fe(III)-oxides) within riverbank sands. The sand with high solid-phase As (>500 mg/kg) was located within the intertidal zone where robust mixing occurs with oxygen-rich river water. Here we present new evidence that upwelling groundwater through a buried silt layer generates the dissolved products of reductive dissolution of Fe(III)-oxides, including As, while mobilization of DOC by upwelling groundwater prevents their reconstitution in the intertidal zone by lowering the redox state. A three end-member conservative mixing model demonstrated mixing between riverbank groundwater above the silt layer, upwelling groundwater through the silt layer, and river water. An electrochemical mass balance model confirmed that Fe(III)-oxides were the primary electron acceptor driving the oxidation of DOC sourced from sediment organic carbon in the silt. Thus, the presence of an intercalating silt layer in the riverbanks of tidal rivers can represent a biogeochemical hotspot of As release while preventing its retention in the hyporheic zone.

Multiple climate change stressors reduce the emergence success of gravel-spawning fish species and alter temporal emergence patterns

Climate change, with its profound effects on stream sediment, hydrological, and temperature dynamics, will exacerbate impacts on habitat conditions for many species, particularly those with vulnerable early life stages relying on the hyporheic zone, such as gravel-spawning fishes. Due to the complex and interactive nature of multiple stressor effects, we employed large-scale outdoor mesocosms to systemically test how the reproductive success of three gravel-spawning fish species brown trout (Salmo trutta), nase, (Chrondrostoma nasus) and Danube salmon (Hucho hucho) was affected by individual and combined effects of warming (+3–4 °C), fine sediment (increase in <0.85 mm by 22 %) and low-flow (eightfold discharge-reduction). Fine sediment had the most detrimental effect on emergence rate and fry length in all three species, reducing the emergence rate to zero in brown trout, 9 % in nase, and 4 % in Danube salmon. The emergence mortality caused by fine sediment surpassed that of hatching distinctly, suggesting that negative effects due to hypoxia were considerably exacerbated by entombment. Warming had only minor effects as a single stressor, but low flow reduced emergence rates of the spring spawning species nase and Danube salmon by 8 and 50 %, respectively. In combined treatments including fine sediment, however, the emergence success of all three species responded strongly negatively, even in the cyprinid species nase, which showed little interactive effects between stressors regarding hatching success. Warming and fine sediment also led to the earlier emergence of fry, implying a risk of asynchrony with available food resources. This study dramatically shows that climate change can have deleterious impacts on the reproductive success of gravel-spawning fish species, irrespective of taxonomic or ecological traits.

Effects of repeated drawdown flushing on riverbed fine sediment dynamics downstream from a dam

Sediment accumulation in reservoirs is frequently a problem, impelling dam managers to implement strategies such as drawdown flushing to limit siltation. Drawdown and other sediment removal methods may induce riverbed clogging downstream of dams, especially in river sections where water is diverted, thereby reducing transport capacity. In this study, we investigate the effects of drawdown flushing downstream of the Plan d’Arem run-of-the-river dam (Upper Garonne, Spain – France border) using an adaptive inter-site comparison strategy to consider a range of spatial and temporal conditions, which allow us to separate the effects of dam storage and flushing from other potential factors. Drawdown flushing has been undertaken three times during the study period over a short span of time (2 months). We couple bed material sampling, which provides direct information on bed composition, with airborne infrared thermal imaging to better interpret whether fine sediment interstitial storage within the bed is associated with clogging. We also measure bed surface grain size and bed mobility in order to investigate their potential role in controlling fine sediment dynamics. We identify surface grain size and water diversion as the main factors controlling fine sediment spatial distribution, with coarse-grained bed-surfaces and by-passing promoting fine sediment enrichment. As a result, sites located upstream and within the by-passed reached of the Plan d’Arem dam show higher fine sediment interstitial storage than sites downstream from the by-passed reach. Results from thermal imagery demonstrate such interstitial storage does not induce clogging effect. The reaches with the most important sediment storage host a high number of cool-water patches, indicating water exchanges with the hyporheic zone. Post-flushing bed composition indicates systematic export of fine sediments from the bed matrix at all sites after the first operation, afterwards fine levels remain low after the second and the third flushing. The low impact of dam flushing in term of clogging is interpreted by the fact that the interstitial material is very sandy.

Numerical simulation of groundwater in hyporheic zone with coupled parameter stochastic scheme

Groundwater numerical modeling is a crucial scientific tool for understanding groundwater circulation and supporting regional water resource planning and management. The effectiveness of these models depends largely on the accuracy of hydrogeological parameters within aquifers, which are often spatially heterogeneous and randomly distributed due to complex geological and tectonic factors. Traditional modeling approaches frequently overlook this randomness, compromising the precision and resolution of groundwater simulations. This study focuses on a section of the Qingshui River in the Huaihe River Basin. Using field and laboratory data, probability distribution functions for key parameters like hydraulic conductivity, specific yield, and specific storage were developed. These functions were integrated into the groundwater model to reflect the inherent stochastic nature of aquifer properties. This integration significantly enhanced model accuracy, reducing the root mean square error of simulated water levels from 0.47–1.43 m to 0.13–0.16 m and improving the Nash-Sutcliffe efficiency coefficients (NSE) from −2.96–0.73 to 0.94–0.98. Additionally, the model facilitated analysis of the interactions between river and groundwater, particularly in the hyporheic zone, under various scenarios. It identified spatial and temporal variations in groundwater recharge dynamics and delay effects at different distances from the river channel. For instance, recharge rates at 50 m and 150 m from the river were 0.295 m/day and 0.015 m/day, respectively, indicating stronger recharge closer to the river. The study also assessed the impact of varying river flows, riverbed permeability, and irrigation practices on water exchanges between the river and groundwater. These factors were found to significantly influence the intensity of water exchange, seepage, and groundwater reserves. This research provides valuable insights for managing river-groundwater interactions and analyzing the ecological environment of surrounding groundwater systems, underscoring the importance of incorporating stochastic characteristics into groundwater modeling.