Open access peer-reviewed chapter

Responses of Community Structure, Productivity and Turnover Traits to Long-Term Grazing Exclusion in a Semiarid Grassland on the Loess Plateau of Northern China

Written By

Jimin Cheng, Wei Li, Jishuai Su, Liang Guo, Jingwei Jin and Chengcheng Gang

Reviewed: 19 February 2019 Published: 23 March 2019

DOI: 10.5772/intechopen.85306

From the Edited Volume

Plant Communities and Their Environment

Edited by Manuel T. Oliveira, Feyza Candan and Anabela Fernandes-Silva

Chapter metrics overview

843 Chapter Downloads

View Full Metrics


Grazing exclusion has been widely used for restoration of degraded grassland all over the world. Based on over a 30-year (from 1982 to 2011) vegetation survey and a 2-year (from 2013 to 2014) field decomposition experiment in Yunwu Mountain Grassland Nature Reserve on the Loess Plateau of China, responses of community structure and productivity and decomposition traits of dominant Stipa species (Stipa bungeana, Stipa grandis and Stipa przewalskyi) litters were determined to reveal the ecosystem cyclic process. Results showed that grassland coverage, plant density, Shannon-Wiener index and aboveground productivity changed in a hump pattern with peaks in 2002. Productivity was significantly positively correlated with mean annual temperature. The direction and magnitude about effects of climatic changes on productivity depended on phonological stages of plant community. Warming in early stage of growing season (April–May) contributed the increase of productivity, while temperature rise after the growing season (September–March in the next following year) was negatively correlated with productivity in the following year. Leaf litters of three Stipa species (S. bungeana, S. grandis and S. przewalskyi) had higher decomposition rates in the growing season than that in the nongrowing season. Nutrient-releasing pattern in litters of three Stipa species followed a different pattern: S. bungeana > S.grandis>S. przewalskyi. Considering productivity and decomposition traits, grazing exclusion promotes carbon sequestration of semiarid grassland, while adjustments in nutrient cycling might explain fluctuations of community structure.


  • yunwu mountain
  • loess plateau
  • grazing exclusion
  • climate variation
  • decomposition
  • Stipa

1. Introduction

As one of the most important and largest terrestrial ecosystems in the world, grasslands cover 30% of the land surface and are mainly distributed in arid and semiarid regions [1]. Due to global climate change and human activity, such as heavy grazing, grasslands in this area have undergone desertification and even virtually disappeared in recent decades [2, 3], making restoration process urgent for degraded grasslands [4, 5]. Current studies about grassland restoration mainly focus on several key components: community composition and structure, species diversity, soil properties and vegetation succession process [6, 7, 8, 9, 10]. Grassland is considered very sensitive to climate changes [11, 12, 13, 14] and also is influenced by soil resource availability [15, 16].

Compared with forest ecosystem and cropland ecosystem, aboveground net primary productivity (ANPP) of grasslands is highly temporally variable [16, 17]. Specifically, climate-driven variability in grassland productivity has important effects on the global carbon balance, ecosystem service delivery, profitability of pastoral livelihoods and the sustainability of grassland resources [11, 18, 19]. Many ecologists have analysed the impacts of annual precipitation and temperature on ANPP at regional and continental scales [17, 20, 21, 22, 23], while numerous site-specific reports have indicated that interannual variability in ANPP is poorly or even not at all correlated with annual climate conditions [19, 24, 25]. Changes in precipitation or temperature during certain parts of the year have been proven to be more relevant drivers of ANPP than annual changes [26, 27, 28, 29], and the impacts on vegetation production varied with seasons [13, 28, 30, 31]. For instance, warming in early spring increased grassland productivity by ameliorating cold temperature constraints on plant growth in northern mid- and high latitudes [32, 33] and advancing spring greening phenology [34, 35, 36]. Temperature increases in summer; however, it can depress productivity by reducing soil moisture and intensifying physiological stress [13].

The Loess Plateau of China has a total area of about 52 million hectares and is widely known for its fragile ecological environment, frequent severe droughts and problems with water runoff and soil erosion [37]. In recent years, the complicated landscape, frequent droughts and severe soil erosion have attracted worldwide attention and caused sustained deterioration of the ecosystem of this region. In contrast to numerous studies in the temperate grasslands of Inner Mongolia and the alpine grasslands of the Tibetan Plateau, very few reports are available on responses of grassland productivity to climate variability on the more arid Loess Plateau in China [3], especially with respect to responses to seasonal climatic variability. Restoration of the natural vegetation is regarded as the most effective method for changing the ecological environment of the Loess Plateau [7, 8, 38].

As a major determinant of nutrient cycling, litter decomposition is a fundamental process of grassland ecosystem functioning [39]. Decomposition traits of plant litters are affected by a number of factors, including litter quality, abiotic environment and soil organisms [40]. In general, plant litters with high C:N ratio and lignin concentration are supposed to have slow decomposition and nutrient immobilisation processes, whereas low C:N ratio and low lignin concentration contribute to fast decomposition and nutrient mineralisation processes. Decomposition traits of plant materials may vary with succession stages. For example, late-seral dominant grasses normally had high tissue N concentrations, low C:N ratios and lignin concentrations, which result into fast decomposition rate and enhanced nutrient mineralisation.

Most previous studies have focused on plant species richness and diversity in abandoned croplands following short-term grazing exclusion in China [8, 41, 42]. Few studies have reported on the restoration succession of typical natural steppe under long-term grazing exclusion [19, 43]. In the present study, Stipa steppe has been fenced from 1982 to the present at Yunwu Mountain National Nature Reserve, and long-term grassland ecological characteristics, productivity and weather records have been collected.

The community in the study area consists of 313 plant species, covering 56 families and 165 genera, with five main families being Compositae, Gramineae, Leguminosae, Rosaceae and Labiatae [44]. The dominant Stipa plants include S. bungeana, S. grandis and S. przewalskyi, and main forbs are Thymus mongolicus, Artemisia sacrorum and Potentilla acaulis [45]. Genus-specific morphological and functional traits contributed the dominance of Stipa plants in temperate, subtropical and tropical steppe in semiarid areas worldwide [46]. Meanwhile, Stipa species showed various adaptabilities to environmental changes, presenting an ecological distribution pattern along the climate gradients [47]. There are 32 species, 1 subspecies and 3 variations in genus Stipa plants in China, mainly distributed in western and northeastern area, and 5 Stipa species are found in our study area. As the constructive species, S. bungeana mainly distributed on the Loess Plateau [48]. Noticeably, replacement of dominant Stipa species occurred during the long-term restoration process, with Stipa bungeana being replaced by S. grandis and S. przewalskyi [49]. The three Stipa species differentiated in their phenotypic traits. In detail, S. grandis owns higher plant height, and S. przewalskyi possesses more tillers. Besides, S. bungeana and S. przewalskyi consistently flower and produce seeds earlier than S. grandis [49].

The temperature and precipitation variability during 1982–2011 were assessed in this study; the ecological characteristics during long-term grazing exclusion were examined; the relationship between grassland productivity and variation in climate variables were explored; and the variations in decomposition traits of three Stipa dominant species (S. bungeana, S. grandis and S. przewalskyi) were determined.


2. Material and method

2.1 Study site

This study was conducted in Yunwu Mountain National Nature Reserve on the Loess Plateau (106°24′–106°28′ E, 36°13′–36°19′ N) (Figure 1) [45, 50]. Grassland in this area was restored from grazing as a long-term ecological monitoring station since 1982. The elevation of this study area is 1800–2180 m and has a total area of 6660 hm2. The mean annual temperature is 7.01°C, and there are on average 137 frost-free days per year [49]. The mean annual precipitation is 425 mm, with 60–75% of rainfall falling during July–September. The mean annual evaporation is 1017–1739 mm. Snow cover depth in winters averaged 1.2 cm during the dormancy period. The vegetation type is typical steppe. Gentianaceae, Stipa and Potentilla are important plant components, and the main dominant species include S. bungeana, Stipa grandis, S. przewalskyi, Thymus mongolicus, Artemisia sacrorum, Potentilla acaulis and Androsace erecta [45]. Soil type is montane grey-cinnamon soil [45].

Figure 1.

Location of experimental site.

2.2 Experimental design and sampling

2.2.1 Grassland ecological survey

The grassland sites have been restored from grazing exclusion since 1982, and consequently goat grazing was excluded [45, 49, 50]. Three equal-sized transect of 300 × 100 m was established at the top, middle and down positions of the same slope, respectively. And, 15 quadrats (1 × 1 m) were established within each transect. The vegetation survey was carried out in mid- or late August each year during 1982–2011. Plant coverage, height, species abundance and plant density in each quadrat were measured. Aboveground parts of grassland plants were clipped and dried at 65°C for 48 h to determine aboveground biomass [43]. Plant roots of 0–120 cm soil layers were collected with a soil auger of 9 cm diameter, then were washed and dried to determine belowground biomass.

Important value (IV) was used to describe the importance of species in grassland community during the restoration process. Shannon-Wiener index was used to indicate diversity and evenness of plant community [50]. All indices were calculated according to 8 and 43.

Important value (IV)


where IV is the important value, RH is the relative height, RC is the relative coverage, RA is the relative abundance and RF is the relative frequency.

Diversity index (H), using Shannon-Wiener index


where S is the total species number of a quadrat and Pi is the relative importance value of species i.

2.2.2 Litter decomposition experiment

Considering the difficulty of gathering sufficient senesced leaves, leaves of three Stipa species (S. bungeana, S. grandis and S. przewalskyi) were collected in August of 2013 and then dried at 40°C as decomposition materials, according to other decomposition studies [51, 52, 53]. Leaf litters were cut into pieces of 10 cm in length and enclosed in nylon bag (15 g bag−1, 15 × 10 cm, 0.15 mm mesh).

In early October of 2013, the leaf litterbags of three Stipa species were transferred to grassland site restored for 23 years. Four plots of 10 × 10 m were established, and seven leaf litterbags of each Stipa species were placed on the soil surface and secured in place with iron nails on each of four plots. Four leaf litterbags of each Stipa species were harvested after 1, 3, 6, 9, 12, 18 and 24 months of incubation.

In the laboratory, leaf litters were removed from bags, cleaned to remove any extraneous material and weighed after drying at 65°C for 48 h. Leaf litters were analysed for carbon (C), nitrogen (N) and phosphorus (P). C was determined by oxidation with potassium dichromate in a heated oil bath. N was determined by the semimicro Kjeldahl method. P was determined by Olsen method [54].

According to [55], decomposition rate (k) of leaf litters was estimated by the negative exponential decay function:


where X is the remaining mass, X0 is the initial mass and t is the decaying time (year).

Based on the nutrient concentration and remaining mass, we further calculated nutrient accumulation index (NAI) for C, N and P of leaf litters during decomposition process [56, 57]:


where X0 and C0 indicate initial leaf litter mass and chemical element concentration, respectively. Xt and Ct indicate remaining leaf litter mass and chemical element concentration after a period of time t (year), respectively.

2.3 Data analyses

All data in the paper are presented as mean ± standard error. A two-way analysis of variance was conducted to determine the effects of decomposition time, species and their interaction on decomposition rate, nutrient concentration and NAI of leaf litters. A linear mixed model was used to examine correlations of vegetative indices with restoration time, productivity with climate variables and remaining mass with decomposition time. Significant differences of all statistical tests were estimated at a significance level of P < 0.05. All statistical analyses were performed using SPSS 18.0 (SPSS Inc., Chicago, IL, USA).

Partial least squares (PLS) regression was used to analyse the responses of grassland productivity to variation in daily temperature and precipitation during all 365 days of the year based on data for 1992–2011 [58, 59]. The two major outputs of PLS analysis are the variable importance in the projection (VIP) and standardised model coefficients. The VIP threshold for considering variables as important is often set to 0.8 [60]. The standardised model coefficients indicate the strength and direction of the impacts of each variable in the PLS model. The root-mean-square errors (RMSE) of the regression analyses were calculated to determine the accuracy of the PLS model. In the PLS analyses, periods with VIP greater than 0.8 and high absolute values of model coefficients represent the relevant phases influencing grassland productivity. Positive model coefficients indicate that increasing temperature or precipitation during the respective period should increase ANPP, while negative model coefficients imply negative impacts on productivity.


3. Results

3.1 Temperature and precipitation changes

The annual mean air temperature had an increasing trend and increased by 1.17°C from 1982 to 2011(Figure 2A). In contrast with mean annual temperature, mean annual precipitation showed a decreasing trend and larger intra- and interannual variations in our study, indicating the warmer and drier climate. The mean annual precipitation from 1982 to 2011 was 425.42 mm, with markedly lower values in 1986, 1991 and 1999 and with higher values in 2003 (Figure 2B).

Figure 2.

Mean annual air temperature (a) and mean annual rainfall (B) of growing season, nongrowing season and entire year at Yunwu Mountain during 1982–2011.

3.2 Vegetative ecological characteristics after grazing exclusion

Grassland coverage, plant species richness (number of plant species), plant density (number of plant individuals) and Shannon-Wiener index had similar variation tendencies during the three-decade restoration process (Figure 3). Initially, the coverage, plant richness, plant density and Shannon-Wiener index significantly increased. After 20 years’ restoration, they reached peak values of 92.47%, 29.33 species m−2, 161.8 individuals m−2 and 2.93, respectively. With grazing exclusion process continuing, the four indices’ values decreased to 88.73%, 28.2 species m−2, 138.7 individuals m−2 and 2.47, respectively (Figure 3).

Figure 3.

Changes of coverage (a), number of plants (B), plant density (C) and Shannon-wiener index (D) of grassland with grazing exclusion time at 5, 10, 15, 20, 25 and 30 years.

3.3 Biomass changes in grassland community after grazing exclusion

There were significant differences in aboveground biomass between four groups and between total aboveground biomass and total belowground biomass (Figure 4). Aboveground biomasses of four plant groups increased with restoration time after grazing exclusion. Aboveground biomass of Gramineae and Compositae peaked at the 20th year, while that of Leguminosae peaked at the 25th year, and that of Weeds families peaked at the 15th year during restoration process after grazing exclusion. Considering the reduced biomasses of weed families, long-term grazing exclusion improved forage quality of grassland. Meanwhile, aboveground and belowground community biomasses were both increased by grazing exclusion. Since grassland mainly consisted of plants belonging to Gramineae and Compositae, peaks of the total above- and belowground community biomass both occurred at the 20th year, with aboveground community biomass of 520.5 g m−2 and belowground community biomass of 3240.2 g m−2 (Figure 4).

Figure 4.

Biomass changes of Gramineae (a), Leguminosae (B), Compositae (C), weeds (D), aboveground community (E) and belowground community (F) of grassland with grazing exclusion time at 5, 10, 15, 20, 25 and 30 years.

3.4 Responses of aboveground productivity to climate variation

Regression analysis showed that ANPP was significantly correlated with MAT (Figure 5b) but was little influenced by AP variations (Figure 5a).

Figure 5.

Correlations between ANPP and annual precipitation (a) and mean annual temperature (b) during 1992–2011 at Yunwushan. AP means annual precipitation and MAT represents mean annual temperature.

The VIP and standardised model coefficients of the PLS analysis showed that impacts of warming on grassland productivity varied with season periods (Figure 6a). Different with the clear-cut impacts of temperature on ANPP, precipitation showed more complex impacts (Figure 6b).

Figure 6.

Results of partial least squares (PLS) regression correlating grassland productivity at Yunwu Mountain during 1992–2011 with 15-day running means of (a) daily mean temperature and (b) daily precipitation previously from September to august. Blue bars in the top row indicate that VIP values are greater than 0.8, the threshold for variable importance. In the middle row, red colour means model coefficients are negative and important, while green colour indicates important positive relationships between grassland productivity and climate variables. The black lines in the bottom panel stand for daily mean temperature and precipitation, while grey, green and red areas represent the standard deviation of daily climate variables.

3.5 Decomposition traits of leaf litters of three dominant Stipa species

The remaining mass of leaf litters decreased with decomposition time and showed significant differences among three Stipa species (Figure 7). At the end of decomposition experiment, the remaining masses of leaf litters of S. bungeana, S. grandis and S. przewalskyi were 64.47%, 61.53% and 65.78%, respectively (Table 1).

Figure 7.

The remaining mass dynamics of leaf litters of three Stipa species during 2 years’ field decomposition process.

SpeciesRemaining massk-Value
First yearSecond yearFirst yearSecond year
S. bungeana70.05 ± 3.91 b64.47 ± 3.66 ab0.3600.236
S. grandis73.97 ± 1.81 ab61.53 ± 5.24 b0.3200.242
S. przewalskyi79.18 ± 1.49 a65.77 ± 1.80 a0.2370.225

Table 1.

Comparisons of litter decomposition traits after 1 and 2 years’ decomposition between three Stipa species.

Different lowercase letters in the same column indicate significant differences (P < 0.05).

During 2 years’ decomposition process, variations of nutrient concentration were affected by the nutrient type (Figure 8). In detail, concentrations of carbon and nitrogen showed species-specific fluctuations with decreasing tendency among three Stipa species. In contrast, phosphorus concentrations in leaf litters were averaged doubled. There were significant differences in C:N ratio and nutrient accumulation index (NAI) of leaf litters among three Stipa species (Table 2).

Figure 8.

Dynamic of carbon (a), nitrogen (b), phosphorus (c), concentrations and C:N ratio (d) of leaf litters of three Stipa species during 2 years’ field decomposition process.

VariablesdfConcentration (g·kg−1)C/NNAI
Time60.575 ns4.701 **39.564 ***3.877**49.738 ***23.944 ***53.070 ***
Species20.613 ns18.860 ***2.991 ns9.074**0.560 ns11.026 ***50.008 ***
121.163 ns1.843 ns1.224 ns1.889 ns1.663 ns1.014 ns1.185 ns

Table 2.

Analysis of variance of decomposition time, species for nutrient concentration, C:N ratio and NAI.

NAI, nutrient accumulation index; ns indicates no significant effects (P > 0.05). ** and *** indicate significant effects at P < 0.01 and P < 0.001 level, respectively.

Different with nutrient concentrations, nutrient accumulation indices in Figure 9 indicated that C, N and P were all mineralised into soils during the decomposition process. There was no significant difference between species for carbon-releasing pattern (Figure 9).

Figure 9.

NAI dynamics for carbon (a), nitrogen (b), phosphorus (c) of leaf litters of three Stipa species during 2 years’ field decomposition process.


4. Discussion

Anthropogenic activities and climate changes have made important impacts on terrestrial ecosystem structure and functions in the last century [30]. Global surface temperatures during the twentieth century was increased by 0.56–0.92°C, while temperatures are predicted to have an increment of 2.0–4.5°C in the twenty-first century [61]. Annual mean air temperature was increased by 1.17°C from 1982 to 2011 in this study, having similar temperature changing trends with study in Xilingol steppe of Inner Mongolia [61]. In detail, temperature rises differentiated with seasons, with temperature rises of 1.01°C and 1.68°C in growing season and nongrowing season, respectively. Thus, the nongrowing season experienced a higher temperature rise than the growing season. In contrast with mean annual temperature, mean annual precipitation showed a decreasing trend and larger intra- and interannual variations in our study, indicating the warmer and drier climate. Previous researches have shown that vegetation characteristics could be improved using grazing exclusion in the degraded sandy grasslands, alpine meadow and wetlands in China [5, 62]. However, many of these restoration studies were based on a relatively short-term scale and the research strategy focusing on the spatial series substitute for temporal series methods [5, 63]. In this study, community coverage, plant species richness, plant density and Shannon-Wiener index had similar variation tendencies during the three-decade restoration process. After 20 years’ restoration, they reached peak values, but these four index values decreased in the following years. These decreases mainly resulted from accumulation of litter, which reduced the access to light for plant seedlings [64, 65]. Overall, 30 years’ restoration made plant species richness increase from 9.5 species m−2 to 28 species m−2 and make grassland coverage increase from 25 to 85%. In addition, plants were categorised into four groups: Gramineae, Leguminosae, Compositae and weeds. Considering the reduced biomasses of weeds, long-term grazing exclusion improved forage quality of grassland. Meanwhile, aboveground and belowground community biomasses were both increased by grazing exclusion.

The rapid recovery due to grazing exclusion played a more important role than climatic variations in regulating grassland ecosystem. Therefore, datasets of aboveground grassland biomass and climate variables during 1992–2011 were used to examine the impacts of climate variations on aboveground net primary productivity (ANPP). Regression analysis showed that ANPP was significantly correlated with MAT and was little influenced by AP variations, while precipitation is regarded as the most important determinant of grassland productivity in arid and semiarid regions [19, 21, 66]. Considering the neglected temporal variation of annual climate variables, more attentions should be paid to studies at higher temporal resolution attributing impacts of climate variation on grassland productivity to seasonal or even daily variation in climatic variables rather than to annual variation [26, 27, 28, 29, 31]. A low root-mean-square error (RMSE) of 8.13 g m−2 indicated a good fit of the data for the resulting PLS model. The VIP and standardised model coefficients of the PLS analysis showed that impacts of warming on grassland productivity varied with season periods. Since model coefficients in April and May were always positive and VIP values mostly exceeded 0.8, warming in this period had a positive impact on grassland productivity. The positive impacts of warming in spring on grassland productivity may result from increased water absorption, N mineralisation, accelerated snowmelt and advanced spring greening for plants, which may lengthen the growing season and increase photosynthesis and carbon acquisition for plants [13, 67, 68, 69].

Warming in summer (June–July) depressed productivity, forming a striking contrast with the impacts of spring warming. The results can be explained by physiological stress for plant growth generated by warming in summer coinciding with drought [70]. Moreover, warming in summer may reduce soil moisture by increasing evapotranspiration [71]. It is believed that climate variations make impacts on grassland productivity through changes of soil moisture [24, 72, 73]. Furthermore, continuous warming and drought in summer reduced productivity by limiting soil resource availability [74, 75]. And, temperature variation in August had no apparent impacts on grassland productivity.

The majority of published studies have focused on productivity responses to climate variability during the growing season. However, the importance of winter climate is getting more and more attentions [76, 77, 78, 79, 80]. Considering the majority of model coefficients during September–March, high temperature at that time was unfavourable for productivity of the following year. Temperature increases during September–October delay the senescence of grassland, which may increase soil nutrient and water depletion, inhibiting biomass production in the following year [36, 69, 81]. Our results were similar with warming experiments in two limestone grasslands in the UK, which showed that winter heating combined with drought reduced the biomass of both communities [11]. Besides, warmer winter can accelerate snowmelt, resulting in declines of snow cover accompanied with increases of frequency of freezing events, which exerted negative impacts on plant growth [76, 82]. Also, warming in winter may delay the fulfilment of chilling requirements of plants for resuming growth in the following spring or even delay onset of spring phenology [58, 59, 77, 78, 79, 83].

Interestingly, some short intervals with positive coefficients during 1 November–29 March were detected during 1992–2011, indicating a complex physiological and ecological process in dormancy period of grassland. Taking a broader view at model coefficients and aiming at consistency with established phenological phases, we interpreted the entire period (November–March) as another relevant period during which temperature increases appeared to reduce grassland productivity. Therefore, we recommend that more scientific attention should be paid to impacts of winter warming on grassland productivity and the timing of spring phenology events.

The daily precipitation values between the previous September and August were also used as independent variables in the PLS analysis. The resulting model still proved to be a good fit for the data, with an RMSE of 6.53 g m−2. Different with the clear-cut impacts of temperature on ANPP, precipitation showed more complex impacts. Precipitation increases in June and July had positive impacts on productivity, while increasing precipitation during the senescence period (September–October) and the early growing season (April–May) was correlated with low productivity. In contrast to studies reporting the positive impacts of precipitation during April–May on grassland productivity [29, 30], results in the present study can be explained by the site hydrology, with frequent winter snow providing sufficient soil water for plant growth, making sporadic precipitation during April–May (with an average of 59.5 mm during 1992–2011) which has less important direct impacts on grassland productivity. Similarly, there was also no significant relationship between grassland ANPP and precipitation in August. Similar results have also been reported for grasslands in Kansas, USA [13]. During the dormancy period, positive impacts of precipitation were almost offset by negative ones; thus, precipitation seemed to have little impacts on grassland productivity during this period.

Investigating the decomposition traits of dominant Stipa species’ (S. bungeana, S. grandis and S. przewalskyi) litters can reveal the ecosystem cyclic process under grazing exclusion and climatic changes. The remaining mass of leaf litters decreased with decomposition time and showed significant differences among three Stipa species. At the end of decomposition experiment, the remaining mass of leaf litters of S. bungeana, S. grandis and S. przewalskyi were 64.47%, 61.53% and 65.78%, respectively. Therefore, S. grandis decomposed fast, and S. przewalskyi had a slow decomposition rate. Additionally, leaf litters decomposed faster in growing season (6–12 month and 18–24 month) than in nongrowing season (0–6 month and 12–18 month). The decomposition rate (k) was calculated based on the regression of negative exponential decay function, with k-values of 0.360, 0.320 and 0.237 after 1 year’s decomposition for S. bungeana, S. grandis and S. przewalskyi, respectively. Similarly, k-values after 2 years’ decomposition of S. bungeana, S. grandis and S. przewalskyi were 0.236, 0.242 and 0.225, respectively. Since higher k-values indicate higher decomposition rates, we concluded that litter’s decaying progress became difficult as decomposition time increases, mainly due to the depletion of soluble compounds and easily decayed parts at the beginning of decomposition process, leaving hard parts such as lignin to decay slowly [39].

The variations of nutrient concentration were affected by nutrient type during 2 years’ decomposition process during 2013–2014(Figure 8 and Table 2). In detail, concentrations of carbon and nitrogen showed species-specific fluctuations with decreasing tendency among three Stipa species (Figure 8a and b). In contrast, phosphorus concentrations in leaf litters were averaged doubled (Figure 8c), indicating immobilisation of P in the leaf litters, possibly due to microbial immobilisation through the uptake of P from soil solution and translocation of P from fungal hyphae [84]. There were significant differences in C:N ratio of leaf litters among three Stipa species (Figure 8d). S. przewalskyi had higher C:N ratio than S. bungeana, which explained the differences of decomposition rates between them. As the dominant species in late succession stage of grassland, C:N ratio of S. przewalskyi litters did not show a lower value as predicted from other studies [51], possibly due to the divergences of climate and species between two regions. C:N ratio has been proven to be negatively correlated with decomposition rate. Besides, the lower k-value after 2 years’ decomposition process could be explained by the increased C:N ratio of leaf litters. Compared with nutrient concentrations, nutrient accumulation indices indicated that C, N and P were all mineralised into soils during the decomposition process. There was no significant difference between species for carbon-releasing pattern. Still, NAI value for C of S. przewalskyi was higher than two other Stipa species after 2 years’ decomposition (Figure 9). The lower NAI values for N and P of S. bungeana indicated that S. bungeana released more N and P to soil than the two other Stipa species. From this perspective, replacement of Stipa species after long-term grazing exclusion might inhibit nutrient cycling of grassland ecosystem, due to the lower nutrient mineralisation in leaf litters of two Stipa species at middle and late succession stage.


5. Conclusion

The present study indicated that grazing exclusion induced positive effects on grassland vegetative characteristics, with peak values in the 20th year (2002), and long-term grazing exclusion led to decreased species diversity and biomass and can inhibit grassland renewing due to the litter accumulation. Besides, nutrient cycling in grassland might be slowed down through replacement of dominant species during long-term grazing exclusion. Grassland productivity was more influenced by temperature than precipitation. Results indicated that analysis of productivity responses should account not only for the magnitude of climate variation but also for its timing. Climate warming might prolong/shorten growing season by advancing/delaying onset of greenness of plant community. Warmer winter further decreases ANPP, and impacts of warming in early spring should also be considered in evaluating ANPP variability. Therefore, more scientific attention should be paid to trends in spring phenology and their impacts on productivity at species and community levels.



This work was supported by the National Key Research and Development Program of China (2016YFC0500700), National Natural Science Foundation of China (41601586, 41230852, 41701606, 31602004, 31601987, 41601053), Deployment Program of the Chinese Academy of Sciences (KJZD-EW-TZ-G10), Technical System of Modern Forage Industry of the Ministry of Agriculture (CARS-35-40) and Key Cultivation Project of the Chinese Academy of Sciences—the promotion and management of ecosystem functions of restored vegetation on the Loess Plateau, China. We thank the staff at the management bureau of Yunwu Mountain National Nature Reserve for fieldworks and weather data collection since the early 1980s.


  1. 1. Shi Y, Wang Y, Ma Y, Ma W, Liang C, Flynn DFB, et al. Field-based observations of regional-scale, temporal variation in net primary production in Tibetan alpine grasslands. Biogeosciences. 2014;11(7):2003-2016
  2. 2. Zhang JT, Ru WM, Li B. Relationships between vegetation and climate on the loess plateau in China. Folia Geobotanica. 2006;41:151-163
  3. 3. Slimani H, Aidoud A, Roze F. 30 years of protection and monitoring of a steppic rangeland undergoing desertification. Journal of Arid Environments. 2010;74:685-691
  4. 4. Simmons MT, Venhaus HC, Windhager S. Exploiting the attributes of regional ecosystems for landscape design: The role of ecological restoration in ecological engineering. Ecological Engineering. 2007;30(3):201-205
  5. 5. Wang XH, Yu JB, Zhou D, Dong HF, Li YZ, Lin QX, et al. Vegetative ecological characteristics of restored reed (Phragmites australis) wetlands in the Yellow River Delta, China. Environmental Management. 2012;49:325-333
  6. 6. Walker K, Stevens P, Stevens D, Mountford J, Manchester S, Pywell R. The restoration and re-creation of species-rich lowland grassland on land formerly managed for intensive agriculture in the UK. Biological Conservation. 2004;119(1):1-18
  7. 7. Zhang JT. Succession analysis of plant communities in abandoned croplands in the eastern loess plateau of China. Journal of Arid Environments. 2005;63(2):458-474
  8. 8. Zhang JT, Dong YR. Factors affecting species diversity of plant communities and the restoration process in the loess area of China. Ecological Engineering. 2010;36(3):345-350
  9. 9. Wei J, Cheng J, Li W, Liu W. Comparing the effect of naturally restored forest and grassland on carbon sequestration and its vertical distribution in the Chinese loess plateau. PLoS One. 2012;7(7):e40123
  10. 10. Li W, Epstein HE, Wen ZM, Zhao J, Jin JW, Jing GH, et al. Community-weighted mean traits but not functional diversity determine the changes in soil properties during wetland drying on the Tibetan plateau. Solid Earth. 2017;8:137-147
  11. 11. Grime JP, Brown VK, Thompson K, Masters GJ, Hillier SH, Clarke IP, et al. The response of two contrasting limestone grasslands to simulated climate change. Science. 2000;289(5480):762-765
  12. 12. Knapp AK, Fay PA, Blair JM, Collins SL, Smith MD, Carlisle JD, et al. Rainfall variability, carbon cycling, and plant species diversity in a Mesic grassland. Science. 2002;298(5601):2202-2205
  13. 13. Craine JM, Nippert JB, Elmore AJ, Skibbe AM, Hutchinson SL, Brunsell NA. Timing of climate variability and grassland productivity. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(9):3401-3405
  14. 14. Hsu JS, Powell J, Adler PB. Sensitivity of mean annual primary production to precipitation. Global Change Biology. 2012;18(7):2246-2255
  15. 15. Hastings A, Byers JE, Crooks JA, Cuddington K, Jones CG, Lambrinos JG, et al. Ecosystem engineering in space and time. Ecology Letters. 2007;10(2):153-164
  16. 16. Fang J, Piao S, Tang Z, Peng C, Ji W. Interannual variability in net primary production and precipitation. Science. 2001;293(5536):1723
  17. 17. Knapp AK, Smith MD. Variation among biomes in temporal dynamics of aboveground primary production. Science. 2001;291(5503):481-484
  18. 18. Guo Q , Hu Z, Li S, Li X, Sun X, Yu G. Spatial variations in aboveground net primary productivity along a climate gradient in Eurasian temperate grassland: Effects of mean annual precipitation and its seasonal distribution. Global Change Biology. 2012;18(12):3624-3631
  19. 19. Sala OE, Gherardi LA, Reichmann L, Jobbágy E, Peters D. Legacies of precipitation fluctuations on primary production: Theory and data synthesis. Philosophical Transactions of the Royal Society B. 2012;367(1606):3135-3144
  20. 20. Rosenzweig ML. Net primary productivity of terrestrial communities: Prediction from climatological data. The American Naturalist. 1968;102(923):67-74
  21. 21. Lauenroth WK, Sala OE. Long-term forage production of north American shortgrass steppe. Ecological Applications. 1992;2(4):397-403
  22. 22. Bai Y, Han X, Wu J, Chen Z, Li L. Ecosystem stability and compensatory effects in the Inner Mongolia grassland. Nature. 2004;431(7005):181-184
  23. 23. Hu Z, Fan J, Zhong H, Yu G. Spatiotemporal dynamics of aboveground primary productivity along a precipitation gradient in Chinese temperate grassland. Science in China Series D: Earth Sciences. 2007;50(5):754-764
  24. 24. Hsu JS, Adler PB. Anticipating changes in variability of grassland production due to increases in inter-annual precipitation variability. Ecosphere. 2014;5(5):1-15
  25. 25. Guo L, Cheng JM, Luedeling E, Koerner SE, He JS, Xu JC, et al. Critical climate periods for grassland productivity on China’s loess plateau. Agricultural and Forest Meteorology. 2017;233:101-109
  26. 26. Milchunas DG, Forwood JR, Lauenroth WK. Productivity of long-term grazing treatments in response to seasonal precipitation. Journal of Range Management. 1994;47(2):133-139
  27. 27. Chou WW, Silver WL, Jackson RD, Thompson AW, Allen-Diaz B. The sensitivity of annual grassland carbon cycling to the quantity and timing of rainfall. Global Change Biology. 2008;14(6):1382-1394
  28. 28. La Pierre KJ, Yuan S, Chang CC, Avolio ML, Hallett LM, Schreck T, et al. Explaining temporal variation in above-ground productivity in a Mesic grassland: The role of climate and flowering. Journal of Ecology. 2011;99(5):1250-1262
  29. 29. Robinson TMP, La Pierre KJ, Vadeboncoeur MA, Byrne KM, Thomey ML, Colby SE. Seasonal, not annual precipitation drives community productivity across ecosystems. Oikos. 2013;122(5):727-738
  30. 30. Ma WH, Liu ZL, Wang ZH, Wang W, Liang CZ, Tang YH, et al. Climate change alters interannual variation of grassland aboveground productivity: Evidence from a 22-year measurement series in the inner Mongolian grassland. Journal of Plant Research. 2010;123:509-517
  31. 31. Hovenden MJ, Newton PCD, Wills KE. Seasonal not annual rainfall determines grassland biomass response to carbon dioxide. Nature. 2014;511(7511):583-586
  32. 32. Chen F, Weber KT. Assessing the impact of seasonal precipitation and temperature on vegetation in a grass-dominated rangeland. Rangeland Journal. 2014;36(2):185-190
  33. 33. Chollet S, Rambal S, Fayolle A, Hubert D, Foulquié D, Garnier E. Combined effects of climate, resource availability, and plant traits on biomass produced in a Mediterranean rangeland. Ecology. 2014;95(3):737-748
  34. 34. Menzel A, Estrella N, Fabian P. Spatial and temporal variability of the phenological seasons in Germany from 1951 to 1996. Global Change Biology. 2001;7(6):657-666
  35. 35. Piao S, Mohammat A, Fang J, Cai Q , Feng J. NDVI-based increase in growth of temperate grasslands and its responses to climate changes in China. Global Environmental Change. 2006;16(4):340-348
  36. 36. Reyes-Fox M, Steltzer H, Trlica MJ, Mcmaster GS, Andales AA, Lecain DR, et al. Elevated CO2 further lengthens growing season under warming conditions. Nature. 2014;510(7504):259-262
  37. 37. Xin Z, Xu J, Zheng W. Spatiotemporal variations of vegetation cover on the Chinese loess plateau (1981-2006): Impacts of climate changes and human activities. Science in China Series D: Earth Sciences. 2008;51(1):67-78
  38. 38. Dang XH, Liu GB, Xue S. Models of soil and water conservation and ecological restoration in the loess hilly region of China. Transactions of the Chinese Society of Agricultural Engineering. 2010;26(9):72-80
  39. 39. Moretto AS, Distel RA. Decomposition of and nutrient dynamics in leaf litter and roots of Poa ligularis and Stipa gyneriodes. Journal of Arid Environments. 2003;55(3):503-514
  40. 40. Callaway RM. Positive interactions among plants. The Botanical Review. 1995;61:306-349
  41. 41. Jiao JY, Tzanopoulos J, Xofis P, Bai WJ, Ma XH, Mitchley J. Can the study of natural vegetation succession assist in the control of soil erosion on abandoned croplands on the loess plateau, China? Restoration Ecology. 2007;15(3):391-399
  42. 42. Jiao JY, Tzanopoulos J, Xofis P, Mitchley J. Factors affecting the distribution of vegetation types on abandoned cropland in the hilly-gullied loess plateau region of China. Pedosphere. 2008;18(1):24-33
  43. 43. Cheng JM, Cheng J, Shao HB, Yang XM. Effects of mowing disturbance on the community succession of typical steppe in the loess plateau of China. African Journal of Agricultural Research. 2012;7(37):5224-5232
  44. 44. Zhu RB, Cheng JM, Liu YJ, Li WJ, Wei L. Floristic study of spermatophyte in Yunwu Mountain natural reserves of China. Acta Agrestia Sinica. 2012;20(3):439-443
  45. 45. Qiu LP, Wei XR, Zhang XC, Cheng JM. Ecosystem carbon and nitrogen accumulation after grazing exclusion in semiarid grassland. PLoS One. 2013;8(1):e55433
  46. 46. Armas C, Kikvidze Z, Pugnaire FI. Abiotic conditions, neighbour interactions, and the distribution of Stipa tenacissima in a semiarid mountain range. Journal of Arid Environments. 2009;73:1084-1089
  47. 47. Lv XM, Zhou GS, Wang YH, Song XL. Sensitive indicators of zonal Stipa species to changing temperature and precipitation in Inner Mongolia grassland, China. Frontiers in Plant Science. 2016;7:73
  48. 48. Lu SL, Wu ZL. On geographical distribution of the genus Stipa L. in China. Acta Phytotaxonomica Sinica. 1996;34(3):242-253
  49. 49. Su JS, Jing GH, Jin JW, Wei L, Liu J, Cheng JM. Identifying drivers of root community compositional changes in semiarid grassland on the loess plateau after long–term grazing exclusion. Ecological Engineering. 2017;99:13-21
  50. 50. Jing ZB, Cheng JM, Su JS, Bai Y, Jin JW. Changes in plant community composition and soil properties under 3-decade grazing exclusion in semiarid grassland. Ecological Engineering. 2014;64:171-178
  51. 51. Moretto AS, Distel RA, Didoné NG. Decomposition and nutrient dynamic of leaf litter and roots from palatable and unpalatable grasses in a semi–arid grassland. Applied Soil Ecology. 2001;18(1):31-37
  52. 52. Mathers NJ, Jalota RK, Dalal RC, Boyd SE. 13C–NMR analysis of decomposing litter and fine roots in the semi–arid Mulga lands of southern Queensland. Soil Biology and Biochemistry. 2007;39(5):993-1006
  53. 53. McLaren JR, Turkington R. Plant functional group identity differentially affects leaf and root decomposition. Global Change Biology. 2010;16(11):3075-3084
  54. 54. Jing ZB, Cheng JM, Chen A. Assessment of vegetative ecological characteristics and the succession process during three decades of grazing exclusion in a continental steppe grassland. Ecological Engineering. 2013;57:162-169
  55. 55. Olson JS. Energy storage and the balance of producers and decomposers in ecological systems. Ecology. 1963;44(2):322-331
  56. 56. Fornara DA, Tilman D, Hobbie SE. Linkages between plant functional composition, fine root processes and potential soil N mineralization rates. Journal of Ecology. 2009;97(1):48-56
  57. 57. Romero LM, Smith TJ, Fourqurean JW. Changes in mass and nutrient content of wood during decomposition in a South Florida mangrove forest. Journal of Ecobiology. 2005;93(3):618-631
  58. 58. Guo L, Dai J, Ranjitkar S, Xu J, Luedeling E. Response of chestnut phenology in China to climate variation and change. Agricultural and Forest Meteorology. 2013;180:164-172
  59. 59. Guo L, Dai J, Wang M, Xu J, Luedeling E. Responses of spring phenology in temperate zone trees to climate warming: A case study of apricot flowering in China. Agricultural and Forest Meteorology. 2015;201:1-7
  60. 60. Wold S. PLS for multivariate linear modeling. In: van der Waterbeemd H, editor. Chemometric Methods in Molecular Design: Methods and Principles in Medicinal Chemistry. Weinheim, Germany: Verlag-Chemie; 1995. pp. 195-218
  61. 61. Han YR, Kun AL, Ma ZG, Masae SYM. Forty-eight-year climatology of air temperature and precipitation changes in Xilinhot, Xilingol steppe (Inner Mongolia), China. Grassland Science. 2011;57:168-172
  62. 62. Yuan JY, Ouyang ZY, Zheng H, Xu WH. Effects of different grassland restoration approaches on soil properties in the southeastern Horqin sandy land, northern China. Applied Soil Ecology. 2012;61:34-39
  63. 63. Zhang C, Xue S, Liu GB, Song ZL. A comparison of soil qualities of different revegetation types in the loess plateau, China. Plant and Soil. 2011;347:163-178
  64. 64. Ungar I. Are biotic factors significant in influencing the distribution of halophytes in saline habitats? The Botanical Review. 1998;64(2):176-199
  65. 65. Lambers H, Mougel C, Jaillard B, Hinsinger P. Plant–microbe–soil interactions in the rhizosphere: An evolutionary perspective. Plant and Soil. 2009;321:83-115
  66. 66. Oesterheld M, Loreti J, Semmartin M, Sala OE. Inter-annual variation in primary production of a semi-arid grassland related to previous-year production. Journal of Vegetation Science. 2001;12(1):137-142
  67. 67. Sierra J. Temperature and soil moisture dependence of N mineralization in intact soil cores. Soil Biology and Biochemistry. 1997;29(9-10):1557-1563
  68. 68. Nemani RR, Keeling CD, Hashimoto H, Jolly WM, Piper SC, Tucker CJ, et al. Climate-driven increases in global terrestrial net primary production from 1982 to 1999. Science. 2003;300(5625):1560-1563
  69. 69. Bradford J, Lauenroth W, Burke I, Paruelo J. The influence of climate, soils, weather, and land use on primary production and biomass seasonality in the US Great Plains. Ecosystems. 2006;9(6):934-950
  70. 70. De Boeck HJ, Bassin S, Verlinden M, Zeiter M, Hiltbrunner E. Simulated heat waves affected alpine grassland only in combination with drought. The New Phytologist. 2016;209(2):531-541
  71. 71. Epstein HE, Lauenroth WK, Burke IC, Coffin DP. Ecological responses of dominant grasses along two climatic gradients in the Great Plains of the United States. Journal of Vegetation Science. 1996;7(6):777-788
  72. 72. Knapp A, Burns C, Fynn RS, Kirkman K, Morris C, Smith M. Convergence and contingency in production–precipitation relationships in north American and south African C4 grasslands. Oecologia. 2006;149(3):456-464
  73. 73. Koerner SE, Collins SL. Interactive effects of grazing, drought, and fire on grassland plant communities in North America and South Africa. Ecology. 2014;95(1):98-109
  74. 74. de Valpine P, Harte J. Plant responses to experimental warming in a montane meadow. Ecology. 2001;82(3):637-648
  75. 75. Dessureault-Rompré J, Zebarth BJ, Georgallas A, Burton DL, Grant CA, Drury CF. Temperature dependence of soil nitrogen mineralization rate: Comparison of mathematical models, reference temperatures and origin of thesoils. Geoderma. 2010;157(3-4):97-108
  76. 76. Bokhorst SF, Bjerke JW, Tømmervik H, Callaghan TV, Phoenix GK. Winter warming events damage sub-Arctic vegetation: Consistent evidence from an experimental manipulation and a natural event. Journal of Ecology. 2009;97(6):1408-1415
  77. 77. Yu H, Luedeling E, Xu J. Winter and spring warming result in delayed spring phenology on the Tibetan plateau. Proceedings of the National Academy of Sciences of the United States of America. 2010;107(51):22151-22156
  78. 78. Cook BI, Wolkovich EM, Parmesan C. Divergent responses to spring and winter warming drive community level flowering trends. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(23):9000-9005
  79. 79. Fu YH, Zhao H, Piao S, Peaucelle M, Peng S, Zhou G, et al. Declining global warming effects on the phenology of spring leaf unfolding. Nature. 2015;526(7571):104-107
  80. 80. Ladwig LM, Ratajczak ZR, Ocheltree TW, Hafich KA, Churchill AC, Frey SJK, et al. Beyond arctic and alpine: The influence of winter climate on temperate ecosystems. Ecology. 2016;97(2):372-382
  81. 81. Qian S, Wang LY, Gong XF. Climate change and its effects on grassland productivity and carrying capacity of livestock in the main grasslands of China. Rangeland Journal. 2012;34(4):341-347
  82. 82. Grimm NB, Chapin FS, Bierwagen B, Gonzalez P, Groffman PM, Luo Y, et al. The impacts of climate change on ecosystem structure and function. Frontiers in Ecology and the Environment. 2013;11(9):474-482
  83. 83. Luedeling E, Guo L, Dai J, Leslie C, Blanke MM. Differential responses of trees to temperature variation during the chilling and forcing phases. Agricultural and Forest Meteorology. 2013;181:33-42
  84. 84. Aerts R, Chapin FS. The mineral nutrition of wild plants revisited: A re-evaluation of processes and patterns. Advances in Ecological Research. 1999;30:1-67

Written By

Jimin Cheng, Wei Li, Jishuai Su, Liang Guo, Jingwei Jin and Chengcheng Gang

Reviewed: 19 February 2019 Published: 23 March 2019