Subscriber access provided by University of Sunderland
Characterization of Natural and Affected Environments
Contrasting warming and ozone effects on denitrifiers dominate soil N2O emissions Yunpeng Qiu, Yu Jiang, Lijin Guo, Kent O. Burkey, Richard W. Zobel, H. David Shew, and Shuijin Hu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b01093 • Publication Date (Web): 29 Aug 2018 Downloaded from http://pubs.acs.org on August 30, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 41
1
Environmental Science & Technology
Contrasting warming and ozone effects on denitrifiers dominate soil N2O emissions
2 3
Yunpeng Qiu†,//, Yu Jiang†,⊥, Lijin Guo†, $, Kent O. Burkey#, &, Richard W. Zobel&, H. David
4
Shew†, and Shuijin Hu†,*
5 6
†
7
27695, USA
8
//
9
210095, China
Department of Entomology & Plant Pathology, North Carolina State University, Raleigh, NC
College of Resources and Environmental Sciences, Nanjing Agricultural University, Nanjing
10
⊥
11
$
12
China
13
#
14
&
15
USA
Institute of Applied Ecology, Nanjing Agricultural University, Nanjing 210095, China College of Plant Science & Technology, Huazhong Agricultural University, Wuhan 430070,
USDA-ARS, Plant Sciences Research Unit, Raleigh, NC 27607, USA Department of Crop and Soil Sciences, North Carolina State University, Raleigh, NC 27695,
16 17 18 19 20 21 22 23
ACS Paragon Plus Environment
Environmental Science & Technology
24
Abstract
25 26
Nitrous oxide (N2O) in the atmosphere is a major greenhouse gas and reacts with volatile organic
27
compounds to create ozone (an air pollutant) in the troposphere. Climate change factors such as
28
warming and elevated ozone (eO3) affect N2O fluxes, but the direction and magnitude of these
29
effects are uncertain and the underlying mechanisms remain unclear. We examined the impact of
30
simulated warming (control + 3.6 oC) and eO3 (control + 45 ppb) on soil N2O fluxes in a soybean
31
agroecosystem. Results obtained showed that warming significantly increased soil labile C,
32
microbial biomass and soil N mineralization, but eO3 reduced these parameters. Warming
33
enhanced N2O-producing denitrifers (nirS- and nirK-type), corresponding to increases in both
34
the rate and sum of N2O emissions. In contrast, eO3 significantly reduced both N2O-producing
35
and N2O-consuming (nosZ-type) denitrifiers but had no impact on N2O emissions. Further, eO3
36
offsets the effects of warming on soil labile C, microbial biomass and the population size of
37
denitrifiers, but still increased N2O emissions, indicating a direct effect of temperature on N2O
38
emissions. Together, these findings suggest that warming may promote N2O production through
39
increasing both the abundance and activities of N2O-producing microbes, positively feeding back
40
to the ongoing climate change.
41 42 43 44 45 46
ACS Paragon Plus Environment
Page 2 of 41
Page 3 of 41
47
Environmental Science & Technology
TOC/Abstract art
48
49 50 51 52 53 54 55 56 57 58 59 60
ACS Paragon Plus Environment
Environmental Science & Technology
61
INTRODUCTION
62 63
N2O is a potent greenhouse gas with a global warming potential 300 times greater than that of
64
carbon dioxide (CO2)1,2 and is the single most important depleter of stratospheric ozone (O3).3 In
65
the troposphere, however, N2O reacts with volatile organic compounds to create O3, an air
66
pollutant that also contributes to global warming.4 Since industrial times, human activities have
67
increased global atmospheric O3 concentration.5,6 Various models predict that global surface
68
temperature will be increased by 1-6 oC by the end of this century.1,7 Many regions of the globe
69
have experienced surface O3 concentration greater than 40 nmol mol-1 (i.e., 40 ppb), a threshold
70
set for injury to O3-sensitive plants.6,8 The current global atmospheric ozone concentration has
71
increased from an estimated preindustrial level of 10-15 nmol mol-1 to current ca. 60 nmol mol-
72
19
73
and chemical N fertilizer applications.1,10 Soil N2O is mainly produced by two microbial
74
processes: nitrification (i.e., oxidation of NH4+ into NO3-) and denitrification (i.e., reduction of
75
NO3- to gaseous forms, N2O, NO and N2).11-13
76
. Globally, agricultural soils account for nearly 70% of the annual N2O budget owing to manure
Climate change factors such as elevated CO2 (eCO2), warming and elevated O3 (eO3) may
77
significantly impact N2O emissions.14-16 While eCO2 tended to increase N2O emissions17,18, the
78
impact of warming was more variable, with positive,19 negative20 or even neutral effects21 having
79
been reported under different conditions.22-24 These inconsistent results may stem from the fact
80
that warming impacts multiple processes, some of which may have contrasting effects on soil
81
N2O production and consumption.14,21 On one hand, warming may increase N2O emissions
82
through directly stimulating the growth and activities of N2O-producing microbes19,25 and/or
83
indirectly enhancing the availability of C via plant growth,21,26 and mineral N via mineralization
ACS Paragon Plus Environment
Page 4 of 41
Page 5 of 41
Environmental Science & Technology
84
(Figure 1).14,27 On the other hand, warming may enhance plant growth, evapotranspiration and N
85
uptake, and reduce soil available N21,28 and soil moisture, thereby constraining nitrifiers and
86
denitrifiers and their activities.29-31 In contrasting to the direct effect of warming, eO3 affects
87
nitrifiers and denitrifiers, and their activities mainly through its indirect effect on plants (Figure
88
1).32,33 This is because O3 itself hardly penetrates into soil to directly affect soil microbes.34
89
Rather, eO3 often reduces plant growth, plant N acquisition and photosynthate allocation
90
belowground for soil microbes,35-37 altering soil N and/or labile C availability for nitrifiers or
91
denitrifiers (Figure 1).38-40 However, only a few experiments have so far examined the effect of
92
eO3 on N2O emissions.32,39 Negative O3 effect on N2O has been reported in a meadow
93
ecosystem41 and in rice fields with the presence42 or the absence of roots,43 suggesting that
94
reduced denitrification resulting from O3-induced decreases in plant-derived C inputs may play a
95
dominant role.
96
An increasing number of experiments have recently examined the impact of warming on
97
microbial activities and N2O emission.19,30,44 However, climate warming will likely concur with
98
other components such as N deposition and eO3 under future climate change scenarios. Yet, few
99
have examined the combined effect of warming and eO3 on N2O emissions. Even fewer have
100
attempted to establish the linkages between climate change-induced alteration in N2O emission
101
and populations of N2O-producing microbes in field, partially due to the difficulty in quantifying
102
N2O–producing microbes.
103
Both nitrifiers and denitrifiers contribute to N2O emissions from soil. Ammonia-
104
oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA) regulate the rate-control step
105
of nitrification that produces N2O under aerobic conditions45,46 and their amoA gene that encodes
106
the ammonia monooxygenase (AMO) has been effectively used to quantify their abundances in
ACS Paragon Plus Environment
Environmental Science & Technology
107
soil.47-49 Although denitrifiers are highly diverse heterotrophs and is currently impossible to
108
quantify their species composition,50 the function genes encoding nitrite reductase (Nir), the key
109
enzyme in the dissimilatory denitrification, have been widely used to characterize the
110
abundances and diversity of denitrifiers.50-52 Genes nirK and nirS encode two structurally
111
different enzymes catalyzing nitrite reduction: the copper and cytochrome cd1-nitrite reductase,
112
respectively.51,53 Another gene, nosZ, which encodes nitrous oxide reductase to convert N2O to
113
N2 has also been widely used to describe denitrifier communities.50,54
114
We initiated a long-term field experiment in a soybean agroecosystem to understand how
115
warming, elevated ozone and their interactions affect plant growth, and soil C and N cycling at
116
Raleigh, NC, USA. Taking advantage of this field study, we examined the impact of warming
117
and ozone on N2O fluxes, soil labile C and N pools, and the abundances of nitrifiers and
118
denitrifiers that are primarily responsible for N2O emissions. Our hypotheses (Figure 1) were that
119
(1) warming increases the abundances of N2O-producing microbes and N2O emissions mainly
120
through enhancing labile C and mineral N sources rather than through its direct impact on the
121
N2O-producing microbes, (2) eO3 reduces the abundances of denitrifiers and N2O emissions by
122
reducing organic C inputs, (3) the effects of warming and eO3 on N2O-producing microbes offset
123
each other.
124 125
MATERIALS AND METHORDS
126 127
The experimental site
128
The long-term experiment was initiated in June 2013 to investigate the response of soybean
129
plants to elevated atmospheric O3 and climate warming using an air exclusion system (AES, 3m
ACS Paragon Plus Environment
Page 6 of 41
Page 7 of 41
Environmental Science & Technology
130
× 10m). The experimental site is located at the Lake Wheeler Experimental Station, 5 km south
131
of North Carolina State University, Raleigh, North Carolina, USA (35。43ʹN, 78。40ʹW; elevation
132
120 m). The soil is an Appling sandy loam (fine, kaolinitic, thermic Typic Kanhapludult), well
133
drained with a pH of 6.2, and contained 12.5 g C and 0.94 g N kg-1 soil when the experiment
134
started. Before O3 and warming treatments were initiated, the soil was repeatedly turned-over
135
using a disc implement and rotovator.
136
Each AES plot consisted of two double-walled open-top chamber panels55 attached to
137
aluminum frames and placed in parallel to establish a 3 m x 10 m footprint. A continuous flow of
138
air was supplied through the holes in the internal wall of each panel by a fan box containing an
139
activated charcoal filter to reduce the O3 concentration of incoming ambient air to ca. 15 nmol
140
mol-1. In elevated temperature plots, each fan box was equipped with two 4000W electrical
141
resistance heaters (a total of 16000W per plot) and an inline water-to-air heat exchanger. The
142
electrical resistance heaters operated 24 hours per day and solar-heated water circulated through
143
the heat exchanger during daylight hours. We designed a warming treatment with an increase on
144
the ambient air temperature by 3.6 oC. We also selected an elevated ozone concentration at 60
145
nmol mol-1, a value close to the average and ozone concentration in our area. We did not increase
146
the O3 concentration because the ambient O3 already has harmful effects on some soybean
147
cultivars. Moisture was added to the heated air using misters inside the double walled panel to
148
elevate the relative humidity to near ambient levels prior to distribution within the experimental
149
plot. O3 was generated from O2 by corona discharge (model TG-20, Ozone Solutions, Hull, IA,
150
USA). The O3 concentration in each plot was monitored at canopy height with a UV photometric
151
O3 analyzer (model 49C, Thermo Environmental Instruments Co., Franklin, MA, USA). Air
152
temperature and relative humidity were monitored using HOBO sensors (model U23-001,
ACS Paragon Plus Environment
Environmental Science & Technology
153
HOBO Pro V2 Temp/RH Datalogger, Onset Computer Corporation, Bourne, MA, USA). Soil
154
temperature was measured with an Enviropro dielectric soil probe (model EP100D, Apcos Pty
155
Ltd, Adelaide, Australia).
156
This experiment was a 2×2 factorial design with four treatments assigned into three
157
blocks (replicates). Four different treatments were: (a) charcoal-filtered air and ambient
158
temperature (CF); (b) charcoal-filtered air plus 3.6 oC increase in temperature (T); (c) charcoal-
159
filtered air plus ambient temperature and 1.4 times ambient ozone (O3);38,39 (d) charcoal-filtered
160
air + warming + 1.4 times ambient ozone (T+O3). The purpose of filtration of ambient air with
161
activated charcoal was to reduce ambient O3 concentrations to levels considered nonphytotoxic
162
to soybeans. The seasonal daily average concentrations (12-h) of O3 and degree of temperature
163
over the experimental duration are shown in Table 1.
164
Soybean (cultivar Jake) was selected based on our preliminary greenhouse experiment in
165
2015 showing that 60 ppb resulted in a moderate foliar ozone response for this cultivar. Soybean
166
seeds were planted in the first week of June, 2015 into AES plots. Each plot consisted of eight
167
plant rows with a row spacing of 25 cm. During soybean growing seasons, plants were irrigated
168
with drip lines to prevent water stress. There were in total three harvesting of aboveground
169
biomass in August, September and October and two harvesting of root biomass in August and
170
September each year. Aboveground biomass was harvested at soil surface. Three soil cores of 5
171
cm in diameter were taken to a 20-cm depth in each plot to harvest three plant roots. Roots were
172
washed with tap water. Aboveground and root biomass were weighed after oven drying at 70 °C
173
for 72 h.
174 175
Field N2O measurements
ACS Paragon Plus Environment
Page 8 of 41
Page 9 of 41
Environmental Science & Technology
176
Each AES plot has 8 rows with two static gas chambers randomly installed into all the plots. Gas
177
samples were taken within 24-48 hours following an irrigation event approximately every 10
178
days. In total, we have 13 gas measurements through the growing season (June to October). The
179
chambers (sampling volume is 6392 ml) were closed with a lid, and fitted with a rubber septum
180
to allow gas sampling via syringe. Gas samples (5 ml) were taken with 20 ml PE syringes
181
(Beceton Dickinson Franklin Lakes, NJ, USA) and immediately injected into N2-preflushed 12
182
mL vials at time 0 and 30 minute after the closure of lids. Air temperature and soil moisture were
183
also recorded during each sampling period. Gas samples were analyzed within 24 h on a gas
184
chromatograph with an autosampler (Shimadzu AOC-5000 Auto-Injector). Two detectors, a
185
flame ionization detector (FID) (Shimadzu GC – 2014, Kyoto, Japan) and an electron capture
186
detector (ECD), were fitted with the GC so that CO2, CH4 and N2O can de simultaneously
187
measured. N2O was detected by ECD with a detection limit of 0.1 ppm. The N2O fluxes and
188
cumulative N2O emissions were calculated using formulas described by Wu et al. (2017).18
189 190
Soil sampling
191
We used a 5-cm diameter soil corer to take six to seven soil cores to 15 cm depth after each
192
harvesting. Soil cores were separated into 0-7.5 cm and 7.5-15 cm depth fractions. All soil
193
samples were sealed in plastic bags, stored in a cooler and transported to the laboratory. In total,
194
there are three soil sampling happened in August, September and October for each plot.
195
Soil samples were mixed thoroughly and sieved through a 2-mm mesh within 24 h of the
196
field sampling and we carefully removed all visible residues and plant roots. Subsamples (~ 20
197
g) were then taken immediately and stored at -20 °C for DNA analysis and the remainder of soils
ACS Paragon Plus Environment
Environmental Science & Technology
198
were stored at 4 °C for other microbial and chemical analyses. A 10-g soil subsample was oven-
199
dried at 105 °C for 48 h and weighed for the determination of the water content.
200 201
Chemical and microbial analyses of soil samples
202
Soil microbial biomass C and N. Soil microbial biomass C (MBC) and biomass N (MBN) were
203
measured using the fumigation-extraction method.56 For each field soil sample, 20-g subsample
204
(dry soil equivalent) was fumigated with ethanol-free chloroform for 48 h and then extracted
205
with 50 mL of 0.5 M K2SO4 by shaking for 30 min. Another 20-g subsample of nonfumigated
206
soil was extracted with 50 mL of 0.5 M K2SO4 after shaking for another 30 min. Soil extractable
207
organic C in both fumigated and non-fumigated extracts was determined using a TOC analyzer
208
(Shimadzu TOC-5050A, Shimadzu Co., Kyoto, Japan). Soluble inorganic N in the extracts of
209
fumigated and non-fumigated soils was quantified on the Lachat flow injection analyzer (Lachat
210
Instruments, Milwaukee, WI, USA) after digestion with alkaline persulfate.57 The differences in
211
extractable organic C and inorganic N between fumigated and non-fumigated soils were used to
212
calculate MBC and MBN using a conversion factor of 0.33 (kEC) and 0.45 (kEN), respectively.56,58
213
Soil extractable C and N. The concentration of soluble organic C in non-fumigated soil extracts
214
was used to represent soil extractable C. The extractable inorganic N referred to the sum of NH4+
215
and NO3- in non-fumigated soil extracts.59
216
Soil microbial respiration. An incubation-alkaline method was used to determine soil
217
heterotrophic respiration.60 Briefly, a subsample of 20.0g soil (dry soil equivalent) was put into a
218
125-mL beaker and adjusted to a moisture level of 60% water holding capacity, and the beaker
219
was then placed in a 1-L Mason jar. 5.0 mL of 0.25 M fresh NaOH were added to trap respired
220
CO2 in a 50-mL beaker that was suspended in each mason jar.61 All mason jars were incubated at
ACS Paragon Plus Environment
Page 10 of 41
Page 11 of 41
Environmental Science & Technology
221
25 °C in the dark for 2 weeks. After the first week incubation, NaOH solutions were replaced
222
with fresh ones. The respired CO2 from the soil was titrated with 0.125 M HCl to determine the
223
soil microbial respiration (SMR). SMR rate was expressed as mg CO2 kg-1 soil d-1 by averaging
224
the data across two1-wk incubations.59
225
Net soil N mineralization. Potential N mineralization was determined following a 4-wk
226
incubation at 25 °C in the dark. After incubation, soil NH4+ and NO3- were extracted with 50 mL
227
of 0.5 M K2SO4 following shaking for 30 min. Inorganic N concentrations were also measured
228
on the Lachat flow injection analyzer. Net mineralized N (NMN) was then calculated by the
229
difference in extractable total inorganic N between un-incubated and incubated soil subsamples.
230
Abundances of functional genes responsible for nitrification and denitrification in soil. Fresh soil
231
(0.50 g) from final harvest was extracted for nucleic acid using FastDNA SPIN kit (MP Bio,
232
Solon, OH, USA) according to the manufacturer’s instructions. Soil DNA quality and size were
233
checked by electrophoresis on a 1% agarose gel. The copy numbers of AOA amoA, AOB amoA,
234
nirS, nirK and nosZ genes of each soil DNA sample were determined using quantitative real-time
235
polymerase chain reaction (PCR) (CFX96 Real-Time PCR Detection System, Bio-Rad,
236
Hercules, CA, USA) with the primers given in SI Table S1. The final qPCR reaction volume was
237
20 µl, with a of 2 µl of template DNA, 14 µl 1× SsoAdvanced TM SYBR Green Supermix (Bio-
238
Rad, Hercules, CA, USA), and 2 µl of each primer. The standard curve for determining the gene
239
copy number was developed with the agarose gel-purified PCR products based upon the method
240
of Chen et al. (2015).62 All the qPCR reactions followed the conditions were shown in SI Table
241
S1.
242 243
Statistical analysis
ACS Paragon Plus Environment
Environmental Science & Technology
244
We examined results for the entire growing season in 2015 and used the linear mixed model63 to
245
test the main effects of temperature, O3 and the interaction of temperature and O3, and whether
246
these changed over time. We employed a set of covariance structures including compound
247
symmetric model (CS), the first-order autoregressive model [AR (1)], and autoregressive with
248
random effect to reduce autocorrelation. The P values for treatments and interaction terms were
249
reported based on the covariance structure that minimized Akaike information criterion (AIC)
250
and Bayesian information criterion (BIC).63 Data for cumulative N2O emissions and gene copy
251
numbers were subjected to the analysis of variance using the mixed model. All statistical
252
analyses were performed using the SAS 9.4 (SAS Institute, Inc., Cary, NC, USA). For all tests,
253
the significant differences were determined at the 95% probability level.
254 255
RESULTS
256 257
Characteristics of climate change treatments
258
During the experiment period, the difference in mean daily temperature between control and
259
elevated temperature treatments (T and T+O3) was 3.6 ± 0.07 oC (Table 1 and SI Figure S1). The
260
mean daily (12-h) difference in O3 concentration between the control and elevated ozone
261
treatments (O3 and T+O3) was 45.6 ± 1.0 ppb (Table 1 and SI Figure S1). The increases in both
262
temperature and O3 well met the respective manipulation targets. Meteorological variables (i.e.,
263
air humidity, soil moisture and soil temperature) recorded on days of sampling gas indicated
264
lower air humidity and higher soil temperature in the warming treatments (T and T+O3)
265
compared with the control (CF). No significant differences in soil moisture were found for the
266
sampling dates during the experimental study (Table 1).
ACS Paragon Plus Environment
Page 12 of 41
Page 13 of 41
Environmental Science & Technology
267 268
Plant biomass
269
Warming increased aboveground plant growth at the early stage but this effect faded then after
270
(Figure 2a). Across the three sampling dates, warming (T) alone had no significant effect on
271
aboveground biomass (Figure 2a). However, it significantly increased root biomass by 22%
272
(Figure 2b). In comparison, eO3 had no effect on aboveground biomass (Figure 2a), but tended to
273
reduce belowground biomass at the final harvest (-11%) (Figure 2b). There were significant
274
effects of the warming × O3 interaction on both aboveground (Figure 2a) and belowground
275
biomass (Figure 2b). Warming and eO3 combination (T+O3) aggravated the warming and ozone
276
stresses on aboveground plant growth and offset the warming effect on root growth, reducing the
277
biomass at the final harvest by 23% and 13%, respectively (Figure 2).
278 279
Soil extractable organic C and N
280
Warming had no significant effect on soil extractable organic C in the surface soil (0-7.5 cm),
281
but significantly enhanced it in the subsoil (7.5-15 cm) (SI Figure S2 and Table S2). eO3 reduced
282
soil extractable organic C (SI Figure S2 and Table S2). Compared to the low ozone plots (CF, T),
283
soil extractable organic C under elevated ozone (O3, T+O3) was 12% (P