Nico Trauth

Hyporheic transport and biogeochemical reactions in pool‐riffle systems under varying ambient groundwater flow conditions

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 O2 and NO3 concentrations in the upwelling groundwater reveal that DN efficiency of the HZ is highest under low upwelling magnitudes accompanied with low concentrations of O2 and NO3. 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.

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

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 × 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.

Single discharge events increase reactive efficiency of the hyporheic zone

In 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.

Robust Optode-Based Method for Measuring in Situ Oxygen Profiles in Gravelly Streambeds

One of the key environmental conditions controlling biogeochemical reactions in aquatic sediments like streambeds is the distribution of dissolved oxygen. We present a novel approach for the in situ measurement of vertical oxygen profiles using a planar luminescence-based optical sensor. The instrument consists of a transparent acrylic tube with the oxygen-sensitive layer mounted on the outside. The luminescence is excited and detected by a moveable piston inside the acrylic tube. Since no moving parts are in contact with the streambed, the disturbance of the subsurface flow field is minimized. The precision of the distributed oxygen sensor (DOS) was assessed by a comparison with spot optodes. Although the precision of the DOS, expressed as standard deviation of calculated oxygen air saturation, is lower (0.2-6.2%) compared to spot optodes (<0.1-0.6%), variations of the oxygen content along the profile can be resolved. The uncertainty of the calculated oxygen is assessed with a Monte Carlo uncertainty assessment. The obtained vertical oxygen profiles of 40 cm in length reveal variations of the oxygen content reaching from 90% to 0% air saturation and are characterized by patches of low oxygen rather than a continuous decrease with depth.

Estimating time‐variable aerobic respiration in the streambed by combining electrical conductivity and dissolved oxygen time series

Aerobic respiration is an important component of in-stream metabolism. The larger part occurs in the streambed, where it is difficult to directly determine actual respiration rates. Existing methods for determining respiration are based on indirect estimates from whole-stream metabolism or provide time invariant results estimated from oxygen consumption measurements in enclosed chambers that do not account for the influence of hydrological changes. In this study we demonstrate a simple method for determining time-variable hyporheic respiration. We use a windowed cross-correlation approach for deriving time-variable travel times from the naturally changing electrical conductivity signal that is transferred into the sediment. By combining the results with continuous in situ dissolved oxygen measurements, variable oxygen consumption rate coefficients in the streambed are obtained. An empirical temperature relationship is derived and used for standardizing the respiration rate coefficients to isothermal conditions. For demonstrating the method, we compare two independent measurement spots in the streambed, which were located upstream and downstream of an in-stream gravel bar and thus exposed strongly diverse travel times. The derived respiration rate results are in accordance with findings of other stream studies. By comparing the travel time and respiration rate coefficient (i.e., Damkohler number) we estimate the contribution of each to the oxygen consumption in the streambed.

River water infiltration enhances denitrification efficiency in riparian groundwater

Nitrate contamination in ground- and surface water is a persistent problem in countries with intense agriculture. The transition zone between rivers and their riparian aquifers, where river water and groundwater interact, may play an important role in mediating nitrate exports, as it can facilitate intensive denitrification, which permanently removes nitrate from the aquatic system. However, the in-situ factors controlling riparian denitrification are not fully understood, as they are often strongly linked and their effects superimpose each other. In this study, we present the evaluation of hydrochemical and isotopic data from a 2-year sampling period of river water and groundwater in the riparian zone along a 3rd order river in Central Germany. Based on bi- and multivariate statistics (Spearmans rank correlation and partial least squares regression) we can show, that highest rates for oxygen consumption and denitrification in the riparian aquifer occur where the fraction of infiltrated river water and at the same time groundwater temperature, are high. River discharge and depth to groundwater are additional explanatory variables for those reaction rates, but of minor importance. Our data and analyses suggest that at locations in the riparian aquifer, which show significant river water infiltration, heterotrophic microbial reactions in the riparian zone may be fueled by bioavailable organic carbon derived from the river water. We conclude that interactions between rivers and riparian groundwater are likely to be a key control of nitrate removal and should be considered as a measure to mitigate high nitrate exports from agricultural catchments.

Evaluating the reliability of time series analysis to estimate variable riparian travel times by numerical groundwater modelling

The transition zones between rivers and adjacent riparian aquifers are locations of high biogeochemical activities that contribute to a removal of potentially hazardous substances in the aquatic system. The potential of the removal processes depends highly on subsurface water travel times, which can be determined by using the propagation of electrical conductivity (EC) signal from the river into the riparian aquifer. Although this method has been applied and verified in many studies, we observe possible limitations for the usage of EC fluctuation analysis. Our findings are based on EC time series analyses during storm events and artificial hydropeaks induced by watermill operations. Travel times derived by cross‐correlation analysis were compared with travel times calculated based on backward particle tracking of a calibrated transient numerical groundwater flow model. The cross‐correlation method produced only reasonable travel times for the artificial hydropeaks. In contrast, cross‐correlation analysis of the EC data during natural storm events resulted in implausibly negative or unrealistically low travel times for the bulk of the data sets. We conclude that the reason for this behaviour is, first, the low EC contrast between river and groundwater in connection with a strong damping of the infiltrating river EC signal into the subsurface during storm events. Second, the existence of old and less‐mineralized riparian water between the river and the monitoring well resulted in bank‐storage‐driven EC breakthrough curves with earlier arrival times and the subsequent estimation of implausible riparian travel times.

Gravel bars are sites of increased CO2 outgassing in stream corridors

Streams are significant sources of CO2 to the atmosphere. Estimates of CO2 evasion fluxes (fCO2) from streams typically relate to the free flowing water but exclude geomorphological structures within the stream corridor. We found that gravel bars (GBs) are important sources of CO2 to the atmosphere, with on average more than twice as high fCO2 as those from the streamwater, affecting fCO2 at the level of entire headwater networks. Vertical temperature gradients resulting from the interplay between advective heat transfer and mixing with groundwater within GBs explained the observed variation in fCO2 from the GBs reasonably well. We propose that increased temperatures and their gradients within GBs exposed to solar radiation stimulate heterotrophic metabolism therein and facilitate the venting of CO2 from external sources (e.g. downwelling streamwater, groundwater) within GBs. Our study shows that GB fCO2 increased fCO2 from stream corridors by [median, (95% confidence interval)] 16.69%, (15.85–18.49%); 30.44%, (30.40–34.68%) and 2.92%, (2.90–3.0%), for 3rd, 4th and 5th order streams, respectively. These findings shed new light on regional estimates of fCO2 from streams, and are relevant given that streamwater thermal regimes change owing to global warming and human alteration of stream corridors.

High resolution monitoring above and below the groundwater table uncovers small – scale hydrochemical gradients

Hydrochemical solute concentrations in the shallow subsurface can be spatially highly variable within small scales, particularly at interfaces. However, most monitoring systems fail to capture these small scale variations. Within this study, we developed a high resolution multilevel well (HR-MLW) with which we monitored water across the interface of the unsaturated and saturated zone with a vertical resolution of 0.05-0.5 m. We installed three of these 4 m deep HR-MLWs in the riparian zone of a third-order river and analyzed for hydrochemical parameters and stable water isotopes. The results showed three distinct vertical zones (unsaturated zone, upper saturated zone, lower saturated zone) within the alluvial aquifer. A 2 m thick layer influenced by river water (upper saturated zone) was not captured by existing monitoring wells with higher screen length. Hydrochemical data (isotopes, total ions) were consistent in all HR-MLWs and showed similar variation over time emphasizing the reliability of the installed monitoring system. Further, the depths zones were also reflected in the NO3-N concentrations; with high spatial variabilities between the three wells. The zonation was constant over time, with seasonal variability in the upper saturated zone due to the influence of river water. This study highlights the use of high resolution monitoring for identifying the spatial and temporal variability of hydrochemical parameters present in many aquifer systems. Possible applications range from riparian zones, agricultural field sites to contaminated site studies, wherever an improved understanding of biogeochemical turnover processes is necessary.

River recharge versus O 2 supply from the unsaturated zone in shallow riparian groundwater: A case study from the Selke River (Germany)

Besides gas-water-exchange in surface waters, respiratory consumption of dissolved oxygen (DO) in adjacent riparian groundwater may trigger the addition of so far hardly explored sources from the unsaturated zone. These processes also systematically influence stable isotope ratios of DO and were investigated together with Cl- as a conservative tracer for water mixing in a near-river riparian groundwater system. The study focused on a losing stream section of the Selke River at the foot of the Harz Mountains (Germany). The study area exposed steep DO gradients between the stream water and riparian groundwater between April 2016 and May 2017. Our results indicated dominant influences of microbial community respiration with observed DO concentration gradients. These observations can be explained by DO from the river that is subject to fractionation by microbial respiration with a typical fractionation factor (αr) of 0.982. However, with such respiration dominance, we expected a simultaneous enrichment of δ18ODO towards values that are more positive than the well-known atmospheric O2 signal of +23.9‰ versus the Vienna Standard Mean Ocean Water standard (VSMOW). Surprisingly, our measurements revealed much lower δ18ODO values between +22‰ and +18‰ in the near-river groundwater. Mass balance calculations revealed that the DO pool in the shallow and unconfined aquifer receives contributions of up to about 80% by diffusion of oxygen from the vadose zone with a distinctly lower isotope value than the one of the atmosphere. This finding about additional oxygen sources from the unsaturated zone has numerous ramifications for oxygen related processes in near-river environments including the oxidation of pollutants, nutrients and ecosystem health.