Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)
Article
Edaphic conditions regulate denitrification directly and indirectly by altering denitrifier abundance in wetlands along the Han River, China Ziqian Xiong, Laodong Guo, Quanfa Zhang, Guihua Liu, and Wenzhi Liu Environ. Sci. Technol., Just Accepted Manuscript • Publication Date (Web): 25 Apr 2017 Downloaded from http://pubs.acs.org on April 25, 2017
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 36
Environmental Science & Technology
1
Edaphic conditions regulate denitrification directly and indirectly by
2
altering denitrifier abundance in wetlands along the Han River, China
3 4
Ziqian Xiong a, b, Laodong Guo c, Quanfa Zhang a, Guihua Liu a*, Wenzhi Liu a, c*
5 6
a
7
Chinese Academy of Sciences, Wuhan 430074, P.R. China
8
b
9
P.R. China
10 11
c
Key Laboratory of Aquatic Botany and Watershed Ecology, Wuhan Botanical Garden,
College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049,
School of Freshwater Sciences, University of Wisconsin-Milwaukee, Milwaukee,
53204, USA
12 13
*Corresponding author
14
Phone: +86 27 87510987
15
Fax: +86 27 87510251
16
Email:
[email protected] (Guihua Liu);
17
[email protected] (Wenzhi Liu)
18 19
Type: Research Articles.
20
Author contributions: QZ, GL, and WL designed the study; ZX, GL, and WL performed
21
the research; ZX, LG, and WL analysed the data; ZX, LG, QZ, GL, and WL wrote the
22
manuscript.
23 24 1
ACS Paragon Plus Environment
Environmental Science & Technology
25
Abstract
26
Riparian wetlands play a critical role in retaining nitrogen (N) from upland runoff
27
and improving river water quality, mainly through biological processes such as soil
28
denitrification. However, the relative contribution of abiotic and biotic factors to
29
riparian denitrification capacity remains elusive. Here we report the spatio-temporal
30
dynamics of potential and unamended soil denitrification rates in 20 wetlands along
31
the Han River, an important water source in central China. We also quantified the
32
abundance of soil denitrifying microorganisms using nirK and nirS genes. Results
33
showed that soil denitrification rates were significantly different between riparian
34
and reservoir shoreline wetlands, but not between mountain and lowland wetlands.
35
In addition, soil denitrification rates showed strong seasonality, with higher values in
36
August (summer) and April (spring) but lower values in January (winter). The
37
potential and unamended denitrification rates were positively correlated with
38
edaphic conditions (moisture and carbon concentration), denitrifier abundance, and
39
plant species richness. Path analysis further revealed that edaphic conditions could
40
regulate denitrification rates both directly and indirectly through their effects on
41
denitrifier abundance. Our findings highlight that not only environmental factors, but
42
also biotic factors including denitrifying microorganisms and standing vegetation,
43
play an important role in regulating denitrification rate and N removal capacity in
44
riparian wetlands.
45 46
Keywords: Denitrifying communities, Microbial abundance, Potential denitrification,
47
Riparian zone, Wetland vegetation
48 2
ACS Paragon Plus Environment
Page 2 of 36
Page 3 of 36
Environmental Science & Technology
TOC/Abstract Art
3
ACS Paragon Plus Environment
Environmental Science & Technology
49
Introduction
50
It is estimated that anthropogenic activities in the world contribute
51
approximately 2.1×1011 kg of reactive nitrogen (N) to terrestrial ecosystems every
52
year.1 Nitrogen mainly comes from industrial and agricultural sources such as
53
wastewater discharge, chemical fertilizers, animal wastes, atmospheric deposition
54
and crop fixation.2‒4 Nitrogen inputs to the terrestrial environments now greatly
55
exceed biological demand, and this excess N can easily leach into rivers, streams, and
56
other aquatic ecosystems.5 Elevated N concentrations in river and stream water
57
commonly result in a number of environmental issues such as deterioration of water
58
quality, eutrophication, toxic algal blooms, and loss of biodiversity.6
59
Riparian wetlands, the transition zone between terrestrial and aquatic systems,
60
play an important role in removing N from runoff and improving river water quality.7,8
61
It has been reported that average N removal efficiency for global riparian zone is
62
approximately 67.5%.9 Nitrogen can be removed or retained by riparian wetlands
63
mainly through soil denitrification, anammox, plant uptake, soil storage, and
64
microbial immobilization.10 Soil denitrification is thought to be the most critical
65
process for N removal, because it can permanently remove the N through reducing
66
nitrate (NO3–) to N gaseous products including nitrous oxide (N2O) and dinitrogen
67
(N2).11 For instance, Kreiling et al. (2011) reported that the denitrification process
68
could account for nearly 82% of the NO3– removal in the upper Mississippi River.12
69
Edaphic characteristics, such as soil moisture, N and carbon (C) availability, can
70
directly and indirectly affect the denitrification rate in riparian wetlands.13,14 For
71
example, Hefting et al. (2006) showed that soil pH, moisture and organic matter were
72
the main determinants of denitrification activity in a forested riparian buffer.15 Hunt 4
ACS Paragon Plus Environment
Page 4 of 36
Page 5 of 36
Environmental Science & Technology
73
et al. (2004) reported that soil NO3– and C levels could account for the majority of
74
variance in soil denitrification rates in a forested riparian wetland.13 Several studies
75
have demonstrated that large-scale riparian features such as geomorphological
76
factors (e.g., slope and altitude) and land use are also significantly correlated with N
77
removal capacity.16,17 Geomorphic character can regulate soil conditions and
78
vegetation characteristics in riparian zones, which may in turn impact the soil
79
denitrification rates.17
80
With the development of molecular technologies, a number of studies have
81
focused on the relationships between soil denitrification rates and denitrifier
82
abundance, but their results are somewhat inconsistent.18‒20 Deslippe et al. (2014)
83
found that the abundance of nirS–type denitrifier was positively correlated to
84
denitrification potential in riparian soils.21 Graham et al. (2014), on the other hand,
85
did not observe a significant increase in explanatory power when functional gene
86
abundances were added to statistical models to predict soil denitrification rates.22
87
Moreover, vegetation has long been recognized as an important biotic factor in
88
determining wetland denitrification process.23 A meta–analysis revealed that wetland
89
macrophytes increased soil denitrification rates by 55% on average, mainly by
90
providing organic matter and introducing oxygen that enhanced generation of NO3–
91
from nitrification.24
92
In this study, we examined seasonal variation in potential and unamended soil
93
denitrification rates in 20 wetlands along the Han River, a water source area of the
94
South-to-North Water Transfer Project in China. We used nirK and nirS genes as
95
markers to quantify the denitrifier abundance in wetland soils. We hypothesized that
96
edaphic conditions such as available C and N could affect denitrification rates both 5
ACS Paragon Plus Environment
Environmental Science & Technology
97
directly and indirectly through their effects on denitrifier abundance. Therefore, our
98
study had several purposes, including (1) investigating the spatio-temporal patterns
99
of soil potential and unamended denitrification rates in wetlands along the Han River,
100
China; (2) quantifying the relative contributions of abiotic and biotic factors to
101
wetland denitrification rates; and (3) elucidating the direct and indirect pathways of
102
the effect of edaphic conditions on wetland denitrification rates.
103 104
Materials and methods
105
Study areas
106
The Han River is the largest tributary of the Yangtze River and drains an area of
107
approximately 159,000 km2 (Fig. 1). It originates from Ningqiang County of Shaanxi
108
Province and flows southeast, emptying into the Yangtze River in Wuhan City. The
109
Han River basin has a north subtropical monsoon climate, where the annual average
110
temperature is 12–16 °C and annual mean precipitation is approximately 804 mm.17
111
Since 1968, eight large reservoirs have been built on the Han River to provide
112
benefits including flood control, freshwater supply, and electric power generation.
113
Among the eight reservoirs, the largest one is the Danjiangkou Reservoir, which has a
114
total area of 745 km2 and a mean water depth of approximately 23 m.14 The majority
115
of these reservoirs have large areas of shoreline wetlands in the transition zone
116
between low and high water levels. In contrast to riparian wetlands, reservoir
117
shoreline wetlands experience large water-level fluctuations (3–15 m annually), and
118
are mainly formed at the former farmlands and forests.14
119
The Han River can be divided into two geomorphically distinct reaches by the
120
Danjiangkou Reservoir.25 The upper reach, from the source to the Danjiangkou 6
ACS Paragon Plus Environment
Page 6 of 36
Page 7 of 36
Environmental Science & Technology
121
Reservoir with a length of 925 km, has typical mountain landscapes (Fig. 1).
122
Therefore, riparian wetlands in the upper Han River basin are commonly restricted to
123
a narrow band (2‒10 m in width) adjacent to the river channel.25 Land cover in this
124
mountain sub-basin is dominated by natural vegetation, especially deciduous and
125
coniferous forests (Table S1). In the lower Han River basin, the mean elevation is
126
about 283 m with a relatively flat topography (Table S1). Agriculture is the dominant
127
landscape in this lowland sub-basin.17
128
In 2002, the Chinese government launched the South-to-North Water Transfer
129
Project to divert freshwater to the northern cities (e.g., Beijing and Tianjin) from the
130
Yangtze River basin.26 The upper Han River is the water source area for the middle
131
rout of this inter-basin water transfer project. Therefore, the water quality and
132
wetland functions of the Han River have been of great concern during recent
133
years.17,27,28 Due to increased nutrients input from industrial wastewater, agricultural
134
fertilizers, and domestic sewage, the Han River has suffered from serious N
135
pollution.27 The concentration of total N (mainly NO3–) in river waters varied between
136
0.07 and 16.73 mg/L, with a mean value of 2.34 mg/L.29 Excess N has resulted in
137
water quality degradation and other attendant ecological issues, including
138
eutrophication and harmful algal blooms.30
139 140
Field sampling
141
In January (representing winter), April (spring), August (summer), and
142
November (autumn) of 2015, field sampling was conducted at 15 riparian wetlands
143
and 5 reservoir shoreline wetlands along the Han River, China (Fig. 1). These
144
wetlands were selected randomly and their coordinates were recorded using a 7
ACS Paragon Plus Environment
Environmental Science & Technology
145
hand-held global positioning system (Unistrong Co., Ltd, Beijing, China). At each
146
sampling wetland, one sampling transect was randomly established perpendicular to
147
river flow and extended from the water's edge to the upland edge of the wetlands
148
(i.e., the highest water-level areas).
149
The transect length varied from less than 5 m to approximately 35 m. One to
150
three plots (1× 1 m) were set up along each transect in 10-m intervals depending on
151
the transect length. Because of the fluctuating water levels of the Han River, a total of
152
35, 39, 46 and 37 plots were established along the 20 transects in January, April,
153
August and November, respectively.
154
At each plot, five soil cores (10 cm deep) were collected randomly using a hand
155
corer and then well mixed immediately to form a composite sample. Approximately
156
300 g of soil from each plot was put into a plastic bag and stored at about 5 °C in a
157
refrigerator (FYL-YS-30L, Fuyilian Co., Beijing, China). In addition, about 10 g of soil
158
was collected at each plot in a centrifuge tube and immediately frozen in liquid N2. At
159
each sampling transect, a 500–ml surface water sample (approximately 0.2 m below
160
surface) was collected from a river location adjacent to this transect.
161 162
Measurements of soil denitrification rates
163
We used the acetylene (C2H2) inhibition method to determine the potential and
164
unamended denitrification rates of wetland soils.31 The C2H2 inhibition technique has
165
a long history of successful use,7,16,23 although it has several drawbacks such as
166
inhibition of nitrification process.32 Potential denitrification rate was an upper-bound
167
estimate of in situ denitrification with C and NO3– amendments, while unamended
168
denitrification (background or basal denitrification) rate provided a conservative 8
ACS Paragon Plus Environment
Page 8 of 36
Page 9 of 36
Environmental Science & Technology
169
estimate and measured without nutrient addition.33
170
Potential denitrification rate was measured according to Xiong et al. (2015).17
171
Specifically, approximately 25 g of fresh soils from each plot were weighed into a
172
250-mL brown glass bottle with 20 mL of incubation solution (final concentrations:
173
0.18 g/L glucose, 0.1 g/L KNO3, and 1 g/L chloramphenicol). Unamended
174
denitrification rate was determined using a similar procedure, but with the addition
175
of 20 mL of unfiltered river water instead of incubation solution. All bottles were
176
sealed with a rubber stopper and purged with 99.999% N2 for 2 min to create anoxic
177
environments. To inhibit the reduction of N2O to N2 during denitrification process, 25
178
mL of C2H2 was then added to each bottle using a syringe. These bottles were
179
incubated in the dark at room temperatures (4.5, 15.1, 25.8 and 11.4 ℃ in January,
180
April, August and November, respectively) for 2 h, as recommended by Opdyke et al.
181
(2006) and Bruesewitz et al. (2012).33,34 At the start and end of incubation, 5 mL of
182
headspace gas samples were collected from each bottle using a syringe after shaking
183
vigorously.
184
The N2O concentrations were determined using a gas chromatograph fitted with
185
an electron capture detector (Agilent 7890, Santa Clara, USA). Potential and
186
unamended soil denitrification rates were calculated as the difference between the
187
start and end N2O concentrations (corrected for N2O dissolved in water using Bunsen
188
coefficient) divided by incubation time (2 h) and soil dry weight, and data were
189
reported in ng N g‒1 h‒1.
190 191 192
Measurements of vegetation, edaphic conditions and water quality In the field, vegetation cover at each plot was visually estimated using a 1×1 m 9
ACS Paragon Plus Environment
Environmental Science & Technology
Page 10 of 36
193
grid frame that was subdivided into 100 cells (10×10 cm).28 Species richness was the
194
number of vascular plant species present in a plot. In the laboratory, soil pH was
195
determined in a soil to water ratio of 1: 5 (v/v) by a pH meter, while soil moisture was
196
measured gravimetrically after drying 30 g of fresh soil at 105°C for 24 h. In order to
197
determine the percentage of fine substrates (i.e., sand, silt and clay), approximately
198
50 g of air-dried soil were sieved through a 10–mesh sieve.17 The soil density was
199
analyzed by weighing 50 cm3 of homogenized fresh soils after drying overnight at
200
105°C. The concentration of soil NO3–‒N was measured by extracting 10 g of fresh
201
soils with 100 mL of KCl for 1 h and using an automatic analyzer (EasyChem plus,
202
Systea, Italy). The soil total carbon (STC) was determined by an elemental analyzer
203
(Vario TOC cube, Elementar, Germany) using air-dried and sieved (100-mesh) samples.
204
Water NH4+‒N and NO3–‒N concentrations were determined by an automatic
205
analyzer (EasyChem plus, Systea, Italy).
206 207
Determination of denitrifier abundances
208
Genomic DNA was extracted from replicate sediment subsamples using the
209
PowerSoil DNA Isolation Kit (MoBio Laboratories, Inc., Carlsbad, USA) following the
210
manufacturer’s instructions. DNA quality was detected in 1% TAE‒agarose gel stained
211
with Goldview and imaged under UV light. The abundance (i.e., copy number) of nirS
212
and nirK genes was determined in triplicate using a Roche LightCycler480 real-time
213
PCR System (Roche Diagnostics, Mannheim, Germany) with the fluorescent dye SYBR
214
green
215
nirSCd3aF/nirSR3cd were applied for nirK and nirS genes, respectively. The 20 μL
216
quantitative PCR mixture contained 10 μL of SybrGreen qPCR Master Mix (2×), 1 μL of
quantitative
PCR
method.
Primer
sets
of
10
ACS Paragon Plus Environment
nirK876/nirK1040
and
Page 11 of 36
Environmental Science & Technology
217
primers (10 μM) and 2 μL of DNA template. Standard curves were constructed with
218
10-fold serial dilutions of a known amount of plasmid DNA containing fragments of
219
nirK and nirS genes. The primer sequence and thermal cycling procedures were listed
220
in Table S2.
221 222
Statistical analyses
223
All data were assessed for normality by using the Shapiro–Wilk test before
224
further statistical analyses. Non-normally distributed data were natural log (ln) or
225
square root (sqrt) transformed to improve their normality. T-test and One-way
226
ANOVA with Tukey post-hoc tests were conducted to examine differences in
227
denitrification rates, denitrifier abundances, vegetation characteristics and edaphic
228
conditions among seasons and between wetland types. The relationships between
229
denitrification rates and biotic and abiotic factors were evaluated using the Pearson
230
correlation coefficient. Stepwise multiple regression analyses were used to identify
231
which factors best explained the variations in denitrification rates. The above
232
statistical analyses were performed using PASW 19.0 statistics software (IBM SPSS
233
Inc., Chicago, USA).
234
To further explore the direct and indirect effects of edaphic conditions on
235
denitrification rates, a path analysis (structural equation modeling without latent
236
variables) consisted of three steps were conducted. Firstly, based on previous
237
research,18 we developed a conceptual model linking abiotic and biotic factors to soil
238
denitrification rates (Figure S1). Secondly, promising explanatory factors were
239
selected to include in path analysis according to results of the above Pearson
240
correlation analysis. To simplify our path model, we used a compound variable 11
ACS Paragon Plus Environment
Environmental Science & Technology
241
(nirK+nirS abundance) instead of nirK abundance and nirS abundance. Finally, path
242
coefficients, R2, and model fit parameters were calculated using AMOS 20.0 software
243
(Amos Development Corporation, Chicago, USA). Maximum likelihood estimation
244
method was used and missing data were handled using the option of “Estimate
245
means and intercepts”. The indirect effects (effects mediated by other factors) refer
246
to the mathematical product of all possible paths.28 A comparative fit index (CFI)
247
value > 0.90 indicated that the final path models provide an acceptable fit.
248 249
Results
250
Wetland denitrification rates
251
Potential denitrification rates showed great variation, ranging from 0.06 to 511
252
ng N g‒1 h‒1 and averaging 51.79 ng N g‒1 h‒1. Unamended denitrification rates varied
253
between 0.03 and 105 ng N g‒1 h‒1, with a mean value of 10.56 ng N g‒1 h‒1.
254
Denitrification rates in lowland wetlands were slightly but not significantly higher
255
than those in mountain wetlands (Fig. 2). However, reservoir shoreline wetlands
256
supported greater denitrification rates than riparian wetlands (Fig. 2). Sampling
257
seasons had a strong influence on both potential and unamended denitrification
258
rates (Fig. 3). The highest denitrification rates occurred in August and April, while the
259
lowest denitrification rates were observed in January.
260 261
Wetland biotic and abiotic factors
262
The abundance of nirK gene in wetland soils ranged from 0.31×105 to 84.50×105
263
gene copies g−1 soil, while the abundance of the nirS gene varied between 0.11×105
264
and 69.20×105 gene copies g−1 soil. The mean abundance of nirK gene (5.22×105 gene 12
ACS Paragon Plus Environment
Page 12 of 36
Page 13 of 36
Environmental Science & Technology
265
copies g−1 soil) was approximately twice that of nirS gene (2.69×105 gene copies g−1
266
soil). There was no significant difference in nirS and nirS gene abundances among
267
seasons (p>0.05; data not shown) and between wetland types (Table S3).
268
Soil moisture was generally higher in lowland wetlands than in mountain
269
wetlands (Table S4). In addition, relative to the riparian wetlands, reservoir shoreline
270
wetlands had significantly higher soil pH and STC (Table S3). There existed a
271
significant difference in all edaphic conditions except fine substrate between
272
sampling seasons (Table S4). The Soil NO3–‒N and STC concentrations in August were
273
significantly higher than those in January and April. Soil moisture was positively
274
related to nirK gene abundance, plant species richness and plant cover, while soil
275
NO3–‒N concentration was negatively correlated with nirS gene abundance and plant
276
species richness (Table 1).
277 278
Determinants of wetland denitrification rates
279
Soil denitrification rates were strongly associated with several abiotic and biotic
280
factors, when all type wetlands were considered (Table 2). Both potential and
281
unamended denitrification rates showed significant positive relationships with
282
denitrifier abundance, plant richness, soil moisture and STC. Stepwise regression
283
analyses demonstrated that nirK abundance was the most important predictor of
284
potential and unamended denitrification rates, explaining 32.5% and 18.6% of their
285
variations, respectively.
286
The CFI values for potential and unamended denitrification models were 0.949
287
and 0.933, respectively, indicating that the path models were acceptable (Fig. 4). For
288
the model of potential denitrification rate (Fig. 4a), denitrifier abundance was the 13
ACS Paragon Plus Environment
Environmental Science & Technology
289
most important factor in explaining the denitrification rate (Table 3). Edaphic
290
conditions, especially moisture, could affect potential denitrification rates both
291
directly and indirectly by altering denitrifier abundance (Table 3). A similar result was
292
obtained from the unamended denitrification model (Fig. 4b). The final models,
293
which included a subset of direct and indirect effects of abiotic and biotic factors,
294
accounted for 61.5% of the variation in potential denitrification rates and 55.4% of
295
the variation in unamended denitrification rates (Fig. 4).
296 297
Discussion
298
Spatial and temporal dynamics of wetland denitrification rates
299
We compared soil denitrification rates in wetlands along the Han River with
300
those in other riparian zones measured by the same C2H2 inhibition method (Table
301
S5). The mean potential denitrification rate in wetlands along the Han Rive was 51.8
302
ng N g‒1 h‒1. This value was considerably higher than values observed in some
303
previous studies, including 25.3 ng N g‒1 h‒1 in a riparian wetland in Maryland, USA
304
and 1.52 ng N g‒1 h‒1 in riparian zones adjacent to the Danjiangkou Reservoir,
305
China.25,35 However, the majority of previous studies reported a potential
306
denitrification rate between 100 and 3000 ng N g‒1 h‒1 (Table S5). This difference
307
could be mainly due to three reasons. First, there was a great difference in
308
concentrations of added N and C in determining potential denitrification rates among
309
studies. For instance, the added N content in our study was approximately 11.09 μg
310
NO3–‒N g-1 soil, which was considerably lower than that in Dhondt et al. (2004) (50
311
μg NO3–‒N g-1 soil) and in McCarty et al. (2007) (600 μg NO3–‒N g-1 soil).35,36 Second,
312
the sampling seasons and incubation temperatures had a significant effect on soil 14
ACS Paragon Plus Environment
Page 14 of 36
Page 15 of 36
Environmental Science & Technology
313
potential denitrification rates.37 Dhondt et al. (2004) found that potential
314
denitrification rates in soils increased considerably with incubation temperatures.36
315
Third, many previous studies used the chloramphenicol to inhibit the microbial
316
growth and block the de novo synthesis of denitrification enzymes.8,13,31,35 In our
317
incubation experiments, the addition of chloramphenicol might significantly reduce
318
the potential and unamended soil denitrification rates.
319
Scaled to an areal basis using the density of surface soils (top 10 cm), the mean
320
unamended denitrification rate in wetlands along the Han Rive was 0.79 mg N m‒2
321
h‒1 or 6920 kg N km‒2 y‒1. Piña-Ochoa and Álvarez-Cobelas (2006) reported that the
322
mean denitrification rate of global rivers was approximately 29540 kg N km‒2 y‒1.38 In
323
addition, Seitzinger et al. (2006) indicated that there were considerable spatial
324
variations in global river denitrification, with a predicted average denitrification rate
325
of 12963 kg N km‒2 y‒1.39 This may suggest that the wetlands along the Han River
326
have relatively low denitrification and N removal rates compared with global river
327
wetlands.
328
Furthermore, reservoir shoreline wetlands along the Han River supported
329
greater denitrification rates than riparian wetlands (Fig. 2). This result can be mainly
330
ascribed to the fact that, in comparison with riparian wetlands, reservoir shoreline
331
wetlands have higher soil C concentrations which are positively correlated to the
332
potential and unamended denitrification rates (Table S3 and Table 2). We also found
333
that sampling season had a strong influence on both potential and unamended
334
denitrification rates in wetlands along the Han River (Fig. 3). The highest
335
denitrification rates were observed in August and April, while the lowest
336
denitrification rates occurred in January. It has been reported that denitrification 15
ACS Paragon Plus Environment
Environmental Science & Technology
337
rates are positively correlated with temperature, with higher values in summer than
338
in winter.40 In the Han River basin, the mean air temperature in August (summer) and
339
April (spring) is 24.8 and 15.2 °C, respectively, while mean temperature in January
340
(winter) is only 2.5 °C.41 Thus, temperature is one of the important determinants of
341
soil denitrification rates in wetlands along the Han River. It should be noted that C
342
and N can also be temporarily adsorbed on soil particles or retained in pore waters.17
343
These effects may somewhat moderate the difference in soil denitrification rates
344
among seasons.
345 346
Direct and indirect effects of edaphic conditions on wetland denitrification rates
347
Edaphic conditions, including moisture, NO3– concentration and C availability,
348
can directly affect the instantaneous rate of soil denitrification.17,18 Consistent with
349
earlier studies,42,43 there existed positive correlations between soil moisture and
350
denitrification rates in wetlands along the Han River (Table 2). The oxygen availability
351
changed rapidly depending on soil moisture and the consequent rate of oxygen
352
diffusion through soils.43 High soil moisture can inhibit oxygen diffusion and thereby
353
create anaerobic or anoxic environments favorable for denitrification. Our study also
354
revealed that soil NO3– content had a positive relationship with the wetland
355
denitrification rate (Table 2). Such a relationship was frequently observed in previous
356
studies,17,44 because NO3– was the terminal electron receptor of the denitrification
357
process. Many studies have examined the influences of C or organic matter contents
358
on soil denitrification rates.7,44 In this study, we observed positive correlations
359
between STC and potential and unamended denitrification rates (Table 2), consistent
360
with many previous studies.17,44 The positive relationships between soil NO3– or C 16
ACS Paragon Plus Environment
Page 16 of 36
Page 17 of 36
Environmental Science & Technology
361
concentrations and denitrification rates suggest that the denitrification process in
362
wetlands along the Han River is probably both C‒ and N‒ limited.
363
Knowledge of how edaphic conditions interactively control denitrifying microbial
364
communities and their ecological functions is relevant to better understanding N
365
removal capacity of wetlands.18 Several studies have assumed that edaphic
366
conditions can affect soil denitrification rates indirectly by altering microbial
367
communities.16,45 However, evidence for this assumption remains limited.18,46 Here, it
368
is evident that the indirect effect of soil moisture, NO3– and STC concentrations on
369
the wetland denitrification rates was mediated through denitrifier abundance (Fig. 4).
370
Soil moisture was positively correlated to nirK gene abundances, while soil NO3– and
371
STC concentrations were negatively correlated with nirS gene abundances (Table 1),
372
implying habitat selection on nirK and nirS denitrifiers and supporting the
373
speculations on niche differentiation.47 In riparian wetlands, soil moisture was
374
negatively correlated with oxygen availability, which could be a critical factor in
375
determining the structure of denitrifying communities.17 Previous studies have
376
reported that nirS–type denitrifiers are found in constantly anoxic conditions and
377
nirK–type denitrifiers are found in environments with fluctuating oxygen levels.16
378
Organic carbon, the primary electron donor for respiratory denitrifying bacteria, has
379
been shown to play an important role in regulating the abundance of denitrifying
380
communities.47 In general, the abundance of nirS gene increased with increasing
381
NO3– concentration in soils, because NO3– served as a terminal electron acceptor.18
382
These inconsistent results demonstrated that the abundance of nirS gene may
383
respond differently to the NO3– loads in different environments. Future work is
384
needed to examine the mechanisms by which soil NO3– and STC regulate the niche 17
ACS Paragon Plus Environment
Environmental Science & Technology
385
differentiation between nirK- and nirS-type denitrifiers.
386 387
Denitrifier abundance and its direct effect on wetland denitrification rates
388
In the present study, the mean abundances of nirK and nirS genes were 5.22 and
389
2.69×105 gene copies g−1 soil, respectively. In some previous studies conducted in a
390
wide range of wetland habitats, nirS gene in soils generally had a greater abundance
391
compared to bacterial nirK gene.48‒50 However, consistent with our results, Dandie et
392
al. (2011) also reported a dominant role of nirK–type denitrifiers in soil denitrifying
393
communities in a riparian zone.51 Therefore, our findings suggest that nirK‒type
394
denitrifiers are important in wetlands along the Han River. Recent studies have
395
indicated that the abundance of nirK and nirS genes differs seasonally.20 However, we
396
found that there was no significant difference in nirK and nirS gene abundances
397
among seasons (p>0.05; data not shown) and between wetland types (Table S4). Our
398
results support earlier observation by Uksa et al. (2014) who also indicated that the
399
abundance of nirS-type denitrifier did not change significantly over time.52 This
400
suggested that denitrifier abundance do not necessarily covary with the
401
denitrification rates at a temporal scale.
402
As soil denitrification is a microbial–mediated process, the denitrifier abundance
403
may directly regulate the denitrification rate.19,53 So far, only a few studies have
404
investigated the relationships between denitrifier abundance and soil denitrification
405
rates in riparian zones, and their results are not consistent.21,51 Dandie et al. (2011)
406
demonstrated that there was no significant correlation between soil denitrification
407
rates and denitrifier abundances in riparian zones.51 However, Deslippe et al. (2014)
408
reported that the abundance of nirS–type denitrifier was positively related to 18
ACS Paragon Plus Environment
Page 18 of 36
Page 19 of 36
Environmental Science & Technology
409
denitrification potential in riparian soils.21 These seemingly contradictory results are
410
not surprising, because soil denitrification rates reflect the readily activated enzymes
411
but the gene pools only partly contribute to the enzyme activities at a given time.54
412
In addition, the presence of nirK and nirS genes in soils does not necessarily mean
413
that such genes will be expressed or that the protein products of denitrifying genes
414
will function equivalently.53 In the present study, soil denitrification rates increased
415
with increasing abundance of nirK and nirS genes (Table 2), suggesting that the
416
denitrifier gene abundance could be an important limiting factor of soil
417
denitrification and N removal rate in wetlands in the Han River basin.
418
In summary, our study indicated that soil denitrification rates were significantly
419
different between riparian and reservoir shoreline wetlands along the Han River. As
420
we hypothesized, edaphic conditions could significantly affect wetland denitrification
421
rates both directly and indirectly. Directly, soil moisture, NO3–‒N and C availability
422
were positively correlated with soil denitrification rates. Indirectly, edaphic
423
conditions could regulate denitrification rate through their effects on denitrifier
424
abundance. Our findings provide evidence that not only environmental factors, but
425
also biotic factors including denitrifying communities and vegetation, play a critical
426
role in regulating soil denitrification rates. Therefore, future ecological restoration
427
such as revegetation and soil C amendment may effectively enhance the N removal
428
capacity of wetlands along the Han River in the Yangtze River basin.
429 430
Acknowledgments
431
We thank Hui Liu, Xiaoliang Jiang, Zhiyong Zhang and Wenyang Li for their
432
assistance with field sampling, sample processing and laboratory analyses, and four 19
ACS Paragon Plus Environment
Environmental Science & Technology
433
anonymous reviewers for their constructive comments which improved the
434
manuscript. This research was funded in part by the National Natural Science
435
Foundation of China (Grant No. 31270583 and 31570463) and the Key Research
436
Program of the Chinese Academy of Sciences (Grant No. ZDRW-ZS-2016-7).
437
438
Notes
439
The authors declare no competing financial interest
440
441
References
442
(1) Fowler, D.; Coyle, M., Skiba, U.; Sutton, M. A.; Cape, J. N.; Reis, S.; Sheppard,
443
L. J.; Jenkins, A.; Grizzetti, B.; Galloway, J. N. The global nitrogen cycle in the
444
twenty-first century. Philos. Trans. R. Soc. B-Biol. Sci. 2013, 368, 20130164.
445
(2) Carpenter, S. R.; Caraco, N. F.; Correll, D. L.; Howarth, R. W.; Sharpley, A. N.; Smith,
446
V. H. Nonpoint pollution of surface waters with phosphorus and nitrogen. Ecol.
447
Appl. 1998, 8, 1–12.
448 449
(3) Gruber, N.; Galloway, J. N. An earth-system perspective of the global nitrogen cycle. Nature. 2008, 451, 293–296.
450
(4) McCrackin, M. L.; Elser, J. J. Atmospheric nitrogen deposition influences
451
denitrification and nitrous oxide production in lakes. Ecology 2010, 91,
452
528-539.
453
(5) Galloway, J. N.; Aber, J. D. ; Erisman, J. W.; Seitzinger, S. P.; Howarth, R. W.;
454
Cowling, E. B.; Cosby, B. J. The nitrogen cascade. Bioscience. 2003, 53, 341–356.
455
(6) Dodds, W. K. Eutrophication and trophic state in rivers and streams. Limnol.
20
ACS Paragon Plus Environment
Page 20 of 36
Page 21 of 36
Environmental Science & Technology
456 457 458
Oceanogr. 2006, 51, 671–680. (7) Groffman, P. M.; Crawford, M. K. Denitrification potential in urban riparian zones. J. Environ. Qual. 2003, 32, 1144–1149.
459
(8) Bettez, N.D.; Groffman, P.M. Denitrification potential in stormwater control
460
structures and natural riparian zones in an urban landscape. Environ. Sci.
461
Technol. 2012, 46, 10909−10917.
462
(9) Mayer, P. M.; Reynolds, S. K.; McCutchen, M. D.; Canfield, T. J. Meta-analysis of
463
nitrogen removal in riparian buffers. J. Environ. Qual. 2007, 36, 1172–1180.
464
(10) Lowrance, R.; Altier, L. S.; Newbold, J. D.; Schnabel, R. R.; Groffman, P. M.; Denver,
465
J. M.; Correll, D. L.; Gilliam, J. W.; Robinson, J. L.; Brinsfield, R. B.; Staver, K. W.;
466
Lucas, W.; Todd, A. H. Water quality functions of riparian forest buffer systems
467
in Chesapeake Bay Watersheds. Environ. Manage. 1997, 21, 687–712.
468
(11) Findlay, S. E. G.; Mulholland, P. J.; Hamilton, S. K.; Tank, J. L.; Bernot, M. J.; Burgin,
469
A. J.; Crenshaw, C. L.; Dodds, W. K.; Grimm, N.B.; McDowell, W. H.; Potter, J. D.;
470
Sobota, D. J. Cross-stream comparison of substrate-specific denitrification
471
potential. Biogeochemistry 2011, 104, 381–392.
472
(12) Kreiling, R. M.; Richardson, W. B.; Cavanaugh, J. C.; Bartsch, L. A. Summer nitrate
473
uptake and denitrification in an upper Mississippi River backwater lake: the
474
role of rooted aquatic vegetation. Biogeochemistry 2011, 104, 309–324.
475
(13) Hunt, P. G.; Matheny, T. A.; Stone, K. C. Denitrification in a coastal plain riparian
476
zone contiguous to a heavily loaded swine wastewater spray field. J. Environ.
477
Qual. 2004, 33, 2367–2374.
478
(14) Liu, W. Z.; Liu, G. H.; Zhang, Q. F. Influence of vegetation characteristics on soil
479
denitrification in shoreline wetlands of the Danjiangkou Reservoir in China. 21
ACS Paragon Plus Environment
Environmental Science & Technology
480
Clean-Soil Air Water 2011, 39, 109–115.
481
(15) Hefting, M. M.; Bobbink, R.; Janssens, M. P. Spatial variation in denitrification
482
and N2O emission in relation to nitrate removal efficiency in a N-stressed
483
riparian buffer zone. Ecosystems. 2006, 9, 550–563.
484
(16) Tatariw, C.; Chapman, E. L.; Sponseller, R. A.; Mortazavi, B.; Edmonds, J. W.
485
Denitrification in a large river: consideration of geomorphic controls on
486
microbial activity and community structure. Ecology 2013, 94, 2249–2262.
487
(17) Xiong, Z. Q.; Li, S. C.; Yao, L.; Liu, G. H.; Zhang, Q. F.; Liu, W. Z. Topography and
488
land use effects on spatial variability of soil denitrification and related soil
489
properties in riparian wetlands. Ecol. Eng. 2015, 83, 437–443.
490
(18) Morrissey, E. M.; Franklin, R. B. Resource effects on denitrification are mediated
491
by community composition in tidal freshwater wetlands soils. Environ.
492
Microbiol. 2014, 17, 1520–1532.
493
(19) Powell, J. R.; Welsh, A.; Hallin, S. Microbial functional diversity enhances
494
predictive models linking environmental parameters to ecosystem properties.
495
Ecology 2015, 96, 1985–1993.
496
(20) Shrewsbury, L. H.; Smith, J. L.; Huggins, D. R.; Carpenter-Boggs, L.; Reardon, C. L.
497
Denitrifier abundance has a greater influence on denitrification rates at larger
498
landscape scales but is a lesser driver than environmental variables. Soil Biol.
499
Biochem. 2016, 103, 221–231.
500
(21) Deslippe, J. R.; Jamali, H.; Jha, N.; Saggar, S. Denitrifier community size, structure
501
and activity along a gradient of pasture to riparian soils. Soil Biol. Biochem.
502
2014, 71, 48–60.
503
(22) Graham, E. B.; Wieder, W. R.; Leff, J. W.; Weintraub, S. R.; Townsend, A. R.; 22
ACS Paragon Plus Environment
Page 22 of 36
Page 23 of 36
Environmental Science & Technology
504
Cleveland, C. C.; Philippot, L.; Nemergut, D. R. Do we need to understand
505
microbial communities to predict ecosystem function? A comparison of
506
statistical models of nitrogen cycling processes. Soil Biol. Biochem. 2014, 68,
507
279–282.
508 509 510 511
(23) Weisner, S. E. B.; Eriksson, P. G.; Granéli, W.; Leonardson, L. Influence of macrophytes on nitrate removal in wetlands. Ambio 1994, 23, 363–366. (24) Alldred, M.; Baines, S. B. Effects of wetland plants on denitrification rates: a meta-analysis. Ecol. Appl. 2016, 26, 676–685.
512
(25) Liu, W. Z.; Xiong, Z. Q.; Liu, H.; Zhang, Q. F.; Liu, G. H. Catchment agriculture and
513
local environment affecting the soil denitrification potential and nitrous oxide
514
production of riparian zones in the Han River Basin, China. Agric. Ecosyst.
515
Environ. 2016, 216, 147–154.
516
(26) Zhang, Q. F. The South-to-North Water Transfer Project of China: environmental
517
implications and monitoring strategy. J. Am. Water Resour. Assoc. 2009, 45,
518
1238–1247.
519
(27) Li, S. Y.; Liu, W. Z.; Gu, S.; Cheng, X. L.; Xu, Z. F.; Zhang, Q. F. Spatio-temporal
520
dynamics of nutrients in the upper Han River basin, China. J. Hazard. Mater.
521
2009, 162, 1340–1346.
522
(28) Liu, W. Z.; Liu, G. H.; Liu, H.; Song, Y.; Zhang, Q, F. Subtropical reservoir shorelines
523
have reduced plant species and functional richness compared with adjacent
524
riparian wetlands. Environ. Res. Lett. 2013, 8, 044007.
525
(29) Wang, S. H.; Wang, W. W.; Jiang, X.; Zhao, L.; Zhang, B.; Wu, D. W.; Chang, L.
526
Spatial and temporal distribution and flux of nitrogen in water of Danjiangkou
527
Reservoir. Res. Environ. Sci. 2016, 29, 995–1005. (In Chinese with English 23
ACS Paragon Plus Environment
Environmental Science & Technology
528
abstract)
529
(30) Liang, K. X.; Wang, X. Y.; Zhang, D. B.; Zhou, Y. H. Ecological conditions of diatom
530
water bloom formulation in the middle and lower reach of the Hanjiang River
531
and strategy for water bloom control. Environ. Sci. Technol. 2012, 35, 113–116.
532
(In Chinese with English abstract)
533 534
(31) Smith, M. S.; Tiedje, J. M. Phases of denitrification following oxygen depletion in soil. Soil Biol. Biochem. 1979, 11, 261–267.
535
(32) Groffman, P. M.; Altabet, M. A.; Böhlke, J.; Butterbach-Bahl, K.; David, M. B.;
536
Firestone, M. K.; Giblin, A. E.; Kana, T. M.; Nielsen, L. P.; Voytek, M. A. Methods
537
for measuring denitrification: diverse approaches to a difficult problem. Ecol.
538
Appl. 2006, 16, 2091–2122.
539
(33) Bruesewitz, D. A.; Tank, J. L.; Hamilton, S. K. Incorporating spatial variation of
540
nitrification and denitrification rates into whole-lake nitrogen dynamics. J.
541
Geophys. Res. 2012, 117, G00N07.
542
(34) Opdyke, M. R.; David, M. B.; Rhoads, B. L. Influence of geomorphological
543
variability in channel characteristics on sediment denitrification in agricultural
544
streams. J. Environ. Qual. 2006, 35, 2103–2112.
545
(35) McCarty, G. W.; Mookherji, S.; Angier, J. T. Characterization of denitrification
546
activity in zones of groundwater exfiltration within a riparian wetland
547
ecosystem. Biol. Fertil. Soils 2007, 43, 691–698.
548
(36) Dhondt, K.; Boeckx. P.; Hofman. G.; Van Cleemput, O. Temporal and spatial
549
patterns of denitrification enzyme activity and nitrous oxide fluxes in three
550
adjacent vegetated riparian buffer zones. Biol. Fertil. Soils. 2004, 40, 243–251.
551
(37) Hopfensperger, K. N.; Kaushal, S. S.; Findlay, S. E. G.; Cornwell, J. C. Influence of 24
ACS Paragon Plus Environment
Page 24 of 36
Page 25 of 36
Environmental Science & Technology
552
plant communities on denitrification in a tidal freshwater marsh of the
553
Potomac River, United States. J. Environ. Qual. 2009, 38, 618–626.
554 555
(38) Piña-Ochoa, E. ; Álvarez-Cobelas, M. Denitrification in aquatic environments: a cross-system analysis. Biogeochemistry 2006, 81, 111–130.
556
(39) Seitzinger, S. P.; Harrison, J. A.; Böhlke, J. K.; Bouwman, A. F.; Lowrance, R.;
557
Peterson, B.; Tobias, C.; Van Drecht, G. Denitrification across landscapes and
558
waterscapes: a synthesis. Ecol. Appl. 2006, 16, 2064–2090.
559
(40) Song, K. Y.; Kang, H.; Zhang, L.; Mitsch, W. J. Seasonal and spatial variations of
560
denitrification and denitrifying bacterial community structure in created
561
riverine wetlands. Ecol. Eng. 2012, 38, 130–134.
562
(41) Cheng, Y.; Chen, H.; Guo, J. L.; Xu, C. Y. Application of regional climate
563
model-RegCM4 in the simulation of precipitation and temperature in Hanjing
564
Basin. J. Water Res. 2013, 2, 294–300.
565 566
(42) Burgin, A. J.; Groffman, P. M.; Lewis, D. N. Factors regulating denitrification in a riparian wetland. Soil Sci. Soc. Am. J. 2010, 74, 1826–1833.
567
(43) Liu, W. Z.; Wang, Z. X.; Zhang, Q. F.; Cheng, X. L.; Lu, J.; Liu, G. H. Sediment
568
denitrification and nitrous oxide production in Chinese plateau lakes with
569
varying watershed land uses. Biogeochemistry 2015, 123, 379–390.
570
(44) Gift, D. M.; Groffman, P. M.; Kaushal, S. S.; Mayer, P. M. Denitrification potential,
571
root biomass, and organic matter in degraded and restored urban riparian
572
zones. Restor. Ecol. 2010, 18, 113–120.
573
(45) Le Roux, X.; Schmid, B.; Poly, F.; Barnard, R. L.; Niklaus, P. A.; Guillaumaud, N.;
574
Habekost, M.; Oelmann, Y.; Philippot, L.; Salles, J. F.; Schloter, M.; Steinbeiss, S.;
575
Weigelt, A. Soil environmental conditions and microbial build-up mediate the 25
ACS Paragon Plus Environment
Environmental Science & Technology
576
effect of plant diversity on soil nitrifying and denitrifying enzyme activities in
577
temperate grasslands. PLoS One 2013, 8, e61069.
578
(46) Petersen, D. G.; Blazewicz, S. J.; Firestone, M.; Herman, D. J.; Turetsky, M.;
579
Waldrop, M. Abundance of microbial genes associated with nitrogen cycling as
580
indices of biogeochemical process rates across a vegetation gradient in Alaska.
581
Environ. Microbiol. 2012, 14, 993–1008.
582
(47) Guo, G. X.; Deng, H.; Qiao, M.; Yao, H. Y.; and Zhu, Y. G. Effect of long-term
583
wastewater irrigation on potential denitrification and denitrifying communities
584
in soils at the watershed scale. Environ. Sci. Technol. 2013, 47, 3105–3113.
585
(48) Abell, G. C. J.; Revill, A. T.; Smith, C.; Bissett, A. P.; Volkman, J. K.; Robert, S. S.
586
Archaeal ammonia oxidizers and nirS-type denitrifiers dominate sediment
587
nitrifying and denitrifying populations in a subtropical macrotidal estuary. ISME
588
J. 2010, 4, 286–300.
589
(49) Chen, Z.; Luo, X. Q.; Hu, R. G.; Wu, M. N.; Wu, J. S.; Wei, W. X. Impact of
590
long-term fertilization on the composition of denitrifier communities based on
591
nitrite reductase analyses in a paddy soil. Microb. Ecol. 2010, 60, 850–861.
592
(50) Baker, B. H.; Kröger, R.; Brooks, J. P.; Smith, R. K.; Prince Czarnecki, J. M.
593
Investigation of denitrifying microbial communities within an agricultural
594
drainage system fitted with low-grade weirs. Water Res. 2015, 87, 193–201.
595
(51) Dandie, C. E.; Wertz, S.; Leclair, C. L.; Goyer, C.; Burton, D. L.; Patten, C. L.;
596
Zebarth, B. J.; Trevors, J. T. Abundance, diversity and functional gene
597
expression of denitrifier communities in adjacent riparian and agricultural
598
zones. FEMS Microbiol. Ecol. 2011, 77, 69–82.
599
(52) Uksa, M.; Fischer, D.; Welzl, G.; Kautz, T.; Köpke, U.; Schloter, M. Community 26
ACS Paragon Plus Environment
Page 26 of 36
Page 27 of 36
Environmental Science & Technology
600
structure of prokaryotes and their functional potential in subsoils is more
601
affected by spatial heterogeneity than by temporal variations. Soil Biol.
602
Biochem. 2014, 75, 197–201.
603
(53) Wallenstein, M. D.; Myrold, D. D.; Firestone, M.; Voytek, M. Environmental
604
controls on denitrifying communities and denitrification rates: Insights from
605
molecular methods. Ecol. Appl. 2006, 16, 2143–2152.
606
(54) Enwall, K.; Throback, I. N.; Stenberg, M.; Soderstrom, M.; Hallin, S. Soil resources
607
influence spatial patterns of denitrifying communities at scales compatible
608
with land management. Appl. Environ. Microbiol. 2010, 76, 2243−2250.
27
ACS Paragon Plus Environment
Environmental Science & Technology
Page 28 of 36
Table 1 Pearson correlation coefficients between biotic factors (denitrifier abundance and vegetation characteristics) and abiotic factors (edaphic conditions and water quality). nirK abundance 5
‒1
nirS abundance 5
Plant richness
‒1
Plant cover
(10 copies g )
(10 copies g )
-0.45**
-0.25
-0.29**
-0.21*
0.33*
0.18
0.19*
0.24**
-0.12
-0.18
-0.05
0.04
-0.22
-0.20
-0.16
-0.16
Soil NO3 ‒N (mg kg )
0.14
-0.29*
-0.29**
-0.12
STC (mg g‒1)
0.01
-0.29*
-0.14
0.08
-0.20
0.23
0.29**
0.13
0.04
-0.24
-0.05
-0.05
Soil pH Moisture (%) Fine substrate (%) ‒3
Density (g cm ) –
‒1
NH4+‒N (mg L‒1) –
‒1
NO3 ‒N (mg L )
(%)
STC, soil total carbon; * p