- Research article
- Open Access
Effects of tree species and topography on soil and microbial biomass stoichiometry in Funiu Mountain, China
BMC Ecology volume 20, Article number: 67 (2020)
Soil and microbial biomass stoichiometry plays an important role in understanding nutrient cycling in terrestrial ecosystems. However, studies on soil and microbial biomass stoichiometry in forests are rare. This study investigated the effect of tree species and topographic factors on the ecological stoichiometry of soil and soil microbial biomass.
Three types of forest stands (Quercus variabilis, Larix principis-ruprechtii, and Cotinus coggygria Scop.) in the Beiru River basin of Funiu Mountain were analyzed in September 2018. Six slope positions (sunny bottom slope, sunny middle slope, sunny top slope, shady bottom slope, shady middle slope, and shady top slope) were selected, and the total number of sampling plots was 108. The stoichiometric indices of soil and microbial biomass were determined.
At a depth of 0–10 cm, the soil organic C contents in different stands followed the order of C. coggygria (27.7 ± 5.2 g/kg) > Q. variabilis (24.5 ± 4.9 g/kg) > L. principis-ruprechtii (20.8 ± 4.3 g/kg) (P < 0.05). The soil organic C contents at depths of 0–10 cm with different slope aspects and at different slope positions also showed significant differences (P < 0.05). The highest MBC content was observed at the slope bottom (1002 ± 157 mg/kg), whereas the lowest was observed at the slope top (641 ± 98.3 mg/kg). Redundancy analysis showed that the contribution of tree species to these differences was 57.1%, whereas that of topographical factors was 36.2%.
Tree species more significantly affected soil nutrients and microbial biomass C, N and P than did topographic factors.
Soil microbes mineralize C, N and P and other nutrient elements in the soil through decomposition and release them into the soil, promoting the material cycle of the ecosystem [1, 2]. Soil microbial community structure and functional characteristics influence the biogeochemical cycle process [3, 4]. Ecological chemometry is a science that studies the balance of C, N, P and other elements in ecosystem interactions and processes . Cleveland and Liptzin  integrated the soil microbial biomass data of the global terrestrial ecosystem and thought that the proportions of C, N and P in the soil microbes was similar to the “Redfield ratio”. The soil microbes also display homeostasis. The soil microbial biomass N:P can be used as a tool to assess nutrient limitation [7, 8].
The C, N and P stoichiometric ratios of the soil microbial biomass can determine the nutritional status and restrictions on microbial growth [7, 9]. Ren et al.  studied a vegetation restoration area in the Loess Plateau, China, and found that nutrient limitation could be determined by soil microbial biomass N:P. In addition, plants display homeostasis, and a “Redfield-like ratio” is present in plants; as influenced by climate change and the biogeographic gradient of soil matrix age, the N:P ratio in plant leaves gradually increases from areas of high latitude to areas of low latitude . The above studies provide new ideas and research methods for the study of nutrient cycling and limitations in forest ecosystems.
Rapid forest vegetation growth consumes excessive soil nutrients [12,13,14]. Stand growth is gradually constrained by soil nutrients, especially by soil N and P. However, forest tree species, slope aspect and position and other conditions are not always consistent. The soil in the bottom slope and that in the shady slope often have better moisture and soil conditions, and the difference in soil nutrients further contributes to the variation in the N:P in tree leaves [15,16,17].
Soil microbial biomass is a critical factor in ecosystem nutrient cycling [16,17,18], and knowledge about its relation with the ecological stoichiometry of soil is crucial for us to gain insight into the nutrient cycling of terrestrial ecosystems [19,20,21,22]. However, there is currently little literature available regarding the variation in the ecological stoichiometry of soil microbes at different forest soil depths. Meanwhile, differences in soil nutrient content, bulk density and moisture content between different soil depths often cause variations in soil microbial biomass at different soil depths [23,24,25,26,27,28]. Little is understood about the effect of tree species and topography on the ecological stoichiometry of soil and soil microbial biomass.
Based on this information, this paper took Funiu Mountain, China, as the research site to study the effects of tree species and topography on soil and microbial biomass stoichiometry. The study aimed to test the following hypotheses: (1) Environmental factors (tree species, slope position and aspect) affect the C, N and P contents in the soil and microbial biomass of different soil depths (0–10 cm, 10–20 cm) and their ecological stoichiometry. (2) Forest soil microbes display homeostasis. (3) Nutrient limitations of the forest ecosystem in the study area can be assessed using the soil microbial biomass N:P.
Materials and methods
The study area, with an elevation of 780–1270 m, is located in the upper reaches of the Beiru River in Funiu Mountain, China. It has a warm temperate continental monsoon climate, with an average annual temperature of 12.4 °C, an annual sun exposure duration of 2103 h, a frost-free period of 198 days, and an average annual precipitation of 670 mm. The soil type is mainly brown soil, according to the soil classification criteria enacted by the General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China . The main tree species are secondary deciduous broadleaf trees, such as Quercus variabilis, Cotinus coggygria Scop. and Acer mono Maxim. In the 1950s, large areas of Larix principis-ruprechtii, Pinus tabuliformis Carr., and other coniferous trees were planted on Funiu Mountain, China. The herbaceous plants mainly include Miscanthus, Cyperus rotundus L., Thalictrum aquilegifolium, Setaria viridis (L.) Beauv., and Spiraea cantoniensis Lour., and the forest coverage percentage is 83.5% in the study area .
In September 2018, the three most representative forest stands in the investigated area, i.e., Q. variabilis, L. principis-ruprechtii and C. coggygria Scop., were selected. The stands were basically consistent in appearance, i.e., tree and crown densities. Three typical shady slopes and three sunny slopes were selected from each forest stand (the definitions for shady slopes and sunny slopes were based on the descriptions in GB 26424-2010-T Technical Specifications for Forest Resources Planning, Design and Investigation of China). Three typical non-catena plots (10 m × 10 m) were set up on each slope from bottom to top according to the elevation, i.e., bottom slope, middle slope and top slope. There were 108 (3 × 3×2 × 2×3) sample plots in total (including three replicates for each combination). In each sample plot, 5–6 sampling points were randomly selected, with a distance between sampling points of about 10 m. The floating leaves were gently removed, and soil samples were collected from depths of 0–10 and 10–20 cm with a soil auger whose bit length was 30 cm. The soil samples from the same depth of each sample plot were fully mixed. One kg of the above soil sample was taken back to the laboratory and screened through a 2-mm sieve to remove materials such as plant roots, stones, and litter. One part of fresh soil was preserved at 4 °C and the contents of soil microbial biomass C, N and P were measured within 10 days. The remaining soil samples were dried naturally and then screened through a 20-mesh sieve to obtain the 1-mm soil sample and then screened through a 0.149-mm sieve to measure the total nutrients of the soil. The specific information about the sampling points is shown in Table 1.
The procedures of this study were approved by the Environmental Protection Bureau of Lushan County, Henan.
The proposed method is based on the complexation of malachite green with phosphomolybdate under acidic conditions.
The total C and N in the soil were determined using a Vario Max C&N analyzer (Elemental Ltd., UK). The contents of C, N and P in soil microbial biomass were determined using chloroform fumigation extraction (CFE) [7, 31, 32]. Each soil sample was made up of subsamples marked “fumigated” and “unfumigated”. Before fumigation, soil samples were cultured at 4 °C for 7 days. The soil samples (equivalent to 3 g in dry weight) were fumigated with ethanol-free chloroform at 25 °C for 24 h, and then C and N were extracted using 0.5-M K2SO4. P was extracted using 0.5M NaHCO3 (pH 8.5) . The unfumigated soil was extracted in the same way, and the extract was filtered through a 0.45-μm syringe-driven filter. Soil organic carbon (SOC) and total nitrogen were measured using a total organic carbon analyzer (Shimadzu TOC-V CPH, Shimadzu Corp., Kyoto, Japan). The transforming factors of microbial biomass C and N are 0.45 and 0.54, respectively [31, 32, 34]. MBC and MBN represent the concentration differences between the fumigated and unfumigated samples . The extractant of P was transferred to a 96 microwell plate, and the content of P was measured by the ammonium molybdate-malachite green method . The transforming factor of microbial biomass P used in this paper was 0.40 [31, 35].
The mass ratio is often used in terrestrial ecosystem research , while the molar ratio is often used in aquatic ecology research; thus, the former was adopted for the stoichiometric ratios of C, N and P of soil and microbial biomass. Microsoft Excel 2010 was used for experimental data sorting and plotting. SPSS 20.0 was used for statistical analysis. Levene’s test was performed to determine the homogeneity of variance. For data with heterogeneous variance, logarithmic transformation was performed. For data with a normal distribution, one-way ANOVA was used for comparisons among groups (α = 0.05). Otherwise, nonparametric Kruskal–Wallis one-factor ANOVA was used. The p values were adjusted by the false discovery rate. Redundancy analysis was performed to determine the contributions of environmental factors to stoichiometry using CANOCO 5 to visually exhibit the relationship between variables and response variables . A general linear model (GLM) was used to assess the effects of tree species, slope aspect, position and their combined effect. Pearson correlation analysis was used to analyze the relationship between the contents of C, N and P in soil and soil microbial biomass and the corresponding C:N, C:P and N:P ratios.
SMATR2.0  was used for the standardized major axis (SMA) analysis of the relationship between the C, N and P contents in soil and those in soil microbial biomass [6, 31] to explore the homeostasis of the soil microbial biomass. SMA analysis presents the “best” binary fitting line between two variables. The relationship between the C, N, and P contents in soil and those in soil microbial biomass was expressed by the equation lgy = a+blgx, where a is the intercept and b is the slope. When the slope b is not significantly greater than 1, the relationship between the two variables is isometric .
Effects of tree species and topography on the C, N and P contents
At depths of 0–10 cm, the SOC contents of the soil of different tree species was in the order of C. coggygria Scop. (27.7 ± 9.8 g/kg) > Q. variabilis (24.5 ± 9.1 g/kg) > L. principis-ruprechtii (20.8 ± 9.8 g/kg) (p < 0.05; Fig. 1). There was no significant difference in the SOC content at 10–20 cm. The SOC content significantly differed (p < 0.05) at different soil depths of the three tree species. The effect of tree species on soil total nitrogen (TN) content was similar to that of SOC, and its effect on soil total phosphorus (TP) content was not significant. The variation in TP content was not significant in relation to tree species and soil depth.
At 0-10 cm, the SOC contents of Q. variabilis, L. principis-ruprechtii and C. coggygria Scop. were in the following order: bottom slope > middle slope > top slope, with the highest value in shady bottom slope and the lowest value in sunny top slope (Fig. 2). The difference between the shady/sunny bottom slope and sunny top slope was significant, and the SOC content in the shady slope was higher than that in the sunny slope overall. At 10–20 cm, the SOC content was greater in the bottom slope than in the top slope. However, there was no significant difference between different slope positions for the sunny slope of Q. variabilis and the shady/sunny slope of C. coggygria Scop. At 0–10 cm, the soil TN contents in the shady and sunny slopes of Q. variabilis, the sunny slope of L. principis-ruprechtii, and the shady slope of C. coggygria Scop. exhibited a similar tendency as that of SOC, that is, the highest content was found in the bottom slope (2.67 g/kg, 2.46 g/kg and 2.05 g/kg, respectively) and the lowest content was found in the top slope (1.6 g/kg, 1.33 g/kg and 1.74 g/kg, respectively), and the difference was significant. There was a significant difference in the soil TN contents between different slope positions for the shady slope of L. principis-ruprechtii and the sunny slope of C. coggygria Scop. However, the variation was not obvious. At 10–20 cm, the soil TN contents of the sunny slopes of Q. variabilis and L. principis-ruprechtii and the shady slope of C. coggygria Scop. were in the following order: bottom slope > middle slope > top slope, and the difference was significant. However, the other slopes did not show regular variation. Overall, the soil C, N and P contents at 1–10 cm were more noticeably influenced by tree species and topography than those at 10–20 cm.
At 0–10 cm, the TP contents of the sunny slopes of Q. variabilis and L. principis-ruprechtii and the shady slope of C. coggygria Scop. were in the following order: bottom slope > middle slope > top slope. However, the difference was not significant (P > 0.05). However, there was no such regular variation in other slopes except that there was a significant difference in the soil TP contents between the shady slope of Q. variabilis and that of L. principis-ruprechtii. At 10–20 cm, the soil TP contents in the sunny slopes of Q. variabilis and L. principis-ruprechtii and the shady slope of C. coggygria Scop. were in the following order: bottom slope > middle slope > top slope. There was a significant difference in the TP content between different slope positions for the sunny slope of Q. variabilis and that of L. principis-ruprechtii, whereas there was no significant difference in the TP content between different slope positions for the shady slope of C. coggygria Scop. There was no obvious regularity for the other slopes. However, there was a significant difference in the soil TP content.
As shown in Fig. 3, at 0–10 cm, the soil microbial biomass C content of different tree species was in the order of C. coggygria Scop. > Q. variabilis > L. principis-ruprechtii, with a significant difference between the three (p < 0.05). At 10–20 cm, there was no significant difference in the soil microbial biomass C contents of the three tree species (p > 0.05). At 0–10 cm, the soil microbial biomass C contents of the three tree species were significantly greater than those at 10–20 cm. The soil microbial biomass N content was similar to the C content. However, there were no significant differences in microbial biomass P between the two soil depths, among the three tree species, or in their interactions.
As shown in Fig. 4, at 0–10 cm, the soil microbial biomass carbon (MBC) contents of Q. variabilis, L. principis-ruprechtii and C. coggygria Scop. were invariably highest in the shady bottom slope (1002 ± 157 mg/kg) and lowest in the sunny top slope (641 ± 98.3 mg/kg). At 10–20 cm, the MBC content was invariably highest in the bottom slope and lowest in the top slope, with the exceptions of the shady slopes of L. principis-ruprechtii. Overall, the soil MBC content in the shady slope was higher than that in the sunny slope.
At 0–10 cm, the changes in the soil microbial biomass nitrogen (MBN) contents of the three tree species did not show regular patterns. The soil MBN content was highest in the shady bottom slope of C. coggygria Scop. and lowest in the sunny top slope of L. principis-ruprechtii. At 10–20 cm, the MBN content in the sunny slope of L. principis-ruprechtii and shady slope of C. coggygria Scop. was in the order of bottom slope > middle slope > top slope. There was no significant difference in the soil MBN between the shady and sunny slopes of the three tree species.
At 0–10 cm, the soil microbial biomass phosphorus (MBP) contents of Q. variabilis, L. principis-ruprechtii and C. coggygria Scop. did not show regular changes. Its contents in the sunny slope of Q. variabilis, sunny slope of L. principis-ruprechtii, and shady slope of C. coggygria Scop. were in the order of bottom slope > middle slope > top slope, and the difference was significant. The highest content of soil MBP was in the bottom slope of C. coggygria Scop (19.0 mg/kg), and the lowest was in the shady top slope of Q. variabilis (16.5 mg/kg). At 10–20 cm, no significant differences in the soil MBP contents were observed among the three slope positions of the sunny slopes of Q. variabilis, L. principis-ruprechtii and C. coggygria Scop.
Effects of tree species and topography on soil and microbial biomass stoichiometry
As shown in Tables 2 and 3, the C:N ratio in Q. variabilis soil was highest on the sunny bottom slope, and the soil C:N ratio was significantly different between the sunny bottom slope and the other slope positions. The C:N ratio in C. coggygria Scop. soil was highest in the sunny bottom slope, and there was no significant difference in the C:N ratio between the sunny bottom slope and the other slope positions. The C:N ratio in the L. principis-ruprechtii soil was highest in the shady top slope, and there was no difference between the different slope positions of L. principis-ruprechtii. The trends of soil C:P and N:P were basically consistent with that of soil C:N.
The soil microbial biomass C:N ratios of Q. variabilis and L. principis-ruprechtii were relatively high on the sunny slope. The soil microbial biomass C:N ratios of C. coggygria Scop. was highest in the shady top slope, and there was a significant difference in the soil microbial biomass C:N ratios between different slope positions. The soil microbial biomass C:P ratios of Q. variabilis and C. coggygria Scop. soil were highest in the sunny bottom slope, making them significantly different from those of the other slope positions. The soil microbial biomass C:P ratio of L. principis-ruprechtii was the highest in the shady top slope, and it was significantly different from that of the other slope positions. The variation trend of the soil microbial biomass N:P was consistent with that of the soil microbial biomass C:P.
Correlations among soil C, N, P, MBC, MBN, MBP, and their ecological stoichiometry
To explore the separate contributions of the environmental factors to soil and microbial biomass stoichiometry, RDA was performed, and the results are shown in Fig. 5. Tree species exhibited a significant impact on soil and microbial biomass stoichiometry (P = 0.002), with a contribution rate of 57.1%. The contribution rate of topographic factors, including slope aspect and position, was only 36.2%. Soil layers did not show a noticeable impact on C:N:P stoichiometry, whose contribution rate was only 6.7%.
As shown in Table 4, there was a strongly significant correlation between soil C, N and P and between soil MBC, MBN and MBP. Additionally, there was a strongly significant correlation between soil C, N and P and soil MBC, MBN and MBP. No significant correlation was found between soil C, N, P and MBC, MBN, and MBP with soil C:N. However, there was a strongly significant correlation between soil C:N and N:P. Additionally, soil N, P and MBC and N:P influenced soil MBC:MBN. Basically, there was no significant correlation between soil C, N, and P and soil MBC, MBN and MBP with MBC:MBP and MBN:MBP. In addition, no significant correlation was found between soil N:P and soil microbial biomass N:P.
Correlation between soil C, N, and P and soil MBC, MBN, and MBP contents and their ecological stoichiometry with environmental factors
Pearson correlation analysis also showed that at 0–10 cm, the soil C, N, and P contents and their ecological stoichiometry were closely related to tree species. The environmental factor analysis showed that the tree species factor could account for 80.5% of the variation, and the slope aspect and position and other topographic factors could account for 8.8% of the variation. The results obtained were basically consistent at both depths of 10–20 and 0–10 cm. Tree species was also the most important influencing factor, being responsible for more than 70% of the variation.
The results of SMA analysis of soil MBC, MBN, and MBP and soil C, N, and P indicated that the slopes between MBC and MBN, between MBC and MBP, and between MBN and MBP were all approximately equal to 1, showing a well-constrained proportion and indicating an isometric relation between soil MBC, MBN and MBP. The slopes of SOC, soil MBC, soil TP and soil MBP were significantly greater than 1 (Table 5). The relationships between soil C, N and P and soil MBC, MBN and MBP presented the characteristics of a nonisometric model, revealing a strong dependence of soil microbes on soil nutrient content.
Different tree species have different litter yields and physicochemical properties, and the amount of nutrient content returned to soil may also vary . The litter layer covering the soil surface can regulate microclimate conditions, such as water moisture and temperature, and impact the nutrient cycling rate of ecosystems [15, 41, 42]. In this study, tree species had a strongly significant impact on the soil C, N and P and soil MBC, MBN and MBP contents (Fig. 5). There is always a thick layer of litter over the soil surface of C. coggygria Scop., making it conducive to accumulating precipitation and reducing water evaporation from the soil mass; the litter layer of L. principis-ruprechtii is relatively thin but difficult to decompose; the differences in nutrient and moisture contents resulted in the difference in soil microbial biomass . For example, Devi and Yadava  showed that soil moisture content had a strong impact on soil microbial biomass, and tree species significantly affected soil and soil MBC, MBN, and MBP contents. Patel et al.  also thought that changes in land use significantly affected soil MBC, MBN and MBP contents. Zhu et al.  also concluded that different vegetation restorations on the Loess Plateau had significant impacts on soil C and N and on soil MBC and MBN.
Slope aspect, position and other topographic factors affect the transport, accumulation and distribution of nutrients. In the Funiu Mountain area, topographic factors affect soil microclimate conditions. Compared with the sunny slope, the shady slope has a relatively low rate of water evaporation and has relatively good moisture conditions, which in turn affect the community composition and species of plants [43,44,46]. In this study, at 0-10 cm, both slope aspect and position significantly influenced MBC, MBN and MBP contents. Soil MBC, MBN and other content indicators, which are sensitive to environmental change, can be used as indicators of soil quality change . The water content is relatively deficient in this study area. Soil MBC, MBN, and MBP are more sensitive to environmental change than are soil C, N, and P. The better water content condition in the shady slope is favorable to microbial growth and reproduction, facilitating nutrient transformation. The significant correlation between the soil C, N and P contents and the soil MBC, MBN, and MBP contents revealed that shady slopes, compared with sunny slopes, have higher soil C, N and P contents and soil MBC, MBN and MBP contents. The better water and nutrient conditions in the shady slope facilitate the growth of tree species. Moreover, the shady slope has a thicker litter layer (Table 1), which is particularly obvious for C. coggygria Scop. (Figure 1). At 10–20 cm, the impact of slope position on the soil C and N contents reached a very significant level, while the impact of slope position on the soil C and N contents reached a significant level. Slope position had a significant impact on soil MBC and MBN, while slope position did not. At 0–10 cm, aspect was a more important influencing factor. In contrast, at 10–20 cm, slope position was a more important factor in terms of its influence on the soil C, N and P and soil MBC, MBN and MBP contents. Zhu et al.  found that the effect of different tree species on soil organic C and N and soil MBC and MBN was greater than the effect of topographic factors (such as slope position) in the Loess Plateau, consistent with the results of this study. At the same time, the SMA analysis results indicate that the soil C, N and P contents and soil MBC, MBN, and MBP contents show a nonisometric relationship (Table 5), the slope between soil MBC and soil C is significantly greater than 1, and the slope between soil MBP and soil P content is even greater. A low availability of soil P and the different activation efficiency of soil P by soil microbes in different locations contribute to this result, indicating that soil microbes are strongly reliant on soil nutrient content, especially soil P, which may be attributed to the fact that P comes only from soil.
The change in plant species is a main factor affecting the soil C:N ratio. The difference in nutrient proportion and the nutrient cycling mechanism within plants impacts litter, which in turn affects soil C:N. At the same time, soil C:N usually reflects the decomposition rate of soil organic matter, and a low C:N ratio represents rapid decomposition of organic matter. In estimating the global soil carbon stock, the soil C:N ratio (10:1) is generally regarded as a constant; generally, the correlation analysis of soil C, N, and P and soil MBC, MBN and MBP with soil C:N suggests that there is no correlation between them, reflecting the stability of soil C:N on a microtopography scale [47, 48]. The results of this study were consistent with those reported in the literature [47, 48]: Slope aspect and position impact only C:P and N:P and have no impact on C:N. Soil ecological stoichiometry also vary according to soil and other microtopography factors .
Tree species also have a significant impact on soil MBC:MBN and MBC:MBP (Fig. 5), and their impact is greater at 0–10 cm than at 10–20 cm, which is consistent with the conclusions reached by Cleveland and Liptzin . Tree species, as a potential factor, influences soil microbial ecological stoichiometric ratios. Different tree species may lead to changes in litter quantity and quality or changes in the main composition of the microbial community. Heuck et al.  discovered that the ecological stoichiometry of soil microbial biomass varied in soil where P was deficient, which might be related to the change in the composition of the soil microbial community. Different microbial communities with different vegetation type conditions may cause changes in soil MBC:MBN and MBC:MBP. Tischer et al.  found that soil and microbial biomass stoichiometry varied according to language use and soil depths and that the soil C:N:P ratio noticeably decreased after forests were interference. However, our study found that tree species had an impact on soil MBC:MBN and MBC:MBP but no impact on soil MBN:MBP. Different from the impact of tree species on soil MBC:MBN and MBC:MBP, slope aspect and position had no impact on the ecological stoichiometry of soil microbial biomass, reflecting its stability on a microtopography scale. Only tree species and slope direction had a combined effect on the ecological stoichiometry of soil microbial biomass. These findings may be related to the unique climate conditions in the Huaihe River basin. Tree species, as the most important influencing factor, impacted soil MBC:MBN and MBC:MBP.
The homeostasis of organism stoichiometry, as a central concept in ecological stoichiometry, is used to reflect the maintenance and circulation of nutrients on scales ranging from subcellular to ecosystem. The homeostasis of ecological chemistry is defined as the extent to which an organism maintains the stability of its own chemical composition when the chemical composition and availability of external resources vary . The slope formed by the element content of organisms and the element content of resources after logarithmic transformation is called homeostasis. Sterner and Elser  defined a slope less than 1 as homeostasis, while McGroddy et al.  defined a slope equal to 1 as homeostasis. Cleveland and Liptzin  conducted an SMA analysis of forest soil and grassland soil microbial biomass and found that the slope of soil microbial biomass C-N, C-P and N-P was approximately 1 or less than 1 and that soil microbial biomass C was significantly related to N and P, indicating that soil microbes had homeostasis.
In this study, Pearson correlation analysis showed significant relationships among MBC, MBN and MBP. In addition, the SMA analysis showed that the slopes formed between soil MBC, MBN and MBP were smaller than—yet close to—1, with a well-constrained ratio. These findings were basically consistent with those reported by Cleveland and Liptzin , which indicated that soil microbes had homeostasis. Fanin et al.  studied the homeostasis of microbes in the decomposition process of plant litter and found that the slope formed between microbial N and P was significantly greater than 1, showing that they did not possess homeostasis. This result is not consistent with the results of this study and those of Cleveland and Liptzin , and this difference may be related to the environmental conditions in which microbes live. There is a relatively stable environment in the soil. When microbes are deficient in elements, they can be obtained from the soil. Unlike soil, a stable environmental condition is absent in the litter. The lack of elements necessary for the formation of microbes during leaf decomposition may cause microbes to rely heavily on N or P and result in a change in homeostasis.
Plant N:P usually varies with changes in tree species or vegetation [6, 51], mainly due to the biogeographical gradient difference in soil matrix age. In this study, the change in tree species did not have a significant impact on soil microbial biomass N:P. There was no relationship between the supply ratio of soil MBN:MBP and soil N:P, which was consistent with the results of Cleveland and Liptzin , indicating that soil microbial biomass N:P conforms to the “Redfield” ratio and that homeostasis controls the ratio of elements. In this study, slope aspect and position had little impact on the ecological stoichiometry of soil microbial biomass (Fig. 5). Correlations showed that the C, N and P contents of soil and microbial biomass and their ecological stoichiometry had no correlation with MBN:MBP. From the perspective of environmental factors, soil microbes exhibited homeostasis. N and P are generally deemed the most important limiting elements to plants in terrestrial ecosystems. The critical ratio of N:P in plant leaves is thought to be an indicator that can be used to assess the nutrient supply of the environment to plant growth [11, 37, 51,52,54]. When N:P is smaller than 14, N is limited; when N:P is greater than 16, P is limited; when N:P is between 14 and 16, the ecosystem is limited by both N and P or no limitations exist due to a sufficient nutrient supply. In areas with complex topography, the N:P ratio of plant leaves is usually used to indicate nutrient limitation in soil.
In tropical ecosystems, soil microbial biomass N:P could also be used to reflect nutrient limitation in soil, and a high soil microbial biomass N:P indicates the presence of P limitation . The N:P ratio of plant leaves has a very wide range . The range of measured leaf ratios included the above N:P breakpoint value ; thus, it was not possible to evaluate the N and P limitations at this site. However, the N:P ratio of soil microbial biomass was limited, and its range was much narrower than that of plant leaves. Therefore, the N:P ratio could better reflect the nutrient limitation of soil. Ren et al.  found that a high N:P ratio in soil microbes in the vegetation restoration area of the Loess Plateau demonstrated that plant growth was limited by P. In this study, the average soil microbial biomass N:P of tree species was 10.9 (equivalent to the atomic ratio of 24.1), which suggested that there was a relative lack of N in this area and that tree species may be mainly limited by N. This result was consistent with the conclusion reached by Bai et al. , i.e., forest region was restricted by N. This study also found that soil microbial biomass N:P was not sensitive to tree species variation, indicating that the homeostasis of soil microbial biomass N:P was better than that of plant leaves , which may be used to reflect the nutrient limitation of ecosystems. In this study area, the leaf N:P ratio of tree species could be combined with the soil microbial biomass N:P ratio to evaluate the nutrient limitations of ecosystems more accurately, providing a basis for ecosystem restoration and management.
Tree species and topography have significant impacts on the ecological stoichiometry of soil and soil microbial biomass, and soil microbes have homeostasis. Our research results are limited to only a small basin in Funiu Mountain, China. Only by a statistical analysis of a large set of sample data covering different climate zones and land uses across the world is it possible to reveal the impact mechanism of tree species and topography on the ecological stoichiometry of the soil and soil microbial biomass.
In conclusion, compared with topographic factors, tree species have a greater impact on soil nutrient content and microbial biomass. Under the same tree species conditions, aspect is the most important factor influencing soil nutrient content and soil microbial biomass at a soil depth of 0–10 cm, whereas slope position has a greater impact at 10-20 cm. There was a very significant correlation between the soil C, N and P contents and the soil MBC, MBN and MBP contents. SMA analysis also revealed that soil microbes heavily rely on soil nutrients, especially soil P, reflecting a close relationship between soil microbes and soil nutrient cycling. However, no correlation was found between the stoichiometric ratios of soil microbial biomass, soil C, N and P contents and their stoichiometric ratios with the slope aspect and position and other topographic factors. According to the SMA analysis of the soil MBC, MBN, and MBP, the soil microbial biomass in the investigated area has a well-constrained ratio, which suggests homeostasis of the soil microbes. The insensitivity of MBN:MBP to tree species and topographic variations and the homeostasis of the microbes may enable the ratio to be used to assess the nutrient limitations of ecosystems.
Availability of data and materials
The datasets analyzed during the current study are available from the corresponding author upon reasonable request.
Chloroform fumigation extraction
General linear model
Microbial biomass phosphorus
Heuck C, Weig A, Spohn M. Soil microbial biomass C: N: P stoichiometry and microbial use of organic phosphorus. Soil Biol Biochem. 2015;85:119–29.
Xia J, Niu S, Wan S. Response of ecosystem carbon exchange to warming and nitrogen addition during two hydrologically contrasting growing seasons in a temperate steppe. Glob Change Biol. 2009;15:544–1556.
An SS, Mentler A, Acosta-Martínez V, Blum WEH. Soil microbial parameters and stability of soil aggregate fractions under different grassland communities on the Loess Plateau, China. Biologia. 2009;64:424–7.
Devi NB, Yadava PS. Seasonal dynamics in soil microbial biomass C, N and P in a mixed-oak forest ecosystem of Manipur, Northeast India. Appl Soil Ecol. 2006;31:220–7.
Sterner RW, Elser JJ. Ecological Stoichiometry: the biology of elements from molecules to the biosphere. New Jersey: Princeton University Press; 2002.
Cleveland CC, Liptzin D. C: N: P stoichiometry in soil: Is there a “Redfield ratio”for the microbial biomass? Biogeochemistry. 2007;85:235–52.
Chen J, Seven J, Zilla T, Dippold MA, Blagodatskaya E, Kuzyakov Y. Microbial C:N: P stoichiometry and turnover depend on nutrients availability in soil: a 14C, 15N and 33P triple labelling study[J]. Soil Biol Biochem. 2019;131:206–16.
Chen XY, Daniell TJ, Neilson R, O’Flaherty V, Griffiths BS. Microbial and microfaunal communities in phosphorus limited, grazed grassland change composition but maintain homeostatic nutrient stoichiometry. Soil Biol Biochem. 2014;75:94–101.
Hartman WH, Richardson CJ. Differential nutrient limitation of soil microbial biomass and metabolic quotients(qCO2): is there a biological stoichiometry of soil microbes? PLoS ONE. 2013;8:e57127.
Ren CJ, Zhao FZ, Kang D, et al. Linkages of C: N: P stoichiometry and bacterial community in soil following afforestation of former farmland. For Ecol Manage. 2016;376:59–66.
Reich PB, Oleksyn J. Global patterns of plant leaf N and P in relation to temperature and latitude. Proc Natl Acad Sci. 2004;101:11001–6.
Shang S, Jiang P, Chang SX, Song Z, Liu J, Sun L. Soil organic carbon in particle size and density fractionations under four forest vegetation-land use types in subtropical China. Forests. 2014;5:1391–408.
Song YL, Wang KQ. Research progress of forestry ecosystem soil carbon storage accounting methods at home and aboard. J Green Sci Technol. 2018;9:1–6 (in Chinese).
Zhu H, Wu J, Guo S, Huang D, Zhu Q, Ge T, Lei T. Land use and topographic position control soil organic C and N accumulation in eroded hilly watershed of the Loess Plateau. CATENA. 2014;120:64–72.
Hu L, Ade L, Wu X, Zi H, Luo X, Wang C. Changes in Soil C:N: P Stoichiometry and microbial structure along soil depth in two forest soils. Forests. 2019;10:1–12.
Hu L, Xiang Z, Wang G, Rafique R, Liu W, Wang C. Changes in soil physicochemical and microbial properties along elevation gradients in two forest soils. Scand J For Res. 2016;31:242–53.
Tomohiro Y, Fujio H, Tokuchi N. Seasonal effects on microbial community structure and nitrogen dynamics in temperate forest soil. Forests. 2018;9:153–70.
Hall EK, Maixner F, Franklin O, Daims H, Richter A, Battin T. Linking microbial and ecosystem ecology using ecological stoichiometry: a synthesis of conceptual and empirical approaches. Ecosystems. 2011;14:261–73.
Bergkemper F, Welzl G, Lang F, Krüger J, Schloter M, Schulz S. The importance of C, N and P as driver for bacterial community structure in German beech dominated forest soils. J Plant Nutr Soil Sci. 2016;179:472–80.
Fanin N, Bertrand I. Aboveground litter quality is a better predictor than belowground microbial communities when estimating carbon mineralization along a land-use gradient. Soil Biol Biochem. 2016;94:48–60.
Heuck C, Spohn M. Carbon, nitrogen and phosphorus net mineralization in organic horizons of temperate forests: stoichiometry and relations to organic matter quality. Biogeochemistry. 2016;131:229–42.
Xu X, Thornton PE, Post WM. A global analysis of soil microbial biomass carbon, nitrogen and phosphorus in terrestrial ecosystems. Glob Ecol Biogeogr. 2013;22:737–49.
Baldrian P, Merhautová V, Petránková M, Cajthaml T, Šnajdr J. Distribution of microbial biomass and activity of extracellular enzymes in a hardwood forest soil reflect soil moisture content. Appl Soil Ecol. 2010;46:177–82.
Brockett BFT, Prescott CE, Grayston SJ. Soil moisture is the major factor influencing microbial community structure and enzyme activities across seven biogeoclimatic zones in western Canada. Soil Biol Biochem. 2012;44:9–20.
Florinsky IV, McMahon S, Burton DL. Topographic control of soil microbial activity: A case study of denitrifies. Geoderma. 2004;119:33–53.
Kirkby CA, Richardson AE, Wade LJ, Batten BGD, Blanchard C, Kirkegaard JA. Carbon-nutrient stoichiometry to increase soil carbon sequestration. Soil Biol Biochem. 2013;60:77–86.
Strickland MS, Rousk J. Considering fungal: bacterial dominance in soils-Methods, controls, and ecosystem implications. Soil Biol Biochem. 2010;42:1385–95.
Tian YW, He CL, Xu SJ, Wang N. Spatial distribution and dynamic changes in soil organic carbon density in a restored grassland converted from farmland. Fresenius Environ Bull. 2019;28:3153–61.
General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China. Soil classification and codes of China in 2019 GBT 17296–2009. Beijing: China Standards Press; 2009.
Tian YW, Li DS, Wang N, Liu J, Wang Y. Changes of runoff, sediment, nitrogen and phosphorus exported from Taowan Watershed, Funiu Mountainous Area. J Northwest For Univ. 2016;31(5):142–7.
Brookes PC, Landman A, Pruden G, Jenkinson DS. Chloroform fumigation and the release of soil nitrogen: a rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biol Biochem. 1985;17:837–42.
Wu J, Joergensen RG, Pommerening B, Chaussod R, Brookes PC. Measurement of soil microbial biomass C by fumigation-extraction: an automated procedure. Soil Biol Biochem. 1990;22:1167–9.
Olsen SR, Cole CV, Wantabe FS, Dean LA. Estimation of available phosphorus in soils by extraction with sodium bicarbonate. USDA Circular No. 939. Washington: US Government Printing Office; 1954.
Joergensen RG, Mueller T. The fumigation-extraction method to estimate soil microbial biomass: calibration of the kEN value. Soil Biol Biochem. 1996;28:33–7.
Bilyera N, Blagodatskaya E, Yevdokimov I, Kuzyakov Y. Towards a conversion factor for soil microbial phosphorus. Eur J Soil Biol. 2018;87:1–8.
D’Angelo E, Crutchfield J, Vandiviere M. Rapid, sensitive, microscale determination of phosphate in water and soil. J Environ Qual. 2001;30:2206–9.
Güsewell S. N: P ratios in terrestrial plants: variation and functional significance. New Phytol. 2004;164:243–66.
Lee SW, Hwang SJ, Lee SB, Hwang HS, Sung HC. Landscape ecological approach to the relationships of land use patterns in watersheds to waterquality characteristics. Landsc Urban Plan. 2009;92:80–9.
Falster DS, Warton DI, Wright IJ (2006) User’s Guide to SMATR: Standardised Major Axis Tests & Routines Version 2.0 (2006-07-02)[2017-06-14]. http://www.biomqeduau/ecology/SMATR2.0. Accessed 14 June 2017
Sigua GC, Coleman S, Albano J, et al. Spatial distribution of soil phosphorus and herbage mass in beef cattle pastures: effects of slope aspect and slope position. Nutr Cycl Agroecosyst. 2010;89:59–70.
Margesin R, Jud M, Tscherko D, Schinner F. Microbial communities and activities in alpine and subalpine soils. FEMS Microbiol Ecol. 2008;67:208–18.
Sayer EJ. Using experimental manipulation to assess the roles of leaf litter in the functioning of forest ecosystems. Biol Rev. 2006;81:1–31.
Patel K, Kumar NJI, Kumar NR, et al. Seasonal and temporal variation in soil microbial biomass C, N and P in different types land uses of dry deciduous forest ecosystem of Udaipur, Rajasthan, Western India. Appl Ecol Environ Res. 2010;8:377–90.
Jobbágy EG, Jackson RB. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol Appl. 2000;10:423–36.
Máren IE, Karki S, Prajapati C, et al. Facing north or south: Does slope aspect impact forest stand characteristics and soil properties in a semiarid trans-Himalayan valley? J Arid Environ. 2015;121:112–23.
Ružek L, Voríšek K, Strnadova S, Nováková M, Barabasz W. Microbial characteristics, carbon and nitrogen content in cambisols and luvisols. Plant Soil Environ. 2004;50:196–204.
Stevenson FJ, Cole MA. Cycles of soils: carbon, nitrogen, phosphorus, sulfur, micronutrients. Hoboken: Wiley; 1999.
Tian HQ, Chen GS, Zhang C, Melillo JM, Hall CAS. Pattern and variation of C:N: P ratios in China’s soils: a synthesis of observational data. Biogeochemistry. 2010;98:139–51.
Tischer A, Potthast K, Hamer U. Land-use and soil depth affect resource and microbial stoichiometry in a tropical mountain rainforest region of southern Ecuador. Oecologia. 2014;175:375–93.
McGroddy ME, Daufresne T, Hedin LO. Scaling of C:N: P stoichiometry in forests worldwide: Implications of terrestrial red field-type ratios. Ecology. 2004;85:2390–401.
Fanin N, Fromin N, Buatois B, et al. An experimental test of the hypothesis of non-homeostatic consumer stoichiometry in a plant litter-microbe system. Ecol Lett. 2013;16:764–72.
Aerts R, Chapin FS III. The mineral nutrition of wild plants revisited: a re-evaluation of processes and patterns. Adv Ecol Res. 2000;30:1–67.
Koerselman W, Meuleman AFM. The vegetation N: Pratio: a new tool to detect the nature of nutrient limitation. J Appl Ecol. 1996;33:1441–50.
Tessier JT, Raynal DJ. Use of nitrogen to phosphorus ratios in plant tissue as an indicator of nutrient limitation and nitrogen saturation. J Appl Ecol. 2003;40:523–34.
Townsend AR, Cleveland CC, Asner GP, Bustamante MMC. Controls of foliar N: P ratios in tropical rain forests. Ecology. 2007;88:107–18.
Bai XJ, Zeng QC, An SS, et al. Ecological stoichiometry characteristics of leaf-litter-soil in different plantations on the Loess Plateau, China. Chin J Appl Ecol. 2016;27:3823–30 (in Chinese).
This research was funded by the National Natural Science Foundation of China (Grant Number 31670616). The funder had no role the design of the study, the collection, analysis and interpretation of data, or in writing the manuscript.
Ethics approval and consent to participate
The procedures of this study were approved by the Environmental Protection Bureau of Lushan County, Henan.
Consent to publish
The authors have no competing interests to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Tian, Y., Qiao, D., Xu, S. et al. Effects of tree species and topography on soil and microbial biomass stoichiometry in Funiu Mountain, China. BMC Ecol 20, 67 (2020). https://doi.org/10.1186/s12898-020-00332-4
- Ecological stoichiometry
- Soil microbes
- Tree species