Next Article in Journal
Assessing the Hydrologic Performance of a Green Roof Retrofitting Scenario for a Small Urban Catchment
Previous Article in Journal
Seasonal Variation of Drainage System in the Lower Ablation Area of a Monsoonal Temperate Debris-Covered Glacier in Mt. Gongga, South-Eastern Tibet
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In Situ Water Quality Improvement Mechanism (Nitrogen Removal) by Water-Lifting Aerators in a Drinking Water Reservoir

1
School of Energy and Environment, Zhongyuan University of Technology, No. 41 Zhongyuan Road (M), Zhengzhou 450007, China
2
Key Laboratory of Northwest Water Resource, Environment and Ecology, MOE, Xi’an University of Architecture and Technology, Xi’an 710055, China
3
Shaanxi Key Laboratory of Environmental Engineering, Xi’an University of Architecture and Technology, Xi’an 710055, China
*
Author to whom correspondence should be addressed.
Water 2018, 10(8), 1051; https://doi.org/10.3390/w10081051
Submission received: 18 July 2018 / Revised: 2 August 2018 / Accepted: 6 August 2018 / Published: 8 August 2018
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
A field scale experiment was performed to explore the nitrogen removal performance of the water and surface sediment in a deep canyon-shaped drinking water reservoir by operating WLAs (water-lifting aerators). Nitrogen removal performance was achieved by increasing the densities and N-removal genes (nirK and nirS) of indigenous aerobic denitrifiers. After the operation of WLAs, the total nitrogen removal rate reached 29.1 ± 0.8% in the enhanced area. Ammonia and nitrate concentrations were reduced by 72.5 ± 2.5% and 40.5 ± 2.1%, respectively. No nitrite accumulation was observed. Biolog results showed improvement of carbon metabolism and carbon source utilization of microbes in the enhanced area. Miseq high-throughput sequencing indicated that the denitrifying bacteria percentage was also higher in the enhanced area than that in the control area. Microbial communities had changed between the enhanced and control areas. Thus, nitrogen removal through enhanced indigenous aerobic denitrifiers by the operation of WLAs was feasible and successful at the field scale.

1. Introduction

Nowadays, more and more reservoirs are being used as water sources by large cities. Excessive nitrogen in these reservoirs results in poor water quality with a severe impact on aquatic ecosystems, and it has given rise to health considerations in humans and in livestock [1,2]. Thus, nitrogen removal from drinking water reservoirs is a critical problem.
Physical, chemical, and biological technologies have been used for nitrogen removal, however, physical and chemical methods incur high costs and are unable to remove nitrogen completely. Bioremediation (biological technology) refers to the use of microorganisms to eliminate or reduce the concentration of hazardous wastes at a contaminated site. Bioremediation has the advantage of no second-pollution or residues [3]. Many researchers have achieved satisfying results for bioremediation in water quality improvement using N-functional bacteria for phosphate, nitrogen, and chemical oxygen demands [4,5,6]. However, bioremediation through bioaugmentation may result in environmental safety problems by changing the composition of indigenous microbial communities [3]. Therefore, N-functional bacteria are rarely reported for use in purifying surface water ecosystems, especially drinking water reservoirs.
In the traditional opinion, bioremediation for nitrogen removal involves autotrophic nitrification and heterotrophic denitrification. Nitrification is achieved under aerobic conditions, whereas denitrification requires anaerobic and anoxic conditions through a sequence of intermediates (nitrate, nitrite, nitric oxide, and nitrous oxide), resulting in nitrogen gas formation [7,8]. Due to inhibition by dissolved oxygen, the traditional process of biological nitrogen removal can seldom be used in reservoirs. However, aerobic denitrifying bacteria that simultaneously utilize oxygen and nitrate as electron acceptors were reported in the early 1980s and provided novel insights in bioremediation for drinking water reservoirs [7,9]. Compared with traditional denitrification, aerobic denitrification has attractive advantages: nitrification and denitrification can exist in one system and denitrification can supply sufficient alkalinity to partially balance the acidity of nitrification [10,11]. The isolated aerobic denitrifying bacteria with satisfying nitrogen removal characteristics included Microvirgula aerodenitrificans, Bacillus methylotrophicus, Pseudomonas stutzeri strains T1, YZN-001, and P. stutzeri ZF31 [2,12,13,14,15]. Moreover, aerobic denitrification also occurred in natural systems, for example, in the Wadden Sea sediment [16] and horizontal subsurface-flow constructed wetlands [17].
Water-lifting aerators (WLAs) which were installed in reservoirs for water quality improvement has been used widely in recent times [18,19,20]. This gives us an opportunity to explore aerobic denitrification in situ. Based on our previous study on a pilot scale [2], an aerobic denitrification biofilm performed well with oxygen supplied by WLAs, and the nitrogen removal performance could meet the requirements of Class 3, based on the Chinese Surface Water Environmental Quality Standard (GB3838-2002). During the WLAs operation in an enclosure system in the Zhoucun Reservoir [8], nitrate in the enhanced system decreased from 1.71 ± 0.01 to 0.8 ± 0.06 mg/L, and the control system showed little removal. The total nitrogen (TN) removal rate of the enhanced system reached 38.33 ± 1.50%. These data suggest that aerobic denitrification is an effective way to decrease the endogenous nitrogen pollution in drinking water reservoirs. Of course, further studies are required to understand the nitrogen removal mechanism [16,17] and water quality improvement mechanism by the WLAs [8].
In this study, we performed a field study in a canyon-shaped drinking water reservoir named Jinpen Reservoir. In the main reservoir (with a surface area of nearly 2 km2), eight WLAs were running during our research. Many studies have shown that under anaerobic conditions, pollutants like N, P, and Fe would be released into the overlying water from the sediment [8,18,19,20]. However, the main function of WLAs is to destroy the stratification and oxygenate the hypolimnion when it is anaerobic. Thus, operation of WLAs provides an aerobic environment that can inhibit sediment release. Simultaneously, the operation can also enhance the indigenous aerobic denitrifiers; as observed, the densities of aerobic denitrifiers and the abundance of denitrification functional genes were both increased by one magnitude compared with that before operation of WLAs. However, in the control area (2 km far away from the WLAs), no obvious increase in aerobic denitrifier densities was observed. Denitrification functional genes (nirS and nirK) were also quantified using real-time polymerase chain reaction (PCR) both in the enhanced and control areas during the operation of WLAs. Biolog technology was used to investigate the microbial metabolic activity and the carbon source utilization. Microbial community variation was also explored to determine which genera of bacterioplankton were dominant in both areas. In the end, we also tried to (1) explore the nitrogen removal performance of the water and surface sediment system in the period of WLAs operation, (2) investigate the microbial metabolic activity and carbon removal efficiency, (3) study the changes in aerobic denitrification bacteria and denitrification functional genes (nirS and nirK), and (4) determine the diversity of the bacterial community and the relationship between the microbial community structure and environment driving factors, especially the bacteria involved in nitrogen cycling.

2. Materials and Methods

2.1. Sampling Sites and Field Work

As shown in Figure S1, a canyon-shaped deep drinking water reservoir named Jinpen Reservoir (at latitudes from 34°13′ N to 34°42′ N and longitudes from 107°43′ E to 108°24′ E) was our research site, and it supplies 80% drinking water to Xi’an City, northwest China. Its main reservoir length is 3.5 km. The total capacity of the reservoir is 2.0 × 108 m3 and its effective capacity is 1.8 × 108 m3. The main function of the Jinpen Reservoir is urban water supply. Agricultural irrigation, power generation, and flood control are its accessory functions. The reservoir had 8 WLAs installed for water quality improvement. The enhanced area was nearly 50 m away from the WLAs, whereas the control area was 2 km away. During the operation of WLAs (15 October to 2 November in 2015), water samples were collected every two days. Every water sample (0.5, 45, and 90 m) was collected in triplicate. The surface sediment samples in both the enhanced and control areas were collected every five days, in triplicate. Water samples from the enhanced and control areas were stored into pre-cleaned high-density polyethylene bottles. The samples were immediately cooled and were stored at 4 °C until analysis.

2.2. Biolog ECO Micro Plate (Carbon Metabolism Characteristics and Diversity of Microbes)

Biolog is an advanced method of investigating the carbon metabolism, community structure, and microbial diversity in environments. Average well color development (AWCD) of the average microbe activity was used to determine the utilization of carbon sources and metabolism characteristics. Absorbance values at wavelengths of 590 nm and 750 nm per hole were subtracted with the absorbance values of the control groups, respectively, and then the absorbance value at 590 nm was subtracted with the value at 750 nm in the same group, respectively (C590–750). Hence, the absorbance value per hole could be obtained. It is notable that a number below 0.06 was recorded as 0 [8]. The formula can be determined as
AWCD590–750 = Σ (C590–750)/31.
where 31 represents the number of carbon source varieties in the Biolog ECO Micro Plate.
Carbon source utilization (non-dimensional) was calculated by the mean values of absorbance for every carbon source in one class. For example, amino acids = (CL-arginine + CL-asparagine + CL-phenylalanine + CL-serine + CL-threonine + CL-glutamic acid)/6.

2.3. Quantification of nirS and nirK Genes

Quantification PCR was used to explore the nitrogen removal gene (nirS and nirK) copy numbers in both the enhanced and control areas (7 times from 15 October to 2 November). DNA was extracted from every 2 L water sample using a Fast DNA kit. Primers for nirS quantification were cd3aF, 5′-GTSAACGTSAAGGARACSGG-3′ and R3cd, 5′-GASTTCGGRTGSGTCTTGA-3′. Primers for nirK quantification were F1aCu, 5′-ATYGGCGGVCAYGGCGA-3′ and R3Cu, 5′-GCCTCGATCAGRTTRTGGTT-3′ [8]. Real-time PCR (qPCR) was performed on the ABI 7500 real-time system (Life Technologies, Carlsbad, CA, USA). SYBR Premix Ex TaqTM (Takara Bio, Dalian, China) was used to perform the amplification reaction. The qPCR amplification of nirS (cd3af/R3cd) and nirK (F1aCu/R3Cu) was performed as described previously [8].

2.4. Number Counts of Indigenous Aerobic Denitrifiers

The numbers of indigenous aerobic denitrifiers in both the enhanced and control areas were determined by plate counts and were tested via gradient dilution. The gradient dilutions were prepared as follows: 10−3, 10−4, 10−5, and 10−6, and 0.2 mL of the diluents was streaked onto a solid screening medium (including (g/L): CH3COONa, 0.10; NaNO3, 0.02; K2HPO4·3H2O, 0.02; CaCl2, 0.01; MgCl2·6H2O, 0.01; and agar, 20; pH 7.2) in triplicate, and incubated at 30 °C for five days. Single colonies were then harvested and counted [8].

2.5. Illumina Miseq Sequencing

The Illumina Miseq Sequencing platform was used to explore the microbial community composition. The DNA extracted from water samples was amplified by PCR using primers 27F (5′-AGAGTTTGATCCTGGCTCAG-3′) and 338R (5′-TGCTGCCTCCCGTAGGAGT-3′) targeting the V2 region of bacterial 16S rRNA genes. The PCR products were then sequenced using an Illumina Miseq Sequencing platform according to the standard protocols specified by the Shanghai Majorbio Bio-pharm Technology Co., Ltd., China. Low quality sequences were removed.

2.6. Physical Analysis

Temperature, DO (dissolved oxygen), pH, Chl-a (chlorophyll-a), turbidity, and ORP (oxidation-reduction potential) were determined in situ every 2 m using a Hach multi-probe water quality analyzer (Hydrolab DS5, Loveland, CO, USA).

2.7. Chemical Analysis in Water

The concentrations of TN (total nitrogen), nitrate, ammonia, and nitrite were determined using a SEAL AA3 HR Auto Analyzer (SEAL, Hamburger, Germany). TOC (total organic carbon) concentrations were measured with a TOC-L total organic carbon analyzer (Shimadzu, Kyoto, Japan).

2.8. Chemical Analysis in Sediments

Sediment samples were air-dried and ground before analysis. Sediment total nitrogen (TN) was determined by Zn-Cd reduction spectrophotometry after digestion with potassium peroxodisulfate and microwaving (121 °C, 1 h). Sediment total phosphorus (TP) was detected using the Mo-Sb colorimetric method. Sediment organic carbon (OC) was determined by the potassium dichromate volumetric method (Mack and D’Antonio, 2012).

2.9. Data Analysis

The fundamental analysis, abundance-based coverage estimators (ACE), Chao richness estimator, coverage, and the Shannon and Simpson index were calculated by MOTHUR. The gradient length of the longest axis explored by detrended correspondence analyses (DCA) was less than 3 SD (standard deviation) units. This indicates that the species exhibited linear responses to environmental gradients and therefore, further analyses were conducted using redundancy analysis (RDA) ordination response models [8].

3. Results and Discussion

3.1. Spatial and Temporal Variation of Environmental Factors in the Enhanced Area

The operation of WLAs decreased the temperature (T) of the surface water and increased that of the bottom water in the enhanced area. As shown in Figure 1, T of the surface water had decreased from 17.1 ± 0.21 °C (15 October) to 14.8 ± 0.12 °C (2 November), and T of the bottom water had increased from 10.1 ± 0.05 °C (15 October) to 14.3 ± 0.09 °C (2 November). T of the vertical water column was homogeneous. The increase in T of the middle and bottom water might help increase the activity of indigenous microbes. DO of the water column was increased, especially that of the bottom water (from 0 mg/L (15 October) to 9.8 ± 0.19 mg/L (2 November)). Thus, on 2 November, the T and DO of the enhanced area were maintained at a nearly homogeneous state. The pH and ORP (oxidation-reduction potential) of the water column increased gradually during the operation. The ORP of the bottom water especially increased to 504 mV from 25 mV. Somehow, the turbidity of the water column showed a decreasing trend from 22.2 ± 0.2 NTU (maximum value on 15 October) to 8.1 ± 0.1 NTU (maximum value on 2 November). CHl-a varied from 2.65 ± 0.18 μg/L to 1.09 ± 0.1 μg/L, and the CHl-a concentration in the water column was also homogeneous.
Compared with the enhanced area, there was little spatial and temporal variation of environmental factors in the control area. As seen in Figure S2, the T, DO, pH, turbidity, ORP, and CHl-a showed little change in the water column in the control area from 15 October to 2 November.
A non-operation year (2011) for WLAs was proposed for some comparison with 2015. As shown in Figure S3, on nearly the same date in 2011, the stratification of Jinpen Reservoir was stable, and the DO concentration in the bottom water was maintained at 0 mg/L until 31 October in 2011. From the results of the control area (2015) and those obtained in 2011, we can conclude that the mixing and oxygenation functions of WLAs were very efficient. Through the operation of WLAs, the DO concentrations of the reservoir—especially in the bottom water—were greatly improved.

3.2. Nitrogen Removal in the Enhanced Area

The changes in TN, nitrate, ammonia, and nitrite at the enhanced and control areas are shown in Figure 2. In the enhanced area, surface TN concentrations decreased from 1.72 ± 0.06 mg/L (15 October) to 1.21 ± 0.05 mg/L (2 November), and the removal rate (see Figure S4) reached 29.1 ± 0.8%. Compared with the control area, the surface TN concentrations decreased from 1.68 ± 0.10 mg/L (15 October) to 1.55 ± 0.03 mg/L (2 November), and the removal rate was only 7.7 ± 1.8%. The TN in the middle and bottom waters showed the same trend as that in the surface water. The TN in the bottom water decreased from 2.36 ± 0.12 mg/L (15 October) to 1.20 ± 0.04 mg/L (2 November), and the removal rate reached as high as 47.5 ± 0.5%.
In the Jinpen Reservoir, nitrate is the main form of nitrogen, accounting for 83.5 ± 5.9%. The nitrate removal rates (shown in Figure S5) of the surface, middle, and bottom water in the enhanced area were 33.6 ± 2.9%, 43.5 ± 2.6%, and 44.4 ± 3.2%, respectively. However, compared with the enhanced area, the nitrate removal rate of the surface, middle, and bottom water in the control area was only 11.1 ± 0.9%, 13.5 ± 0.6%, and 10.4 ± 0.2%, respectively. During the operation, no nitrite accumulation was observed both in the enhanced and control areas (Figure S6). Nitrite in both areas was maintained at low concentrations, and the maximum concentration in the enhanced area reached 0.04 ± 0.005 mg/L, decreasing gradually from 15 October to 2 November, and the nitrite in the control area was maintained at 0.03 ± 0.004 mg/L. As seen in Figure S7, the ammonia concentration in all three water layers showed a decreasing trend, especially in the bottom water, decreasing from 0.48 ± 0.06 mg/L to 0.12 ± 0.02 mg/L, and its removal rate reaching 72.5 ± 2.5%. No ammonia accumulation was observed during the nitrate removal period in the enhanced area. Compared with the enhanced area, the ammonia removal rate in the control area was quite small, 19.2 ± 0.5%. In order to explore where the nitrogen had gone, the nitrogen content in the sediments was also determined.
During the operation of WLAs, the surface sediments in both the enhanced and control areas were determined on 15 October, 21 October, 28 October, and 2 November, and the results are shown in Table 1. The OC and TN concentrations in the enhanced area had all decreased from 3.42 ± 0.20% and 1522 ± 66 mg/kg to 2.66 ± 0.12% and 1242 ± 61 mg/kg, and their removal rates were 22.2 ± 6.5% and 18.4 ± 3.6%, respectively. In contrast, the OC and TN concentrations in the control area were slightly increased. The reason of this increase may be inflow sedimentation as described in our previous study [8].
Based on these results, the nitrogen content in both the water and the sediment of the enhanced area was significantly removed. Thus, the activity of indigenous aerobic denitrifiers may be enhanced in the enhanced area in situ. As is already known, aerobic denitrification bacteria have been found in natural environments such marine sediments, constructed wetlands, and reservoirs [21,22,23], and have been reported only by these three studies. In this study, we also observe the phenomenon of nitrogen removal in a drinking water reservoir in situ with nitrite and ammonia accumulation. We consider that the indigenous aerobic denitrifiers were activated by the in situ operation of WLAs. Deviating from our previous research [23], the present study was performed in the field with the in situ operation of eight WLAs, and not in an enclosure system. The results obtained in the study of the enhanced area were consistent with those described before. However, a lower nitrogen removal rate was also found in the control area, which was inconsistent with previous results [23]. In our study, ammonia and nitrite were simultaneously removed via nitrification and denitrification, and similar results were obtained in previous studies [17,23]. In order to clearly understand whether the activity of indigenous aerobic denitrifiers had been improved and to determine the growth of indigenous aerobic denitrifiers, the Biolog method and plate count method were used for both the enhanced and control areas.

3.3. Carbon Metabolism Characteristics and Carbon Sources Utilization in the Enhanced and Control Areas

The AWCD value reflects the microbial microflora with different carbon sources and the change rate showed the metabolic activity of the microbial microflora [17]. A higher AWCD value and growth rate indicates a higher metabolic activity. During the operation of the WLAs, water samples (0.5, 45, and 90 m) were tested using the Biolog ECO plate on 15 October, 25 October, and 2 November, respectively. The results are shown in Figure 3A–C. The increase of AWCD mainly occurred in the first 120 h, and then stabilized. This means that the carbon source utilization mainly occurred in the first 120 h. According to the 120 h data point, for the bottom in the enhanced area, the AWCD value increased from 0.53 ± 0.08 (15 October) to 0.75 ± 0.11 (25 October), and then to 0.87 ± 0.12 (2 November). Compared with the enhanced area, these values in the control area were 0.60 ± 0.04 (15 October), 0.54 ± 0.05 (25 October), and 0.49 ± 0.02 (2 November), respectively, demonstrating a decreasing trend. As shown in Figure S8, the increase trend line of enhanced area of 0.5 m water was y = 0.0668x, y = 0.1029x, and y = 0.1318x, respectively. The increase rate was obviously improved. Compared with the control area of 0.5m water was y = 0.0852x, y = 0.0814x, and y = 0.0770x, respectively. It is clear that the in situ operation of WLAs resulted in improved metabolic activity of the indigenous microbes in the enhanced area, whereas that in the control area had reduced naturally.
Figure 3D–F showed the results of carbon source utilization. All the carbon sources were classified into six types: carbohydrates, amino acids, esters, alcohols, amines, and carboxylic acids. Absorbance was used to represent a carbon source utilization quantity. The carbon source utilization quantity was generally improved by the operation of WLAs. The total carbon sources utilization quantities in the enhanced area on 15 October were 4.49 ± 0.12 (0.5 m), 3.74 ± 0.15 (45 m), and 3.86 ± 0.12 (90 m). On 25 October, these were increased to 5.12 ± 0.21 (0.5 m), 4.41 ± 0.18 (45 m), and 4.87 ± 0.21 (90 m). On 2 November, the total carbon source utilization quantities in the enhanced area had greatly increased compared to those on 15 October, which were 6.60 ± 0.32 (0.5 m), 5.50 ± 0.25 (45 m), and 5.75 ± 0.12 (90 m). Compared with the enhanced area, these values in the control area showed little increase. The total carbon source utilization quantities showed the same trend as TOC variation. As seen in Figure S5, TOC concentrations in the enhanced area decreased from 3.69 ± 0.12 (15 October, 0.5 m), 3.42 ± 0.20 (15 October, 45 m), and 4.02 ± 0.18 (15 October, 90 m) to 2.72 ± 0.15 (2 November, 0.5 m), 2.65 ± 0.18 (2 November, 45 m), and 2.85 ± 0.20 (2 November, 90 m). The TOC removal rates of the enhanced area were 25.5 ± 0.8% (0.5 m), 23.8 ± 1.2% (45 m), and 28.2 ± 2.1% (90 m), and those of the control area were only 3.4 ± 0.4% (0.5 m), 6.5 ± 1.2% (45 m), and 5.6 ± 0.5% (90 m).
The results of the Biolog ECO Micro Plate and TOC (Figure S9) determination showed that the metabolic activity and the total carbon source utilization quantities of the indigenous aerobic denitrifiers were improved in the enhanced area compared with those in the control area. However, the amount of aerobic denitrifiers and the number of nitrogen removal functioning genes were unknown. Therefore, it was necessary to quantify the aerobic denitrifiers and nitrogen removal functioning genes in these areas.

3.4. Density Changes for Aerobic Denitrifiers and Denitrification Functional Genes nirK and nirS

As seen in Figure 4A, the densities of the indigenous aerobic denitrifiers were determined through the plate count method in both the enhanced and the control areas during the operation of WLAs. The number of aerobic denitrifiers in the enhanced area showed an increasing trend from 2.14 ± 0.08 × 104 cfu/mL (0.5 m, 15 October), 1.34 ± 0.05 × 104 cfu/mL (45 m, 15 October), and 0.78 ± 0.02 × 104 cfu/mL (90 m, 15 October) to 6.06 ± 0.18 × 104 cfu/mL (0.5 m, 2 November), 5.84 ± 0.09 × 104 cfu/mL (45 m, 2 November), and 4.67 ± 0.10 × 104 cfu/mL (90 m, 2 November). The number of aerobic denitrifiers was increased by 2–5 times after the operation of WLAs. No obvious increase in aerobic denitrifiers was observed in the control area and their number was maintained at 0.98 ± 0.06 × 104 cfu/mL to 1.85 ± 0.12 × 104 cfu/mL during the operation in all water layers.
Deviating from our previous study [23] in the Zhoucun Reservoir with the enclosure system, the number of aerobic denitrifiers in the enhanced area (0.5, 45, and 90 m) all showed an increasing trend throughout the monitoring period, and no decreasing trend was observed. The reason for this may be that in the field scale, there is no lack of TOC due to the severe exchange between water bodies (Figure S8). Moreover, the number of aerobic denitrifiers in the enhanced area was less than that in the enclosure system [23], and so the nitrogen removal rate in this study was also a slightly lower than that in the enclosure system.
Nitrite reductase occurs in two structurally different but functionally equivalent forms: nirK, containing copper, and nirS containing iron [23]. The nirK and nirS genes have been used to examine denitrifying bacteria in many environments [22,23,24,25,26,27]. Thus, quantification of denitrifying genes is of great importance for better understanding of denitrifying activity in the reservoir. As is known, few studies have focused on quantifying nirK and nirS. Similar to our previous study, we quantified nirK and nirS by real-time quantitative in both the enhanced and control areas.
As seen in Figure 4B,C, both nirK and nirS in the enhanced area were increased during the operation of WLAs. On 15 October, nirK in the enhanced area showed (4.82 ± 0.12) × 104 copies/μL (0.5 m), (3.85 ± 0.08) × 104 copies/μL (45 m), and (2.45 ± 0.13) × 104 copies/μL (0.5 m). After the operation, on 2 November, these values increased to (6.35 ± 0.24) × 104 copies/μL (0.5 m), (5.65 ± 0.15) × 104 copies/μL (45 m), and (5.52 ± 0.18) × 104 copies/μL (0.5 m). The nirK copy number in the control area showed little increase in all the water layers and was maintained at (2.41–4.76) × 104 copies/μL. The nirS copy number showed the same trend as that of nirK. The nirS copy number in the enhanced area increased by 0.24–0.63 times via the in situ operation. In the control area, the nirS genes showed a slight decrease.
From the above results, we can conclude that both aerobic denitrifier density and the denitrifying gene copy numbers were increased via the in situ operation of WLAs in the enhanced area compared with the control area. The Biolog ECO plate results showed that microbial metabolic activity had also improved greatly. Thus, the nitrogen was removed by the indigenous aerobic denitrifiers via the in situ operation of WLAs.

3.5. Spatial and Temporal Changes of Microbial Community

Miseq high-throughput sequencing was used to explore the microbial community changes in the enhanced and control areas. As shown in Table 2, of the 800,156 sequenced reads, 15,627 OTUs (operational taxonomic units) with 97% similarity were obtained. The number of OTUs obtained from the 12 water samples in the enhanced area wad 8264, which was higher than that obtained from samples of the control area (7363 OTUs). Thus, operation of the WLAs improved the microbial diversity in the enhanced area compared with that in the control area. The highest OTU number was obtained in the enhanced area at 45 m (15 October). The lowest number was obtained in the control area at 0.5 m (15 October).
The AEC and Chao diversity estimators, as shown in Table 2, varied mainly across depths. AEC diversity estimators ranged from 612 to 1205 and Chao diversity estimators varied from 499 to 1122 across depths. The average coverage of the enhanced and control areas was 0.9938 and 0.9943, respectively. This showed that the real structure of microbial communities was wonderfully reflected [28].
As shown in Figure 5, 15,627 OTUs with 97% similarity were affiliated with 11 main phyla including Proteobacteria, Bacteroidetes, Cyanobacteria, Planctomycetes, Nitrospirae, Actinobacteria, Firmicutes, Chloroflexi, Acidobacteria, Verrucomicrobia, and Deinococcus-thermus. Both in the enhanced and control areas, the phyla varied significantly across different periods. Similar to the results from the Zhoucun Reservoir study, Proteobacteria was the largest dominant population in the Jinpen Reservoir ranging from 29% to 50%. In the control area, on 15 October, Cyanobacteria were the dominant population, but it decreased during the monitoring period. In the enhanced area, Cyanobacteria decreased from 14.9% (15 October) to 1.7% (2 November). Bacteroidetes was also a dominant population in both the enhanced and control areas ranging from 6.1% to 26.4% varying from depth to depth. The Actinobacteria phylum varied significantly from 2.5% to 39.1%, and with the percentage in the enhanced area usually being higher than that in the control.
Microbial compositions varied greatly across different depths and dates. The abundance of the dominant aerobic denitrifying bacterial genera in the enhanced and the control areas during the operation of WLAs are shown in Figure 6. In the enhanced area, the hgcI_clade (6.0–22.5%), LD12_freshwater_group_norank (5.5–14.5%), Bacillus (1.0–11.2%), CL500-29_marine_group (4.5–13.8%), and Limnohabitans (1.9–13.5%) were the dominant genera. During the operation of WLAs, the hgcI_clade (mean value of three water layers) increased from 11.0% (15 October) to 17.9% (2 November), the LD12_freshwater_group_norank increased from 6.2% (15 October) to 11.3% (2 November), whereas the Bacillus decreased from 7.6% (15 October) to 3.1% (2 November).

3.6. Variation of Nitrogen Removal Bacteria

Considering the nitrogen removal performance of indigenous aerobic denitrifiers in the enhanced area via the in situ operation of WLAs, analyzing the variation in nitrogen functional bacteria is of great importance.
A previous study demonstrated that Proteobacteria (especially Acidovorax and Hydrogenophaga) were involved in nitrogen cycling [29]. In the enhanced area, the percentage of Acidovorax increased by 15 times from less than 0.01% (15 October) to 0.15% (2 November). Novosphingobium including the aerobic denitrifiers (ZHF2, GenBank no. KP717095) and ZHF8, GenBank no. KP717087) isolated from the Zhoucun Reservoir [8] increased from less than 0.01% (15 October) to 0.10% (2 November). The Hydrogenophaga in the enhanced area increased from 0% (15 October) to 0.05% (2 November). The aerobic denitrifiers Pseudomonas and Methylophaga increased from less than 0.01% and 0.02% to 0.06% and 0.05%, respectively. The aerobic denitrifiers Sphingomonadaceae and Bacteroidetes (ZMF6, GenBank no. KP717084) were also increased via the operation of WLAs.
Thus, it is clear that the activities and densities of the indigenous aerobic denitrifiers listed above were greatly improved via the WLAs operation. This is the main reason for nitrogen removal in the enhanced area. Future research should be focused on the variation in the microbial community of the surface sediment in both the enhanced and control areas. Moreover, the cultivation, enrichment, and utilization of indigenous aerobic denitrifiers in other water bodies should also be explored.

3.7. PCA of Microbial Communities

The principal component analysis methodology was used to explore the spatial and temporal variation in the bacterial community across water samples. As shown in Figure S9, the first two principle components (PC1 and PC2) closely approached 65.87% and the accumulated contributions of PC1 and PC2 were 42.90% and 22.97%, respectively, which was sufficient to explain the variation of the sampling sites. All the samples (enhanced and control areas) were displayed in the plot. The distribution of different samples is likely affected by a variety of factors, such as spatial and temporal differences and environment parameters. As shown in Figure 7, the samples from the same period showed tighter clustering, whereas the samples from different areas were distributed relatively more widely. Moreover, the samples were primarily separated by PC1. The samples of the control area were placed on the positive side of PC1, whereas those of the enhanced area were scattered on the negative side of PC1. Deviating from our previous study in the enclosure system in the Zhoucun Reservoir, microbial communities were not clearly separated between the enhanced and control areas in this field experiment, and the influence of WLAs in the control area (although 2 km far away) might be the possible reason for this effect.

3.8. Relationship of Microbial Community and Environmental Variables

Multiple statistical analyses (RDA) were used to determine the relationship between the water quality and the bacterial community. The different bacterial communities present in the enhanced and control water samples were well discriminated at the gene level (Figure 5). The results in Figure 8 showed than the first two RDA dimensions with 12 parameters (VIF < 30, except TN = 46) explained 35.23% of the microbial community variation (F = 1.0, p = 0.502) and the accumulated contribution ratios of RDA1 and RDA2 achieved 28.77% and 9.27%, respectively. The water samples were separated by RDA2 (samples of the enhanced area were located in the positive values of RDA2, whereas those of the enhanced area were scattered on the negative values of RDA2). Meanwhile, the parameters including the physical environment factors (CHl-a, pH, T, ORP, and DO), chemical parameters (TOC, TN, and TP), and the biological parameters (AD and denitrification functional genes) significantly influenced the microbial community composition. This was consistent with the operation of WLAs in the Jinpen reservoir ecosystem. It is clear that different sampling areas exhibited different bacterial communities. Thus, nitrogen removal by in situ enhancement of oxygen and indigenous aerobic denitrification bacteria using WLAs in the Jinpen reservoir was effective. However, the mechanism of nitrogen removal was not clear, and therefore, further exploration involving denitrifying enzymes, denitrification activities, and nitrogen metabolic genes (haoA, amoA, nosZ, nirK, and nirS) needs to be carried out in the future.

4. Conclusions

Field tests of nitrogen removal were performed through operation of WLAs. In the enhanced area, nitrogen removal rate reached 29.1 ± 0.8%. Microbial carbon metabolism activity and carbon source utilization were improved compared to those in the control area. Aerobic denitrifier densities and denitrifying gene (nirK and nirS) copies were also increased by the operation of WLAs. The number of aerobic denitrifiers was increased by 2–5 times and nirS copy number in the enhanced area increased by 0.24–0.63 times via the in situ operation. Miseq high-throughput sequencing indicated that denitrifying bacteria percentage was higher in the enhanced area, for example, in the enhanced area, the percentage of Acidovorax increased by 15 times from less than 0.01% (15 October) to 0.15% (2 November). In conclusion, WLAs can be used to improve water quality especially through nitrogen removal by the active indigenous aerobic denitrifiers.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/10/8/1051/s1.

Author Contributions

T.H. conceived the idea. Z.Z. designed the experiments and performed the experiments. Y.L., W.G., Y.L., Y.D., F.Z., S.Z., and W.M. analyzed the data and wrote the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 51478378 and 50830303.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (nos. 51478378 and 50830303). The authors would like to express their gratitude to the agencies involved and the study participants.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fenech, C.; Rock, L. The potential for a suite of isotope and chemical markers to differentiate sources of nitrate contamination: A review. Water Res. 2012, 46, 2023–2041. [Google Scholar] [CrossRef] [PubMed]
  2. Huang, T.; Lin, G. Nitrogen-removal efficiency of a novel aerobic denitrifying bacterium, Pseudomonas stutzeri strain ZF31, isolated from a drinking-water reservoir. Bioresour. Technol. 2015, 196, 209–216. [Google Scholar] [CrossRef] [PubMed]
  3. Song, Z.F.; Yang, X.L. Isolation and characterization of an aerobic denitrifying Bacillus sp. YX-6 from shrimp culture ponds. Aquaculture 2011, 319, 188–193. [Google Scholar] [CrossRef]
  4. Naeem, A. Enhanced biological nutrient removal by the alliance of a heterotrophic nitrifying strain with a nitrogen removing ecosystem. J. Environ. Sci. 2008, 20, 216–223. [Google Scholar]
  5. Patureau, D.; Helloin, E. Combined phosphate and nitrogen removal in a sequencing batch reactor using the aerobic denitrifier, Microvirgula aerodenitrificans. Water Res. 2001, 35, 189–197. [Google Scholar] [CrossRef]
  6. Qian, C.; Ni, J. Bioaugmentation treatment of municipal wastewater with heterotrophic-aerobic nitrogen removal bacteria in a pilot-scale SBR. Bioresour. Technol. 2015, 183, 25–32. [Google Scholar] [CrossRef]
  7. Robertson, L.A.; Kuenen, J.G. Thiosphaera pantotropha gen. nov. sp. nov., a Facultatively Anaerobic, Facultatively Autotrophic Sulphur Bacterium. Microbiology 1983, 129, 2847–2855. [Google Scholar] [CrossRef]
  8. Zhou, S.; Huang, T. Nitrogen removal characteristics of indigenous aerobic denitrifiers and changes in the microbial community of a reservoir enclosure system via in situ oxygen enhancement using water lifting and aeration technology. Bioresour. Technol. 2016, 214, 63–73. [Google Scholar] [CrossRef] [PubMed]
  9. Yao, S.; Ni, J. Enrichment and characterization of a bacteria consortium capable of heterotrophic nitrification and aerobic denitrification at low temperature. Bioresour. Technol. 2013, 127, 151–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Carter, J.P.; Hsaio, Y.H. Soil and sediment bacteria capable of aerobic nitrate respiration. Appl. Environ. Microbiol. 1995, 61, 2852–2858. [Google Scholar] [PubMed]
  11. Joo, H.S.; Hirai, M. Piggery wastewater treatment using Alcaligenes faecalis strain No. 4 with heterotrophic nitrification and aerobic denitrification. Water Res. 2006, 40, 3029–3036. [Google Scholar] [CrossRef] [PubMed]
  12. Guo, L.; Chen, Q. Application potential of a newly isolated indigenous aerobic denitrifier for nitrate and ammonium removal of eutrophic lake water. Bioresour. Technol. 2013, 142C, 45–51. [Google Scholar] [CrossRef] [PubMed]
  13. Patureau, D.; Godon, J.J. Microvirgula aerodenitrificans gen. nov., sp. nov., a new gram-negative bacterium exhibiting co-respiration of oxygen and nitrogen oxides up to oxygen-saturated conditions. Int. J. Syst. Bacteriol. 1998, 48 Pt 3, 775–782. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, J.; Wu, P. Heterotrophic nitrification and aerobic denitrification by the bacterium Pseudomonas stutzeri YZN-001. Bioresour. Technol. 2011, 102, 9866–9869. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, Q.L.; Liu, Y. The characteristics of a novel heterotrophic nitrification–aerobic denitrification bacterium, Bacillus methylotrophicus strain L7. Bioresour. Technol. 2012, 108, 35–44. [Google Scholar] [CrossRef] [PubMed]
  16. Gao, H.; Schreiber, F. Aerobic denitrification in permeable Wadden Sea sediments. ISME J. 2011, 5, 417–426. [Google Scholar] [CrossRef]
  17. Coban, O.; Kuschk, P. Nitrogen transforming community in a horizontal subsurface-flow constructed wetland. Water Res. 2015, 74, 203–212. [Google Scholar] [CrossRef] [PubMed]
  18. Chai, B.; Huang, T. A new method of inhibiting pollutant release from source water reservoir sediment by adding chemical stabilization agents combined with water-lifting aerator. J. Environ. Sci. 2011, 23, 1977–1982. [Google Scholar] [CrossRef]
  19. Gantzer, P.A.; Bryant, L.D. Controlling soluble iron and manganese in a water-supply reservoir using hypolimnetic oxygenation. Water Res. 2009, 43, 1285–1294. [Google Scholar] [CrossRef] [PubMed]
  20. Zhou, Z.Z.; Huang, T.L. Impacts of water quality variation and rainfall runoff on Jinpen Reservoir, in Northwest China. Water Sci. Eng. 2015, 8, 301–308. [Google Scholar] [CrossRef]
  21. Gao, H.; Schreiber, F.; Collins, G.; Jensen, M.M.; Lavik, G.; Kuypers, M.M.M. Aerobic Denitrification in Permeable Intertidal Sediments from the Wadden Sea; EGU General Assembly: Vienna, Austria, 2009. [Google Scholar]
  22. Huang, T.; Xuan, L. Effects of storm runoff on the thermal regime and water quality of a deep, stratified reservoir in a temperate monsoon zone, in Northwest China. Sci. Total Environ. 2014, 485, 820–827. [Google Scholar] [CrossRef] [PubMed]
  23. Zhou, S.; Huang, T. Nitrogen removal characteristics of enhanced in situ indigenous aerobic denitrification bacteria for micro-polluted reservoir source water. Bioresour. Technol. 2015, 201, 195–207. [Google Scholar] [CrossRef] [PubMed]
  24. Herbert, R.B.; Winbjörk, H. Nitrogen removal and spatial distribution ofdenitrifier and anammox communities in a bioreactor for mine drainage treatment. Water Res. 2014, 66, 350–360. [Google Scholar] [CrossRef] [PubMed]
  25. Ligi, T.; Truu, M. Effects of soil chemical characteristics and water regime on denitrification genes (nirS, nirK, and nosZ) abundances in a created riverine wetland complex. Ecol. Eng. 2013, 72, 47–55. [Google Scholar] [CrossRef]
  26. Wang, L.; Zheng, B. Diversity of bacterial community and detection of nirS- and nirK-encoding denitrifying bacteria in sandy intertidal sediments along Laizhou Bay of Bohai Sea, China. Mar. Pollut. Bull. 2014, 1, 215–223. [Google Scholar] [CrossRef] [PubMed]
  27. Zeng, W.; Zhang, J. Denitrifying phosphorus removal from municipal wastewater and dynamics of “Candidatus Accumulibacter” and denitrifying bacteria based on genes of ppk1, narG, nirS and nirK. Bioresour. Technol. 2016, 207, 322–331. [Google Scholar] [CrossRef] [PubMed]
  28. Luo, J.; Liang, H. Microbial community structures in a closed raw water distribution system biofilm as revealed by 454-pyrosequencing analysis and the effect of microbial biofilm communities on raw water quality. Bioresour. Technol. 2013, 148, 189–195. [Google Scholar] [CrossRef] [PubMed]
  29. Yan, Q.; Bi, Y. Impacts of the Three Gorges Dam on microbial structure and potential function. Sci. Rep. 2015, 5, 8605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Variation of T, DO, pH, turbidity, ORP, and CHl-a in the enhanced area during the operation of WLAs (E—indicates “enhanced area”).
Figure 1. Variation of T, DO, pH, turbidity, ORP, and CHl-a in the enhanced area during the operation of WLAs (E—indicates “enhanced area”).
Water 10 01051 g001
Figure 2. Variations in the concentration of total nitrogen, nitrate, ammonia, and nitrite during the operation of WLAs in the enhanced area.
Figure 2. Variations in the concentration of total nitrogen, nitrate, ammonia, and nitrite during the operation of WLAs in the enhanced area.
Water 10 01051 g002
Figure 3. AWCD and carbon source utilization variation of microbes in the enhanced and control areas during the operation of WLAs ((A,D): 15 October, (B,E): 25 October, (C,F): 2 November; E—indicates “enhanced area”; C—indicates “control area”).
Figure 3. AWCD and carbon source utilization variation of microbes in the enhanced and control areas during the operation of WLAs ((A,D): 15 October, (B,E): 25 October, (C,F): 2 November; E—indicates “enhanced area”; C—indicates “control area”).
Water 10 01051 g003
Figure 4. Variation of aerobic denitrifier density and the genes nirS and nirK in the enhanced and control areas ((A): 15 October, (B): 25 October, (C): 2 November; E—indicates “enhanced area”; C—indicates “control area”).
Figure 4. Variation of aerobic denitrifier density and the genes nirS and nirK in the enhanced and control areas ((A): 15 October, (B): 25 October, (C): 2 November; E—indicates “enhanced area”; C—indicates “control area”).
Water 10 01051 g004
Figure 5. Variation of microbial phyla in the enhanced and control areas (E—indicates “enhanced area”; C—indicates “control area”).
Figure 5. Variation of microbial phyla in the enhanced and control areas (E—indicates “enhanced area”; C—indicates “control area”).
Water 10 01051 g005
Figure 6. Abundance of the dominating aerobic denitrifying genus in the enhanced and control water layers during the operation of WLAs (“E-01-0m” indicates “enhanced area, the first day of WLAs operation, surface water 0.5 m”; “C-09-90” indicates “control area, the ninth day of WLAs operation, bottom water 90 m”).
Figure 6. Abundance of the dominating aerobic denitrifying genus in the enhanced and control water layers during the operation of WLAs (“E-01-0m” indicates “enhanced area, the first day of WLAs operation, surface water 0.5 m”; “C-09-90” indicates “control area, the ninth day of WLAs operation, bottom water 90 m”).
Water 10 01051 g006aWater 10 01051 g006b
Figure 7. Principal component analysis (PCA) of water samples in the enhanced and control areas.
Figure 7. Principal component analysis (PCA) of water samples in the enhanced and control areas.
Water 10 01051 g007
Figure 8. Redundancy analysis of water bacterial communities in the enhanced and control areas.
Figure 8. Redundancy analysis of water bacterial communities in the enhanced and control areas.
Water 10 01051 g008
Table 1. Changes in OC and TN concentrations in the surface sediment (0–1 cm) at the enhanced and control areas during the operation of WLAs, in triplicate
Table 1. Changes in OC and TN concentrations in the surface sediment (0–1 cm) at the enhanced and control areas during the operation of WLAs, in triplicate
Sampling SiteDateOC (%)TN (mg/kg)
Enhanced Area15 October3.42 ± 0.201522 ± 66
21 October2.92 ± 0.111471 ± 43
28 October2.78 ± 0.111313 ± 56
2 November2.66 ± 0.121242 ± 61
Control Area15 October3.33 ± 0.121644 ± 72
21 October3.48 ± 0.091687 ± 54
28 October3.43 ± 0.081744 ± 36
2 November3.64 ± 0.121751 ± 44
Note: “enhanced area” was nearly 50 m away from the WLAs, “control area” was 2 km away from the WLAs. “OC” means organic carbon in sediment; “TN” means total nitrogen in sediment.
Table 2. Spatial and temporal distribution of the microbial community diversity and richness estimators in the enhanced and control areas
Table 2. Spatial and temporal distribution of the microbial community diversity and richness estimators in the enhanced and control areas
PeriodWater DepthEnhanced AreaControl Area
Reads Number0.97 LevelReads Number0.97 Level
OTUsDiversityCoverageRichnessOTUsDiversityCoverageRichness
ACEChao1ShannonSimpsonACEChao1ShannonSimpson
15 October0.5 m24,6397129319230.99164.550.027538,0613636124990.99703.260.1035
45 m39,804852103410210.99514.850.019830,5525447087380.99493.730.0921
90 m40,082765101910290.99434.380.031329,84162910269050.99324.280.0346
21 October0.5 m27,963650110510040.99204.320.02926,9866859189080.99244.410.0327
45 m29,8175797407800.99474.220.033834,2687029149470.99414.370.034
90 m30,4186548588770.99384.450.025930,516846111111220.99184.910.0171
28 October0.5 m33,4316499538430.99464.470.025627,0364656006410.99524.040.0441
45 m31,3027159609810.99304.470.026133,1725957587730.99534.230.0379
90 m38,056768100710180.99424.330.039739,4546548208390.99584.380.0333
2 November0.5 m19,6965669058100.99054.150.040242,4394616176160.99683.590.0566
45 m46,08269910459520.99574.340.030531,7837029489340.99324.090.0473
90 m45,37965510418820.99553.950.050229,379717120510650.99144.150.0405

Share and Cite

MDPI and ACS Style

Zhou, Z.; Huang, T.; Gong, W.; Li, Y.; Liu, Y.; Zhao, F.; Dou, Y.; Zhou, S.; Ma, W. In Situ Water Quality Improvement Mechanism (Nitrogen Removal) by Water-Lifting Aerators in a Drinking Water Reservoir. Water 2018, 10, 1051. https://doi.org/10.3390/w10081051

AMA Style

Zhou Z, Huang T, Gong W, Li Y, Liu Y, Zhao F, Dou Y, Zhou S, Ma W. In Situ Water Quality Improvement Mechanism (Nitrogen Removal) by Water-Lifting Aerators in a Drinking Water Reservoir. Water. 2018; 10(8):1051. https://doi.org/10.3390/w10081051

Chicago/Turabian Style

Zhou, Zizhen, Tinglin Huang, Weijin Gong, Yang Li, Yue Liu, Fuwang Zhao, Yanyan Dou, Shilei Zhou, and Weixing Ma. 2018. "In Situ Water Quality Improvement Mechanism (Nitrogen Removal) by Water-Lifting Aerators in a Drinking Water Reservoir" Water 10, no. 8: 1051. https://doi.org/10.3390/w10081051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop