Greywater reuse is a strategy to address water scarcity, necessitating the selection of treatment processes that balance cost-efficiency and human health risks. A key aspect in evaluating these risks is understanding pathogen contamination levels in greywater, a complex task due to intermittent pathogen occurrences. To address this, faecal indicator organisms like E. coli are often monitored as proxies to evaluate faecal contamination levels and infer pathogen concentrations. However, the wide variability in faecal indicator concentrations poses challenges in their modelling for quantitative microbial risk assessment (QMRA). Our study critically assesses the adequacy of parametric models in predicting the variability in E. coli concentrations in greywater. We found that models that build on summary statistics, like medians and standard deviations, can substantially underestimate the variability in E. coli concentrations. More appropriate models may provide more accurate estimations of, and uncertainty around, peak E. coli concentrations. To demonstrate this, a Poisson lognormal distribution model is fit to a data set of E. coli concentrations measured in shower and laundry greywater sources. This model estimated arithmetic mean E. coli concentrations in laundry waters at approximately 1.0E + 06 MPN 100 mL−1. These results are around 2.0 log10 units higher than estimations from a previously used hierarchical lognormal model based on aggregated summary data from multiple studies. Such differences are considerable when assessing human health risks and setting pathogen reduction targets for greywater reuse. This research highlights the importance of making raw monitoring data available for more accurate statistical evaluations than those based on summary statistics. It also emphasizes the crucial role of model comparison, selection, and validation to inform policy-relevant outcomes.@en

In light of increasingly diverse greywater reuse applications, this study proposes risk-based log-removal targets (LRTs) to aid the selection of treatment trains for greywater recycling at different collection scales, including appliance-scale reuse of individual greywater streams. An epidemiology-based model was used to simulate the concentrations of prevalent and treatment-resistant reference pathogens (protozoa: Giardia and Cryptosporidium spp., bacteria: Salmonella and Campylobacter spp., viruses: rotavirus, norovirus, adenovirus, and Coxsackievirus B5) in the greywater streams for collection scales of 5-, 100-, and a 1000-people. Using quantitative microbial risk assessment (QMRA), we calculated LRTs to meet a health benchmark of 10^–4 infections per person per year over 10′000 Monte Carlo iterations. LRTs were highest for norovirus at the 5-people scale and for adenovirus at the 100- and 1000-people scales. Example treatment trains were designed to meet the 95 % quantiles of LRTs. Treatment trains consisted of an aerated membrane bioreactor, chlorination, and, if required, UV disinfection. In most cases, rotavirus, norovirus, adenovirus and Cryptosporidium spp. determined the overall treatment train requirements. Norovirus was most often critical to dimension the chlorination (concentration × time values) and adenovirus determined the required UV dose. Smaller collection scales did not generally allow for simpler treatment trains due to the high LRTs associated with viruses, with the exception of recirculating washing machines and handwashing stations. Similarly, treating greywater sources individually resulted in lower LRTs, but the lower required LRTs nevertheless did not generally allow for simpler treatment trains. For instance, LRTs for a recirculating washing machine were around 3-log units lower compared to LRTs for indoor reuse of combined greywater (1000-people scale), but both scenarios necessitated treatment with a membrane bioreactor, chlorination and UV disinfection. However, simpler treatment trains may be feasible for small-scale and application-scale reuse if: (i) less conservative health benchmarks are used for household-based systems, considering the reduced relative importance of treated greywater in pathogen transmission in households, and (ii) higher log-removal values (LRVs) can be validated for unit processes, enabling simpler treatment trains for a larger number of appliance-scale reuse systems.

@en

Intense rainfall and snowmelt events may affect the safety of drinking water, as large quantities of fecal material can be discharged from storm or sewage overflows or washed from the catchment into drinking water sources. This study used β-d-glucuronidase activity (GLUC) with microbial source tracking (MST) markers: human, bovine, porcine mitochondrial DNA markers (mtDNA) and human-associated Bacteroidales HF183 and chemical source tracking (CST) markers including caffeine, carbamazepine, theophylline and acetaminophen, pathogens (Giardia, Cryptosporidium, adenovirus, rotavirus and enterovirus), water quality indicators (Escherichia coli, turbidity) and hydrometeorological data (flowrate, precipitation) to assess the vulnerability of 3 drinking water intakes (DWIs) and identify sources of fecal contamination. Water samples were collected under baseline, snow and rain events conditions in urban and agricultural catchments (Québec, Canada). Dynamics of E. coli, HF183 and WWMPs were similar during contamination events, and concentrations generally varied over 1 order of magnitude during each event. Elevated human-associated marker levels during events demonstrated that urban DWIs were impacted by recent contamination from an upstream municipal water resource recovery facility (WRRF). In the agricultural catchment, mixed fecal pollution was observed with the occurrences and increases of enteric viruses, human bovine and porcine mtDNA during peak contaminating events. Bovine mtDNA qPCR concentrations were indicative of runoff of cattle-derived fecal pollutants to the DWI from diffuse sources following rain events. This study demonstrated that the suitability of a given MST or CST indicator depend on river and catchment characteristics. The sampling strategy using continuous online GLUC activity coupled with MST and CST markers analysis was a more reliable source indicator than turbidity to identify peak events at drinking water intakes.@en

Antimicrobial resistance (AMR) poses a global health threat, causing millions of deaths annually, with expectations of increased impact in the future. Wastewater surveillance offers a cost-effective, non-invasive tool to understand AMR carriage trends within a population. We monitored extended-spectrum β-lactamase producing Escherichia coli (ESBL-E. coli) weekly in influent wastewater from six wastewater treatment plants (WWTPs) in Switzerland (November 2021 to November 2022) to investigate spatio-temporal variations, explore correlations with environmental variables, develop a predictive model for ESBL-E. coli carriage in the community, and detect the most prevalent ESBL-genes. We cultured total and ESBL-E. coli in 300 wastewater samples to quantify daily loads and percentage of ESBL-E. coli. Additionally, we screened 234 ESBL-E. coli isolates using molecular methods for the presence of 18 ESBL-gene families. We found a population-weighted mean percentage of ESBL-E. coli of 1.9% (95% confidence interval: 1.8–2%) across all sites and weeks, which can inform ESBL-E. coli carriage. Concentrations of ESBL-E. coli varied across WWTPs and time, with higher values observed in WWTPs serving larger populations. Recent precipitations (previous 24/96 h) showed no significant association with ESBL-E. coli, while temperature occasionally had a moderate impact (P < 0.05, correlation coefficients approximately 0.40) in some locations. We identified blaCTX-M-1, blaCTX-M-9, and blaTEM as the predominant ESBL-gene families. Our study demonstrates that wastewater-based surveillance of culturable ESBL-E. coli provides insights into AMR trends in Switzerland and may also inform resistance. These findings establish a foundation for long term, nationally established monitoring protocols and provide information that may help inform targeted public health interventions.@en

Enteric pathogen infections are a leading cause of morbidity and mortality globally, with the highest disease burden in low-income countries. Hands act as intermediaries in enteric pathogen transmission, transferring enteric pathogens between people and the environment through contact with fomites, food, water, and soil. In this study, we conducted a systematic review of prevalence and concentrations of fecal indicator microorganisms (i.e., E. coli, fecal coliform) and enteric pathogens on hands. We identified 84 studies, reporting 35,440 observations of hand contamination of people in community or household settings. The studies investigated 44 unique microorganisms, of which the most commonly reported indicators were E. coli and fecal coliforms. Hand contamination with 12 unique enteric pathogens was reported, with adenovirus and norovirus as the most frequent. Mean E. coli prevalence on hands was 62% [95% CI 40%–82%] and mean fecal coliform prevalence was 66% [95% CI 22%–100%]. Hands were more likely to be contaminated with E. coli in low/lower-middle-income countries (prevalence: 69% [95% CI 48%–88%]) than in upper-middle/high-income countries (6% [95% CI 2%–12%]). The Review also highlights the importance of standardizing hand sampling methods, as hand rinsing was associated with greater fecal contamination compared to other sampling methods.@en

Waterborne enteric viruses in lakes, especially at recreational water sites, may have a negative impact on human health. However, their fate and transport in lakes are poorly understood. In this study, we propose a coupled water quality and quantitative microbial risk assessment (QMRA) model to study the transport, fate and infection risk of four common waterborne viruses (adenovirus, enterovirus, norovirus and rotavirus), using Lake Geneva as a study site. The measured virus load in raw sewage entering the lake was used as the source term in the water quality simulations for a hypothetical scenario of discharging raw wastewater at the lake surface. After discharge into the lake, virus inactivation was modeled as a function of water temperature and solar irradiance that varied both spatially and temporally during transport throughout the lake. Finally, the probability of infection, while swimming at a popular beach, was quantified and compared among the four viruses. Norovirus was found to be the most abundant virus that causes an infection probability that is at least 10 times greater than the other viruses studied. Furthermore, environmental inactivation was found to be an essential determinant in the infection risks posed by viruses to recreational water users. We determined that infection risks by enterovirus and rotavirus could be up to 1000 times lower when virus inactivation by environmental stressors was accounted for compared with the scenarios considering hydrodynamic transport only. Finally, the model highlighted the role of the wind field in conveying the contamination plume and hence in determining infection probability. Our simulations revealed that for beaches located west of the sewage discharge, the infection probability under eastward wind was 43% lower than that under westward wind conditions. This study highlights the potential of combining water quality simulation and virus-specific risk assessment for a safe water resources usage and management.@en

Preventing failures of water treatment barriers can play an important role in meeting the increasing demand for microbiologically safe water. The development and integration of failure prevention strategies into quantitative microbial risk assessment (QMRA) offer opportunities to support the design and operation of treatment trains. This study presents existing failure models and extends them to guide the development of risk-based operational monitoring strategies. For barriers with rapid performance loss, results show that a failure of 15 s should be reliably detected to verify a log reduction value (LRV) of 6.0; thus, detecting and remediating these failures may be beyond current technology. For chemical disinfection with a residual, failure durations in order of minutes should be reliably detected to verify a LRV of 6.0. Short-term failures are buffered because the disinfectant residual concentration sustains a partial reduction performance. Therefore, increasing the contact time and hydraulic mixing reduces the impact of failures. These findings demonstrate the importance of defining precise frequencies to monitor barrier performances during operation. Overall, this study highlights the utility of process-specific models for developing failure prevention strategies for water safety management.@en

Legionella are natural inhabitants of building plumbing biofilms, where interactions with other microorganisms influence their survival, proliferation, and death. Here, we investigated the associations of Legionella with bacterial and eukaryotic microbiomes in biofilm samples extracted from 85 shower hoses of a multiunit residential building. Legionella spp. relative abundance in the biofilms ranged between 0–7.8%, of which only 0–0.46% was L. pneumophila. Our data suggest that some microbiome members were associated with high (e.g. Chthonomonas, Vrihiamoeba) or low (e.g. Aquabacterium, Vannella) Legionella relative abundance. The correlations of the different Legionella variants (30 Zero-Radius OTUs detected) showed distinct patterns, suggesting separate ecological niches occupied by different Legionella species. This study provides insights into the ecology of Legionella with respect to: (i) the colonization of a high number of real shower hoses biofilm samples; (ii) the ecological meaning of associations between Legionella and co-occurring bacterial/eukaryotic organisms; (iii) critical points and future directions of microbial-interaction-based-ecological-investigations.@en

Minimum treatment requirements are set in response to established or anticipated levels of enteric pathogens in the source water of drinking water treatment plants (DWTPs). For surface water, contamination can be determined directly by monitoring reference pathogens or indirectly by measuring fecal indicators such as Escherichia coli (E. coli). In the latter case, a quantitative interpretation of E. coli for estimating reference pathogen concentrations could be used to define treatment requirements. This study presents the statistical analysis of paired E. coli and reference protozoa (Cryptosporidium, Giardia) data collected monthly for two years in source water from 27 DWTPs supplied by rivers in Canada. E. coli/Cryptosporidium and E. coli/Giardia ratios in source water were modeled as the ratio of two correlated lognormal variables. To evaluate the potential of E. coli for defining protozoa treatment requirements, risk-based critical mean protozoa concentrations in source water were determined with a reverse quantitative microbial risk assessment (QMRA) model. Model assumptions were selected to be consistent with the World Health Organization (WHO) Guidelines for drinking-water quality. The sensitivity of mean E. coli concentration trigger levels to identify these critical concentrations in source water was then evaluated. Results showed no proportionalities between the log of mean E. coli concentrations and the log of mean protozoa concentrations. E. coli/protozoa ratios at DWTPs supplied by small rivers in agricultural and forested areas were typically 1.0 to 2.0-log lower than at DWTPs supplied by large rivers in urban areas. The seasonal variations analysis revealed that these differences were related to low mean E. coli concentrations during winter in small rivers. To achieve the WHO target of 10⁻⁶ disability-adjusted life year (DALY) per person per year, a minimum reduction of 4.0-log of Cryptosporidium would be required for 20 DWTPs, and a minimum reduction of 4.0-log of Giardia would be needed for all DWTPs. A mean E. coli trigger level of 50 CFU 100 mL⁻¹ would be a sensitive threshold to identify critical mean concentrations for Cryptosporidium but not for Giardia. Treatment requirements higher than 3.0-log would be needed at DWTPs with mean E. coli concentrations as low as 30 CFU 100 mL⁻¹ for Cryptosporidium and 3 CFU 100 mL⁻¹ for Giardia. Therefore, an E. coli trigger level would have limited value for defining health-based treatment requirements for protozoa at DWTPs supplied by small rivers in rural areas.

@en

A monitoring strategy was implemented at two drinking water treatment plants in Quebec, Canada, to evaluate microbial reduction performances of full-scale treatment processes under different source water conditions. β-D-glucuronidase activity in source water was automatically monitored in near-real-time to establish baseline and event conditions at each location. High-volume water samples (50–1,500 L) were collected at the inflow and the outflow of coagulation/flocculation, filtration, and UV disinfection processes and were analysed for two naturally occurring surrogate organisms: Escherichia coli and Clostridium perfringens. Source water Cryptosporidium data and full-scale C. perfringens reduction data were entered into a quantitative microbial risk assessment (QMRA) model to estimate daily infection risks associated with exposures to Cryptosporidium via consumption of treated drinking water. Daily mean E. coli and Cryptosporidium concentrations in source water under event conditions were in the top 5% (agricultural site) or in the top 15% (urban site) of what occurs through the year at these drinking water treatment plants. Reduction performances of up to 6.0-log for E. coli and 5.6-log for C. perfringens were measured by concentrating high-volume water samples throughout the treatment train. For both drinking water treatment plants, removal performances by coagulation/flocculation/sedimentation processes were at the high end of the range of those reported in the literature for bacteria and bacterial spores. Reductions of E. coli and C. perfringens by floc blanket clarification, ballasted clarification and rapid sand filtration did not deteriorate during two snowmelt/rainfall events. QMRA results suggested that daily infection risks were similar during two rainfall/snowmelt events than during baseline conditions. Additional studies investigating full-scale reductions would be desirable to improve the evaluation of differences in treatment performances under various source water conditions.

@en

This study investigates short-term fluctuations in virus concentrations in source water and their removal by full-scale drinking water treatment processes under different source water conditions. Transient peaks in raw water faecal contamination were identified using in situ online β-D-glucuronidase activity monitoring at two urban drinking water treatment plants. During these peaks, sequential grab samples were collected at the source and throughout the treatment train to evaluate concentrations of rotavirus, adenovirus, norovirus, enterovirus, JC virus, reovirus, astrovirus and sapovirus by reverse transcription and real-time quantitative PCR. Virus infectivity was assessed through viral culture by measurement of cytopathic effect and integrated cell culture qPCR. Virus concentrations increased by approximately 0.5-log during two snowmelt/rainfall episodes and approximately 1.0-log following a planned wastewater discharge upstream of the drinking water intake and during a β-D-glucuronidase activity peak in dry weather conditions. Increases in the removal of adenovirus and rotavirus by coagulation/flocculation processes were observed during peak virus concentrations in source water, suggesting that these processes do not operate under steady-state conditions but dynamic conditions in response to source water conditions. Rotavirus and enterovirus detected in raw and treated water samples were predominantly negative in viral culture. At one site, infectious adenoviruses were detected in raw water and water treated by a combination of ballasted clarification, ozonation, GAC filtration, and UV disinfection operated at a dose of 40 mJ cm⁻². The proposed sampling strategy can inform the understanding of the dynamics associated with virus concentrations at drinking water treatment plants susceptible to de facto wastewater reuse.

@en

Waterborne pathogens are heterogeneously distributed across various spatiotemporal scales in water resources, and representative sampling is therefore crucial for accurate risk assess-ment. Since regulatory monitoring of microbiological water quality is usually conducted at fixed time intervals, it can miss short-term fecal contamination episodes and underestimate underlying microbial risks. In the present paper, we developed a new automated sampling methodology based on near real-time measurement of a biochemical indicator of fecal pollution. Online monitoring of β-D-glucuronidase (GLUC) activity was used to trigger an automated sampler during fecal contamination events in a drinking water supply and at an urban beach. Significant increases in protozoan parasites, microbial source tracking markers and E. coli were measured during short-term (<24 h) fecal pollution episodes, emphasizing the intermittent nature of their occurrence in water. Synchro-nous triggering of the automated sampler with online GLUC activity measurements further re-vealed a tight association between the biochemical indicator and culturable E. coli. The proposed event sampling methodology is versatile and in addition to the two triggering modes validated here, others can be designed based on specific needs and local settings. In support to regulatory monitoring schemes, it should ultimately help gathering crucial data on waterborne pathogens more efficiently during episodic fecal pollution events.

@en

In several jurisdictions, the arithmetic mean of Escherichia coli concentrations in raw water serves as the metric to set minimal treatment requirements by drinking water treatment plants (DWTPs). An accurate and precise estimation of this mean is therefore critical to define adequate requirements. Distributions of E. coli concentrations in surface water can be heavily skewed and require statistical methods capable of characterizing uncertainty. We present four simple parametric models with different upper tail behaviors (gamma, log-normal, Lomax, mixture of two log-normal distributions) to explicitly account for the influence of peak events on the mean concentration. The performance of these models was tested using large E. coli data sets (200–1800 samples) from raw water regulatory monitoring at six DWTPs located in urban and agricultural catchments. Critical seasons of contamination and hydrometeorological factors leading to peak events were identified. Event-based samples were collected at an urban DWTP intake during two hydrometeorological events using online β-D-glucuronidase activity monitoring as a trigger. Results from event-based sampling were used to verify whether selected parametric distributions predicted targeted peak events. We found that the upper tail of the log-normal and the Lomax distributions better predicted large concentrations than the upper tail of the gamma distribution. Weekly sampling for two years in urban catchments and for four years in agricultural catchments generated reasonable estimates of the average raw water E. coli concentrations. The proposed methodology can be easily used to inform the development of sampling strategies and statistical indices to set site-specific treatment requirements.

@en

The identification of appropriately conservative statistical distributions is needed to predict microbial peak events in drinking water sources explicitly. In this study, Poisson and mixed Poisson distributions with different upper tail behaviors were used for modeling source water Cryptosporidium and Giardia data from 30 drinking water treatment plants. Small differences (<0.5-log) were found between the “best” estimates of the mean Cryptosporidium and Giardia concentrations with the Poisson–gamma and Poisson–log-normal models. However, the upper bound of the 95% credibility interval on the mean Cryptosporidium concentrations of the Poisson–log-normal model was considerably higher (>0.5-log) than that of the Poisson–gamma model at four sites. The improper choice of a model may, therefore, mislead the assessment of treatment requirements and health risks associated with the water supply. Discrimination between models using the marginal deviance information criterion (mDIC) was unachievable because differences in upper tail behaviors were not well characterized with available data sets ((Formula presented.)). Therefore, the gamma and the log-normal distributions fit the data equally well but may predict different risk estimates when they are used as an input distribution in an exposure assessment. The collection of event-based monitoring data and the modeling of larger routine monitoring data sets are recommended to identify appropriately conservative distributions to predict microbial peak events.

@en

Temporal variations in concentrations of pathogenic microorganisms in surface waters are well known to be influenced by hydrometeorological events. Reasonable methods for accounting for microbial peaks in the quantification of drinking water treatment requirements need to be addressed. Here, we applied a novel method for data collection and model validation to explicitly account for weather events (rainfall, snowmelt) when concentrations of pathogens are estimated in source water. Online in situ β-d-glucuronidase activity measurements were used to trigger sequential grab sampling of source water to quantify Cryptosporidium and Giardia concentrations during rainfall and snowmelt events at an urban and an agricultural drinking water treatment plant in Quebec, Canada. We then evaluate if mixed Poisson distributions fitted to monthly sampling data ((Formula presented.) = 30 samples) could accurately predict daily mean concentrations during these events. We found that using the gamma distribution underestimated high Cryptosporidium and Giardia concentrations measured with routine or event-based monitoring. However, the log-normal distribution accurately predicted these high concentrations. The selection of a log-normal distribution in preference to a gamma distribution increased the annual mean concentration by less than 0.1-log but increased the upper bound of the 95% credibility interval on the annual mean by about 0.5-log. Therefore, considering parametric uncertainty in an exposure assessment is essential to account for microbial peaks in risk assessment.

@en

Urban source water protection planning requires the characterization of sources of contamination upstream of drinking water intakes. Elevated pathogen concentrations following Combined Sewer Overflows (CSOs) represent a threat to human health. Quantifying peak pathogen concentrations at the intakes of drinking water plants is a challenge due to the variability of CSO occurrences and uncertainties with regards to the fate and transport mechanisms from discharge points to source water supplies. Here, a two-dimensional deterministic hydrodynamic and water quality model is used to study the fluvial contaminant transport and the impacts of the upstream CSO discharges on the downstream concentrations of Escherichia coli in the raw water supply of two drinking water plants, located on a large river. CSO dynamic loading characteristics were considered for a variety of discharges. As a result of limited Cryptosporidium data, a probability distribution of the ratio of E. coli to Cryptosporidium based on historical data was used to estimate microbial risk from simulated CSO-induced E. coli concentrations. During optimal operational performance of the plants, the daily risk target was met (based on the mean concentration during the peak) for 80% to 90% of CSO events. For suboptimal performance of the plants, these values dropped to 40% to 55%. Mean annual microbial risk following CSO discharge events was more dependent on treatment performance rather than the number of CSO occurrences. The effect of CSO-associated short term risk on the mean annual risk is largely dependent on the treatment performance as well as representativeness of the baseline condition at the intakes, demonstrating the need for assessment of treatment efficacy. The results of this study will enable water utilities and managers with a tool to investigate the potential alternatives in reducing the microbial risk associated with CSOs.@en

Past waterborne outbreaks have demonstrated that informed vulnerability assessment of drinking water supplies is paramount for the provision of safe drinking water. Although current monitoring frameworks are not designed to account for short-term peak concentrations of fecal microorganisms in source waters, the recent development of online microbial monitoring technologies is expected to fill this knowledge gap. In this study, online near real-time monitoring of β-d-glucuronidase (GLUC) activity was conducted for 1.5 years at an urban drinking water intake impacted by multiple point sources of fecal pollution. Parallel routine and event-based monitoring of E. coli and online measurement of physico-chemistry were performed at the intake and their dynamics compared over time. GLUC activity fluctuations ranged from seasonal to hourly time scales. All peak contamination episodes occurred between late fall and early spring following intense rainfall and/or snowmelt. In the absence of rainfall, recurrent daily fluctuations in GLUC activity and culturable E. coli were observed at the intake, a pattern otherwise ignored by regulatory monitoring. Cross-correlation analysis of time series retrieved from the drinking water intake and an upstream Water Resource Recovery Facility (WRRF) demonstrated a hydraulic connection between the two sites. Sewage by-passes from the same WRRF were the main drivers of intermittent GLUC activity and E. coli peaks at the drinking water intake following intense precipitation and/or snowmelt. Near real-time monitoring of fecal pollution through GLUC activity enabled a thorough characterization of the frequency, duration and amplitude of peak contamination periods at the urban drinking water intake while providing crucial information for the identification of the dominant upstream fecal pollution sources. To the best of our knowledge, this is the first characterization of a hydraulic connection between a WRRF and a downstream drinking water intake across hourly to seasonal timescales using high frequency microbial monitoring data. Ultimately, this should help improve source water protection through catchment mitigation actions, especially in a context of de facto wastewater reuse.@en