Patrick M.
D'Aoust
a,
Robert
Delatolla
*a,
Alexandre
Poulain
b,
Galen
Guo
b,
Ru
Wang
a,
Colin
Rennie
a,
Liyu
Chen
a and
Frances R.
Pick
b
aDepartment of Civil Engineering, Faculty of Engineering, University of Ottawa, Room A108, 161 Louis Pasteur Private, Ottawa, K1N 6N5 Canada. E-mail: Robert.Delatolla@uottawa.ca; Tel: 1 613 562 5800 ext. 2677
bDepartment of Biology, Faculty of Science, University of Ottawa, Ottawa, K1N 9B4 Canada
First published on 25th May 2017
Stormwater retention ponds have become an integral component of stormwater management across the world. Under prolonged hypoxia, these ponds are capable of releasing large quantities of hydrogen sulfide (H2S) gas. In this study, water quality constituents and bacterial communities in sediment were analyzed in two stormwater retention ponds, RSP1 (reference pond) and RSP2 (problematic pond) over a period of two years, to identify the factors driving H2S production and understand the microbial community associated with H2S production in stormwater ponds. It was found that the background total sulfide concentrations were not statistically different between the two ponds during summer (RSP2: 0.012 ± 0.001 mg L-S−1; RSP1: 0.010 ± 0.001 mg L-S−1) and were statistically different during ice covered winter operation (RSP2: 6.375 ± 1.135 mg L-S−1; RSP1: 0.016 ± 0.009 mg L-S−1). The study showed a lack of correlation between total sulfide concentrations in RSP2 and soluble chemical oxygen demand, sulfate, soluble total phosphorus, total ammonia nitrogen, nitrate, nitrite and pH. However, DO concentrations demonstrated a strong negative correlation with total sulfides concentrations in RSP2 (p < 0.006, r = −0.58, n = 26), which confirmed DO as the critical water quality parameter linked to H2S production in stormwater ponds. Finally, it was found that seasonal change, ice covered versus non-ice covered operation and a comparison between a H2S emitting pond and non-emitting pond all did not promote a measurable proliferation of sulfate-reducing bacteria nor a community shift in the sulfate-reducing bacterial population. Hence, the study demonstrates that sulfide production is a result of increased ubiquitous SRB activity in stormwater retention ponds and the emission of H2S gas is not indicative of SRB proliferation or a population shift towards specific SRB taxa.
Water impactIce covered operation of stormwater retention ponds exacerbates hypoxic conditions in northern climates, which can lead to hydrogen sulfide (H2S) production. Water constituent and molecular microbial analyses show that H2S production was correlated to low dissolved oxygen concentrations and that sulfate-reducing bacterial populations were not observed to vary between a H2S emitting pond and non-emitting pond, during ice covered operation as comparted to non-ice covered operation or seasonally in stormwater ponds. |
The European Water Framework Directive (EU-WFD) and the Canadian Water Act provide guidelines for the design of stormwater retention ponds for all member-states/provinces,9 where every individual member-state/province also hold their own set of regulations and guidelines.9 As municipalities attempt to mitigate flooding risks, infrastructure damage, land washouts and negative quality impacts on the receiving waters, many have now adopted stormwater retention ponds as one of their main tools to mitigate the environmental impact of increased urbanization. The International Stormwater Best Management Practices (BMP) Database (ISBMPD), which collects and repertories government-submitted data and studies of stormwater facilities, has currently logged over 530 BMP studies, investigating over 16000 stormwater management facilities. These studies investigated performance and treatment efficiency of stormwater facilities; with approximately 57 of the facilities being stormwater retention ponds. An analysis of the 2014 ISBMPD summary report10 reveals that many stormwater retention ponds failed to meet their treatment objectives for total dissolved solids, and failed to reduce the loading of certain dissolved metals, such as nickel.11–13 In addition, climate change is projected to change global weather patterns,14,15 which will have potentially adverse consequences on the hydrological cycle and these facilities. Northern temperate areas, such as Southern Canada,16 the Northern United States17 and Northern Europe18 are expected to experience increased, more intense or more frequent precipitations19,20 and warmer daily minimum temperatures.16,21 In response, governments are striving to improve urban stormwater management planning, practices and policies. New strategies include increases to the number and size of stormwater retention ponds to accommodate the anticipated heightened precipitations.22 Although larger ponds will improve the retention capacity for large rain events, this large retention capacity may also impact the water quality of the ponds and increase the occurrence of hydrogen sulfide (H2S) production at these facilities due to the potential of lower dissolved oxygen conditions in larger designed retention basins.
H2S is a noxious and toxic gas produced by sulfate-reducing bacteria (SRB). SRBs are anaerobic microorganisms which utilize sulfate (SO42−) to obtain energy.23 A by-product of this reaction is the release of H2S. The occurrence of H2S gas in stormwater retention ponds is an indicator of sub-optimal facility design or operational problems as it is produced during periods of significant hypoxia.24 Although Makepeace25 recognized H2S production as a potential problem in stormwater retention ponds, the factors leading to H2S gas in these systems have not been clearly identified. There is presently a fundamental lack of knowledge and understanding of the processes and factors affecting the initiation and sustained production of H2S in stormwater retention ponds.
The aim of the study is to identify and quantify the key parameters influencing H2S production and total sulfide concentrations in stormwater retention ponds during various seasons of operation, including ice covered operation during winter months. In particular, water quality constituents and the microbial structure of pond sediment collected from two residential stormwater retention pond facilities in Ottawa, Canada were studied and compared across a period of a year. Water samples were collected in a manner permitting the analysis of spatial and depth variations throughout the facilities. The two ponds studies in this research include a problematic pond, Riverside South Pond #2 (RSP2), that has a history of H2S emission and a reference pond located in close proximity to RSP2, Riverside South Pond #1 (RSP1), that does not emit H2S.
![]() | ||
Fig. 1 Stormwater ponds, a) RSP1 with sampling locations and flow direction, b) RSP2 with sampling location and flow direction, c) bathymetry of RSP2 with depth in meters. |
Water samples were collected at these sites at 0.20 m and 1.50 m depths. Samples were collected using a Wildco 1520 C25 Kemmerer 2.2 L TT water sampler (Yulee, FL). A simple modification was performed by adding a 0.40 m long, 12.7 mm inner diameter piece of silicone tubing on the decanting valve of the sampler; the added tubing restricted air entrainment into the sample containers during the collection of water at depth in the pond. The minimization of air entrainment into the water samples reduced the effect of oxygen on the total sulfide concentration of the water samples and ultimately increased the precision of the total sulfide measurements. Total sulfide samples were immediately preserved on-site with the addition of zinc acetate and sodium hydroxide solutions, as per standard method 4500-S2 D.26
The following water quality variables were measured in accordance with standard methods26 and US EPA methods:27 i) total sulfides (SM 4500-S2-D), ii) total ammonia nitrogen (SM 4500-NH3 B), iii) sulfate (US EPA 375.4 US), iv) nitrate (SM 4500-NO3− B), v) nitrite (SM 4500-NO2− B), vi) soluble chemical oxygen demand (SM 5220 D), vii) soluble total phosphorus (SM 4500-P E) and viii) pH (SM 4500-H + B).
Molecular amplicon sequence datasets were analysed using the Quantitative Insights Into Microbial Ecology (QIIME) pipeline.28 The data was quality filtered using QIIME default parameters (quality score = 25, min length = 200, max length = 1000). Additional quality filtering and operational taxonomic unit (OTU) clustering was performed with the USEARCH V9.0,29 which utilizes the UCHIME algorithm to identify chimera sequences for removal against the Gold databases. De novo OTU picking with UCLUST was used to form the representative OTU dataset.29 Bacterial taxonomy was assigned using the RDP classification algorithm against the Greengenes databases. Subsequent community structure analyses were conducted using R with the phyloseq package.30 OTUs were analyzed in the context of relative abundance and UniFrac distances. Principal coordinate analysis (PCoA) using the Analyses of Phylogenetics and Evolution (ape) R statistical package31 and permutational multivariate analysis of variance using distance matrices (Adonis) was performed and using the vegan statistical package.32
ddPCR evaluation was conducted using a BIO-RAD QX200TM ddPCR system (Hercules, CA). Count data from the ddPCR was acquired using the Quantasoft software, developed by BIO-RAD (Hercules, CA). The expected amplicon size for SRB and methanogens were of 221 bp and 491 bp, respectively. Table 1 outlines the primers utilized for ddPCR. All primers were purchased from Integrated DNA Technologies (Coralville, IA).
Microbial population | Primer | Sequence |
---|---|---|
Sulfate-reducing bacteria | dsr1-F RT | 5′-ACS CAC TGG AAG CAC GGC GG-3′ |
dsr-R RT | 5′-GTG GMR CCG TGC AKR TTG G-3′ | |
Methanogens | mcrA R | 5′-CGT TCA TBG CGT AGT TVG GRT AGT-3′ |
mlas F | 5′ GGT GGT GTM GGD TTC ACM CAR TA-3′ |
High total sulfides concentrations are often a result of sulfate-reduction.34 The concentrations of total sulfides in this study are within the range found in other cold, deep, and/or strongly stratified aquatic systems. These includes a stormwater retention pond in Edmonton, Canada that experienced 1.4–3.6 mg L-S−1;24 Onondaga Lake NY, US that demonstrated a maximum of 56.23 mg L-S−1 (ref. 35) and the Torquay Canal, US that demonstrated ≥40.90 mg L-S−1.36 Additionally, as seen in Fig. 2, during summer between June 12th and 25th, 2015 an H2S production event lasting approximately two weeks was measured at the outlet of RSP2 (i.e. RSP2-4). This summer event showed a maximum concentration of total sulfides at the RSP2 outlet of 0.628 ± 0.007 mg L-S−1, while concentrations during the same period in RSP1 were measured at 0.025 ± 0.002 mg L-S−1. Hence, the concentration of total sulfides in RSP2 was significantly greater than the concentrations found in RSP1 during this event. The production of sulfides found at the RSP2 outlet during the summer production event is similar to other reported cases of sulfide production in warm lakes, such as in Lake Brooker FL, US where concentrations were recorded at 0.176 ± 0.069 mg L-S−1.37
The background daily total sulfide concentrations at all sampling locations in the two ponds (RSP2-1 to 4: 0.012 ± 0.001 mg L-S−1; RSP1-1: 0.010 ± 0.001 mg L-S−1) were not statistically different to each other during periods with low total sulfides production. Further, the calculated average background total sulfides concentrations at depths of 0.20 m and 1.50 m along with the maximum concentrations of total sulfides measured in RSP2 were shown to not be statistically different at the two sampling depths and were shown not be different spatially throughout the year of operation. As such, both RSP1 and RSP2 were not strongly chemically stratified throughout the year, with the exception of during ice cover (December 15th, 2014 to April 8th, 2015) and during the sulfide production event at the outlet of RSP2 (June 12th to June 25th 2015).
Although the daily, average and maximum total sulfide concentrations did not show differences spatially across the pond or with depth, two differences were observed between RSP2-4 and other sampling locations of RSP2: first, the summer sulfide production event was unique to the outlet of RSP2 and hence location RSP2-4 and secondly, there was a lack of statistically validated stratification of H2S with depth during the ice covered event at RSP2-4. The sampling location RSP2-4 is located in close proximity to the outlet of the pond, at the location with the greatest depth and was qualitatively observed in the field to accumulate the greatest quantity of sediment as compared to other locations in the pond. Based on the summer sulfide production event isolated to RSP2-4 and the saturated water column with hydrogen sulfide during ice covered conditions, it is likely that the deepest portion of RSP2 with the greatest accumulated quantity of sediment initiated sulfide production in the pond.
Further, sediment samples collected from the two ponds demonstrated that iron concentrations were slightly lower in RSP1 compared to RSP2 throughout the study period. Iron concentrations of 20.06 ± 0.33 and 23.17 ± 0.33 mg g−1 dry sediment were measured for RSP1 and RSP2 respectively. Hence, the slightly lower iron concentrations in RSP1 compared to RSP2 throughout the study period indicate that iron sequestering in the ponds investigated in this study is likely not a dominant factor affecting sulfide concentrations in the pond water columns. This finding is further supported by the fact that both facilities were dug from the same native clay.
There was no observed lag period between the onset of hypoxic conditions and a significant increase in total sulfides at warmer temperatures between June 12-25, 2015 in RSP2-4 or under ice cover at all locations in RSP2 or RSP1. Low DO concentrations (<2.0 mg L−1) at depth (1.50 m) in RSP2 were first observed at RSP2-4 on January 7th, 2015, followed by RSP2-2 and RSP2-3 on January 9th, 2015, and finally at RSP2-1 on January 21st, 2015. Low DO concentrations at RSP1-1 were only first observed approximately a month later, on February 12th, 2015. DO concentrations <2.0 mg L−1 near the surface (0.20 m) were first observed at all locations (in RSP2 on February 12th, 2015). DO concentrations <2.0 mg L−1 in RSP1-1 near the surface occurred, again at a later date compared to RSP2, on March 3rd, 2015. Additionally, there was periodic stratification of DO concentrations during the ice covered period at all locations with stratification occurring at a later date at RSP1-1, as shown in Fig. 2, starting at the end of December 2015 and continuing during January and February of 2016.
During the ice covered period, NH3/NH4+ concentrations increased, peaking at all locations on March 20th to 30th 2015. Similarly to sulfide, there was a slow progression of high ammonium/ammonia concentrations at depth which progressed from the bottom to 0.20 m below the surface. Initially, concentrations were determined to be statistically different at depth versus near the surface (apparent stratification), but towards the end of the ice covered period (March 2015), concentrations were similar at all locations and at all depths within RSP2 (1.59 ± 0.52 mg L−1 NH3–N). At the same time (March 2015), concentrations in RSP1-1 were slightly lower at (1.23 ± 0.48 mg L−1 NH3–N).
During summer, total ammonia nitrogen concentrations were low and similar at all locations. The average concentrations measured in RSP1 and RSP2 at 1.50 m of depth were 0.32 ± 0.27 mg L−1 NH3–N and 0.25 ± 0.28 mg L−1 NH3–N, respectively. It is hypothesized that low ammonium/ammonia concentrations measured during the summer period are due to the uptake of ammonium by algae and aquatic plants.
During the ice covered period, the average pH value in RSP2 decreased slightly (low of 7.27 ± 0.28), while values in RSP1-1 remained stable (7.48 ± 0.20). The decrease in pH is possibly the result of greater anaerobic activity of the benthic sediment of RSP2,46 and coincides with a decrease in DO and production of sulfides (Fig. 2 and 4). During summer periods, the average pH values in RSP2 were 7.90 ± 0.41, while values in RSP1-1 were 7.35 ± 0.44. It is hypothesized that the higher pH values in RSP2 were the result of higher primary production rates in RSP2 than at RSP1-1,47 due to the higher uptake of carbon dioxide (and removal of carbonic acid).
Water quality constituent | Averaged across RSP2 | RSP2-4 | RSP1-1 |
---|---|---|---|
Sulfate (mg SO42− L−1) | 50.1 ± 10.9 | 49.51 ± 12.6 | 46.5 ± 8.5 |
Soluble chemical oxygen demand (mg L−1) | 20.2 ± 12.4 | 21.2 ± 14.6 | 16.6 ± 8.0 |
Soluble total phosphorus (mg-P L−1) | 0.13 ± 0.13 | 0.15 ± 0.15 | 0.09 ± 0.08 |
Nitrate (mg-N L−1) | 1.04 ± 0.20 | 1.04 ± 0.23 | 0.92 ± 0.38 |
Nitrite (mg-N L−1) | <0.012 | <0.012 | <0.012 |
Sulfate, soluble chemical oxygen demand, nitrate and soluble total phosphorous concentrations were all stable spatially and at various depths in RSP2 throughout the entire study period. The measured sulfate concentrations indicate that sufficient sulfate was present for SRB activity.48 The soluble chemical oxygen demand concentrations are also within the required ranges for SRB activity.48 Soluble nitrate concentrations were stable spatially and at various depths throughout the entire study period, with average values of 0.92 ± 0.38 mg N L−1 in RSP1-1 and 1.04 ± 0.20 mg N L−1 in RSP2. Soluble total phosphorus concentrations were indicative of non-limited phosphorus conditions for microbial activity and hence SRB activity. Nitrite concentrations were below the practical quantification limit (PQL) of 0.012 mg N L−1 for the majority of the tested samples throughout the study, with all quantifiable concentrations being measured below 0.090 mg N L−1.
The percent abundance of the SRB population of the outlet sediments of RSP1-1 and RSP2-4 were not shown to differ significantly (p = 0.78) (Fig. 6), with the percent abundance of the microbial community being approximately 5.01 ± 0.79 and 6.22 ± 2.11 at RSP2-4 and RSP1-1 respectively. The top 10 dominant SRBs and their relative abundance identified in the outlet sediment of both ponds were also not significantly different (Table 3). Hence, these findings indicate that periodic hypoxia and hydrogen sulfide production events, as seen at RSP2-4, do not cause statistical distinction in the sediment microbial structure and the SRB population of stormwater ponds.
RSP1-1 | RSP2-4 | |
---|---|---|
Organism (genus) | Percent organisms (%) | Percent organisms (%) |
Family Desulfobulbaceae, unclassified genus | 2.39 ± 1.58 | 1.98 ± 0.49 |
Desulfococcus | 1.45 ± 0.92 | 1.14 ± 0.26 |
Family Desulfobacteraceae, unclassified genus | 1.25 ± 0.65 | 1.04 ± 0.26 |
Geobacter | 0.38 ± 0.07 | 0.29 ± 0.11 |
Desulfobulbus | 0.16 ± 0.12 | 0.17 ± 0.02 |
Desulfomonile | 0.15 ± 0.08 | 0.12 ± 0.10 |
Synthrophobacter | 0.11 ± 0.06 | 0.06 ± 0.01 |
Family Desulfuromonadales, unclassified genus | 0.07 ± 0.03 | 0.05 ± 0.03 |
Desulfobacca | 0.17 ± 0.23 | 0.05 ± 0.04 |
Desulfomicrobium | 0.04 ± 0.02 | 0.04 ± 0.01 |
The ddPCR counts of the SRB and methanogens in the sediment of RSP1-1 and RSP2-4 normalized per gram of sediment also showed no statistical differences in the quantity of SRB or methanogenic populations (Fig. 7). SRB counts were shown to be higher than methanogen bacterial counts, at all locations, regardless of season or temperature (Fig. 8). Both SRB and methanogenic populations showed no statistical correlation to season of operation, ice covered versus non ice covered conditions or to periods of elevated and background H2S concentrations at RSP2-1, RSP2-2, RSP2-3 or RSP2-4. Furthermore, no statistically significant correlations were found to exist between SRB counts and DO, total sulfide concentrations or temperature at RSP2-1, RSP2-2, RSP2-3 or RSP2-4. The lack of distinction between the microbial communities of the two ponds at various operational seasons and the lack of measured population shift or change in SRB count numbers during H2S production events indicates that hydrogen sulfide production in stormwater ponds results from a simple increase in the activity of a ubiquitous SRB sediment population that is initiated by hypoxic conditions.
![]() | ||
Fig. 7 Sulfate-reducing and methanogenic bacterial copies per g sediment in benthic sediment at RSP2-4 and RSP1-1. |
![]() | ||
Fig. 8 Total sulfides, SRB and methanogen population counts, water temperature and DO at RSP1-1, RSP2-1, RSP2-2, RSP2-3 and RSP2-4. |
During H2S production events it was observed that the problematic pond simply initiated an earlier and more significant decrease in DO concentration in the pond compared to the reference pond, which subsequently led to an earlier onset of H2S production in this pond. The concentration of total sulfides was shown to increase at shallower depths in the water column of the problematic pond and reach elevated total sulfide concentrations across a period of approximately 14 days. Furthermore, the deepest portion of the problematic pond, which corresponded to locations with the greatest quantity of accumulated sediment, were observed to show the highest propensity for the production of H2S.
Microbial structural analyses of the pond sediment shows that the microbial communities and in particular the SRB populations of the problematic and reference stormwater ponds are not statistically distinct. Furthermore, the microbial community of the two ponds did not undergo a shift and the SRB counts did not demonstrate a statistically significant change during varying environmental conditions of operation, or during H2S production. The lack of distinction between the microbial communities of the two ponds and the lack of measured population shift or change in SRB count numbers during H2S production events indicates that these events are a results of an increase in the activity of a ubiquitous SRB sediment population during the on-set of hypoxic conditions. These findings are supported by the rapid increase in total sulfide concentrations at depth during hypoxia at depth.
This journal is © The Royal Society of Chemistry 2017 |