The hyporheic zone is the ecotone between stream and river channel flow and groundwater that can process nutrients and improve water quality. Transient hyporheic zones occur in the riparian zone (bank storage or “lung model” exchange) during channel stage fluctuations. Recent studies show that soil pipes are widespread in stream banks and beneath floodplains, creating highly preferential flow between channel and riparian groundwater such that the traditional Darcy model of flow does not apply. We used MODFLOW with the conduit flow package to model a series of stream bank soil pipes and examined soil pipe density (number per m), length, diameter, height above baseflow water surface, connectivity, and matrix hydraulic conductivity on transient particle flow paths and total hyporheic exchange volume (i.e., bank storage) over the course of a peak flow (e.g., storm) event. We found that adding five soil pipes per meter more than doubled hyporheic volume. Soil pipe length was the most important control; adding one 1.5‐m‐long soil pipe caused a 73.4% increase in hyporheic volume. The effect of increasing soil pipe diameter on hyporheic volume leveled off at ~1 cm, as flow limitation switched from pipe flow to pipe‐matrix exchange. To validate our approach, we used the model to successfully reproduce trends from field studies. Our results highlight the need to consider soil pipes when modeling, monitoring, or managing bank storage, floodplain connectivity, or hyporheic exchange.
Peak flow events in gaining stream/river channels cause lung model hyporheic exchange with the banks (bank storage), which fosters beneficial reactions as polluted channel water cycles through riparian groundwater. Soil pipes are common along stream/riverbanks, and enhance exchange, yet their effect on reactions such as denitrification is unknown. We used MODFLOW with the Conduit Flow Package to simulate lung model exchange during a peak flow event in a section of stream bank/riparian soil with soil pipes, and MT3D‐USGS to estimate nitrate transport and denitrification. We varied soil matrix hydraulic conductivity (
- NSF-PAR ID:
- 10368407
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Water Resources Research
- Volume:
- 58
- Issue:
- 4
- ISSN:
- 0043-1397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Dam removals are on the increase across the US with Pennsylvania currently leading the nation. While most dam removals are driven by aquatic habitat and public safety considerations, we know little about how dam removals impact water quality and riparian zone processes. Dam removals decrease the stream base level, which results in dewatering of the riparian zone. We hypothesized that this dewatering of the riparian zone would increase nitrification and decrease denitrification, and thus result in nitrogen (N) leakage from riparian zones. This hypothesis was tested for a 1.5 m high milldam removal. Stream, soil water, and groundwater N concentrations were monitored over 2 years. Soil N concentrations and process rates and
δ 15N values were also determined. Denitrification rates and soilδ 15N values in riparian sediments decreased supporting our hypothesis but no significant changes in nitrification were observed. While surficial soil water nitrate‐N concentrations were high (median 4.5 mg N L−1), riparian groundwater nitrate‐N values were low (median 0.09 mg N L−1), indicating that nitrate‐N leakage was minimal. We attribute the low groundwater nitrate‐N to denitrification losses at the lower, more dynamic, groundwater interface and/or dissimilatory nitrate reduction to ammonium (DNRA). Stream water nitrate‐N concentrations were high (median 7.6 mg N L−1) and contrary to our dam‐removal hypothesis displayed a watershed‐wide decline that was attributed to regional hydrologic changes. This study provided important first insights on how dam removals could affect N cycle processes in riparian zones and its implications for water quality and watershed management. -
more » « less
Abstract Log jams alter gradients in hydraulic head, increase the area available for hyporheic exchange by creating backwater areas, and lead to the formation of multiple channel branches and bars that drive additional exchange. Here, we numerically simulated stream‐groundwater interactions for two constructed flume systems—one without jams and one with a series of three jams—to quantify the effects of interacting jam structures and channel branches on hyporheic exchange at three stream flow rates. In simulations without jams, average hyporheic exchange rates ranged from 2.1 × 10−4to 2.9 × 10−4 m/s for various stream discharge scenarios, but with jams, exchange rates increased to a range of 1.3 × 10−3–3.5 × 10−3 m/s. Largely due to these increased hyporheic exchange rates, jams increased stream‐groundwater connectivity or decreased the turnover length that stream water travels before interacting with the hyporheic zone, by an order of magnitude, and drove long flow paths that connected multiple jams and channel threads. Decreased turnover lengths corresponded with greater reaction significance per km, a measure of the potential for the hyporheic zone to influence stream water chemistry. For low‐flow conditions, log jams increased reaction significance per km five‐fold, from 0.07 to 0.35. Jams with larger volumes led to longer hyporheic residence times and path lengths that exhibited multiple scales of exchange. Additionally, the longest flow paths connecting multiple jams occurred in the reach with multiple channel branches. These findings suggest that large gains in hydrologic connectivity can be achieved by promoting in‐stream wood accumulation and the natural formation of both jams and branching channels.
-
more » « less
Abstract In streams where water temperatures stress native biota, management of riparian shade or hyporheic exchange are both considered viable management strategies for reducing the peaks of daily and seasonal stream channel temperature cycles. Although shade and hyporheic exchange may have similar effects on stream temperatures, their mechanisms differ. Improved understanding of the heat‐exchange mechanisms influenced by shade and hyporheic exchange will aid in the appropriate application of either stream temperature management strategy. To illustrate a conceptual model highlighting shade as ‘thermal insulation’ and hyporheic exchange imparting ‘thermal capacitance’ to a stream reach, we conducted an in‐silico simulation modelling experiment increasing shade or hyporheic exchange parameters on an idealized, hypothetical stream. We assessed the potential effects of increasing shade or hyporheic exchange on a stream reach using an established process‐based heat‐energy budget model of stream‐atmosphere heat exchange and incorporated an advection‐driven hyporheic heat exchange routine. The model tracked heat transport through the hyporheic zone and exchange with the stream channel, while including the effects of hyporheic water age distribution on upwelling hyporheic temperatures. Results showed that shade and hyporheic exchange similarly damped diurnal temperature cycles and differentially altered seasonal cycles of our theoretical stream. In winter, hyporheic exchange warmed simulated channel temperatures whereas shade had little effect. In summer, both shade and hyporheic exchange cooled channel temperatures, though the effects of shade were more pronounced. Our simple‐to‐grasp analogies of ‘thermal insulation’ for shade effects and ‘thermal capacitance’ for hyporheic exchange effects on stream temperature encourage more accurate conceptualization of complex, dynamic heat exchange processes among the atmosphere, stream channel, and alluvial aquifer.
-
more » « less
Abstract Nitrous oxide (N2O) is a potent ozone‐depleting greenhouse gas produced by incomplete denitrification. Recent works on riverine N2O emissions focus mainly on contributions from in‐channel, benthic, and fluvial hyporheic environments under assumptions of steady‐state conditions and homogeneous sediment hydraulic conductivity (
K ). However, riparian floodplains are also a potentially important N2O source characterized by complex sediment heterogeneity and dynamic surface and groundwater interactions. We use numerical flow and reactive transport models to investigate the influence of complex sedimentary architecture and high‐flow events (e.g., storms) on N2O production. We interpret the correlation between flow and solute fields with the flow topological Okubo‐Weiss metric (OW) and the scalar dissipation rate weighted by soil organic matter (OM) fraction and soil saturation. We model a heterogeneous riparian floodplain based on field observations from the Theis Environmental Monitoring and Modeling Site, Ohio, USA. N2O production is greatest within intermediate‐K sediments (e.g., sands) where denitrification rates are highest, and emissions increase by more than an order of magnitude during storms. Sensitivity analysis reveals that the denitrification rate is most influential for N2O flux, accounting for nearly 46% of the variance in production rates. Denitrification rates adapt to spatial changes in the flow topology (measured by OW) related to sediment heterogeneity and are strongly influenced by subsurface mixing dynamics. Mixing is greatest in shear strain‐dominated regions, while vorticity promotes OM dissolution and prolongs residence times. Accurate lithologic representation is imperative for characterizing subsurface N2O production dynamics, especially given growing concern regarding climate change driven hydrologic changes within watersheds worldwide.
