Magnitude and Pathways of Increased Nitrous Oxide Emissions from

Natali, S. M.; Schuur, E. A. G.; Mauritz, M.; Schade, J. D.; Celis, G.; Crummer, K. G.; Johnston, C.; Krapek, J.; Pegoraro, E.; Salmon, V. G.; Webb, E...
0 downloads 0 Views 1MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Environmental Processes

Magnitude and pathways of increased nitrous oxide emissions from uplands following permafrost thaw Guibiao Yang, Yunfeng Peng, Maija Marushchak, Yongliang Chen, Guanqin Wang, Fei Li, Dianye Zhang, Jun Wang, Jianchun Yu, Li Liu, Shuqi Qin, Dan Kou, and Yuanhe Yang Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b02271 • Publication Date (Web): 09 Jul 2018 Downloaded from http://pubs.acs.org on July 10, 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 34

Environmental Science & Technology

1

Magnitude and pathways of increased nitrous oxide emissions from uplands

2

following permafrost thaw

3 4

Guibiao Yang1,2, Yunfeng Peng1, Maija E. Marushchak3, Yongliang Chen1, Guanqin

5

Wang1,2, Fei Li1,2, Dianye Zhang1,2, Jun Wang1,2, Jianchun Yu1,2, Li Liu1,2, Shuqi

6

Qin1,2, Dan Kou1,2 and Yuanhe Yang1,2*

7 8

1.

State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China

9 10

2.

University of Chinese Academy of Sciences, Beijing 100049, China

11

3.

Department of Environmental and Biological Sciences, University of Eastern Finland, Kuopio 70211, Finland

12 13 14

*

15

10-62836632, E-mail: [email protected]

Corresponding author: Dr. Yuanhe Yang, tel.: + 86 10-62836638, fax: + 86

16

1

ACS Paragon Plus Environment

Environmental Science & Technology

17

ABSTRACT: Permafrost thawing may release nitrous oxide (N2O) due to large N

18

storage in cold environments. However, N2O emissions from permafrost regions have

19

received little attention so far, particularly with respect to the underlying microbial

20

mechanisms. We examined the magnitude of N2O fluxes following upland

21

thermokarst formation along a 20-year thaw sequence within a thermo-erosion gully

22

in a Tibetan swamp meadow. We also determined the importance of environmental

23

factors and the related-microbial functional gene abundance. Our results showed that

24

permafrost thawing led to a mass release of N2O in recently collapsed sites (3 years

25

ago), particularly in exposed soil patches, which presented post-thaw emission rates

26

equivalent to those from agricultural and tropical soils. In addition to abiotic factors,

27

soil microorganisms exerted significant effects on the variability in the N2O emissions

28

along the thaw sequence and between vegetated and exposed patches. Overall, our

29

results demonstrate that upland thermokarst formation can lead to enhanced N2O

30

emissions, and that the global warming potential (GWP) of N2O at the thermokarst

31

sites can reach 60% of the GWP of CH4 (vs. ~6% in control sites), highlighting the

32

potentially strong non-carbon (C) feedback to climate warming in permafrost regions.

33 34

TOC Art

35 36

Key-words: climate warming; denitrifiers; microbial functional genes; N2O emission;

37

non-carbon (C) feedback; thermokarst; Tibetan Plateau

38

2

ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

Environmental Science & Technology

39

Introduction

40

Permafrost regions account for a large proportion of the global carbon (C) reservoirs,

41

due to the slow microbial decomposition of soil organic matter under cold

42

conditions1-3. The large amount of organic C stored in permafrost soils may represent

43

increasingly relevant sources of methane (CH4) and carbon dioxide (CO2) to the

44

atmosphere under global warming, which is more pronounced at high-latitudes than

45

the global average4-7. The effects of permafrost thaw on gaseous C release have been

46

extensively studied to evaluate the magnitude of permafrost C climate feedback4, 7-12.

47

However, the large amount of nitrogen (N) stored in permafrost soils, which is a

48

potential source for another greenhouse gas, nitrous oxide (N2O)8,

49

ignored. Understanding N2O emissions from thawing permafrost is particularly

50

critical because it has a high global warming potential (GWP), which is ~300-fold

51

greater per unit mass than that of CO2 over a 100-year horizon22, and N2O is also the

52

most important substance leading to ozone depletion in the stratosphere23.

13-21

, is often

53 54

Permafrost thaw, leading to either gradual active layer deepening or often abrupt

55

development of thermokarst landforms24, may increase N2O emissions due to the

56

unlocking of vast N stocks stored in frozen soils. However, current experimental

57

evidence are mainly derived from lowland permafrost ecosystems14,

58

observations are available for upland ecosystems. Upland permafrost landscapes

59

prone to thermokarst formation cover as much as 0.91×106 km2, or ~5% of the

60

permafrost regions24 and belong to an important category of permafrost thawing 3

ACS Paragon Plus Environment

17

, and limited

Environmental Science & Technology

Page 4 of 34

61

around the world. More importantly, compared with lowlands which tend to be

62

waterlogged following thermokarst formation25, thermokarst subsidence in upland

63

permafrost usually results in intermediate or even lowered soil moisture26, 27 and the

64

formation of adjacent aerobic and anaerobic soil environments8. These intermediate

65

soil moisture conditions may promote coupled nitrification (oxidation of ammonia to

66

NO3-) and denitrification (anaerobic respiration using NO3- as an electron acceptor),

67

which may lead to high N2O production and emissions15, 28. Moreover, disturbances of

68

the vegetation cover may also promote N2O emissions via improved nutrient

69

availability for soil microbes in the absence of plant uptake15, 17. In support of these

70

deductions, elevated N2O concentrations in the soil have been observed in upland

71

thermokarst features8. However, direct field evidence of post-thaw N2O emissions in

72

upland thermokarst remains unavailable.

73 74

The escape of N2O from soils to the atmosphere after permafrost thaw is the result of

75

the balance of multistep processes related to nitrification and denitrification29, and

76

these processes are primarily driven by microbial pathways30. Multiple specialized

77

microbial groups are responsible for nitrification-mediated (ammonia oxidizing

78

archaea

79

microorganisms). Further, the denitrification pathway is dependent on nitrification for

80

the supply of nitrate, the initial electron acceptor for denitrification. Soil abiotic

81

properties indirectly influence N2O flux via regulating microbial abundance and their

82

activities31. Although the importance of plant and abiotic factors in the regulation of

and

bacteria)

and

denitrification-mediated

4

ACS Paragon Plus Environment

pathways

(denitrifying

Page 5 of 34

Environmental Science & Technology

83

N2O fluxes has been widely examined in permafrost areas (e.g., absence of plant N

84

uptake, favorable moisture content and pH, low soil C:N ratios and relatively high

85

inorganic N concentrations)13-15,

86

directly associated with N2O production under permafrost thaw remain lacking19.

17

, investigations of the microbial communities

87 88

The Tibetan permafrost region, covering 67%, or ~1,500,000 km2, of the plateau

89

area32, has experienced significant climate warming and consequential permafrost

90

thawing over the past decades. The alpine swamp meadows, which contains a high

91

soil N density similar to that in the northern circumpolar permafrost region (Fig. S1),

92

offer an opportunity to examine the changes in N2O emissions following permafrost

93

collapse together with the underlying mechanisms. Here, we measured the in situ N2O

94

flux along a thermo-erosion gully in an alpine swamp meadow underlain by

95

permafrost. The abiotic factors and microbial functional gene abundance were

96

determined to gain a mechanistic understanding of how thermokarst formation affects

97

N2O emissions. We hypothesize that 1) upland thermokarst formation enhances N2O

98

emissions because of the mobilization of permafrost N and reduced plant N uptake,

99

co-occurring with favorable moisture content, and that 2) thaw-induced changes in

100

N2O flux are controlled by interactions between abiotic factors and soil

101

microorganisms.

102 103

Materials and Methods

104

Experimental site 5

ACS Paragon Plus Environment

Environmental Science & Technology

105

The study site is an alpine swamp meadow underlain by discontinuous permafrost,

106

located on the northeastern Tibetan Plateau, China (N 37°28′, E 100°17′, ~3900 m

107

a.s.l; Fig. S2)33. The region is characterized by a wet alpine climate with a mean

108

annual air temperature of -3.3°C and a mean annual precipitation of 460 mm. The

109

active layer, which is affected by seasonal thawing, is ca. 90 cm thick as measured by

110

a GSSI georadar (SIR-20, Laurel, Santa Clara, USA) and validated by manual

111

measurements by a thaw probe. Vegetation at the site is mainly dominated by

112

graminoids, such as Kobresia tibetica, K. royleana and Carex atrofuscoides. The soil

113

in this area has a silty loam texture with 2.0% clay, 53.7% silt, and 44.3% sand. The

114

chemical properties of the 0–15 cm soil layer are as follows: soil organic C, 19.4%;

115

total N, 1.5%; bulk density, 0.27 g/cm3; and pH, 5.6.

116 117

We measured the N2O flux along a thermo-erosion gully during two summers in 2015

118

and 2016. The gully occurred in the early 1990s according to local residents and was

119

~240 m long and approximately 2 m at the maximum depth at the time of

120

measurements. In May 2015, four sites with an area of 15×20 m2 were established

121

outside and inside the gully. The control site was an undisturbed area outside the

122

thermo-erosion gully. The other three sites were dispersed at intervals of ~80 m from

123

the head to the end along the gully, representing different thawing stages (early stage:

124

collapse occurred 3 years ago; mid-stage: collapse occurred 12 years ago; and late

125

stage: collapse occurred 20 years ago).

126 6

ACS Paragon Plus Environment

Page 6 of 34

Page 7 of 34

Environmental Science & Technology

127

N2O flux measurements

128

Fluxes of N2O were determined using static chamber and gas chromatography

129

techniques34. In 2015, ten vegetated patches were randomly chosen at the control site

130

and at each of the three sites within the gully. A PVC ring (26 cm in diameter, 12 cm

131

in height) with a collar was inserted into the soil to a depth of 10 cm. Because

132

thermokarst formation often creates a mosaic of vegetated and exposed patches8 with

133

potentially different N2O emission rate13,

134

exposed patches inside the gully in 2016. The N2O flux was measured from 5

135

positions in the exposed patches in each site due to their smaller proportion of the

136

feature area. During gas collection, an opaque cylindrical chamber (25 cm tall) was

137

placed on the collar. The chamber was covered with insulation materials to avoid a

138

rapid increase in chamber temperature and equipped with an electric fan for proper

139

mixing of the headspace air. Gas samples (50 ml) were collected from the headspace

140

at 10-min intervals (0, 10, 20, 30, and 40 min) and immediately injected into

141

pre-evacuated air bags (Delin Inc., Dalian, Liaoning, China). Gas samples were then

142

transported to the laboratory and analyzed for N2O concentration using a gas

143

chromatograph (Agilent 7890A, Agilent Technologies Inc., Santa Clara, California,

144

USA; for analytical methodologies, see Wang & Wang, 2003). The rate of N2O flux

145

was calculated from the slope of the N2O concentration–time function. Gas samples

146

were collected between 9:00 am and 12:00 pm; in 2015, ca. every 10 days in July and

147

August, and in 2016 once in early July, thrice during consecutive days at the end of

148

July, once at the end of August. In total, 2275 gas samples were analyzed to obtain the

17

, we also measured the N2O flux in

7

ACS Paragon Plus Environment

Environmental Science & Technology

149

N2O fluxes. At the time of gas collection, the soil temperature and volumetric water

150

content (VWC) at a 0–10 cm soil depth were measured in each patch adjacent to the

151

PVC ring using a digital thermometer and a portable TDR-100 soil moisture probe

152

(Spectrum Technologies Inc., Plainfield, IL, USA).

153 154

Soil sampling and biochemical analyses

155

To discern the mechanism underlying the thermokarst effects on N2O emissions, we

156

collected soil samples from the four sites and measured the soil physico-chemical

157

characteristics and functional gene abundance of denitrifying microbes. Six soil cores

158

(2.5 cm in diameter and 15 cm in depth) were collected near each PVC ring and

159

homogenized to gain one composite sample during the intensive flux measurement

160

campaign at the end of July 2016. In total, 240 soil cores were collected (40 soil

161

samples). The soil samples were passed through a 2-mm-mesh sieve after removal of

162

roots and divided into three parts. The first part was stored at 4°C to determine the soil

163

nitrate (NO3 ) concentration. The second part was air-dried to measure the soil pH, and

164

the third part was stored at -80°C to quantify the key microbial functional gene

165

abundances. The soil NO3- concentration was determined using a flow injection

166

analyzer (Autoanalyzer 3 SEAL; Bran and Luebbe, Norderstedt, Germany) after

167

extracting fresh soil with 2 M KCl at a soil to solution ratio of 1:5 (v/v). Soil pH was

168

measured in a soil water mixture (dry soil to deionized water ratio of 1:5 (v/v)) with a

169

pH electrode (PB-10, Sartorius, Germany).

-

170 8

ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

Environmental Science & Technology

171

For the soil microbial functional gene quantification, the present study only focused on

172

denitrifying microorganisms because the soil water pore space (WFPS) was between

173

62.6% and 84.1% across treatments and growing seasons. Under this condition,

174

microbial denitrification is considered the predominant pathway of N2O production35.

175

Here, we quantified the abundance of genes that encoded nitrite reductase (nirK and

176

nirS) and nitrous oxide reductase (nosZ). For the analysis, DNA was extracted from

177

wet soil (0.4 g) using the Power Soil® DNA Isolation Kit (MoBio Laboratories,

178

Carlsbad, CA, USA) according to the manufacturers’ instructions. The DNA quality

179

was checked with a NanoDrop-2000 (Thermo Fisher Scientific, Madison, WI, USA),

180

which was assessed based on the absorbance values at 260/280 nm and 260/230 nm. A

181

quantitative real-time polymerase chain reaction (qPCR) analysis was performed in

182

triplicate on a StepOne Plus Real-Time PCR system (Applied Biosystems, Inc., CA,

183

USA) in 20 µl of reaction mixture containing 10 µl of SYBR premix Ex Taq (Tli

184

RNaseH Plus), 0.4 µl forward and reverse primers, 0.4 µl of ROX Reference Dye

185

(50×), 6.8 µl of DNA-free water and 2 µl of fivefold diluted template DNA. The

186

amplification conditions, reading temperature, primer pairs, number of cycles and

187

references used to quantify the abundance of the nirK, nirS and nosZ genes are

188

summarized in Table S1. For quantification, standard curves were constructed with a

189

known concentration of purified plasmid DNA carrying the cloned portion of the

190

corresponding gene with r2 > 0.95. The abundances of the nirK, nirS and nosZ genes

191

are presented as a copy number per gram of dry weight (copies g-1 dw).

192 9

ACS Paragon Plus Environment

Environmental Science & Technology

193

Statistical analysis

194

All data were checked for normality before the analysis. The N2O flux data were

195

log-transformed to meet the assumption of a normal distribution. To assess the effects

196

of thermokarst formation on N2O fluxes, a repeated-measures ANOVA was used to

197

compare the differences in the N2O flux among the thawing stages and between

198

vegetated and exposed patches, with the thawing stage or patch type as the

199

between-subject factor, and the sampling date as the within-subject factor. One-way

200

ANOVA were employed to evaluate the differences in seasonal means of the N2O flux,

201

soil temperature and VWC as well as the soil pH, NO3 concentration, and the

202

microbial functional gene abundance associated with denitrification (nirK, nirS, nosZ)

203

among the thawing stages and between vegetated and exposed patches. We used post

204

hoc tests (Tukey’s HSD) to test whether the above mentioned parameters exhibited

205

significance at the level of P