Zongfei Wang · Fenli Zheng

Abstract Long-term natural vegetation succession plays a substantial role in the accumulation and distribution of plant and soil C:N:P stoichiometry. However, how plant and soil C:N:P relationships or ratios change along with successional stages over a century in the severely eroded areas remain unclear. These were measured over a 100-year natural succession in five successional stages from annual grasses to climax forests. The results show that natural succession had significant effects on carbon (C), nitrogen (N) and phosphorous (P) concentrations in leaf-litter-soil and their ratios in severely eroded areas. Nitrogen concentrations and N:P ratios in leaf and litter increased from annual grasses to the shrub stage and then decreased in the late successional forest stages. Leaf P levels decreased from annual grasses to shrub stages and did not significantly change during late successional stages. Litter P concentration decreased in the early successional stages and increased during late successional stages, with no overall significant change. Soil C and N concentrations and C:N, C:P and N:P ratios increased with successional stages. Soil C and N concentrations decreased with the increasing soil depth. Both were significantly different between any successional stages and controls (cropland) in the upper 10 cm and 10–20 cm soil layers. Leaf N:P ratios may be used to indicate nutrient limitations and this study suggests that plant growth during the grass stages was limited by N, during the shrub stage, by P, and during the forest stages, by both of N and P. In addition, there were close correlations between litter and leaf C:N:P ratios, soil and litter C and N levels, and C:P and N:P ratios. These results show that long-term natural vegetation succession is effective in restoring degraded soil properties and improving soil fertility, and provide insights into C:N:P relationships of leaf, litter and soil influenced by vegetation succession stage.

Keywords Ecological stoichiometry · Leaf-litter-soil · Loess Platea · Natural succession · Nutrient elements

Introduction

Degradation of the environment seriously affects human survival and the sustainable development of social economies (Roberts et al. 2009), and is of concern worldwide. Site restoration, including natural successional and afforestation approaches, is an effective measure to control soil erosion, address ecological degradation, improve environmental quality, and achieve sustainability in ecosystems (Garcia et al. 2002; Zheng 2006; Frouz et al. 2008; Jiao et al. 2012; Deng et al. 2013; Zhao et al. 2015). Secondary successional stages, which include grass stages, shrub stages and forest stages, enhance plant diversity and community structure (Letcher 2010), alter plant traits in species (Navas et al. 2010), improve soil nutrient concentrations (Davidson et al. 2007), recover soil caused by erosion (Frouz et al. 2008), increase soil carbon storage (Zhang and Shangguan 2016), and maintain soil fertility (Deng et al. 2013). Understanding the impacts of natural succession on plant and soil nutrients in severely eroded soils is of considerable importance for improving ecological function and providing guidance for reconstruction and rehabilitation of degraded ecological systems.

Ecological stoichiometry concerns the balance ofenergy and multiple chemical elements in organismal-based systems (Elser et al. 2000; Tian et al. 2010). It is also an effective and useful approach for analyzing the influence of the mass balance of multiple chemical elements in ecological process, for studying the limitations of plant growth (He et al. 2008), and interactions between soil and plants and nutrient cycling (Elser 2006). Most studies have suggested that leaf N:P ratios may be used to indicate nutritional limitations for growth (Koerselman and Meuleman 1996; Güsewell 2004). These have been widely accepted and applied to the study of terrestrial ecosystems (Reich and Oleksyn 2004; Han et al. 2005; Zhang et al. 2015). Schreeg et al. ( 2014) reported that different species have different effects on chemical relationships. Zhang et al. ( 2015) found that the successional stage was the main factor affecting leaf N:P ratios in a Karst region of China. Richardson et al. ( 2004) observed that leaf N:P ratios varied between 9.2 and 16.0, and increased along a soil chronosequence in temperate rainforests and that N and P limitations changed over the succession. Other studies also documented nutrient limitation changes with succession (Davidson et al. 2007; Selmants and Hart 2008; Laliberté et al. 2012; Zhang et al. 2015). In addition, Broadbent et al. ( 2014) found that P limitation increased during forest succession following slash-and-burn agriculture in the Bolivian Amazon, and this is similar to studies by Du et al. ( 2011) and Huang et al. ( 2013). However, McGroddy et al. ( 2008) and Rastetter et al. ( 2013) reported that, over succession, plant growth was nitrogen- limited in early stages, phosphorous-limited in middle stages and co-limited by N and P in late stages. Therefore, with ecological stoichiometry as a research tool, studies on nutrient distribution, circulation and restriction of C, N and P is of considerable value to understanding the influencing factors and mechanism of ecosystem processes.

The Loess Plateau in China, with an area of 62.4 × 10 4 km 2 , is considered as one of the most seriously eroded areas in the world (Fu et al. 2010). Historically, most of the forests were removed to meet the increasing food needs of a growing population, resulting in severe soil erosion and land degradation. The “Grain-to-Green-Program” (GTGP) was implemented in 1999 to control soil erosion and restore the ecosystems. The main goal of the GTGP was to transfer low-yielding, steep sloping croplands into grasslands, shrublands or forests (Zheng 2006; Jiao et al. 2012, 2013; Deng et al. 2013). With the development of a site restoration and succession program, plant diversity and community productivity have improved (Wang et al. 2010), soil nutrient levels have increased (Deng et al. 2013), degraded soil properties have recovered (Zhou et al. 2012), soil C storage has increased (Zhang and Shangguan 2016), and soil fertility has been maintained (Zhang et al. 2016). Several researchers have applied methods of plant or soil ecological stoichiometry to investigate nutrient limitations over vegetation succession and have shown that ecological stoichiometry can be used to assess the effects of site restoration in degraded ecosystems (Davidson et al. 2007; An and Shangguan 2010; Du et al. 2011; Laliberté et al. 2012; Huang et al. 2013; Jiao et al. 2013; Rastetter et al. 2013; Chai et al. 2015; Zhang et al. 2015; Zhang and Shangguan 2018). However, plant components and soil ecological relationship changes during secondary successional stages are not particularly clear (Zeng et al. 2017), which limits understanding of nutrient cycling in terrestrial ecosystems.

In recent decades, a wide range of secondary forest landscapes have formed in the Ziwuling area of the Loess Plateau, which is the only natural secondary forest on the Loess Plateau. Studying the long-term natural site succession (100 years or more) can enhance understanding of succession effects on nutrient relationships in degraded ecosystems. In this study, leaf-litter-soil nutrient stoichiometry was investigated across a century of natural succession with five successional stages from annual grasses to climax forest communities. The specific objectives were: (1) to investigate changes in C, N and P concentrations and their stoichiometry in leaf-litter-soil by successional stages; (2) to clarify the limiting nutrient factors of plant growth at each stage; (3) to determine soil profile distribution of nutrients at each stage; and, (4) to analyze the correlations of C, N and P concentrations and their stoichiometry between leaf, litter and soil across the successional stages. The results can provide theoretical basis for planning vegetation restoration and improving the ecological stoichiometry of degraded ecosystems on Chinese Loess Plateau.

Materials and methods

Study site

Fig.1 Locations of the sampling plots

This study was conducted at the Fuxian Observatory for Soil Erosion and Eco-environment (35° 44.1′–36° 23.4′ N, 108° 29.5′–109° 42.9′ E) located in the center of the Ziwuling natural secondary forest of Chinese Loess Plateau (Fig.1). The study area is situated in a loess hilly gully region with elevations ranging from 920 to 1683 m (Zheng 2006). The mean annual temperature is 7.4 °C and average annual rainfall is 576.7 mm, of which 60% falls in the rainy season which lasts from June to September. The air relative humidity is 65–70%. The dominant soil is a fine-grained loess, which can be classified as a Calcic Cambisol (USDA NRCS 1999). Soil texture is 28.3% sand (> 50 μm), 58.1% silt (50–2 μm), 13.6% clay (< 2 μm), and the thickness of loess in this region is more than 50 m. Forest canopy closure is more than 0.6 in this area. According to historical records and previous research (Zheng et al. 1997; Kang et al. 2014), most forests were removed and the land cultivated over 100 years ago. Croplands occupied 30% of the area, and the estimate of soil loss was about 8000–10,000 t km 2 a?1. The ephemeral or transitory gully density on hillsides was 20–50 km km?2, which reflects severe soil erosion over time. Moreover, the organic horizon was completely removed and loessal parent material exposed to the surface. In 1866–1870, a civil war in this region dispersed populations, and vegetation began to recover and eventually formed a local climax forest largely of Liaotung oak ( Quercus liaotungensis Koidz). However, since 1949, there have been several periods of vegetation disturbance, such as the late 1950 (timber and fuelwood needs), 1960–1962 (food demands due to three year’s drought disaster), 1966–1976 (the 10-Year Cultural Revolution), and the 1980s (implementing household contract responsibility system) (Wang and Tang 1992). Therefore, at present there are different stages of successions especially, since more than a century has passed; large areas of low-yielding crop filed on steep slopes have converted into grasslands, shrubs and pioneer tree species, and climax forest communities. Furthermore, main pioneer forest is Betula platyphylla Sukaczev; there are two dominant shrub species Sophora viciifolia (Franch.) Skeels and Hippophae rhamnoides Linn., but both does not concur in the same location (Zheng 2006). The principal herb species are Artemisia scoparia Waldst.et (Kit.), Artemisia sacrorum (Ledeb.) and Bothriochloa ischaemun (Linn.) (An and Shangguan 2010; Chai et al. 2015). Accordingly, these dominate species were selected to study the impact of vegetation succession on leaf-littersoil C:N:P stoichiometry and their intrinsic relationships.

Soil and plant sampling

In this study, the approach of substituting “space” for “time” was used, which has been adopted by several researchers (Sparling et al. 2003; An and Shangguan 2010; Jiao et al. 2013; Zhang et al. 2015; Zeng et al. 2017). To select communities representing different successional stages and to determine their ages included the following steps. First, based on previous research (An and Shangguan 2010; Chai et al. 2015) and further field investigation, five stages were selected: annual grass stage (G1), perennial grass stage (G2), shrub stage (S), pioneer forest stage (F1), and climax forest stage (F2). Second, ages ofeach stage were identified. The early stage (< 20 years) was easily determined because the “Grain-to-Green” program started in 1999. Identifying the intermediate stage (20–75 years) and later stage (> 75 years) included the following steps: (1) the successional age for each site was determined by field investigation and vegetation type; local farmers provided land contracts between themselves and local communities in order to estimate how long farmlands were abandoned to vegetation restoration; (2) past and current village heads were interviewed for information on land use history ofeach sampling site and to further confirm the successional ages; and, (3) all information was combined to determine the age ofeach successional stage. The five vegetation successional stages were approximately 2-, 20-, 55-, 75- and 100-years-old, respectively. Cropland near the sample plots was selected as a control. For each stage or sere, three experimental sites with similar topographical conditions were chosen for sample collection. In addition, the distance between all experimental sites was approximately 15 km, which reduced the impacts of previous site conditions on the successional process. Leaves and litter ofonly the dominant species in each stage were collected. Artemisia scoparia Waldst.et (Kit.), was dominant in G1, Artemisia sacrorum (Ledeb.) in G2, Sophora viciifolia (Franch.) Skeels in S, Betula platyphylla Sukaczev in F1 and Quercus liaotungensis Koidz in F2. All shrub and tree species are deciduous. Leaves and soil samples were obtained in late July 2016, when the plants were in a vigorous growth period, and litter samples that lay on the soil surface over multiple years but had not decomposed were collected in late October 2016. Table 1 shows the characteristics of the sample plots of the different successional stages.

For F1 and F2 studies, plots were 20 m × 20 m; for S, plots were 5 m × 5 m, and for G1 and G2, plots were 1 m × 1 m. Each plot had three replicates. Fully expanded leaves were randomly collected from five to ten healthy plants in each plot for the shrub and tree species, and then 10–20 leaves were collected from each shrub or tree, for a total of 50–200 leaves from each plot. For the grasses, stems and leaves were cut from three 0.25-m 2 sampling sites. All leaf samples were evenly mixed and bagged. Litter samples for each plot were collected from three 1 m × 1 m squares along diagonal lines, and all litter from these squares was mixed and stored. All samples were transported to the laboratory. Soil samples from the 100-cm profile were collected using a 5-cm diameter soil sampler along an S-shaped route in each plot. Before each sample was collected, the sampler was sterilized with alcohol to ensure samples were free from cross-infection. In addition, the 100-cm profile was divided into six layers (0–10, 10–20, 20–40, 40–60, 60–80, and 80–100 cm) and samples taken from each layer at five points. Soil samples ofeach layer were then mixed evenly, stored in a plastic bag and taken to the laboratory.

Table 1 Descriptions of the sample plots at different successional stages

Sample analysis

Plant samples were oven- dried at 70 °C for 48 h to a constant mass and then weighed. The samples were ground to a uniform fine powder using a plant-sample mill (1093 Sample Mill, Foss, Sweden). Soil samples were air- dried and sieved through a 0.25 mm sieve. The Walkley–Black modified acid–dichromate FeSO4titration method measured organic C concentrations of leaf, litter and soil samples (Bao 2000). The Kjeldahl method (KJELTE2300, Sweden) determined total N concentrations of the samples and total P levels were quantified by the standard ammonium molybdate method (Bao 2000). Each sample analysis was repeated three times, and the absolute error indicator (less than 5%) was applied for quality control. C, N and P concentrations were expressed as g kg?1on a dry-weight basis. C: N: P ratios in leaf, litter and soil were calculated as mass ratios.

Statistical analysis

All results were reported as the means ± standard errors (SE) and tested for normality of distributions and homogeneity of variances before analysis. One-way ANOVA determined the effects of successional stages on C, N and P concentrations and their stoichiometry of leaves, litter and soil. Two-way ANOVAs analyzed the effects of successional stage, soil depth and their interactions on soil C, N, and P concentrations and their stoichiometry. Pearson correlation was used to test the relationship between C, N and P concentrations in leaf, litter and soil. A linear regression analysis tested the relationship between leaf, litter and soil stoichiometric characteristics. All significant differences were shown at P < 0.05. SPSS 19.0 software (SPSS, Inc., Chicago, IL, USA) was used to perform the statistical analysis.

Results

Response of leaf and litter C, N and P concentrations to successional stages

C, N, and P levels differed significantly with different successional stages (Fig.2). Concentrations in leaves were higher than in litter (Fig.2 a–c). In general, leaf C concentrations increased with successional stage, and in the 75-year-old pioneer forest stage, leaf C was significantly higher than in other stages. Litter C levels decreased with successional stages and in the 55-year-old shrub stage was significantly lower than in other stages (Fig.2 a). Leaf N levels firstly increased and then later decreased, and was highest in the shrub community (55-year-old) and lowest in the annual grass stage (2-year-old) (Fig.2 b), but it was not significantly different among G2, F1 and F2 stages (Fig.2 b). The dynamic changes of N levels in litter was similar to that of leaf N concentrations (Fig.2 b). P levels in leaves significantly decreased over the early seral stages and did not significantly changed in the late successional stages (Fig.2 c). Litter P concentrations decreased in early grass stages and increased in subsequent stages. In the 20-year-old perennial grass community stage, P concentrations were significantly lower than in other stages (Fig.2 c).

Changes of leaf and litter C, N and P stoichiometric characteristics with successional stages

The stage of succession had considerable effects on C:N, C:P and N:P ratios in leaves and litter (Fig.2 d–e). C:N ratios in leaf and litter significantly decreased with successional stages, reaching the lowest value in the 55-year-old shrub stage and subsequently increased; were highest in the 2-yearold annual stage (Fig.2 d). Litter C:P ratios were higher than that leaf C:P ratios (Fig.2 e). These increased with the stage of successional, and no significant difference was observed in the 55-year-old (S), 75-year-old (F1) and 100-year-old ages (F2) stages, while litter C:P ratios increased from G1 to G2 and then decreased (Fig.2 e). N:P ratios in leaves and were opposite to C:N ratios in leaves and litter, first increasing in early succession, becoming significantly higher in the shrub stage and then decreasing (Fig.2 f). In general, patterns of leaf C:N:P stoichiometry increased along stage of succession (Table 2).

Distribution of soil C, N and P concentrations in soil profiles

Figure 3 shows changes of soil C, N and P concentrations with soil depths. Over five successional stages and the cropland control, C levels were highest in the upper 10 cm layer, then decreased with increasing depth, remaining stable below 40 cm (Fig.3 a). The distribution of soil N was similar to soil C levels (Fig.3 b). Concentrations of soil P were relatively stable with increasing soil depth (Fig.3 c). Compared with the control treatment, over five successional stages soil C and N concentrations in both the 0–10 and 10–20 cm layers were significantly higher and increased with stage of succession. At the 20–40 cm depth, C and N levels were not significantly different between the controls and any succession stage, but increased with stage of succession. Bellow 40 cm, C and N levels were stable. Consequently, changes of soil nutrient concentrations and their stoichiometry with successional stages at 0–10, 0–20 and 0–40 cm soil depths were analyzed in late section.

Fig.2 Change of stoichiometry characteristics of C, N and P in leaf and litter over successional stages. All data are in the form of mean ± SE. Different lowercase letters indicate significant differences at different vegetation successional stages whether in leafor in litter ( P < 0.05)

Table 2 C:N:P ratios in leaves, litter and soil at different successional stages

Fig.3 Vertical profile distribution of soil C, N and P concentrations at different successional stages. CL: crop land; all data are in the form of mean ± SE

Change of soil C, N and P concentrations with successional stages

Soil C concentrations in the 0–10, 0–20 and 0–40 cm layers significantly increased with stage of succession (Fig.4 a). Soil C concentrations over successional stages decreased with increasing of soil depth, ranging from 9.6 to 21.9 g/kg in the 0–10 cm layer, 7.6 to 16.8 g/kg in the 0–20 cm layer and 5.6 to 10.8 g/kg in the 0–40 cm layer. Soil C concentration was highest in the F2 (climax forest community) stage and lowest the controls in all three depths (Fig.4 a). Changes in soil N with successional stages were similar to soil C (Fig.4 b), ranging from 0.8 to 1.8 g/kg in the upper 10 cm layer, 0.8 to 1.4 g/kg in 0–20 cm soil depth and 0.6 to 1.0 g/kg in 0–40 cm soil depth. It was highest in the F2 (climax forest) stage and lowest for the controls in these three soil depths. However, soil P concentrations in the different soil depths (0–10, 0–20 and 0–40 cm) first decreased from the G1 and G2 grass stages to the shrub stages, increased from the S stage to the F1 (pioneer forest stage), and decreased in the climax forest (Fig.4 c).

Changes of soil C, N and P stoichiometric characteristics with successional stages

The dynamic soil C, N and P stoichiometry changes with vegetation succession in different soil layers are shown in Fig.4 d and f. The C:N ratios in three soil depths (0–10, 0–20 and 0–40 cm) numerically increased with successional stages, and were highest in the climax forest community stage (Fig.4 d). Soil C:P ratios in the 0–10 cm soil depth was significantly higher than that in the 0–20 and 0–40 cm soil depths (Fig.4 e). Soil C:P ratios increased with successional stages and were highest in the climax forest stage in these three soil depths. Changes of N:P ratios by soil depth were similar to changes in soil C:P ratios. In the 0–10 cm depth, N:P ratios in F2 stage were significantly higher than in other stages, but in the 0–20 and 0–40 cm soil depths, there were no significant differences between S and F2 stages (Fig.4 f). Soil C:N:P stoichiometry increased over successional stages and were highest in F2 stage and lowest in the controls (Table 2).

The results of two-way ANOVA indicated that both successional stage and soil depth significantly affected soil C, N and P concentrations and their ratios. The interactions between successional stage and soil depth significantly affected soil C and N concentrations, C:N, C:P and N:P ratios but not P concentration (Table 3).

Relationships between plant and soil nutrients and their stoichiometry

There were significant correlations between leaf and litter C, N and P concentrations (Fig.5). Meanwhile, C:N ratios in leaves increased linearly with an increase of litter C:N ratios (R 2 = 0.82 and P < 0.001) (Fig.5 d). Litter C:P ratios were linearly correlated with leaf C:P ratios (R 2 = 0.24 and P < 0.05) (Fig.5 e). In addition, litter N:P ratios increased linearly with the increase of N:P ratios (R 2 = 0.72 and P < 0.001) (Fig.5 f).

Fig.4 Stoichiometry characteristics of soil C, N and P at 0–10, 0–20 and 0–40 cm soil depths by successional stages. CL: crop land; All data are in the form of mean ± SE; Different lowercase letters indicate significant differences at different vegetation successional stages in the same soil depth ( P < 0.05)

Based on the above results of soil nutrients profile distribution and changes of soil nutrients with successional stages, the relationship between soil nutrients in 0–20 cm soil depth and leaf and litter nutrients were analyzed. There were significant correlations between litter and soil C concentration in 0–20 cm soil depth (R 2 = 0.20 and P < 0.05) (Fig.6 a). Litter N concentrations were correlated with soil N concentration (R 2 = 0.20 and P < 0.05) (Fig.6 b). However, the relationship between litter and soil P concentration was not significant (R 2 = 0.061 and P = 0.189) (Fig.6 c). Meanwhile, litter C:N ratios were not correlated with soil C:N ratios (R 2 = 0.05 and P = 0.247) (Fig.6 d), and litter C:P ratios were linearly correlated with soil C:P ratios (R 2 = 0.26 and P < 0.05) (Fig.6 e). Additionally, litter N:P ratios increased linearly with an increase in soil N:P ratios (R 2 = 0.63 and P < 0.001) (Fig.6 f). There was a significant positive correlation between leaf and soil C concentrations in the 0–20 cm soil depth (R 2 = 0.62 and P < 0.001) (Fig.7 a), but the relationships between soil andleaf N and P concentrations were insignificant ( P = 0.054, Fig.7 b) ( P = 0.727, Fig.7 c). Leaf C:N ratios were not correlated with soil C:N ratios ( P = 0.785) (Fig.7 d). In addition, soil C:P ratios increased linearly with an increase in leaf C:P ratios (R 2 = 0.75 and P < 0.001) (Fig.7 e) and leaf N:P ratios were positively correlated with soil N:P ratios (R 2 = 0.45 and P < 0.001) (Fig.7 f).

Discussion

The effect of successional stages on plant C, N and P concentrations and C:N:P ratios

As a major component terrestrial ecosystem, plants have a vital function in stabilizing the whole system. C, N, P and other elements are essential for plant development (John et al. 2007). Leaf stoichiometric characteristics affect plant growth and reflect their internal characteristics and growth state (Koerselman and Meuleman 1996). Additionally, leaf stoichiometry can reveal the adaptation of plants to the external environment (McIntyre et al. 1999). Plant community succession is a result of the substitution of species leading to changes in community structure, composition and function (Raevel et al. 2012). Due to species composition in various successional stages, nutrient utilization strategies are different, resulting in significant differences in nutrient concentrations and N:P ratios in leaves (Navas et al. 2010). In this study, the differences of dominant species in the five successional stages resulted in differences in C, N and P concentrations and C:N:P stoichiometry in leaves. Leaf P concentrations in the G1 grass stage was higher than that in the S (shrub) and F1 and F2 forest stages. These results are consistent with studies by Chai et al. ( 2015) and Zhang et al. ( 2015). Previous studies showed that the rapid growth of plants requires the synthesis of large amounts of protein and RNA, and this requires a high P concentration (Elser et al. 2003). Kazakou et al. ( 2006) also observed that in the succession of abandoned farmlands in the Mediterranean area, early successional species had high P concentration and relatively high growth rate, and late successional species had longer leaf life. He et al. ( 2008) also reported that P concentrations in woody plants were significantly lower than in herbaceous plants. Reich and Oleksyn ( 2004) and Sardans and Pe?uelas ( 2013) found that plants with long life cycles and slow growth rates had lower P concentrations than those with shorter life cycles and rapid growth rates. Leaf N concentrations in the shrub stage were significantly higher than that in other stages, because this stage is characterized by Sophora viciifolia (Franch.) Skeels, a nitrogen-fixing species that utilizes atmospheric N2. Our study has reported that the concentrations of C, N and P in litter were lower than those of plant leaves over all successional stages. This is consistent with the results of Zeng et al. ( 2017). A possible reason is that leaves transfer nutrients to other plant tissues before senescence (Schreeg et al. 2014). In addition, this disbursement improves nutrient utilization rates and adaptation to the environment (Drenovsky and Richards 2004). In this study, litter N concentrations first increased and then decreased during the succession process. The change in litter N was similar to leaf N concentrations. Takashima et al. ( 2004) reported that rapid plant growth requires more nitrogen consumption to participate in photosynthesis, and higher N transfer rates leads to the reduction of N in litter. In this study, higher growth rates by annual grasses may be the reason that litter N concentrations were lower than that in other stages.

Fig.5 Correlations between leaf and litter C:N:P stoichiometric characteristics by successional stages

Research by Güsewell ( 2004) indicated that deficiencies in N and P can be reflected in changes in plant growth. Therefore, leaf N:P ratios can be used as one of the indices to estimate nutrient limitations (Koerselman and Meuleman 1996; Güsewell 2004). Güsewell ( 2004) reported that the threshold of leaf N:P ratios were 10 and 20: when the N:P ratio was > 20, plant growth was P-limited, and when the ratio was < 10, plant growth was N-limited. When the N:P ratio was between 10 and 20, plant growth was limited by both N and P. Thus, based on the Güsewell ( 2004), the limiting nutrient factor in each successional stage has been identified in this research. In the annual and perennial grass stages (G1 and G2), plant growth was N-limited, whereas it was P-limited in the shrub stage (s). In the pioneer and climax forest stages (F1 and F2), growth was limited by both N and P. This is consistent with the results reported by Rastetter et al. ( 2013). There are several possible reasons for this. The first is related to utilization approaches to obtaining N and P elements. Soil P is mainly derived from the weathering of rocks, which is a long-term biogeochemical cycle, and its utilization is influenced by soil pH value (Turrión et al. 2000). There are various ways to acquire N elements (Reich and Oleksyn 2004), and N in litter is easier to mineralize and release than that of P (Güsewell 2004). In early successional stages, the release of P element from soil parent material is a result of surface water erosion, making it easier accumulate P than later successional stages due to increased vegetation cover. As succession proceeds, the demand for nutrients also increase, resulting in insuffi cient P, which is caused by the lack rate of P to be significantly greater than that of N. Secondly, the shrub stage is characterized by Sophora viciifolia (Franch.) Skeels, a nitrogen-fixing plant and can utilize atmospheric N2. Thus, leaf N:P ratios in mid successional stage was higher than that in early grass stage and forest stage. Although leaf N:P ratios may effectively reflect N or P limitations, the importance of the N:P ratios is mainly in its function as an indicator (Güsewell 2004). In addition, a fertilization experiment should be conducted in the near future to identify limiting nutrient factors and verify the conclusion from this study.

Fig.6 Correlations between litter and soil C:N:P stoichiometric characteristics at 0–20 cm soil depth by successional stages

The effect of successional stages on soil C, N and P concentrations and C:N:P stoichiometry

Fig.7 Correlations between leaf and soil C:N:P stoichiometric characteristics at 0–20 cm soil depth by successional stages

Site restoration by vegetation changes physical, chemical and biological soil properties, especially the circulation of the three important nutrients (C, N and P) in the soil (Wei et al. 2009). Converting farmland to grassland and ultimately forestland had considerable effects in reducing erosion, improving ecological environments and enhancing soil quality (Garcia et al. 2002; Zheng 2006; Jiao et al. 2012). Long-term site restoration will increase soil nutrients, and significantly improve soil quality (Deng et al. 2013; Jiao et al. 2013; Zhao et al. 2015). Fu et al. ( 2010) observed that the restoration of vegetation increased C and N fixation because long-term site restoration reduced the loss of these elements caused by soil erosion. In this study, the results indicated that the soil C and N concentrations increased significantly over successional stages, this is consistent with the results of Deng et al. ( 2013). The main reason is that the composition, structure and level of the community gradually became more complex and diversified, and the enhanced activities of soil animals and microorganisms accelerated litter decomposition which increased the input oforganic compounds, which in turn, led to high C and N concentrations. This study showed that P concentrations were different at different successional stages, decreasing from G1 and G2 stage and slowly increasing to the F1 stage, which is not in agreement with Zhang et al. ( 2015). In addition, C and N concentrations decreased with increasing depth, and in the top soil layer (0–10 cm) were higher than those in the lower layers, while soil P concentrations were relatively stable with increasing depth. This is mainly due to different sources of C, N and P. Soil C and N concentrations were affected by decomposition of the litter, and the nutrient return first occurs at the soil surface and leaches downward, while P concentration mainly influenced by weathering of parent materials and its low mobility, making it less variable (Wei et al. 2009).

C, N and P are the three major elements of living organisms, and their ratios are an important trait oforganisms and ofecosystems. Soil C:N:P ratio are an important index oforganic matter composition, soil quality and nutrient supply capacity. A well-constrained global C:N:P ratio is 186:13:1 at a depth of 10 cm (0–10 cm organic-rich soil) (Cleveland and Liptzin, 2007). In China, the average soil C:N:P ratio is 134:9:1 for surface soil (0–10 cm mineral soil) (Tian et al. 2010). In this study, the highest C:N:P ratio was 31.9:2.7:1.0 for the upper 10 cm organic-rich soil (Table 2), lower than the average for China and the global mean. This is due to the serious erosion in the study area which results in considerable surface carbon loss and low C:N:P ratios. Soil C:N:P ratios increased with successional stages. A possible reason is that the absorption and release of root system secretions changes C:N:P ratio; additionally, N and P released by litter decomposition in forestland affects the C:N:P ratio. Among these ratios, the C:N ratio was a sensitive indicator of soil quality and could affect C and N cycling. In this study, there was a significant positive correlation between C and N concentrations (R 2 = 0.93), and the C:N ratio was approximately 10.2 on the Loess Plateau, similar to the average value in China (10.0–12.0) (Tian et al. 2010) but lower than the global average (13.3) (Cleveland and Liptzin 2007). However, the result was higher than the reported by Liu et al. ( 2010) (C:N = 7.62). Tian et al. ( 2010) observed that C and N concentrations had wide spatial variability, whereas C:N ratios were relatively stable. Our results are consistent with studies that show that C and N were closely related to plant structural components and that the C:N ratio was relatively stable ratio in the processes of accumulation and consumption (Agren 2008). The soil N:P ratio may be used as an effective index of nutrient limitation (Tessier and Raynal 2003). The maximum value in this study of the N:P ratio was 2.7, lower than the 5.2 average in China (Tian et al. 2010). This is because the soil was alkaline and concentrations of soil N are low on the Loess Plateau. Therefore, plant growth was limited by N or by both N and P. In addition, the N:P ratio was lower than the global value (5.9) (Cleveland and Liptzin 2007). The low ratio also shows that N concentrations were relatively low in the region. In addition, the N:P ratio in the climax forest stage was higher than in other stages; the reason is that from the annual grass to climax forest stage, the increase of N concentration was greater than that of P concentration.

Relationship between soil and plant stoichiometry

Plants sequester C through photosynthesis and gradually add it to the soil in the form of litter, while soil, as the medium of plant growth, provides all plant nutrients. The availability of soil nutrients determines the concentrations of chemical elements in plants. A better understanding was required to determine how nutrient cycles interact among plants, litter and soil during the successional stages in terrestrial ecosystems. The data presented here show a significant correlation between leaves and litter for C, N and P concentrations and their ratios, Lambers et al. ( 2010) reported correlations between N and P concentrations in soils and leaves in south-western Australia. Parfitt et al. ( 2005) also reported that soil N concentrations were positively correlated with leaf N concentration in temperate rain forests over three vegetation successional stages in New Zealand. Research shows that P concentration in leaves was closely related to P concentration in soil (Hedin 2004). However, other studies have reported that there were no relationships between soil and forest plants in terms of C, N and P concentrations (Ladanai et al. 2010; Yu et al. 2010). Significant positive correlations were also found between soil and leaves for C:N:P ratios in Australia’s major native vegetation ecosystems (Bui and Henderson 2013), and Bell et al. ( 2014) found a strong correlation between soil and leaves for C:N:P ratios in semi-arid grassland ecosystems. In this study, the data indicate that leaf N and P concentrations have no correlations with soil N and P concentrations, a finding consistent with Ladanai et al. ( 2010). This is because of the absorption of mineral elements by plants and their transport, storage, and utilization area affected by many factors, including soil temperatures, water, microbial activities, and interspecific competition (Güsewell 2004; Reich and Oleksyn 2004). In addition, the species in the study region have long adapted to their habitat, and concentrations of N and P in leaves may be more influenced by the characteristics of the species themselves than by nutrients in the soil. However, there were close relations between leaf N:P and C:P ratios and soil N:P and C:P ratios. The results are consistent with those previous studies (Bui and Henderson 2013; Bell et al. 2014; Zeng et al. 2017).

In general, correlations among leaf, litter and soil C, N and P concentrations and C:N:P stoichiometry reflect, to some degree, the inheritance and co-invariance of nutrient concentrations. However, factors such as vegetation type and age, and soil depth result in differences in the correlations between leaf, litter and soil in terms of the C, N and P concentrations and in the C:N:P stoichiometric characteristics, reflecting a flexible adaptation mechanism of plants to environmental change.

Conclusion

The results show that changes of leaf-litter-soil C, N and P concentrations and their stoichiometric characteristics were different with the successional stages. Leaf N:P ratios can be used to indicate nutritional limits for plant growth and this study found that plant growth in the grass stage was limited by N, while in the shrub stage, growth was limited by P, and in the forest stages, it was restricted by both of N and P. There were significant correlations between leaf and litter nutrients ( P < 0.05) and leaf N concentrations were significantly correlated with litter N concentrations ( P < 0.05). However, there were no significant correlations between N and P concentrations in leaves and soil ( P = 0.054 and P = 0.727). The results provide a basis for the management of degraded ecosystems on the Chinese Loess Plateau.

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons.org/licen ses/by/4.0/.

References

Agren GI (2008) Stoichiometry and nutrition of plant growth in natural communities. Annu Rev Ecol Evol Syst 39(1):153–170

An H, Shangguan ZP (2010) Leaf stoichiometric trait and specific leaf area of dominant species in the secondary succession of the Loess Plateau. Pol J Ecol 58:103–113

Bao SD (2000) Soil and agricultural chemistry analysis. China Agriculture Press, Beijing, pp 30–78 (in Chinese)

Bell C, Carrillo Y, Boot CM, Rocca JD, Pendall E, Wallenstein MD (2014) Rhizosphere stoichiometry: are C:N:P ratios of plants, soils, and enzymes conserved at the plant species-level? New Phytol 201:505–517

Broadbent EN, Zambrano AMA, Asner GP, Soriano M, Field CB, de Souza HR, Pena-Claros M, Adams RI, Dirzo R, Giles L (2014) Integrating stand and soil properties to understand foliar nutrient dynamics during forest succession following slash-and-burn agriculture in the Bolivian Amazon. Plos One 9(2):e86042

Bui EN, Henderson BL (2013) C:N: P stoichiometry in Australian soils with respect to vegetation and environmental factors. Plant Soil 373:553–568

Chai YF, Liu X, Yue M, Guo JC, Wang M, Wan PC, Zhang XF, Zhang CG (2015) Leaf traits in dominant species from different secondary successional stages of deciduous forest on the Loess Plateau of northern China. Appl Veg Sci 18:50–63

Cleveland CC, Liptzin D (2007) C:N:P stoichiometry in soil: is there a “Redfield ratio” for the microbial biomass? Biogeochemistry 85:235–252

Davidson EA, Carvalho CJ, Figueira AM, Ishida FY, Ometto JP, Nardoto GB, Saba RT, Hayashi SN, Leal EC, Vieira IC, Martinelli LA (2007) Recuperation of nitrogen cycling in Amazonian forests following agricultural abandonment. Nature 447:995–998

Deng L, Wang KB, Chen ML, Shangguan ZP, Sweeney S (2013) Soil organic carbon storage capacity positively related to forest succession on the Loess Plateau, China. CATENA 110:1–7

Drenovsky RE, Richards JH (2004) Critical N:P values: predicting nutrient deficiencies in desert shrublands. Plant Soil 259(1–2):59–69

Du Y, Pan G, Li L, Hu Z, Wang X (2011) Leaf N/P ratio and nutrient reuse between dominant species and stands: predicting phosphorus deficiencies in Karst ecosystems, southwestern China. Environ Earth Sci 64:299–309

Elser J (2006) Biological stoichiometry: a chemical bridge between ecosystem ecology and evolutionary biology. Am Nat 168:S25–S35

Elser JJ, Sterner RW, Gorokhova E, Fagan WF, Markow TA, Cotner JB, Harrison JF, Hobbie SE, Odell GM, Weider LJ (2000) Biologi

cal stoichiometry from genes to ecosystems. Ecol Lett 3:540–550

Elser JJ, Acharya K, Kyle M, Cotner W, Makion W, Makion T, Wstts T, Hobbie S, Fagan W, Schade J, Hood J, Sterner RW (2003) Growth rate-stoichiometry couplings in diverse biota. Ecol Lett 6(10):936–943

Frouz J, Prach K, Pi?l V, Háněl L, Stary J, Tajovsky K, Materna J, Balík V, Kal?ík J, ?ehounková K (2008) Interactions between soil development, vegetation and soil fauna during spontaneous succession in post mining sites. Eur J Soil Biol 44(1):109–121

Fu XL, Shao MA, Wei XR, Horton R (2010) Soil organic carbon and total nitrogen as affected by vegetation types in Northern Loess Plateau of China. Geoderma 15:31–35

Garcia C, Hernandez T, Roldan A, Martin A (2002) Effect of plant cover decline on chemical and microbiological parameters under Mediterranean climate. Soil Biol Biochem 34:635–642

Güsewell S (2004) N:P ratios in terrestrial plants: variation and functional significance. New Phytol 164:243–266

Han WX, Fang JY, Guo DL, Zhang Y (2005) Leaf nitrogen and phosphorus stoichiometry across 753 terrestrial plant species in China. New Phytol 168:377–385

He JS, Wang L, Flynn DFB, Wang X, Ma W, Fang J (2008) Leaf nitrogen: phosphorus stoichiometry across Chinese grassland biomes. Oecologia 155:301–310

Hedin LO (2004) Global organization of terrestrial plant-nutrient interactions. 2004. PNAS 101:10849–10850

Huang W, Liu J, Wang Y, Zhou G, Han T, Li Y (2013) Increasing phosphorus limitation along three successional forests in southern China. Plant Soil 364:181–191

Jiao JY, Zhang ZG, Bai WJ, Jia YF, Wang N (2012) Assessing the ecological success of restoration by afforestation on the Chinese Loess Plateau. Restor Ecol 20:240–249

Jiao F, Zhong ZM, An SS, Yuan Z (2013) Successional changes in soil stoichiometry after land abandonment in Loess Plateau, China. Ecol Eng 58:249–254

John R, Dalling JW, Harms KE, Yavitt JB, Stallard RF, Mirabello M, Hubbell SP, Valencia R, Navarrete H, Vallejo M (2007) Soil nutrients influence spatial distributions of tropical tree species. PNAS 104:864–869

Kang D, Guo Y, Ren C, Zhao F, Feng Y, Han X, Yang G (2014) Population structure and spatial pattern of main tree species in secondary Betula platyphylla forest in Ziwuling Mountains, China. Sci Rep 4:1–8

Kazakou E, Vile D, Shipley B, Gallet C, Garnier E (2006) Co-variations in litter decomposition, leaf traits and plant growth in species from a Mediterranean old-field succession. Funct Ecol 20(1):21–30

Koerselman W, Meuleman AFM (1996) The vegetation N:P ratio: a new tool to detect the nature of nutrient limitation. J Appl Ecol 33:441–1450

Ladanai S, Agren GI, Olsson BA (2010) Relationships between tree and soil properties in Picea abies and Pinus sylvestris forests in Sweden. Ecosystems 13:302–316

Laliberté E, Turner BL, Costes T, Pearse SJ, Wyrwoll K-H, Zemunik G, Lambers H (2012) Experimental assessment of nutrient limitation along a 2-million-year dune chronosequence in the south-western Australia biodiversity hotspot. J Ecol 100:631–642

Lambers H, Brundrett MC, Raven JA, Hopper SD (2010) Plant mineral nutrition in ancient landscapes: high plant species diversity on infertile soils is linked to functional diversity for nutritional strategies. Plant Soil 334:11–31

Letcher SG (2010) Phylogenetic structure of angiosperm communities during tropical forest succession. Proc R Soc B Biol Sci 277(1678):97–104

Liu X, Li FM, Liu DQ, Sun GJ (2010) Soil organic carbon, carbon fractions and nutrients as affected by land use in semi–arid region of Loess Plateau of China. Pedosphere 20:146–152

McGroddy ME, Baisden WT, Hedin LO (2008) Stoichiometry of hydrological C, N, and P losses across climate and geology: an environmental matrix approach across New Zealand primary forests. Global Biogeochem Cycles 22(1):GB1026

McIntyre S, Lavorel S, Landsberg J, Forbes TDA (1999) Disturbance response in vegetation—towards a global perspective on functional traits. J Veg Sci 10:621–630

Navas ML, Roumet C, Bellmann A, Laurent G, Garnier E (2010) Suites of plant traits in species from different stages of a Mediterranean secondary succession. Plant Biol 12:183–196

Nrcs USDA (1999) Soil taxonomy: a basic system of soil classification for making and interpreting soil surveys. Agricultural handbook 436, 2nd edn. U.S. Government Printing Offi ce, Washington DC

Parfitt RL, Ross DJ, Coomes DA, Richardson SJ, Smale MC, Dahlgren RA (2005) N and P in New Zealand soil chronosequences and relationships with foliar N and P. Biogeochemistry 75:305–328

Raevel V, Violle C, Munoz F (2012) Mechanisms ofecological succession: insights from plant functional strategies. Oikos 121:1761–1770

Rastetter EBR, Yanai RD, Thomas RQ, Vadeboncoeur MA, Fahey TJ, Fiak MC, Kwiatkowshi BL, Hamburg SP (2013) Recovery from disturbance requires resynchronization ofecosystem nutrient cycles. Ecol Appl 23(23):621–642

Reich PB, Oleksyn J (2004) Global patterns of plant leaf N and P in relation to temperature and latitude. PNAS 101:11001–11006

Richardson S, Peltzer D, Allen R, McGlone M, Parfitt R (2004) Rapid development of phosphorus limitation in temperate rainforest along the Franz Josef soil chronosequence. Oecologia 139:267–276

Roberts L, Stone R, Sugden A (2009) The rise of restoration ecology. Science 325(5940):555

Sardans J, Pe?uelas J (2013) Tree growth changes with climate and forest type are associated with relative allocation of nutrients, especially phosphorus, to leaves and wood. Global Ecol Biogeogr 22(4):494–507

Schreeg LA, Santiago LS, Wright SJ, Turner BL (2014) Stem, root, and older leaf N:P ratios are more responsive indicators of soil nutrient availability than new foliage. Ecology 95:2062–2068

Selmants PC, Hart SC (2008) Substrate age and tree islands influence carbon and nitrogen dynamics across a retrogressive semiarid chronosequence. Global Biogeochem Cycles 22(1):GB1021

Sparling GP, Schipper LA, Bettjeman W, Hill R (2003) Soil quality monitoring in New Zealand: practical lessons from a 6-year trial. Agric Ecosyst Environ 104(3):523–534

Takashima T, Hikosaka K, Hirose T (2004) Photosynthesis or persistence: nitrogen allocation in leaves ofevergreen and deciduous Quercusspecies. Plant Cell Environ 27(8):1047–1054

Tessier JT, Raynal DJ (2003) Use of nitrogen to phosphorus ratios in plant tissue as an indicator of nutrient limitation and nitrogen saturation. J Appl Ecol 40(3):523–534

Tian HQ, Chen GS, Zhang C, Charles ASH, Jerry MM (2010) Pattern and variation of C:N:P ratios in China’s soils: a synthesis of observational data. Biogeochemistry 98:139–151

Turrión MB, Glaser B, Solomon D, Ni A, Zech W (2000) Effects of deforestation on phosphorus poos in mountain soils of the Alay Range, Khyrgyzia. Biol Fertil Soils 31:134–142

Wang BK, Tang KL (1992) A study on temporal variation of reclamation on the Loess Plateau. J Soil Water Conserv 6(2):63–67 (in Chinese)

Wang KB, Shao RX, Shangguan ZP (2010) Changes in species richness and community productivity during succession on the Loess Plateau of China. Pol J Ecol 58(3):549–558

Wei XR, Shao MA, Fu XL, Horton R, Li Y, Zhang XC (2009) Distribution of soil organic C, N and P in three adjacent land use patterns in the northern Loess Plateau, China. Biogeochemistry 96:149–162

Yu Q, Chen QS, Elser JJ, He NP, Wu HH, Zhang GM, Wu JG, Bai YF, Han XG (2010) Linking stoichiometric homoeostasis with ecosystem structure, functioning and stability. Ecol Lett 13:1390–1399

Zeng QC, Liu Y, Fang Y, Rentian M, Rattan L, An SS, Huang YM (2017) Impact of vegetation restoration on plants and soil C:N:P stoichiometry on the Yunwu Mountain Reserve of China. Ecol Eng 109:92–100

Zhang YW, Shangguan ZP (2016) The coupling interaction of soil water and organic carbon storage in the long vegetation restoration on the Loess Plateau. Ecol Eng 91:574–581

Zhang YW, Shangguan ZP (2018) Interaction of soil water storage and stoichiometrical characteristics in the long-term natural vegetation restoration on the Loess Plateau. Ecol Eng 116:7–13

Zhang W, Zhao J, Pan F, Li D, Chen H, Wang K (2015) Changes in nitrogen and phosphorus limitation during secondary succession in a karst region in southwest China. Plant Soil 391:77–91

Zhang YW, Deng L, Yan WM, Shangguan ZP (2016) Interaction of soil water storage dynamics and long-term natural vegetation succession on the Loess Plateau, China. CATENA 137:52–60

Zhao FZ, Di K, Han XH, Yang GH, Feng YZ, Ren GX (2015) Soil stoichiometry and carbon storage in long-term afforestation soil affected by understory vegetation diversity. Ecol Eng 74:415–422

Zheng FL (2006) Effect of vegetation changes on soil erosion on the Loess Plateau. Pedosphere 16:420–427

Zheng FL, Tang KL, Zhang K, Cha X, Bai H (1997) Relationship ofeco-environmental change with natural erosion and artificially accelerated erosion. J China Geol 7(2):75–84 (in Chinese)

Zhou D, Zhao S, Zhu C (2012) The grain for green project induced land cover change in the Loess Plateau: a case study with Ansai County Shanxi Province, China. Ecol Indic 23:88–94