Model-based evaluation of reduction strategies for micropollutants from wastewater treatment plants in complex river networks

Predicting environmental risks of pharmaceutical residues by wastewater surveillance: An analysis based on pharmaceutical sales and their excretion data

Pharmaceutical residues are increasingly becoming a significant source of environmental water pollution and ecological risk. This study, leveraging official national pharmaceutical sales statistics, predicts the environmental concentrations of 33 typical pharmaceuticals in the Tianjin area. The results show that 52 % of the drugs have a PEC/MEC (Predicted Environmental Concentration/Measured Environmental Concentration) ratio within the acceptable range of 0.5-2, including atenolol (1.21), carbamazepine (1.22), and sulfamethoxazole (0.91). Among the selected drugs, tetracycline, ciprofloxacin, and acetaminophen had the highest predicted concentrations. The EPI (Estimation Programs Interface) biodegradation model, a tool from the US Environmental Protection Agency, is used to predict the removal efficiency of compounds in wastewater treatment plants. The results indicate that the EPI predictions are acceptable for macrolide antibiotics and β-blockers, with removal rates of roxithromycin, spiramycin, acetaminophen, and carbamazepine being 14.1 %, 61.2 %, 75.1 %, and 44.5 %, respectively. However, the model proved to be less effective for fluoroquinolone antibiotics. The ECOSAR (Ecological Structure-Activity Relationships) model was used to supplement the assessment of the potential impacts of drugs on aquatic ecosystems, further refining the analysis of pharmaceutical environmental risks. By combining the concentration and detection frequency of pharmaceutical wastewater, this study identified 9 drugs with significant toxicological risks and marked another 24 drugs as substances of potential concern. Additionally, this study provides data support for addressing pharmaceutical residues of priority concern in subsequent research.

A graph-based modeling framework for tracing hydrological pollutant transport in surface waters

Anthropogenic pollution of hydrological systems affects diverse communities and ecosystems around the world. Data analytics and modeling tools play a key role in fighting this challenge, as they can help identify key sources as well as trace transport and quantify impact within complex hydrological systems. Several tools exist for simulating and tracing pollutant transport throughout surface waters using detailed physical models; these tools are powerful, but can be computationally intensive, require significant amounts of data to be developed, and require expert knowledge for their use (ultimately limiting application scope). In this work, we present a graph modeling framework - which we call HydroGraphs - for understanding pollutant transport and fate across waterbodies, rivers, and watersheds. This framework uses a simplified representation of hydrological systems that can be constructed based purely on open-source data (National Hydrography Dataset and Watershed Boundary Dataset). The graph representation provides a flexible intuitive approach for capturing connectivity and for identifying upstream pollutant sources and for tracing downstream impacts within small and large hydrological systems. Moreover, the graph representation can facilitate the use of advanced algorithms and tools of graph theory, topology, optimization, and machine learning to aid data analytics and decision-making. We demonstrate the capabilities of our framework by using case studies in the State of Wisconsin; here, we aim to identify upstream nutrient pollutant sources that arise from agricultural practices and trace downstream impacts to waterbodies, rivers, and streams. Our tool ultimately seeks to help stakeholders design effective pollution prevention/mitigation practices and evaluate how surface waters respond to such practices.

Campus node-based wastewater surveillance enables COVID-19 case localization and confirms lower SARS-CoV-2 burden relative to the surrounding community

Wastewater-based surveillance (WBS) has been established as a powerful tool that can guide health policy at multiple levels of government. However, this approach has not been well assessed at more granular scales, including large work sites such as University campuses. Between August 2021 and April 2022, we explored the occurrence of SARS-CoV-2 RNA in wastewater using qPCR assays from multiple complimentary sewer catchments and residential buildings spanning the University of Calgary's campus and how this compared to levels from the municipal wastewater treatment plant servicing the campus. Real-time contact tracing data was used to evaluate an association between wastewater SARS-CoV-2 burden and clinically confirmed cases and to assess the potential of WBS as a tool for disease monitoring across worksites. Concentrations of wastewater SARS-CoV-2 N1 and N2 RNA varied significantly across six sampling sites - regardless of several normalization strategies - with certain catchments consistently demonstrating values 1-2 orders higher than the others. Relative to clinical cases identified in specific sewersheds, WBS provided one-week leading indicator. Additionally, our comprehensive monitoring strategy enabled an estimation of the total burden of SARS-CoV-2 for the campus per capita, which was significantly lower than the surrounding community (p≤0.001). Allele-specific qPCR assays confirmed that variants across campus were representative of the community at large, and at no time did emerging variants first debut on campus. This study demonstrates how WBS can be efficiently applied to locate hotspots of disease activity at a very granular scale, and predict disease burden across large, complex worksites.

Spatiotemporal distribution of veterinary and human drugs and its predictability in Japanese catchments

Little is known about the predictability of mass flows of veterinary drugs in Asian catchments, where effluent from livestock farms is a major source. We therefore conducted this study to understand the applicability and limitations of a population-based emission model, which assumed usage of veterinary and human drugs to be evenly distributed over the national livestock or human population throughout the year, and sources to be effluent discharges at livestock farms, households, and sewage treatment plants in Japanese catchments. We monitored five veterinary drugs (lincomycin, sulfamonomethoxine, tiamulin, tylosin, and tilmicosin), two human and livestock drugs (sulfamethoxazole and trimethoprim), two human drugs (carbamazepine and clarithromycin), and a metabolite (sulfapyridine) of a human drug once a month over 2 years in eight Japanese rivers which have active livestock farming in their catchments. Mass flows of carbamazepine and sulfapyridine were stable, while those of veterinary drugs fluctuated widely, especially sulfamonomethoxine and tilmicosin, whose 25 %-100 % ranges averaged 1.5 and 1.2 log units, respectively, attributable mainly to their usage patterns. The model accurately predicted mean mass flows of carbamazepine in the rivers with errors of <±0.3 log unit. Although it slightly to moderately overestimated those of the other four human-related compounds, the incorporation of an empirical correction factor, determined to minimize mean absolute error (MAE) among the rivers, substantially lowered their MAEs to <0.23 log units. However, the MAEs of the five veterinary drugs were as high as 0.42 (sulfamonomethoxine) to 0.60 (tiamulin) log units even with the coefficient, likely due mainly to the spatial distribution of their usage per capita. So as not to overlook spatiotemporal elevation of risks of veterinary drugs, a stochastic method should be applied in their management. This is the first study to assess the use of spatiotemporal homogeneity in usage per capita of veterinary drugs in Asian catchments.

Fungi-enabled pore channel regulation and defect engineering of a novel micro-reactor for treating complex effluents

Defect engineering has become a significant research area in recent years; however, little has been reported on the biological method for modulating the intrinsic carbon defects of the biochar framework. Herein, a fungi-enabled method for the fabrication of porous carbon/Fe₃O₄/Ag (PC/Fe₃O₄/Ag) composites was developed, and the mechanism underlying the hierarchical structure is elucidated for the first time. By regulating the cultivation process of fungi on water hyacinth biomass, a well-developed interconnected structure and carbon defects acting as potential catalytic active sites were formed. This new material with antibacterial, adsorption and photodegradation properties could be an excellent choice for treating the mixed dyestuff effluents with oils and bacteria, also guiding pore channel regulation and defect engineering in materials science. Numerical simulations were carried out to demonstrate the remarkable catalytic activity.

Occurrence of pharmaceuticals in the wastewater of a Greek hospital: Combining consumption data collection and LC-QTOF-MS analysis

In this article we applied drug consumption approach and chemical analysis in parallel to investigate the concentrations of a large number of pharmaceuticals in different streams of a General Hospital. Drugs consumption data was collected during two periods (Period 1, 2) and the predicted environmental concentrations (PECs) were estimated for the wastewater of a building housing specific medical services (Point A) and for the entire hospital (Point B). Hospital wastewater samples (HWW) samples were also collected from these points and periods and the measured environmental concentrations (MEC) were determined using UHPLC-ESI-QTOF-MS/MS. According to consumption data, the highest number of drugs was consumed in the departments of Hematology, Intensive Care Unit, Cardiology, Internal Medicine, and Oncology, while the number of active substances used in the hospital was 413 (Period 1) and 362 (Period 2). For most substances, much higher PEC and MEC values were found at the HWW of Point A indicating that on-site treatment of this stream could be examined in the future. The application of wide-scope target analysis allowed the quantification of 122 compounds, while 21 additional substances were identified using suspect screening. The highest mean concentrations in Period 1 were found for acetaminophen (1100 μg/L) and rifaximin (723 μg/L), while in Period 2 for iopromide (458 μg/L) and acyclovir (408 μg/L). Among the detected compounds, 19 metabolites were determined. Atenolol acid, 1-hydroxy-midazolam and clopidogrel carboxylic acid were quantified at concentrations much higher than parent compounds indicating the importance of metabolites' monitoring in HWW. Calculation of PEC/MEC ratio for 36 pharmaceuticals showed sufficient correlation of these values for 19 % to 33 % of the substances depending on the examined period and sampling point. The parallel collection of drugs consumption data and chemical analysis give a thorough picture of the substances present in HWW and their main sources, facilitating decision-making for their better management.

An advection-reaction model of emerging contaminants in ground and surface waterbodies with known topology: a data-based approach

A simple data-based advection-reaction (reactive transport) model applicable to both rivers and aquifers monitoring networks is proposed. It is built on (a) available monitoring data, and (b) graph-theoretical concepts, specifically making use of the Laplacian matrix to capture the network topology and the advection process. The method yields useful information regarding the dynamic spatial behavior of the variables monitored, expressed in terms of quantitative parameters like characteristic length, entropy, first-order decay constants, synchronization between sites, and the external inputs/outputs to the system. The model was tested in an unconfined shallow aquifer located in the lower Besòs River (Spain), in which 37 pharmaceutical compounds were monitored at 7 sites, alongside two campaigns (February and May 2021). Characteristic lengths were, on average, of the same order (24.5 m) as the mean distance between consecutive monitoring sites (33.6 m), thus reflecting an adequate monitoring network design. From an estimated mean advection velocity (0.24 m·h^-1), first-order decay constants were calculated for each compound and campaign, with mean values of 0.025 h^-1 (February) and 0.005 h^-1 (May). Whereas entropy was generally slightly larger values in February than in May (mean values of 1.02 and 0.9 entropy units respectively), synchronization showed the opposite trend (mean values of 62.4% and 68.8% respectively). The input/output profiles were generally site-dependent, regardless of the compound, and campaign considered. • A new advection-reaction modeling approach directly based on experimental data obtained from monitoring campaigns together with the network topology is proposed. • The method yields new quantitative information regarding the dynamic behavior of the variables monitored, useful for both research and management purposes.

European scale assessment of the potential of ozonation and activated carbon treatment to reduce micropollutant emissions with wastewater

Micropollutants (MPs) in wastewater pose a growing concern for their potential adverse effects on the receiving aquatic environment, and some countries have started requiring that wastewater treatment plants remove them to a certain extent. Broad spectrum advanced treatment processes, such as ozonation, activated carbon or their combination, are expected to yield a significant reduction in the toxicity of effluents. Here we quantify the reduction of effluent toxicity potentially achieved by implementing these advanced treatment solutions in a selection of European wastewater treatment plants. To this end, we refer to a list of "total pollution proxy substances" (TPPS) composed of 1337 chemicals commonly found in wastewater effluents according to a compilation of datasets of measured concentrations. We consider these substances as an approximation of the "chemical universe" impinging on the European wastewater system. We evaluate the fate of the TPPS in conventional and advanced treatment plants using a compilation of experimental physicochemical properties that describe their sorption, volatilization and biodegradation during activated sludge treatment, as well as known removal efficiency in ozonation and activated carbon treatment, while filling the gaps through in silico prediction models. We estimate that the discharge of micropollutants with wastewater effluents in the European Union has a cumulative MP toxicity to the environment equal to the discharge of untreated wastewater of ca. 160 million population equivalents (PE), i.e. about 30 % of the generated wastewater in the EU. If all plants above a capacity of 100,000 PE were equipped with advanced treatment, we show that this load would be reduced to about 95 million PE. In addition, implementing advanced treatment in wastewater plants above 10,000 PE discharging to water bodies with an average dilution ratio smaller than 10 would yield a widespread improvement in terms of exposure of freshwater ecosystems to micropollutants, almost halving the part of the stream network exposed to the highest toxic risks. Our analysis provides background for a cost-effectiveness appraisal of advanced treatment "at the end of the pipe", which could lead to optimized interventions. This should not be regarded as a stand-alone solution, but as a complement to policies for the control of emissions at the source for the most problematic MPs.

Biomarkers selection for population normalization in SARS-CoV-2 wastewater-based epidemiology

Wastewater-based epidemiology (WBE) has been one of the most cost-effective approaches to track the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) levels in the communities since the coronavirus disease 2019 (COVID-19) outbreak in 2020. Normalizing SARS-CoV-2 concentrations by the population biomarkers in wastewater is critical for interpreting the viral loads, comparing the epidemiological trends among the sewersheds, and identifying the vulnerable communities. In this study, five population biomarkers, pepper mild mottle virus (PMMoV), creatinine (CRE), 5-hydroxyindoleacetic acid (5-HIAA), caffeine (CAF) and its metabolite paraxanthine (PARA) were investigated and validated for their utility in normalizing the SARS-CoV-2 loads through two normalizing approaches using the data from 64 wastewater treatment plants (WWTPs) in Missouri. Their utility in assessing the real-time population contributing to the wastewater was also evaluated. The best performing candidate was further tested for its capacity for improving correlation between normalized SARS-CoV-2 loads and the clinical cases reported in the City of Columbia, Missouri, a university town with a constantly fluctuating population. Our results showed that, except CRE, the direct and indirect normalization approaches using biomarkers allow accounting for the changes in wastewater dilution and differences in relative human waste input over time regardless flow volume and population of the given WWTP. Among selected biomarkers, PARA is the most reliable population biomarker in determining the SARS-CoV-2 load per capita due to its high accuracy, low variability, and high temporal consistency to reflect the change in population dynamics and dilution in wastewater. It also demonstrated its excellent utility for real-time assessment of the population contributing to the wastewater. In addition, the viral loads normalized by the PARA-estimated population significantly improved the correlation (rho=0.5878, p < 0.05) between SARS-CoV-2 load per capita and case numbers per capita. This chemical biomarker complements the current normalization scheme recommended by CDC and helps us understand the size, distribution, and dynamics of local populations for forecasting the prevalence of SARS-CoV2 within each sewershed.

Biomarkers Selection for Population Normalization in SARS-CoV-2 Wastewater-based Epidemiology

Wastewater-based epidemiology (WBE) has been one of the most cost-effective approaches to track the SARS-CoV-2 levels in the communities since the COVID-19 outbreak in 2020. Normalizing SARS-CoV-2 concentrations by the population biomarkers in wastewater can be critical for interpreting the viral loads, comparing the epidemiological trends among the sewersheds, and identifying the vulnerable communities. In this study, five population biomarkers, pepper mild mottle virus (pMMoV), creatinine (CRE), 5-hydroxyindoleacetic acid (5-HIAA), caffeine (CAF) and its metabolite paraxanthine (PARA) were investigated for their utility in normalizing the SARS-CoV-2 loads through developed direct and indirect approaches. Their utility in assessing the real-time population contributing to the wastewater was also evaluated. The best performed candidate was further tested for its capacity for improving correlation between normalized SARS-CoV-2 loads and the clinical cases reported in the City of Columbia, Missouri, a university town with a constantly fluctuated population. Our results showed that, except CRE, the direct and indirect normalization approaches using biomarkers allow accounting for the changes in wastewater dilution and differences in relative human waste input over time regardless flow volume and population at any given WWTP. Among selected biomarkers, PARA is the most reliable population biomarker in determining the SARS-CoV-2 load per capita due to its high accuracy, low variability, and high temporal consistency to reflect the change in population dynamics and dilution in wastewater. It also demonstrated its excellent utility for real-time assessment of the population contributing to the wastewater. In addition, the viral loads normalized by the PARA-estimated population significantly improved the correlation ( rho =0.5878, p <0.05) between SARS-CoV-2 load per capita and case numbers per capita. This chemical biomarker offers an excellent alternative to the currently CDC-recommended pMMoV genetic biomarker to help us understand the size, distribution, and dynamics of local populations for forecasting the prevalence of SARS-CoV2 within each sewershed.

HIGHLIGHT bullet points

The paraxanthine (PARA), the metabolite of the caffeine, is a more reliable population biomarker in SARS-CoV-2 wastewater-based epidemiology studies than the currently recommended pMMoV genetic marker.SARS-CoV-2 load per capita could be directly normalized using the regression functions derived from correlation between paraxanthine and population without flowrate and population data.Normalizing SARS-CoV-2 levels with the chemical marker PARA significantly improved the correlation between viral loads per capita and case numbers per capita.The chemical marker PARA demonstrated its excellent utility for real-time assessment of the population contributing to the wastewater.

Functionalized Multi-Walled Carbon Nanotube–Based Aptasensors for Diclofenac Detection

Driven by the increasing concern about the risk of diclofenac (DCF) residues as water pollutants in the aqueous environment and the growing need for its trace determination, a simple but sensitive electrochemical aptasensor for the trace detection of DCF was developed. To construct the aptasensor, the amine-terminated DCF aptamer was covalently immobilized on the surface of the carboxylic acid–functionalized multi-walled carbon nanotube (f-MWCNT)–modified glassy carbon electrode (GCE) through EDC/NHS chemistry. The f-MWCNTs provide a reliable matrix for aptamer immobilization with high grafting density, while the aptamer serves as a biorecognition probe for DCF. The obtained aptasensor was incubated with DCF solutions at different concentrations and was then investigated by electrochemical impedance spectroscopy (EIS). It displays two linear ranges of concentration for DCF detection, from 250 fM to 1pM and from 1 pM to 500 nM with an extremely low detection limit of 162 fM. Also, the developed biosensor shows great reproducibility, acceptable stability, and reliable selectivity. Therefore, it offers a simple but effective aptasensor construction strategy for trace detection of DCF and is anticipated to show great potential for environmental applications.

Effects of diclofenac on the swimming behavior and antioxidant enzyme activities of the freshwater interstitial crustacean Bryocamptus pygmaeus (Crustacea, Harpacticoida)

Diclofenac (DCF) is one of the most widespread pharmaceutical compounds found in freshwaters as a pseudo-persistent pollutant due to its continuous release from point and diffuse sources, being its removal in Wastewater Treatment Plants incomplete. Moreover, DCF is particularly persistent in interstitial habitats and potentially toxic for the species that spend their whole life cycle among the same sediment grains. This study is aimed at offering a first contribution to the assessment of DCF effects on freshwater invertebrate species living in the interstitial habitats of springs, rivers, lakes and groundwaters. The Crustacea Copepoda are one of the main components of the freshwater interstitial communities, with the primacy taken by the worm-like and small-sized harpacticoids. A sub-lethal concentration of 50 μg L^-1 DCF significantly affected six out of the eight behavior parameters of the burrower/interstitial crustacean harpacticoid Bryocamptus pygmaeus recorded by video tracking analysis. DCF exposure reduced swimming speed, swimming activity, exploration ability and thigmotaxis, and increased swimming path tortuosity. The biochemical approach revealed a reduced level of the mitochondrial superoxide dismutase 2 in individuals exposed to DCF. It could be explained by a decline in mitochondrial performance or by a reduced number of functional mitochondria. Since mitochondrial dysfunction may determine ATP reduction, it comes that less energy is produced for maintaining the cell functions of the DCF-exposed individuals. In addition, the increasing energy demand for the detoxification process further contributes to decrease the total energetic budget allocated for other physiological activities. These observations can explain the changes we have observed in the swimming behavior of the copepod B. pygmaeus.

Shifted mineral ions transport in the mollusk Unio pictorum exposed to environmental concentrations of diclofenac

Previous studies showed that diclofenac (DCF), when released in the environment, can be toxic to aquatic animals (fish and mollusks), affecting gills, which are the main organ of ionic regulation. This study focuses on detecting the effects of relevant environmental concentrations of DCF (0.1-1 μg L^-1) on the transport of main mineral cations, i.e. sodium (Na), potassium (K), calcium (Ca), and magnesium (Mg), by widely distributed freshwater bivalve mollusks Unio pictorum. After 96-h exposure to river aerated water at 25 °C with DCF concentrations of 0 (control), 0.1 (treatment I), and 1 μg L^-1 (treatment II), the mollusks were transferred to deionized water, and daily (for 7 days) concentrations of these cations in the medium have been measured. Animals exposed to 1 μg L^-1 DCF maintained the ionic balance between the organism and the diluted medium at a significantly higher level of Na, K, and Mg ions in water compared to the control and animals exposed to 0.1 μg L^-1 DCF. At 0.1 μg L^-1 DCF, the greater loss concerning the control (p < 0.05) was found only for Na ion. There were no differences in the dynamics of Ca ions between control and both treatments. This study showed that detectable environmental concentrations of DCF in natural waters can influence the transport of main cations required by freshwater animals to maintain their ionic balance, and the observed effect (elevated ion loss) is ion-specific and also dose-dependent.

Geo-referenced simulation of pharmaceuticals in whole watersheds: application of GREAT-ER 4.1 in Germany

The geo-referenced regional exposure assessment tool for European rivers (GREAT-ER) is designed to support river basin management or the implementation process within the EU Water Framework Directive by predicting spatially resolved exposure concentrations in whole watersheds. The usefulness of the complimentary application of targeted monitoring and GREAT-ER simulations is demonstrated with case studies for three pharmaceuticals in selected German watersheds. Comparison with monitoring data corroborates the capability of the probabilistic model approach to predict the expected range of spatial surface water concentrations. Explicit consideration of local pharmaceutical emissions from hospitals or private doctor's offices (e.g., for X-ray contrast agents) can improve predictions on the local scale without compromising regional exposure assessment. Pharmaceuticals exhibiting low concentrations hardly detectable with established analytical methods (e.g., EE2) can be evaluated with model simulations. Management scenarios allow for a priori assessment of risk reduction measures. In combination with targeted monitoring approaches, the GREAT-ER model can serve as valuable support tool for exposure and risk assessment of pharmaceuticals in whole watersheds.

Can source control of pharmaceuticals decrease the investment needs in urban wastewater infrastructure?

The source control of pharmaceuticals involves influencing the everyday consumption volume and compound choice. This paper evaluates how source control contributes to protecting the environmental health and decreasing the investment needs in urban wastewater infrastructure. Different levels of reduction in diclofenac consumption (as recommended by the European Medicines Agency) compensated by equivalent increases in naproxen consumption (a less environmentally harmful compound) are evaluated. The different loads of compounds are fed into a microcontaminant fate and transport model of the Llobregat river basin (Spain) to assess the investment needs in tertiary treatment to reach diclofenac and naproxen concentrations below environmental quality standards. The results show that, despite the implementation of source control measures, tertiary treatment upgrades are still required in every scenario evaluated. Even though source control of pharmaceuticals decreases the investment needs in urban wastewater infrastructure, apparent concentrations reductions (i.e. statistically significant differences relative to the reference situation) are only observed in drastic substitutions of diclofenac by naproxen (a reduction in the total diclofenac consumption by 73% and a corresponding increase in naproxen consumption). The results also show that Spain is on good track with regards to the substitution of diclofenac by naproxen (among the top 5 in Europe), and this paper shows how positive this substitution can be for the environment.

Tool for selecting indicator substances to evaluate the impact of wastewater treatment plants on receiving water bodies

The elimination of organic micropollutants (OMPs) from wastewater could in future become mandatory for operators of wastewater treatment plants (WWTPs). Indicator substances are a great help and a cost-efficient way in monitoring the pollution of water bodies with OMPs caused by the discharge of WWTPs. However, with the still increasing number of OMPs in our environment, the selection of suitable indicator substances presents a challenge. A concept was developed to help identify representative indicator substances. The derived indicator substances are not only used to assess water pollution, but can also be used to calculate elimination efficiencies of WWTPs. In the present investigations, the indicator substances were used to evaluate the reduction of OMPs in the water body on the basis of the expansion of a WWTP with an ozonation plant. The transferability of the tool was verified with a second WWTP. Furthermore, the impact of the number of measurements was analysed via statistical combinatorics. With the tool, 36 substances were classified, leading to the identification of 9 suggested indicator substances. Among them ibuprofen and diclofenac attracted attention due to their ecotoxicological relevance. Detailed data analyses were carried out using principal component analysis (PCA) and loads.

Management actions to mitigate the occurrence of pharmaceuticals in river networks in a global change context

Human consumption of pharmaceuticals leads to high concentrations of pharmaceuticals in wastewater, which is usually not or insufficiently collected and treated before release into freshwater ecosystems. There, pharmaceuticals may pose a threat to aquatic biota. Unfortunately, occurrence data of pharmaceuticals in freshwaters at the global scale is scarce and unevenly distributed, thus preventing the identification of hotspots, the prediction of the impact of Global Change (particularly streamflow and population changes) on their occurrence, and the design of appropriate mitigation actions. Here, we use diclofenac (DCL) as a typical pharmaceutical contaminant, and a global model of DCL chemical fate based on wastewater sanitation, population density and hydrology to estimate current concentrations in the river network, the impact of future changes in runoff and population, and potential mitigation actions in line with the Sustainable Development Goals. Our model is calibrated against measurements available in the literature. We estimate that 2.74 ± 0.63% of global river network length has DCL concentrations exceeding the proposed EU Watch list limit (100 ng L^-1). Furthermore, many rivers downstream from highly populated areas show values beyond 1000 ng L^-1, particularly those associated to megacities in Asia lacking sufficient wastewater treatment. This situation will worsen with Global Change, as streamflow changes and human population growth will increase the proportion of the river network above 100 ng L^-1 up to 3.10 ± 0.72%. Given this background, we assessed feasible source and end-of-pipe mitigation actions, including per capita consumption reduction through eco-directed sustainable prescribing (EDSP), the implementation of the United Nations Sustainable Development Goal (SDG) 6 of halving the proportion of population without access to safely managed sanitation services, and improvement of wastewater treatment plants up to the Swiss standards. Among the considered end-of-pipe mitigation actions, implementation of SDG 6 was the most effective, reducing the proportion of the river network above 100 ng L^-1 down to 2.95 ± 0.68%. However, EDSP brought this proportion down to 2.80 ± 0.64%. Overall, our findings indicate that the sole implementation of technological improvements will be insufficient to prevent the expected increase in pharmaceuticals concentration, and that technological solution need to be combined with source mitigation actions.

Spatial Differences among Micropollutants in Sewer Overflows: A Multisite Analysis Using Passive Samplers

Untreated sewer overflows can contaminate receiving waters with micropollutants. Although concentrations of discharged micropollutants can be ecotoxicologically relevant, only limited data is available to assess occurrence and spatial differences among sewer overflow catchments. Therefore, we present an innovative type of data obtained with passive samplers at 20 combined sewer overflow sites (2-7 events per site; 95 events in total). The data sheds light on concentration ranges for 13 representative polar organic micropollutants and shows that micropollutants in both municipal wastewater and stormwater can be relevant sources of contaminants. We identify indicator micropollutants for further studies: benzotriazole (80% interquantile of time-weighted average concentration: 250-4800 ng/L), carbamazepine (33-910 ng/L), diclofenac (78-1000 ng/L), carbendazim (21-900 ng/L), diazinon (2.1-53 ng/L), diuron (22-1100 ng/L), mecoprop (98-5300 ng/L), metolachlor (6-230 ng/L), and terbutryn (29-810 ng/L). These concentration estimates are assumed to be on the safe side for comparison with environmental quality standards (EQS). A majority of sewer overflow sites (13 of 20) show discharge concentrations above acute EQS for at least one micropollutant and thus would have to rely on dilution by receiving waters to not exceed any EQS. The intersite variability among sewer overflows exceed the within-site variability. Hence, future monitoring studies should cover more sewer overflow sites. No correlation could be found with event durations, specific storage volume or land use data, thus showing the complexity of micropollutant occurrence and indicating that other factors led to the observed high spatial variability. In conclusion, our results clearly show the potential relevance of micropollutants in sewer overflows and the need to assess site-specific measures.

Sale-based estimation of pharmaceutical concentrations and associated environmental risk in the Japanese wastewater system

Information on sales and emission of selected pharmaceuticals were used to predict their concentrations in Japanese wastewater influent through a >300 of pharmaceuticals data sink. A combined wastewater-based epidemiology and environmental risk analysis follow was established. By comparing predicted environmental concentrations (PECs) of pharmaceuticals in wastewater influent against measured environmental concentrations (MECs) reported in previous studies, it was found that the model gave accurate results for 17 pharmaceuticals (0.5 < PEC/MEC < 2), and acceptable results for 32 out of 40 pharmaceuticals (0.1 < PEC/MEC < 10). Although the majority of pharmaceuticals considered in the model were antibiotics and analgesics, pranlukast, a receptor antagonist, was predicted to have the highest concentration in wastewater influent. With regard to the composition of wastewater effluent, the Estimation Program Interface (EPI) suite was used to predict pharmaceutical removal through activated sludge treatment. Although the performance of the EPI suite was variable in terms of accurate prediction of the removal of different pharmaceuticals, it could be an efficient tool in practice for predicting removal under extreme scenarios. By using the EPI suite with input data on PEC in the wastewater influent, the PEC values of pharmaceuticals in wastewater effluent were predicted. The concentrations of 26 pharmaceuticals were relatively high (>1 μg/L), and the PECs of 6 pharmaceuticals were extremely high (>10 μg/L) in wastewater effluent, which could be attributed to their high usage rates by consumers and poor removal rates in wastewater treatment plants (WWTPs). Furthermore, environmental risk assessment (ERA) was carried out by calculating the ratio of predicted no effect concentration (PNEC) to PEC of different pharmaceuticals, and it was found that 9 pharmaceuticals were likely to have high toxicity, and 54 pharmaceuticals were likely to have potential toxicity. It is recommended that this is further investigated in detail. The priority screening and environmental risk assessment results on pharmaceuticals can provide reliable basis for policy-making and environmental management.

Biotransformation Changes Bioaccumulation and Toxicity of Diclofenac in Aquatic Organisms

Biotransformation plays a crucial role in regulating the bioaccumulation potential and toxicity of organic compounds in organisms but is, in general, poorly understood for emerging contaminants. Here, we have used diclofenac as a model compound to study the impact of biotransformation on the bioaccumulation potential and toxicity in two keystone aquatic invertebrates: Gammarus pulex and Hyalella azteca. In both species, diclofenac was transformed into several oxidation products and conjugates, including two novel products, that is, diclofenac taurine conjugate (DCF-M403) and unexpected diclofenac methyl ester (DCF-M310.03). The ratios of biotransformation products to parent compound were 12-17 for DCF-M403 and 0.01-0.7 for DCF-M310.03 after 24 h exposure. Bioconcentration factors (BCFs) of diclofenac were 0.5 and 3.2 L kg_ww^-1 in H. azteca and G. pulex, respectively, whereas BCFs of DCF-M310.03 was 164.5 and 104.7 L kg_ww^-1, respectively, representing a 25- to 110-fold increase. Acute toxicity of DCF-M310.03 was also higher than the parent compound in both species, which correlated well with the increased bioconcentration potential. The LC₅₀ of diclofenac in H. azteca was 216 mg L^-1, while that of metabolite DCF-M310.03 was reduced to only 0.53 mg L^-1, representing a 430-fold increase in acute toxicity compared to diclofenac. DCF-M403 is less toxic than its parent compound toward H. azteca, which may be linked to its slightly lower hydrophobicity. Furthermore, the transformation of diclofenac to its methyl ester derivative was explored in crude invertebrate extracts spiked with an S-adenosylmethionine cofactor, revealing possible catalysis by an S-adenosylmethionine-dependent carboxylic acid methyltransferase. Methylation of diclofenac was further detected in fish hepatocytes and human urine, indicating a broader relevance. Therefore, potentially methylated metabolites of polar contaminants should be considered for a comprehensive risk assessment in the future.

Coupling River Concentration Simulations with a Toxicokinetic Model Effectively Predicts the Internal Concentrations of Wastewater-Derived Micropollutants in Field Gammarids

Although the exposure assessment of wastewater-derived micropollutants via chemical, bioanalytical, and modeling methods in environmental compartments is becoming more frequent, the whole-body burden (i.e., internal concentrations) in nontarget organisms is rarely assessed. An understanding of the internal concentration fluctuation is especially important when exploring the mechanistic linkage between exposure and effects. In this study, we coupled a simple river model with a first-order toxicokinetic (TK) model to predict the concentrations of wastewater-derived micropollutants in freshwater invertebrates (Gammarus spp.). We applied Monte Carlo simulations and conducted laboratory experiments to account for the uncertain input data and the lack of uptake/depuration rate constants required for the TK model. The internal concentrations in field gammarids were predicted well, and the estimates varied only by a factor of 0.1-1.9. Fast equilibrium may also be assumed such that bioconcentration factors (BCFs) are used together with the daily river dilution patterns to predict internal concentrations. While this assumption is suitable for compounds observed in our experiment to reach the steady state within 48 h in gammarids, the model overpredicted the concentrations of substances that reach this condition after longer periods. Nevertheless, this approach provides conservative estimates and simplifies the coupling of models as BCFs are slightly more accessible than the rate constants. However, if one is interested in a more detailed exposure information (e.g., peak concentration and the whole-body burden recovery after a spill), then the nonsteady-state formulation should be employed.

Occurrence, interactive effects and ecological risk of diclofenac in environmental compartments and biota - a review

Diclofenac, a nonsteroidal anti-inflammatory drug has turned into a contaminant of emerging concern; hence, it was included in the previous Watch List of the EU Water Framework Directive. This review paper aims to highlight the metabolism of diclofenac at different trophic levels, its occurrence, ecological risks, and interactive effects in the water cycle and biota over the past two decades. Increased exposure to diclofenac not only raises health concerns for vultures, aquatic organisms, and higher plants but also causes serious threats to mammals. The ubiquitous nature of diclofenac in surface water (river, lake canal, estuary, and sea) is compared with drinking water, groundwater, and wastewater effluent in the environment. This comprehensive survey from previous studies suggests the fate of diclofenac in wastewater treatment plants (WWTPs) and may predict its persistence in the environment. This review offers evidence of fragmentary available data for the water environment, soil, sediment, and biota worldwide and supports the need for further data to address the risks associated with the presence of diclofenac in the environment. Finally, we suggest that the presence of diclofenac and its metabolites in the environment may represent a high risk because of their synergistic interactions with existing contaminants, leading to the development of drug-resistant strains and the formation of newly emerging pollutants.

Linking microbial diversity and population dynamics to the removal efficiency of pharmaceutically active compounds (PhACs) in an anaerobic/anoxic/aerobic (A2O) system

The removal efficiencies (REs) of nineteen pharmaceutically active compounds (PhACs) (six antibiotics-clarithromycin, ofloxacin, sulfadiazine, sulfamethazine, sulfamethoxazole and trimethoprim -, four β-blockers -atenolol, metoprolol, propranolol and sotalol-, two antihypertensives/diuretics -furosemide and hydrochlorothiazide-, three lipid regulators -bezafibrate, fenofibrate and gemfibrozil-, and four psychiatric medications -carbamazepine, diazepam, lorazepam and paroxetine) were ascertained in a pilot-scale anaerobic/anoxic/aerobic (A²O) system treating urban wastewater, long term operated during two experimental phases using different sets of environmental conditions and operating parameters. Illumina MiSeq sequencing was used to investigate the structure, diversity and population dynamics of bacteria, archaea and fungi communities in the activated sludge. The results showed that mixed liquor suspended solids (MLSS) and food-to-microorganisms ratio (F/M) were operational parameters significantly influencing the REs of five of the analyzed PhACs in the A²O system. Biota-environment (BIO-ENV) analysis revealed strong correlations between population shifts of the activated sludge community and the REs of PhACs of the different pharmaceutical families. Increased REs of clarithromycin, furosemide, bezafibrate and gemfibrozil were concomitant to higher relative abundances of bacterial phylotypes classified within the Rhodobacteraceae and Sphingomonadaceae (Alphaproteobacteria), while those of Betaproteobacteria, Chloroflexi and Methanomethylovorans (Euryarchaea) correlated positively with the REs of up to seven PhACs belonging to different therapeutic groups.

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.

Modelling micropollutant fate in sewer systems - A new systematic approach to support conceptual model construction based on in-sewer hydraulic retention time

Conceptual sewer models are useful tools to assess the fate of micropollutants (MPs) in integrated wastewater systems. However, the definition of their model structure is highly subjective, and obtaining a realistic simulation of the in-sewer hydraulic retention time (HRT) is a major challenge without detailed hydrodynamic information or with limited measurements from the sewer network. This study presents an objective approach for defining the structure of conceptual sewer models in view of modelling MP fate in large urban catchments. The proposed approach relies on GIS-based information and a Gaussian mixture model to identify the model optimal structure, providing a multi-catchment conceptual model that accounts for HRT variability across urban catchment. This approach was tested in a catchment located in a highly urbanized Italian city and it was compared against a traditional single-catchment conceptual model (using a single average HRT) for the fate assessment of reactive MPs. Results showed that the multi-catchment model allows for a successful simulation of dry weather flow patterns and for an improved simulation of MP fate compared to the classical single-catchment model. Specifically, results suggested that a multi-catchment model should be preferred for (i) degradable MPs with half-life lower than the average HRT of the catchment and (ii) MPs undergoing formation from other compounds (e.g. human metabolites); or (iii) assessing MP loads entering the wastewater treatment plant from point sources, depending on their location in the catchment. Overall, the proposed approach is expected to ease the building of conceptual sewer models, allowing to properly account for HRT distribution and consequently improving MP fate estimation.

Large-scale determination of micropollutant elimination from municipal wastewater by passive sampling gives new insights in governing parameters and degradation patterns

A simple balancing method using passive samplers over a week's period has been developed and tested successfully to determine elimination rates of 22 common micropollutants of household and industrial sources in 18 full-scale wastewater treatment plants of different design and performance. Independent reactor tests to delineate elimination rates with native sludge of the treatment plants correlated very well with the full-scale elimination rate determinations. As opposed to common assumptions, this large dataset indicated that shorter sludge retention times - read: higher active biomass - showed higher micropollutant elimination rates in many cases. Multivariate statistical analysis of the elimination rates over the 18 treatment plants was able to group compounds according to common degradation pathways and showed that sensitivity to SRT drove the grouping. The dataset also allowed to determine population equivalent normalized loads of the investigated micropollutants. The application of WWTP balancing with passive sampling makes it relatively easy to gather elimination rates and inlet loads on a much broader basis than before and gives orientation for more in-depth analysis of degradation pathways.

Estimating the eco-toxicological risk of estrogens in China's rivers using a high-resolution contaminant fate model

The contamination of freshwater systems arises in many river basins due to industrialization and population growth, posing risks to ecosystems and human health. Despite these concerns, the fate and potential impact of many emerging pollutants are currently unknown, especially when the compounds are released into surface waters from populations distributed across large spatial scales. In order to address this shortcoming, a spatially-explicit contaminant fate model was developed as an extension of the global, vector-based river routing model HydroROUT. HydroROUT operates at very high spatial resolution (∼500 m), simulates river and stream chemical transport with in-stream removal, and contains links to a set of lakes and reservoirs, which act as a partial sink during the transport. The chemical fate model was applied to China and includes a consumption and release module based on county-level population demographics, considers point-source contributions from 2547 wastewater treatment plants, and accumulates contributions of rural and urban populations not connected to sewage treatment plants. As a case study, the sources and fates of the estrogens estrone (E₁), 17β-estradiol (E₂), estriol (E₃), as well as the synthetic estrogenic steroid hormone 17α-ethinylestradiol (EE₂) were modelled in Chinese surface water bodies. Preliminary validation of the results showed predictions to be within the ranges of concentrations reported in literature, with exception of EE₂. The total estrogenic mass in the entire river and lake system amounted to 370 tonnes of estrogens, with about 1.3 tonnes per year discharged to the ocean, neighboring countries or to inland sinks. Under a selected baseline scenario, eco-toxicological risk-i.e., contaminant concentrations in excess of the predicted no effect concentration (PNEC)-is found in 23.6% of all analyzed rivers of China with an average flow > 0.1 m³/s. Out of these, about 4% of rivers showed a high level of risk of 10 times or more above PNEC. Medium-to-large rivers are disproportionally affected; for example, 23.6%, 37.3%, 29.0% and 21.6% of river length are at risk in rivers of 1-10, 10-100, 100-1,000, and 1,000-10,000 m³/s of discharge, respectively, whereas no risk was predicted in the largest rivers (i.e., >10,000 m³/s) of China. Wastewater treatment plants process 22.5% of the total hormone load and thus play an important role in water quality control by reducing the risk in substantial portions of the river network, which would otherwise show elevated risk. Releases from untreated population dominate by far the overall contribution to risk.

Balancing environmental quality standards and infrastructure upgrade costs for the reduction of microcontaminant loads in rivers

Investments for upgrading wastewater treatment plants (WWTPs) with tertiary treatment to reduce microcontaminant loads in surface waters at a catchment scale can be daunting. These investments are highly sensitive to the selection of environmental quality standards (EQSs) for the target microcontaminants. Our hypothesis is that there is a balance between EQS selection and investment that needs to be considered in decision-making. We used a customized microcontaminant fate and transport model coupled to an optimization algorithm to validate this hypothesis in the Llobregat river basin and for the pharmaceutical compound diclofenac. The algorithm optimizes the number of WWTPs in this catchment requiring an upgrade to minimize the total amount of diclofenac that exceeds the EQS in every river section and the total cost. We simulated and optimized 40 scenarios representing a combination of 4 potential EQSs (10, 30, 50 and 100 ng L^-1), 5 levels of uncertainty bounds in the predictions of river concentrations and 2 hydrological scenarios (average flows, flows annually exceeding 30% of the days; and environmental flows, flows annually exceeding 99% of the days). The results showthat there is a nonlinear relationship between the EQS and the required investment. The investment increases by 100% from an EQS of 100 ng L^-1 to 10 ng L^-1, significantly increasing (by 60%) from 30 to 10 ng L^-1. Thus, establishing an EQS of 30 ng L^-1 would balance environmental protection and costs. The selection of the hydrological conditions also plays a key role in the upgrade analysis because the costs for environmental flows are 50% higher than for average flows. Finally, we highlight that the investment in research would allow the reduction of uncertainties, hence allowing more qualified decisions to be made and a reduction in the WWTP upgrade costs (up to 4 €·household^-1·year^-1).

Research progress on distribution, migration, transformation of antibiotics and antibiotic resistance genes (ARGs) in aquatic environment

Antimicrobial and antibiotics resistance caused by misuse or overuse of antibiotics exposure is a growing and significant threat to global public health. The spread and horizontal transfer of antibiotic resistant bacteria (ARB) and antibiotic resistance genes (ARGs) by the selective pressure of antibiotics in an aquatic environment is a major public health issue. To develop a better understanding of potential ecological risks die to antibiotics and ARGs, this study mainly summarizes research progress about: (i) the occurrence, concentration, fate, and potential ecological effects of antibiotics and ARGs in various aquatic environments, (ii) the threat, spread, and horizontal gene transfer (HGT) of ARGs, and (iii) the relationship between antibiotics, ARGs, and ARB. Finally, this review also proposes future research direction on antibiotics and ARGs.

Oxidation Processes in Water Treatment: Are We on Track?

Chemical oxidants have been applied in water treatment for more than a century, first as disinfectants and later to abate inorganic and organic contaminants. The challenge of oxidative abatement of organic micropollutants is the formation of transformation products with unknown (eco)toxicological consequences. Four aspects need to be considered for oxidative micropollutant abatement: (i) Reaction kinetics, controlling the efficiency of the process, (ii) mechanisms of transformation product formation, (iii) extent of formation of disinfection byproducts from the matrix, (iv) oxidation induced biological effects, resulting from transformation products and/or disinfection byproducts. It is impossible to test all the thousands of organic micropollutants in the urban water cycle experimentally to assess potential adverse outcomes of an oxidation. Rather, we need multidisciplinary and automated knowledge-based systems, which couple predictions of kinetics, transformation and disinfection byproducts and their toxicological consequences to assess the overall benefits of oxidation processes. A wide range of oxidation processes has been developed in the last decades with a recent focus on novel electricity-driven oxidation processes. To evaluate these processes, they have to be compared to established benchmark ozone- and UV-based oxidation processes by considering the energy demands, economics, the feasibilty, and the integration into future water treatment systems.

Management strategies for trace organic chemicals in water - A review of international approaches

To ensure an appropriate management of potential health risks and uncertainties from the release of trace organic chemicals (TOrCs) into the aqueous environment, many countries have evaluated and implemented strategies to manage TOrCs. The aim of this study was to evaluate existing management strategies for TOrCs in different countries to derive and compare underlying core principles and paradigms and to develop suggestions for more holistic management strategies to protect the environment and drinking water supplies from the discharge of undesired TOrCs. The strategies in different industrial countries were summarized and subsequently compared with regards to three particular questions: 1) Do the approaches different countries have implemented manage all or only specific portions of the universe of chemicals; 2) What implementation and compliance strategies are used to manage aquatic and human health risk and what are their pros and cons; and 3) How are site-specific watershed differences being addressed? While management strategies of the different countries target similar TOrCs, the programs differ in several important aspects, including underlying principles, the balance between aquatic or human health protection, implementation methods, and financing mechanisms used to fund regulatory programs.

Source estimation of pharmaceuticals based on catchment population and in-stream attenuation in Yodo River watershed, Japan

Fifty-five pharmaceuticals were monitored at four rivers and inlets and/or outlets of three sewage treatment plants (STPs) in Yodo River watershed, Japan over 17 sampling events. Twenty-six quantified pharmaceuticals were classified by source and fate. The load per person (LPP) of nine pharmaceuticals, including six with observed mass balance in studied river stretch of <80%, was appreciably lower in river water (RW) than in the effluent (EF) of STPs (RW/EF <0.5), indicating that they were susceptible to in-stream attenuation in the study area, while the others were relatively conservative. The LPP of 12 pharmaceuticals in RW were within ±50% of that in EF. Because their mass loadings in rivers were correlated with human population in the catchment and most people use the sewer system, the major source of the 12 pharmaceuticals was considered to be STPs. The LPP of the three most labile pharmaceuticals in STPs (caffeine, theophylline, and acetaminophen) was >1.5 in RW/EF and <1.0 in RW/influent (IF) of STPs. Poorly treated sewage discharged from households without using the sewer system was considered to be influential source of the three pharmaceuticals. The LPP (RW/EF) of caffeine, a pharmaceutical contained in food and beverage, was considerably higher than that of the other two, and this is attributable to untreated gray water discharged at households using the night-soil treatment system. The LPP of two veterinary drugs (sulfamonomethoxine and lincomycin) were >1.5 (RW/EF) and >1.0 (RW/IF). Their mass loadings in rivers showed a positive correlation with swine population in the catchment, although sulfamonomethoxine is equally used in both cattle and swine farming. This was attributable to application of cattle excrement as manure, and lability of sulfamonomethoxine during composting processes. The major source of the two veterinary drugs was considered to be on-site treatment facilities of swine urine.

Assessing watercourse quality: challenges in implementing European and Swiss legal frameworks

Hydro-climatic changes and the increasing release of pollutants into rivers by human activities tend to affect the quality of watercourses, to alter aquatic ecosystems and to reduce the amount of useable water. The ecological and chemical states of rivers and their evolution is thus of growing concern. In Europe and Switzerland, water policies are progressively shifting towards a holistic approach of river systems. The European Commission notably established a framework to highlight rivers' ecological deficits and to enhance regional or local water management plans. In Switzerland, a similar framework is currently under development. In this paper, both procedures are compared and implemented in a Swiss catchment dominated by agricultural activities. The aim is to define the challenges that still need to be addressed to assess and sustain river health. The hydromorphological, ecological, and ecotoxicological quality of the river was evaluated. Both frameworks highlighted the fact that no section of the river can currently be classified as being in a good environmental state and that the state deteriorates as tributaries and wastewater discharge flow into the main riverbed. Chemical issues and water quality changes due to hydro-climatic variations and management strategies were also pinpointed. Both frameworks are thus useful tools to survey changes in rivers quality in space and over time. However, challenges remain regarding the appropriate strategies to monitor and analyze chemicals, the definition of target values and conditions, the evaluation and integration of human-induced pressures, and the overall evaluation of the state of a river. The development of integrated indicators or of ecosystem services approaches is considered as a potential solution to explore river health and to define efficient restoration measures by water managers.

Incorporating model uncertainty into the evaluation of interventions to reduce microcontaminant loads in rivers

Models of microcontaminant fate and transport in wastewater treatment plants (WWTPs) and rivers have been developed and used to assist decision-making in the field of water management. These models come with parameter uncertainties that must be properly incorporated in the decision-making process. The main goal of this study is to evaluate how the magnitudes of key model parameter uncertainties influence the selection of end-of-pipe interventions (at WWTPs) designed to reduce the microcontaminant loads in rivers. We developed a model that describes the fate and removal of pharmaceuticals in WWTPs and the river network based on 3 key parameters: human pharmaceutical consumption and excretion (F) and the pharmaceutical degradation constants in WWTPs (k_WWTP) and rivers (k_river). We modelled the fate and transport of diclofenac in the Llobregat River basin (NE Spain). We calibrated the model using a Bayesian approach, which resulted in an accurate prediction of measured diclofenac loads at 9 locations along the Llobregat River and at the influents and effluents of 2 WWTPs (R² = 0.95). Using different scenarios, we evaluated three levels of uncertainty in the key model parameters. The first level of uncertainty corresponded to the reference distributions obtained from the Bayesian calibration. Then, for each parameter, we generated a narrower PDF (decreased uncertainty with respect to the reference) and a wider PDF (increased uncertainty). For each level of uncertainty, we evaluated increasing removal efficiencies of diclofenac at the WWTPs, from 38% to 98%. We assumed that removal efficiencies of up to 75% can be achieved by upgrading secondary treatment; beyond 75%, tertiary treatment is needed. The scenario analysis showed that achieving diclofenac removal efficiencies corresponding to tertiary treatment results in apparent concentration reductions (statistically significant differences relative to the reference situation), regardless of the level of uncertainty applied to the model parameters. However, upgrades in the secondary treatment resulted in apparent reductions only in the case of reduced uncertainty. We concluded that model uncertainty greatly influences the decisions that river basin authorities must make to reduce the microcontaminant loads released by WWTPs into rivers. In addition, we discussed research priorities to help reduce model uncertainty and thereby make more appropriate decisions.

Predictive framework for estimating exposure of birds to pharmaceuticals

We present and evaluate a framework for estimating concentrations of pharmaceuticals over time in wildlife feeding at wastewater treatment plants (WWTPs). The framework is composed of a series of predictive steps involving the estimation of pharmaceutical concentration in wastewater, accumulation into wildlife food items, and uptake by wildlife with subsequent distribution into, and elimination from, tissues. Because many pharmacokinetic parameters for wildlife are unavailable for the majority of drugs in use, a read-across approach was employed using either rodent or human data on absorption, distribution, metabolism, and excretion. Comparison of the different steps in the framework against experimental data for the scenario where birds are feeding on a WWTP contaminated with fluoxetine showed that estimated concentrations in wastewater treatment works were lower than measured concentrations; concentrations in food could be reasonably estimated if experimental bioaccumulation data are available; and read-across from rodent data worked better than human to bird read-across. The framework provides adequate predictions of plasma concentrations and of elimination behavior in birds but yields poor predictions of distribution in tissues. The approach holds promise, but it is important that we improve our understanding of the physiological similarities and differences between wild birds and domesticated laboratory mammals used in pharmaceutical efficacy/safety trials, so that the wealth of data available can be applied more effectively in ecological risk assessments. Environ Toxicol Chem 2017;36:2335-2344. © 2017 SETAC.

Assessing the population equivalent and performance of wastewater treatment through the ratios of pharmaceuticals and personal care products present in a river basin: Application to the River Thames basin, UK

The quality of surface waters in lowland rivers is largely dependent on the efficiency of wastewater treatment. Even in the developed countries, there have been difficulties in evaluating the effectiveness of wastewater management and the proportion of wastewater content (WWC) in the river, as well as in estimating the contributing human population. This study aimed to develop a wastewater quality and quantity assessment based on the occurrence of pharmaceuticals in the receiving waters. A survey of 53 pharmaceuticals in 324 samples (river water and influent and effluent of sewage (wastewater) treatment plants) was carried out in southern England in the River Thames catchment over four years. Carbamazepine was selected as stable marker and from its concentration WWC in the rivers and cumulative human populations along the catchment were estimated. The estimated population had a strong relationship (R²=0.94) with that reported by the local water company. The concentration ratio of the labile marker caffeine to carbamazepine indicated the efficiency of wastewater treatment in the different treatment systems (i.e. trickling filter or activated sludge) and in the receiving waters. The ratio in some river samples revealed unexpected discharges of untreated or poorly treated wastewater, with a total concentration of the analytes (up to 20μg/L) five times higher than that in treated wastewater. Such information could be valuable to estimate the discharge or occurrence of not only non-targeted chemicals, but also pathogens within the basin.

Removal of Antibiotics in Biological Wastewater Treatment Systems-A Critical Assessment Using the Activated Sludge Modeling Framework for Xenobiotics (ASM-X)

Many scientific studies present removal efficiencies for pharmaceuticals in laboratory-, pilot-, and full-scale wastewater treatment plants, based on observations that may be impacted by theoretical and methodological approaches used. In this Critical Review, we evaluated factors influencing observed removal efficiencies of three antibiotics (sulfamethoxazole, ciprofloxacin, tetracycline) in pilot- and full-scale biological treatment systems. Factors assessed include (i) retransformation to parent pharmaceuticals from e.g., conjugated metabolites and analogues, (ii) solid retention time (SRT), (iii) fractions sorbed onto solids, and (iv) dynamics in influent and effluent loading. A recently developed methodology was used, relying on the comparison of removal efficiency predictions (obtained with the Activated Sludge Model for Xenobiotics (ASM-X)) with representative measured data from literature. By applying this methodology, we demonstrated that (a) the elimination of sulfamethoxazole may be significantly underestimated when not considering retransformation from conjugated metabolites, depending on the type (urban or hospital) and size of upstream catchments; (b) operation at extended SRT may enhance antibiotic removal, as shown for sulfamethoxazole; (c) not accounting for fractions sorbed in influent and effluent solids may cause slight underestimation of ciprofloxacin removal efficiency. Using tetracycline as example substance, we ultimately evaluated implications of effluent dynamics and retransformation on environmental exposure and risk prediction.

Bromide Sources and Loads in Swiss Surface Waters and Their Relevance for Bromate Formation during Wastewater Ozonation

Bromide measurements and mass balances in the catchments of major Swiss rivers revealed that chemical industry and municipal waste incinerators are the most important bromide sources and account for ∼50% and ∼20%, respectively, of the ∼2000 tons of bromide discharged in the Rhine river in 2014 in Switzerland. About 100 wastewater treatment plants (WWTPs) will upgrade their treatment for micropollutant abatement in the future to comply with Swiss regulations. An upgrade with ozonation may lead to unintended bromate formation in bromide-containing wastewaters. Measured bromide concentrations were <0.05 mg L(-1) in ∼75% of 69 WWTPs, while they ranged from 0.4 to ∼50 mg L(-1) in WWTPs with specific bromide sources (e.g., municipal waste incinerators, landfill leachate, and chemical industry). Wastewater ozonation formed little bromate at specific ozone doses of ≤0.4 mg O3/mg DOC, while the bromate yields were almost linearly correlated to the specific ozone dose for higher ozone doses. Molar bromate yields for typical specific ozone doses in wastewater treatment (0.4-0.6 mg O3/mg DOC) are ≤3%. In a modeled extreme scenario (in which all upgraded WWTPs release 10 μg L(-1) of bromate), bromate concentrations increased by <0.4 μg L(-1) in major Swiss rivers and by several micrograms per liter in receiving water bodies with a high fraction of municipal wastewater.

Predicted and measured concentrations of pharmaceuticals in hospital effluents. Examination of the strengths and weaknesses of the two approaches through the analysis of a case study

This study deals with the chemical characterization of hospital effluents in terms of the predicted and measured concentrations of 38 pharmaceuticals belonging to 11 different therapeutic classes. The paper outlines the strengths and weaknesses of the two approaches through an analysis of a case study referring to a large hospital. It highlights the observed (and expected) ranges of variability for the parameters of the adopted model, presents the results of an uncertainty analysis of direct measurements (due to sampling mode and frequency and chemical analysis) and a sensitivity analysis of predicted concentrations (based on the annual consumption of pharmaceuticals, their excretion rate and annual wastewater volume generated by the hospital). Measured concentrations refer to two sampling campaigns carried out in summer and winter in order to investigate seasonal variability of the selected compounds. Predicted concentrations are compared to measured ones in the three scenarios: summer, winter and the whole year. It was found that predicted and measured concentrations are in agreement for a limited number of compounds (namely atenolol, atorvastatin and hydrochlorothiazide), and for most compounds the adoption of the model leads to a large overestimation in all three periods. Uncertainties in predictions are mainly due to the wastewater volume and excretion factor, whereas for measured concentrations, uncertainties are mainly due to sampling mode.

Management of pharmaceutical substances in the environment: Lithuanian case study

Investigation on the sources, discharges and related risks for the environment of the pharmaceutical substance (PhS) diclofenac (DCF) was performed in Lithuania, a country of the Baltic Sea region, for the first time. The investigation only refers to DCF as a PhS for human use; emissions from animal husbandry were not considered. In the first stage of the research, the main sources and pathways of DCF via substance flow analysis were identified within the country. During the second stage, DCF flows along the wastewater treatment plants (WWTPs) in two different cities were measured in order to assess the current levels of pharmaceutical residues in the environment. Furthermore, environmental risk assessment was carried out by taking into account the parameters of consumption data and elimination rate in WWTPs. Then, the assessment of different technical and managerial removal approaches was accomplished in an environmental management model of wastewater containing PhS, based on the framework of environmental systems theory.

Hospital-Use Pharmaceuticals in Swiss Waters Modeled at High Spatial Resolution

A model to predict the mass flows and concentrations of pharmaceuticals predominantly used in hospitals across a large number of sewage treatment plant (STP) effluents and river waters was developed at high spatial resolution. It comprised 427 geo-referenced hospitals and 742 STPs serving 98% of the general population in Switzerland. In the modeled base scenario, domestic, pharmaceutical use was geographically distributed according to the population size served by the respective STPs. Distinct hospital scenarios were set up to evaluate how the predicted results were modified when pharmaceutical use in hospitals was allocated differently; for example, in proportion to number of beds or number of treatments in hospitals. The hospital scenarios predicted the mass flows and concentrations up to 3.9 times greater than in the domestic scenario for iodinated X-ray contrast media (ICM) used in computed tomography (CT), and up to 6.7 times greater for gadolinium, a contrast medium used in magnetic resonance imaging (MRI). Field measurements showed that ICM and gadolinium were predicted best by the scenarios using number of beds or treatments in hospitals with the specific facilities (i.e., CT and/or MRI). Pharmaceuticals used both in hospitals and by the general population (e.g., cyclophosphamide, sulfamethoxazole, carbamazepine, diclofenac) were predicted best by the scenario using the number of beds in all hospitals, but the deviation from the domestic scenario values was only small. Our study demonstrated that the bed number-based hospital scenarios were effective in predicting the geographical distribution of a diverse range of pharmaceuticals in STP effluents and rivers, while the domestic scenario was similarly effective on the scale of large river-catchments.

Pharmaceutical concentrations in screened municipal wastewaters in Victoria, British Columbia: A comparison with prescription rates and predicted concentrations

Pharmaceuticals and personal care products (PPCPs) are emerging chemicals of concern detected in surface waters globally. Recent reviews advocate that PPCP occurrence, fate, and exposure need to be better predicted and characterized. The use of pharmaceutical prescription rates to estimate PPCP concentrations in the environment has been suggested. Concentrations of 7 pharmaceuticals (acetylsalicylic acid, diclofenac, fenoprofen, gemfibrozil, ibuprofen, ketoprofen, and naproxen) were measured in municipal wastewater using gas chromatography/ion trap-tandem mass spectroscopy (GC/IT-MS/MS). Subregional pharmaceutical prescription data were investigated to determine whether they could predict measured effluent concentrations (MECs) in wastewaters. Predicted effluent concentrations (PECs) for 5 of the 7 pharmaceuticals were within 2-fold agreement of the MECs when the fraction of parent pharmaceutical excreted was not considered. When the fraction of parent pharmaceutical excreted was considered, the respective PECs decreased, and most were within an order of magnitude of the MECs. Regression relationships of monthly PECs versus MECs were statistically significant (p < 0.05) but weak (R(2) = 0.18-0.56) for all pharmaceuticals except ketoprofen. This suggests high variability in the data and may be the result of factors influencing MECs such as the analytical methods used, wastewater sampling frequency, and methodology. The PECs were based solely on prescription rates and did not account for inputs of pharmaceuticals that had a significant over-the-counter component or were from other sources (e.g., hospitals).

Risk assessment of down-the-drain chemicals at large spatial scales: Model development and application to contaminants originating from urban areas in the Saint Lawrence River Basin

Chemicals released into freshwater systems threaten ecological functioning and may put aquatic life and the health of humans at risk. We developed a new contaminant fate model (CFM) that follows simple, well-established methodologies and is unique in its cross-border, seamless hydrological and geospatial framework, including lake routing, a critical component in northern environments. We validated the model using the pharmaceutical Carbamazepine and predicted eco-toxicological risk for 15 pharmaceuticals in the Saint-Lawrence River Basin, Canada. The results indicated negligible to low environmental risk for the majority of tested chemicals, while two pharmaceuticals showed elevated risk in up to 13% of rivers affected by municipal effluents. As an integrated model, our CFM is designed for application at very large scales with the primary goal of detecting high risk zones. In regulatory frameworks, it can help screen existing or new chemicals entering the market regarding their potential impact on human and environmental health. Due to its high geospatial resolution, our CFM can also facilitate the prioritization of actions, such as identifying regions where reducing contamination sources or upgrading treatment plants is most pertinent to achieve targeted pollutant removal or to protect drinking water resources.

Contaminants of Emerging Concern: Mass Balance and Comparison of Wastewater Effluent and Upstream Sources in a Mixed-Use Watershed

Understanding the sources, transport, and spatiotemporal variability of contaminants of emerging concern (CECs) is important for understanding risks and developing monitoring and mitigation strategies. This study used mass balances to compare wastewater treatment plant (WWTP) and upstream sources of 16 CECs to a mixed-use watershed in Minnesota, under different seasonal and hydrological conditions. Three distinct CEC groups emerged with respect to their source proportionality and instream behavior. Agricultural herbicides and daidzein inputs were primarily via upstream routes with the greatest loadings and concentrations during high flows. Trimethoprim, mecoprop, nonprescription pharmaceuticals, and personal care products entered the system via balanced/mixed pathways with peak loadings and concentrations in high flows. Carbaryl, 4-nonylphenol, and the remaining prescription pharmaceuticals entered the system via WWTP effluent with relatively stable loadings across sampling events. Mass balance analysis based on multiple sampling events and sites facilitated CEC source comparisons and may therefore prove to be a powerful tool for apportioning sources and exploring mitigation strategies.

Wastewater treatment--adsorption of organic micropollutants on activated HTC-carbon derived from sewage sludge

Organic micropollutants (MPs), in particular xenobiotics and their transformation products, have been detected in the aquatic environment and the main sources of these MPs are wastewater treatment plants. Therefore, an additional cleaning step is necessary. The use of activated carbon (AC) is one approach to providing this additional cleaning. Industrial AC derived from different carbonaceous materials is predominantly produced in low-income countries by polluting processes. In contrast, AC derived from sewage sludge by hydrothermal carbonization (HTC) is a regional and sustainable alternative, based on waste material. Our experiments demonstrate that the HTC-AC from sewage sludge was able to remove most of the applied MPs. In fact more than 50% of sulfamethoxazole, diclofenac and bezafibrate were removed from artificial water samples. With the same approach carbamazepine was eliminated to nearly 70% and atrazine more than 80%. In addition a pre-treated (phosphorus-reduced) HTC-AC was able to eliminate 80% of carbamazepine and diclofenac. Atrazine, sulfamethoxazole and bezafibrate were removed to more than 90%. Experiments using real wastewater samples with high organic content (11.1 g m(-3)) succeeded in proving the adsorption capability of phosphorus-reduced HTC-AC.

A review of hydrological and chemical stressors in the Adige catchment and its ecological status

Quantifying the effects of multiple stressors on Alpine freshwater ecosystems is challenging, due to the lack of tailored field campaigns for the contemporaneous measurement of hydrological, chemical and ecological parameters. Conducting exhaustive field campaigns is costly and hence most of the activities so far have been performed addressing specific environmental issues. An accurate analysis of existing information is therefore useful and necessary, to identify stressors that may act in synergy and to design new field campaigns. We present an extended review of available studies and datasets concerning the hydrological, chemical and ecological status of the Adige, which is the second longest river and the third largest river basin in Italy. The most relevant stressors are discussed in the light of the information extracted from a large number of studies. The detailed analysis of these studies identified that hydrological alterations caused by hydropower production are the main source of stress for the freshwater ecosystems in the Adige catchment. However, concurrent effects with other stressors, such as the release of pollutants from waste water treatment plants or from agricultural and industrial activities, have not been explored at depth, so far. A wealth of available studies address a single stressor separately without exploring their concurrent effect. It is concluded that a combination of extended experimental field campaigns, focusing on the coupled effects of multiple stressors, and modeling activities is highly needed in order to quantify the impact of the multifaceted human pressures on freshwater ecosystems in the Adige river.

Balancing the health benefits and environmental risks of pharmaceuticals: Diclofenac as an example

Pharmaceuticals are designed to improve human and animal health, but even the most beneficial pharmaceuticals might raise some questions concerning the consequences of exposure to non-target organisms. To illustrate this situation and using diclofenac as a case-study, we analyze global consumption and occurrence data to identify hot spots of consumption without occurrence data, review the scientific literature on the harmful environmental effects to determine whether the observed concentrations in freshwater are of environmental concern, summarize the current pharmaceutical and environmental policies to highlight policy gaps, and suggest a series of research and policy recommendations, which can be summarized as follows: we need to improve the current knowledge on occurrence in freshwaters to properly implement environmental policies (i), diclofenac might pose a risk to non-target organisms in freshwater (ii); the harmful effects that some pharmaceuticals may have on the environment are not always addressed by environmental policies (iii).

Risk of endocrine disruption to fish in the Yellow River catchment in China assessed using a spatially explicit model

The global water availability assessment (GWAVA) model, incorporating regional water abstractions and reservoir information, was used to model the human-sourced steroid estrogens estrone (E1) and estradiol (E2) in the Yellow River catchment (China). The river flows in the main stem were calibrated using gauged flows. Following a review of Chinese data on estrogen discharge from a range of sewage treatment plants, low, median, and high discharge rates were identified and used as best, expected, and worst-case scenarios, respectively. For any given location, the temporal variation of modeled estrogen levels was summarized using the mean and upper 90th percentile, which is where the model predicts 90% of values would be below this concentration. The predicted means and 90th percentiles for E1 were comparable to previous E1 measurements reported in the river. For the whole catchment, only 19% (mean value) of the river system by length was predicted to exceed 1 ng/L E2 equivalents (EEQs) using expected estrogen sewage discharge. Only 3% of the network by length was predicted to exceed the dangerously high 10 ng/L EEQ when considering 90th percentile concentrations. The highest exposures were in the Fen and Wei tributaries. Endocrine disruption risk from estrogens was predicted to be minimal in the main stem. Only in the worst-case discharge scenario and 90th percentile predicted concentrations were the most downstream river reaches of the main stem predicted to be at risk. Reservoirs appeared to be helpful in reducing estrogen concentrations thanks to longer water residence facilitating biodegradation.