[1] Exchange of water and solutes across the stream-sediment interface is an important control for biogeochemical transformations in the hyporheic zone (HZ). In this paper, we investigate the interplay between turbulent stream flow and HZ flow in pool-riffle streams under various ambient groundwater flow conditions. Streambed pressures, derived from a computational fluid dynamics (CFD) model, are assigned at the top of the groundwater model, and fluxes at the bottom of the groundwater model domain represent losing and gaining conditions. Simulations for different Reynolds numbers (Re) and pool-riffle morphologies are performed. Results show increasing hyporheic exchange flows (m 3 /d) for larger Re and a concurrent decrease in residence time (RT). Losing and gaining conditions were found to significantly affect the hyporheic flow field and diminish its spatial extent as well as rates of hyporheic exchange flow. The fraction of stream water circulating through the hyporheic zone is in the range of 1 Â 10 À5 to 1 Â 10 À6 per meter stream length, decreasing with increasing discharge. Complex distributions of pressure across the streambed cause significant lateral hyporheic flow components with a mean lateral travel distance of 20% of the longitudinal flow paths length. We found that the relationship between pool-riffle height and hyporheic exchange flow is characterized by a threshold in pool-riffle amplitude, beyond which hyporheic exchange flow becomes independent of riffle height. Hyporheic residence time distributions (RTD) are log-normally distributed with medians ranging between 0.7 and 19 h.

At the interface between stream water, groundwater, and the hyporheic zone (HZ), important biogeochemical processes that play a crucial role in fluvial ecology occur. Solutes that infiltrate into the HZ can react with each other and possibly also with upwelling solutes from the groundwater. In this study, we systematically evaluate how variations of gaining and losing conditions, stream discharge, and pool-riffle morphology affect aerobic respiration (AR) and denitrification (DN) in the HZ. For this purpose, a computational fluid dynamics model of stream water flow is coupled to a reactive transport model. Scenarios of variations of the solute concentration in the upwelling groundwater were conducted. Our results show that solute influx, residence time, and the size of reactive zones strongly depend on presence, magnitude, and direction of ambient groundwater flow. High magnitudes of ambient groundwater flow lower AR efficiency by up to 4 times and DN by up to 3 orders of magnitude, compared to neutral conditions. The influence of stream discharge and morphology on the efficiency of AR and DN are minor, in comparison to that of ambient groundwater flow. Different scenarios of O 2 and NO 3 concentrations in the upwelling groundwater reveal that DN efficiency of the HZ is highest under low upwelling magnitudes accompanied with low concentrations of O 2 and NO 3 . Our results demonstrate how ambient groundwater flow influences solute transport, AR, and DN in the HZ. Neglecting groundwater flow in stream-groundwater interactions would lead to a significant overestimation of the efficiency of biogeochemical reactions in fluvial systems.

Hyporheic 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 3 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.

No abstract