Monitoring Effect of SO2 Emission ... - ACS Publications

Aug 3, 2017 - the capital city) reached the highest value in 2006 (0.86 Mt year. −1. ) ...... controlling the chemical response of streams of wester...
0 downloads 0 Views 1MB Size
Subscriber access provided by Warwick University Library

Article

Monitoring effect of SO2 emission abatement on recovery of acidified soil and stream water in southwest China Qian Yu, Ting Zhang, Xiaoxiao Ma, Ronghua Kang, Jan Mulder, Thorjørn Larssen, and Lei Duan Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b01147 • Publication Date (Web): 03 Aug 2017 Downloaded from http://pubs.acs.org on August 4, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37

Environmental Science & Technology

1

Monitoring effect of SO2 emission abatement on recovery of acidified

2

soil and stream water in southwest China

3 4 5

Qian Yu †, Ting Zhang†, Xiaoxiao Ma‡, Ronghua Kang§, Jan Mulder§, Thorjørn Larssen#, Lei Duan†,∇,*

6 7



Environment, Tsinghua University, Beijing 100084, China;

8 9 10

State Key laboratory of Environmental Simulation and Pollution Control, School of



State Grid Xingyuan Company Limited, Beijing 100761, China;

§

Department of Environmental Sciences, Norwegian University of Life Sciences, Box 5003, NO-1432 Ås, Norway.

11 12

#

Norwegian Institute for Water Research, Gaustadalleen 21, 0349 Oslo, Norway;

13



Collaborative Innovation Centre for Regional Environmental Quality, Tsinghua

14

University, Beijing 100084, China;

15 16

ABSTRACT: Following Europe and North America, East Asia has become a global

17

hotspot of acid deposition, with very high deposition of both sulfur (S) and nitrogen (N)

18

occurring in large areas in southwest and southeast China. Great efforts have been made

19

in reducing national emission of sulfur dioxide (SO2) since 2006 in China. However, the

20

total emission of nitrogen oxides (NOx) continued to increase until 2011. In order to

21

evaluate the effects of SO2 and NOx emission abatement on acid deposition and

22

acidification of soil and water, we monitored the chemical composition of throughfall,

23

soil water, and stream water from 2001 to 2013 in a small, forested catchment near

24

Chongqing city in Southwestern China. The deposition of S decreased significantly,

25

whereas N deposition increased in the recent years. This clearly showed the effect of

ACS Paragon Plus Environment

Environmental Science & Technology

26

SO2 abatement but not of NOx. Overall the rate of acid deposition was reduced.

27

However, there was delay in the recovery of soil and surface water from acidification,

28

probably due to desorption of previously stored sulfate (SO42-) and increase in nitrate

29

(NO3-) leaching from soil. The average acid input by N transformations has greatly

30

exceeded the H+ input directly by atmospheric deposition. The reversal of acidification

31

with an increase in pH of soil water, requires additional abatement of emissions of both

32

SO2 and NOx.

33

KEYWORDS: sulfur deposition, nitrogen deposition; acidification; sulfur adsorption;

34

nitrogen saturation; Chongqing

35 36

INTRODUCTION

37

Atmospheric acid deposition has been shown to affect the chemistry of soils and

38

surface waters and accelerate their acidification. Trends of soil water in forested

39

catchments and other ecosystems sensitive to acid deposition in the late 1980s and in the

40

1990s, were analyzed in many countries in Europe and North America.1~4 After the

41

Convention on Long-Range Transboundary Air Pollution to reduce emissions of

42

acidifying pollutants such as sulfur dioxide (SO2) and nitrogen oxides (NOx) was

43

implemented, significant recovery of acidified soil water and surface water in forest

44

ecosystems was reported.5~8

45

Since emission control in Europe and North America has taken effect, East Asia has

46

become the global hotspot of acid deposition. Elevated deposition of both sulfur (S) and

47

nitrogen (N) occurs in southwest and southeast China.9~11 Acid deposition was first

ACS Paragon Plus Environment

Page 2 of 37

Page 3 of 37

Environmental Science & Technology

48

reported in China in the late 1970s, with the lowest annual average pH of rainwater and

49

highest S deposition in southwest China, especially around Chongqing city.12~14 After a

50

30-year continuous increase, the emissions of SO2 and NOx have been reduced at the

51

national scale in recent years, following the implementation of the National SO2 Total

52

Emission Control in 2006 and the National NOx Total Emission Control in 2011,

53

respectively.15 However, emissions of NH3 have continued to increase because of

54

technical and economic difficulties in controlling these emissions. So far, there have

55

been few long-term monitoring reports of the effects of acid deposition on soils and

56

waters in China.16

57

Much of our understanding of the effect of acid deposition on soils and waters is

58

based on work in temperate regions of Europe and North America, where mineral

59

weathering is the main acid-buffering process in soils.17 Further extension of acid

60

deposition to Eastern and Southern Asia, and South America in recent years,13,18,19

61

warrants the collection of long-term monitoring data on soil and water acidification in

62

these regions. Although acid forest soils in tropical and subtropical regions,

63

characterized by slow mineral weathering and thus limited acid-buffering capacity, are

64

potentially sensitive to acidification,20,21 both atmospheric base cation deposition,22,23

65

denitrification,16 and S accumulation and adsorption24 may enhance the acid neutralizing

66

capacity of these soils.

67

Here, we use long-term monitoring data from a forested catchment near Chongqing

68

city in southwest China during 2001-2013 to examine the trends in soil and water

ACS Paragon Plus Environment

Environmental Science & Technology

69

acidification, and to assess the effect of recent reductions in national S and N emissions.

70

Accompanying the increasing importance of nitrogen oxides (NOx) and ammonia (NH3)

71

emissions, the acidifying effect of nitrogen (N) deposition has become a research and

72

management focus in recent years.25~27 The main objectives of this study were: (1) to

73

assess the effect of SO2 emission abatement on acid deposition, (2) to investigate the

74

response of soil water and stream water chemistry to changes in acid deposition, and (3)

75

to evaluate the contribution of NOx and NH3 emissions to soil and water acidification.

76 77

METHODOLOGY

78

Study site. Tieshanping (29°38′N, 106°43′E), one of the five long-term monitoring

79

catchments involved in the IMPACTS (Integrated Monitoring Program on Acidification

80

of Chinese Terrestrial Systems) project,9 is a 16-ha subtropical headwater catchment

81

about 25 km northeast from the center of Chongqing City in southwestern China (Fig.1).

82

As one of the biggest cities in southwest China, Chongqing has high emissions of both

83

SO2 and NOx. The SO2 emissions of Chongqing province (with an area of 8.24 Mha;

84

Chongqing as the capital city) reached the highest value in 2006 (0.86 Mt yr-1) and

85

decreased afterwards in response to the national SO2 emission abatement implemented

86

in 2006 (Fig. 2). The NOx emission increased until 2011, when the national NOx

87

abatement was implemented. Due to the relatively constant N fertilizer consumption in

88

Chongqing province since 2001, it is unlikely that the emission of NH3 changed much

89

during 2001-2013.

ACS Paragon Plus Environment

Page 4 of 37

Page 5 of 37

Environmental Science & Technology

90

Annual mean temperature and precipitation (1971 to 2000) at the nearest

91

meteorological station in Chongqing was 18.2°C and 1105 mm, respectively. The

92

bedrock of the catchment is sandstone. The forest soil is locally called Yellow Earth

93

(corresponding to a Haplic Acrisol in FAO, and representative for this part of China),

94

with a thin (< 1 cm) O horizon and low pH value (3.7~4.1 from the O horizon to the

95

lower B horizon).24 Based on earlier studies, the Acrisol at Tieshanping has a dense

96

argic B horizon at about 30 cm soil depth which reduced hydraulic conductivity, thus

97

creating interflow along the hill slope.28 The major biogeochemical processes

98

determining water transport and S and N transformations in soil are thus limited to the

99

upper soil horizons (0-30cm). The vegetation is Masson pine (Pinus massoniana Lamb.)

100

dominated forest, which is widely distributed in the humid subtropical areas of south and

101

southwest China.29 The annual average pH of precipitation at the site ranged from 4.0 to

102

4.2 during 2001–2004, with elevated deposition of S (10 keq·ha−1·yr−1), N (3

103

keq·ha−1·yr−1, about two-thirds as NH4+ and the rest as NO3-) and calcium (Ca) (6

104

keq·ha−1·yr−1).17 The high acid deposition and ambient air SO2 concentration caused

105

acidification of soil and stream water, and impacted forest health.30 In general, the study

106

site could present sensitive forested catchment to acidification in this region.

107

Sampling and chemical analysis. In the catchment of Tieshanping we studied two 10

108

m by 10 m sampling plots named C, and L (Fig. 1), each of which had four throughfall

109

collectors placed in the corners and a set of ceramic suction lysimeters placed close to

110

the center of the plots, at two different depths (about 3 and 30 cm) for soil water

ACS Paragon Plus Environment

Environmental Science & Technology

111

collection. The self-made throughfall collectors consisted of a 10-cm-diameter funnel to

112

collect throughfall and an opaque bottle to store water. Plot C is located at the foot of the

113

hillslope with an elevation of 560 m above sea level, near the stream, whereas plot L

114

(548 m in elevation) is located near the summit of the hillslope. Continuous sampling of

115

throughfall and soil water was performed weekly at both plots from 2001 to 2013 (Note

116

that in 2007, only throughfall was sampled and analyzed). Stream water was also

117

collected by fixed volume (500 mL) weekly (2001-2004, 2009-2013) at a dam (The dam

118

was moved from downstream of a pool to upstream in 2008), where a flow-measuring

119

flume was installed to record the stream water flux.

120

Every four weeks, the samples from all throughfall collectors, each lysimeter, and

121

stream were pooled into a bulk sample (monthly sample with volume proportion)

122

respectively for chemical analysis. All water samples were filtered (before filtration a

123

sub-sample was used for pH measurement) and analysed for major cations (Ca2+, Mg2+,

124

K+, Na+, and NH4+) and anions (NO3−, SO42−, Cl−, and F−) by ion chromatography.

125

Monomeric (mainly labile) aluminium (Ala) fraction of soil water and surface water was

126

determined by measuring the amount of 8-hydroxy-quinoline aluminium complex that

127

formed within 20 s.31 The organic fraction (Alo) of Ala was determined by passing the

128

sample through a pH-adjusted cation exchange column before 8-hydroxy-quinoline

129

complex formation. Inorganic monomeric Al (Ali) was calculated as the difference

130

between Ala and Alo.32 Al speciation, especially Al3+ concentration, was calculated by

131

the ALCHEMI program.33

ACS Paragon Plus Environment

Page 6 of 37

Page 7 of 37

Environmental Science & Technology

132

Data analysis. The flux of each ion was calculated from ion concentration and water

133

flux in each compartment. Since the water flux in soil was not directly measured in this

134

study, it was estimated on the basis of the Na+ balance assumption (The ratio of soil

135

water flux to throughfall water flux is assumed to equal the ratio of Na+ concentration in

136

throughfall to that in soil water). 24 It should be noted that Cl- is generally assumed to be

137

conservative, and thus can be used as a tracer to estimate soil water flux.34 However,

138

long-term monitoring in this catchment showed significant sink of Cl- in soil, as

139

indicated by similar or even lower Cl- concentration in soil water and stream water than

140

in throughfall water (Fig. S1 in Supplementary Information). In comparison with Cl-,

141

Na+ seemed more conservative in this catchment (Fig. S2). Considering the existence of

142

potential sources/sinks of Na+ in soil such as mineral weathering and cation exchange,

143

even extremely low in this catchment,24,35 the estimation of soil water flux based on the

144

real ratios of Na+ concentrations between throughfall and soil water might cause

145

uncertainty. Therefore, a fixed ratio of soil water flux to throughfall water flux was

146

assumed, 0.79 for S1 and 0.43 for S2, according to a fertilizing experiment by NaNO3

147

addition in the same catchment,24 where Na+ was more conservative than the untreated

148

soil because soil Na+ weathering or exchange could be ignored due to huge amount of

149

Na+ added.

150

The acid neutralizing capacity (ANC) was calculated as the sum of base cations (Ca2+

151

+ Mg2+ + K+ + Na+ + NH4+, in equivalent) minus the sum of acid anions (SO42- + NO3- +

152

Cl- + F-, in equivalent) based on the charge balance approach.36 The molar ratio of base

ACS Paragon Plus Environment

Environmental Science & Technology

153

cation to Al3+, Bc/Al, was calculated as the sum of base cation (Ca2+ + Mg2+ + K+, in

154

molar concentration) divided by the molar concentration of Al3+. In addition to H+

155

deposition, N transformations may contribute to soil acidification (H+N):37

156

H+N = (NH4+i –NH4+o) + (NO3-o – NO3-i)

157

where NH4+i and NO3-i were throughfall fluxes of NH4+ and NO3- respectively, and

158

NH4+o and NO3-o were the leaching fluxed of NH4+ and NO3- in soil water.

159

Statistically significant temporal changes in the chemical composition of throughfall

160

and soil water was analyzed by the flow-adjusted, seasonal Kendall test (SKT) and the

161

Sen slope.38-40 The analysis was done for two periods, 2001-2006 and 2008-2013,

162

respect to the SO2 emission control since 2006. Monthly data of pH, SO42-, NO3-, Ca2+,

163

and Al3+ concentration were used for analysis. In addition, post hoc analysis using

164

one-way ANOVA was applied to analyze differences in water chemistry in different

165

periods (2001-2006 and 2008-2013). Difference were considered significant at the p ≤

166

0.05 level.

167 168

RESULTS

169

Throughfall. The monthly concentration of major ions, such as SO42-, NO3-, NH4+,

170

Ca2+, and Al3+ in throughfall water, as well as the water amount and pH of throughfall,

171

show clear seasonal variation (Fig. 3). The pH (annual average ranging from 3.35 to

172

4.53) was high in summer but low in winter, which coincided with the seasonal variation

173

of precipitation (shown by throughfall water amount). Throughfall pH was only greater

ACS Paragon Plus Environment

Page 8 of 37

Page 9 of 37

Environmental Science & Technology

174

than 5.6 during a few months in summer. The concentrations of SO42-, NO3-, and Ca2+

175

were all higher in winter than in summer. In general, lower concentrations of SO42- and

176

NO3- due to the diluting effect of abundant rain in summer might lead to higher pH. The

177

NH4+ concentration showed not only peaks in winter, but also some peaks in summer,

178

which reflected high NH3 emission from agricultural field in the growing season.

179

The pH of throughfall showed a significant trend of decrease during 2001-2006

180

(negative Sen slope by SKT; at a rate of -0.0953 yr-1), but a significant increase during

181

2008-2013 (at a rate of 0.1792 yr-1) (Table 1). Several ions had significant

182

increasing/decreasing trends in concentration values. For example, until 2006, the SO42-

183

concentration in throughfall increased significantly (at 97.25 µeq⋅L-1⋅yr-1). However,

184

after 2008 the concentration of SO42- decreased (at a rate of -110.5 µeq⋅L-1⋅yr-1). In

185

contrast, the NO3- concentration showed an increasing trend during the whole period (at

186

12.59 µeq⋅L-1⋅yr-1); This rate became even higher in recent years. The concentration of

187

NH4+ and Ca2+ also showed increasing trends during the whole period. The faster rate of

188

SO42- concentration decreasing than the increasing rate of NO3- concentration after 2008

189

(while no significant trends in NH4+ and Ca2+) indicated decreasing trend of throughfall

190

water acidity, coincided with the increasing trend of pH.

191

Soil water and surface water. Temporal variations of monthly concentration of

192

SO42- and Ca2+ in soil water and stream water coincided with those in throughfall (with

193

some delay), but with smaller amplitude (especially in S2 and stream water) (Fig. 3).

194

Due to nitrification or uptake, little NH4+ occurred in soil water (especially in S2). As a

ACS Paragon Plus Environment

Environmental Science & Technology

195

result of NO3- formation by nitrification, more frequent variation of NO3- concentration

196

was found in S1 than in throughfall, according to the peaks of DIN (NH4+ + NO3-) in

197

throughfall. The pH of soil water showed opposite trend to throughfall, i.e., lower in

198

summer but higher in winter. The stream water chemistry was even more stable than that

199

of soil water in the bottom layer (S2).

200

Unlike throughfall, soil water and stream water showed a decreasing trend of pH in

201

the whole period (at -0.0171, -0.0369 and -0.0877 yr-1 in S1, S2 and SW respectively)

202

(Table 1). The annual average pH of stream water was in the range of 3.96-4.85, which

203

was greater than that of soil water (S1, 3.39-4.03;S2, 3.59-4.30). Unlike the decreasing

204

trend of SO42- concentration in throughfall after 2008, the SO42- concentration of soil

205

water and stream water showed no significant trend during the same period, even

206

increasing trend during the whole study period (at 23.53 and 40.64 µeq⋅L-1⋅yr-1 in S2 and

207

SW respectively). The concentration of NO3- and NH4+ in soil water and stream water

208

showed also an increasing trend, similar as in throughfall. The increase in Ca2+

209

concentration in stream water might neutralize the increase in SO42- concentration (i.e.,

210

the changing rates of Ca2+ and SO22- were similar), while decreasing pH coincided with

211

increasing NO3-. It should be noted that more variable stream water chemistry during

212

2009-2013 occurred than during 2001-2004 period (Fig. 3), because the dam for stream

213

water monitoring was moved from the downstream to upstream of a small pond.

214

Accompanying the decrease in pH, the Al3+ concentration in soil water showed an

215

increasing trend during the whole period (at 20.80 and 32.26 µeq⋅L-1⋅yr-1 in S1 and S2

ACS Paragon Plus Environment

Page 10 of 37

Page 11 of 37

Environmental Science & Technology

216

respectively) (Table 1). However, there was significant decline in the saturation of the

217

soil solution with respect to Al(OH)3, as can be seen from the significant decline in log

218

Q (defined as log Al + 3 pH)41 (Fig. 3).

219

Fluxes. Yearly fluxes of S, N, Ca + Mg, and H+ in throughfall, soil water, and stream

220

water showed also variation during 2001-2013 (Fig.4). The fluxes of S and Ca + Mg

221

showed generally a decreasing trend along the hydrological continuums (i.e., from

222

throughfall to soil water, then to stream water). In contrast, the fluxes of total inorganic

223

nitrogen (TIN) in top soil water (S1) were larger than that in throughfall. The average H+

224

depostion based on throughfall was significantly lower during 2008-2013 than during

225

2001-2006 (one-way ANOVA; Table S1), which coincided with the significantly higher

226

average pH of throughfall during 2008-2013 than during 2001-2006. However, the

227

average H+ flux in soil water (S2) was significantly higher during 2008-2013 than

228

during 2001-2006, while there is no significant difference in H+ flux in stream water.

229

There were no significant differences in fluxes of SO42-, NO3- and Ca2+ between the two

230

periods in stream water, with only exception of NO3- depostion, which was significantly

231

higher during 2008-2013 than during 2001-2006 (one-way ANOVA; Table S1).

232 233

DISCUSSION

234

Change of acid deposition. The significant declining trend of S depositon after 2008

235

found in this study is among the first to document a successful SO2 emission abatement

236

in China. During the monitoring period, the annual S deposition (based on throughfall)

ACS Paragon Plus Environment

Environmental Science & Technology

237

increased to maximum values in 2004 and 2005 (>14 keq·ha-1·yr-1; Fig. 4). After 2008,

238

the S deposition showed a decreasing trend (significant negative Sen slope; Table 1),

239

reaching 7.7 keq·ha-1·yr-1 in 2013, which is smaller than in 2001. This temporal pattern

240

in throughfall deposition of S, with its decline after 2008, coincides with a decrease in

241

SO2 emission in Chongqing province (Fig. 2). In contrast, the annual N deposition

242

continued to increase, reaching 5.3 keq·ha-1·yr-1 in 2010 - 2012, which is twice as much

243

as in 2001 (2.5 keq·ha-1·yr-1) (Fig. 4). This trend also coincides with the fast increasing

244

NOx emission before 2011 (Fig. 2), which resulted in larger increase in NO3- than NH4+

245

in N deposition, and decrease in NH4+/NO3- ratio (Fig. 4). The net effect of the two N

246

species on throughfall was thus decrease in pH. The effect of decreasing NOx emission

247

after 2011 has not been seen yet.

248

Although S and N deposition at the site were large compared to values in Europe and

249

North America in 1980s,9 the precipitation was not exceptionally acidic, with annual

250

average pH of throughfall of 3.79, which is within the range of values reported for

251

typically acidified regions in North America and Europe (3.5~4.5).1~3,42 The annual

252

average pH of deposition would have been around 3.0 without the high deposition of

253

Ca2+ (and other base cations), which reflected the importance of the alkaline dust,

254

mainly from anthropogenic sources such as cement production and iron steel industry in

255

China.23 The annual Ca2+ deposition showed an increasing trend (due to increasing

256

anthropogenic emission of Ca-content dust),43 reaching a value as high as 8.2

257

keq·ha-1·yr-1. The annual Mg2+ deposition fluxes (see Fig. S3 for monthly variation in

ACS Paragon Plus Environment

Page 12 of 37

Page 13 of 37

Environmental Science & Technology

258

concentration) amounted to about 23% of Ca2+ deposition (Fig. 4), and also contributed

259

greatly to ANC. As a result of decreasing S deposition and increasing Ca deposition

260

since 2008, the H+ deposition decreased after reaching maximum values of about 4

261

keq·ha-1·yr-1 (Fig. 4).

262

S and N budet. The overall budget (i.e., the difference between throughfall

263

deposition flux and stream water output) indicated considerable retention of S, N and Ca

264

in the catchment (Table 2). The S output with soil water in the top layer (S1) was lower

265

than input by deposition, especially in the earlier years. This highlights a significant

266

SO42- sink in the upper soil. In a recent study we found that SO42- adsorption in the acid

267

soils was the most likely mechanism removing SO42- from soil water.24 The sulfate

268

adsorbed pool was 16.3-28.5 keq ha-1 in the upper 30 cm soil layer, and there might be

269

larger adsorption capacity in deeper soil layer.16 In addition, incorporation to organic

270

matter might be another important sulfate rentention mechanism, especially in the topper

271

soil with higer organic carbon content. The organic S pool was estimated to be 11 keq

272

ha-1 according to the ratio of organic S to organic C

273

soil.24 Since the S deposition was averaged to 11.5 keq ha-1 yr-1 during 2001-2006 (Table

274

S1), this SO42- sink in top soil (3.33 keq ha-1 yr-1 in Table2) amounted to about 29% of

275

the S deposition. As SO42- adsorption is a process of ligand exchange with OH-,44 it may

276

account for an increase in the ANC with about 3 keq ha-1 yr-1, which is larger than the

277

annual H+ depostion (2.60 keq ha-1 yr-1). Adsorption and desorption may also be

278

important to control short-term fluctuations of the SO42- flux, such as the observed

48

and the organic C content in

ACS Paragon Plus Environment

Environmental Science & Technology

279

year-to-year variation. When S deposition decreased after 2008, the S sink in S1

280

decreased to 1.30 keq ha-1 yr-1, which accounted for only 12% of the S deposition

281

(averagely 11.0 keq ha-1 yr-1) during 2008-2013 (Table S1). In addition, much more S

282

retention occurred below the upper 30 cm soil layer (Table 2). Together with SO42-

283

adsorption in the deeper soil, SO42- reduction to sulfides in the wet soils of the

284

groundwater discharge zone (GDZ) may contribute significantly to S sink in the

285

catchment, which need be tested in the future.

286

In soil water, the flux of NH4+-N was very low, while that of NO3--N was large and

287

even exceeded the sum of the input of NH4+ and NO3- in throughfall (Fig. 4). This

288

indicates that most NH4+ derived from throughfall deposition and mineralization of soil

289

organic matter was nitrified. The large leaching of NO3--N from the top soil layer (S1)

290

indicates large N internal cycle, as illustrated by the rapid mineralization of annual

291

litterfall.45

292

decreased due to immobilization, nutrient uptake, and denitrification16. As a result, the

293

upper 30cm soil layer showed a slight N sink (averaged 1.38 keq ha-1 yr-1; according to

294

Table 2), which indicated little less N flux in soil solution (S2) than in throughfall

295

deposition. In addition, there were larger variations of N fluxes from year to year in the

296

soil solution; the fluxes were larger in the wet years, which can be attributed to

297

mineralized N from litterfall and soil organic matter at high moisture conditions (driven

298

by rainfall). It may be concluded that N saturation46 is extreme in hillslope soils at

299

Tieshanping with all N deposition leached. This is also shown in the studies of soil N

From the top layer (S1) to the lower layer (S2), the N flux in soil water

ACS Paragon Plus Environment

Page 14 of 37

Page 15 of 37

Environmental Science & Technology

300

status,45 N mass balances,24 natural 15N abundances47,48 and 15N tracer49 in the catchment.

301

Similar to S, a high N retention occurs in the groundwater discharge zone, that is, after

302

percolation water has left the upper 30 cm of soil but before the water leaves the

303

catchment as stream water at the weir (Table 2). Denitrification, especially in GDZ, may

304

be the key reason of N loss.16 Since loss of nitrate by denitrification provides ANC,

305

stream water was much less acidified than soil water, by nitrogen deposition.

306

Like S, Ca2+ sink mainly occurred in the upper 30 cm of the soil (Table 2). In China,

307

Ca2+ deposition is a major source of Ca2+ in soil water.16,23 The soil Ca2+ sink was too

308

large to be explained by plant uptake, since the Ca2+ uptake by forest trees at the site was

309

estimated to only 0.25-1.10 keq ha-1 yr-1.50 A possible mechanism might be the building

310

up of cation exchange capacity (CEC),16 possibly due to the increasing content of soil

311

organic matter caused by N and Ca deposition.27,51

312

Acidification and recovery. The concentration of SO42- and NO3- in soil water in the

313

catchment was very high, even one order higher than that seen in Europe in 1980s and

314

1990s.1~4 Although the concentration of base cations (especially Ca2+) was also very

315

high, it was not high enough to compensate for the negative charge associated with

316

SO42- and NO3-. Thus, the pH of soil water and stream water was relatively low.

317

The pH and the concentration (and deposition) of SO42- in throughfall showed

318

acidification of rainwater during 2001-2006, and recovery (not only alkalinisation

319

because of significant decrease in S deposition) after 2008. However, pH of soil water

320

did not show a significant increase in recent years, even with significant increase in pH

ACS Paragon Plus Environment

Environmental Science & Technology

321

of stream water (Table 1 and Table S1). The delayed recovery of forest soil and surface

322

water from acidification was also widely found in Europe and North America, and

323

attributed to the release of previously stored SO42- from soil, coupled to the leaching of

324

Al3+ and base cations.3,4,52~54 The large adsorption of SO42- by soil on oxide surfaces

325

buffered soil acidity thus limiting the concentration of H+ in soil water, especially in the

326

earlier years (2001-2006). This occurred even with much NH4+ nitrified to produce

327

acidity. In recent years, accompanying the decrease in SO42- concentration in soil water

328

due to the decrease in S deposition, SO42- adsorption largely decreased and the

329

previously adsorbed sulfate outflowed from the upper 30 cm soil layer, which led to an

330

increase in ANC sink and thus a H+ source in this layer (Table 2). However, the total

331

SO42- sink remained unchanged because of the increase in SO42- adsorption in deeper

332

soil layer (Table 2).

333

Associated with the increasing N deposition, N transformations played an

334

increasingly important role in the acidification of soils and waters. The average acid

335

input by N transformation was 3.03 keq ha-1 yr-1, exceeding the H+ input directly by

336

atmospheric deposition (1.91 keq ha-1 yr-1) (Fig. 5). The exceedance was larger than the

337

previous findings.25 Since the S flux in soil water and stream water did not show

338

significantly increasing trend in recent years, the decrease in pH in stream water (Table

339

1) might be caused by increase in NO3- leaching. Actually, since NOx emission

340

abatement only began in 2011, its effect on acidification was not expected to see.

ACS Paragon Plus Environment

Page 16 of 37

Page 17 of 37

Environmental Science & Technology

341

The Al3+ concentration of soil water did not change with the decrease of pH,

342

providing further evidence that gibbsite equilibrium does not present a good model to

343

describe the Al solubility.24,55,56 There was significant decline in log Q (defined as log

344

Al + 3 pH), while Al3+ concentration increasing (Table 1). The results may support the

345

earlier hypothesis that in strongly acidified forest soils complexion by solid phase

346

organics controls the solubility of Al, even in mineral soil layers, relatively low in

347

organic C.41,57

348

The Bc/Al in soil water is commonly used as a risk indicator for forest decline due to

349

acidification,58 and 1.0 is widely used as a critical limit in critical load assessments in

350

Europe.59 In this study, the molar ratios of Bc/Al in soil water of S2 were averaged to

351

2.42 at Tieshanping, much higher than the critical limit used. However, possible damage

352

to forest by acidification had occurred at this site due to high Al concentration. 32 This

353

implies that a higher critical limit should be applied in subtropical China.60 In addition,

354

no significant trend of increasing Bc/Al was found after 2008 than before 2006.

355 356

POLICY IMPLICATION

357

The long-term monitoring on throughfall, soil water and stream water during 2001 to

358

2013 in Tieshanping catchment, representing catchments sensitive to acidification in

359

southwest China, showed clear effect of emission abatement on recovery of acidified

360

soil and stream water. Coincided with the decrease in SO2 emission in Chongqing

361

province after 2006, a significantly decreasing trend of S deposition was detected during

ACS Paragon Plus Environment

Environmental Science & Technology

362

2008-2013. In contrast, N deposition continued increasing in accordance with the

363

increase in NOx (until 2011) and NH3 emission. As a net effect, acid deposition had

364

been reduced, and the acidification of soil water and surface water had leveled off. The

365

trends of deposition in Tieshanping from 2001-2013 could represent the regional or even

366

national trends. It should be noted that high SO42- adsorption in the acid soil, together

367

very high deposition of base cations, especially Ca2+, provided large ANC to the

368

ecosystem. Those might be the main reasons that the soil water acidification in China

369

was not very serious in comparison with that in Europe and North America, even with

370

much higher S and N deposition. However, the recovery of acidified soil water and

371

surface water was delayed, due probably to desorption of previously stored sulfate

372

(SO42-) and enhanced leaching of nitrate (NO3-) from soil. Since soil acidification is still

373

very serious, shown by quite low pH in soil water, future emission abatement of both

374

SO2 and NOx should be strengthened.

375 376

ASSOCIATED CONTENT

377

Supporting Information

378

Additional information on the water chemistry during the periods and the calculation of

379

soil water flux are presented in a supplementary document.

380

AUTHOR INFORMATION

381

Corresponding Author

382

* Phone: +86 10 6278 3758; fax: +86 10 6277 3597; e-mail: [email protected]

383

Notes

ACS Paragon Plus Environment

Page 18 of 37

Page 19 of 37

Environmental Science & Technology

384

The authors declare no competing financial interest.

385 386

ACKNOWLEDGEMENT

387

The authors are grateful for the financial support of the National Natural Science

388

Foundation of China (21221004), the State Environmental Protection Public Welfare

389

Project of China (201209001) and the Research Council of Norway (209696/E10) to

390

carry out this study.

391 392

REFERENCE

393

(1) Meesenburg, H.; Meiwes, K. J.; Rademacher P. Long term trends in atmospheric

394

deposition and seepage output in northwest German forest ecosystems. Water, Air,

395

Soil Pollut. 1995, 85, 611-616.

396

(2) Kvaalen, H.; Solberg, S.; Clarke, N.; Torp, T.; Aamlid, D. Time series study of

397

concentrations of SO42- and H+ in precipitation and soil waters in Norway. Environ

398

Pollut. 2002, 117, 215-224.

399

(3) Moffat, A.J.; Kvaalen, H.; Solberg, S.; Clarke, N. Temporal trends in throughfall

400

and soil water chemistryat three Norwegian forests, 1986-1997. Forest Eco Manag.

401

2002, 168, 15-28.

402

(4) Fölster, J.; Bringmark, L.; Lundin, L. Temporal and spatial variations in soilwater

403

chemistry at three acid forest sites. Water, Air, Soil Pollut. 2003, 146, 171-195.

ACS Paragon Plus Environment

Environmental Science & Technology

404

(5) Lawrence, G. B.; Hazlett, P. W.; Fernandez, I. J.; Ouimet, R.; Bailey, S. W.;

405

Shortle, W. C.; Smith, K. T.; Antidormi, M. R., Declining Acidic Deposition

406

Begins Reversal of Forest-Soil Acidification in the Northeastern US and Eastern

407

Canada. Environ. Sci. Technol. 2015, 49, (22), 13103-13111

408

(6) Garmo, Ø. A.; Skjelkvåle, B. L.; Wit, H. A. d.; Colombo, L.; Curtis, C.; Fölster, J.;

409

Hoffmann, A.; Hruška, J.; Høgåsen, T.; Jeffries, D. S., Trends in Surface Water

410

Chemistry in Acidified Areas in Europe and North America from 1990 to 2008.

411

Water, Air, Soil Pollut. 2014, 225, (3), 1-14.

412

(7) Strock, K. E.; Nelson, S. J.; Kahl, J. S.; Saros, J. E.; Mcdowell, W. H., Decadal

413

trends reveal recent acceleration in the rate of recovery from acidification in the

414

northeastern U.S. Environ. Sci. Technol. 2014, 48, (9), 4681-9.

415

(8) Sutherland, J. W.; Acker, F. W.; Bloomfield, J. A.; Boylen, C. W.; Charles, D. F.;

416

Daniels, R. A.; Eichler, L. W.; Farrell, J. L.; Feranec, R. S.; Hare, M. P., Brooktrout

417

Lake case study: biotic recovery from acid deposition 20 years after the 1990 Clean

418

Air Act Amendments. Environ. Sci. Technol. 2015, 49, (5), 2665-74.

419 420

(9) Larssen, T.; Lydersen, E.; Tang, D. G.; et al. Acid rain in China. Environ. Sci. Technol. 2006, 40, 418-425.

421

(10) Zhao, Y.; Duan, L.; Xing, J.; Larssen, T.; Nielsen, C.; Hao, J. M. Soil acidification

422

in China: Is controlling SO2 emissions enough? Environ. Sci. Technol. 2009, 43,

423

8021-8026.

ACS Paragon Plus Environment

Page 20 of 37

Page 21 of 37

Environmental Science & Technology

424

(11) Zhao, B.; Wang, S. X.; Dong, X. Y.; Wang, J. D.; Duan, L.; Fu, X.; Hao, J. M.; Fu,

425

J. Environmental effects of the recent emission changes in China: Implications for

426

particulate matter pollution and soil acidification. Environ. Res Lett. 2013, 8,

427

024031.

428

(12) Zhao, D. W.; Sun, B. Z. Air pollution and acid rain in China. Ambio 1986, 15, 2-5.

429

(13) Galloway, J. N.; Zhao, D. W.; Xiong, J. L.; Likens. G. E. Acid rain: a comparison

430

of China, United States and a remote area. Science 1987, 236, 1559-1962.

431

(14) Hao, J. M.; Xie, S. D.; Duan, L.; Ye X. M. Acid deposition and ecosystem

432

sensitivity in China. In: Bashkin V., Park S. (eds.). Acid Deposition and Ecosystem

433

Sensitivity in East Asia. Nova Science Publishers, New York, 1998, pp267-311.

434

(15) MEP (Ministry of Environmental Protection) of China. Notice of the State Council

435

on Issuing the Twelfth Five-Year Plan for National Environmental Protection. 2011.

436

http://zfs.mep.gov.cn/fg/gwyw/ 201112/t20111221_221570.htm. (in Chinese)

437

(16) Larssen, T.; Duan, L.; Mulder, J. Deposition and leaching of sulfur, nitrogen and

438

calcium in four forested catchments in China: implications for acidification.

439

Environ. Sci. Technol. 2011, 45, 1192-1198.

440

(17) Sverdrup, H.; Warfvinge, P. Weathering of primary silicate minerals in the natural

441

soil environment in relation to a chemical weathering model. Water, Air, Soil Pollut.

442

1988, 38, 387-408.

443 444

(18) Figueredo, D. V. Influence of calcareous soil particulates on acid rain: Belo Horizonte Metropolitan Region, Brazil. Ambio 1999, 28, 514-518.

ACS Paragon Plus Environment

Environmental Science & Technology

445 446 447 448

(19) Bhattacharya, A.; Mudgal, R.; Taneja, A. Acid deposition and critical load analysis in Agra, India. J Hazard Mater 2004, 106, 157-160. (20) Duan, L.; Xie, S. D.; Zhou, Z. P.; Hao, J. M. Critical loads of acid deposition on soil in China. Water, Air, Soil Pollut. 2000, 118, 35-51.

449

(21) Lu, X. K.; Mao, Q. G.; Mo, J. M.; Gilliam, F. S.; Zhou, G. Y.; Luo, Y. Q.; Zhang,

450

W.; Huang, J. Divergent responses of soil buffering capacity to long-term N

451

deposition in three typical tropical forests with different land-use history. Environ.

452

Sci. Technol. 2015, 49, 4072-4080.

453

(22) Zhao, Y.; Duan, L.; Larssen, T.; Hu L. H.; Hao, J. M. Simultaneous assessment of

454

depositions of base cations, sulfur and nitrogen using an extended critical load

455

function for acidification. Environ. Sci. Technol. 2007, 41, 1815-1820.

456

(23) Zhao, Y.; Duan, L.; Lei, Y.; Xing, J.; Nielsen, C. P.; Hao, J. M. Will PM control

457

undermine China’s efforts to reduce soil acidification? Environ. Pollut. 2011, 159,

458

2726-2732.

459

(24) Huang, Y. M.; Kang, R. H.; Mulder, J.; Zhang, T.; Duan, L. Nitrogen saturation,

460

soil acidification and ecological effects in a subtropical pine forest on acid soil in

461

Southwest China. J Geophys Res. 2015, 120, 2457–2472.

462

(25) Chen, X. Y.; Mulder, J.; Wang, Y. H.; Zhao, D. W.; Xiang, R. J. Atmospheric

463

deposition, mineralization and leaching of nitrogen in subtropical forested

464

catchments, south China. Environ Geochem Hlth 2004, 26, 179-186.

ACS Paragon Plus Environment

Page 22 of 37

Page 23 of 37

Environmental Science & Technology

465 466

(26) Liu, X. J.; Duan, L.; Mo, J. M.; et al. Nitrogen deposition and its ecological impact in China: an overview. Environ. Pollut. 2011, 159, 2251-2264.

467

(27) Lu, X. K.; Mao, Q. G.; Gilliam, F. S.; Luo, Y. Q.; Mo, J. M. Nitrogen deposition

468

contributes to soil acidification in tropical ecosystems. Glob. Change Biol. 2014, 20,

469

3790-3801.

470

(28) Sørbotten, L. Hill slope unsaturated flowpaths and soil moisture variability in a

471

forested catchment in southwest China. Master Dissertation, Norwegian University

472

of Life Sciences, Aas, 2011.

473

(29) Wu, Z. Y. China Vegetation. Science Press, Beijing, 1980. (in Chinese)

474

(30) Wang, Y. H.; Solberg, S.; Yu, P. T. Assessments of tree crown condition of two

475

Masson pine forests in the acid rain region in south China. Forest Ecol Manag,

476

2007, 242, 530-540.

477 478 479 480 481 482 483 484

(31) Barnes, R. B. The determination of specific forms of aluminium in natural water. Chem Geol 1975, 15, 177-191. (32) Driscoll, C. T. A procedure for the fractionation of aqueous aluminium in dilute acidic waters. Int J Environ An Ch 1984, 16, 267-283. (33) Schecher, W. D.; Driscoll, C. T. An evaluation of uncertainty associated with aluminium equilibrium calculations. Water Resour Res 1987, 23, 525-534. (34) Drever, J. I.; Clow, D. W. Weathering rates in catchment. Rev Mineral 1995, 31, 463-483.

ACS Paragon Plus Environment

Environmental Science & Technology

485

(35) Duan, L.; Liu, J.; Xin, Y.; Larssen, T., Air-pollution emission control in China:

486

Impacts on soil acidification recovery and constraints due to drought. Sci. Total

487

Environ. 2013, 463, 1031-1041.

488 489 490 491 492 493 494 495

(36) Reuss, J. O.; Christophersen, N.; Seip, H. M., A critique of models for freshwater and soil acidification. Water, Air, Soil Pollut. 1986, 30, (30), 909-930. (37) Van Breemen, N.; Mulder, J.; Driscoll, C. T. Acidification and alkalinization of soil. Plant Soil 1983, 75, 283–308. (38) Hirsch, R. M.; Slack, J. R. A nonparametric test for seasonal data with serial dependence. Water Resour Res 1984, 20, 727-732. (39) Sen, P. K. Estimations of the regression coefficient base on Kendall’s Tau. J Am Stat Assoc 1968, 63, 1379-1389.

496

(40) Hirsch, R. M.; Slack, J. R.; Smith, R. A. Techniques of trend analysis for monthly

497

water quality data. Water Resour Res 1982, 18, 107-121.Mulder, J.; Stein, A. The

498

solubility of aluminum in acidic forest soils: Long-term changes due to acid

499

deposition. Geochim Cosmochim Ac 1994, 58, 85-94.

500 501

(41) Likens, G. E.; Driscoll, C. T.; Buso, D. C. Long-term effects of acid rain: Response and recovery of a forest ecosystem. Science 1996, 272, 244-246.

502

(42) Lei, Y.; Zhang, Q.; He, K. B.; Streets, D. G. Primary anthropogenic aerosol

503

emission trends for China, 1990–2005. Atmos Chemd Phys 2011, 11, 931-954.

504

(43) McBride, M. B. Environmental Chemistry of Soils. Oxford University Press, New

505

York, 1994.

ACS Paragon Plus Environment

Page 24 of 37

Page 25 of 37

Environmental Science & Technology

506 507 508 509

(44) Chen, X. Y.; Mulder, J. Indicators for nitrogen status and leaching in subtropical forest ecosystems, South China. Biogeochemistry 2007, 82, 165-180. (45) Aber, J. D.; Nadelhoffer K. J.; Steudler P.; Melillo J. M. Nitrogen saturation in northern forest ecosystems. BioScience 1989, 39, 378-386.

510

(46) Zhu, J.; Mulder, J.; Bakken, L.; Dorsch, P., The importance of denitrification for

511

N2O emissions from an N-saturated forest in SW China: results from in situ N-15

512

labeling experiments. Biogeochemistry 2013, 116, (1-3), 103-117.

513

(47) Yu, L.; Zhu, J.; Mulder, J.; Dörsch, P., Multiyear dual nitrate isotope signatures

514

suggest that N‐saturated subtropical forested catchments can act as robust N sinks.

515

Glob. Change Biol. 2016, 22, (11), 3662-3674.

516

(48) Yu, L.; Kang, R.; Mulder, J.; Zhu, J.; Dorsch, P., Distinct fates of atmogenic NH4+

517

and NO3- in subtropical, N-saturated forest soils. Biogeochemistry 2017, 133, (3),

518

279-294.

519

(49) Duan, L.; Huang, Y. M.; Hao, J. M.; Xie, S. D.; Hou, M. Vegetation uptake of

520

nitrogen and base cations in China and its role in soil acidification. Sci. Total

521

Environ. 2004, 330, 187-198.

522

(50) Luo, Y.; Duan, L.; Xu, G. Y.; Hao, J. M. Inhibition of mercury release from forest

523

soil by high atmospheric deposition of Ca2+ and SO42-. Chemosphere 2015, 134,

524

113-119.

ACS Paragon Plus Environment

Environmental Science & Technology

525

(51) Kopacek, J.; Hejzlar, J.; Porcal, P.; Posch, M., Sulphate leaching from diffuse

526

agricultural and forest sources in a large central European catchment during

527

1900-2010. Sci. Total Environ. 2014, 470, 543-550.

528

(52) Prechtel, A.; Alewell, C.; Armbruster, M.; et al. Response of sulphur dynamics in

529

European catchments to decreasing sulphate deposition. Hydrol Earth Syst Sc 2001,

530

5, 311-325.

531

(53) Robison, A. L.; Scanlon, T. M.; Cosby, B. J.; Webb, J. R.; Galloway, J. N., Roles

532

of sulfate adsorption and base cation supply in controlling the chemical response of

533

streams of western Virginia to reduced acid deposition. Biogeochemistry 2013, 116,

534

(1-3), 119-130.

535

(54) Gruba, P.; Mulder J. Relationship between aluminum in soils and soil water in

536

mineral horizons of a range of acid forest soils. Soil Sci Soc Am J 2008, 72,1150–

537

1157.

538

(55) Gruba, P.; Mulder J.; Brożek S. Modelling the pH dependency of dissolved calcium

539

and aluminium in O, A and B horizons of acid forest soils. Geoderma, 2013, 206,

540

85–91.

541

(56) Berggren, D.; Mulder, J. The role of organic matter in controlling aluminium

542

solubility in acidic mineral soils horizons. Geochim Cosmochim Ac 1995, 59,

543

4167-4180.

544 545

(57) Cronan, C. S.; Grigal, D. F. Use of calcium/aluminum ratios as indicators of stress in forest ecosystems. J Environ Qual 1995, 24, 209-226.

ACS Paragon Plus Environment

Page 26 of 37

Page 27 of 37

Environmental Science & Technology

546 547

(58) Sverdrup, H.; Warfvinge, P.; Frogner, T.; Håøya, A. O.; Johansson, M.; Andersen, B. Critical loads for forest soils in the Nordic countries. Ambio 1992, 21, 348-355.

548

(59) Sverdrup, H., Warfvinge, P. Effect of Soil Acidification on Growth of Trees and

549

Plants as Expression by the (Ca+Mg+K)/Al Ratio. Report 2: 1993, Lund University,

550

Lund, 1993.

551

ACS Paragon Plus Environment

Environmental Science & Technology

552

Page 28 of 37

Tables

553 554

Table 1 Sen slopes of throughfall (CTF), soil water (S1 and S2), and stream water (SW)

555

concentrationa Component

Period

pH

SO42-

NO3-

NH4+

Ca2+

Al3+

logQ

CTF

2001~2013

0.0073**

n.s.

12.59**

6.232*

16.69*

-

-

2001~2006

-0.0953**

97.25**

9.756**

n.s.

n.s.

-

-

2008~2013

0.1792**

-110.5**

12.78**

n.s.

n.s.

-

-

2001~2013

-0.0171**

n.s.

30.94**

2.681**

n.s.

20.80**

-0.0798**

2001~2006

-0.0073**

n.s

29.70*

n.s

n.s

43.80**

n.s.

2008~2013

0.0042*

n.s

n.s

n.s

n.s

n.s.

n.s.

2001~2013

-0.0369**

23.53**

26.18**

n.s

n.s.

32.26**

-0.0963**

2001~2006

-0.0263**

44.70**

35.54**

n.s.

n.s

32.57**

n.s.

2008~2013

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

n.s.

2001~2013

-0.0877**

40.64**

9.452**

2.580**

40.82**

-

-

2001~2006

-

-

-

-

-

-

-

2008~2013

-0.0713*

n.s.

n.s.

n.s.

n.s.

-

-

S1

S2

SW

556

a

Analyzed by seasonal Kendall test (SKT), with change of concentration in µeq⋅L-1⋅yr-1, and change of pH in yr-1;

557

* Trend significant at P