Climate factors such as aridity significantly influence soil carbon (C) and
nitrogen (N) stocks in terrestrial ecosystems. Further, soil texture plays
an important role in driving changes of soil C and N contents at regional
scale. However, it remains uncertain whether such changes resulted from
the variation of different soil particle–size factions and/or the C and N
concentrations in those fractions. We examined the distribution of total C
and N in both bulk soil and different soil particle–size fractions,
including sand (53–2000
Grasslands, which cover nearly 40 % of the world's land area, store approximately one-third of the total carbon (C) in terrestrial ecosystems, and more than 70 % of C in grassland ecosystems is stored in the top 1m soil layer (He et al., 2012; White et al., 2000). Consequently, C turnover in grassland soils is considered to be a critical component of the global C cycle (Fisher et al., 1994; Wang et al., 2009). In China, grasslands account for 41.7 % of terrestrial land area and are mainly distributed in arid and semi-arid regions (NSBC, 2002). Carbon and nitrogen (N) cycling in the grasslands of northern China are relatively sensitive to global change factors, such as increases in drought, more extreme precipitation events, and global warming (Wang et al., 2014; Song et al., 2012). However, this sensitivity may show regional variations due to variability in both climate and soil characteristics. Better understanding of the controls on regional variation of soil C and N stocks in Chinese grasslands would facilitate projections of regional C and N cycling under global change scenarios.
Soil C and N stocks in grassland ecosystems are closely correlated with climatic conditions. In arid and semi-arid ecosystems, soil C and N stocks are positively correlated with mean annual precipitation (MAP) at the regional scale (He et al., 2014; Nichols, 1984). This positive relationship is driven by the fact that water availability is the dominant limiting factor for plant growth (and thus soil organic matter inputs) in these ecosystems. In contrast, mean annual temperature (MAT) is negatively correlated with soil C and N stocks, as higher temperature generally enhances microbial decomposition more than detrital production (Homann et al., 2007; Miller et al., 2004; Schimel et al., 1994).
Aridity, which is intensified by decreasing MAP and increasing MAT, is projected to increase in drylands worldwide during this century (Dai, 2013; Delgado-Baquerizo et al., 2013); this change may significantly diminish soil C and N stocks in those regions (Delgado-Baquerizo et al., 2013; Sanaullah et al., 2014). Variation in soil fraction composition and the C and N concentrations in different soil particle sizes may strongly influence the pattern of decreasing bulk soil C and N pools observed with increasing aridity (Amelung et al., 1998; He et al., 2014). However, compared with these well-described patterns of variation in soil C and N stocks along climate and soil particle size gradients, the relative influence of these factors on terrestrial C and N pools under climate change scenarios such as increasing aridity is less clear (Delgado-Baquerizo et al., 2013). Previous studies found that soil C and N tends to decrease with increasing aridity (that is the degree of dryness of the climate at a given location) largely due to decreased primary productivity, as well as other ecosystem processes, such as reduction of total plant cover, shifts of species composition, changes of litter quality, altered litter decomposition rates (Carrera and Bertiller, 2010; Delgado-Baquerizo et al., 2013; Sanaullah et al., 2014).
Variation in soil particle–size fractions also exerts significant controls on the stock and turnover of soil organic matter (SOM; Chen et al., 2010; Christensen, 2001; Qin et al., 2010) and increasing attention has focused upon the responses of C and N pools in different soil particle–size fractions to climate change (Amelung et al., 1998; He et al., 2014). Clay and silt fractions in soil usually have higher C and N concentrations and stocks than that of sand fraction, and thus soils with higher clay and silt contents generally have higher soil organic C (SOC) and N stocks (Amelung et al., 1998; Feller and Beare, 1997; Follett et al., 2012; Hassink, 1997). This pattern reflects that organic materials are preferably decayed from pools of coarse soil particles; these relatively C- and N-rich decomposition products tend to accumulate in finer clay and silt particles (Amelung et al., 1998). Moreover, clay and silt may physically protect organic materials from decomposition and promote the accumulation of recalcitrant material in the fine particle–size fractions of soils (Hassink, 1997; Zhao et al., 2006; Chen and Chiu, 2003).
To increase our understanding of the variation in both the components of different soil particle–size fractions and their C and N concentrations with increasing aridity at the regional scale, we collected soil samples from 58 sites along a 3000 km transect in northern China which covered a wide range of grassland ecosystems locating at the eastern end of the contiguous Eurasian steppe with distinct aridity gradients. The objectives of this study were to (1) examine the distribution of C and N in various particle–size fractions of soils across the aridity gradient, and (2) evaluate the relative contributions of soil fraction components and their element concentrations to the changes of soil C and N stocks across the aridity gradient. We hypothesized the following: (1) that concentrations and stocks of C and N in soil particle–size fractions would be negatively correlated with increasing aridity; and that (2) soil C and N stocks would decline with increasing aridity due to both an increase in the relative proportion of sand to silt and clay and a decrease of C and N concentrations in the three soil fractions.
This study was carried out along a 3000 km west–east transect of arid and
semi-arid grasslands across Xinjiang, Gansu, and Inner Mongolia in northern
China. The distinctive features of this transect include complete
meteorological records, relatively gentle geographical relief, and light
human disturbances (Luo et al., 2015). The longitude of the transect ranged
from 87
Fifty-eight sites were set up along the transect at an interval of 50–100 km. The location and elevation of the sampling sites were measured by GPS
(eTrex Venture, Garmin, USA). For each site, one large plot (50 m
Particle–size fractionation was completed by disrupting soil aggregates of
bulk soil samples using ultrasonic energy and separating the particle–size
fractions by a combination of wet sieving and continuous flow centrifugation
(Chen and Chiu, 2003; He et al., 2009). Briefly, 40 g of sieved soil
(< 2 mm) was dispersed in 200 mL of deionized water (the floating
visible debris was removed) using a probe-type ultrasonic cell disrupter
system (scientz-IID) operating for 15 min in the continuous mode at 361 W.
We used a sieve to separate sand (particle size, 53–2000
Carbon and N stocks (Mg C ha
Similarly, C and N stocks (Mg C ha
All of the relationships between variables were explored by using simple linear regression analyses (58 sites with five subplots as replications in each site). We observed that the relationships were best-fitted by either a first-order equation or a second-order equation. As the contents of sand, silt and clay in soils are not independent of each other, stepwise multiple regression analyses, which are highly conservative (Fornara and Tilman, 2008), were used to determine the simultaneous effects of soil fraction composition and C and N concentrations in soil particle–size fractions on soil C and N stocks. All analyses were performed using SPSS V13.0 (SPSS, Chicago, IL, USA).
The relationships of soil particle–size fractions contents
with aridity differed, with positive correlation between sand content and
aridity, negative correlation between silt content and aridity, and no
significant relationship between clay content and aridity. Data are
presented as mean
Sand was the most abundant fraction for most sites, accounting for
21.62–90.65 % of the total soil weight along the transect. The content of
sand was positively correlated with increasing aridity (Fig. 1). The silt
content, which accounted for 4.19–49.29 % of the total soil weight,
decreased with increasing aridity (Fig. 1). The content of clay was
relatively low across the transect, ranging from 1.36–33.7 %. There were
no significant relationships between clay content and aridity (Fig. 1). Bulk
density ranged from 0.90 to 1.72 g cm
Soil C and N concentrations in bulk soils and different
soil particle–size fractions (sand, silt and clay) at 58 sampling sites in
arid and semi-arid grasslands of northern China (data are represented as
means
Carbon
Total C (Fig. 2a) and N concentrations (Fig. 2b) in the bulk soil
significantly decreased with increasing aridity. Soil C concentration ranged
from 2.71 to 50.33 g C kg
The C
Across the whole transect, C stock in bulk soils (0–10 cm) ranged from 4.36
to 46.16 Mg C ha
Soil C and N stocks in bulk soils and different soil
particle–size fractions (sand, silt, and clay) at 58 sampling sites in arid
and semi-arid grasslands of northern China (data are represented as means
Across the transect, the concentrations and stocks of C and N in bulk soils were negatively correlated with the content of sand and positively correlated with the contents of silt and clay (Fig. 4). The concentrations and stocks of C and N in bulk soils were positively correlated with their concentrations in sand, silt, and clay (Fig. 5).
Across the transect, the concentrations and stocks of C
Stepwise multiple regression analyses allowed us to quantify the
simultaneous effects of soil fraction composition, element concentrations in
soil particle–size fractions, and BD on bulk soil C and N stocks. The
multiple regression model for C stocks in bulk soils included the variables
(Table 3): clay C concentration (with the value of normalized regression
coefficient for this variable
Results of the multiple regressions refer to the final accepted model which just included the effects of the significant variables for C stocks in bulk soils of arid and semi-arid grasslands.
Results of the multiple regressions refer to the final accepted model which just included the effects of the significant variables for N stocks in bulk soils of arid and semi-arid grasslands.
Across this 3000 km aridity gradient, sand, which accounted for 21.62–90.65 % of the total soil weight, was the most abundant fraction; the contents of silt (4.19–49.29 %) and clay (1.36–33.7 %) were much lower, especially in soils from the extremely arid sites. We suspect that this pattern is partly caused by the wind erosion and dust storms, which can be exacerbated by increasing aridity and frequently occur in higher aridity areas of northern China (Wang et al., 2013; Yan et al., 2013; Zhang and Liu, 2010). Wind erosion favors losses of fine soil particles and consequently leads to changes of the soil texture (Feng et al., 2001; Wang et al., 2006; Yan et al., 2013). In arid and semi-arid ecosystems experiencing increasing aridity, soils become more vulnerable to wind erosion because vegetation coverage declines (Zhang and Liu, 2010). Similar to our results, Liu et al. (2008) found that sand fractions in soils of steppe and meadow were negatively correlated with MAP and positively correlated with MAT due to drought-driven vegetation cover decline in the semi-arid East Asian steppe.
The concentrations and stocks of C
Other factors may also contribute to the pattern observed here. Parent materials, land use (e.g. grazing), and topography, can largely influence soil formation process and the contents of soil fractions (Barthold et al., 2013; Deng et al., 2015). For example, soils derived from limestone and quartzite have a lower content of sand fractions and a higher content of silt, compared to soils derived from granite (Belnap et al., 2014). In the grasslands of Inner Mongolia, climate and land use (e.g. intensified grazing leads to soil degradation by diminishing the fine soil fraction) are of greater importance than parent material and topography in controlling soil type distribution (Barthold et al., 2013). Therefore, we suspect that those factors associated with climate would be more important than other factors in structuring the pattern of soil particle–size distribution in our study sites.
Our results showed that C and N concentrations were highest in clay, followed by silt, and much lower in sand across a 3000 km aridity gradient. This pattern, which may be caused because fine fractions in soil have high surface area which can enhance formation of organo-mineral complexes that protect SOM from microbial degradation (Hassink, 1997; Zhang and Liu, 2010), supports previous findings that soil fractionation is a useful tool for examining different C and N pools in soil (Amelung et al., 1998; Gerzabek et al., 2001; Stemmer et al., 1999). Across the transect, we observed that C and N concentrations and stocks in bulk soils were negatively correlated with sand content and positively correlated with silt and clay contents. Similarly, Bai et al. (2007) demonstrated that there was a negative correlation between SOC content and sand content, and there were positive correlations between SOC content and clay and silt contents based in wetland soils in northeastern China. Positive correlations between soil C and N concentrations and silt and clay contents were also found in Inner Mongolian grasslands (He et al., 2014).
Supporting our first hypothesis, we found that C and N concentrations and stocks in soil particle–size fractions tended to be negatively correlated with increasing aridity. The higher aridity sites have lower primary productivity (Wang et al., 2014) and thus a lower input of plant detritus into soil. Lower litter input is correlated with lower C and N concentrations and stocks in soil fractions (Yang et al., 2011; He et al., 2014). Additionally, aridity has been identified to be a major factor affecting bacterial diversity, community composition and taxon abundance in this system (Wang et al., 2015). Therefore, with increasing aridity, microbially mediated litter decomposition may also change due to altered microbial community composition, which may further influence soil C and N (Carrera and Bertiller, 2010). Paralleling our results, He et al. (2014) found that C and N concentrations in soil particle–size fractions were positively correlated with MAP; moreover, they considered that MAP was better than MAT to model the variation of soil C stock in an Inner Mongolia grassland. In the present study, we quantified the relationship between soil C and N stocks and aridity, which combines MAP and MAT, across different sampling sites. Our results suggest that aridity is a robust predictor for the regional variation of C and N stocks in soil fractions.
We found that C stock in sand was first decreased and then increased along the aridity gradient, which seems paradoxical given the results that the C concentrations in sand linearly declined with increasing aridity, while the content of sand linearly increased with increasing aridity across the transect. The observed variation of C stock in sand across the transect may be due to the shifts of dominant controller for the C stock in sand across the aridity transect. Sand C concentration appears to be more important than sand content in driving the variation of sand C stock in the ecosystem with aridity value is less than 0.8 (where the C concentration in sand was relatively higher and sand content was relatively lower). In contrast, sand content appears to be more important than C concentration in determining sand C stock when the aridity value exceeds 0.8. Our results highlight the importance of considering both soil particle size and the C concentration of different particles in order to better understand the influence of aridity on soil C pools.
We found that total C and N concentrations and stocks in bulk soils generally decreased with increasing aridity across the whole transect. Previous studies have reported that soil C and N stocks in the upper soil layers were positively correlated with MAP and negatively correlated with MAT; these findings are similar to our observations along a large aridity gradient (Follett et al., 2012; He et al., 2014; Liu et al., 2012; Miller et al., 2004). The depletion of fine soil particles due to the intensified wind erosion with increasing aridity could further deplete C and nutrients in arid systems because these particles have disproportionately greater amounts of C and nutrients than larger particles (Yan et al., 2013). Furthermore, the decline of plant coverage and aboveground biomass under higher aridity would also contribute to the decreased C and N content along this aridity gradient. Actually, aboveground primary productivity was significantly decreased with increasing aridity along this transect (Wang et al., 2014).
Our results suggest that the decreases of bulk soil C and N stocks along the aridity gradient were resulted not only from the changes of composition of different soil fractions but also from the decreases of C and N concentrations in each of those fractions. While both soil C and N stocks decreased with increasing aridity, the stepwise multiple regression analyses indicated that the simultaneous influences of variation of different soil fractions and the element concentrations were different for C and N. For bulk soil C stock, the most robust regression model did not include sand content, sand C concentration, and BD, whereas for bulk soil N stock, silt content and silt N concentration were excluded from the model. Our results thus demonstrate that sand content is less important than silt content for controlling variation of soil C stock, whereas silt is less important for the variation of soil N stock at regional scale in the arid and semi-arid grasslands of northern China.
These findings are somewhat in agreement with previous findings that C is readily mineralized from un-complexed organic matter in sand-sized aggregates whereas N is not, while silt tends to be more enriched in C than N (Christensen, 2001). We found that clay content and clay element concentrations were the most important factors for predicting the variation of both the soil C and N stocks across this aridity gradient. Similarly, Burke et al. (1989) observed that clay was an important predictor of soil C for American grassland soils. Together, these results indicate differences in the relative importance of different soil particle–size fractions in driving soil C and N stocks, although it is generally accepted that the dynamics of those two elements in soils are closely correlated (Finzi et al., 2011).
This large-scale field investigation provides strong evidence that increasing aridity would reduce the soil C and N stocks in arid and semi-arid ecosystems due both to the changes of particle–sized fractions in soils (i.e. relatively more coarse fraction content, but less fine fraction content with increasing aridity) and to the decline of C and N concentrations in each fraction. This study provides novel insights into the patterns underlying regional changes of soil C and N from a soil particle–size fractions perspective. Given the predicted increases in aridity in this century for the global drylands (Dai, 2013), this study indicates that the soil C and N pools in those arid ecosystems may decline in the future. Because wind erosion would lead to greater loss of relatively fine silt and clay particles (Yan et al., 2013), our results suggest that land use practices which reduce wind erosion (e.g. reducing the intensity of grazing) will play an important role in sustaining soil C sequestration in dryland regions globally.
Along the transect, aridity was an important factor driving the changes of soil C and N concentrations and stocks in the arid and semi-arid grasslands. Both of the C and N concentrations and stocks in the three particle–size fractions as well as in bulk soils tended to be negatively correlated with aridity. The concentrations and stocks of C and N in bulk soils were negatively correlated with sand content but positively correlated with both silt and clay contents, suggesting that fine soil fractions can protect SOM from microbial degradation. There were positive correlations between the concentrations and stocks of C or N in bulk soils and the C or N concentrations in the three soil particle–size fractions. Our results have significant implications for better understanding soil C and N cycles under scenarios of increasing aridity in global drylands that are predicted to occur this century.
We thank all the members in the Shenyang Sampling Campaign Team from the Institute of Applied Ecology, Chinese Academy of Sciences for their assistance in field sampling. We appreciate the constructive comments from two anonymous reviewers and the editor. This work was supported by National Natural Science Foundation of China (31470505), the National Basic Research Program of China (2015CB150802), Strategic Priority Research Program of the Chinese Academy of Sciences (XDB15010403 and XDB15010401), the Key Research Program from CAS (KFZD-SW-305-002) and Youth Innovation Promotion Association CAS (2014174). Edited by: Z. Jia