Subscriber access provided by Kaohsiung Medical University
Characterization of Natural and Affected Environments
Can we swim yet? Systematic review, meta-analysis, and risk assessment of aging sewage in surface waters Alexandria B Boehm, Katherine E Graham, and Wiley Jennings Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01948 • Publication Date (Web): 06 Aug 2018 Downloaded from http://pubs.acs.org on August 6, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 39
1 2 3 4 5 6 7 8 9 10 11 12
Environmental Science & Technology
Can we swim yet? Systematic review, meta-analysis, and risk assessment of aging sewage in surface waters Alexandria B. Boehm*, Katherine E. Graham, Wiley C. Jennings Department of Civil and Environmental Engineering, Stanford University, Stanford, California, USA 94305 Submitted to Environmental Science & Technology
*
Corresponding author: 473 Via Ortega, Room 189, Dept. of Civil & Environmental Engineering , Stanford University, Stanford, CA 994305.
[email protected], (650) 724-9128 (tel)
ACS Paragon Plus Environment
1
Environmental Science & Technology
13 14
Page 2 of 39
TOC Art
15 16
ACS Paragon Plus Environment
2
Page 3 of 39
Environmental Science & Technology
17
Abstract
18
This study investigated the risk of gastrointestinal illness associated with swimming in surface
19
waters with aged sewage contamination. First, a systematic review compiled 333 first order
20
decay rate constants (k) for human norovirus and its surrogates feline calicivirus and murine
21
norovirus, Salmonella, Campylobacter, Escherichia coli O157:H7, Giardia, and
22
Cryptosporidium, and human-associated indicators in surface water. A meta-analysis
23
investigated effects of sunlight, temperature, and water matrix on k. There was a relatively large
24
number of k for bacterial pathogens and some human-associated indicators (n>40), fewer for
25
protozoan and (n=14-22), and few for human norovirus and its Caliciviridae surrogates (n=2-4).
26
Average k ranked: Campylobacter>human-associated markers>Salmonella> E. coli O157:H7 >
27
norovirus and its surrogates >Giardia>Cryptosporidium. Compiled k values were used in a
28
quantitative microbial risk assessment (QMRA) to simulate gastrointestinal illness risk
29
associated with swimming in water with aged sewage contamination. The QMRA used human-
30
associated fecal indicator HF183 as an index for the amount of sewage present and thereby
31
provided insight into how risk relates to HF183 concentrations in surface water. Because
32
exposure to norovirus contributed the majority of risk, and HF183 k is greater than norovirus k,
33
the risk associated with exposure to a fixed HF183 concentration increases with the age of
34
contamination. Swimmer exposure to sewage after it has aged ~3 days results in median risks
35
less than 30/1000. A risk-based water quality threshold for HF183 in surface waters that takes
36
into account uncertainty in contamination age is derived.
37
ACS Paragon Plus Environment
3
Environmental Science & Technology
Page 4 of 39
38
Introduction. Surface waters, including fresh, estuarine, and marine waters, serve as drinking
39
water sources, sites for recreation, sources of food, and essential organism habitat. They are
40
susceptible to microbial pollution1–5 from runoff6, animal feces7, and sewage discharges8.
41
Microbial pollutants include pathogenic protozoa, viruses, and bacteria. Microbial pollution
42
around the world is routinely assessed using concentrations of fecal indicator bacteria (FIB)
43
including Escherichia coli and enterococci as proxies for pathogenic organisms9. FIB presence in
44
bathing waters is linked quantitatively to gastrointestinal illness risk in swimmers when the FIB
45
source is runoff, sewage, or wastewater effluent10–14. However, FIB can be found in a variety of
46
sources including animal feces and environmental reservoirs15–17 making them non-ideal
47
indicators. Therefore, new indicators that are fecal-source associated have been recently
48
developed. For example, human-associated fecal indicators like HF183 and HumM218–21 are
49
highly abundant in, and specific to, human feces.
50 51
Once introduced to surface waters, microbial pollutants are transported and dispersed, and
52
removed from the water column by a variety of processes including inactivation22. Inactivation
53
of microbial pollutants in surface waters can occur via light and dark pathways. Inactivation via
54
dark pathways is caused by lack of nutrients, stress, senescence, and exposure to extracellular
55
biocidal compounds as well as grazing by bacterivorous organisms22,23. Inactivation can also
56
occur via light pathways where photons from the sun induce die-off directly (for example UVB
57
damaging genomic DNA) or indirectly via reactive species (generated when photons interact
58
with sensitizers like humic acids)24,25.
59
ACS Paragon Plus Environment
4
Page 5 of 39
Environmental Science & Technology
60
Several recent studies have compiled inactivation data on Escherichia coli and enterococci, as
61
well as human-associated indicators and pathogens in surface waters26–29. However, a systematic
62
screening of the literature on pathogen and human-associated indicators suggested that available
63
data are considerably more extensive than reported in these studies. A quantitative synthesis of
64
available inactivation data of waterborne pathogens and human indicators would facilitate the
65
modeling of these organisms in surface waters by providing a range of relevant decay rate
66
constants30. Such data would also be useful in quantitative microbial risk assessment (QMRA)
67
scenarios where fecal contamination may be aged31.
68 69
QMRAs have been used extensively to explore how risk of enteric illness varies with bathing
70
water exposure to microbial pollution. For example, QMRA was used to show that swimming in
71
water containing enterococci from different fecal sources, including gull feces, cow feces, and
72
sewage, is associated with different risks32 owing to the fact that different pathogens, with
73
different pathogenicity, are present in the different fecal sources. Other bathing water QMRA
74
applications allowed calculation of risk-based water quality thresholds for human-33 and gull-
75
associated indicators34 for recreational waters; when concentrations of the human- and gull-
76
associated indicators are below the risk-based thresholds, then water is safe for swimming. In
77
these QMRAs, however, the fecal source was assumed to be “unaged”, meaning that the ratio of
78
indicators to pathogens in the fecal source material, an important model input, was assumed to be
79
the same as that in the surface water. In reality, once introduced to the environment, the indicator
80
and pathogens may decay at different rates, thus altering their ratios.
81
ACS Paragon Plus Environment
5
Environmental Science & Technology
Page 6 of 39
82
In this study, a novel systematic review and meta-analysis of decay rate constants for key
83
pathogens and human-associated indicators was conducted. Target-specific decay rate constants
84
were summarized as statistical distributions so they may be readily used in a variety of modeling
85
applications. We subsequently used the decay rate constants in a QMRA to investigate how
86
exposure to fecal contamination from untreated sewage of different ages during swimming
87
affects the resultant gastrointestinal illness risk. The QMRA provides a framework for
88
identifying a risk-based water quality threshold for the human-associated fecal marker HF183.
89 90
Materials and Methods.
91
Systematic review. The systematic review and meta-analysis followed PRISMA guidelines35.
92
The goal of the review was to compile from the peer reviewed literature quantitative information
93
on the decay of waterborne pathogens and human-associated indicators of fecal pollution in
94
surface waters under environmentally-relevant conditions. Pathogens included in the review
95
were human norovirus, Campylobacter, Salmonella, Giardia, Cryptosporidium, and E. coli
96
O157:H7. Human norovirus surrogates in the Caliciviridae family, feline calicivirus and murine
97
norovirus, were also included. Human-associated fecal indicators included were BacHum-UCD,
98
HumM2, and HF183. These organisms were included as they were central to our planned QMRA
99
application.
100 101
Web of Science core collection (search field = topic), Scopus (search field = article title, abstract,
102
keyword), and PubMed (search field = all fields) were searched in August 2017 (Tables 1 and
103
S1). The search terms were “(X) AND (water OR seawater OR stormwater) AND (die-off OR
104
persistence OR survival OR inactivat* OR decay)” where X is the target-specific text (Table 1).
ACS Paragon Plus Environment
6
Page 7 of 39
Environmental Science & Technology
105
Identified articles were assembled and duplicates were removed. Details of the review process,
106
which involved two independent full text reviews of papers, are provided in the supporting
107
information (SI). The inclusion criteria were the paper: (1) contains quantitative data on the
108
decay of the target of interest in raw (unaltered) surface water, (2) is in English, (3) is not a
109
review, presents primary data, and is peer reviewed, (4) does not contain data solely on
110
disinfection treatments such as addition of oxidants or SODIS, (5) included data from decay
111
experiments where the temperature is greater than 4°C and less than 30°C, and (6) describes
112
methods to enumerate the target that are logical and justifiable.
113 114
Decay rate constants were extracted from papers by a single reviewer. First-order decay rate
115
constants (k), in units per day (d-1), calculated from natural log (ln)-transformed concentration
116
data as used in Chick’s law36, were sought. If a study presented k values, then they were
117
extracted from the paper along with any reported errors and model fit values (R2 and/or root
118
mean square error (RMSE)), and unit conversions were applied where appropriate. If a study
119
reported decay model parameters from a model that was not a first-order model (for example a
120
shoulder log-linear model, or biphasic model), then we extracted those reported model
121
parameters and any associated errors and model fit values. If a study only reported T90 or T99
122
values (time to 90% of 99% reduction in concentration, respectively), then they were converted
123
to first order decay rate constants assuming Chick’s law applied. If no first order decay rate
124
constant was reported by the study authors, then Plot Digitizer
125
(http://plotdigitizer.sourceforge.net) was used to digitize the concentration times series appearing
126
in graphs within the publication. To be clear, this included data from studies that only reported
127
decay model parameters from other types (not first-order log-linear) of decay models. k was then
ACS Paragon Plus Environment
7
Environmental Science & Technology
Page 8 of 39
128
calculated as the regression slope of ln(C/Co) versus time (in days) using linear least-squares
129
regression in R. In this formulation C is the concentration at time t, and Co is the concentration at
130
the start of the experiment at t=0. k and its associated error as well as model fit parameters were
131
recorded. In carrying out the linear regression, values reported at or below the detection limit
132
were included if and only if they were not preceded by other consecutive values at or below the
133
detection limit; the value directly reported by the author was used in these cases.
134 135
Once all data were compiled, datasets and model parameters were examined to assess whether a
136
non-linear model was needed to describe decay. The goodness of log-linear model fit to the data
137
(R2 and RMSE), and the number of data points that appeared to “deviate” from the log-linear
138
model were considered. In general, if R2 values were greater than 0.7 and RMSE was relatively
139
small (~1 ln unit), only one data point visually deviated from a straight line fit between time and
140
ln(C/Co), or the non-log-linear model fit was no better than the log-linear fit, then a log-linear
141
curve fit was deemed acceptable.
142 143
In addition to extracting information on the decay of the target of interest, whether the
144
experiment was conducted in (1) freshwater, estuarine water, or seawater and (2) natural sunlight
145
or the dark was also noted. If an experiment was reportedly carried out in sunlight, but at a depth
146
in the water column greater than ~25 cm or in a container that was opaque to UVA and UVB, the
147
experiment was categorized as carried out in the dark. This is justified by previous work
148
suggesting that such experiments do not receive sufficient photons to be deemed affected by
149
sunlight37. The temperature at which the experiment was conducted was also recorded. If a range
150
of temperatures was provided, the mean of the reported range was used. Finally, the method used
ACS Paragon Plus Environment
8
Page 9 of 39
Environmental Science & Technology
151
for target enumeration was noted (culture, microscopy, quantitative PCR (QPCR), reverse-
152
transcription QPCR (RT-QPCR), or propidium monoazide QPCR (PMA-QPCR)).
153 154
Twenty percent of the papers from which data were extracted by a single reviewer were
155
randomly chosen for a second round of data extraction by a different reviewer. Data extracted by
156
the two reviewers were compared to ensure consistency. A single author conducted detailed
157
review of all datasets to identify missing data, data outliers, and data entry mistakes.
158 159
Meta-analysis. Statistical distributions were fit to target-specific k values. Goodness of fit was
160
assessed by visual inspection of residual and Q-Q plots. This yielded log-normal fits for all
161
organisms/targets with the number of k values n>14. For congruity, log-normal distributions
162
were also used when n