1. Introduction
Detection of Escherichia coli (E. coli) in natural freshwater bodies is used as an indicator of fecal contamination in assessing environmental water quality worldwide. Such indicators positively correlate with fecally-derived pathogens (e.g., Salmonella, Campylobacter, Giardia, enteric viruses, etc.) and have been used by both public and private environmental monitoring agencies for the last 100 years.
1.1. Studies of Indicator Distribution in Natural Waters
To date, many reports have been published on the vertical and horizontal variation of E. coli density and other indicator bacteria within standing water bodies such as freshwater lakes, reservoirs, and water/sewage treatment ponds. Environmental parameters such as temperature, oxygen, and/or nutrient availability [1] -[4] , sediment deposits [1] [5] -[7] , solar inactivation and insolation [8] -[10] , predation [3] , presence of biological host(s) [11] -[13] , input from watershed soils and sediments [14] [15] , nucleic acid content [16] , and variety of bacterial indicator assessed [17] have been examined with respect to depth in water column in these water bodies. Still other reports have examined stratification of viral populations in freshwater lakes [18] and estuaries [19] .
Although there exists information about bacteria stratification in lakes [4] and a recent study of estuarine water reports variation in E. coli numbers with respect to lateral location and bank side [20] , distributions of E. coli and/or other indicator bacteria in freshwater streams and rivers (i.e., moving freshwater bodies) have received little attention [4] . Perhaps the main reason for the lack of such data is due to the early adoption of standard field techniques for the capture of stream samples that mirrored protocols used for capture of standing water. Furthermore, because of the association of bacteria and suspended sediment, standardized techniques for the collection of water quality data have been published for sampling streams and/or flowing water sites according to rigorous protocols for suspended sediment analysis [21] -[23] . The purpose of such standards is to: 1) provide methods to minimize bias within and between data sets; 2) provide reproducible data; 3) encourage consistent field sampling methodology for data comparisons; and 4) provide citable documentation of data collection protocols to reinforce all the above. The philosophy of standardized water capture is to achieve consistency from one sampling event/date to the next.
1.2. Historical Perspective on Sampling for Indicator Bacteria
In 1965, early microbial ecologist A.G. Rodina stated “the selection of sampling locations in water mass is an important task which must be performed with all possible care” [24] . Sampling flowing water for microbial presence requires sufficient sample size to adequately represent the material sampled, yet reduced volume to allow both efficient transport and preservation of sample quality. Bacterial presence in flowing waters is measured from selected points or sampling sites along the drainage. Higher numbers of bacteria tend to be associated with warm months, proximity to point source contributors (e.g., wastewater plants, agricultural operations, storm drains, etc.), and storm runoff [5] [25] . Depth- and width-integrated approaches to sampling moving waters are difficult, especially in streams with low depth and velocities [23] . Therefore, in smaller streams, catch- or grab- samples can be obtained by hand from the main channel at mid-depth [24] with the understanding that any departure from standardized technique should be based on sound field judgment and be quality-assured, and documented.
1.3. The Dilemma of Obtaining a “Representative” Sample
Scientists and technicians within the public and private water monitoring community grapple with the fact that streams and other moving water bodies present a dynamic medium from which to obtain a representative sample for bacterial assessment for either surveillance or regulatory purpose. In spite of published standard methodologies, stream hydrography can change between sampling periods, reducing sampling precision as vertical (depth) and horizontal position of the main stream channel varies. Accordingly, in an attempt to obtain representative bacterial counts for a sample site, environmental regulatory agencies suggest taking multiple samples during a prescribed time to obtain an average or geometric mean for the location.
So the question remains: As care is taken to sample for bacteria in a stream, the question of variability in bacterial concentration for any one sampling site, regardless of depth or lateral location, remains to be examined. To our knowledge, these data are the first to address this question within a moving freshwater body.
1.4. The Experimental Question
This study examines the differences between counts of E. coli indicator bacteria from samples obtained at different points within the vertical and lateral water columns of four freshwater streams in the upper Appomattox River watershed in the Piedmont region of central Virginia. It represents one of the first reports in the literature examining bacterial numbers within the potentially stratified water column of a moving water body. The results provide critical information and implications relative to standard stream sampling protocols for water quality technicians and scientific professionals within the field of water quality assessment, specifically, and for environmental microbiologists, in general.
1.5. Objectives
The primary objective of this study was to determine whether or not the standard bacterial sampling protocol utilized and approved by most environmental agencies (e.g., Virginia Department of Environmental Quality [VA-DEQ] or United States Environmental Protection Agency [USEPA]) provided a representative count of bacterial concentrations within peripheral locations (both vertical and lateral) at a designated stream collection site. A secondary objective was to determine if numbers of E. coli indicator bacteria varied within the vertical and lateral water column of a moving freshwater stream at a given time and place.
2. Material and Methods
2.1. Sampling Locations
Weekly water samples were taken during the spring and fall months for two successive years from four different stream sites (labeled as APP1, APP2, BUF1, and VAU14) in the upper Appomattox River watershed within Buckingham, Cumberland, and Prince Edward counties in central Virginia (n = 936 samples). The sampling sites are located within a 15-mile radius of the town of Farmville, Virginia (Figure 1). The streams reveal aver-
Figure 1. Map of water sampling sites within central Virginia, USA.
age depths between 0.3 m and 1.5 m at baseline flow and 10-year average discharges of ≤245 cfs. Geographically, all sample sites are classified as second order streams with the exception of the APP2 location, which is on a third order section of stream. Stream site locations with GPS coordinates are presented in Table 1. For all sampling locations, the predominant land use of the sub-watersheds is agricultural (beef, dairy, poultry, row crops, and hay fields) interspersed with mixed hardwood forest and low housing density.
All four stream sampling sites are located at highway bridge access points. At each stream site for each sample date, three mid-column water samples were obtained within the deepest channel of the stream from the bridge perch that spanned the stream (identified as “bridge” samples). In addition, three replicate water sample sets were obtained within a 50 m distance upstream of the bridge from randomly chosen locations both within or peripheral to the main channel and assessed for bacterial numbers from each of three depth locations (high, mid, and low water column).
2.2. Sampling Protocol
At each site for each sample date, three separate samples were obtained from a bridge perch via a weighted, surface sterilized, tethered stainless steel cylinder (similar to a dip sampler) holding a sterile 120 mL polystyrene clinical sample container (NCS Diagnostics, Etobicoke, Ont.) and lowered to mid-column in the deepest channel of the stream. In addition, three separate sets of stratified water samples were taken within 50 m upstream of the bridge via hand grab method including: upstream-high column (water sample taken at 2.5 cm below stream surface), upstream-mid column (water sample taken at mid-column in stream), and upstream-low column (water sample taken at 2.5 cm above stream bottom). The stratified water samples were obtained facing upstream at three random but separate locations in succession, moving upstream of the bridge perch (Figure 2). Care was taken to obtain stratified samples at the specified depth using a meter stick placed in the water as a guide (Figure 3). In each case, a closed sample bottle was submerged to the measured depth, pointed into the current and
Table 1. Water sampling sites with coordinates within the upper Appomattox River watershed in central Virginia.
Figure 2. Photo showing three random sampling locations within 50 m upstream of representative bridge perch (APP2 site) showing three random upstream locations where separate sets of stratified water samples were taken (A = high column, B = mid column, C = low column) within the stream water column. Three mid-column samples from the main stream channel were also obtained from the bridge perch.
Figure 3. Capture of stratified water samples showing the use of a meter stick while opening sample bottle facing upstream and towards the water current.
opened briefly to allow entry of water to ensure clean capture of water. All water samples were captured in separate sterile sample containers, capped tightly, and immediately submerged in ice for transport back to the laboratory to be assessed for bacterial concentration using the Colilert® defined substrates (IDEXX, Westbrook, ME) medium for enumeration of total coliforms and E. coli.
2.3. Laboratory Analysis
All samples were aseptically processed within the laboratory for bacterial analysis within three hours of sampling in accordance with APHA guidelines [26] . Coliform bacteria and E. coli were assayed using Colilert Quanti-tray 2000® envelopes (IDEXX, Westbrook, ME) with a 25% sample dilution (25 mL test sample: 75 mL sterile phosphate buffered dilution water) according to manufacturer’s directions, and incubated at 44.5˚C ± 0.2˚C for 24 ± 2 hrs. For this analysis, only counts of E. coli were used for comparison. One sample blank was pro- cessed with each sampling date cohort to ensure sterility of dilution water. After incubating, wells showing chromogenicity (yellow) and fluorescing under long UV (365 nm) illumination were counted positive for E. coli. All enumerations were performed using a most probable number (MPN)-based system with a quantification range between <1 and 9676 CFU per 100 mL when using a 25 mL sample dilution. Weekly aliquots of several randomly selected Quanti-tray wells showing both chromogenic and fluorogenic responses were aseptically transferred using sterile, disposable 3 cc medical syringes to confirm lactose fermentation into lactose broth (EC broth, BDL, Sparks, MD) and to check IMViC responses (broth/agar media-DIFCO Labs, Detroit, MI). Confirmation of E. coli is enhanced by IMViC testing accompanied by the elevated (44.5˚C) incubation temperature [27] . All wells presumed to contain E. coli were confirmed as positive. Recent data in this and other labs [28] [29] indicate confirmatory tests are not needed when using the Colilert medium on freshwater samples. Colilert test results were recorded as number of total coliform bacteria (not reported in this study) and as number of E. coli per 100 mL for all data points.
2.4. Statistical Analysis
Counts of bacterial “colony forming units” per 100 mL (CFU/dL) of ambient water sample are generally believed to be variable among both field duplicate samples and sample or laboratory quality assurance (QA) duplicates (“splits”). The experimental design of our study violates a fundamental assumption of parametric analysis of variance (ANOVA)―the independence of measurements [30] . Samples that are collected within minutes of one another from various locations at the same field site are not independent; the samples are expected to be related, if not identical in a statistical sense. However, there are no published data that test this assumption involving bacterial assays, and most microbiologists and field technicians believe bacterial counts from such samples to be different (personal communications).
2.4.1. Experimental Treatments
The range of treatments for each sampling event in this study include: 1) Bridge-Mid samples (obtained from the bridge perch via weighted sampling device submerged to mid water column), and within-stream samples composed of 2) Upstream-High or near-surface channel samples (obtained by hand at 2.5 cm depth), 3) Upstream- Mid or mid-depth channel samples (obtained by hand at mid-depth in water column), and 4) Upstream-Low or near-bottom channel samples (obtained by hand at 2.5 cm from bottom substrate).
2.4.2. Experimental Blocks
Three replicate samples were collected for each of the four treatments at one sampling site (n = 12) for each sampling event or date. This statistical analysis groups data by site and sampling event to ensure comparison of homogenous data, and for logical reasons, disregards data comparisons between sites or sampling events.
2.5. Nested ANOVA
After the first year of data was recorded, variance components of bacterial counts were evaluated across 37 sam- pling events (site-date combinations), across four locations per event (bridge, high, mid, low), with three replicates per location (n = 444 samples) using a two-level nested parametric analysis of variance. All bacterial counts were log-transformed (X’ = log10 [X]) prior to the analysis. Because no zero bacterial counts were observed, the standard modified log-transformation (X’ = log10 [X+1]) was not required.
2.6. Friedman One-Way ANOVA
Subsequent to the results of the parametric analysis, we continued sampling with the same experimental design through the second year to increase the sample size at each site (i.e., to increase the number of blocks). A second phase analysis was then applied individually to each study site using the non-parametric Friedman one-way ANOVA based on ranks. The Friedman test is a non-parametric analysis of variance designed for making comparisons among several related samples, such as the results of different “treatments” (in this case, Bridge-Mid, Upstream-High, -Mid, and -Low samples) applied within “blocks” (date/site combination) or groups of experimental units that are similar in type. Its outcome is based on comparison of ranks of measurements within each complete block and is especially appropriate when comparisons between blocks are not of interest or are not possible [31] . A slight modification of the Friedman procedure, as described by Conover [31] , was applied to the replicated within-block measurements of the current study. Again, each sampling event (site/date combination) is considered as a block, and the various sampling depths and/or locations within a site are analogous to treatments.
All calculations were carried out on Excel spreadsheets using formulas from Conover [31] and standard Excel functions. In the case of tied measurements, each of the tied observations was assigned the mean rank of the tied group; the “RANK” (number, ref, order) function in Excel assigned the lowest rank to all members of tied groups.
3. Results
3.1. Pooled Data Comparisons
Pooled means (±std. dev.) of E. coli counts obtained from each of the treatment depths for each sampling location are presented in Figure 4. Although two of the sampling locations (APP1 and VAU14) contained comparatively lower bacterial concentrations, both the mean and variability of E. coli counts remained uniform across sample depths (p > 0.5). Likewise, the sampling sites that revealed higher concentrations of E. coli (APP2 and BUF1) showed no differences (p > 0.5) with respect to mean or variability across sample depths. The relatively high variability in bacterial counts from each sample site (Figure 4) was due to sampling dates occurring on/ near precipitation events as well as during dry periods which generally corresponded to higher and lower bacterial counts, respectively.
Figure 4. Pooled mean data. Mean ± std. dev. of pooled E. coli counts at different depths within the water column from each of the four water sampling sites revealed no differences (p > 0.5) either within the cleaner (APP 1 and VAU 14) or more contaminated (APP 2 and BUF 1) sampling sites. Note: “Bridge” = Bridge-Mid, “Hi” = Upstream-High, “Mid” = Upstream-Mid, and “Lo” = Upstream-Low sample sources.
3.2. Nested ANOVA of Early Data
The results of the early analysis showed that, across 37 sampling events at four sites on 18 separate dates for the first year of the study, variation in E. coli counts between site-date sampling events accounted for 91.22% of the total observed variation among bacterial counts (Table 2). For each individual site (APP1, APP2, BUF1, and VAU14) and sampling date, variation in E. coli counts among the four sample locations (replicates: 8.66%) was higher than the variation between the water column depths (0.12%).
The results of the nested ANOVA comparing differences between counts of E. coli at different sample depths (Bridge-Mid; Upstream-High, -Mid, -Low) were not significant (Table 2, p = 0.3921). Only the differences between site-date combinations (blocks) were significant (F36, 443 = 122.4865, p < 0.0001); however, this study is not examining differences between the sites or sampling dates, so this finding has no relevance to the purpose of our study.
These results provided an early indication that there were no significant differences in E. coli counts among the four sampling locations/depths at each site (i.e., no indication was observed that stratification of bacteria occurred within the water column). Two lingering doubts, however, led us to continue sampling and to carry out an additional statistical analysis; 1) we could not be certain that log-transformation of the bacterial counts met all the assumptions of the parametric ANOVA carried out in the first phase, and 2) perhaps more importantly, we could not be certain whether stratification at one or more site-date combinations may have been masked by the lack of stratification at other site-date combinations in the overall analysis. We believed that local variations in
(a) (b)
Table 2. Two-level nested ANOVA of bacterial counts obtained from the water columns at each of the four stream sampling sites.
hydrography might cause differences in stratification of bacterial presence within a particular stream site.
3.3. Friedman Analysis of Complete Data Set
Table 3 shows the results of the Friedman non-parametric analyses comparing bacterial counts obtained from the sampling depths for each of the four sampling locations for the entire data set. No differences were observed among bacterial counts within each block of data from the four sampling locations at each site. None of the Chi Square values of the analysis were significant, showing p values greater than 0.25 (Table 3). Accordingly, one data point (i.e., one sample) appears to be representative and suitable to describe the point distribution (bacterial counts) within any one block (date/site combination) for each of the four streams of this study.
3.4 Relationship of Pooled Means and Variances
Our data closely adhere to a Poisson distribution in that the variance values closely approximate the means. Figure 5 shows the relationship between log-transformed pooled means and variances for each of the 78 sampling events across the study. The data reveal a high, positive correlation (r = 0.8680) between mean and variance which relates to generally higher variability in E. coli counts when high concentrations of these bacteria are assessed and lower, more restricted variation when low E. coli concentrations are measured. In the context of water quality assessment, higher concentrations of E. coli can be associated with increased water flow in the drainage area that would occur during precipitation events and/or periods of increased surface run-off.
4. Discussion
Our data address both objectives for this study as they positively support our hypothesis that concentrations of the commonly used E. coli bacterial indicator are not stratified, but are homogenous within a generalized sampling area (of at least 50 m) and throughout the accompanying water column for any particular sampling date for small, freshwater streams. A representative sample for a bacterial indicator such as E. coli for each of the streams of this study can be provided either from a bridge perch sampling position obtained by mechanical sampler device or from an in-stream position obtained manually if care is taken to insure for clean capture of water (e.g., with no detritus, benthic material, or surface film). Concentrations of E. coli do not differ meaningfully within the adjacent upstream water column of a bridge access point nor do E. coli numbers differ when obtained by either hand grab or mechanical method for any individual sample date/time. Current field sampling methods stress the importance of obtaining “clean” samples of moving environmental waters taken at or near the
Table 3. Friedman one-way non-parametric ANOVA based on bacterial count rankings.
Figure 5. Correlation between pooled log-transformed means and variances for each of the individual site-date sampling events (n = 78). The computed correlation coefficient (r = 0.8680) reveals a strong, positive correlation between means and variation of E. coli counts per sampling site and date.
midpoint of the water column [32] . Although written mainly for physical and chemical analyses, the USGS National Field Manual [23] suggests several samples for bacterial assay should be taken within the cross-sectional area of a stream in order to be sufficiently representative and to lessen sample bias. Additionally, sampling protocols suggest water samples be taken from reaches upstream of bridges and other structures to avoid contaminating and/or toxic effects of surface runoff, etc. With the caveat of obtaining clean water samples near the main channel of the stream, these data reveal homogenous counts of indicator bacteria across multiple samples for a given sample date and homogeneity of counts either at or within 50 m upstream of a bridge perch.
To our knowledge, this is the first study to examine stratification of a bacterial water quality indicator within a small freshwater stream. These data reveal consistent homogeneity of E. coli concentrations at a variety of depths within the water column on any particular sampling date provided they are collected within or peripheral to the main channel. Whereas these data appear to contend the need for depth- and width-integrated approaches to sampling freshwaters proposed for decades [21] [26] [33] for bacterial samples, the authors of this report stress the importance of standards of sample quality (e.g., without disruption of stream bottom or inclusion of surface film) as reported in the literature [34] [35] .
Although not reported here, our data also reveal within-event consistency of counts with respect to season and precipitation effects. Bacterial counts during cold or warm seasons varied to the extent reported by others in a variety of aquatic environments [36] [37] , but the data did not reveal differences in count with respect to position in the water column for any specific sample date within any season. Precipitation events cause concomitant rises in counts of indicator bacteria by resuspension of bottom substrate particles [38] -[40] and increased runoff from adjacent land surfaces. Whenever precipitation occurred during or just prior to a sampling date for this study, a rise in E. coli indicators was observed across the cohort of samples for any one sample site, yet the variability of concentrations with respect to stream depth and location was insignificant.
It is imperative that environmental and public health agencies obtain reliable estimates of the bacteriologic quality of their nearby waterways to instill confidence in their constituent audiences and to insure for public safety. Our data suggest that a water sample for indicator bacterial assay taken with care on the same day/time from several different depths or from an adjacent or peripheral region along a moving water body is similar to one taken at the exact midpoint of the water column of the main channel. This finding allows some latitude with respect to choice of sampling locations for technicians, scientists, or other professionals involved in assessments of E. coli indicators of surface water quality for both urban as well as rural water monitoring. Our findings also suggest that rivers and streams showing high concentrations of indicator bacteria such as E. coli reveal such high counts homogenously throughout the cross-sectional area of the moving water body. Accordingly, the potential for human contact with fecally-derived pathogens (e.g., Salmonella, Campylobacter, Giardia, enteric viruses, etc.) for which these bacteria highly correlate/indicate may be equal throughout the water column.
Acknowledgements
This work was funded, in part, by grants from the Cook-Cole College of Arts and Sciences Undergraduate Research Support, Longwood University, and from the Virginia Department of Environmental Quality (VADEQ Contract numbers 14,309 and 14,733). Many thanks to W. Witschey for providing the map of our study sites. This work is dedicated to R. K. Riley and W. J. Vail, Emeritus Professors of Biology, Department of Biology, Frostburg State University.
NOTES
*Corresponding author.