Ambient Groundwater Flow Research Articles

A Monte‐Carlo Chemical Budget Approach to Assess Ambient Groundwater Flow in Bedrock Open Boreholes

AbstractIn low‐permeability rocks, ambient groundwater flow in open boreholes may go undetected using conventional borehole‐flowmeter tools and alternative approaches may be needed to identify flow. Understanding ambient flow in open boreholes is important for tracking of cross contamination in groundwater. Chlorinated volatile organic compound (CVOC) concentrations from three open boreholes set in a crystalline‐rock aquifer (two of three open boreholes) and a siltstone aquifer (one of three open boreholes) were examined using a new approach and associated software program called the AFCE (Aqueous‐Flow‐Concentration‐Estimator). The program allows comparison of coupled chemical datasets through a Monte‐Carlo simulation and a chemical‐budget approach to assess ambient groundwater flow in open boreholes. The coupled datasets required for the comparison include aqueous CVOC concentrations from groundwater samples from (1) discrete fractures, such as those measured from temporary deployment of straddle‐borehole packer assemblies; and (2) the concentration of the open borehole (wellbore) water, as measured by a vertical profile of passive samplers from within the same open borehole. Because results from the passive samplers represent a composite mixture of the results from the discrete samples under ambient groundwater‐flow conditions, potentially at unknown proportions, the comparison between coupled datasets affords the ability to discern likely water contributions of CVOC from discrete fractures (or fracture zones), and which fractures may be dominating the water chemistry of the open borehole.

Scaling of groundwater flow subject to managed aquifer recharge using injection boreholes: Physical experiments and upscaled numerical models

Managed aquifer recharge (MAR) holds great promise for solving the issue of groundwater resource depletion and quality deterioration in urban areas, but the scaling of groundwater flow patterns under MAR influence from controlled laboratory experiments to complex field scale conditions has not been fully resolved. This study applied physical and numerical groundwater flow experiments at the laboratory scale (100 m) along with ten numerical models at increasing spatial scales (100–103 m) for porous media representative of unconsolidated sand aquifers. We considered several spatial and temporal features that characterize groundwater flow patterns, such as the water table elevation, the thickness of the capillary fringe, the MAR lens geometry (defined as formed by the plume of artificially recharged water), the locations of stagnation points (defined as the points of convergence between ambient groundwater flow and injected water flow), and the groundwater travel times, which were examined at relative scales for comparison. The analysis showed that small-scale sandbox experiments overall underestimate the magnitude of the effects of MAR on groundwater flow pathways. The trends with increasing scales for these parameters were non-linear and could be fitted to power law relationships. These trends were more prominent at scales <20 times the size of the laboratory sand tank. Importantly, the heterogeneity of the sand material, as reflected by non-uniform spatial distributions of hydraulic conductivities (K) also had a major effect on the scaling of groundwater flow patterns, even for the relatively low degree of heterogeneity considered. Specifically, an increase in the average K and its standard deviation further exacerbates the scale dependency. This study provides insights into the appropriate scaling techniques to be used when applying small scale experimental results (e.g., from laboratory) to predict MAR-induced groundwater flow dynamics at larger field scales in moderately heterogeneous, unconsolidated sands.

Experimental investigations on spreading and displacement of fluid plumes around an injection well in a contaminated aquifer

ｄe concentration of the injected fluid decays smoothly around the location of the corresponding sharp interface, which is beyond the theoretical range. This concentration spread should be considered in redefining outer mixing boundaries around injection well between the injected and ambient contaminated solutions. In this study, a radial transport experiment was performed (naphthol green B was first injected as a contaminant plume followed by potassium permanganate as the remediation agent) to mimic the interface mixing in more realistic conditions. The radial advection–dispersion model (e.g., dominated by interfacial driving, R (C/C0 = 0.5), combined with dispersion mixing) was used to obtain the position of the “boundary”, rf, at the mixing fringe of the injected fluid. The results show that the injected oxidant front ran after the contaminated plume at a migration velocity of more than 1.29 times, increasing the mechanical dispersion difference between potassium permanganate and naphthol green B. The increase in the injection rate and ambient groundwater flow velocity improved the mixing of oxidant and contaminated plume (the area difference between the contaminated plume and the oxidant decreased) but also accelerated the contaminated plume coverage. Changing the injection mode (such as linear injection) was an effective method, which enhanced the mixing of the oxidant-contamination plume and weakly expanded the contaminated plume (The swept area, 254 cm2, was the smallest among the three injection modes). The relationship between injection and mixing needs to be further strengthened.

Dipole Flow Tracer Test Models for Identifying Skin Properties and Ambient Groundwater Flow Effects.

As a special single-well test with dual-screened configuration, dipole flow test with a tracer (DFTT) has been an efficient tool for aquifer parameter interpretation. The tool is excellent for characterizing the vertical distribution of hydraulic conductivity, dispersivity and in situ reaction rate in aquifers. It is noteworthy that the effects of skin and ambient groundwater flow (AGF) on the solute transport, which have been widely investigated in terms of conventional tracer tests, have not been comprehensively studied with respect to DFTT. Meanwhile, it is important to elucidate the potentials of DFTT in identifying the primary skin properties as well as the methods to avoid the AGF effects on DFTT. Consequently, we developed a mathematical model for DFTT accounting for skin and AGF. The results show that a negative skin acts as a preferential pathway for tracer transport and the AGF effects on breakthrough curves (BTCs) can be easily mistaken for reaction. Through extensive mathematical model-based experiments, the BTCs calculated from the DFTT models with skin effects were categorized into four types, each of which is able to imply a specific kind of skin. Additionally, a new dimensionless parameter , composed of injection/extraction rate, hydraulic conductivity, hydraulic gradient of AGF and screen parameters, was defined to quantitatively indicate the AGF effects on the overestimation degree of the in situ first-order reaction, implying that one can adjust the field parameters of DFTT to keep less than 3.7% to achieve an overestimation error of reaction rate below 50%.

Numerical simulation of overbank hyporheic transport and biogeochemical reactions in a compound channel

AbstractTransverse hyporheic transport and biogeochemical reactions in a compound channel remain poorly understood. As floodplains are submerged, solutes in the surface water enter the hyporheic zone of the compound channel and react with the solutes upwelling below the floodplain. Aerobic respiration (AR) and denitrification (DN) processes in the hyporheic zone of the compound channel still need to be clarified. In this paper, a 3D hydrodynamic model coupled to a 2D groundwater flow and biogeochemical model was applied to investigate such processes. The model results were verified using laboratory experimental measurements. The effects of bank slope angle, ambient groundwater flow, and ambient groundwater solute concentration on aerobic and anaerobic biogeochemical processes were studied. The denitrification zone was found below the interface between the main channel and the floodplain. The lower bank slope angle favours nitrogen removal. Small losing groundwater discharges promote aerobic respiration and denitrification reactions. As the ambient groundwater oxygen concentration rises, O2 consumption increases while N removal decreases. As the NO3 concentration in ambient groundwater increases, NO3 consumption increases and proportionally NO3 consumption decreases. These findings not only provide a better understanding of biogeochemical reactions in a compound channel but can also be applied to river restoration to efficiently remove nitrate by modifying bank slope angles.

Aerobic respiration in riparian exchange zones of regulated river corridors

AbstractRiver stage fluctuations drive surface water‐groundwater exchanges within river corridors. This study evaluates how repeated daily stage fluctuations, representative of hydropeaking conditions, influence aerobic respiration of river‐sourced dissolved organic carbon (DOC) in the riparian exchange zone using reactive flow and transport simulations. Over 50 hypothetical scenarios were modelled to evaluate how the duration of the daily flood signal, river DOC concentration, aquifer hydraulic conductivity and ambient groundwater flow condition affect the fate and transport of DOC and DO in the riparian aquifer. Time series subsurface snapshots highlight how the various factors influence the subsurface distribution of DOC and DO. The total mass of DOC respired per meter of river had a wide range depending on the parameters, spanning from 1.4 to 71 g over 24‐h, with high hydraulic conductivity and losing ambient groundwater flow conditions favouring the largest amount of DOC respired. The ratio of DOC mass entering the riparian zone with the mass returning to the river showed that as little as 5% to as much as 76% of the DOC that enters the bank during stage fluctuations returns to the river. This return ratio is dependent on river DOC concentration, hydraulic conductivity and ambient groundwater flow. The results illustrate that stage variations due to river regulation can be a significant control on aerobic respiration in riparian exchange zones.

Effects of Decaying Hydraulic Conductivity on the Groundwater Flow Processes in a Managed Aquifer Recharge Area in an Alluvial Fan

Groundwater artificial recharge and medium characteristics represent the major factors in controlling the groundwater flow processes in managed aquifer recharge areas. According to the depositional features of alluvial fans, an analogous homogeneous phreatic sand tank aquifer and the corresponding inhomogeneous scale numerical models were established to investigate the groundwater flow under the combined influence of artificial recharge (human activities) and decaying hydraulic conductivity (medium characteristics). In this study, groundwater flow through a managed aquifer recharge area in an alluvial fan was analyzed under the conditions of decaying hydraulic conductivity (K) with depth or length from apex to apron. The results showed that groundwater flow processes induced by artificial recharge were significantly controlled by the increasing decay exponents of K. The decaying K with depth or length in alluvial fan areas expanded the degree of influence of artificial recharge on groundwater flow. With the increase of decay exponents, the flow directions gradually changed from a horizontal to vertical direction. Groundwater age and spatial variability could also be increased by the increasing decay exponents. The residence time distributions (RTDs) of ambient groundwater and artificially recharged water exhibited logarithmic, exponential, and power law behavior. Penetration depth and travel times of ambient groundwater flow could be affected by artificial recharge and decay exponents. Furthermore, with the increase of decay exponents, the thickness of the artificially recharged water lens and travel times of artificially recharged water were increased. These findings have important implications for the performance of managed aquifer recharge in alluvial fan areas as well as the importance of considering the gradual decrease of K with depth and length.

High-Resolution Integrated Transport Model for Studying Surface Water-Groundwater Interaction.

Transport processes that lead to exchange of mass between surface water and groundwater play a significant role for the ecological functioning of aquatic systems, for hydrological processes and for biogeochemical transformations. In this study, we present a novel integral modeling approach for flow and transport at the sediment-water interface. The model allows us to simultaneously simulate turbulent surface and subsurface flow and transport with the same conceptual approach. For this purpose, a conservative transport equation was implemented to an existing approach that uses an extended version of the Navier-Stokes equations. Based on previous flume studies which investigated the spreading of a dye tracer under neutral, losing and gaining flow conditions the new solver is validated. Tracer distributions of the experiments are in close agreement with the simulations. The simulated flow paths are significantly affected by in- and outflowing groundwater flow. The highest velocities within the sediment are found for losing condition, which leads to shorter residence times compared to neutral and gaining conditions. The largest extent of the hyporheic exchange flow is observed under neutral condition. The new solver can be used for further examinations of cases that are not suitable for the conventional coupled models, for example, if Reynolds numbers are larger than 10. Moreover, results gained with the integral solver provide high-resolution information on pressure and velocity distributions at the rippled streambed, which can be used to improve flow predictions. This includes the extent of hyporheic exchange under varying ambient groundwater flow conditions.

Analysis of the Effects of Dam Release Properties and Ambient Groundwater Flow on Surface Water‐Groundwater Exchange Over a 100‐km‐Long Reach

AbstractHydroelectric dams often create highly dynamic downstream flows that promote surface water‐groundwater (SW‐GW) interactions including bank storage, the temporary storage of river water in the riverbank. Previous research on SW‐GW exchanges in dammed rivers has primarily been at single study sites, which has limited the understanding of how these exchanges evolve as dam releases travel downstream. This study evaluates how dam releases affect SW‐GW exchange continuously over a 100‐km distance. This is accomplished by longitudinally routing water releases through a synthetic river and modeling bed and bank fluid and solute exchange across transverse transects spaced along the reach. Peak and square dam release hydrograph shapes with three magnitudes (0.5, 1.0, and 1.5 m) were considered. The effect of four ambient groundwater flow conditions (very slightly losing, neutral, and two gaining from the perspective of the river) was evaluated for each dam release scenario. Both types of dam release shapes cause SW‐GW interaction over the entire 100‐km distance, and our results show that square type releases cause bank storage exchange well beyond this distance. Strongly gaining conditions reduce the amount of exchange and allow flushing of river‐sourced solute out of the bank after the dam pulse has passed. Both neutral and losing conditions have larger fluid and solute flux into the bank and limit the amount of solute that returns to the river. Our results support that river corridors downstream of dams have increased river‐aquifer connectivity and that this enhanced connectivity can extend at least 100 km downstream.

Assessment of short-term aquifer thermal energy storage for demand-side management perspectives: Experimental and numerical developments

In the context of demand-side management and geothermal energy production, our proposal is to store thermal energy in shallow alluvial aquifers at shorter frequencies than classical seasonal aquifer thermal energy storage. We first conducted a one-week experiment in a shallow alluvial aquifer, which is characterized by a slow ambient groundwater flow, to assess its potential for thermal energy storage and recovery. This experiment has shown that up to 90% of the stored thermal energy can be recovered and would therefore suggest that aquifer thermal energy storage could be considered for demand-side management applications. We then conceptualized, developed, and calibrated a deterministic 3D groundwater flow and heat transport numerical model representing our study site, and we simulated 77 different scenarios to further assess this potential. This has allowed us to demonstrate that low-temperature aquifer thermal energy storage (temperature differences of −4 K for precooling and 3, 6, and 11 K for preheating) is efficient with energy recovery rates ranging from 78 to 87%, in a single aquifer thermal energy storage cycle. High-temperature aquifer thermal energy storage (temperature differences between 35 and 65 K) presents lower energy recovery rates, from 53 to 71%, with all other parameters remaining equals. Energy recovery rates decrease with increasing storage duration and this decrease is faster for higher temperatures. Retrieving directly useful heat (without upgrading with a groundwater heat pump) using only a single storage and recovery cycle appears to be complicated. Nevertheless, there is room for aquifer thermal energy storage optimization in space and time with regard to improving both the energy recovery rates and the recovered absolute temperatures.

Dynamic Hyporheic Zones: Exploring the Role of Peak Flow Events on Bedform‐Induced Hyporheic Exchange

AbstractDischarge varies in space and time, driving hyporheic exchange processes in river corridors that affect biogeochemical cycling and ultimately control the dynamics of biogeochemical hot spots and hot moments. Herein, we use a reduced‐order model to conduct the systematic analysis of the interplay between discharge variability (peak flow intensities, duration, and skewness) and streambed topography (bedform aspect ratios and channel slopes) and their role in the flow and transport characteristics of hyporheic zones (HZs). We use a simple and robust conceptualization of single peak flow events for a series of periodic sinusoidal bedforms. Using the model, we estimate the spatial extent of the HZ, the total amount of exchange, and the residence time of water and solutes within the reactive environment and its duration relative to typical time scales for oxygen consumption (i.e., a measure of the denitrification potential). Our results demonstrate that HZ expansion and contraction is controlled by events yet modulated by ambient groundwater flow. Even though the change in hyporheic exchange flux (%) relative to baseflow conditions is invariant for different values of channel slopes, absolute magnitudes varied substantially. Primarily, peak flow events cause more discharge of older water for the higher aspect ratios (i.e., for dunes and ripples) and lower channel slopes. Variations in residence times during peak flow events lead to the development of larger areas of potential nitrification and denitrification in the HZ for longer durations. These findings have potential implications for river management and restoration, particularly the need for (re)consideration of the importance of hyporheic exchange under dynamic flow conditions.

Analytical solutions for predicting thermal plumes of groundwater heat pump systems

Groundwater heat pump (GWHP) systems have gained attention for space heating and cooling due to their efficiency and low installation costs. Their number is growing in many countries, and therefore in some areas, dense installations are expected. This might lead to thermal interferences between neighbouring groundwater wells and a decrease in efficiency. In the presented study, three analytical formulations are inspected for the prediction of the thermal plume around such open-loop systems under various hydrogeological conditions. A thermal radial transport scenario without background groundwater flow and two advective scenarios with moderate to significant ambient flow velocities (1 and 10 m d−1) are analytically simulated and compared with numerical simulations. Two-dimensional (2D) numerical models are used to estimate the validity of analytical results for a homogeneous confined aquifer, without considering heat transfer in upper and lower layers of the aquifer. In order to represent more realistic aquifer conditions of limited vertical extension, an additional three-dimensional numerical model (3D) is deployed to account for vertical heat losses. The estimated relative errors indicate that the analytical solution of the radial heat transport is in good agreement with both numerical model results. For the advective scenarios, the suitability of the linear and planar advective heat transport models strongly depend on ambient groundwater flow velocity and well injection rate. For low groundwater velocities (1 m d−1), the planar model fits both numerical model results better than the linear advective model. However, the planar model's ability to estimate thermal plumes considerably decreases for high injection rates (>0.6 l s−1). In contrast, the linear advective model shows a good agreement with the two-dimensional numerical results for high groundwater flow conditions (≥10 m d−1). The comparison with the three-dimensional numerical models indicates that the vertical heat transfer is challenging for all of the selected analytical solutions. Despite this, there is a wide range of applicability for the provided analytical solutions in studying the thermal impact of GWHP systems. Hence, the inspected solutions prove to be useful candidates for first-tier impact assessment in crowded areas with potential thermal interferences.

ATES systems in aquifers with high ambient groundwater flow velocity

Aquifer thermal energy storage (ATES) is a technology with worldwide potential to provide sustainable space heating and cooling using groundwater stored at different temperatures. In areas with high ambient groundwater flow velocity (>25 m/y) thermal energy losses by displacement of groundwater may be prevented by application of multiple doublets. In such configurations two or more warm and two or more cold wells are aligned in the direction of the ambient groundwater flow. By controlling the infiltration and extraction rates of the upstream and downstream wells, the advection by ambient groundwater flow can be compensated by storing thermal energy through the upstream well, while re-extracting it from the downstream well.This study uses analytical and numerical tools and a case study to analyze the relevant processes, and provides guidelines for well placement and an operation strategy for ATES wells in aquifers with considerable groundwater flow. The size of the thermal radius relative to ambient groundwater flow velocity is an important metric. With multiple wells to counteract groundwater flow, this ratio affects the pumping scheme of these wells. The optimal distance between them is around 0.4 times the distance traveled by the groundwater in one year. A larger distance negatively affects the efficiency during the first years of operation.

Analysis of the impact of storage conditions on the thermal recovery efficiency of low-temperature ATES systems

Aquifer thermal energy storage (ATES) is a technology with worldwide potential to provide sustainable space heating and cooling using groundwater stored at different temperatures. The thermal recovery efficiency is one of the main parameters that determines the overall energy savings of ATES systems and is affected by storage specifics and site-specific hydrogeological conditions. Although beneficial for the optimization of ATES design, thus far a systematic analysis of how different principal factors affect thermal recovery efficiency is lacking. Therefore, analytical approaches were developed, extended and tested numerically to evaluate how the loss of stored thermal energy by conduction, dispersion and displacement by ambient groundwater flow affect thermal recovery efficiency under different storage conditions. The practical framework provided in this study is valid for the wide range of practical conditions as derived from 331 low-temperature (<25°C) ATES systems in practice.Results show that thermal energy losses from the stored volume by conduction across the boundaries of the stored volume dominate those by dispersion for all practical storage conditions evaluated. In addition to conduction, the displacement of stored thermal volumes by ambient groundwater flow is also an important process controlling the thermal recovery efficiencies of ATES systems. An analytical expression was derived to describe the thermal recovery efficiency as a function of the ratio of the thermal radius of the stored volume over ambient groundwater flow velocity (Rth/u). For the heat losses by conduction, simulation results showed that the thermal recovery efficiency decreases linearly with increasing surface area over volume ratios for the stored volume (A/V), as was confirmed by the derivation of A/V-ratios for previous ATES studies. In the presence of ambient groundwater flow, the simulations showed that for Rth/u<1year, displacement losses dominated conduction losses. Finally, for the optimization of overall thermal recovery efficiency as affected by these two main processes, the optimal design value for the ratio of well screen length over thermal radius (L/Rth) was shown to decrease with increasing ambient flow velocities while the sensitivity for this value increased. While in the absence of ambient flow a relatively broad optimum exists around an L/Rth-ratio of 0.5–3, at 40m/year of ambient groundwater flow the optimal L/Rth-value ranges from 0.25 to 0.75. With the insights from this study, the consideration of storage volumes, the selection of suitable aquifer sections and well screen lengths can be supported in the optimization of ATES systems world-wide.

Flow and Residence Times of Dynamic River Bank Storage and Sinuosity‐Driven Hyporheic Exchange

AbstractHydrologic exchange fluxes (HEFs) vary significantly along river corridors due to spatiotemporal changes in discharge and geomorphology. This variability results in the emergence of biogeochemical hot‐spots and hot‐moments that ultimately control solute and energy transport and ecosystem services from the local to the watershed scales. In this work, we use a reduced‐order model to gain mechanistic understanding of river bank storage and sinuosity‐driven hyporheic exchange induced by transient river discharge. This is the first time that a systematic analysis of both processes is presented and serves as an initial step to propose parsimonious, physics‐based models for better predictions of water quality at the large watershed scale. The effects of channel sinuosity, alluvial valley slope, hydraulic conductivity, and river stage forcing intensity and duration are encapsulated in dimensionless variables that can be easily estimated or constrained. We find that the importance of perturbations in the hyporheic zone's flux, residence times, and geometry is mainly explained by two‐dimensionless variables representing the ratio of the hydraulic time constant of the aquifer and the duration of the event (Γd) and the importance of the ambient groundwater flow (). Our model additionally shows that even systems with small sensitivity, resulting in small changes in the hyporheic zone extent, are characterized by highly variable exchange fluxes and residence times. These findings highlight the importance of including dynamic changes in hyporheic zones for typical HEF models such as the transient storage model.

Denitrification in the banks of fluctuating rivers: The effects of river stage amplitude, sediment hydraulic conductivity and dispersivity, and ambient groundwater flow

AbstractHyporheic exchange induced by periodic river fluctuations leads to important biogeochemical processes, particularly nitrogen cycling, in riparian zones (RZs) where chemically distinct surface water and groundwater mix. We developed a two‐dimensional coupled flow, reactive transport model to study the role of bank storage induced by river fluctuations on removing river‐borne nitrate. Sensitivity analyses were conducted to quantify the effects of river amplitude, sediment hydraulic conductivity and dispersivity, and ambient groundwater flow on nitrate removal rate. The simulations showed that nitrification occurred in the shallower zone adjacent to the bank where oxic river water and groundwater interacted while denitrification occurred deeper into the aquifer and in the riverbed sediments where oxygen was depleted. River fluctuations greatly increased the amount of nitrate being removed; the nitrate removal rate increased as river amplitude increased. Similarly, increasing hydraulic conductivity increased overall nitrate removal since it expanded the denitrifying zone but decreased efficiency. In contrast, increasing sediment dispersivity increased the removal efficiency of nitrate because it promoted mixing between electron acceptors and donors. The presence and direction of ambient groundwater flow had a significant impact on nitrate removal rate when compared to neutral conditions. A losing river showed a larger nitrate removal rate, whereas a gaining river showed a smaller nitrate removal rate. Our results demonstrated that daily river fluctuations created denitrification hot spots within the RZ that would not otherwise exist under naturally neutral or gaining conditions.

Ambient groundwater flow diminishes nitrate processing in the hyporheic zone of streams

AbstractModeling and experimental studies demonstrate that ambient groundwater reduces hyporheic exchange, but the implications of this observation for stream N‐cycling is not yet clear. Here we utilize a simple process‐based model (the Pumping and Streamline Segregation or PASS model) to evaluate N‐cycling over two scales of hyporheic exchange (fluvial ripples and riffle‐pool sequences), ten ambient groundwater and stream flow scenarios (five gaining and losing conditions and two stream discharges), and three biogeochemical settings (identified based on a principal component analysis of previously published measurements in streams throughout the United States). Model‐data comparisons indicate that our model provides realistic estimates for direct denitrification of stream nitrate, but overpredicts nitrification and coupled nitrification‐denitrification. Riffle‐pool sequences are responsible for most of the N‐processing, despite the fact that fluvial ripples generate 3–11 times more hyporheic exchange flux. Across all scenarios, hyporheic exchange flux and the Damköhler Number emerge as primary controls on stream N‐cycling; the former regulates trafficking of nutrients and oxygen across the sediment‐water interface, while the latter quantifies the relative rates of organic carbon mineralization and advective transport in streambed sediments. Vertical groundwater flux modulates both of these master variables in ways that tend to diminish stream N‐cycling. Thus, anthropogenic perturbations of ambient groundwater flows (e.g., by urbanization, agricultural activities, groundwater mining, and/or climate change) may compromise some of the key ecosystem services provided by streams.

Single discharge events increase reactive efficiency of the hyporheic zone

AbstractIn this study, we investigate the impact of single stream discharge events on water exchange, solute transport, and reactions in the hyporheic zone below a natural in‐stream gravel bar. We set up a reactive transport groundwater model with streamflow scenarios that vary by event duration and peak discharge. A steady ambient groundwater flow field is assumed that results in losing, neutral, or gaining stream conditions depending on the stream stage. Across the streambed dissolved oxygen, organic carbon, and nitrate are transported into the subsurface. Additional nitrate is received from upwelling groundwater. Aerobic respiration and denitrification are simulated for scenarios with different stream solute concentrations. Results show that hyporheic exchange flux, solute transport, and consumption increase during events. However, their intensities depend highly on the interplay between event characteristics and ambient groundwater conditions. During events where reversals in the hydraulic gradient occur stream water and solutes infiltrate deeper into the aquifer where they have more time to react. For those events, the reactive efficiency of the hyporheic zone (solute consumption as fraction of influx) for aerobic respiration and denitrification is up to 2.7 and 10 times higher compared to base flow conditions. The fraction of stream nitrate load consumed in the hyporheic zone increases with stream discharge (up to 150 mg/m2/h), but remains below the value under base flow conditions for weak events. Events also increase denitrification of groundwater borne nitrate, but groundwater nitrate flux to the stream decreases by up to 33% due to temporary gradient reversals.

Mixing interfaces, fluxes, residence times and redox conditions of the hyporheic zones induced by dune-like bedforms and ambient groundwater flow

Recent studies highlighted the importance of the interface between streams and their surrounding sediment, known as the hyporheic zone, where stream waters flow through the alluvium. These pore water fluxes stem from the interaction among streambed morphology, stream hydraulics and surrounding groundwater flow. We analytically model the hyporheic hydraulics induced by a spatially uniform ambient groundwater flow made of a horizontal, underflow, and a vertical, basal, component, which mimics gaining and losing stream conditions. The proposed analytical solution allows to investigate the control of simple hydromorphological quantities on the extent, residence time and redox conditions of the hyporheic zone, and the thickness of the mixing interface between hyporheic and groundwater cells. Our analysis shows that the location of the mixing zone shallows or deepens in the sediment as a function of bedform geometry, surface hydraulic and groundwater flow. The point of stagnation, where hyporheic flow velocities vanish and where the separation surface passes through, is shallower than or coincides with the deepest point of the hyporheic zone only due to underflow. An increase of the ambient flow causes a reduction of the hyporheic zone volume similarly in both losing and gaining conditions. The hyporheic residence time is lognormally distributed under neutral, losing and gaining conditions, with the residence time moments depending on the same set of parameters describing dune morphology and stream flow.

Hydraulic controls of in‐stream gravel bar hyporheic exchange and reactions

AbstractHyporheic exchange transports solutes into the subsurface where they can undergo biogeochemical transformations, affecting fluvial water quality and ecology. A three‐dimensional numerical model of a natural in‐stream gravel bar (20 m × 6 m) is presented. Multiple steady state streamflow is simulated with a computational fluid dynamics code that is sequentially coupled to a reactive transport groundwater model via the hydraulic head distribution at the streambed. Ambient groundwater flow is considered by scenarios of neutral, gaining, and losing conditions. The transformation of oxygen, nitrate, and dissolved organic carbon by aerobic respiration and denitrification in the hyporheic zone are modeled, as is the denitrification of groundwater‐borne nitrate when mixed with stream‐sourced carbon. In contrast to fully submerged structures, hyporheic exchange flux decreases with increasing stream discharge, due to decreasing hydraulic head gradients across the partially submerged structure. Hyporheic residence time distributions are skewed in the log‐space with medians of up to 8 h and shift to symmetric distributions with increasing level of submergence. Solute turnover is mainly controlled by residence times and the extent of the hyporheic exchange flow, which defines the potential reaction area. Although streamflow is the primary driver of hyporheic exchange, its impact on hyporheic exchange flux, residence times, and solute turnover is small, as these quantities exponentially decrease under losing and gaining conditions. Hence, highest reaction potential exists under neutral conditions, when the capacity for denitrification in the partially submerged structure can be orders of magnitude higher than in fully submerged structures.