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

Highly sensitive passive sampling of new pollutants in urban reclaimed water using hydrophilic-lipophilic balance sorbent-embedded cellulose acetate membrane

Urban reclaimed water is important water resource in China, while the residual new pollutants in the water largely challenge their safety for further use. China's action program for the management of new pollutants (also known as emerging contaminants) requires effective method for monitoring diverse new pollutants in water. This work proposed a highly sensitive passive sampling method for monitoring diverse new pollutants in urban reclaimed water. Hydrophilic-lipophilic balance sorbent-embedded cellulose acetate membranes (HECAMs) were dynamically deployed in self-developed continuous flow integrative sampling device (CFISD) for sampling four types of new pollutants with wide polarity ranges (1.11 < log Kow < 9.49) in a reclaimed water network for landscape irrigation in Beijing, China. The estimated equilibrium partition ratios of most chemicals between HECAM and water were over 104, which would provide low detection limits. In the 7-d high-efficiency deployment, thirty new pollutants were detected, which indicated incompletely removal of various new pollutants in wastewater treatment process. The dynamical accumulation data were successfully fitted by first-order kinetic model and different contaminants reached different accumulation phases in the HECAMs during 7 d. Acceptable and steady uptake rate constants and sampling rates were obtained with the use of CFISD in field sampling. The estimated time-weighted average concentrations for contaminants had wide concentration range and were from 0.03 ng L-1 (pendimethalin) to 3,394 ng L-1 (dibutyl phthalate) and this may lead to environmental risk for further use the reclaimed water. Dynamical deployment results also provided sensitive information of concentration fluctuations and twelve pollutants showed concentration fluctuations during the sampling period. In summary, HECAM showed high sensitivities and credible performance of monitoring diverse new pollutants in the urban reclaimed water.

Mass variations and transfer process of shrimp farming pollutants in aquaculture drainage systems: Effects of DOM features and physicochemical properties

The expansion of aquaculture produces increasing pollutant loads, necessitating the use of drainage systems to discharge wastewater into surface water. To assess the mass variations and transfer process of aquaculture wastewater, an entire aquaculture drainage investigation lasting for 48 h was conducted, focusing on the nutrients, heavy metals, dissolved organic matter (DOM), and physicochemical properties of drainage in a commercial shrimp farm. The findings revealed that early drainage produced more heavy metals, total nitrogen (TN), dissolved organic nitrogen (DON), and feed-like proteins from aquaculture floating feed and additives, whereas late drainage produced more PO43--P and total dissolved phosphorus (TP). A few pollutants, including DON, Cu, and feed-like proteins, were effectively removed, whereas the contents of TN, dissolved inorganic nitrogen, and Zn increased in the multi-level aquaculture drainage system. Limited dilution indicated that in-stream transfer was the main process shaping pollutant concentrations within the drainage system. In the lower ditches, NO3--N, heavy metals, and feed-like proteins exhibited evident in-stream attenuation, while TN and NH4+-N underwent significant in-stream enrichment processes, especially in ditch C, with the transfer coefficient values (vf) of -1.74E-5 and -2.04E-5. This indicates that traditional aquaculture drainage systems serve as nitrogen sinks, rather than efficient nutrient purge facilitators. Notably, DOM was identified as a more influential factor in shaping the in-stream transfer process in aquaculture drainage systems, with an interpretation rate 40.79% higher than that of the physiochemical properties. Consequently, it is necessary to eliminate the obstacles posed by DOM to pollutant absorption and net zero emissions in aquaculture drainage systems in the future. ENVIRONMENTAL IMPLICATIONS: Nutrients, heavy metals, and dissolved organic matter are hazardous pollutants originating from high-density aquaculture. As the sole conduit to natural waters, aquaculture drainage systems have pivotal functions in receiving and purifying wastewater, in which the in-stream transfer process is affected by ambient conditions. This field study investigated the spatial variations, stage distinctions, effects of physicochemical properties, and dissolved organic matter (DOM) features. This finding suggests that the aquaculture drainage system as a nitrogen sink and DOM source. While the DOM is the key factor in shaping the in-stream transfer process, and obstacles for pollutant elimination. This study helps in understanding the fate of aquaculture pollutants and reveals the drawbacks of traditional aquaculture drainage systems.

Relevance of photocatalytic redox transformations of selected pharmaceuticals in a copper- and iron-rich Mediterranean intermittent river

This work aimed at investigating specific attenuation pathways of pharmaceuticals in copper- and iron-rich Mediterranean intermittent and sunlit rivers by combining lab- and field-scale studies. Poorly photodegradable and biodegradable compounds such as fluconazole, oxazepam and venlafaxine attenuated in two river stretches with short hydraulic residence times (<3 h). This result was assumed to be related to their capacity to interact with photoreactive free Cu2+ and Fe3+ or their associated oxides. Lab-scale photodegradation experiments under simulated solar irradiation revealed the beneficial impact of a mixture Cu2+ and colloidal iron hydroxides at environmental concentrations and at neutral pH on the pharmaceuticals photodegradation kinetic rate constants. These latter were consistent with the in-stream attenuation rate constants of targeted contaminants which ranged from 0.104 to 0.154 h-1. Further identification of phototransformation products by LC-HRMS highlighted reductive transformation pathways including reductive dehalogenation and hydrogenation reactions. Several TPs were found to be stable under irradiation and were detected in field monitoring, accordingly. This was ascribed to the formation of a Cu/Fe composite material under solar irradiation with photocatalytic properties. The role of Cu was to trap the electron in the conduction band of the iron-based photocatalyst, which promoted separation efficiency of electron-hole pairs as well as enhanced photoreduction processes at the expense of oxidation ones. Even though, these mechanisms have been reported in water treatment field for organic micropollutants removal, their significance was demonstrated for the first time in natural settings.

Analysis on the attenuation characteristics of PPCPs in surface water and their influencing factors based on a compilation of literature data

The attenuation characteristics of PPCPs play an important part in predicting their environmental concentrations. However, considerable uncertainty remains in reported laboratory data on the attenuation characteristics of PPCPs. In this analysis, we compile information on laboratory-observed photodegradation half-lives (t1/2), biodegradation t1/2, the organic carbon normalized adsorption constant (KOC) and field-observed overall attenuation t1/2 for PPCPs in water bodies from more than 200 peer-reviewed studies. To mitigate the effects of such uncertainty, we derive representative values (RV) for PPCP degradability from these records to better compare the characteristics of different PPCPs. We further examine the influence of experimental conditions and environmental drivers on the determination of t1/2 using difference analysis and correlation analysis. The results indicate that for laboratory photodegradation tests, different light sources, initial concentration and volume significantly affect t1/2, whereas there is no significant difference between values obtained from tests conducted in pure water and natural water. For biodegradation, laboratory-measured t1/2 values in batch, flume and column studies gradually decrease, marking the controlling role of experimental setup. Redox condition, initial concentration and volume are also recognized as important influencing factors. For adsorption, water-sediment ratio is the primary reaction parameter. As two frequently investigated factors, however, pH and temperature are not significant factors in almost all cases. In field observations, the persistence of carbamazepine, typically used as a tracer, is in doubt. Water depth and latitude are the most correlated drivers of t1/2, indicating the predominant status of photodegradation in the overall attenuation rates. These findings call for caution when selecting experimental parameters and environmental drivers in determining PPCP's attenuation rates and establishing PPCP fate models in the field.

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.

Do biotransformation data from laboratory experiments reflect micropollutant degradation in a large river basin?

Identifying a chemical's potential for biotransformation in the aquatic environment is crucial to predict its fate and manage its potential hazards. Due to the complexity of natural water bodies, especially river networks, biotransformation is often studied in laboratory experiments, assuming that study outcomes can be extrapolated to compound behavior in the field. Here, we investigated to what extent outcomes of laboratory simulation studies indeed reflect biotransformation kinetics observed in riverine systems. To determine in-field biotransformation, we measured loads of 27 wastewater treatment plant effluent-borne compounds along the Rhine and its major tributaries during two seasons. Up to 21 compounds were detected at each sampling location. Measured compound loads were used in an inverse model framework of the Rhine river basin to derive k'_bio,field values - a compound-specific parameter describing the compounds' average biotransformation potential during the field studies. To support model calibration, we performed phototransformation and sorption experiments with all the study compounds, identifying 5 compounds that are susceptible towards direct phototransformation and determining K_oc values covering four orders of magnitude. On the laboratory side, we used a similar inverse model framework to derive k'_bio,lab values from water-sediment experiments run according to a modified OECD 308-type protocol. The comparison of k'_bio,lab and k'_bio,field revealed that their absolute values differed, pointing towards faster transformation in the Rhine river basin. Yet, we could demonstrate that relative rankings of biotransformation potential and groups of compounds with low, moderate and high persistence agree reasonably well between laboratory and field outcomes. Overall, our results provide evidence that laboratory-based biotransformation studies using the modified OECD 308 protocol and k'_bio values derived thereof bear considerable potential to reflect biotransformation of micropollutants in one of the largest European river basins.

Persistence of N-oxides transformation products of tertiary amine drugs at lab and field studies

This work aimed at studying the formation and persistence of N-oxides transformation products (TPs) of tertiary amine drugs by combining laboratory and field studies relevant for surface water. A monitoring study using passive samplers was first achieved for assessing attenuation of selected pharmaceuticals and their related N-oxides and N-, O-dealkylated TPs (i.e., venlafaxine, tramadol, amisulpride and sulpiride) along a 1.7 km river stretch between two sampling sites. This study revealed the stability of tramadol-N-oxide, amisulpride-N-oxide and the fast dissipation of O-desmethylvenlafaxine-N-oxide, as well as the significance of N-oxidized TPs in comparison to N-dealkylated TPs and parent compounds in river. Lab-scale experiments were then implemented for a better understanding of their mechanisms of formation and degradation under aerobic water/sediment testing and under simulated solar photochemistry. N-oxidation reactions were always a minor transformation pathway under both degradation conditions with respect to N-and O-dealkylation reactions. The amount of generated N-oxides were similar for venlafaxine, tramadol and sulpiride and peaked in the 8.4-12.8% and <4% of their initial concentration (100 μg/L), during photodegradation and biodegradation experiments, respectively. Other transformation pathways such as hydroxylation and α-C-hydroxylation followed by oxidation to amide or dehydration were also identified. Investigated N-oxides TPs (except O-desmethylvenlafaxine-N-oxide) were found stable under solar photolysis and aerobic biodegradation with a very slight reverse reaction to parent compound observed for tramadol-N-oxide and amisulpride-N-oxide. Lab-scale degradation experiments were not able to anticipate the high occurrence levels of N-oxide compounds in the environment. This was most likely due to faster degradation kinetics and/or higher sorption to sediment of parent compounds and dealkylated TPs over N-oxide TPs, resulting in higher relative accumulation of the latter.

Modeling the Dynamics of Mixture Toxicity and Effects of Organic Micropollutants in a Small River under Unsteady Flow Conditions

The presence of anthropogenic organic micropollutants in rivers poses a long-term threat to surface water quality. To describe and quantify the in-stream fate of single micropollutants, the advection-dispersion-reaction (ADR) equation has been used previously. Understanding the dynamics of the mixture effects and cytotoxicity that are cumulatively caused by micropollutant mixtures along their flow path in rivers requires a new concept. Thus, we extended the ADR equation from single micropollutants to defined mixtures and then to the measured mixture effects of micropollutants extracted from the same river water samples. Effects (single and mixture) are expressed as effect units and toxic units, the inverse of effect concentrations and inhibitory concentrations, respectively, quantified using a panel of in vitro bioassays. We performed a Lagrangian sampling campaign under unsteady flow, collecting river water that was impacted by a wastewater treatment plant (WWTP) effluent. To reduce the computational time, the solution of the ADR equation was expressed by a convolution-based reactive transport approach, which was used to simulate the dynamics of the effects. The dissipation dynamics of the individual micropollutants were reproduced by the deterministic model following first-order kinetics. The dynamics of experimental mixture effects without known compositions were captured by the model ensemble obtained through Bayesian calibration. The highly fluctuating WWTP effluent discharge dominated the temporal patterns of the effect fluxes in the river. Minor inputs likely from surface runoff and pesticide diffusion might contribute to the general effect and cytotoxicity pattern but could not be confirmed by the model-based analysis of the available effect and chemical data.

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.

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.

Modeling Risk Dynamics of Contaminants of Emerging Concern in a Temperate-region Wastewater Effluent-dominated Stream

Wastewater effluent-dominated streams are becoming increasingly common worldwide, including in temperate regions, with potential impacts on ecological systems and drinking water sources. We recently quantified the occurrence/ spatiotemporal dynamics of pharmaceutical mixtures in a representative temperate-region wastewater effluent-dominated stream (Muddy Creek, Iowa) under baseflow conditions and characterized relevant fate processes. Herein, we quantified the ecological risk quotients (RQs) of 19 effluent-derived contaminants of emerging concern (CECs; including: 14 pharmaceuticals, 2 industrial chemicals, and 3 neonicotinoid insecticides) and 1 run-off-derived compound (atrazine) in the stream under baseflow conditions, and estimated the probabilistic risks of effluent-derived CECs under all-flow conditions (i.e., including runoff events) using stochastic risk modeling. We determined that 11 out of 20 CECs pose medium-to-high risks to local ecological systems (i.e., algae, invertebrates, fish) based on literature-derived acute effects under measured baseflow conditions. Stochastic risk modeling indicated decreased, but still problematic, risk of effluent-derived CECs (i.e., RQ≥0.1) under all-flow conditions when runoff events were included. Dilution of effluent-derived chemicals from storm flows thus only minimally decreased risk to aquatic biota in the effluent-dominated stream. We also modeled in-stream transport. Thirteen out of 14 pharmaceuticals persisted along the stream reach (median attenuation rate constant k<0.1 h^-1) and entered the Iowa River at elevated concentrations. Predicted and measured concentrations in the drinking water treatment plant were below the human health benchmarks. This study demonstrates the application of probabilistic risk assessments for effluent-derived CECs in a representative effluent-dominated stream under variable flow conditions (when measurements are less practical) and provides an enhanced prediction tool transferable to other effluent-dominated systems.

Lignin peroxidase-catalyzed direct oxidation of trace organic pollutants through a long-range electron transfer mechanism: Using propranolol as an example

In this work, lignin peroxidase (LiP) was extracted for the in vitro degradation of a persistent compound (propranolol, PPN). The results showed that 94.2% of PPN was degraded at 30 U L^-1 LiP activity and 10 mg L^-1 PPN. The PPN degradation rate increased from 33.5% to 94.2% when the veratryl alcohol (VA) concentration varied from 0 to 180 µM, but decreased to 73.1% with further VA addition. This phenomenon confirmed that VA was indispensable, however, it also acted as a competitive inhibitor of PPN oxidation. Computational analysis revealed that the Trp171…iron porphyrin (TRP-FeP) path was responsible for specific substrate (e.g., VA) transformation, and another long-range electron transfer (LRET) path through His-Asp…FeP (HSP-FeP) was discovered for non-specific substrate (e.g., PPN) degradation. These two electron-transfer routes shared one catalytic center, and VA protected the enzyme from H₂O₂-dependent inactivation. The HSP-FeP path transformed PPN through single electron transfer or H abstraction mechanisms. In addition, hydroxyl radicals generated in the LiP/H₂O₂ system were involved in the hydroxylation of the PPN intermediates. Possible degradation pathways were deduced using these degradation mechanisms and mass-spectrometry analysis. The multipath degradation mechanism endowed LiP with a remarkable capacity for removing various recalcitrant pollutants in environmental remediation.

Sedimentary spectrum and potential ecological risks of residual pharmaceuticals in relation to sediment-water partitioning and land uses in a watershed

Pharmaceutical residues in river surficial sediment are prone to anthropogenic impacts and environmental factors in watershed, but the mechanisms remain unclear. This study attempted to reveal surficial sediment-water pseudo-partitioning and anthropogenic (land use) patterns of pharmaceutical residues in surficial sediment among 23 subwatersheds of Jiulong River, southeast China with a gradient of urban land use percentile in dry and wet seasons. Thirty-eight out of target 86 compounds from six-category pharmaceuticals were quantified and ranged from below the quantification limits (0.001 mg kg^-1 dry mass) up to 8.19 mg kg^-1 dry mass (chlortetracycline) using a developed SPE-HPLC-MS/MS protocol. Antibiotics and non-steroidal anti-inflammatory drugs (NSAIDs) collectively dominated sedimentary pharmaceutical residues for 34.5-99.8% of the total quantified compounds (median at 92%). Land uses in subwatersheds showed high consistency with sedimentary pharmaceutical residues in the dry season rather than the wet season, especially for human use only and veterinary use only compounds. Surficial sediment-water partitioning of pharmaceutical compounds influenced their sedimentary residues regardless of season, which were determined by properties of compound and surficial sediment interactively. All tetracycline compounds, trimethoprim (sulfonamides synergist), caffeine (central nervous system drug), and oxfendazole (antiparasitic drug) were quantified to pose high potential ecological risks to aquatics. Findings of this study suggest that pseudo-persistent legacy of human and veterinary pharmaceuticals requires a wider coverage of pharmaceutical compounds for a comprehensive ecological assessment in the environment and more involvement of anthropogenic impacts and socioeconomic factors in the future studies.

Analytical method to monitor contaminants of emerging concern in water and soil samples from a non-conventional wastewater treatment system

Nonconventional wastewater treatments, such as vegetation filters (VFs), are propitious systems to attenuate contaminants of emerging concern (CECs) in small municipalities. The development of standardised multiresidue and multimatrix methods suitable for measuring a reliable number of CEC in environmental samples is crucial for monitoring infiltrating concentrations and for ensuring these systems' treatment capacity. The objective of this study is to develop and validate an analytical method for the simultaneous determination of CECs, including transformation products (TPs), with diverse physico-chemical properties, in environmental samples. The optimised method is based on sample clean-up and preconcentration by solid-phase extraction (SPE), followed by liquid chromatography electrospray ionization tandem mass spectrometry (LC-MS/MS). The method is able to detect and quantify 40 target CECs, including pharmaceuticals of different classes (analgesics, antibiotics, antihypertensives, lipid regulators, anticonvulsants, antidepressants, antiarrhythmics, beta-blockers, amongst others), hormones and lifestyle products with good reproducibility (variations below 23%), in different water matrices, and 28 CECs, in soil samples. Acceptable recoveries (65-120%) were obtained for most of the CECs in all the matrices. However in the soil samples, as complexity required a prior extraction treatment, the recovery of some analytes was affected, which reduced the number of target CECs. The achieved methodological quantification limits (0.05-5 ng/L and 0.04-1.1 ng/g levels for the water and the soil matrices, respectively) were reasonably low for most CECs. The proposed method was successfully applied to monitor CECs in a VF. The CECs detected at higher concentrations are some of the world's most widely used products (e.g. acetaminophen or caffeine and its main TP, paraxanthine). The results showed an almost 70% reduction in CEC concentrations during infiltration. The groundwater data indicated that the VF treatment operation did not affect the underlying aquifer (Cmax found in GW <1 µg/L).

Coupling Suspect and Nontarget Screening with Mass Balance Modeling to Characterize Organic Micropollutants in the Onondaga Lake-Three Rivers System

Characterizing the occurrence, sources, and fate of organic micropollutants (OMPs) in lake-river systems serves as an important foundation for constraining the potential impacts of OMPs on the ecosystem functions of these critical landscape features. In this work, we combined suspect and nontarget screening with mass balance modeling to investigate OMP contamination in the Onondaga Lake-Three Rivers system of New York. Suspect and nontarget screening enabled by liquid chromatography-high-resolution mass spectrometry led to the confirmation and quantification of 105 OMPs in water samples collected throughout the lake-river system, which were grouped by their concentration patterns into wastewater-derived and mixed-source clusters via hierarchical cluster analysis. Four of these OMPs (i.e., galaxolidone, diphenylphosphinic acid, N-butylbenzenesulfonamide, and triisopropanolamine) were prioritized and identified by nontarget screening based on their characteristic vertical distribution patterns during thermal stratification in Onondaga Lake. Mass balance modeling performed using the concentration and discharge data highlighted the export of OMPs from Onondaga Lake to the Three Rivers as a major contributor to the OMP budget in this lake-river system. Overall, this work demonstrated the utility of an integrated screening and modeling framework that can be adapted for OMP characterization, fate assessment, and load apportionment in similar surface water systems.

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.

Tandem field and laboratory approaches to quantify attenuation mechanisms of pharmaceutical and pharmaceutical transformation products in a wastewater effluent-dominated stream

Evolving complex mixtures of pharmaceuticals and transformation products in effluent-dominated streams pose potential impacts to aquatic species; thus, understanding the attenuation dynamics in the field and characterizing the prominent attenuation mechanisms of pharmaceuticals and their transformation products (TPs) is critical for hazard assessments. Herein, we determined the attenuation dynamics and the associated prominent mechanisms of pharmaceuticals and their corresponding TPs via a combined long-term field study and controlled laboratory experiments. For the field study, we quantified spatiotemporal exposure concentrations of five pharmaceuticals and six associated TPs in a small, temperate-region effluent-dominated stream during baseflow conditions where the wastewater plant was the main source of pharmaceuticals. We selected four sites (upstream, at, and two progressively downstream from effluent discharge) and collected water samples at 16 time points (64 samples in total, approximately twice monthly, depending on flows) for 1 year. Concurrently, we conducted photolysis, sorption, and biodegradation batch tests under controlled conditions to determine the major attenuation mechanisms. We observed 10-fold greater attenuation rates in the field compared to batch tests, demonstrating that connecting laboratory batch tests with field measurements to enhance predictive power is a critical need. Batch systems alone, often used for assessment, are useful for determining fate processes but poorly approximate in-stream attenuation kinetics. Sorption was the dominant attenuation process (t_1/2<7.7 d) for 5 of 11 compounds in the batch tests, while the other compounds (n = 6) persisted in the batch tests and along the 5.1 km stream reach. In-stream parent-to-product transformation was minimal. Differential attenuation contributed to the evolving pharmaceutical mixture and created changing exposure conditions with concomitant implications for aquatic and terrestrial biota. Tandem field and laboratory characterization can better inform modeling efforts for transport and risk assessments.

Dosing the Coast: Leaking Sewage Infrastructure Delivers Large Annual Doses and Dynamic Mixtures of Pharmaceuticals to Urban Rivers

Pharmaceuticals are commonly detected at low concentrations in surface waters, where they disrupt biological and ecological processes. Despite their ubiquity, the annual mass of pharmaceuticals exported from watersheds is rarely quantified. We used liquid chromatography-mass spectroscopy to screen for 92 pharmaceuticals in weekly samples from an urban stream network in Baltimore, MD, USA, that lacks wastewater treatment effluents. Across the network, we detected 37 unique compounds, with higher concentrations and more compounds in streams with higher population densities. We also used concentrations and stream discharge to calculate annual pharmaceutical loads at the watershed outlet, which range from less than 1 kg to ∼15 kg and are equivalent to tens of thousands of human doses. By calculating annual watershed mass balances for eight compounds, we show that ∼0.05 to ∼42% of the pharmaceuticals consumed by humans in this watershed are released to surface waters, with the importance of different pathways (leaking sewage vs treated wastewater effluent) differing among compounds. These results demonstrate the importance of developing, maintaining, and improving sewage infrastructure to protect water resources from pharmaceutical contamination.

Manganese peroxidase mediated oxidation of sulfamethoxazole: Integrating the computational analysis to reveal the reaction kinetics, mechanistic insights, and oxidation pathway

In this study, manganese peroxidase (MnP) was applied to induce the in vitro oxidation of sulfamethoxazole (SMX). The results indicated that 87.04% of the SMX was transformed and followed first-order kinetics (k_obs=0.438 h^-1) within 6 h when 40 U L^-1 of MnP was added. The reaction kinetics were investigated under different conditions, including pH, MnP activity, and H₂O₂ concentration. The active species Mn³⁺ was responsible for the oxidation of SMX, and the Mn³⁺ production rate was monitored to reveal the interaction among MnP, Mn³⁺, and SMX. By integrating the characterizations analysis of the MnP/H₂O₂ system with the density functional theory (DFT) calculations, the proton-coupled electron transfer (PCET) process dominated the catalytic circle of MnP and the transformation of Mn³⁺. Additionally, possible oxidation pathways of SMX were proposed based on single-electron transfer mechanism, which primarily included the S-N bond cleavage, the C-S bond cleavage, and one electron loss without bond breakage. It was then transformed to hydrolysis, N-H oxidation, self-coupling, and carboxylic acid coupling products. This study provides insights into the atomic-level mechanism of MnP and the transformation pathways of sulfamethoxazole, which lays a significant foundation for the potential of MnP in wastewater treatment applications.

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.

Separation of Photochemical and Non-Photochemical Diurnal In-Stream Attenuation of Micropollutants

For a better process understanding of in-stream attenuation of trace organic contaminants (TrOCs), quantitative comparisons between field studies under different environmental conditions and controlled laboratory experiments are important to separate different processes. However, this is hampered by the challenge to transfer kinetics from the laboratory to different field conditions due to the lack of good quantitative measures to account for different boundary conditions. For phototransformation, in situ light conditions in a river are difficult to determine because light is reduced, for instance, by absorption, scattering on suspended particles, and shading effects. In this study, we present an approach to separate photochemical from non-photochemical diurnal in-stream attenuation based on rate constants relative to diclofenac, as a reference compound, to account for the difference in the in situ light conditions combined with laboratory experiments. 12 out of 45 detected target TrOCs showed a diurnal attenuation at a selected river stretch. A non-photochemical process, potentially biotransformation, was responsible for the diurnal attenuation of bisoprolol, metoprolol, O-desmethylvenlafaxine, tramadol, and venlafaxine. Attenuation of amisulpride, flufenamic acid, hydrochlorothiazide, naproxen, and xipamide can be quantitatively explained by phototransformation, partially for sotalol. Attenuation rate constants of hydrochlorothiazide at different field sites from this study and from published data range over 2 orders of magnitude. Differences can be quantitatively explained by different light exposures but not by water chemical parameters.

Transformation of organic micropollutants along hyporheic flow in bedforms of river-simulating flumes

Urban streams receive increasing loads of organic micropollutants from treated wastewaters. A comprehensive understanding of the in-stream fate of micropollutants is thus of high interest for water quality management. Bedforms induce pumping effects considerably contributing to whole stream hyporheic exchange and are hotspots of biogeochemical turnover processes. However, little is known about the transformation of micropollutants in such structures. In the present study, we set up recirculating flumes to examine the transformation of a set of micropollutants along single flowpaths in two triangular bedforms. We sampled porewater from four locations in the bedforms over 78 days and analysed the resulting concentration curves using the results of a hydrodynamic model in combination with a reactive transport model accounting for advection, dispersion, first-order removal and retardation. The four porewater sampling locations were positioned on individual flowpaths with median solute travel times ranging from 11.5 to 43.3 h as shown in a hydrodynamic model previously. Highest stability was estimated for hydrochlorothiazide on all flowpaths. Lowest detectable half-lives were estimated for sotalol (0.7 h) and sitagliptin (0.2 h) along the shortest flowpath. Also, venlafaxine, acesulfame, bezafibrate, irbesartan, valsartan, ibuprofen and naproxen displayed lower half-lives at shorter flowpaths in the first bedform. However, the behavior of many compounds in the second bedform deviated from expectations, where particularly transformation products, e.g. valsartan acid, showed high concentrations. Flowpath-specific behavior as observed for metformin or flume-specific behavior as observed for metoprolol acid, for instance, was attributed to potential small-scale or flume-scale heterogeneity of microbial community compositions, respectively. The results of the study indicate that the shallow hyporheic flow field and the small-scale heterogeneity of the microbial community are major controlling factors for the transformation of relevant micropollutants in river sediments.

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.

Assessment of UV combined with free chlorine for removal of valsartan acid from water samples

Valsartan acid (VALA) is a persistent and mobile pollutant, ubiquitously distributed in the aquatic environment. Herein, we assessed the efficiency of UV/free chlorine for the removal of this pollutant. Degradation experiments were performed using different water samples, considering several pH values and concentrations of inorganic anions. Time-course of VALA was measured by injection of different reaction time aliquots in a liquid chromatography (LC) triple quadrupole (QqQ) mass spectrometry (MS) system, whilst the study of potential transformation products (TPs) was evaluated by LC combined with a hybrid quadrupole time-of-flight (QTOF) MS system. Formation of volatile disinfection by-products (DBPs) was investigated by gas chromatography (GC) with TOF-MS detection. Compared to free chlorine treatment and UV photolysis, the combination of both parameters significantly enhanced the degradability of VALA. At neutral pH, UV/free chlorine was also more effective than UV/H₂O₂ to remove VALA from spiked water solutions. Three TPs of VALA were tentatively identified by LC-QTOF-MS, although only one was stable in the UV/free chlorine media. As regards volatile DBPs, the formation of chloroform, dichloroacetonitrile, di- and trichloroacetic acid was noticed. The mass yield of DBPs formation from VALA varied from 0.3% (dichloroacetonitrile) to 1.1% (chloroform). The efficiency of UV/free chlorine was first investigated in spiked solutions with increasing complexities: ultrapure, river and treated wastewater. Thereafter, the feasibility of reducing VALA levels in polluted river water was demonstrated.

Revealing the hydrological transport and attenuation of 14 antibiotics in a low-flow stream

The fate of antibiotics has mostly been studied in lab-scale systems. Few studies have precisely evaluated the attenuation of antibiotics in natural streams. We used Lagrangian sampling combined with a tracer test and one-dimensional transport with inflow and storage model (OTIS) to reveal the effects of dilution and self-attenuation processes, and diurnal variation of light on the degradation of 14 antibiotics in a 3.6 km low-flow stream. The results showed that the order of in-stream attenuation rates were macrolides (0.18-0.25 h^-1) > tetracyclines (0.16-0.18 h^-1) > fluoroquinolones (0.094-0.13 h^-1) > sulfonamides (0.056-0.082 h^-1); half of the mass of antibiotics were lost within 0.44-1.96 km. The dilution effect, including longitudinal dispersion and transient storage, accounted for 20.70%-91.60% of the total attenuation while self-attenuation processes accounted for 8.40%-79.30%. Over 60% of sulfonamides were dissipated by dilution, while over 68% of the removal of macrolides and tetracyclines was attributed to self-attenuation. A comparison of the attenuation rates between day and night demonstrated that photo-dependent attenuation played a dominant role, especially for sulfonamides, accounting for more than 50% of their self-attenuation. Photo-independent attenuation reduced most macrolides and tetracyclines. This in situ experiment increased our understanding of antibiotic attenuation in natural streams.

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.

Fate of Trace Organic Compounds in Hyporheic Zone Sediments of Contrasting Organic Carbon Content and Impact on the Microbiome

The organic carbon in streambed sediments drives multiple biogeochemical reactions, including the attenuation of organic micropollutants. An attenuation assay using sediment microcosms differing in the initial total organic carbon (TOC) revealed higher microbiome and sorption associated removal efficiencies of trace organic compounds (TrOCs) in the high-TOC compared to the low-TOC sediments. Overall, the combined microbial and sorption associated removal efficiencies of the micropollutants were generally higher than by sorption alone for all compounds tested except propranolol whose removal efficiency was similar via both mechanisms. Quantitative real-time PCR and time-resolved 16S rRNA gene amplicon sequencing revealed that higher bacterial abundance and diversity in the high-TOC sediments correlated with higher microbial removal efficiencies of most TrOCs. The bacterial community in the high-TOC sediment samples remained relatively stable against the stressor effects of TrOC amendment compared to the low-TOC sediment community that was characterized by a decline in the relative abundance of most phyla except Proteobacteria. Bacterial genera that were significantly more abundant in amended relative to unamended sediment samples and thus associated with biodegradation of the TrOCs included Xanthobacter, Hyphomicrobium, Novosphingobium, Reyranella and Terrimonas. The collective results indicated that the TOC content influences the microbial community dynamics and associated biotransformation of TrOCs as well as the sorption potential of the hyporheic zone sediments.

Association between Aquatic Micropollutant Dissipation and River Sediment Bacterial Communities

Assessment of micropollutant biodegradation is essential to determine the persistence of potentially hazardous chemicals in aquatic ecosystems. We studied the dissipation half-lives of 10 micropollutants in sediment-water incubations (based on the OECD 308 standard) with sediment from two European rivers sampled upstream and downstream of wastewater treatment plant (WWTP) discharge. Dissipation half-lives (DT50s) were highly variable between the tested compounds, ranging from 1.5 to 772 days. Sediment from one river sampled downstream from the WWTP showed the fastest dissipation of all micropollutants after sediment RNA normalization. By characterizing sediment bacteria using 16S rRNA sequences, bacterial community composition of a sediment was associated with its capacity for dissipating micropollutants. Bacterial amplicon sequence variants of the genera Ralstonia, Pseudomonas, Hyphomicrobium, and Novosphingobium, which are known degraders of contaminants, were significantly more abundant in the sediment incubations where fast dissipation was observed. Our study illuminates the limitations of the OECD 308 standard to account for variation of dissipation rates of micropollutants due to differences in bacterial community composition. This limitation is problematic particularly for those compounds with DT50s close to regulatory persistence criteria. Thus, it is essential to consider bacterial community composition as a source of variability in regulatory biodegradation and persistence assessments.

Temporal and spatial variable in-stream attenuation of selected pharmaceuticals

Organic micropollutants enter rivers mainly with discharges of wastewater treatment plants (WWTP) and pose a risk to aquatic ecosystems and water quality. A considerable knowledge gap exists for disentangling overlapping processes and driving conditions that control the fate of these pollutants. Thus, the aim of this study was to identify the driving parameters for attenuation of selected pharmaceuticals (carbamazepine, diclofenac, tramadol and venlafaxine) under field conditions. The presented study was performed at a small river (Ammer River, mean discharge 0.87 m³ s^-1) which is hydrologically complex due to karstification, numerous artificial discharges, and engineered modifications of the channel. We applied a Lagrangian sampling scheme at two sequential river reaches. In general, for the investigated compounds and over the length of the tested reaches, the absolute net attenuation representative for 24 h was low (≤ 23% net attenuation), yet calculated half-lives were comparable to literature. Photodegradation is specifically relevant for the first river reach characterized by a higher net attenuation of the photosensitive compound diclofenac (14.5% ±11.3%) compared to the second section (9.8% ±13.7%). This is likely due to a spatial difference in canopy shading, which is supported by significant correlations (R² ≥ 0.8) of the temporally changing 'temperature' and 'solar radiation' with time-specific degradation rate constants of photosensitive compounds for consecutive hourly water parcels. In general, the presented spatially and temporally resolved approach is a suitable tool to determine the attenuation of organic micropollutants and to narrow down the interpretation of net attenuation to a few reasonable processes.

Enantiomeric Fractionation during Biotransformation of Chiral Pharmaceuticals in Recirculating Water-Sediment Test Flumes

Many organic contaminants entering the aquatic environment feature stereogenic structural elements that give rise to enantiomerism. While abiotic processes usually act identical on enantiomers, biotic processes, such as biodegradation often result in enantiomeric fractionation (EFr), i.e., the change of the relative abundance of enantiomers. Therefore, EFr offers the opportunity to differentiate biodegradation in complex environmental systems from abiotic processes. In this study, an achiral-chiral two-dimensional liquid chromatographic method for the enantioseparation of selected pharmaceuticals was developed. This method was then applied to determine the enantiomeric compositions of eight chiral pharmaceuticals in 20 water-sediment test flumes and test EFr as an indicator of biodegradation. While all eight substances were attenuated by at least 60%, five (atenolol, metoprolol, celiprolol, propranolol, and flecainide) displayed EFr. No EFr was observed for citalopram, fluoxetine, and venlafaxine despite almost complete attenuation (80 to 100%). Celiprolol, a barely studied β-blocker, revealed the most distinct EFr among all investigated substances; however, EFr varied considerably with biodiversity. Celiprolol-H2 was identified as a biological transformation product possibly formed by reduction of the celiprolol keto group through a highly regio- and enantioselective carbonyl reductase. While celiprolol-H2 was observed across all flumes, as expected, its formation was faster in flumes with high bacterial diversity where also EFr was highest. Overall, EFr and transformation product formation together served as good indicators of biological processes; however, the strong dependence of EFr on biodiversity limits its usefulness in complex environmental systems.

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.

Comparing non-targeted chemical persistence assessed using an unspiked OECD 309 test to field measurements

Previous research has shown that unspiked OECD 309 tests can be used to quantify chemical biodegradation in surface waters, relying on chemical residues already present in the water. Here we test the hypothesis that unspiked OECD 309 tests can quantitatively predict chemical persistence in the environment by comparing chemical half-lives assessed in the laboratory against those measured in the field. The study object was a Swedish lake heavily impacted by treated municipal wastewater. Half-lives in the field were measured by mass balance over 12 weeks. In parallel, half-lives in the lab were determined with an unspiked OECD 309 test run for 60 days. Chemical analysis was conducted using a non-target screening approach. The field study yielded a half-life <100 days for 38 chemicals for which the dominant source was wastewater; 32 of these were also detected in the lab test, whereby 18 had half-lives with a well-constrained uncertainty that did not intersect infinity. For 14 of the 18 chemicals, the field and lab half-lives agreed within a factor 3. In summary, the lab test predicted chemical attenuation in the field well. Limitations of the approach include the need for measurable chemical concentrations in the water body and failure to account for some attenuation mechanisms like phototransformation.

Analysis of mobile chemicals in the aquatic environment-current capabilities, limitations and future perspectives

Persistent and mobile water contaminants are rapidly developing into a focal point of environmental chemistry and chemical regulation. Their defining parameter that sets them apart from the majority of regularly monitored and regulated contaminants is their mobility in the aquatic environment, which is intrinsically tied to a high polarity. This high polarity, however, may have severe implications in the analytical process and thus the most polar of these mobile contaminants may not be covered by widely utilized trace-analytical methods, and thus, alternatives are required. In this review, we infer the physical and chemical properties of mobile water contaminants from a set of almost 1800 prioritized REACH chemicals and discuss the implications these substance properties may have on four integral steps of the analytical process: sampling and sample storage, sample pre-treatment, separation and detection. We discuss alternatives to widely utilized trace-analytical methods, examine their application range and limitations, highlight potential analytical techniques on the horizon and emphasize research areas we believe still offer the most room for further improvement. While we have a comprehensive set of analytical methods to cover a large portion of the known mobile chemicals, these methods are still only infrequently utilized. Graphical abstract.

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.

Impact of Bed Form Celerity on Oxygen Dynamics in the Hyporheic Zone

Oxygen distribution and uptake in the hyporheic zone regulate various redox-sensitive reactions and influence habitat conditions. Despite the fact that fine-grain sediments in streams and rivers are commonly in motion, most studies on biogeochemistry have focused on stagnant sediments. In order to evaluate the effect of bed form celerity on oxygen dynamics and uptake in sandy beds, we conducted experiments in a recirculating indoor flume. Oxygen distribution in the bed was measured under various celerities using 2D planar optodes. Bed morphodynamics were measured by a surface elevation sensor and time-lapse photography. Oxygenated zones in stationary beds had a conchoidal shape due to influx through the stoss side of the bed form, and upwelling anoxic water at the lee side. Increasing bed celerity resulted in the gradual disappearance of the upwelling anoxic zone and flattening of the interface between the oxic (moving fraction of the bed) and the anoxic zone (stationary fraction of the bed), as well as in a reduction of the volumetric oxygen uptake rates due shortened residence times in the hyporheic zone. These results suggest that including processes related to bed form migration are important for understanding the biogeochemistry of hyporheic zones.

Mapping Micro-Pollutants and Their Impacts on the Size Structure of Streambed Communities

Recently there has been increasing concern over the vast array of emerging organic contaminants (EOCs) detected in streams and rivers worldwide. Understanding of the ecological implications of these compounds is limited to local scale case studies, partly as a result of technical limitations and a lack of integrative analyses. Here, we apply state-of-the-art instrumentation to analyze a complex suite of EOCs in the streambed of 30 UK streams and their effect on streambed communities. We apply the abundance–body mass (N–M) relationship approach as an integrative metric of the deviation of natural communities from reference status as a result of EOC pollution. Our analysis includes information regarding the N and M for individual prokaryotes, unicellular flagellates and ciliates, meiofauna, and macroinvertebrates. We detect a strong significant dependence of the N–M relationship coefficients with the presence of EOCs in the system, to the point of shielding the effect of other important environmental factors such as temperature, pH, and productivity. However, contrary to other stressors, EOC pollution showed a positive effect on the N–M coefficient in our work. This phenomenon can be largely explained by the increase in large-size tolerant taxa under polluted conditions. We discuss the potential implications of these results in relation to bioaccumulation and biomagnification processes. Our findings shed light on the impact of EOCs on the organization and ecology of the whole streambed community for the first time.

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.

An immission perspective of emerging micropollutant pressure in Luxembourgish surface waters: A simple evaluation scheme for wastewater impact assessment

While wastewater treatment plants have been identified as the most prominent source of emerging micropollutants in surface waters, prediction of their ambient concentrations remains a challenge. This is due to the variability of loads entering individual treatment plants and of the elimination capacity by the latter as well as potential attenuation in the river network. Although geospatially detailed models exist, they suffer from the same data input uncertainties. Here, we investigated the concentration profiles of 20 emerging pollutants in different river stretches in Luxembourg with variable sanitary pressures. Using carbamazepine as a recalcitrant wastewater indicator, the correlation of the compounds to the latter revealed source and fate variability as well as specific emitters. Relating carbamazepine to sanitary pressure, expressed as the sum of population equivalents in a catchment divided by its surface [PE ha^-1] allowed predicting the impact of emerging pollutants on the entire river network. The limited variability of the pollutant profiles allowed for prioritization of impacted stretches depending on the different sanitary pressures at risk quotient exceedance. The main drivers of impact were triclosan, diclofenac, clarithromycine and diuron.

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.