High methane emissions largely attributed to ebullitive fluxes from a

Mar 13, 2019 - Over the two-year period, annual CH4 emissions averaged 29.52 mmol m-2 d-1, amounting to 10.78 mol m-2 yr-1, making the river a strong ...
1 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF LOUISIANA

Environmental Processes

High methane emissions largely attributed to ebullitive fluxes from a subtropical river draining a rice paddy watershed in China Shuang Wu, Shuqing Li, Ziheng Zou, Tao Hu, Zhiqiang Hu, Shuwei Liu, and Jianwen Zou Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.8b05286 • Publication Date (Web): 13 Mar 2019 Downloaded from http://pubs.acs.org on March 16, 2019

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 31

Environmental Science & Technology

1

Research Article

2

Number of tables: 1

3

Number of figures: 4

4

Number of text pages: 30

5

Number of references: 51

6 7

High methane emissions largely attributed to ebullitive fluxes from a

8

subtropical river draining a rice paddy watershed in China

9

Author list:

10

Shuang Wu1,a, Shuqing Li1,b, Ziheng Zouc, Tao Hua, Zhiqiang Hud, Shuwei Liua,b* and Jianwen Zoua,b

11

Affiliations:

12

a

13

Environmental Sciences, Nanjing Agricultural University, Nanjing, China; bJiangsu Key Lab and

14

Engineering Center for Solid Organic Waste Utilization, Jiangsu Collaborative Innovation Center for

15

Solid Organic Waste Resource Utilization, Nanjing Agricultural University, Nanjing, China; cCollege

16

of Overseas Education, Nanjing Tech University, Nanjing, China; dTaizhou University, Taizhou, China.

Jiangsu Key Laboratory of Low Carbon Agriculture and GHGs Mitigation, College of Resources and

17 18

1

These authors contributed equally to this work

19 20

Keywords: Agricultural impacted watershed; Methane; Diffusive flux; Ebullitive flux;

21

Greenhous gas; Subtropical river

22 23

To be submitted to Environmental Science & Technology

24

(Revised version, manuscript ID: es-2018-052869) 1

ACS Paragon Plus Environment

Environmental Science & Technology

25

ABSTRACT

26

Rivers are of increasing concern as sources of atmospheric methane (CH4), while

27

estimates of global CH4 emissions from rivers are poorly constrained due to a lack of

28

representative measurements in tropical and subtropical latitudes. Measurements of

29

complete CH4 flux components from subtropical rivers draining agricultural watersheds

30

are particularly important since these rivers are subject to large organic and nutrient loads.

31

Two-year measurements of CH4 fluxes were taken to assess the magnitude of CH4

32

emissions from the Lixiahe River (a tributary of the Grand Canal) draining a subtropical

33

rice paddy watershed in China. Over the two-year period, annual CH4 emissions averaged

34

29.52 mmol m-2 d-1, amounting to 10.78 mol m-2 yr-1, making the river a strong source of

35

atmospheric CH4. The CH4 emissions from rivers during the rice-growing season (June–

36

October) accounted for approximately 70% of the annual total, with flux rates at one to

37

two orders of magnitude greater than those for rice paddies in this area. Ebullition

38

contributed to 44–56% of the overall CH4 emissions from the rivers and dominated the

39

emission pathways during the summer months. Our data highlight that rivers draining

40

agricultural watersheds may constitute a larger component of anthropogenic CH4

41

emissions than is currently documented in China.

42 43 44 45 46 2

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

Environmental Science & Technology

47

TOC (Table of Contents)

48 49 50

 INTRODUCTION

51

Methane (CH4) is a potent greenhouse gas (GHG) that exhibits a relative flux-sustained

52

global warming potential (SGWP) of 45 times that of carbon dioxide (CO2) on a 100-year

53

time horizon.1 The attribution of observed rising biogenic CH4 emissions to anthropogenic

54

and natural sources is difficult since the current network is insufficient to characterize

55

emission rates by region and source process.1 While natural wetlands and rice paddies

56

have been well documented as major sources of atmospheric CH4, freshwaters impacted

57

by agriculture as an anthropogenic source of CH4 are poorly represented in global CH4

58

budgets.2–6

59

While several studies have dedicated to quantify global flux strength of CH4 from

60

freshwaters, the current estimates of global CH4 emissions from freshwater bodies are still

61

poorly constrained, primarily due to the small amount of data and limited geographic

62

distribution of measurements. Global CH4 emissions from freshwater are estimated to be 3

ACS Paragon Plus Environment

Environmental Science & Technology

63

40-231 Tg yr-1, showing high uncertainties across studies.3,4,7 The current estimate of CH4

64

emissions from rivers even spans more than one order of magnitude (1.5-28 Tg yr-1).3,5,8,9

65

Global CH4 emissions from rivers were estimated to be 1.5 Tg yr-1 in an earlier study but

66

were only based on 13 measurements from nine studies, dramatically underestimated for

67

the tropics.3,9 By summarizing a larger dataset, Stanley et al.5 recently updated the

68

estimate rate to be 26.8-28.0 Tg yr-1. As the authors cautioned, however, the geographic

69

distribution of studies is still clustered, with heavy representation from Europe and North

70

America, while it is extremely limited in vast continental areas, such as in Asia, the

71

Middle East and Australia.3,5,8

72

Given that CH4 emissions from waters are closely related to organic material and

73

nutrients inputs,6,10,11 field measurements of CH4 fluxes from rivers draining agricultural

74

watersheds would therefore be highly needed due to high nutrient and sediment loading

75

rates through cropland drainage, leaching, runoff and soil erosion.7 High nutrient and

76

sediment loads can stimulate CH4 production by providing microbial communities with

77

organic substrates, and algal blooms from excessive nutrient loads can further add to this

78

with more enriched labile C. These factors may together result in higher CH4 fluxes from

79

agricultural waters relative to other types of freshwater or freshwaters in more natural

80

landscape areas.12–14 However, low rates of methanogenesis and the coldest environment

81

unfavorable for methanogenesis during non-cropping months may balance high CH4

82

emission rates during cropping duration, complicating efforts to predict the annual

83

magnitude of CH4 emissions from rivers draining agricultural watersheds. Unfortunately,

84

few measurements of CH4 fluxes from rivers have been taken on an annual basis, which 4

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Environmental Science & Technology

85

would limit our ability to assess the magnitude of CH4 emissions from waters impacted by

86

agriculture.6,13,15–17

87

Methane produced in sediments can be transported to the water surface via gas bubble

88

ebullition, molecular diffusion, or advection through plant stems.13,16,18,19 Diffusive fluxes

89

are determined from the water-air CH4 concentration gradient and piston velocity, which

90

have been most frequently employed in previous studies. Although eddy-covariance

91

method allows to capture ebullition of CH4 with fine resolution, ebullition fluxes have

92

been poorly represented for being episodic and not representatively captured by the usual

93

short-term chamber method.3,5,20 The previous efforts in tropical and subtropical rivers21-24

94

have highlighted the potential importance of rivers as a source of atmospheric CH4, and

95

recent studies have emphasized the importance of ebullition in shallow flowing

96

waters,8,13,25 while few flux measurements have included ebullitive fluxes in the recently

97

updated global database.5,26 Therefore, more studies on CH4 emissions from rivers are

98

highly expected to include the evaluation of both ebullition and diffusive components.

99

China is now the largest producer of rice and the second largest irrigator in the

100

world.27 Frequent irrigation and drainage episodes are currently practiced in rice paddies,

101

which have been threatening the rivers draining them in China with increased sediment

102

loading and nutrient enrichment.28,29 While midseason drainage and moist irrigation can

103

reduce seasonal CH4 emissions by as much as 50% in rice paddies,28 this reduction may

104

be offset by increased CH4 emissions from rivers draining rice paddy watersheds in China.

105

However, little is known about CH4 fluxes from rivers that receive high C and N loads

106

through cropland drainage, leaching, runoff and soil erosion.14 A model study estimated 5

ACS Paragon Plus Environment

Environmental Science & Technology

107

total GHG emissions from agricultural irrigation to be 36.72-54.16 Tg yr-1 in China,30

108

while this model estimate needs to be further examined by flux measurements.

109

Here, we present the results from a two-year measurement of CH4 fluxes from rivers

110

draining a rice paddy watershed in southeast China. The main objectives of this study

111

were first to gain insights into the magnitude of annual CH4 fluxes separated into diffusive

112

and ebullitive components from rivers draining rice paddy watersheds. We also aimed to

113

examine the dependence of seasonal CH4 fluxes on water and soil/sediment parameters.

114

Overall, we attempted to shed light on the role of rivers draining agricultural watersheds

115

in the national anthropogenic CH4 budget in China.

116 117

 MATERIALS AND METHODS

118

Study Site. Lixiahe River, a tributary of the Grand Canal, represents a typical river

119

impacted by lowland agriculture in southeast China. It is located in Xinghua, Jiangsu

120

Province, China (32°52’N, 119°50’E), where it has been dominated by paddy

121

rice-cropping systems over hundreds of years. This rice paddy watershed is characterized

122

by a subtropical monsoon climate, with an annual mean temperature of 17.8 °C and

123

precipitation of 1090 mm. High organic material and N fertilizer inputs (typically applied

124

at a seasonal rate of ~300 kg N ha-1) as well as different irrigation regimes (e.g.,

125

waterlogging, midseason drainage and moist irrigation episodes) are currently practiced in

126

rice paddies.29 The sampling area is 5 km downstream of the headwater stream of the river.

127

The main river of the water body is 8 m wide, with an average water depth of 2.5 m. The

128

average velocity of the water was approximately 0.5 m s-1 over the sampling days. Fall 6

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Environmental Science & Technology

129

turnover occurs during the month of October to November, when thermal stratification

130

begins to degrade and isothermal conditions are achieved. The details of water and

131

sediment physicochemical properties are summarized in Table S1. Over the two-year

132

experimental period, the chemical oxygen demand (COD) concentration in surface water

133

(0-20 cm depth) averaged 55.62 mg L-1, ranging from 15.9 to 132.6 mg L-1.

134

Experimental Design. Three typical measurement sites, including main-river (MP),

135

intersection (IP) and side-river points (SP), were located to simultaneously collect both

136

the gas and surface water samples. Gas flux measurements and water samples were taken

137

4–6 times a month, and sediment samples were collected twice a month. The samples

138

were taken mainly between 9:00 and 11:00 h to minimize potential discrepancies due to

139

diurnal variations.30 All the measurements at each site were taken with three replicates.

140

CH4 Flux Measurements. The fluxes of CH4 were measured using the floating

141

chamber method from September 15, 2014 to September 15, 2016. The chamber

142

deployment, sampling procedures and configuration of the GC for measuring CH4 fluxes

143

are extensively included in the Supporting Information (see Supporting Information, SI,

144

Materials and Methods-Extensive), as detailed in our previous studies.28,31-34

145

A water-air gas transfer equation was used to calculated the diffusive CH4 fluxes from

146

surface water to the atmosphere following previously documented methods, whereby a

147

gas-transfer coefficient of velocity across the water-air interface was adopted in addition

148

to the water dissolved CH4 concentrations in equilibrium with the atmosphere.3,8,20

149 150

=

(



) × 1000

(1)

Where F is the flux of CH4 (mmol m-2 d-1), k is the piston velocity of CH4 across the 7

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 31

151

water-air interface (m d-1), Cm is the concentration (mol m-3) of CH4 measured in the

152

surface water, and Ce is the concentration (mol m-3) of CH4 in the surface water at

153

equilibrium with the CH4 partial pressure in the floating chamber (typically calculated

154

from Henrys law).

155

As Sawakuchi et al.8 noted, however, the flux is partially driven by the change in

156

concentration in Eqn. (1), which shows that the flux will decrease with time as the

157

concentration increases inside the chambers. Therefore, this simple calculation will

158

underestimate the instantaneous flux rate. To reduce this error, we solve for k by using the

159

common gas law (PV = nRT) and make the equation continuous (dP/dt instead of Pt – P0).

160

Hence, the floating chamber-based total CH4 fluxes were determined by the Eqn. (2)

161

instead of the Eqn. (1)8,34:

162

×

=

× 1000

(2)

163

Where dP/dt is the slope of CH4 accumulation in the chamber (Pa d-1) over the

164

sampling period, V is the chamber volume (m3), A is the chamber enclosed surface area

165

(m2), R is the universal gas constant (8.314 Pa K-1 mol-1), and T is the temperature (K).

166

Diffusive and Ebullition Components. To determine which chambers captured

167

ebullition, we used the distributions and variance of the apparent k values as described in

168

Bastviken et al.20,36 and Sawakuchi et al.8. According to Henrys Law (C = KhP), the CH4

169

concentration numbers are also expressed as corresponding gas pressure. Thus, Eqn. (1)

170

could be rewritten as:

171 172

=

(



) × 1000

(3)

where Kh is the Henry’s Law constant for CH4 (mol m-3 Pa-1), Pm is the partial 8

ACS Paragon Plus Environment

Page 9 of 31

Environmental Science & Technology

173

pressure of CH4 in the chamber at equilibrium with Cm (Pa), and P0 is the partial pressure

174

of CH4 in the chamber at time 0 (approximately the same as atmosphere). Combining

175

Eqns (2) and (3), thus, k is calculated by the equation:

176

(

×

=

)

(4)

177

Further, the calculated apparent k values for each chamber were normalized to

178

Schmidt numbers of 600 (k600), allowing k values to be compared for any gas and

179

temperature.

180

=

(

)

.

(5)

181

where Sc is the Schmidt number of CH4 at a given temperature.35

182

First, we normalized the calculated apparent k values for each chamber to Schmidt

183

numbers of 600 (k600), allowing k values to be compared for any gas and temperature.

184

Assuming that the CH4 emission consists of ebullition and diffusive fluxes, ebullition

185

makes the calculated apparent k600 values substantially greater than those receiving CH4

186

only by diffusive flux, allowing the separation of the two flux components. Given that the

187

diffusive flux is assumed to have a constant rate at a given sampling site and time, the

188

chambers receiving only diffusive flux would have lower and similar flux rates and could

189

be distinguished from the chambers receiving ebullition. However, once all chambers at

190

one sampling site showed a large discrepancy in flux rates, suggesting that all chambers

191

could have received ebullition. In this case, the lowest value similar with the diffusion

192

measured in nearby sites was assumed as diffusion. In each measurement, the minimum

193

k600 of CH4 (minimum k600-CH4) found in the same set of chambers was attributed solely to

194

diffusive flux and we calculated the ratio of apparent k600 of CH4 (k600-CH4) and minimum 9

ACS Paragon Plus Environment

Environmental Science & Technology

195

k600 of CH4 (k600-CH4/minimun k600-CH4) for each chamber.

196

In this study, the frequency distribution of this ratio (k600-CH4/minimun k600-CH4) for all

197

chambers clearly showed two distinct groups of ratios, one between 1.0 and 2.5 and

198

another > 2.5 (Figure 1). Thereby, a ratio of 2.5 was accepted as threshold for significant

199

contribution of ebullition to total CH4 flux. For the chambers that received ebullition flux,

200

we first calculated the diffusive flux by Eqn. (1) using the averaged k600 from the other

201

chambers that received diffusive flux only, and the remaining CH4 flux into the chambers

202

was attributed to ebullition. Meanwhile, we also calculated k600 values of CO2 (k600-CO2)

203

based on CO2 data in comparison with CH4 for each chamber. If the CH4 emissions were

204

purely diffusive, then the values of k600-CH4 and k600-CO2 should tend to be identical. The

205

ratio of k600-CH4 to k600-CO2 (k600-CH4/k600-CO2) was used as an indicator of ebullition when

206

values of k600-CH4 far exceed those of k600-CO2 for the chambers capturing CH4 rich bubbles.

Measurement numbers

80

60

40

20

0

Calculated k 600-CH4/Minimum k 600-CH4

207 208

Figure 1. Frequency distribution of the calculated k600-CH4/minimum k600-CH4 ratio for CH4

209

measurements using the floating chamber method, used for distinction between fluxes

210

consisting of both ebullition and diffusion and those of only diffusion. 10

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

Environmental Science & Technology

211 212

While gas sampling, surface water samples (0–30 cm depth) were collected to

213

measure the dissolved CH4 concentrations. Water samples were collected using PTE vials

214

(100 mL) and the PTE vials were injected by Helium (He, 99.999%) gas to stop biological

215

alteration of water samples. After water sampling, the vials were immediately stored at

216

4 °C during transportation prior to analysis within two hours. The wind speed at

217

approximately 1 m above the water surface was automatically measured with a previously

218

installed weather information recorder during gas and water sampling. Gas and water

219

samples were simultaneously taken with wind speed ranging from 0.41 to 6.4 m s-1. The

220

headspace equilibrium technique was used to measure dissolved CH4 concentrations.36,37

221

Other data measurements are included in the Supporting Information.

222

Statistical Analyses. Differences in CH4 emissions determined by the floating

223

chamber method among the three sampling sites over the two-year period were examined

224

by a two-way analysis of variance (ANOVA). Linear regressions were used to identify the

225

key factors influencing CH4 fluxes from agricultural irrigation waters. All statistical

226

analyses were performed using SPSS version 19.0 (SPSS Inc., USA), and statistical

227

significance was determined at the 0.05 probability level.

11

ACS Paragon Plus Environment

Environmental Science & Technology

228 12

ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31

Environmental Science & Technology

229

Figure 2. Seasonal dynamics of dissolved CH4 concentrations (a), floating chamber-based

230

CH4 fluxes (b) and the water parameters (c-e) in the river draining a rice paddy watershed.

231

Bars represent the mean ± SE (n = 3). The sampling sites of MP, IP and SP refer to

232

main-river, intersection and side-river points, respectively.

233 234

 RESULTS

235

Dissolved CH4 Concentrations, CH4 Emissions and Water Parameters. Seasonal

236

dissolved CH4 concentrations and CH4 emissions showed similar variation patterns

237

(Figure 2a and 2b). These patterns did not differ among sampling sites or between the two

238

experimental years (Figure 2a and 2b). In general, both dissolved CH4 concentrations and

239

CH4 fluxes ascended steadily during the spring (March-May) until they attained their

240

highest values in August over the summer months (June-August). Thereafter, they

241

decreased dramatically and remained relatively lower during the autumn until they

242

reached their lowest point in the winter (December-February). Over the two-year period,

243

the mean and median of dissolved CH4 concentrations of surface water were 0.42 and

244

0.24 µmol L-1, respectively, ranging from a seasonal mean of 0.22 µmol L-1 in the winter

245

to 0.81 µmol L-1 in the summer. Among sampling sites, mean annual dissolved CH4

246

concentrations were higher at the MP site (0.46 µmol L-1) relative to the IP and SP sites

247

(0.39 µmol L-1).

248

The seasonal pattern of CH4 fluxes was generally similar to those of water

249

temperature and water dissolved organic carbon (DOC), while it was contrary to that of

250

water dissolved oxygen (DO). Over the two experimental years, single CH4 measurements 13

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 31

251

greatly differed with time and site (ANOVA, p < 0.001), ranging from 0.07 to 297.86

252

mmol m-2 d-1 at the MP site, 0.05 to 155.19 mmol m-2 d-1 at the IP site, and 0.04 to 171.60

253

mmol m-2 d-1 at the SP site. The standard error of a single measurement ranged from 0.01

254

to 83.10 mmol m-2 d-1 at the MP site, 0.01 to 35.74 mmol m-2 d-1 at the IP site, and 0.01 to

255

36.59 mmol m-2 d-1 at the SP site (Figure 2b). Similar to CH4 fluxes, the water

256

temperature was lowest in December to January and highest in July to August. When

257

averaged across sites, the mean annual water temperature was 18.98 °C and 19.85 °C in

258

the first and second years, respectively (Figure 2c). Over the two-year period, the mean

259

(or median) water DO content of 258 measurements was 13.12 mg L-1 over the three

260

sampling sites, with a variation from 0.34 to 30.67 mg L-1 (Figure 2d). Across sampling

261

sites, water DOC concentrations averaged 78.1 mg L-1 and 92.6 mg L−1, ranging from

262

25.2 to 261.3 mg L-1 and from 17.2 to 283.6 mg L-1 in the first and second years,

263

respectively (Figure 2e).

Annual CH4 emissions (mol m-2)

30 Ebullition 25

Site (S): p < 0.01 Year (Y): p < 0.01

Diffusion

S×Y: non-significant

20

15

10

5

0 SP

IP

MP

First-year

SP

IP

MP

Second-year

264 265

Figure 3. Annual total of CH4 emissions and diffusive and ebullitive components at the 14

ACS Paragon Plus Environment

Page 15 of 31

Environmental Science & Technology

266

three sampling sites in agricultural irrigation waters over the two experimental years. Bars

267

represent the mean ± SE (n = 3).

268 269

Annual total of CH4 fluxes greatly differed with the sampling site and experimental

270

year, while they were not significantly affected by their interaction (Figure 3). At an

271

identical sampling site, average annual CH4 fluxes were significantly greater in the second

272

experimental year (September 2015–September 2016) than in the first year (Figure 3).

273

Among the three sampling sites, CH4 fluxes were the greatest at the MP site, with an

274

average of 29.33 and 64.80 mmol m-2 d-1 in the first and second experimental years,

275

respectively. Annual CH4 fluxes at the SP site averaged 13.13 mmol m-2 d-1, amounting to

276

4.79 mol m-2 in the first year, which was 44% lower than the average fluxes in the second

277

year (Figure 3). Over the two experimental years, annual CH4 emissions at the IP site

278

totaled 6.52–10.33 mol m-2, with an annual average flux of 17.85–28.31 mmol m-2 d-1,

279

which was 16–27% greater than that at the SP site (Figure 3). In contrast, annual CH4

280

fluxes averaged 20.10 and 38.95 mmol m-2 d-1, amounting to 7.34 and 14.22 mol m-2 in

281

the first and second years, respectively. Across the two experimental years, annual CH4

282

emissions totaled 10.78 mol m-2, with an average flux of 29.52 mmol m-2 d-1. Substantial

283

CH4 emissions occurred during the rice-growing season (June–October), accounting for

284

68-73% of the annual total.

15

ACS Paragon Plus Environment

Environmental Science & Technology

Page 16 of 31

285

Table 1. Seasonal variations in CH4 fluxes (Mean ± SE, mmol m-2 d-1), diffusive and ebullitive components, and k600-CH4/k600-CO2 for the river

286

draining rice paddy watersheds. First year

Second year

Autumn

Winter

Total fluxes

5.40 ± 1.53

0.62 ± 0.02

1.05 ± 0.24

68.66 ± 10.28

12.62 ± 4.08

3.71 ± 0.69

38.94 ± 6.89

96.05 ± 13.26

Diffusive fluxes

2.88 ± 0.35

0.62 ± 0.02

0.80 ± 0.08

8.72 ± 0.71

8.64 ± 1.76

2.54 ± 0.56

8.03 ± 1.01

9.72 ± 0.74

Ebullition

2.52 ± 1.52

0

0.26 ± 0.12

59.94 ± 9.87

3.98 ± 0.96

1.17 ± 0.11

30.92 ± 6.75

86.33 ± 12.87

47

0

24

87

32

31

79

89

0.94 ± 0.09

4.59 ± 1.38

33.52 ± 5.16

13.30 ± 3.88

5.27 ± 1.48

28.50 ± 4.16

58.80 ± 11.70

Ebullition (%) k600-CH4/k600-CO2

18.95 ± 6.32

Spring

Summer

Autumn

16

ACS Paragon Plus Environment

Winter

Spring

Summer

Page 17 of 31

Environmental Science & Technology

120

80

k600-CH4/k600-CO2

(a)-Spring

SP

IP

(b)-Summer

MP

60

90

40

60

20

30

0

0 Mar

Apr

May

Jun

120

Jul

Aug

Jan

Feb

20 (c)-Autumn

(d)-Winter

100

16

80 12 60 8 40 4

20

0

0 Sept

Oct

Nov

Dec Month

287 288

Figure 4. Ratios of k600-CH4 to k600-CO2 in four seasons from September 2014 to September

289

2016 in agricultural irrigation waters (error bars SE, n =3).

290 291

Diffusive and Ebullitive Fluxes. Similar to the total CH4 fluxes, both diffusive and

292

ebullitive fluxes showed parallel seasonal variations, with the lowest rates in the winter

293

and the highest in the summer (Table 1). The seasonal average of diffusive fluxes ranged

294

from 0.62 to 8.72 mmol m-2 d-1 in the first year, which was lower than that with a range of

295

2.54–9.72 mmol m-2 d-1 in the second year (Table 1). Diffusion was the main component

296

of CH4 emissions when total CH4 fluxes were under 12.50 mmol m-2 d-1, particularly in

297

the winter and late autumn. With the exception that ebullition was not pronounced in the

298

winter of the first year, seasonal ebullitive fluxes showed large variation over the seasons, 17

ACS Paragon Plus Environment

Environmental Science & Technology

299

with the lowest rate of 0.26 mmol m-2 d-1 in the spring of the first year and the highest rate

300

of 86.33 mmol m-2 d-1 in the summer of the second year (Table 1). On average, diffusion

301

contributed approximately 56% and 42% to the total CH4 emissions, while almost 44%

302

and 58% of total CH4 emissions could be attributed to ebullition in the first and second

303

years, respectively (Figure 3). Among the three sampling sites, ebullitive fluxes were the

304

highest at the MP site, while they were lowest at the SP site (Figure 3).

305

Piston Velocities for CO2 and CH4. We estimated the median and mean of gas piston

306

velocities (k600) to be 0.27 and 0.28 m d-1 for CO2 and 2.32 and 6.74 m d-1 for CH4,

307

respectively. The ratio of k600-CH4 to k600-CO2, referring to the difference between the CH4

308

and CO2 k600 values, was much higher in the summer and early autumn, while it was

309

relatively lower in the winter (Figure 4). Over the two experimental years, the mean ratio

310

of k600-CH4 to k600-CO2 was estimated to be 2.9 in the winter, 16.5–21.8 in the spring and

311

autumn, and 45.1 in the summer. Overall, greater k600 values for CH4 relative to CO2

312

suggested that floating chambers intercepted CH4-rich bubbles rising to the water surface,

313

particularly in the summer and early autumn. On the other hand, the difference between

314

the k600-CH4 and k600-CO2 values was most extreme at the MP site, except for a relatively

315

higher ratio at the IP site in September (Figure 4).

316

Dependence of CH4 Fluxes on Water/Sediment Parameters. By pooling

317

measurements over the two years, the CH4 fluxes were positively related to water

318

temperature, DOC and NO3--N, while they were negatively related to water DO

319

concentrations (Figure S1). Relative to water parameters, CH4 fluxes were positively

320

related to sediment pH, electrical conductivity, organic matter and DOC (Figure S2).

321 322 323

 DISCUSSION CH4 Emissions in Comparison with Previous Estimates. The CH4 fluxes in his 18

ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

Environmental Science & Technology

324

study were generally within the range previously reported in tropical and subtropical lakes

325

and rivers, but had relatively higher release levels.12,38-41 Several explanations may be

326

given for higher CH4 emissions in this study. First, we believe that this may be associated

327

with our high frequency of sampling measurements, which was typically once a week,

328

allowing us to capture more flux peaks (Figure 2). Second, the use of floating chambers at

329

different sampling sites enabled us to capture a larger number of events with ebullitive

330

fluxes.8 Annual total CH4 emissions, including both diffusive and ebullitive fluxes from

331

the rivers, averaged 29.52 mmol m-2 d-1 across the two years, amounting to 10.78 mol m-2

332

yr-1 in this study. This rate is generally greater than in previous studies in rivers impacted

333

by agriculture in this subtropical watershed using only dissolved CH4 concentration data

334

to calculate passive diffusive fluxes of CH4.22,23 In this watershed, previous studies

335

estimated annual diffusive CH4 emissions to be 0.53 mol m-2 in hypereutrophic Donghu

336

Lake and 0.04–6.50 mol m-2 in Taihu Lake,22,23 which is almost comparable to our

337

diffusive CH4 emission rates (2.88–7.25 mol m-2). Third, high nutrient loading,

338

sedimentation and algae blooms in the summer months can stimulate methanogenesis in

339

rivers draining rice paddies, resulting in CH4 emissions that are higher than previously

340

reported in rivers and lakes not impacted by agriculture8 but that are fairly comparable to

341

the estimates of studies on subtropical Mexican eutrophic rivers and streams draining

342

human-dominated landscapes in Wisconsin, USA.12,14 In addition, as a consequence of

343

cropland irrigation and drainage episodes, outgassing to the atmosphere of dissolved CH4

344

delivered by drainage flows may also contribute to CH4 emissions from rivers draining

345

rice paddies, as similarly documented in reservoirs.12

346

Spatiotemporal Variation in CH4 Fluxes. Our results showed that CH4 fluxes were

347

the greatest at the MP site among sampling sites. The greatest CH4 fluxes at the MP site

348

could be primarily ascribed to the nutrient and sediment loads through drainage water, 19

ACS Paragon Plus Environment

Environmental Science & Technology

349

leaching, runoff and soil erosion directly into the main-river location, rather than the

350

intersection and side-river locations. Over the two-year period, the mean and median of

351

sediment DOC contents at the MP site were 27–48% and 24–41% greater than those at the

352

IP and SP sites, respectively. In the river water, mean NH4+-N and NO3 --N contents at the

353

MP site were also 18–27% and 10–14% greater than those at the IP and SP sites,

354

respectively. These higher sediment organic C and nutrient loads would have contributed

355

to greater CH4 emissions at the MP site relative to the IP and SP sites,41,42 as supported by

356

the evidence that CH4 emissions were highly related to sediment C and water nutrient

357

loads in this study (Supporting Information, Figures S1 and S2). In addition, greater CH4

358

fluxes could also be associated with higher water-dissolved CH4 concentrations and a

359

larger ebullition contribution at the MP site relative to at the other two sites. Among

360

sampling sites, mean annual dissolved CH4 concentrations and the contribution of

361

ebullitive CH4 fluxes to the overall total CH4 fluxes at the MP site were 14–16% and 20–

362

36% greater than those at IP and SP sites, respectively (Table 1, Figure 3). Therefore, both

363

diffusive and ebullitive fluxes contributed to greater CH4 fluxes at the MP site than at the

364

other two sites.

365

Seasonal dynamics of CH4 fluxes from the rivers were typically characterized by a

366

pattern in which CH4 fluxes were the lowest in the winter and the highest in the summer,

367

and small peaks frequently happened during the late spring and early autumn (Figure 2).

368

This seasonal pattern was in line with those reported in other waters impacted by

369

agriculture.10,18 The seasonal pattern of CH4 fluxes was generally in agreement with those

370

of water temperature and DOC but in contrast to that of water DO. On average, mean

371

monthly fluxes of CH4 showed two orders of magnitude differences between the winter

372

and summer seasons (Table 1). Low CH4 fluxes in the winter were primarily ascribed to

373

low C and nutrient loads, low water temperature and high O2 availability in bottom water 20

ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

Environmental Science & Technology

374

and sediments, which reduces the substrate availability and metabolic activity of

375

methanogens.43 The higher dissolved CH4 during the 2016 autumn season was not

376

accompanied by a higher CH4 flux, which could be associated with high surface water DO

377

(Figure 2) and strong thermal stratification that prevented the advective transport of CH4

378

from the hypolimnion to the epilimnion11 and then reduced ebullitive CH4 emissions

379

(Table 1). In contrast, surface CH4 emissions from the waters reached a maximum in the

380

summer months. Intensive nutrient and sediments loads, high water temperature and low

381

O2 availability in bottom waters and sediments would facilitate sediment CH4 production

382

and diffusive CH4 fluxes in the summer. Greater CH4 emissions could also be ascribed to

383

higher ebullitive CH4 fluxes in the summer (Table 1), as supported by higher ratios of

384

k600-CH4 to k600-CO2 in the summer than in other seasons (Figure 4). In particular, frequent

385

water withdrawal and recharge courses due to frequent irrigation and drainage in rice

386

paddies would create the environments where algal blooms from excessive nutrient loads

387

could further enrich water and sediment labile C and deplete water O2. Therefore, our

388

results suggested that the seasonal variation in CH4 emissions could be larger in rivers

389

draining agricultural watersheds than in watersheds not impacted by agriculture.

390

Moreover, CH4 emissions also showed large interannual variations in this study. On

391

average, annual CH4 emissions were significantly greater in the second year than in the

392

first year (Figure 2), which is attributed to a combination of water and sediment variables,

393

including water temperature (mean annual: 18.98 ºC in the first year vs. 19.85 ºC in the

394

second year), DOC (78.1 vs. 92.6 mg L-1), DO (15.76 vs. 10.58 mg L-1), and sediment

395

SOC (77.18 vs. 90.89 mg kg-1). All these factors together would have consistently

396

contributed to greater CH4 emissions in the second year than in the first year. In addition

397

to higher mean annual dissolved CH4 concentrations in water (0.34 vs. 0.49 µmol L-1),

398

mean seasonal ebullitive CH4 fluxes were much greater in the second year than in the first 21

ACS Paragon Plus Environment

Environmental Science & Technology

399

year (Table 1). Therefore, both diffusive and ebullitive fluxes were responsible for

400

interannual variations of CH4 emissions between the two experimental years.

401

Ebullitive CH4 Fluxes. The current estimates of freshwater CH4 emissions are still

402

constrained with a poor assessment of ebullitive fluxes.3,5 In the present study, the floating

403

chambers captured both diffusive and ebullitive CH4 fluxes from the rivers draining rice

404

paddy watersheds. The ebullitive CH4 fluxes were responsible for more than 50% of the

405

total emissions over the two years. This value fell well within the range of 10–80%, as

406

reported on the contribution of ebullition to the total CH4 emissions from streams and

407

rivers.8,25,44,45 For example, Sawakuchi et al.8 found that ebullition contributed to more

408

than 50% of the total CH4 emissions for some rivers in the Amazon Basin. In particular,

409

ebullitive fluxes showed a large seasonal variation, and summer ebullitive fluxes

410

contributed to almost 90% of the total CH4 emissions. High ebullitive CH4 fluxes are

411

primarily ascribed to large organic material loads fueling CH4-containing bubble

412

formation in sediments and high water temperature facilitating bubble transportation to

413

the atmosphere in shallow freshwaters.46,47 This may be intensified in watersheds

414

dominated by rice paddies and wetlands subject to high soil erosion rates and sediment

415

loading.

416

The k600 values for CH4 relative to CO2 also provided insight into the relative

417

importance of diffusive and ebullitive pathways. The k600-CO2 values suggested that

418

diffusion was the dominant CO2 emission pathway, and the high value of k600-CH4 may

419

imply that the CH4 emission pathway attributed more to ebullition rather than diffusion.48

420

The mean and median of k600-CO2 values were comparable to values reported for diffusive

421

emissions, suggesting that the floating chambers did not artificially disturb the gas

422

exchange between the air and water interface.13,18,49 The mean k600-CH4 values were in good

423

agreement with those reported in a large river and a reservoir by Beaulieu et al.13,49, but 22

ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

Environmental Science & Technology

424

greater than the values reported in boreal lakes.18 Greater piston velocities calculated from

425

CH4 relative to CO2 were found for the majority of measurements in this study, suggesting

426

the presence of surfacing CH4-rich bubbles in addition to gas diffusive fluxes.50,51 In

427

particular, the k600-CH4 values and ratios of k600-CH4 to k600-CO2 were greatest and extremely

428

high in the summer months and at the MP site (Figure 4), suggesting that wind speed and

429

water velocity at the main-river location relative to the intersection and side-river location

430

may favor greater sediment formation of CH4-containing bubbles and their transportation,

431

especially in summer months.

432

In the present study, the rivers draining rice paddy watersheds showed a consistent

433

source of atmospheric CH4 throughout the years. On average, annual CH4 emissions from

434

the rivers draining rice paddy watersheds in this study were one to 2 orders of magnitude

435

greater than previously reported for rice paddies in this watershed.28,31,33,39 Given field

436

drainage has been fairly approved to be an effective practice to reduce CH4 emissions

437

from rice paddies, this reduction could be offset by high CH4 emissions as a consequence

438

of nutrient enrichment and sediment loading in rivers draining rice paddy watersheds.

439

While rice paddies have been well recognized as one of major sources of CH4 emissions,

440

this study highlights the importance of including rivers draining agricultural watersheds in

441

national anthropogenic CH4 budget in China given that China is the largest paddy rice

442

production and the second agricultural irrigation countries in the world.

443 444

 AUTHOR INFORMATION

445

Corresponding author

446

*Phone: +86 25 8439 6286; Fax: +86 25 8439 5210; Email: [email protected]

447

Author Contributions

448

1

S.W. and S.L. contributed equally to this work 23

ACS Paragon Plus Environment

Environmental Science & Technology

449

Notes

450

The authors declare no competing financial interest.

451 452

 ACKNOWLEDGMENTS

453

This work was supported by the Natural Science Foundation of China (41301244,

454

41771268), the National Key Research and Development Program of China

455

(2016YFD0201200), the Fundamental Research Funds for the Central Universities

456

(KYZ201621, KJQN201516) and PADA.

457 458

 SUPPORTING INFORMATION

459

Supplementary Texts

460

The details on CH4 flux measurements and the other data measurements are integral parts

461

of the ‘Materials and Methods’ and the extensive “Discussion” on the dependence of CH4

462

emissions on water and sediment parameters and mitigation strategies.

463

Supplementary Tables

464

Table S1. Physicochemical properties of water and sediment (0–20 cm depth) in the river.

465

Supplementary Figures

466

Figure S1. Dependence of CH4 fluxes on water parameters in agricultural irrigation water.

467

(a) water temperature; (b) water dissolved oxygen (DO); (c) water dissolved organic

468

carbon (DOC); (d) water nitrate concentrations (NO3--N).

469

Figure S2. Dependence of CH4 fluxes on sediment parameters in rivers that drain rice

470

paddy watersheds. (a) sediment pH; (b) sediment electrical conductivity; (c) sediment

471

organic matters; (d) sediment dissolved organic carbon (DOC).

472 473



REFERENCES 24

ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

Environmental Science & Technology

474

(1) Intergovernmental Panel on Climate Change (IPCC) Climate Change 2013: The

475

Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report

476

of the Intergovernmental Panel on Climate Change; Cambridge University Press,

477

Cambridge, 2013.

478

(2) Cole, J. J.; Prairie, Y. T.; Caraco, N. F.; McDowell, W. H.; Tranvik, L. J.; Striegl, R. G.;

479

Duarte, C. M.; Kortelainen, P.; Downing, J. A.; Middelburg, J. J. Plumbing the global

480

carbon cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems

481

2007, 10 (1), 171–184.

482

(3) Bastviken, D.; Tranvik, L. J.; Downing, J. A.; Crill, P. M.; Enrich-Prast, A. Freshwater

483

methane emissions offset the continental carbon sink. Science 2011, 331 (6013), 50.

484

(4) Kirschke, S.; Bousquet, P.; Ciais, P.; Saunois, M.; Canadell, J. G.; Dlugokencky, E. J.;

485

Bergamaschi, P.; Bergmann, D.; Blake, D. R.; Bruhwiler, L. Three decades of global

486

methane sources and sinks. Nat. Geosci. 2013, 6 (10), 813–823.

487

(5) Stanley, E. H.; Casson, N. J.; Christel, S. T.; Crawford, J. T.; Loken, L. C.; Oliver, S. K.

488

The ecology of methane in streams and rivers: Patterns, controls, and global significance.

489

Ecol. Monogr. 2016, 86 (2), 146–171.

490

(6) Borges, A. V.; Darchambeau, F.; Lambert, T.; Bouillon, S.; Morana, C.; Brouyère, S.;

491

Hakoun, V.; Jurado, A.; Tseng, H.-C.; Descy, J.-P.; Roland, F. A. E. Effects of agricultural

492

land use on fluvial carbon dioxide, methane and nitrous oxide concentrations in a large

493

European river, the Meuse (Belgium). Sci. Total Environ. 2018, 610-611: 342–355.

494

(7) O’Connor, F. M.; Boucher, O.; Gedney, N.; Jones, C. D.; Folberth, G. A.; Coppell, R.;

495

Friedlingstein, P.; Collins, W. J.; Chappellaz, J.; Ridley, J. Possible role of wetlands,

496

permafrost, and methane hydrates in the methane cycle under future climate change: A

497

review. Rev. Geophys. 2010, 48 (4), 1–33.

498

(8) Sawakuchi, H. O.; Bastviken, D.; Sawakuchi, A. O.; Krusche, A. V.; Ballester, M. V. 25

ACS Paragon Plus Environment

Environmental Science & Technology

499

R.; Richey, J. E. Methane emissions from amazonian rivers and their contribution to the

500

global methane budget. Global Chang. Biol. 2014, 20 (9), 2829–2840.

501

(9) Borges, A. V.; Darchambeau, F.; Teodoru, C. R.; Marwick, T. R.; Tamooh, F.; Geeraert,

502

N.; Omengo, F. O.; Guérin, F.; Lambert, T.; Morana, C. Globally significant

503

greenhouse-gas emissions from African inland waters. Nat. Geosci. 2015, 8 (8), 637–642.

504

(10) Schrier-Uijl, A. P.; Veraart, A. J.; Leffelaar, P. A.; Berendse, F.; Veenendaal, E. M.

505

Release of CO2 and CH4 from lakes and drainage ditches in temperate wetlands.

506

Biogeochemistry 2011, 102 (1–3), 265–279.

507

(11) Beaulieu, J. J.; Smolenski, R. L.; Nietch, C. T.; Townsend-Small, A.; Elovitz, M. S.

508

High methane emissions from a midlatitude reservoir draining an agricultural watershed.

509

Environ. Sci. Technol. 2014, 48 (19), 11100–11108.

510

(12) Gonzalez-Valencia, R.; Sepulveda-Jauregui, A.; Martinez-Cruz, K.; Hoyos-Santillan,

511

J.; Dendooven, L.; Thalasso, F. Methane emissions from Mexican freshwater bodies:

512

Correlations with water pollution. Hydrobiologia 2014, 721 (1), 9–22.

513

(13) Tranvik, L. J.; Downing, J. A.; Cotner, J. B.; Loiselle, S. A.; Striegl, R. G.; Ballatore,

514

T. J.; Dillon, P.; Finlay, K.; Fortino, K.; Knoll, L. B. Lakes and reservoirs as regulators of

515

carbon cycling and climate. limnol. Oceanogr. 2009, 54 (6 PART 2), 2298–2314.

516

(14) Crawford, J. T.; Stanley, E. H. Controls on methane concentrations and fluxes in

517

streams draining human-dominated landscapes. Ecol. Appl. 2016, 26 (5), 1581–1591.

518

(15) Whalen, S. C.; Reeburgh, W. S. A methane flux time series for tundra environments.

519

Global Biogeochem. Cycles 1988, 2 (4), 399–409.

520

(16) Laini, A.; Bartoli, M.; Castaldi, S.; Viaroli, P.; Capri, E.; Trevisan, M. Greenhouse

521

gases (CO2, CH4 and N2O) in lowland springs within an agricultural impacted watershed

522

(Po River Plain, Northern Italy). Chem. Ecol. 2011, 27 (2), 177–187.

523

(17) Rantakari, M.; Heiskanen, J.; Mammarella, I.; Tulonen, T.; Linnaluoma, J.; Kankaala, 26

ACS Paragon Plus Environment

Page 26 of 31

Page 27 of 31

Environmental Science & Technology

524

P.; Ojala, A. Different apparent gas exchange coefficients for CO2 and CH4: comparing a

525

brown-water and a clear-water lake in the boreal zone during the whole growing season.

526

Environ. Sci. Technol. 2015, 49 (19), 11388–11394.

527

(18) Bastviken, D.; Cole, J. J.; Pace, M. L.; Van de-Bogert, M. C. Fates of methane from

528

different lake habitats: Connecting whole-lake budgets and CH4 emissions. J. Geophys.

529

Res. Biogeosciences 2008, 113 (2), 1–13.

530

(19) Wilcock, R. J.; Sorrell, B. K. Emissions of greenhouse gases CH4 and N2O from

531

low-gradient streams in agriculturally developed catchments. Water Air Soil Pollut. 2008,

532

188 (1–4), 155–170.

533

(20) Bastviken, D.; Cole, J.; Pace, M.; Tranvik, L. Methane emissions from lakes:

534

Dependence of lake characteristics, two regional assessments, and a global estimate.

535

Global Biogeochem. Cycles 2004, 18 (4), 1–12.

536

(21) Richey, J. E.; Devol, A. H.; Wofsy, S. C.; Victoria, R.; Riberio, M. N. G. Biogenic

537

gases and the oxidation and reduction of carbon in amazon river and floodplain waters.

538

Limnol. Oceanogr. 1988, 33 (4), 551–561.

539

(22) Xing, Y.; Xie, P.; Yang, H.; Ni, L.; Wang, Y.; Rong, K. Methane and carbon dioxide

540

fluxes from a shallow hypereutrophic subtropical lake in China. Atmos. Environ. 2005, 39

541

(30), 5532–5540.

542

(23) Wang, H.; Lu, J.; Wang, W.; Yang, L.; Yin, C. Methane fluxes from the littoral zone

543

of hypereutrophic Taihu lake, China. J. Geophys. Res. Atmos. 2006, 111 (17), 1–8.

544

(24) Koné, Y. J. M.; Abril, G.; Delille, B.; Borges, A. V. Seasonal variability of methane in

545

the rivers and lagoons of ivory coast (West Africa). Biogeochemistry 2010, 100 (1), 21–

546

37.

547

(25) Baulch, H. M.; Dillon, P. J.; Maranger, R.; Schiff, S. L. Diffusive and ebullitive

548

transport of methane and nitrous oxide from streams: Are bubble-mediated fluxes 27

ACS Paragon Plus Environment

Environmental Science & Technology

549

important? J. Geophys. Res. Biogeosciences 2011, 116 (G4).

550

(26) Stanley, E.; Loken, L.; Crawford, J.; Casson, N.; Oliver, S.; Gries, C.; Christel, S. A

551

global database of methane concentrations and atmospheric fluxes for streams and rivers.

552

Long Term Ecol. Res. Netw. 2015, https://portal.lternet.edu/nis/mapbrowse.

553

(27) Wu, K.; Lu, B.; Yuan, Z. The recent developments and the contribution of farmland

554

irrigation to national grain safeness in China. J. Irrig. Drain. 2006, 25 (4), 7–10 (in

555

Chinese).

556

(28) Zou, J.; Huang, Y.; Jiang, J.; Zheng, X.; Sass, R. L. A 3-year field measurement of

557

methane and nitrous oxide emissions from rice paddies in China: Effects of water regime,

558

crop residue, and fertilizer application. Global Biogeochem. Cycles 2005, 19 (2), 1–9.

559

(29) Zou, J.; Huang, Y.; Qin, Y.; Liu, S.; Shen, Q.; Pan, G.; Lu, Y.; Liu, Q. Changes in

560

fertilizer-induced direct N2O emissions from paddy fields during rice-growing season in

561

China between 1950s and 1990s. Global Change Biol. 2009, 15 (1), 229–242.

562

(30) Zou, X.; Li, Y.; Li, K.; Cremades, R.; Gao, Q.; Wan, Y.; Qin, X. Greenhouse gas

563

emissions from agricultural irrigation in China. Mitig. Adapt. Strateg. Global Change

564

2015, 20 (2), 295–315.

565

(31) Wu, S.; Hu, Z.; Hu, T.; Chen, J.; Yu, K.; Zou, J.; Liu, S. Annual methane and nitrous

566

oxide emissions from rice paddies and inland fish aquaculture wetlands in southeast China.

567

Atmos. Environ. 2018, 175, 135–144.

568

(32) Hu, Z.; Wu, S.; Ji, C.; Zou, J.; Zhou, Q.; Liu, S. A comparison of methane emissions

569

following rice paddies conversion to crab-fish farming wetlands in southeast China.

570

Environ. Sci. Pollut. Res. 2016, 23 (2), 1505–1515.

571

(33) Liu, S.; Hu, Z.; Wu, S.; Li, S.; Li, Z.; Zou, J. Methane and nitrous oxide emissions

572

reduced following conversion of rice paddies to inland crab-fish aquaculture in southeast

573

China. Environ. Sci. Technol. 2016, 50, 633–642. 28

ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

Environmental Science & Technology

574

(34) Gao, J.; Zheng, X.; Wang, R.; Zou, J. Preliminary comparison of the static floating

575

chamber and the diffusion model methods for measuring water−atmosphere exchanges of

576

methane and nitrous oxide from inland water bodies. Clim. Environ. Res. 2014, 19 (3):

577

290–302 (in Chinese).

578

(35) Wanninkhof, R. H. Relationship between wind speed and gas exchange over the

579

ocean. J. Geophys. Res. Ocean 1992, 97, 7373–7382.

580

(36) Bastviken, D.; Santoro, A. L.; Marotta, H.; Pinho, L. Q.; Calheiros, D. F.; Crill, P.;

581

Enrich-Prast, A. Methane emissions from pantanal, south America during the low water.

582

Environ. Sci. Technol. 2010, 44 (14), 5450–5455.

583

(37) Amouroux, D.; Roberts, G.; Rapsomanikis, S.; Andreae, M. O. Biogenic gas (CH4 ,

584

N2O , DMS ) emission to the atmosphere from near-shore and shelf waters of the

585

north-western Black Sea . Estuar. Coast. Shelf Sci. 2002, 54, 575–587.

586

(38) Huttunen, J. T.; Väisänen, T. S.; Hellsten, S. K.; Heikkinen, M.; Nykänen, H.;

587

Jungner, H.; Niskanen, A.; Virtanen, M. O.; Lindqvist, O. V.; Nenonen, O. S. Fluxes of

588

CH4 , CO2 , and N2O in hydroelectric reservoirs lokka and porttipahta in the northern

589

boreal zone in Finland. Global Biogeochem. Cycles 2002, 16 (1), 3–1–3–17.

590

(39) Liu, S.; Zhang, L.; Liu, Q.; Zou, J. Fe(III) fertilization mitigating net global warming

591

potential and greenhouse gas intensity in paddy rice-wheat rotation systems in China.

592

Environ. Pollut. 2012, 164, 73–80.

593

(40) Xing, Y.; Xie, P.; Yang, H.; Ni, L.; Wang, Y.; Tang, W. Diel variation of methane

594

fluxes in summer in a eutrophic subtropical lake in China. J. Freshw. Ecol. 2004, 19 (4),

595

639–644.

596

(41) Delsontro, T.; Kunz, M. J.; Kempter, T.; Wüest, A.; Wehrli, B.; Senn, D. B. Spatial

597

heterogeneity of methane ebullition in a large tropical reservoir. Environ. Sci. Technol.

598

2011, 45 (23), 9866–9873. 29

ACS Paragon Plus Environment

Environmental Science & Technology

599

(42) Heathcote, A. J.; Downing, J. A. Impacts of eutrophication on carbon burial in

600

freshwater lakes in an intensively agricultural landscape. Ecosystems 2012, 15 (1), 60–70.

601

(43) Kelly, C. A.; Chynoweth, P. The Contributions of temperature and of the input of

602

organic matter in controlling rates of sediment methanogenesis. Limnol. Oceanogr. 1981,

603

26 (5), 891–897.

604

(44) Beaulieu, J. J.; Balz, D. A.; Birchfield, M. K.; Harrison, J. A.; Nietch, C. T.; Platz, M.

605

C.; Squier, W. C.; Waldo, S.; Walker, J. T.; White, K. M. Effects of an experimental

606

water-level drawdown on methane emissions from a eutrophic reservoir. Ecosystems 2018,

607

21 (4), 657–674.

608

(45) Attermeyer, K.; Flury, S.; Jayakumar, R.; Fiener, P.; Steger, K.; Arya, V.; Wilken, F.;

609

Van Geldern, R.; Premke, K. Invasive Floating Macrophytes reduce greenhouse gas

610

emissions from a small tropical lake. Sci. Rep. 2016, 6, 1–10.

611

(46) Zhu, D.; Wu, Y.; Chen, H.; He, Y.; Wu, N. Intense methane ebullition from open

612

water area of a shallow peatland lake on the eastern Tibetan Plateau. Sci. Total Environ.

613

2016, 542, 57–64.

614

(47) Schmid, M.; Ostrovsky, I.; McGinnis, D. F. Role of gas ebullition in the methane

615

budget of a deep subtropical lake: What can we learn from process-based modeling?

616

Limnol. Oceanogr. 2017, 62 (6), 2674–2698.

617

(48) Matthews, C. J.; St Louis, V. L.; Hesslein, R. H. Comparison of three techniques used

618

to measure diffusive gas exchange from sheltered aquatic surfaces. Environ. Sci. Technol.

619

2003, 37(4), 772.

620

(49) Beaulieu, J. J.; Shuster, W. D.; Rebholz, J. A. Controls on gas transfer velocities in a

621

large river. J. Geophys. Res. Biogeosciences 2012, 117 (2), 1–13.

622

(50) Prairie, Y. T.; del Giorgio, P. A. A new pathway of freshwater methane emissions and

623

the putative importance of microbubbles. Inl. Waters 2013, 3 (3), 311–320. 30

ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31

Environmental Science & Technology

624

(51) McGinnis, D. F.; Kirillin, G.; Tang, K. W.; Flury, S.; Bodmer, P.; Engelhardt, C.;

625

Casper, P.; Grossart, H. P. Enhancing surface methane fluxes from an oligotrophic lake:

626

Exploring the microbubble hypothesis. Environ. Sci. Technol. 2015, 49 (2), 873–880.

627

31

ACS Paragon Plus Environment