The fate of polar trace organic compounds in the hyporheic zone

How redox gradient potentially influences nitrate reduction coupled with sulfur cycling: A new insight into nitrogen cycling in the hyporheic zone of effluent-dominated rivers

The coupled N and S cycling in variable redox gradients in the hyporheic zone (HZ) of the rivers receiving effluents from wastewater treatment plants is unclear. Using two representative effluent-dominated rivers as model systems, a metagenome approach was employed to explore the spatiotemporal redox zonation of the HZ and the N/S cycling processes within the system. The results manifested that nitrate reduction represented the fundamental nitrogen pathway in the HZ. Interestingly, DNRA coupled with sulfur reduction, and denitrification coupled with sulfur oxidation were respectively abundant in the oxic and anoxic zone. Lower nitrate concentration (0-2.72 mg-N/L) and more abundant genes involved in denitrification (napB, NarGHI) and sulfur oxidation (sseA, glpE) were detected in the anoxic zone. Contrarily, the nitrate concentration (0.07-4.87 mg-N/L) and the abundance of genes involved in sulfur reduction (ttrB, sudA) and DNRA (nirBD) were observed more abundant in the oxic zone. Therefore, the results verified the oxygen-limited condition did not suppress but rather facilitated the denitrification process in the presence of active S cycling. The high relative abundances of nosZ gene encoding sequence (3-5 % relative to all nitrogen-cycling processes) in both the effluent-discharging area and downstream area highly confirmed that HZ was capable of alleviating the N2O emission in the region. The functional keystone taxa were revealed through co-occurrence network analysis. The structural equation model shows that the genes of N/S cycling were positively impacted by functional keystone taxa, especially the N cycling genes. Functional keystone taxa were proven driven by the redox gradient, demonstrating their positive roles in mediating N/S cycling processes. The promoting effect on nitrate reduction coupled with sulfur cycling was clarified when redox conditions oscillated, providing a new perspective on mitigating nitrogen pollution and greenhouse gas emissions in effluent-receiving rivers.

Deciphering solute transport, microbiota assembly patterns and metabolic functions in the hyporheic zone of an effluent-dominated river

We lack a clear understanding of how anthropogenic pressures, exemplified by effluent discharge from wastewater treatment plants, destabilize microbial communities in the hyporheic zone (HZ) of receiving rivers. In this study, the spatiotemporal characteristics of hydrological parameters, and the physicochemical properties of surface and subsurface water in a representative effluent-dominated river were monitored. Sequencing of 16S rRNA amplicons and metagenomes revealed the microbial community structure in the HZ of both effluent discharge area and downstream region. The keystone taxa (taxa vital in determining the composition of each microbial cluster) and the keystone functions they controlled were subsequently identified. Effluent discharge amplified the depth of the oxic/suboxic zone and the hyporheic exchange fluxes in the effluent discharge area, which was 50-120% and 40-300% higher than in the downstream region, respectively. Microbial community structure pattern analysis demonstrated an enhancement in the rate of dispersal, an increase in microbial diversity, and an improved community network complexity in the effluent discharge area. By contrast, the number of keystone taxa in the effluent discharge area was 50-70% lower than that of the downstream region, resulting in reduced community network stability and functionality. The keystone taxa controlling metabolic functions in the networks categorized to effluent discharge area were comprised of more genera related to nitrogen and sulfur cycling, e.g., Dechloromonas, Desulfobacter, Flavobacterium, Nitrosomonas, etc., highlighting a research need in monitoring species associated with nutrient element cycling in the HZ of receiving waterbodies. The results showed that the keystone taxa could contribute positively to network stability, which was negatively correlated to hyporheic exchange fluxes and redox gradients. This study provides valuable insights that will improve our understanding of how river ecosystems respond to changes in anthropogenic pressures.

Redox gradients drive microbial community assembly patterns and molecular ecological networks in the hyporheic zone of effluent-dominated rivers

The impacts of effluent discharge on receiving waterbodies have been a research hotspot. Nonetheless, limited information is available on the microbial community assembly patterns in the hyporheic zones (HZ) responding to the changes in the microenvironments, e.g., solute transport and redox gradient variations. Using two representative effluent-dominated rivers as model systems, the spatio-temporal bacterial community dynamics and assembly patterns in oxic and suboxic zones in the shallow riverbed sediments were disentangled via null model- and neutral model-based approaches. Bacterial dynamics in community composition were observed driven by environmental filtering, i.e., impacts of environmental variables, more than geographic distances, i.e., the depths of sediments. The communities in samples collected in summer were largely shaped by stochasticity, in which homogeneous selection occupied a higher proportion in oxic (∼39%) than in suboxic zone (∼23%). Deterministic processes contributed to a more complex community structure for samples from oxic zones, whereas weakened the interspecies interactions in suboxic zones. The richness and abundances of non-neutral community were confirmed governing the deterministic assembly in oxic zones. Key species ascribed to 'connectors' and 'network hubs' dominated the community assembly variations in samples collected in winter, and in oxic zones, respectively. Significant positive relationships between β-nearest taxon index and dissolved organic nitrogen (DON) and nitrate highlighted their vital roles in community assembly via deterministic selective pressures in oxic zones. The significance thresholds of nitrogen species for community transition in winter (ΔDON: 2.81 mg-N/L, ΔNO3-: 1.09 mg-N/L) were lower than in summer, probably implying that stricter effluent quality standards should be established in colder seasons. Combined, our work poses first insights on the roles of redox zonation in driving microbial community assembly in HZ, which is of significance in guiding ecological remediation processes in effluent-dominated rivers.

Biodegradation of trace sulfonamide antibiotics accelerated by substrates across oxic to anoxic conditions during column infiltration experiments

Frequent occurrence of trace organic contaminants in aquatic environments, such as sulfonamide antibiotics in rivers receiving reclaimed water, is concerning. Natural attenuation by soil and sediment is increasingly relied upon. In the case of riverbank filtration for water purification, the reliability of antibiotic attenuation has been called into question due to incomplete understanding of their degradation processes. This study investigated influence of substrates and redox evolution along infiltration path on biotransformation of sulfonamides. Eight sand columns (length: 28 cm) with a riverbed sediment layer at 3-8 cm were fed by groundwater-sourced tap water spiked with 1 μg/L of sulfadiazine (SDZ), sulfamethazine (SMZ), and sulfamethoxazole (SMX) each, with or without amendments of dissolved organic carbon (5 mg-C/L of 1:1 yeast and humics) or ammonium (5 mg-N/L). Two flow rates were tested over 120 days (0.5 mL/min and 0.1 mL/min). Iron-reducing conditions persisted in all columns for 27 days during the initial high flow period due to respiration of sediment organics, evolving to less reducing conditions until the subsequent low flow period to resume more reducing conditions. With surplus substrates, the spatial and temporal patterns of redox conditions differentiated among columns. The removal of SDZ and SMZ in effluents was usually low (15 ± 11%) even with carbon addition (14 ± 9%), increasing to 33 ± 23% with ammonium addition. By contrast, SMX removal was higher and more consistent among columns (46 ± 21%), with the maximum of 64 ± 9% under iron-reducing conditions. When sulfonamide removal was compared between columns for the same redox zones during infiltration, their enhancements were always associated with the availability of dissolved or particulate substrates, suggesting co-metabolism. Manipulation of the exposure time to optimal redox conditions with substrate amendments, rather than to simply prolong the overall residence time, is recommended for nature-based solutions to tackle target antibiotics.

Prediction models and major controlling factors of antibiotics bioavailability in hyporheic zone

Recently, antibiotics have been frequently detected in the hyporheic zone (HZ) as a novel contaminant. Bioavailability assessment has gradually attracted more attention in order to provide a more realistic assessment of human health risks. In this study, two typical antibiotics, oxytetracycline (OTC) and sulfamethoxazole (SMZ), were used as target pollutants in the HZ of the Zaohe-Weihe River, and the polar organics integrated sampler was used to analyze the variation of antibiotics bioavailability. According to the characteristics of the HZ, the total concentration of pollutants, pH, and dissolved oxygen (DO) were selected as major predictive factors to analyze their correlation with the antibiotics bioavailability. Then the predictive antibiotic bioavailability models were constructed by stepwise multiple linear regression method. The results showed that there was a highly significant negative correlation between OTC bioavailability and DO (P < 0.001), while SMZ bioavailability showed a highly significant negative correlation with total concentration of pollutants (P < 0.001) and a significant negative correlation with DO (P < 0.01). The results of correlation analysis were further verified by Principal Component Analysis. Based on the experimental data, we constructed eight prediction models for the bioavailability of two antibiotics and verified them. The data points of the six prediction models were distributed in the 95% prediction band, indicating that the models were more reliable and accurate. The prediction models in this study provide reference for the accurate ecological risk assessment of the bioavailability of pollutants in the HZ, and also provide a new idea for predicting the bioavailability of pollutants in practical applications.

Experimental and numerical elucidation of the fate and transport of antibiotics in aquatic environment: A review

This review highlights various experimental and mathematical modeling strategies to investigate the fate and transport of antibiotics that elucidate antimicrobial selective pressure in aquatic environments. Globally, the residual antibiotic concentrations in effluents from bulk drug manufacturing industries were 30- and 1500-fold greater than values reported in municipal and hospital effluents, respectively. The antibiotic concentration from different effluents enters the waterbodies that usually get diluted as they go downstream and undergo various abiotic and biotic reactive processes. In aquatic systems, photolysis is the predominant process for antibiotic reduction in the water matrix, while hydrolysis and sorption are frequently reported in the sediment compartment. The rate of antibiotic reduction varies widely with influencing factors such as the chemical properties of the antibiotics and hydrodynamic conditions of river streams. Among all, tetracycline was found to more unstable (log Kow =  - 0.62 to - 1.12) that can readily undergo photolysis and hydrolysis; whereas macrolides were more stable (log Kow = 3.06 to 4.02) that are prone to biodegradation. The processes like photolysis, hydrolysis, and biodegradation followed first-order reaction kinetics while the sorption followed a second-order kinetics for most antibiotic classes with reaction rates occurring in the decreasing order of Fluoroquinolones and Sulphonamides. The reports from various experiments on abiotic and biotic processes serve as input parameters for an integrated mathematical modeling to predict the fate of the antibiotics in the aquatic environment. Various mathematical models viz. Fugacity level IV, RSEMM, OTIS, GREAT-ER, SWAT, QWASI, and STREAM-EU are discussed for their potential capabilities. However, these models do not account for microscale interactions of the antibiotics and microbial community under real-field conditions. Also, the seasonal variations for contaminant concentrations that exert selective pressure for antimicrobial resistance has not been accounted. Addressing these aspects collectively is the key to exploring the emergence of antimicrobial resistance. Therefore, a comprehensive model involving antimicrobial resistance parameters like fitness cost, bacterial population dynamics, conjugation transfer efficiency, etc. is required to predict the fate of antibiotics.

Metformin Contamination in Global Waters: Biotic and Abiotic Transformation, Byproduct Generation and Toxicity, and Evaluation as a Pharmaceutical Indicator

Metformin is the first-line antidiabetic drug and one of the most prescribed medications worldwide. Because of its ubiquitous occurrence in global waters and demonstrated ecotoxicity, metformin, as with other pharmaceuticals, has become a concerning emerging contaminant. Metformin is subject to transformation, producing numerous problematic transformation byproducts (TPs). The occurrence, removal, and toxicity of metformin have been continually reviewed; yet, a comprehensive analysis of its transformation pathways, byproduct generation, and the associated change in adverse effects is lacking. In this review, we provide a critical overview of the transformation fate of metformin during water treatments and natural processes and compile the 32 organic TPs generated from biotic and abiotic pathways. These TPs occur in aquatic systems worldwide along with metformin. Enhanced toxicity of several TPs compared to metformin has been demonstrated through organism tests and necessitates the development of complete mineralization techniques for metformin and more attention on TP monitoring. We also assess the potential of metformin to indicate overall contamination of pharmaceuticals in aquatic environments, and compared to the previously acknowledged ones, metformin is found to be a more robust or comparable indicator of such overall pharmaceutical contamination. In addition, we provide insightful avenues for future research on metformin.

Fate of trace organic compounds in the hyporheic zone: Influence of microbial metabolism

The hyporheic zone (HZ) is considered a hydrodynamically-driven bioreactor with significant pollutant removal capacities and can therefore not only improve wholestream water quality but also preserve human and ecosystem health. Microbial metabolism is hypothesized to play a key role in pollutant transformation in hyporheic sediments of natural streams. However, previous work investigating the influence of microbial metabolism on pollutant transformation has been predominantly laboratory studies. The key challenge for field studies is the appropriate determination of net microbial metabolism, i.e. information on the actual exposure times to specific microbial processes in the investigated system. The present study uses reactive fluorescent tracers to determine microbial metabolism and ultimately its influence on pollutant transformation, e.g. for trace organic compounds, in hyporheic sediments under natural conditions. In particular, the reactive fluorescent tracers resazurin and its main transformation product resorufin were used to determine the microbial metabolism of facultative or obligate aerobes. The influence of the derived microbial metabolism on the transformation of 20 trace organic compounds, such as pharmaceuticals, including 3 parent-daughter pairs, was examined. The present findings validate laboratory results on the microbially-mediated transformation of the anticonvulsant gabapentin to its main transformation product gabapentin lactam under natural conditions. All other TrOCs investigated did not show a clear link between TrOC reactivity to the microbial metabolism informed by the resazurin-resorufin-system. Overall, the present study not only demonstrates the use of the fluorescent tracer-system resazurin and resorufin for determining microbial metabolism of facultative or obligate aerobes but also generally highlights the potential of reactive fluorescent tracers to disentangle specific reactive properties and ultimately their influence on the fate of pollutants in natural HZs.

In-stream sorption of azithromycin and levofloxacin in a river receiving sewage treatment plant effluent

Modelling natural attenuation is crucial to managing pharmaceuticals. However, little is known about the mechanism behind their in-stream sorption. To better understand the in-stream attenuation of the highly sorptive antibiotics azithromycin (AZM) and levofloxacin (LVF), we monitored them in a 2.1-km stretch of the Asano River under diverse flow conditions. This stretch receives effluent directly from a sewage treatment plant (STP), which was a dominant source of the pharmaceuticals. Average distribution coefficients between dissolved and particulate phases (K_d,SPM) in the outflow river water were 6.3×10⁵ L/kg for AZM and 7.5×10⁴ L/kg for LVF, while those in the STP effluent were 1-2 orders of magnitude lower. Mass balances in the river stretch calculated by considering only dissolved phase (MB_w) and both dissolved and particulate phases (MB_s) were 8%-52% and 58%-102%, respectively, for AZM, and 58%-71% and 60%-105% for LVF. MB_w<MB_s is attributed to an increase in suspended particulate matter (SPM)-mediated mass flows in the river stretch, i.e., in-stream sorption to SPM, which was caused mainly by their much higher river K_d,SPM values than those in the effluent. Their river K_d,SPM values increased on higher-flow days with decreasing effluent content in the river water, resulting in the increase of their in-stream SPM sorption. Their in-stream loss from the entire water column (i.e., 100-MB_s), which was attributable to their mass transfer from the overlying water to sediment through sorption, was decreased on higher-flow days by hydrological factors. A key finding is that AZM and LVF mostly entered the river stretch in the dissolved phase of STP effluent, whereas they existed substantially in the particulate phase in the outflow river water, especially on high-flow days.

The biogeochemical responses of hyporheic groundwater to the long-run managed aquifer recharge: Linking microbial communities to hydrochemistry and micropollutants

Interactions of surface water and groundwater (SW-GW) in hyporheic zones produce biogeochemical hotspots. However, response patterns of hyporheic groundwater to external influences remain unclear. In this study, three datasets (hydrochemistry, antibiotics, and microbiome) were collected over a hydrological year to explore the influence of a 12-year managed aquifer recharge (MAR) project. We observed that the long-term MAR practice elevated nutrient and antibiotic levels while reduced redox potential in hyporheic groundwater, and these impacts depended on decreasing SW-GW interaction intensity with aquifer depth. In contrast, the long-term MAR practice increased community dissimilarity of 30-m groundwater but had little impact on 50-m or 80-m groundwater. Moreover, hyporheic community assembly was dominated by dispersal limitation, and thereby co-varied hydrochemistry and antibiotics only attributed to small community variability. The long-term MAR practice decreased species-interaction intensity and changed the abundance of metabolic functions in hyporheic groundwater. Furthermore, predicted community functions involving carbon, nitrogen, sulfur, and manganese cycles for 30-m groundwater showed higher abundances than those for 50- and 80-m groundwater. Collectively, we showed that hyporheic groundwater was sensitive to the SW-GW interaction and human activities, with the interactions of hydrochemistry, contaminants, and microbiome linking to hyporheic groundwater quality and ecosystem functioning.

Effects of ecohydrological interfaces on migrations and transformations of pollutants: A critical review

With the rapid development of society, the soil and water environments in many countries are suffering from severe pollution. Pollutants in different phases will eventually gather into the soil and water environments, and a series of migrations and transformations will take place at ecohydrological interfaces with water flow. However, it is still not clear how ecohydrological interfaces affect the migration and the transformation of pollutants. Therefore, this paper summarizes the physical, ecological, and biogeochemical characteristics of ecohydrological interfaces on the basis of introducing the development history of ecohydrology and the concept of ecohydrological interfaces. The effects of ecohydrological interfaces on the migration and transformation of heavy metals, organic pollutants, and carbon‑nitrogen‑phosphorus (C-N-P) pollutants are emphasized. Lastly, the prospects of applying ecohydrological interfaces for the removal of pollutants from the soil and water environment are put forward, including strengthening the ability to monitor and simulate ecohydrological systems at micro and macro scales, enhancing interdisciplinary research, and identifying main influencing factors that can provide theoretical basis and technical support.

Determining hyporheic removal rates of trace organic compounds using non-parametric conservative transport with multiple sorption models

Assessing the transport and reactive processes of contaminants in freshwater streams is crucial in managing water resources sustainably. Particularly the hyporheic zone, the sediment-water interface where surface water and groundwater mix, may possess significant contaminant removal capacities due to its myriad physical, chemical, and microbiological processes. However, modelling approaches aiming at assessing the hyporheic zone's reactivity are either based on simple assumptions, such as, predefining the shape of the residence times distribution (RTD) function, or are computationally not feasible due to a too detailed system characterisation. In addition, parent-daughter reactions of contaminants are barely investigated. The present study introduces a numerical modelling framework for assessing hyporheic reactions of contaminant transformation reactions based on a non-parametric residence time approach combined with multiple sorption models and first-order removal reactions. The proposed framework uses natural electrical conductivity fluctuations to determine conservative transport properties and is demonstrated by interpreting time series of hyporheic point measurements of trace organic compounds, such as pharmaceuticals, and their transformation products using two commonly-used sorption models, namely the simple retardation and the first-order kinetic sorption model. The developed approach gives similar reaction rate coefficient estimates for all contaminants considered for both sorption models tested. The findings highlight that (i) the accurate shape of the RTD is most certainly important for reactive parameter determination and (ii) the daughter reaction rate coefficient may be underestimated if its parent transformation is ignored. The model provides reactive parameter estimates of contaminant transformation reactions with high parameter identifiability and informs which specific parent-daughter-pathway has occurred.

Biodegradation of metoprolol in oxic and anoxic hyporheic zone sediments: unexpected effects on microbial communities

Abstract

Metoprolol is widely used as a beta-blocker and considered an emerging contaminant of environmental concern due to pseudo persistence in wastewater effluents that poses a potential ecotoxicological threat to aquatic ecosystems. Microbial removal of metoprolol in the redox-delineated hyporheic zone (HZ) was investigated using streambed sediments supplemented with 15 or 150 μM metoprolol in a laboratory microcosm incubation under oxic and anoxic conditions. Metoprolol disappeared from the aqueous phase under oxic and anoxic conditions within 65 and 72 days, respectively. Metoprolol was refed twice after initial depletion resulting in accelerated disappearance under both conditions. Metoprolol disappearance was marginal in sterile control microcosms with autoclaved sediment. Metoprolol was transformed mainly to metoprolol acid in oxic microcosms, while metoprolol acid and α-hydroxymetoprolol were formed in anoxic microcosms. Transformation products were transient and disappeared within 30 days under both conditions. Effects of metoprolol on the HZ bacterial community were evaluated using DNA- and RNA-based time-resolved amplicon Illumina MiSeq sequencing targeting the 16S rRNA gene and 16S rRNA, respectively, and were prominent on 16S rRNA rather than 16S rRNA gene level suggesting moderate metoprolol-induced activity-level changes. A positive impact of metoprolol on Sphingomonadaceae and Enterobacteriaceae under oxic and anoxic conditions, respectively, was observed. Nitrifiers were impaired by metoprolol under oxic and anoxic conditions. Collectively, our findings revealed high metoprolol biodegradation potentials in the hyporheic zone under contrasting redox conditions associated with changes in the active microbial communities, thus contributing to the attenuation of micropollutants.

Key points

• High biotic oxic and anoxic metoprolol degradation potentials in the hyporheic zone.

• Key metoprolol-associated taxa included Sphingomonadaceae, Enterobacteraceae, and Promicromonosporaceae.

• Negative impact of metoprolol on nitrifiers.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11466-w.

Macrophyte Controls on Urban Stream Microbial Metabolic Activity

Urban rivers worldwide are affected directly by macrophyte growth, causing reduced flow velocity and risks of flooding. Therefore, cutting macrophytes is a common management practice to ensure free drainage. The impacts of macrophyte removal on transient storage dynamics and microbial metabolic activity of wastewater-fed urban streams are unknown, preventing any assessment of the hydrodynamic and biogeochemical consequences of this management practice. Slug tracer injections were performed with the conservative tracer uranine and the reactive tracer resazurin to quantify the implications of macrophyte cutting on stream flow dynamics and metabolism. Macrophyte cutting reduced mean tracer arrival times in managed stream reaches but did not significantly decrease whole-stream microbial metabolic activity. In fact, transient storage indices were found to have increased after cutting, suggesting that macrophyte removal and the resulting increase in flow velocity may have enhanced hyporheic exchange flow through streambed sediments. Our results evidence that macrophyte cutting in nutrient-rich urban streams does not necessarily lead to lower in-stream storage and metabolism but that the gain in hyporheic exchange and streambed microbial metabolic activity can compensate for reduced in-stream storage. Increased stream flow resulting from macrophyte removal may therefore even enhance nutrient and pollutant attenuation capacity of streambed sediments.

Anaerobic degradation of xenobiotic organic contaminants (XOCs): The role of electron flow and potential enhancing strategies

In groundwater, deep soil layer, sediment, the widespread of xenobiotic organic contaminants (XOCs) have been leading to the concern of human health and eco-environment safety, which calls for a better understanding on the fate and remediation of XOCs in anoxic matrices. In the absence of oxygen, bacteria utilize various oxidized substances, e.g. nitrate, sulphate, metallic (hydr)oxides, humic substance, as terminal electron acceptors (TEAs) to fuel anaerobic XOCs degradation. Although there have been increasing anaerobic biodegradation studies focusing on species identification, degrading pathways, community dynamics, systematic reviews on the underlying mechanism of anaerobic contaminants removal from the perspective of electron flow are limited. In this review, we provide the insight on anaerobic biodegradation from electrons aspect - electron production, transport, and consumption. The mechanism of the coupling between TEAs reduction and pollutants degradation is deconstructed in the level of community, pure culture, and cellular biochemistry. Hereby, relevant strategies to promote anaerobic biodegradation are proposed for guiding to an efficient XOCs bioremediation.

Simultaneous attenuation of trace organics and change in organic matter composition in the hyporheic zone of urban streams

Trace organic compounds (TrOCs) enter rivers with discharge of treated wastewater. These effluents can contain high loads of dissolved organic matter (DOM). In a 48 h field study, we investigated changes in molecular composition of seven DOM compound classes (FTICR-MS) and attenuation of 17 polar TrOCs in a small urban stream receiving treated wastewater. Correlations between TrOCs and DOM were used to identify simultaneous changes in surface water and the hyporheic zone. Changes in TrOC concentrations in surface water ranged between a decrease of 29.2% for methylbenzotriazole and an increase of 152.2% for the transformation product gabapentin-lactam. In the hyporheic zone, only decreasing TrOC concentrations were observed, ranging from 4.9% for primidone to 93.8% for venlafaxine . TrOC attenuation coincided with a decline of molecular diversity of easily biodegradable DOM compound classes while molecular diversity of poorly biodegradable DOM compound classes increased. This concurrence indicates similar or linked attenuation pathways for biodegradable DOM and TrOCs. Strong correlations between TrOCs and DOM compound classes as well as high attenuation of TrOCs primarily occurred in the hyporheic zone. This suggests high potential for DOM turnover and TrOC mitigation in rivers if hyporheic exchange is sufficient.

Bacterial Diversity Controls Transformation of Wastewater-Derived Organic Contaminants in River-Simulating Flumes

Hyporheic zones are the water-saturated flow-through subsurfaces of rivers which are characterized by the simultaneous occurrence of multiple physical, biological, and chemical processes. Two factors playing a role in the hyporheic attenuation of organic contaminants are sediment bedforms (a major driver of hyporheic exchange) and the composition of the sediment microbial community. How these factors act on the diverse range of organic contaminants encountered downstream from wastewater treatment plants is not well understood. To address this knowledge gap, we investigated dissipation half-lives (DT50s) of 31 substances (mainly pharmaceuticals) under different combinations of bacterial diversity and bedform-induced hyporheic flow using 20 recirculating flumes in a central composite face factorial design. By combining small-volume pore water sampling, targeted analysis, and suspect screening, along with quantitative real-time PCR and time-resolved amplicon Illumina MiSeq sequencing, we determined a comprehensive set of DT50s, associated bacterial communities, and microbial transformation products. The resulting DT50s of parent compounds ranged from 0.5 (fluoxetine) to 306 days (carbamazepine), with 20 substances responding significantly to bacterial diversity and four to both diversity and hyporheic flow. Bacterial taxa that were associated with biodegradation included Acidobacteria (groups 6, 17, and 22), Actinobacteria (Nocardioides and Illumatobacter), Bacteroidetes (Terrimonas and Flavobacterium) and diverse Proteobacteria (Pseudomonadaceae, Sphingomonadaceae, and Xanthomonadaceae). Notable were the formation of valsartan acid from irbesartan and valsartan, the persistence of N-desmethylvenlafaxine across all treatments, and the identification of biuret as a novel transformation product of metformin. Twelve additional target transformation products were identified, which were persistent in either pore or surface water of at least one treatment, indicating their environmental relevance.

The effect of unsteady streamflow and stream-groundwater interactions on oxygen consumption in a sandy streambed

Streamflow dynamics are often ignored when studying biogeochemical processes in the hyporheic zone. We explored the interactive effects of unsteady streamflow and groundwater fluxes on the delivery and consumption of oxygen within the hyporheic zone using a recirculating flume packed with natural sandy sediments. The flume was equipped with a programmable streamflow control and drainage system that was used to impose losing and gaining fluxes. Tracer tests were used to measure hyporheic exchange flux and a planar optode was used to measure subsurface oxygen concentration patterns. It was found that the volume of the oxic zone decreased when the losing flux declined, and was drastically decreased when gaining conditions were applied. It was also found that unsteady streamflow led to a slight increase in the average volume of the oxic zone, compared to the average volume of the oxic zone under steady streamflow. However, the average oxygen consumption rates were significantly higher under unsteady streamflow compared to steady streamflow under all groundwater conditions with the exception of the highest losing flux. The present study provides the first insight into the interactions between streamflow unsteadiness and losing/gaining fluxes and improve understanding of their impact on microbial metabolism in the hyporheic zone.

Using recirculating flumes and a response surface model to investigate the role of hyporheic exchange and bacterial diversity on micropollutant half-lives

Enhancing the understanding of the fate of wastewater-derived organic micropollutants in rivers is crucial to improve risk assessment, regulatory decision making and river management. Hyporheic exchange and sediment bacterial diversity are two factors gaining increasing importance as drivers for micropollutant degradation, but are complex to study in field experiments and usually ignored in laboratory tests aimed to estimate environmental half-lives. Flume mesocosms are useful to investigate micropollutant degradation processes, bridging the gap between the field and batch experiments. However, few studies have used flumes in this context. We present a novel experimental setup using 20 recirculating flumes and a response surface model to study the influence of hyporheic exchange and sediment bacterial diversity on half-lives of the anti-epileptic drug carbamazepine (CBZ) and the artificial sweetener acesulfame (ACS). The effect of bedform-induced hyporheic exchange was tested by three treatment levels differing in number of bedforms (0, 3 and 6). Three levels of sediment bacterial diversity were obtained by diluting sediment from the River Erpe in Berlin, Germany, with sand (1 : 10, 1 : 1000 and 1 : 100 000). Our results show that ACS half-lives were significantly influenced by sediment dilution and number of bedforms. Half-lives of CBZ were higher than ACS, and were significantly affected only by the sediment dilution variable, and thus by bacterial diversity. Our results show that (1) the flume-setup is a useful tool to study the fate of micropollutants in rivers, and that (2) higher hyporheic exchange and bacterial diversity in the sediment can increase the degradation of micropollutants in rivers.

Passive sampling of organic contaminants across the water-sediment interface of an urban stream

Passive sampling is a well-established tool for monitoring time-weighted average concentrations of polar and semi-polar organic contaminants in streams at flow velocities between 0.1 and 0.4 m s^-1. However, its application under low-flow conditions (10^-5 to 0.01 m s^-1) - as encountered in hyporheic zones - has been scarcely reported. In this study, 3 novel passive sampler configurations were developed for the monitoring of (semi-)polar organic pollutants and related transformation products across the water-sediment interface and thus across varying hydrodynamic conditions. Their design was inspired by Chemcatcher and diffusive gradients in thin films for organics. To determine the most optimal sampler design, an uptake experiment was completed involving the 3 novel passive sampler configurations and a reference Chemcatcher in polar configuration. The experiments consisted of a circular flume that simulated the main channel of a stream and an aquarium with stagnant water that represented the underlying hyporheic zone. The systems were exposed to 192 organic pollutants at environmental concentrations, and the samplers were then collected, extracted and analyzed using liquid chromatography high-resolution mass spectrometry after 2, 6 and 14 days. The configuration that was most insensitive to different hydrodynamic conditions consisted of a reversed-phase sulfonated styrenedivinylbenzene disk as the receiving phase that was covered by an agarose diffusion gel and topped with a polyethersulfone membrane filter. To further evaluate its environmental application, samplers were installed downstream of a sewage treatment plant located at an urban stream in Berlin, Germany (Erpe). The samplers were mounted on custom-made holders which were subsequently embedded in the stream bed to position samplers above (0.30 m) and within the sediment (-0.15/-0.30/-0.45 m) for 11 days. Target and suspect screening workflows were then applied to identify common concentration patterns and link parent attenuation to transformation product formation. A total of 104 concentration profiles were determined, suggesting the efficiency of the proposed sampling strategy in the water-sediment interface. Valsartan acid was the only known transformation product indicative of hyporheic zone-driven attenuation as its concentration in porewater by far exceeded its concentration in surface water. Similar patterns were observed for a larger list of suspected transformation products, of which a sotalol transformation product was tentatively identified. Overall, the established sampling methodology can be effectively used to quantify organic contaminants during low-flow conditions and is suitable for the characterization of attenuation patterns of organic pollutants in hyporheic zones.

Is the Hyporheic Zone Relevant beyond the Scientific Community?

Rivers are important ecosystems under continuous anthropogenic stresses. The hyporheic zone is a ubiquitous, reactive interface between the main channel and its surrounding sediments along the river network. We elaborate on the main physical, biological, and biogeochemical drivers and processes within the hyporheic zone that have been studied by multiple scientific disciplines for almost half a century. These previous efforts have shown that the hyporheic zone is a modulator for most metabolic stream processes and serves as a refuge and habitat for a diverse range of aquatic organisms. It also exerts a major control on river water quality by increasing the contact time with reactive environments, which in turn results in retention and transformation of nutrients, trace organic compounds, fine suspended particles, and microplastics, among others. The paper showcases the critical importance of hyporheic zones, both from a scientific and an applied perspective, and their role in ecosystem services to answer the question of the manuscript title. It identifies major research gaps in our understanding of hyporheic processes. In conclusion, we highlight the potential of hyporheic restoration to efficiently manage and reactivate ecosystem functions and services in river corridors.

Qualitative and quantitative use of micropollutants as source and process indicators. A review

Nowadays, micropollutants such as pharmaceuticals, pesticides and personal care products can be found ubiquitously in the anthropogenically influenced water cycle. As micropollutants have virtually no natural background concentrations they are significantly more sensitive in detecting processes and flow paths than classic inorganic tracers and indicators and at the same time they are often highly source specific. Therefore, using micropollutants as environmental indicators for anthropogenic activities is a common and frequently applied method today. As they interact in many ways with environmental matrices they can be used for source apportionment as well as to estimate flow paths and residence times in waterbodies. This review gives a systematic overview over the large variety of micropollutants used as indicators in the aquatic environment over the last decades together with the prerequisites on their use. Their application is subdivided into their qualitative (compound presence or absence) and quantitative (volume flows) use and shows the numerous possibilities from gaining basic information on the water regime up to advanced applications such as wastewater-based epidemiology.

Temperature-dependent redox zonation, nitrate removal and attenuation of organic micropollutants during bank filtration

River bank filtration (RBF) is considered to efficiently remove nitrate and trace organic micropollutants (OMP) from polluted surface waters. This is essential for maintaining good groundwater quality and providing high quality drinking water. Predicting the fate of OMP during RBF is difficult as the biogeochemical factors controlling the removal efficiency are not fully understood. To determine in-situ removal efficiency and degradation rates of nitrate and OMP indicator substances we conducted a field study in a RBF system during a period of one and a half years incorporating temporally and spatially varying redox conditions and temperature changes typically occurring in temperate climates. RBF was analyzed by means of mixing ratios between infiltrated river water and groundwater as well as average residence times of surface water towards the individual groundwater observation wells. These results were used to calculate temperature dependent first order degradation rates of redox sensitive species and several OMP. Five out of ten investigated OMP were completely removed along RBF pathways. We demonstrate that degradation rates of several OMP during bank filtration were controlled by redox conditions and temperature whereby temperature itself also had a significant influence on the extent of the most reactive oxic zone. The seasonal variations in temperature alone could explain a considerable percentage of the variance in dissolved oxygen (34%), nitrate (81%) as well as the OMPs diclofenac (44%) and sulfamethoxazole (76%). Estimated in-situ degradation rates roughly varied within one order of magnitude for temperature changes between 5 °C and 20 °C. This study highlights that temporal variability in temperature and redox zonation is a significant factor for migration and degradation of nitrate and several OMPs.

Revealing chlorinated ethene transformation hotspots in a nitrate-impacted hyporheic zone

Hyporheic zones are increasingly thought of as natural bioreactors, capable of transforming and attenuating groundwater pollutants present in diffuse baseflow. An underappreciated scenario in the understanding of contaminant fate in hyporheic zones is the interaction between point-source trichloroethene (TCE) plumes and ubiquitous, non-point source pollutants such as nitrate. This study aims to conceptualise critical biogeochemical gradients in the hyporheic zone which govern the export potential of these redox-sensitive pollutants from carbon-poor, oxic aquifers. Within the TCE plume discharge zone, discrete vertical profiling of the upper 100 cm of sediment pore water chemistry revealed an 80% increase in dissolved organic carbon (DOC) concentrations and 20-60 cm thick hypoxic zones (<2 mg O₂ L^-1) within which most reactive transport was observed. A 33% reduction of nitrate concentrations coincided with elevated pore water nitrous oxide concentrations as well as the appearance of manganese and the TCE metabolite cis-1,2-dichloroethene (cDCE). Elevated groundwater nitrate concentrations (>50 mg L^-1) create a large stoichiometric demand for bioavailable DOC in discharging groundwater. With the benefit of a high-resolution grid of pore water samplers investigating the shallowest 30 cm of hypoxic groundwater flow paths, we identified DOC-rich hotspots associated with submerged vegetation (Ranunculus spp.), where low-energy metabolic processes such as mineral dissolution/reduction, methanogenesis and ammonification dominate. Using a chlorine index metric, we show that enhanced TCE to cDCE transformation takes place within these biogeochemical hotspots, highlighting their relevance for natural plume attenuation.

In situ removal of four organic micropollutants in a small river determined by monitoring and modelling

Organic micropollutants (OMPs) are widely detected in surface waters. So far, the removal processes of these compounds in situ in river systems are not yet totally revealed. In this study, a combined monitoring and modelling approach was applied to determine the behaviour of 1-H benzotriazole, carbamazepine, diclofenac and galaxolide in a small river system. Sewage treatment plant effluents and the receiving waters of the river Swist were monitored in 9 dry weather sampling campaigns (precipitation < 1 mm on the sampling day itself and <5 mm total precipitation two days before the sampling) during different seasons over a period of 3 years. With the results gained through monitoring, mass balances have been calculated to assess fate in the river. With the DWA Water Quality Model, OMP concentrations in the river were successfully simulated with OMP characteristics gained through literature studies. No removal was determined for 1-H benzotriazole and carbamazepine, whereas diclofenac showed removal that coincided with light intensity. Moreover, modelling based on light sensitivity of diclofenac also suggested relevant degradation at natural light conditions. These two approaches suggest removal by photodegradation. The highest removal in the river was detected for galaxolide, presumably due to volatilisation, sorption and biodegradation. Furthermore, short-term concentration variability in the river was determined, showing that daily concentration patterns are influenced by dynamics of sewage treatment plant effluent volumes and removal processes in the river.

Nontarget Screening Reveals Time Trends of Polar Micropollutants in a Riverbank Filtration System

The historic emissions of polar micropollutants in a natural drinking water source were investigated by nontarget screening with high-resolution mass spectrometry and open cheminformatics tools. The study area consisted of a riverbank filtration transect fed by the river Lek, a branch of the lower Rhine, and exhibiting up to 60-year travel time. More than 18,000 profiles were detected. Hierarchical clustering revealed that 43% of the 15 most populated clusters were characterized by intensity trends with maxima in the 1990s, reflecting intensified human activities, wastewater treatment plant upgrades and regulation in the Rhine riparian countries. Tentative structure annotation was performed using automated in silico fragmentation. Candidate structures retrieved from ChemSpider were scored based on the fit of the in silico fragments to the experimental tandem mass spectra, similarity to openly accessible accurate mass spectra, associated metadata, and presence in a suspect list. Sixty-seven unique structures (72 over both ionization modes) were tentatively identified, 25 of which were confirmed and included contaminants so far unknown to occur in bank filtrate or in natural waters at all, such as tetramethylsulfamide. This study demonstrates that many classes of hydrophilic organics enter riverbank filtration systems, persisting and migrating for decades if biogeochemical conditions are stable.

Fate of Trace Organic Compounds in the Hyporheic Zone: Influence of Retardation, the Benthic Biolayer, and Organic Carbon

The fate of 28 trace organic compounds (TrOCs) was investigated in the hyporheic zone (HZ) of an urban lowland river in Berlin, Germany. Water samples were collected hourly over 17 h in the river and in three depths in the HZ using minipoint samplers. The four relatively variable time series were subsequently used to calculate first-order removal rates and retardation coefficients via a one-dimensional reactive transport model. Reversible sorption processes led to substantial retardation of many TrOCs along the investigated hyporheic flow path. Some TrOCs, such as dihydroxy-carbamazepine, O-desmethylvenlafaxine, and venlafaxine, were found to be stable in the HZ. Others were readily removed with half-lives in the first 10 cm of the HZ ranging from 0.1 ± 0.01 h for iopromide to 3.3 ± 0.3 h for tramadol. Removal rate constants of the majority of reactive TrOCs were highest in the first 10 cm of the HZ, where removal of biodegradable dissolved organic matter was also the highest. Because conditions were oxic along the top 30 cm of the investigated flow path, we attribute this finding to the high microbial activity typically associated with the shallow HZ. Frequent and short vertical hyporheic exchange flows could therefore be more important for reach-scale TrOC removal than long, lateral hyporheic flow paths.

Spatial and Temporal Variability in Attenuation of Polar Organic Micropollutants in an Urban Lowland Stream

Contamination of rivers by trace organic compounds (TrOCs) poses a risk for aquatic ecosystems and drinking water quality. Spatially- and temporally varying environmental conditions are expected to play a major role in controlling in-stream attenuation of TrOCs. This variability is rarely captured by in situ studies of TrOC attenuation. Instead, snap-shots or time-weighted average conditions and corresponding attenuation rates are reported. The present work sought to investigate this variability and factors controlling it by analysis of 24 TrOCs over a 4.7 km reach of the River Erpe (Berlin, Germany). The factors investigated included sunlight and water temperature as well as the presence of macrophytes. Attenuation rate constants in 48 consecutive hourly water parcels were tracked along two contiguous river sections of different characteristics. Section 1 was less shaded and more densely covered with submerged macrophytes compared to section 2. The sampling campaign was repeated after macrophyte removal from section 1. The findings show, that section 1 generally provided more favorable conditions for both photo- and biodegradation. Macrophyte removal enhanced photolysis of some compounds (e.g., hydrochlorothiazide and diclofenac) while reducing the biodegradation of metoprolol. The transformation products metoprolol acid and valsartan acid were formed along the reach under all conditions.

Evaluating emerging organic contaminant removal in an engineered hyporheic zone using high resolution mass spectrometry

The hyporheic zone (HZ), located at the interface of surface and groundwater, is a natural bioreactor for attenuation of chemical contaminants. Engineered HZs can be incorporated into stream restoration projects to enhance hyporheic exchange, with flowpaths optimized to promote biological habitat, water quantity, and water quality improvements. Designing HZs for in-stream treatment of stormwater, a significant source of flow and contaminant loads to urban creeks, requires assessment of both the hydrology and biogeochemical capacity for water quality improvement. Here, we applied tracer tests and high resolution mass spectrometry (HRMS) to characterize an engineered hyporheic zone unit process, called a hyporheic design element (HDE), in the Thornton Creek Watershed in Seattle, WA. Dye, NaCl, and bromide were used to hydrologically link downwelling and upwelling zones and estimate the hydraulic retention time (HRT) of hyporheic flowpaths. We then compared water quality improvements across hydrologically-linked surface and hyporheic flowpaths (3-5 m length; ∼30 min to >3 h) during baseflow and stormflow conditions. We evaluated fate outcomes for 83 identified contaminants during stormflow, including those correlated with an urban runoff mortality syndrome in coho salmon. Non-target HRMS analysis was used to assess holistic water quality improvements and evaluate attenuation mechanisms. The data indicated substantial water quality improvement in hyporheic flowpaths relative to surface flow and improved contaminant removal with longer hyporheic HRT (for ∼1900 non-target compounds detected during stormflow, <17% were attenuated >50% via surface flow vs. 59% and 78% via short and long hyporheic residence times, respectively), and strong contributions of hydrophobic sorption towards observed contaminant attenuation.

Determination of polar organic micropollutants in surface and pore water by high-resolution sampling-direct injection-ultra high performance liquid chromatography-tandem mass spectrometry

Hyporheic zones (HZs) are dynamic and complex transition regions between rivers and aquifers which are thought to play an important role in the attenuation of environmental micropollutants. Non-steady state and small-scale hyporheic processes which affect micropollutants in the HZ are poorly characterized due to limitations in existing analytical methodologies. In this work we developed a method for high spatio-temporal resolution analysis of polar organic micropollutants (POMs) in hyporheic pore- and surface waters by combining (semi-) automatic low volume sampling techniques with direct-injection ultra-high performance liquid chromatography tandem mass spectrometry. The method is capable of quantifying 25 parent compounds and 18 transformation products (TPs) using only 0.4 mL of water and few preparation steps. Application of the method to both surface and pore water revealed significant (i.e. > an order of magnitude) differences in POM concentrations over small time and spatial scales (i.e. < a few hours and tens of cm, respectively). Guanylurea, a TP of the antidiabetic drug metformin was detected at unprecedentedly high concentrations. Collectively, this method is suitable for in situ characterization of POMs at high spatial and temporal resolution and with minimal disturbance of natural flow paths and infiltration of surface water.

Hyporheic Exchange Controls Fate of Trace Organic Compounds in an Urban Stream

First-order half-lives for 26 trace organic compounds (TrOCs) were determined in the hyporheic zone (HZ) and along a 3 km reach of a first-order stream in South Australia during both dry and wet seasons. Two salt tracer experiments were conducted and evaluated using a transient storage model to characterize seasonal differences in stream residence time and transient storage. Lagrangian and time-integrated surface water sampling were conducted to calculate half-lives in the surface water. Half-lives in the HZ were calculated using porewater samples obtained from a modified mini-point sampler and hyporheic residence times measured via active heat-pulse sensing. Half of the investigated TrOCs (e.g., oxazepam, olmesartan, candesartan) were not significantly removed along both the investigated river stretch and the sampled hyporheic flow paths. The remaining TrOCs (e.g., metformin, guanylurea, valsartan) were found to be significantly removed in the HZ and along the river stretch with relative removals in the HZ correlating to reach-scale relative removals. Using the modeled transport parameters, it was estimated that wet season reach-scale removal of TrOCs was predominately caused by removal in the HZ when the intensity of hyporheic exchange was also higher. Factors that increase HZ exchange are thus likely to promote in-stream reactivity of TrOCs.