1
Decreased soil cation exchange capacity across northern China’s
2
grasslands over the last three decades
3 4 5 6
Kai Fang1,2, Dan Kou1,2, Guanqin Wang1,2, Leiyi Chen1, Jinzhi Ding1,2, Fei Li1,2, Guibiao Yang1,2, Shuqi Qin1,2, Li Liu1,2, Qiwen Zhang1,2 and Yuanhe Yang1,2*
7
1
8
Chinese Academy of Sciences, Beijing 100093, China
9
2
University of Chinese Academy of Sciences, Beijing 100049, China
11
*
Corresponding author: Dr. Yuanhe Yang, tel.: + 86 10-6283 6638, fax: + 86
12
10-6283 6632, E-mail:
[email protected]
State Key Laboratory of Vegetation and Environmental Change, Institute of Botany,
10
13 14 15 16
Revised manuscript submitted to Journal of Geophysical Research:
17
Biogeosciences
18
21th- Sep.-2017
19 20 21 22 23
Manuscript information: 31 pages, 1 table, 5 figures and 1 supplementary dataset.
24
Key Points:
25
Decreases in topsoil cation exchange capacity were observed across northern
26 27 28 29 30
China’s grasslands over the last three decades. The decreases were detected on both the Tibetan Plateau and the Inner Mongolian Plateau, among various grassland types. The decreases were attributed to soil carbon loss, soil desertification, soil acidification, and extreme precipitation.
31 1
32
Abstract
33
Cation exchange capacity (CEC) helps soils hold nutrients and buffer pH, making it
34
vital for maintaining basic function of terrestrial ecosystems. However, little is known
35
about the temporal dynamics of CEC over broad geographical scales. In this study, we
36
used Random Forest method to compare historical CEC data from the 1980s with new
37
data from the 2010s across northern China’s grasslands. We found that topsoil CEC in
38
the 2010s was significantly lower than in the 1980s, with an overall decline of about
39
14%. Topsoil CEC decreased significantly in alpine meadow, alpine steppe, meadow
40
steppe, and typical steppe by 11%, 20%, 27% and 9% respectively. Desert steppe was
41
the only ecosystem type which experienced no significant change. CEC was
42
positively related to soil carbon content, silt content, and mean annual precipitation,
43
suggesting that the decline was potentially associated with soil organic carbon loss,
44
soil degradation, soil acidification, and extreme precipitation across northern China’s
45
grasslands since the 1980s. Overall, our results demonstrate topsoil CEC loss due to
46
environmental changes, which may alter the vegetation community composition and
47
its productivity and thus trigger grassland dynamics under a changing environment.
48 49
Keywords: environmental change, geochemistry, grassland ecosystems, soil carbon
50
loss, soil desertification, soil acidification, extreme precipitation
2
51
1. Introduction
52
Cation exchange capacity (CEC) is usually expressed as the number of moles of
53
exchangeable cations adsorbed by electrostatic force per unit mass soil. CEC plays an
54
important role in holding soil nutrients [Brady and Weil, 2002] and buffering soil pH
55
changes [Xu et al., 2012; Luo et al., 2015], which makes it vital for maintaining the
56
structure and function of terrestrial ecosystems. CEC can affect ecosystem structure
57
and function by controlling the supply of exchangeable cations (i.e., Ca2+, Mg2+, K+,
58
Na+, Al3+, and Fe3+) in the soil [Lucas et al., 2011; Mueller et al., 2012]. Given that
59
different plant species require different concentrations and ratios of Ca2+, Mg2+, K+
60
and Na+, changes in these cations can alter plant biodiversity [Chen et al., 2013].
61
Besides being essential to plant growth [Likens et al., 1998; McLaughlin and Wimmer,
62
1999], the availability of these nutrients affects plant tolerance to drought, frost, and
63
pathogens [DeHayes et al., 1999; Demchik and Sharpe, 2000], therefore affecting
64
primary productivity and overall structure of terrestrial ecosystems. In contrast,
65
increased Al3+ and Fe3+ could limit vegetation productivity because of its toxicity
66
[Lieb et al., 2011]. CEC can also maintain the stability of terrestrial ecosystems by
67
buffering changes in soil acidity [Bowman et al., 2008], since soil pH decreases could
68
alter biological activity and cation supply and have negative impacts on terrestrial
69
ecosystems [Kirk et al., 2010]. Overall, in the context of environmental changes, a
70
deep understanding of the temporal dynamics of CEC is crucial for predicting the
71
responses of ecosystem structure and function, and in guiding policies that maintain
72
the stability of terrestrial ecosystems. 3
73 74
Given the importance of CEC, studies about the effects of environmental changes on
75
CEC have gradually accumulated in recent decades [Watmough and Dillon, 2003;
76
Högberg et al., 2006; Fissore et al., 2007; Ruiz Sinoga et al., 2012; Xu et al., 2012; Lu
77
et al., 2015]. It has been reported that CEC exhibited positive correlations with
78
2:1-type clay minerals [Xu et al., 2012] and organic matter content [Ruiz Sinoga et al.,
79
2012; Xu et al., 2012]. Due to the close relationships between soil organic matter and
80
climatic parameters, CEC has also been observed to be positively related with
81
precipitation [Ruiz Sinoga et al., 2012] but negatively related with mean annual
82
temperature [Fissore et al., 2007]. These studies have generated broad understanding
83
of spatial variation and environmental drivers of CEC over broad geographic scales,
84
but little is known about temporal CEC dynamics at large spatial scales. Long-term
85
site monitoring has demonstrated cation loss from soils due to soil acidification,
86
harvesting, and land-use change [Watmough and Dillon, 2003; Högberg et al., 2006;
87
Lu et al., 2015]. However, these observations are limited to the site scale and it is
88
unknown if they are representative of larger landscape or continental trends.
89
Considering that the driving factors for CEC dynamics (e.g., climatic factors, edaphic
90
factors, and human activity) exhibit substantial variations in space and time, a
91
comprehensive study about the temporal dynamics of CEC along broad environmental
92
gradients is needed.
93 94
China’s grasslands are distributed across large geographic scales, making them an 4
95
ideal ecosystem to explore the temporal dynamics of CEC for three reasons. First,
96
along with the geographic gradient, grassland types are diverse (including both alpine
97
and temperate grasslands) and environmental parameters such as precipitation and
98
edaphic variables exhibit a wide range [Liao and Jia, 1996]. Second, since the 1980s,
99
northern China has experienced significant environmental changes, which mainly
100
include the following: (1) Overgrazing and cultivation [Bridges and Oldeman, 1999;
101
Yu et al., 2012] have led to soil desertification [Yang et al., 2007; Mu et al., 2013] and
102
pasture degradation [Dai et al., 2011], thereby decreasing grassland productivity
103
[Babel et al., 2014] and organic carbon storage [Xie et al., 2007; Dai et al., 2011]
104
across the study area; (2) Elevated atmospheric acid deposition has induced soil
105
acidification across northern China, with a decrease in soil pH of 0.63 units over the
106
past two decades [Yang et al., 2012b]; (3) The frequency and intensity of extreme
107
precipitation has increased across northern China [Zhai et al., 2005], specifically,
108
precipitation in the desert region of China has increased more than 40% since the
109
beginning of the 21st century [Xu et al., 2010]. Third, soil physical and chemical
110
properties during the 1980s were well documented in historical soil inventory data
111
[Office of Soil Census of Inner Mongolia, 1994; Tibet Land Management Bureau,
112
1994; Office of Agriculture Resource Layouting of Qinghai, 1997]. Together, these
113
geographic conditions, environmental changes, and historical inventory make possible
114
an investigation of the temporal dynamics of CEC over the last three decades.
115 116
In this study, we investigated changes in CEC across northern China’s grasslands 5
117
from the 1980s to 2010s. We collected surface soil samples from 251 sites on the
118
Tibetan Plateau and Inner Mongolian Plateau from 2013-2015, comparing these data
119
with 125 soil profiles measured during the 1980s, obtained from China’s Second
120
National Soil Inventory. We combined these two datasets of CEC with the related
121
environmental parameters, and used Random Forest method to explore the changes of
122
CEC between the two periods. We hypothesized that topsoil CEC across northern
123
China’s grasslands would show a significant decrease due to the substantial
124
environmental changes over the past three decades.
125 126
2. Materials and Methods
127
2.1. Study Area
128
Our study area mainly contained alpine grasslands on the Tibetan Plateau and
129
temperate grasslands on the Inner Mongolian Plateau (Figure 1), whose longitudes
130
and latitudes are 80.10°-121.93° E and 28.07°-49.62° N, respectively. During the
131
2010s, the mean annual temperature (MAT) varied from -3.06 °C to 8.13 °C, and the
132
mean annual precipitation (MAP) varied from 111.51 mm to 719.28 mm. The
133
dominant grassland types are alpine meadow (AM), alpine steppe (AS), meadow
134
steppe (MS), typical steppe (TS) and desert steppe (DS) [Chinese Academy of
135
Sciences, 2001]. The dominant species are Kobresia pygmaea and K. tibetica in the
136
AM, Stipa purpurea and Festuca ovina in the AS, S. baicalensis and Leymus
137
chinensis in the MS, S. grandis and S. krylovii in the TS, and S. klemenzii, and S.
138
breviflora in the DS [Ma et al., 2010]. According to the classification of the World 6
139
Reference Base for Soil Resources 2006, the soil types related to the five grassland
140
types were defined as Cambisols (in AM and AS), Chernozems (in MS), Kastanozems
141
(in TS) and Calcisols (in DS) [IUSS Working Group WRB, 2007].
142 143
2.2. Historical observations during the 1980s
144
China’s Second National Soil Inventory during 1979-1989 recorded 74 and 51
145
representative grassland soil profiles on the Tibetan Plateau and Inner Mongolian
146
Plateau, respectively (Figure 1). From the 1980s data, we obtained CEC and other
147
edaphic properties (e.g. SOC content, pH, and silt content) in the A horizon of 125
148
typical profiles, with a median depth of 15 cm [Office of Soil Census of Inner
149
Mongolia, 1994; Tibet Land Management Bureau, 1994; Office of Agriculture
150
Resource Layouting of Qinghai, 1997]. According to the historical records, CEC was
151
determined by hydrochloric acid neutralization titration, after extraction with a mixed
152
solution consisting of 1 M ammonium acetate and 0.005 M EDTA (adjusted to pH 7.0
153
for acidic soil and pH 8.5 for alkaline soil, respectively), and subsequent distillation of
154
ammonium [Chang, 1979; Sparks et al., 1996]. SOC content was determined by
155
ferrous sulfate titration, after potassium dichromate oxidation [Nelson and Sommers,
156
1982]. Soil pH was measured by potentiometry, using a water soil ratio of 2.5:1
157
[Koncki et al., 1992]. Silt content (soil texture) was measured by the sieve-hydrometer
158
method, based on the speed of particle sedimentation at multiple time steps [Di
159
Stefano et al., 2010]. Grassland types were determined by the dominant species
160
related to the local soil types. Meanwhile, based on the available detailed 7
161
administrative address, longitude and latitude at these sites were determined from a
162
digital map developed by the National Administration of Surveying, Mapping, and
163
Geoinformation (http://map.tianditu.com/map/index.html). Historical MAT and MAP
164
at these sites were extracted using Kriging interpolation based on atmospheric
165
temperature and precipitation data from 120 meteorological stations on the Tibetan
166
Plateau and the Inner Mongolian Plateau from 1982-1989 (http://data.cma.cn/).
167 168
2.3. Field inventory during the 2010s
169
To characterize the current status of grassland soil CEC, we investigated 251 sites
170
(173 on the Tibetan Plateau [Ding et al., 2016] and 78 on the Inner Mongolian Plateau)
171
throughout northern China’s grasslands during July and August of 2013-2015. The
172
sampling sites during the 2010s covered extensive climatic gradients and major
173
grassland types (Figure 1). Specifically, latitude and longitude were acquired using
174
GPS (M-GIS T20, China), and grassland types were determined by the dominant
175
species present. At each site, we set up five 1 m×1 m quadrats at every corner and in
176
the center of a 10 m×10 m square. We further selected three quadrats along a
177
diagonal line of the square to take topsoil (0-10 cm) samples. Soil samples were
178
collected in plastic bags and sent to the laboratory.
179 180
We analyzed soil CEC and other edaphic properties in the laboratory. First, composite
181
soil samples (equally mixed subsamples from three samples from each site) were
182
air-dried and sieved (< 2 mm), and fine roots were manually removed. To avoid the 8
183
potential uncertainties caused by measuring method, soil CEC, SOC content, pH, and
184
silt content were all measured by the same methods as in the 1980s. MAT and MAP at
185
these sites were extracted in the same method but from meteorological data of
186
2010-2015.
187 188
2.4. Data Analyses
189
To examine the dynamics of soil CEC from 1980s to 2010s, it is necessary to have
190
time-paired data for each sampling site. This method has been used to examine soil
191
organic [Yu et al., 2009; Yang et al., 2014] and inorganic carbon [Yang et al., 2012a]
192
dynamics, as well as soil pH changes [Yang et al., 2012b]. Because the site locations
193
in the two sampling periods were inconsistent, we used the Random Forest model to
194
estimate CEC values for the missing points from each sampling period. Random
195
Forest extends Classification and Regression Tree analyses to improve prediction
196
accuracy without excessive computation [Breiman, 2001a, 2001b; Wiesmeier et al.,
197
2011; Ziegler and Konig, 2014; Were et al., 2015]. To reduce the uncertainties caused
198
by the uneven distribution of sampling sites during the two periods, we selected a
199
grassland transect (including 91 among the 125 sites in the 1980s and 189 among the
200
251 sites in the 2010s) from the central-southern Tibetan Plateau to the eastern Inner
201
Mongolian Plateau (Figure 1) to examine the CEC dynamics. Briefly, the dataset of
202
measured CEC was split into two parts: a training dataset used to construct the
203
Random Forest model (90% of the data points) and a test dataset used to assess the
204
performance of the model (10% of the data points). Based on the historical soil 9
205
inventory training dataset from the 1980s and related environmental variables (i.e.,
206
longitude, latitude, grassland type, MAT, MAP, SOC, pH, and silt content), we
207
constructed a Random Forest model to predict CEC for those sites without
208
measurement in the 1980s (189 sites that surveyed during the 2010s). By doing so, we
209
obtained a dataset consisting of 189 data points including measured soil CEC from the
210
2010s and corresponding estimated data for the 1980s (detailed procedure in Figure 2).
211
Similarly, using the training dataset from the current survey during the 2010s and the
212
corresponding environmental data, we created another Random Forest model to
213
simulate CEC during the 2010s for those sites without measurements from the 2010s
214
(91 sites measured in the 1980s). Thereafter, we obtained another dataset consisting of
215
91 data points including measured CEC from 1980s and estimated data from the
216
2010s (Figure 2). By combining these data sets, we finally generated a time-paired
217
dataset across the 280 sites with either measured or estimated data for the two
218
sampling periods for each sampling site. We used Pearson correlation coefficient (r)
219
and root-mean-squared error (RMSE) to test the accuracy of the model predictions.
220
We repeated the Random Forest training 1000 times, and the averages were treated as
221
the predicted CEC values. All the Random Forest Model computations were
222
conducted using R software [R Development Core Team, 2012]. Based on the
223
time-paired dataset obtained through the Random Forest method, we used a paired
224
t-test to explore the temporal dynamics of CEC. We also examined the relationships
225
between CEC dynamics and environmental variables using linear regression model.
226
All the statistical analyses were performed using SPSS 20.0 (IBM SPSS, Chicago, 10
227
Illinois, USA), with a significance level of 0.05.
228 229
3. Results
230
3.1. Reliable Estimation by the Random Forest Model
231
To test the accuracy of the model predictions, we examined the correlation between
232
the predicted and measured CEC values. Our results showed that during the 1980s, r
233
and RMSE between the testing data set and the measured values were 0.90 and 5.41
234
cmol kg-1, respectively (Figure 3a). Similarly, r and RMSE during the 2010s were
235
0.93 and 3.48 cmol kg-1, respectively (Figure 3b). These results showed high
236
correlation and low estimation error, indicating a reliable prediction ability for the
237
Random Forest model.
238 239
3.2. Temporal Dynamics of CEC across Northern China’s Grasslands
240
Topsoil CEC across northern China’s grasslands during the 2010s was significantly
241
lower than the 1980s (P < 0.001; Figure 4a). According to the results of 280 sampling
242
sites within the transect, topsoil CEC significantly declined by about 14% (Table 1).
243
However, substantial variations existed among the various grassland types (i.e., AM,
244
AS, MS, TS and DS). Topsoil CEC did not exhibit a significant change in DS perhaps
245
because of its lower initial values, but decreased by 11%, 20%, 27% and 9% in AM,
246
AS, MS and TS, respectively (Figure 4b; Table 1).
247 248
4. Discussion
249
There has been significant loss of topsoil CEC across northern China’s grasslands 11
250
over the last three decades. A decline in base cations over a longer time period was
251
also detected in North America [Watmough and Dillon, 2003], Europe [Högberg et al.,
252
2006], Oceania [Hartemink, 1998; Curtin et al., 2015], Africa [Jaiyeoba, 2003] and
253
Asia [Ali et al., 1997; Lu et al., 2015], indicating that this phenomenon is widespread.
254
The decrease of CEC observed in this study is within the range of previous studies,
255
varying from 8% to 50% [Ali et al., 1997; Hartemink, 1998; Jaiyeoba, 2003; Curtin et
256
al., 2015]. Such a pattern could be due to the following four phenomena. First, a
257
reduction in SOC content can lead to a decrease in CEC. Overgrazing and cultivation
258
in past decades [Bridges and Oldeman, 1999; Yang et al., 2007; Yu et al., 2012; Mu et
259
al., 2013] has increasingly degraded grasslands [Dai et al., 2011] and causing declines
260
in grassland productivity [Babel et al., 2014]. Consequently, the SOC stock of China’s
261
grasslands decreased by 3.56 Pg C (1 Pg = 1015 g), dropping from approximately
262
37.71 to 34.15 Pg C from the 1980s to the 2000s, with particularly marked changes in
263
Tibet, Qinghai, and Inner Mongolia [Xie et al., 2007]. The decrease in SOC has a
264
profound influence on CEC for two reasons. One is that the dissociation of some
265
functional groups (e.g., carboxyl and phenolic hydroxyl) decreases with the decrease
266
in SOC. As a result, the amount of negative charge possessed by humus is reduced,
267
weakening the complexation of the humus with exchangeable base cations [Brady and
268
Weil, 2002; Chapin et al., 2011]. The other is that cation-bridging (e.g., Ca2+, Fe3+,
269
Al3+) between clay minerals and humus can decline due to a reduction in SOC,
270
making humus-clay mineral complexes unstable [Brady and Weil, 2002; von Lutzow
271
et al., 2006; Mueller et al., 2012]. Consequently, the protection for SOM against 12
272
decomposition is weakened and causes the decrease in CEC. Consistent with this
273
inference, our results showed that CEC was positively correlated with SOC, and its
274
variations were mostly explained by SOC in both the 1980s (TP: r2 = 0.72, P < 0.001;
275
IM: r2 = 0.71, P < 0.001; Figure 5a) and 2010s (TP: r2 = 0.89, P < 0.001; IM: r2 =
276
0.84, P < 0.001; Figure 5e). Many studies also have found a positive relationship
277
between CEC and SOC [Tůma et al., 2011; Lu et al., 2014; Gruba and Mulder, 2015].
278 279
Second, the CEC reduction can also be caused by soil desertification. It has been
280
reported that soils in arid and semi-arid regions were more likely to be degraded
281
[Bridges and Oldeman, 1999]. Both the Tibetan Plateau and Inner Mongolian Plateau,
282
which are located in arid and semi-arid regions, have been suffering the effects of
283
climate warming, human disturbance, overgrazing, cultivation, and pikas, all of which
284
can result in soil degradation and desertification [Yang et al., 2007; Xue et al., 2009;
285
Xue et al., 2011; Yu et al., 2012; Qin et al., 2015]. It has been reported that lots of
286
Inner Mongolian grasslands have undergone degradation or desertification from 1985
287
to 2000 [Mu et al., 2013]. Once soil desertification occurs, wind erosion alters soil
288
texture, blowing away fine particles (e.g. clay and silt) and leaving coarse and barren
289
sand [Chepil, 1957; Hoffmann et al., 2008; Amundson et al., 2015]. Moreover, with
290
the decrease in particle contents, the surface areas provided for exchangeable cations
291
become smaller [Brady and Weil, 2002]. Furthermore, soil desertification causes a
292
decrease in soil productivity and SOC content [Larney et al., 1998; Yan et al., 2005].
293
All of these processes will reduce CEC, along with soil desertification. Our data 13
294
support the above conclusions, with CEC positively correlated with silt content
295
(1980s TP: r2 = 0.12, P < 0.01; IM: r2 = 0.71, P < 0.001; Figure 5b; 2010s TP: r2 =
296
0.25, P < 0.001; IM: r2 = 0.68, P < 0.001; Figure 5f), but negatively associated with
297
sand content (1980s TP: r2 = 0.12, P < 0.01; IM: r2 = 0.62, P < 0.001; Figure 5c;
298
2010s TP: r2 = 0.14, P < 0.001; IM: r2 = 0.67, P < 0.001; Figure 5g).
299 300
Third, the CEC reduction can still be induced by soil acidification. It has been
301
reported that enhanced atmospheric nitrogen deposition [Liu et al., 2013] and sulfur
302
deposition [Zhao et al., 2009] resulted in soil acidification across China’s grasslands
303
[Yang et al., 2012b]. Sustained acidification leads to a lower soil pH, which can be
304
buffered by carbonate and base cations. It has been hypothesized that when the soil
305
pH is over 7.5, the change is mainly buffered by carbonate [Bowman et al., 2008;
306
Yang et al., 2012a, b], while the exchangeable base cations play a lesser role. When
307
the soil pH is in the range of 4.5-7.5, it is primarily the exchangeable cations that
308
buffer pH change [Bowman et al., 2008]. Exchangeable ions including Ca2+, Mg2+, K+
309
and Na+ are adsorbed on the surface of the soil colloid, making CEC value higher.
310
With the increases in acidity, H+ gradually occupies the exchangeable sites. At the
311
same time, base cations disassociate and are more likely to be leached by rainfall,
312
resulting in decreased CEC [Watmough and Dillon, 2003; Lapenis et al., 2004; Lieb et
313
al., 2011; Chen et al., 2013; Lu et al., 2014]. This classic hypothesis is supported by
314
experimental studies on long-term dynamics of CEC and base cations at the local
315
scale [Watmough and Dillon, 2003; Högberg et al., 2006; Lu et al., 2015]. To be 14
316
specific, results from a manipulative experiment also showed that soil pH and base
317
saturation significantly increased when acid addition stopped [Högberg et al., 2006].
318
Therefore, we argue that acidification mainly affects CEC when the pH was low (e.g.
319
in AM and MS).
320 321
Fourth, the decrease in CEC is possibly related to extreme precipitation. In recent
322
decades, extreme precipitation has been observed more frequently [IPCC, 2007;
323
Moldan et al., 2012]. Recent studies indicated that the frequency and intensity of
324
extreme precipitation were gradually increasing across western China and the Tibetan
325
Plateau [Zhai et al., 2005; Yang et al., 2008; Xu et al., 2010], making soil base cations
326
leached. Normally, precipitation is positively correlated with SOC [Ruiz Sinoga et al.,
327
2012], which is strongly associated with CEC [Gruba and Mulder, 2015]. Normally,
328
precipitation is positively correlated with SOC [Ruiz Sinoga et al., 2012], which is
329
strongly associated with CEC [Gruba and Mulder, 2015]. In both periods, our results
330
showed that CEC increased along the precipitation gradient (1980s TP: r2 = 0.12, P
Mg2+ > K+ > Na+, and so, if CEC is reduced, base cations
347
will be lost in different degrees. As a result, species that rely on a certain cation may
348
decrease or even become extinct on long time scales. In addition, in alpine meadow
349
and meadow steppe environments with a low pH, soil acidification caused by nitrogen
350
and sulfur deposition leads to the decrease in CEC, eventually releasing Al3+ and Fe3+
351
into the soil [Bowman et al., 2008]. Consequently, species that are sensitive to Al3+ or
352
Fe3+ may decrease. Second, ecosystem function across northern China’s grasslands
353
may be weakened when CEC is reduced. The decrease in CEC and the increase in
354
Al3+ and Fe3+ weaken the abilities of the nutrient (i.e., Ca2+, Mg2+, K+ and Na+) supply
355
and plant tolerance to drought, frost or pathogen attacks [DeHayes et al., 1999;
356
Demchik and Sharpe, 2000]. Thus, to some degree, it may potentially decrease
357
vegetation productivity [Likens et al., 1998; McLaughlin and Wimmer, 1999]. To
358
avoid the ecological consequences mentioned above, it is vital to reduce the
359
interference of human activities in these ecosystems and strictly control the intensity 16
360
of grazing and cultivation [Mu et al., 2013]. Reducing the emissions of SO2 and NOx
361
is also essential to relieve soil acidification and environment aggravation caused by
362
continuous CEC loss [Lucas et al., 2011].
363 364
Acknowledgments
365
Associated data are available in supporting information. We are grateful to the
366
members of the IBCAS Sampling Campaign Teams for their assistance in field
367
investigations. This work was supported by the National Natural Science Foundation
368
of China (31670482 and 41371213), Key Research Program of Frontier Sciences,
369
Chinese Academy of Sciences (QYZDB-SSW-SMC049), Chinese Academy of
370
Sciences-Peking University Pioneer Cooperation Team, and the Thousand Young
371
Talents Program.
372 373
References
374
Ali, M. M., S. M. Saheed, D. Kubota, T. Masunaga, and T. Wakatsuki (1997), Soil
375
degradation during the period 1961-1995 in Bangladesh .2. Selected chemical
376
characters, Soil Sci. Plant Nutr., 43(4), 879-890.
377
Amundson, R., A. A. Berhe, J. W. Hopmans, C. Olson, A. E. Sztein, and D. L. Sparks
378
(2015), Soil and human security in the 21st century, Science, 348(6235), doi:
379
10.1126/science.1261071.
380 381
Babel, W., et al. (2014), Pasture degradation modifies the water and carbon cycles of the
Tibetan
highlands,
Biogeosciences,
11(23),
6633-6656, 17
382
doi:10.5194/bg-11-6633-2014.
383
Bowman, W. D., C. C. Cleveland, Ĺ. Halada, J. Hreško, and J. S. Baron (2008),
384
Negative impact of nitrogen deposition on soil buffering capacity, Nat. Geosci.,
385
1(11), 767-770, doi:10.1038/ngeo339.
386 387 388 389 390 391 392
Brady, N. C., and R. R. Weil (2002), The Nature and Properties of Soils, Prentice Hall, New Jersey. Breiman,
L.
(2001a),
Random
5-32,
Bridges, E. M., and L. R. Oldeman (1999), Global assessment of human-induced soil
doi:10.1080/089030699263212.
398
45(1),
199-215, doi:10.1214/ss/1009213726.
394
397
Learn.,
Breiman, L. (2001b), Statistical modeling: The two cultures, Statistical Science, 16(3),
degradation,
396
Mach.
doi:10.1023/a:1010933404324.
393
395
forests,
Arid
Soil
Res.
Rehabil.,
13(4),
319-325,
Chang, Z. (1979), Routine analysis programs and methods in the soil census, Xinjiang Agricultural Sciences, 2, 31-33. Chapin, F. S., P. A. Matson, and P. M. Vitousek (2011), Principles of Terrestrial Ecosystem Ecology, Springer, New York.
399
Chen, D., Z. Lan, X. Bai, J. B. Grace, and Y. Bai (2013), Evidence that
400
acidification-induced declines in plant diversity and productivity are mediated by
401
changes in below-ground communities and soil properties in a semi-arid steppe,
402
J. Ecol., 101(5), 1322-1334, doi:10.1111/1365-2745.12119.
403
Chepil, W. S. (1957), Sedimentary characteristics of dust storms .1. sorting of 18
404 405 406
wind-eroded soil material, Am. J. Sci., 255(1), 12-22. Chinese Academy of Science (2001), Vegetation Atlas of China, Science Press, Beijing.
407
Curtin, D., P. M. Fraser, and M. H. Beare (2015), Loss of soil organic matter
408
following cultivation of long-term pasture: effects on major exchangeable
409
cations
410
doi:10.1071/sr14173.
and
cation
exchange
capacity,
Soil
Res.,
53(4),
377-385,
411
Dai, F. Q., Z. A. Su, S. Z. Liu, and G. C. Liu (2011), Temporal variation of soil
412
organic matter content and potential determinants in Tibet, China, Catena, 85(3),
413
288-294, doi:10.1016/j.catena.2011.01.015.
414
DeHayes, D. H., P. G. Schaberg, G. J. Hawley, and G. R. Strimbeck (1999), Acid rain
415
impacts
on
calcium
nutrition
and
forest
health
-
Alteration
of
416
membrane-associated calcium leads to membrane destabilization and foliar
417
injury in red spruce, Bioscience, 49(10), 789-800, doi:10.2307/1313570.
418
Demchik, M. C., and W. E. Sharpe (2000), The effect of soil nutrition, soil acidity and
419
drought on northern red oak (Quercus rubra L.) growth and nutrition on
420
Pennsylvania sites with high and low red oak mortality, For. Ecol. Manage.,
421
136(1-3), 199-207, doi:10.1016/s0378-1127(99)00307-2.
422
Di Stefano, C., V. Ferro, and S. Mirabile (2010), Comparison between grain-size
423
analyses using laser diffraction and sedimentation methods, Biosyst. Eng., 106(2),
424
205-215, doi:10.1016/j.biosystemseng.2010.03.013.
425
Ding, J., et al. (2016), The permafrost carbon inventory on the Tibetan Plateau: a new 19
426
evaluation using deep sediment cores, Glob. Change Biol., 22(8), 2688-2701,
427
doi:10.1111/gcb.13257.
428
Fissore, C., C. P. Giardina, R. K. Kolka, C. C. Trettin, G. M. King, M. F. Jurgensen, C.
429
D. Barton, and S. D. McDowell (2007), Temperature and vegetation effects on
430
soil organic carbon quality along a forested mean annual temperature gradient in
431
North
432
doi:10.1111/j.1365-2486.2007.01478.x.
America,
Glob.
Change
Biol.,
14(0),
193-205,
433
Gruba, P., and J. Mulder (2015), Tree species affect cation exchange capacity (CEC)
434
and cation binding properties of organic matter in acid forest soils, Sci. Total
435
Environ., 511, 655-662, doi:10.1016/j.scitotenv.2015.01.013.
436
Hartemink, A. E. (1998), Soil chemical and physical properties as indicators of
437
sustainable land management under sugar cane in Papua New Guinea, Geoderma,
438
85(4), 283-306, doi:10.1016/s0016-7061(98)00048-2.
439
Hoffmann, C., R. Funk, Y. Li, and M. Sommer (2008), Effect of grazing on wind
440
driven carbon and nitrogen ratios in the grasslands of Inner Mongolia, Catena,
441
75(2), 182-190, doi:10.1016/j.catena.2008.06.003.
442
Högberg, P., H. B. Fan, M. Quist, D. Binkley, and C. O. Tamm (2006), Tree growth
443
and soil acidification in response to 30 years of experimental nitrogen loading on
444
boreal
445
doi:10.1111/j.1365-2486.2005.01102.x.
446 447
forest,
Glob.
Change
Biol.,
12(3),
489-499,
IPCC (2007), Climate Change 2007: Synthesis Report. Intergovernment Panel on Climate Change, pp73. 20
448 449 450
IUSS Working Group WRB (2007), World reference base for soil resources 2006. First update 2007, World Soil Resources Report No. 103. (FAO: Rome). Jaiyeoba, I. A. (2003), Changes in soil properties due to continuous cultivation in
451
Nigerian
semiarid
Savannah,
452
doi:10.1016/s0167-1987(02)00138-1.
Soil
Tillage
Res.,
70(1),
91-98,
453
Kirk, G. J. D., P. H. Bellamy, and R. M. Lark (2010), Changes in soil pH across
454
England and Wales in response to decreased acid deposition, Glob. Change Biol.,
455
16(11), 3111-3119, doi:10.1111/j.1365-2486.2009.02135.x.
456 457
Knapp, A. K., et al. (2008), Consequences of more extreme precipitation regimes for terrestrial ecosystems, Bioscience, 58(9), 811-821, doi:10.1641/b580908.
458
Koncki, R., P. Leszcyński, A. Hulanicki, and S. Gła̧b (1992), Urea sensors based on
459
glass pH electrodes with physically immobilized urease, Anal. Chim. Acta,
460
257(1), 67-72, doi:http://dx.doi.org/10.1016/0003-2670(92)80151-V.
461
Lapenis, A. G., G. B. Lawrence, A. A. Andreev, A. A. Bobrov, M. S. Torn, and J. W.
462
Harden (2004), Acidification of forest soil in Russia: From 1893 to present, Glob.
463
Biogeochem. Cycle, 18(1), doi:10.1029/2003gb002107.
464
Larney, F. J., M. S. Bullock, H. H. Janzen, B. H. Ellert, and E. C. S. Olson (1998),
465
Wind erosion effects on nutrient redistribution and soil productivity, J. Soil Water
466
Conserv., 53(2), 133-140.
467 468 469
Liao, G. F., and Y. L. Jia (1996), Rangeland Resources of China, China Science and Technology Press, Beijing. Lieb, A. M., A. Darrouzet-Nardi, and W. D. Bowman (2011), Nitrogen deposition 21
470
decreases acid buffering capacity of alpine soils in the southern Rocky
471
Mountains, Geoderma, 164(3-4), 220-224, doi:10.1016/j.geoderma.2011.06.013.
472
Likens, G. E., et al. (1998), The biogeochemistry of calcium at Hubbard Brook,
473 474 475
Biogeochemistry, 41(2), 89-173, doi:10.1023/a:1005984620681. Liu, X., et al. (2013), Enhanced nitrogen deposition over China, Nature, 494(7438), 459-462, doi:10.1038/nature11917.
476
Lu, X., Q. Mao, F. S. Gilliam, Y. Luo, and J. Mo (2014), Nitrogen deposition
477
contributes to soil acidification in tropical ecosystems, Glob. Change Biol.,
478
20(12), 3790-3801, doi:10.1111/gcb.12665.
479
Lu, X., Q. Mao, J. Mo, F. S. Gilliam, G. Zhou, Y. Luo, W. Zhang, and J. Huang (2015),
480
Divergent responses of soil buffering capacity to long-term N deposition in three
481
typical tropical forests with different land-use history, Environ. Sci. Technol.,
482
49(7), 4072-4080, doi:10.1021/es5047233.
483
Lucas, R. W., J. Klaminder, M. N. Futter, K. H. Bishop, G. Egnell, H. Laudon, and P.
484
Högberg (2011), A meta-analysis of the effects of nitrogen additions on base
485
cations: Implications for plants, soils, and streams, For. Ecol. Manage., 262(2),
486
95-104, doi:10.1016/j.foreco.2011.03.018.
487
Luo, W. T., et al. (2015), Contrasting pH buffering patterns in neutral-alkaline soils
488
along a 3600 km transect in northern China, Biogeosciences, 12(23), 7047-7056,
489
doi:10.5194/bg-12-7047-2015.
490
Ma, W., J.-S. He, Y. Yang, X. Wang, C. Liang, M. Anwar, H. Zeng, J. Fang, and B.
491
Schmid (2010), Environmental factors covary with plant diversity-productivity 22
492
relationships among Chinese grassland sites, Glob. Ecol. Biogeogr., 19(2),
493
233-243, doi:10.1111/j.1466-8238.2009.00508.x.
494
McLaughlin, S. B., and R. Wimmer (1999), Tansley Review No. 104 - Calcium
495
physiology and terrestrial ecosystem processes, New Phytol., 142(3), 373-417,
496
doi:10.1046/j.1469-8137.1999.00420.x.
497
Moldan, F., J. Hruska, C. D. Evans, and M. Hauhs (2012), Experimental simulation of
498
the effects of extreme climatic events on major ions, acidity and dissolved
499
organic carbon leaching from a forested catchment, GAyenrdsjon, Sweden,
500
Biogeochemistry, 107(1-3), 455-469, doi:10.1007/s10533-010-9567-6.
501
Mu, S. J., Y. Z. Chen, J. L. Li, W. M. Ju, I. O. A. Odeh, and X. L. Zou (2013),
502
Grassland dynamics in response to climate change and human activities in Inner
503
Mongolia, China between 1985 and 2009, Rangeland J., 35(3), 315-329,
504
doi:10.1071/rj12042.
505
Mueller, K. E., D. M. Eissenstat, S. E. Hobbie, J. Oleksyn, A. M. Jagodzinski, P. B.
506
Reich, O. A. Chadwick, and J. Chorover (2012), Tree species effects on coupled
507
cycles of carbon, nitrogen, and acidity in mineral soils at a common garden
508
experiment,
509
doi:10.1007/s10533-011-9695-7.
510 511 512 513
Biogeochemistry,
111(1-3),
601-614,
Nelson, D. W., and L. E. Sommers (1982), Methods of Soil Analysis II, American Society of Agronomy, Madison. Office of Agriculture Resource Layouting of Qinghai (1997), Soil Species of Qinghai, China Agricultural Press, Beijing. 23
514 515
Office of Soil Census of Inner Mongolia (1994), Soil Species of Inner Mongolia Autonomous Region, China Agricultural Press, Beijing.
516
Qin, Y., J. Chen, and S. Yi (2015), Plateau pikas burrowing activity accelerates
517
ecosystem carbon emission from alpine grassland on the Qinghai-Tibetan Plateau,
518
Ecol. Eng., 84, 287-291, doi:10.1016/j.ecoleng.2015.09.012.
519
Quinto-Mosquera, H., and F. Moreno-Hurtado (2016), Precipitation effects on soil
520
characteristics in tropical rain forests of the Chocó biogeographical region,
521
Revista Facultad Nacional de Agronomía, Medellín, 69(1), 7813-7823,
522
doi:10.15446/rfna.v69n1.54749.
523 524
R Development Core Team (2012), R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna.
525
Ruiz Sinoga, J. D., S. Pariente, A. Romero Diaz, and J. F. Martinez Murillo (2012),
526
Variability of relationships between soil organic carbon and some soil properties
527
in Mediterranean rangelands under different climatic conditions (South of Spain),
528
Catena, 94, 17-25, doi:10.1016/j.catena.2011.06.004.
529
Sparks, D. L., A. L. Page, P. A. Helmke, R. H. Loeppert, P. N. Soltanpour, M. A.
530
Tabatabai, C. T. Johnston, and M. E. Sumner (1996), Methods of Soil Analysis.
531
Part III. Chemical Methods.
532 533
Tibet Land Management Bureau (1994), Soil Species of Tibet Autonomous Region, Science Press, Beijing.
534
Tůma, I., K. Fiala, J. Záhora, and P. Holub (2011), The role of Athyrium
535
distentifolium in reduction of soil acidification and base cation losses due to acid 24
536
deposition in a deforested mountain area, Plant Soil, 354(1-2), 107-120,
537
doi:10.1007/s11104-011-1048-8.
538
von Lutzow, M., I. Kogel-Knabner, K. Ekschmitt, E. Matzner, G. Guggenberger, B.
539
Marschner, and H. Flessa (2006), Stabilization of organic matter in temperate
540
soils: mechanisms and their relevance under different soil conditions - a review,
541
Eur. J. Soil Sci., 57(4), 426-445, doi:10.1111/j.1365-2389.2006.00809.x.
542
Watmough, S. A., and P. J. Dillon (2003), Base cation and nitrogen budgets for a
543
mixed hardwood catchment in South-Central Ontario, Ecosystems, 6(7), 675-693,
544
doi:10.1007/s10021-002-0164-y.
545
Were, K., D. T. Bui, O. B. Dick, and B. R. Singh (2015), A comparative assessment of
546
support vector regression, artificial neural networks, and random forests for
547
predicting and mapping soil organic carbon stocks across an Afromontane
548
landscape,
549
doi:10.1016/j.ecolind.2014.12.028.
Ecological
Indicators,
52,
394-403,
550
Wiesmeier, M., F. Barthold, B. Blank, and I. Koegel-Knabner (2011), Digital mapping
551
of soil organic matter stocks using Random Forest modeling in a semi-arid
552
steppe ecosystem, Plant Soil, 340(1-2), 7-24, doi:10.1007/s11104-010-0425-z.
553
Xie, Z., J. Zhu, G. Liu, G. Cadisch, T. Hasegawa, C. Chen, H. Sun, H. Tang, and Q.
554
Zeng (2007), Soil organic carbon stocks in China and changes from 1980s to
555
2000s,
556
doi:10.1111/j.1365-2486.2007.01409.x.
557
Glob.
Change
Biol.,
13(9),
1989-2007,
Xu, L., H. Zhou, C. Liang, L. Du, and H. Li (2010), Spatial and temporal variability 25
558
of annual and seasonal precipitation over the desert region of China during
559
1951-2005, Hydrol. Process., 24(20), 2947-2959, doi:10.1002/hyp.7708.
560
Xu, R., A. Zhao, J. Yuan, and J. Jiang (2012), pH buffering capacity of acid soils from
561
tropical and subtropical regions of China as influenced by incorporation of crop
562
straw
563
doi:10.1007/s11368-012-0483-3.
biochars,
J.
Soils
Sediments,
12(4),
494-502,
564
Xue, X., J. Guo, B. S. A. Han, Q. W. Sun, and L. C. Liu (2009), The effect of climate
565
warming and permafrost thaw on desertification in the Qinghai-Tibetan Plateau,
566
Geomorphology, 108(3-4), 182-190, doi:10.1016/j.geomorph.2009.01.004.
567
Xue, X., Y. Q. Luo, X. H. Zhou, R. Sherry, and X. H. Jia (2011), Climate warming
568
increases soil erosion, carbon and nitrogen loss with biofuel feedstock harvest in
569
tallgrass
570
doi:10.1111/j.1757-1707.2010.01071.x.
prairie,
Global
Change
Biology
Bioenergy,
3(3),
198-207,
571
Yan, H., S. Q. Wang, C. Y. Wang, G. P. Zhang, and N. Patel (2005), Losses of soil
572
organic carbon under wind erosion in China, Glob. Change Biol., 11(5), 828-840,
573
doi:10.1111/j.1365-2486.2005.00950.x.
574
Yang, J., Z. Jiang, P. Wang, and Y. Chen (2008), Temporal and spatial characteristic of
575
extreme precipitation event in China, Climatic and Environmental Research,
576
13(1), 75-83.
577
Yang, X., Z. Ding, X. Fan, Z. Zhou, and N. Ma (2007), Processes and mechanisms of
578
desertification in northern China during the last 30 years, with a special reference
579
to the Hunshandake Sandy Land, eastern Inner Mongolia, Catena, 71(1), 2-12, 26
580
doi:10.1016/j.catena.2006.10.002.
581
Yang, Y., J. Fang, C. Ji, W. Ma, A. Mohammat, S. Wang, S. Wang, A. Datta, D.
582
Robinson, and P. Smith (2012a), Widespread decreases in topsoil inorganic
583
carbon stocks across China's grasslands during 1980s-2000s, Glob. Change Biol.,
584
18(12), 3672-3680, doi:10.1111/gcb.12025.
585
Yang, Y., C. Ji, W. Ma, S. Wang, S. Wang, W. Han, A. Mohammat, D. Robinson, and P.
586
Smith (2012b), Significant soil acidification across northern China's grasslands
587
during
588
doi:10.1111/j.1365-2486.2012.02694.x.
1980s-2000s,
Glob.
Change
Biol.,
18(7),
2292-2300,
589
Yang, Y., P. Li, J. Ding, X. Zhao, W. Ma, C. Ji, and J. Fang (2014), Increased topsoil
590
carbon stock across China's forests, Glob. Change Biol., 20(8), 2687-2696,
591
doi:10.1111/gcb.12536.
592
Yu, C., Y. Zhang, H. Claus, R. Zeng, X. Zhang, and J. Wang (2012), Ecological and
593
environmental issues faced by a developing Tibet, Environ. Sci. Technol., 46(4),
594
1979-1980, doi:10.1021/es2047188.
595
Yu, Y. Y., Z. T. Guo, H. B. Wu, J. A. Kahmann, and F. Oldfield (2009), Spatial
596
changes in soil organic carbon density and storage of cultivated soils in China
597
from 1980 to 2000, Glob. Biogeochem. Cycle, 23, doi:10.1029/2008gb003428.
598
Zhai, P. M., X. B. Zhang, H. Wan, and X. H. Pan (2005), Trends in total precipitation
599
and frequency of daily precipitation extremes over China, J. Clim., 18(7),
600
1096-1108, doi:10.1175/jcli-3318.1.
601
Zhao, Y., L. Duan, J. Xing, T. Larssen, C. P. Nielsen, and J. Hao (2009), Soil 27
602
acidification in China: Is controlling SO2 emissions enough?, Environ. Sci.
603
Technol., 43(21), 8021-8026, doi:10.1021/es901430n.
604
Ziegler, A., and I. R. Konig (2014), Mining data with random forests: current options
605
for real-world applications, Wiley Interdiscip. Rev.-Data Mining Knowl. Discov.,
606
4(1), 55-63, doi:10.1002/widm.1114.
607
28
608
Table 1. Means of CEC and its changes among five grassland types within a grassland
609
transect across northern China from the 1980s-2010s. The values in
610
parenthesis show the 95% confidence interval calculated using the Random
611
Forest Model. Notes: n, sample size; AM, alpine meadow; AS, alpine steppe;
612
MS, meadow steppe; TS, typical steppe; DS, desert steppe. CEC during the 1980s
CEC during the 2010s
Change in CEC
(cmol kg-1)
(cmol kg-1)
(cmol kg-1)
Grassland type
n
Total
280
17.06 [16.04, 18.08] 14.69 [13.76, 15.61]
-2.37 [-2.97, -1.78]
AM
97
22.57 [20.92, 24.23] 20.13 [18.35, 21.91]
-2.45 [-3.71, -1.18]
AS
90
13.84 [12.30, 15.39]
11.07 [9.85, 12.29]
-2.77 [-3.60, -1.94]
MS
16
26.41 [22.14, 30.67] 19.39 [17.24, 21.55]
-7.01 [-10.32, -3.70]
TS
52
13.44 [12.19, 14.69] 12.22 [11.26, 13.18]
-1.22 [-2.17, -0.26]
DS
25
8.83 [7.99, 9.66]
8.72 [7.26, 10.18]
-0.11 [-1.54, 1.32]
613
29
614
Figure Legends
615
Figure 1. The distribution of sampling sites during the 1980s and the 2010s. The
616
background colors show grassland extent and type across northern China
617
[Chinese Academy of Sciences, 2001]. The rectangle represents the
618
grassland transect selected for exploring the temporal dynamics of CEC.
619
Blue points represent the sites during the 1980s and green points represent
620
the 2010s.
621 622 623
Figure 2. Procedure of constructing Random Forest (RF) models to predict data for those sampling sites during the 2000s (a) and during the 1980s (b).
624 625 626
Figure 3. Comparison between predicted and measured CEC for the period of (a) the 1980s and (b) the 2010s.
627 628
Figure 4. Changes in CEC (a) along a grassland transect across northern China’s
629
grasslands and (b) among various grassland types between the 1980s and
630
the 2010s. The inset in panel (a) shows the comparison of CEC between
631
the two periods, with a 1:1 line for reference. Different letters above
632
box-and-whisker plots express significant differences between the two
633
periods. The gap and square, lower and upper edges and bars in the boxes
634
represent median and mean values, 25th and 75th percentiles, and standard
635
deviations of all data, respectively. AM: alpine meadow, AS: alpine steppe, 30
636
MS: meadow steppe, TS: typical steppe, DS: desert steppe.
637 638
Figure 5. Relationships of CEC with (a, e) SOC, (b, f) silt content, (c, g) sand content,
639
and (d, h) MAP. Panels (a-d) and (e-h) illustrate the relationships between
640
CEC and environmental factors during the 1980s and the 2010s,
641
respectively. Data and trends on the Tibetan Plateau are represented by red
642
points and lines, while those on the Inner Mongolian Plateau are shown by
643
blue points and lines.
31
Figure 1.
Figure 2.
Figure 3.
Figure 4.
Figure 5.