Hyporheic Research Articles

Spatial distribution of antibiotic resistance genes of the Zaohe-Weihe Rivers, China: exerting a bottleneck in the hyporheic zone.

The hyporheic zone (HZ) is an active biogeochemical region where groundwater and surface water mix and a potential reservoir for antibiotic resistance genes (ARGs). In this paper, the relative abundance and spatial distribution of ARGs in the HZ media were investigated, taking into consideration both the five speciation of six metals and the local characteristics. The samples of surface water, groundwater, and sediment were collected from Zaohe-Weihe Rivers of Xi'an City, which is a representative city with characteristics of the northwest region of China. Of 271 ARGs associated with 9 antibiotics, 228ARGs were detected, with a total detection rate of 84%. Sulfonamide and aminoglycoside ARGs were the dominant types of ARGs. The top 6 ARGs and mobile genetic elements (MGEs) in terms of abundance were tnpA-04, cepA, sul1, aadA2-03, sul2and intI1. The results of principal component analysis (PCA) showed that the distribution characteristics of ARGs were not associated with the sampling sites but with the environmental medias. Similarity in the water phases and significant differences in the water and sediment phases were found. The redundancy analysis (RDA) identified the key factors controlling ARG pollution, including dissolved oxygen (DO) in surface water, total nitrogen (TN) in groundwater, and total organic carbon (TOC) in sediment. In terms of the speciation of heavy metals, we further revealed the promotion effect between ARGs and heavy metals, especially the residual fraction of Ni. In terms of horizontal transfer mechanism, ARGs were significantly correlated with tnpA-03 in water phase and tnpA-04 in sediment. In the three media, intI1 and ARGs all show a significant correlation. These findings showed that hyporheic zone exerted a bottleneck effect on the distribution and transfer of ARGs.

Hyporheic hydraulic geometry: Conceptualizing relationships among hyporheic exchange, storage, and water age.

Hyporheic exchange is now widely acknowledged as a key driver of ecosystem processes in many streams. Yet stream ecologists have been slow to adopt nuanced hydrologic frameworks developed and applied by engineers and hydrologists to describe the relationship between water storage, water age, and water balance in finite hydrosystems such as hyporheic zones. Here, in the context of hyporheic hydrology, we summarize a well-established mathematical framework useful for describing hyporheic hydrology, while also applying the framework heuristically to visualize the relationships between water age, rates of hyporheic exchange, and water volume within hyporheic zones. Building on this heuristic application, we discuss how improved accuracy in the conceptualization of hyporheic exchange can yield a deeper understanding of the role of the hyporheic zone in stream ecosystems. Although the equations presented here have been well-described for decades, our aim is to make the mathematical basis as accessible as possible and to encourage broader understanding among aquatic ecologists of the implications of tailed age distributions commonly observed in water discharged from and stored within hyporheic zones. Our quantitative description of "hyporheic hydraulic geometry," associated visualizations, and discussion offer a nuanced and realistic understanding of hyporheic hydrology to aid in considering hyporheic exchange in the context of river and stream ecosystem science and management.

Production of hydrogen peroxide in an intra-meander hyporheic zone at East River, Colorado

The traditionally held assumption that photo-dependent processes are the predominant source of H2O2 in natural waters has been recently questioned by an increrasing body of evidence showing the ubiquitiousness of H2O2 in dark water bodies and in groundwater. In this study, we conducted field measurement of H2O2 in an intra-meander hyporheic zone and in surface water at East River, CO. On-site detection using a sensitive chemiluminescence method suggests H2O2 concentrations in groundwater ranging from 6 nM (at the most reduced region) to ~ 80 nM (in a locally oxygen-rich area) along the intra-meander transect with a maxima of 186 nM detected in the surface water in an early afternoon, lagging the maximum solar irradiance by ∼ 1.5 h. Our results suggest that the dark profile of H2O2 in the hyporheic zone is closely correlated to local redox gradients, indicating that interactions between various redox sensitive elements could play an essential role. Due to its transient nature, the widespread presence of H2O2 in the hyporheic zone indicates the existence of a sustained balance between H2O2 production and consumption, which potentially involves a relatively rapid succession of various biogeochemically important processes (such as organic matter turnover, metal cycling and contaminant mobilization). More importantly, this study confirmed the occurrence of reactive oxygen species at a subsurface redox transition zone and further support our understanding of redox boundaries on reactive oxygen species generation and as key locations of biogeochemical activity.

Local flow convergence, bed scour, and aquatic habitat formation during floods around wooden training structures placed on sand-gravel bars

Riverine ponds, which are formed and sustained through sediment erosion and deposition, are key habitats for enhancing biodiversity in river reaches. The objective of this study was to understand the roles of traditional river-training wooden structures called “seigyu” on the formation of ponds on nonvegetated bars. Here, the spatial and temporal patterns of the flow and bedform coupled with pond formations for several flood events were assessed. The surface flow patterns were monitored by an unmanned aerial vehicle (UAV) and evaluated by large-scale particle image velocimetry (LSPIV); the maximum flow velocities were 1.3 m s−1 and 1.9 m s−1 during floods when seigyu units were partially and fully submerged, respectively. Although the overall mean flow velocity was greater for the latter events, the spatial variation in flow velocity and dissipation rate of turbulent kinetic energy (TKE) around seigyu was greater for the former events. Such flow patterns affected both bed formation and ecological habitats; ponds were formed at locations beside and behind seigyu, where the flow converged and bed scouring occurred during floods. The frequency and size of ponds around seigyu increased in the early half of the season, and they decreased in the other half when floods were greater in magnitude, which suggests that the bed scouring effect of seigyu was greater in flood stages with partial than fully submergence. Although the bar ponds lack shade to temper the effects of incident light and atmospheric conditions, the ponds displayed smaller daily oscillations in temperature than did the main river, probably due to hyporheic water supply to the ponds. Because many aquatic species cannot tolerate extremely high temperatures in summer, the generation and maintenance of deeper ponds by scouring with sufficient water exchange with the hyporheic zone can be key to enhancing colonization by various aquatic species.

Isotopic evidence for nitrate sources and controls on denitrification in groundwater beneath an irrigated agricultural district

The application of N fertilisers to enhance crop yield is common throughout the world. Many crops have historically been, or are still, fertilised with N in excess of the crop requirements. A portion of the excess N is transported into underlying aquifers in the form of NO3−, which is potentially discharged to surface waters. Denitrification can reduce the severity of NO3− export from groundwater. We sought to understand the occurrence and hydrogeochemical controls on denitrification in NO3−-rich aquifers beneath the Emerald Irrigation Area (EIA), Queensland, Australia, a region of extensive cotton and cereal production. Multiple stable isotope (in H2O, NO3−, DIC, DOC and SO42−) and radioactive isotope (3H and 36Cl) tracers were used to develop a conceptual N process model. Fertiliser-derived N is likely incorporated and retained in the soil organic N pool prior to its mineralisation, nitrification, and migration into aquifers. This process, alongside the near absence of other anthropogenic N sources, results in a homogenised groundwater NO3− isotopic signature that allows for denitrification trends to be distinguished. Regional-scale denitrification manifests as groundwater becomes increasingly anaerobic during flow from an upgradient basalt aquifer to a downgradient alluvial aquifer. Dilution and denitrification occurs in localised electron donor-rich suboxic hyporheic zones beneath leaking irrigation channels. Using approximated isotope enrichment factors, estimates of regional-scale NO3− removal ranges from 22 to 93% (average: 63%), and from 57 to 91% (average: 79%) beneath leaking irrigation channels. In the predominantly oxic upgradient basalt aquifer, raised groundwater tables create pathways for NO3− to be transported to adjacent surface waters. In the alluvial aquifer, the transfer of NO3− is limited both physically (through groundwater-surface water disconnection) and chemically (through denitrification). These observations underscore the need to understand regional- and local-scale hydrogeological processes when assessing the impacts of groundwater NO3− on adjacent and end of system ecosystems.

Effects of In-Channel Structure on Chinook Salmon Spawning Habitat and Embryo Production

Adult salmonids are frequently observed building redds adjacent to in-channel structure, including boulders and large woody debris. These areas are thought to be preferentially selected for a variety of reasons, including energy and/or predation refugia for spawners, and increased hyporheic exchange for incubating embryos. This research sought to quantify in-channel structure effects on local hydraulics and hyporheic flow and provide a mechanistic link between these changes and the survival, development, and growth of Chinook salmon Oncorhynchus tshawytscha embryos. Data were collected in an eight-kilometer reach, on the regulated lower Mokelumne River, in the California Central Valley. Nine paired sites, consisting of an area containing in-channel structure paired with an adjacent area lacking in-channel structure, were evaluated. Results indicated that in-channel structure disrupts surface water velocity patterns, creating pressure differences that significantly increase vertical hydraulic gradients within the subsurface. Overall, in-channel structure did not significantly increase survival, development, and growth of Chinook salmon embryos. However, at several low gradient downstream sites containing in-channel structure, embryo survival, development, and growth were significantly higher relative to paired sites lacking such features. Preliminary data indicate that adding or maintaining in-channel structure, including woody material, in suboptimal spawning reaches improves the incubation environment for salmonid embryos in regulated reaches of a lowland stream. More research examining temporal variation and a full range of incubation depths is needed to further assess these findings.

Zasolenie strefy hyporeicznej rzeki Moszczenicy w rejonie wysadu solnego Rogóźno

Zasolenie strefy hyporeicznej rzeki Moszczenicy w rejonie wysadu solnego Rogóźno

Hydrobiogechemical interactions in the hyporheic zone of a sulfate-impacted, freshwater stream and riparian wetland ecosystem.

Coupled abiotic and biotic processes in the hyporheic zone, where surface water and groundwater mix, play a critical role in the biogeochemical cycling of carbon, nutrients, and trace elements in streams and wetlands. Dynamic hydrologic conditions and anthropogenic pollution can impact redox gradients and biogeochemical response, although few studies examine the resulting hydrobiogeochemical interactions generated within the hyporheic zone. This study examines the effect of hyporheic flux dynamics and anthropogenic sulfate loading on the biogeochemistry of a riparian wetland and stream system. The hydrologic gradient as well as sediment, surface water, and porewater geochemistry chemistry was characterized at multiple points throughout the 2017 spring-summer-fall season at a sulfate-impacted stream flanked by wetlands in northern Minnesota. Results show that organic-rich sediments largely buffer the geochemical responses to brief or low magnitude changes in hydrologic gradient, but sustained or higher magnitude fluxes may variably alter the redox regime and, ultimately, the environmental geochemistry. This has implications for a changing climate that is expected to dramatically alter the hydrological cycle. Further, increased sulfate loading and dissolved or adsorbed ferric iron complexes in the hyporheic zone may induce a cryptic sulfur cycle linked to iron and carbon cycling, as indicated by the abundance of intermediate valence sulfur compounds (e.g., polysulfide, elemental sulfur, thiosulfate) throughout the anoxic wetland and stream-channel sediment column. The observed deviation from a classical redox tower coupled with potential changes in hydraulic gradient in these organic-rich wetland and stream hyporheic zones has implications for nutrient, trace element, and greenhouse gas fluxes into surface water and groundwater, ultimately influencing water quality and global climate.

Hyporheic Exchange Due to Cobbles on Sandy Beds

AbstractHyporheic exchange induced by riverbed topography and roughness provides important ecosystem functions. We investigated how cobble clusters embedded in a finer‐grained sand bed affect the near‐bed channel flow and the exchange of surface and subsurface water. We tested how the spacing and embeddedness of cobbles altered hyporheic exchange through a three‐dimensional fully coupled surface‐subsurface model. The 3D modeling framework captured the full physics of the exchange process locating the lateral position of upwelling zones on the side of the cobble. As the cobbles protrude more into the channel, eddies appear downstream, with smaller vortices where cobbles are closer together. In our simulations, hyporheic exchange increased with the spacing and protrusion ratio. The travel time in the hyporheic zone also increased with the protrusion ratio because the hyporheic flow paths are longer when the cobble is more protruded. With different spacings at the highest protrusion ratio, the average slope of early time breakthrough curve was steeper for large spacing settings: the travel time increases as the spacing decreases.

Seasonal Differences in the Attenuation of Polar Trace Organics in the Hyporheic Zone of an Urban Stream

AbstractAttenuation of trace organic compounds (TrOCs) in a river occurs to a large extent in its hyporheic zone. A major part of the attenuation of polar TrOCs is of microbial origin. As microbial activity depends on temperature and redox conditions, seasonal differences in TrOC attenuation are likely. We investigated TrOC attenuation at a river influenced by treated wastewater during two sampling campaigns, one in summer and one in winter. In addition to redox conditions and temperature, we also determined residence times of porewater in sediment using three methods: (a) non‐parametric deconvolution of electrical conductivity time series, (b) the model VFLUX 2.0 based on temperature time series (only summer), and (c) applying Darcy's law to differences in hydraulic heads (only summer). Contrary to our expectations, we found higher attenuation for 12 out of 18 TrOCs in winter, while three TrOCs were better attenuated in summer. Sediment conditions varied between seasons as more of the top sandy layer with a higher hydraulic permeability accumulated on the river bed in summer. As a result, residence times in the sediment were shorter in summer. In winter, longer residence times, lower temperatures, and a steeper oxygen gradient in sediment coincided with higher TrOC attenuation. Further research is needed to understand our unexpected findings and underlying mechanisms.

Removal of pollution accumulated in the process of wastewater discharge from the bottom layer of river systems

Three-dimensional numerical simulation has been performed to investigate the transport of impurities in a two-layer system “liquid-porous medium” for different hydrodynamic regimes of river systems. The dynamics of waste water containing heavy impurities in areas of industrial activity is modeled. The parameters of the numerical experiment correspond to the data of full-scale flow experiments in natural water bodies, part of the reservoir where the Solikamsk-Berezniki industrial hub is located. Due to the influx of brines, pollution accumulates in the hyporheic zone of the river, the fluid flow carries this pollution away through the porous medium, thus the problem is solved of washing out the brine distributed in the porous medium. The computational domain consists of a fluid mass (river) propagating over a porous medium (river bottom, hyporheic zone). It is shown that the concentration of polluting ingredients in the soil of the river bottom near the discharge points increases in the presence of "heavy" impurities. Such problems arise acutely when excess highly mineralized brines are discharged into water bodies located in the industrial areas with “heavy” wastewater, for which, due to significant intra-annual fluctuations in the water level and ice cover, the option with near-bottom wastewater disposal is considered preferable.

Structure and Function Response of Microbial Communities Towards Antibiotic Contamination in Hyporheic Zone Sediments

Structure and Function Response of Microbial Communities Towards Antibiotic Contamination in Hyporheic Zone Sediments

Degradation of Mcpa, Metolachlor and Propiconazole in the Hyporheic Zone of an Agriculturally Impacted River

Degradation of Mcpa, Metolachlor and Propiconazole in the Hyporheic Zone of an Agriculturally Impacted River

Drivers of Functional Diversity in the Hyporheic Zone of a Large River

Drivers of Functional Diversity in the Hyporheic Zone of a Large River

Environmental context mediates the interaction between fine sediment and hyporheic macroinvertebrate

AbstractExcessive deposition of fine‐sediment filling interstitial spaces in riverbed substrates is widely acknowledged to negatively impact the aquatic ecosystems. Hyporheic invertebrates, sensitive to ecological status, respond to the clogged substratum through more intense bioturbation. To better understand the interaction between hyporheic macroinvertebrates and fine sediment, a continuum of a naturally increasing fine‐sediment content in the Weihe River was investigated through field sampling in four seasons. Nonmetric multidimensional scaling (NMDS) and cluster analysis were employed to assess the distribution and the homogeneity in the composition variation of macroinvertebrate assemblages. Canonical correspondence analysis (CCA) was performed to reveal the context‐dependent relevance of bioturbation activity with bioturbation potential value. Our results suggest that the impact of fine‐sediment clogging on community composition and bioturbation activity was seasonally dependent due to the variability of environmental factors, such as flood events and water quality. The highest abundance and bioturbation activity occurred at moderately clogged sampling sites, while reductions in richness and bioturbation activity were observed at heavily clogged sampling sites, and the samples collected at some of the heavily clogged sites downstream did not even contain macroinvertebrates. The results indicated that the interaction between organisms and fine sediments is limited due to the complexity of the benthic and hyporheic zone. The study highlights the feasibility of hyporheic macroinvertebrates as indicators for the effects of fine‐sediment loading pressures on lotic ecosystems and the necessity of reach‐specific conservation and management strategies, especially in the downstream.

Dynamics of microbiotic patterns reveal surface water groundwater interactions in intermittent and perennial streams

Exchange between groundwater (GW), hyporheic zone waters (HZ) and surface waters (SW) is critical for water quality, quantity, and the ecological health and functioning of all three ecosystems. Hydrological exchange is particularly important in intermittent creeks, such as in the Murray Darling Basin, Australia, where stream reaches shift from losing to gaining depending on the volume of surface flows. In this study we used hydrochemistry to identify SW-GW exchange and combined this with eDNA data to analyse the response of eukaryote and prokaryote communities to differing flow conditions within intermittent and perennial stream reaches. Our study suggested that SW and GW microbial communities were only around 30% similar. Differences in microbiota between SW, HZ and GW habitats were driven by changes in relative abundances of surface water dominant organisms (such as those capable of photosynthesis) as well as anaerobic taxa typical of GW environments (e.g., methanogens), with GW and HZ microbial communities becoming increasingly different to those in SW as flow ceased in intermittent creeks. Fine-scale hydrologic changes were identified through microbial communities in the perennial Maules Creek, indicating the importance of GW-SW exchange to biotic communities. This study highlights the importance of flow in shaping microbial communities and biogeochemical cycling within intermittent creeks and their connected alluvial aquifers. Our results suggest that microbiota may prove a useful indicator of SW-GW exchange, and in some circumstances, may be more sensitive in demonstrating fine-scale changes in SW-GW interactions than water chemistry. This knowledge furthers our understanding of GW-SW exchange and its impacts on ecological health.

Role of River–Lake System Sediments and Microbial Activity in the Hyporheic Zone

The effect of river–lake systems on the surface water self-purification process is a significant and not fully recognised scientific issue. The conditions prevailing in the hyporheic zone of these ecosystems are of great importance in the process of component exchange between water and sediments. The aim of this study was to investigate the influence of the type of sediments located at the bottom of the riverbed being part of a river–lake system on microbial activity in the hyporheic zone. An ex situ experiment was used to study the microbiological activity and the transformation of components in the collected river sediments. It was found that the specific properties of sediments varied depending on their location in the riverbed between the lakes comprising the system and that the prevailing meteorological conditions can also have an effect on microbial activity in the hyporheic zone, e.g., aerobic conditions. These conditions determined the intensity of component conversion in the sediments due to microbial metabolism. A closer understanding of the processes occurring in the hyporheic zone may allow the processes of water self-purification within river–lake systems to be supported in the future, which will contribute to the improvement of surface water quality.

A systematic model-based evaluation of the influence of hydraulic conductivity, heterogeneity and domain depth on hyporheic nutrient transformation

The hyporheic zone provides opportunities for the turnover and removal of nutrients such as ammonium and nitrate from stream water. Hyporheic flow patterns and dynamics govern the transport of nutrients and other solutes through the hyporheic zone and in turn control their transformation and removal. To better understand the links between flow and turnover, we systematically studied nutrient turnover rates with simulated sets of idealized streambed domains which varied in domain depth, hydraulic conductivity, its heterogeneity and anisotropy. We show that both aerobic (e.g. nitrification) and anaerobic (e.g. denitrification) turnover rates generally decrease monotonically with decreasing sediment conductivity at sufficiently large hyporheic zone depth. However, a reduction of the domain depth leads to an optimum of anaerobic turnover rates at intermediate conductivity due to the rate or transit time limitations of the reaction at the far ends of the conductivity range. Furthermore, conductivity heterogeneity generally lowers the turnover rates of both aerobic and anaerobic reactions. It is demonstrated how to reproduce this reduced nutrient turnover in heterogeneous sediments via a simplified anisotropy, which facilitates the prediction of turnover rates in heterogeneous sediments in certain cases. Our results foster the understanding of the interplay between hyporheic exchange and nutrient transformations and they aim to facilitate the conceptualization of large scale numerical simulations by demonstrating the effect of study design decisions like domain extent, heterogeneity and anisotropy.

Identification of Characteristic Spatial Scales to Improve the Performance of Analytical Spectral Solutions to the Groundwater Flow Equation

AbstractAnalytical solutions for the three‐dimensional groundwater flow equation have been widely used to gain insight about subsurface flow structure and as an alternative to computationally expensive numerical models. Of particular interest are solutions that decompose prescribed hydraulic head boundaries (e.g., Dirichlet boundary condition) into a collection of harmonic functions. Previous studies estimate the frequencies and amplitudes of these harmonics with a least‐square approach where the amplitudes are fitted given a pre‐assigned set of frequencies. In these studies, an ad hoc and structured discretization of the frequency domain is typically used, excluding dominant frequencies while assigning importance to spurious frequencies, with significant consequences for estimating the fluxes and residence times. This study demonstrates the advantages of using a pre‐assigned frequency spectrum that targets the dominant frequencies based on rigorous statistical analysis with predefined significance levels. The new approach is tested for three hydrologic conceptualizations: (a) a synthetic periodic basin, (b) synthetic bedforms, and (c) a natural mountainous watershed. The performance of the frequency spectrum selection is compared with exact analytical or approximate numerical solutions. We found that the new approach better describes the fluxes and residence times for Dirichlet boundaries with well‐defined characteristics spatial scales (e.g., periodic basins and bedforms). For more complex scenarios, such as natural mountainous watersheds, both pre‐assigned frequency spectrums present similar performance. The spectral solutions presented here can play a central role in developing reduced‐complexity models for assessing regional water and solute fluxes within mountain watersheds and hyporheic zones.

Improving Predictions of Stream CO2 Concentrations and Fluxes Using a Stream Network Model: A Case Study in the East River Watershed, CO, USA

AbstractInland waters are an important component of the global carbon budget. However, our ability to predict carbon fluxes from stream systems remains uncertain, as pCO2 varies within streams at scales of 1–100 m. This makes direct monitoring of large‐scale CO2 fluxes impractical. We incorporate CO2 input and output fluxes into a stream network advection‐reaction model, representing the first process‐based representation of stream CO2 dynamics at watershed scales. This model includes groundwater (GW) CO2 inputs, water column (WC), benthic hyporheic zone (BHZ) respiration, downstream advection, and atmospheric exchange. We evaluate this model against existing statistical methods including upscaling and multiple linear regressions through comparisons to high‐resolution stream pCO2 data collected across the East River Watershed in the Colorado Rocky Mountains (USA). The stream network model accurately captures GW, evasion, and respiration‐driven pCO2 variability and significantly outperforms multiple linear regressions for predicting pCO2. Further, the model provides estimates of CO2 contributions from internal versus external sources suggesting that streams transition from GW‐ to BHZ‐dominated sources between 3rd and 4th Strahler orders, with GW, BHZ, and WC accounting for 49.3%, 50.6%, and 0.1% of CO2 fluxes from the watershed, respectively. Lastly, stream network model atmospheric CO2 fluxes are 4‐12x times smaller than upscaling technique predictions, largely due to relationships between stream pCO2 and gas exchange velocities. Taken together, this stream network model improves our ability to predict stream CO2 dynamics and efflux. Furthermore, future applications to regional and global scales may result in a significant downward revision of global flux estimates.