Abstract
Forest soils contain a large amount of organic carbon and contribute to terrestrial carbon sequestration. However, we still have a poor understanding of what nutrients limit soil microbial metabolism that drives soil carbon release across the range of boreal to tropical forests. Here we used ecoenzymatic stoichiometry methods to investigate the patterns of microbial nutrient limitations within soil profiles (organic, eluvial and parent material horizons) across 181 forest sites throughout China. Results show that, in 80% of these forests, soil microbes were limited by phosphorus availability. Microbial phosphorus limitation increased with soil depth and from boreal to tropical forests as ecosystems become wetter, warmer, more productive, and is affected by anthropogenic nitrogen deposition. We also observed an unexpected shift in the latitudinal pattern of microbial phosphorus limitation with the lowest phosphorus limitation in the warm temperate zone (41-42°N). Our study highlights the importance of soil phosphorus limitation to restoring forests and predicting their carbon sinks.
Similar content being viewed by others
Introduction
Forests have been proposed as a natural-solution to enhance carbon (C) sequestration in terrestrial ecosystems1,2 because they cover ~4200 million hectares of the world’s land surface3, and 50–90% of the total annual C flux of terrestrial ecosystems occurs at the forest-atmosphere interface4. However, primary productivity and other biological processes in forest ecosystems are widely constrained by soil nutrient availabilities especially by soil nitrogen (N) and phosphorus (P)5,6,7,8. Microbial metabolism serves as the base of detrital food-webs, mediates soil nutrient cycles and releases soil carbon (C) to atmosphere9. Biogeochemical constraints on catabolic processes, relative to primary production, determine the magnitude of soil respiration and the retention of detrital organic C in soils9,10. Therefore, nutrient limitation to microbial metabolism is a potential control in soil C turnover and release in forest ecosystems9,11. It is especially relevant to quantify the potential of forests to capture C in their soils as restoration efforts are expected to plant millions of trees over the next few decades.
Previous studies suggest that microbial metabolism is commonly limited by P in tropical forests12,13, while limited by N in temperate forests14,15. However, recent studies suggest that soil microbial metabolism can be limited by P in non-tropical zones16,17, and by N in subtropical soils6. In general, forest N and P availability vary with distance from the equator, plant productivity, climate, and soil depth7,18,19. Accurately predicting the role of soil nutrients in controlling microbial activity must account for all these environmental factors simultaneously. For example, C inputs and nutrient availability are generally lower in subsoil than topsoil so that limiting factors for microbial metabolism often vary with depth within a single profile17,19. Therefore, exploring the patterns and drivers of microbial metabolic limitation in soil profiles across large geographical scales could help decipher the effects of both local and geographical controls.
The forested areas of China account for 5% of the world’s forests and currently store over 4.75 Pg of C20. They are experiencing rapid revegetation in recent decades and their soils are considered significant C sinks with large C sequestration potential21,22. Although several studies have estimated pervasive soil P limitation to aboveground plant production in this region7,8, comprehensive estimates of belowground microbial nutrient limitations are unavailable. Ecoenzymatic stoichiometry based on the activities of extracellular enzymes (ecoenzymes), the proximate agents of microbial metabolism, connects the metabolic theory of ecology and ecological stoichiometry theory to predict microbial nutrient assimilation and growth9, and thus is an appropriate method to predict microbial metabolic limitation at a community level in soils17,23,24,25,26.
Here, we used ecoenzymatic stoichiometry methods to investigate patterns of microbial metabolic limitations in contrasting forests across a wide latitudinal gradient, (from tropical to cold regions) and at different soil depths (organic (O), eluvial (A), and parent material (C) layers) in China. Our results revealed widespread soil P limitation to microbial metabolism across Chinese forest ecosystems with the lowest limitation in warm-temperate regions. Significant correlations between microbial and plant P limitations over a wide geographic scale indicate synchronous P constraints in both microbial metabolism and primary productivity. These results imply that soil P availability is a long-term and inherent constraint on the turnover of soil organic C and the potential of ecosystem C sink in Chinese forests.
Results and discussion
Patterns of microbial metabolic limitation in soil profiles and latitudes across Chinese forests
Soil samples from 181 sites of Chinese forests showed microbial P limitation (microbial N/P limitation value > 0) at most locations (Fig. 1a), from tropical to boreal forests, with significant variability across soil horizons, vegetation types and climate zones (P < 0.001; Fig. 1, Supplementary Tables 4 and 5). Although 31% of all sampling sites from warm-temperate to cold-temperate zones were N limited (microbial N/P limitation value <0) (Supplementary Table 6), P limitation was pervasive across these Chinese forests and had much higher magnitude than N limitation. We thus focused on microbial P rather than N limitations in subsequent discussions.
Within soil profiles, microbial P limitation significantly increased with soil depth across climate zones and vegetation types, with the O-layer having the lowest microbial P limitation (P < 0.001; Fig. 1a, f, g), consistent with a recent study by Jing et al.17. In contrast, relative C limitation was significantly higher in the surface soils than in the subsoils for all vegetation types and climate zones (P < 0.001; Fig. 1a, d, e). Because relative C limitation indicates the hydrolysis potential of organic C rather than the metric of absolute microbial C limitation25, this result suggests that microorganisms in surface soils have stronger potential to drive organic C decomposition than in subsoils. This profile pattern of relative C limitation can be attributed mainly to higher organic C content in surface than subsoil layers providing more organic C to support microbial metabolism (Supplementary Fig. 4).
Geographically, piecewise linear regression analyses identified significant latitudinal breakpoints at about latitude 27°N and 41°N for relative C limitation and N/P limitation, respectively (P < 0.001; Fig. 2a, b), although these breakpoints have obvious variations among soil layers (Fig. 2c). Relative C limitation significantly increased from 19°N to 27°N and then decreased from 27°N to 54°N, suggesting the highest potential for organic C decomposition in the subtropical zone (latitude 24–30°N) (Fig. 1a). An opposite latitude pattern in N/P limitation indicated higher microbial P limitation in tropical and cold-temperature zones, with the lowest microbial P limitation in the warm-temperate zone (latitude 41–42°N) (Fig. 1b). These results suggested that soil P limitation to microbial metabolism remains unexpectedly high at relatively high latitudes in Chinese forests. It provides novel evidence that P limitation to soil microorganisms might be far more important than expected for boreal forests, compared to forests in other regions of the world usually thought to be limited by N27,28.
Drivers of microbial P limitation in Chinese forests
Increased microbial P limitation with soil depth in our study seems to contradict the prevailing idea of soil parent material being the main source of P29,30. However, two mechanisms could generate this vertical pattern of microbial P limitation in our study. First, soil organic matter is a potential pool of available P9,31. Higher organic matter decomposition reflected by higher relative C limitation in the O-layer than in A- and C-layers (Fig. 1a, d, e) suggests higher P release from organic matter in the O-layer compared to A- and C-layers. Thus, organic matter decomposition could help alleviate microbial P limitation in the O-layer in a P-deficient environment. Second, differential uptake of P by root systems with depth could be another mechanism driving vertical patterns of microbial P limitation. Vegetation index showed positive and increasing effects on microbial P limitation from O- to C-layers (Fig. 3d–j), suggesting an increased competition in P-acquisition between plants and microbes with increasing roots densities from O- to C-layers. A consistent pattern of P limitation between plants and microbes across this large geographical scale and soil layers (Fig. 4a–c) supports this argument and suggests stronger nutrient competition of plants with microbes in subsoils with abundant root systems than in topsoils32,33. Our previous study on rhizospheres also found that plants and microbes show strong competition for P when microbes are P constrained34. Therefore, P competition between roots and microbes could contribute to shaping these vertical patterns of microbial P limitation.
Overall, we found that the combined effects of climate, vegetation and soil factors played the most important roles in latitudinal patterns of microbial P limitation across soil profiles (Fig. 3a–c). Random-forest models and partial least squares path modeling further found that vegetation index driven by climate conditions (including aridity index and mean annual precipitation) has the largest positive effects on microbial P limitation (Fig. 3d–g). These results suggested that primary productivity has a central role in contributing to the geographic pattern of microbial P limitation. Many processes relevant to plant production can limit the availability of soil P in forests including biomass accumulation and its effect on soil pH6,30. The most direct process is that the fixation of P in biomass by plant growth reduces the availability of soil P, which is supported by the consistent pattern we found between P limitation in plants and soil microbial communities (Fig. 4a–c).
Highly productive ecosystems often develop on acidic soils, which can further decrease P availability for plants and microbes. This is especially important in warm, wet ecosystems in the tropics12. However, high P limitation at high latitudes can result from cold and dry conditions that suppress microbial activities and decomposition19,35 (Fig. 2a), subsequently slowing P release from organic matter9. In addition, P limitation should be exacerbated by anthropogenic N deposition (Fig. 3 and Supplementary Figs. 1–3), because it can create an imbalance in soil N:P stoichiometry and increase soil acidification36. Indeed, some low latitude forests around the world, especially in China, are receiving high N deposition due to geographic distributions of industry and human populations37,38,39, which implies an even greater increase in microbial P limitation at lower latitudes in the future.
Implications and uncertainties
Our recent incubation experiments found that P addition significantly alleviated microbial P limitation, and resulted in an 11.5–29.1% increase in C emissions from organic C decomposition40. This suggested that constraining the decomposition of organic matter under P limitations can potentially decrease soil C loss. In plant-microbe systems, however, soil P limitation constrains not only microbial metabolism but also plant growth18,41. In contrast, soils with sufficient P supply could maintain both high primary productivity and high soil C emissions via microbial decomposition. This argument is also supported in part by the counterbalance hypothesis of plant-microbe interaction that purports a homeostatic mechanism stabilizing microbial metabolism in relation to gross primary production wherein the greatest release of C by soil microorganisms occurs in ecosystems having the highest levels of plant production41. Therefore, soil P supply may be a key driver of C turnover in forest ecosystems, i.e., low turnover rate under P limitation and higher rate with sufficient P, whereas the potential soil C sink depends on relative availabilities of soil P to plants versus microbes.
The consistent effects of soil P limitation on both microbial metabolism and primary productivity in our study provide cues to explaining the “soil C saturation” hypothesis (i.e., soils have a limited capacity to store C as organic inputs increase)42 (Fig. 4a–c). Recent studies found that tree plantings in organic soils did not increase net C sequestration after 12 to 39 years43, and CO2 enrichment did not increase C sequestration at the ecosystem level in mature forests44. An inherent constraint suggested by our study is that soil P availability simultaneously limits both plant growth and microbial catabolism, thus constraining both ecosystem C balance and soil C sink. Consequently, forest ecosystems may not sequester as much C from vegetation restoration or afforestation as expected due to inherent P constraints.
Although our study reveals geographical and vertical profile patterns of microbial metabolic limitations in the forest soils of China, the present study has uncertainties. First, the ecoenzymatic stoichiometry method we used is based on the activities of relatively few enzymes to represent microbial nutrient acquisition, and is bound to include uncertainties as suggested by several recent studies45,46. Second, we used plant nutrient limitation data based on geographic coordinates that corresponded to the microbial nutrient limitations derived in our study to explore their spatial couplings, but we did not account for likely temporal and spatial differences in both sets of data. The lack of precise space-time pairing of sites limits our knowledge of how nutrient limitation affects plant-microbial interaction on micro-scales such as in the rhizosphere32. Third, global change and human activity are increasing uncertainties in the coupling of soil nutrient constraints to ecosystem C cycles43,47. For example, different adaptabilities of plants and microbes to environmental changes41,47 means that N deposition may affect plant growth and microbial catabolism in different ways and magnitudes through changes soil nutrient limitations. This may result in a dynamic non-equilibrium state of the C cycle in microbial-plant systems, resulting in a large C source or sink effect within any terrestrial ecosystem. Regardless of these uncertainties, our study is of direct relevance to current prediction of forest C sinks, and has significant implications for future study and policy decisions in forest restoration.
Materials and methods
Study area and soil sampling
We sampled almost all accessible mountains having Chinese forest ecosystems (Fig. 1a and Supplementary Data), which span a wide geographic range (18.9°–53.5°N, 101.0°–129.7°E, 210-3830 m above sea level). These forests are found in national parks or nature reserves with limited human interference, and include broadleaf forests (BF), coniferous-broadleaf mixed forests (CBF) and coniferous forests (CF). Climates include tropical, subtropical, warm-temperate, temperate, and cold-temperate zones. The geographical coordinates of each sample site were recorded using a global positioning system (GPS) device (eTrex Venture, USA). Mean annual precipitation (MAP) and mean annual temperature (MAT) ranged from −5.91 °C to 23.18 °C and 242 mm to 2667 mm, respectively.
Field sampling was conducted from July 2012 to March 2013. At each selected forest, one to fifteen sampling sites were located along an altitudinal transect based on the range of community types available (Supplementary Data). At each site, three 10 m × 10 m plots located 50 m away from each other were established. Five samples of each O, A and C-horizon were collected from each plot and pooled as one composite sample. In total, 504 soil samples were collected from 181 sites (some sites had only 1 or 2 soil horizons). After removing roots, litter, debris, and stones, composite soil samples were divided into two subsamples. One subsample was stored at 4 °C for later analysis of enzymatic activity. Because extracellular enzymes (i.e., ecoenzymes) are the products of long-term rather than instantaneous metabolism of microorganisms and they are relatively stable in soils9, their seasonal variability could be much smaller than the differences between environmental gradients. Thus, the enzymatic activities from our study would be representative of the large-scale sampling across Chinese forests. The other subsample was passed through a 2 mm sieve and air-dried for physicochemical analysis. Soil bulk density (BD) was measured using a soil bulk sampler and a stainless-steel cutting ring (5 cm in diameter by 5 cm long), adjacent to the soil sampling points in each plot.
Environmental factors
Our study included 14 environmental variables (climate, vegetation, and soil), obtained either from field investigations or from satellites and databases. MAT, MAP and both temperature and precipitation seasonality (TSEA and PSEA) were obtained from the Worldclim database (https://www.worldclim.org; ~1 km resolution). The aridity index of each site was estimated by the ratio of MAP to potential evapotranspiration (PET) (aridity index = MAP/PET). The MAT, MAP and aridity index values ranged from 5.91 °C to 23.18 °C, 242 mm to 2667 mm, and 0.24 to 2.42, respectively. We obtained bulk N deposition rates and acid rain data for each region in China from 1980 to 2015 from Yu et al.39 and the National Earth System Science Data Center (http://www.geodata.cn), respectively. Vegetation index (NDVI, 2001-2015) was obtained from the Advanced Very High Resolution Radiometer (AVHRR) developed by the Global Inventory Modeling and Mapping Studies (GIMMS) group (https://ecocast.arc.nasa.gov/data/pub/gimms/3g.v1/) at a 1/12° resolution.
Soil properties of samples (soil pH, soil organic C (SOC), total N and total P) were measured for each horizon (O, A, C) using standardized protocols according to methods described in a previous study34. Each sample was measured in triplicate.
Measurements of microbial metabolic parameters
The activities of two extracellular C-acquiring enzymes (β−1,4-glucosidase (BG) and β-d-cellobiosidase (CBH)), two extracellular N-acquiring enzymes (β−1,4-N-acetylglucosaminidase (NAG) and L-leucine aminopeptidase (LAP)), and one extracellular organic P-acquiring enzyme (alkaline or acid phosphatase (AP)) (Supplementary Table 1) were determined using the method of Saiya-Cork et al.48 and German et al.49. Enzyme activities were expressed as nanomoles of substrate released per hour per gram of dry soil (nmol g−1 h−1). The experimental procedures are described in our previous study16.
Calculating microbial metabolic limitations
Microbial metabolic limitations were quantified by calculating the vector lengths and angles of enzymatic activity based on untransformed proportional activities (e.g. [BG + CBH]/[BG + CBH + NAG + LAP]). Vector length, representing relative C limitation of microbial metabolism, was calculated as the square root of the sum of x2 and y2 (Eq. 1), where x is the relative activity of C vs. P-acquiring enzymes, and y is the relative activity of C vs. N-acquiring enzymes25. Vector angle, representing microbial N/P limitation, was calculated as the arctangent of the line extending from the plot origin to point (x, y). The relative C limitation to microbes increases with vector length whereas a vector angle of 45° roughly represents the boundary between microbial P limitation (> 45°) and N limitation (< 45°). This definition is based on the global ecoenzymatic stoichiometry law for C:N:P activities near 1:1:1 at a global scale, i.e., the slopes of their regression relationships are close to 45°23,24. In the present study, the slopes of the relationships between enzyme activities for C vs. P and C vs. N acquisition were very close to 1.0 for all sites (Fig. 1b, c), indicating that the ecoenzymatic stoichiometry method and the 45° boundary of N/P limitation using in Chinese forests is reasonable. The relative strength of microbial P limitation increases with vector angle and strength of microbial N limitation increases as the angle decreases.
Considering potential competition of microbes with plants for soil N and P33, we further linked our predictions of microbial metabolic limitation with the plant nutrient limitations from Du et al.8 to evaluate the relationships of microbial N/P limitation with plant N/P limitation. We extracted the ln-transformed NREdom/PREdom (ln NREdom/PREdom) from the database of Du et al.8 on the same or adjacent sites to our sampling sites. In total, values of ln NREdom/PREdom from 55 matched sites were obtained with a wide latitudinal range across Chinese Forests (21.9°–51.6°N). Du et al.8 defines ln NREdom/PREdom > 0 indicating plant N limitation and the value < 0 indicating plant P limitation. To facilitate interpretations in the present study, we redefined the plant nutrient limitation as follows: plant N/P limitation = −(ln NREdom/PREdom). Thus, in our study, the value of plant N/P limitation > 0 indicates plant P limitation, while the value < 0 indicates plant N limitation. As a result, consistent with soil microbial P and N limitations, the strength of plant P limitation increases with value of plant N/P limitation and strength of plant N limitation increases as the value decreases.
Statistical analysis
We first tested the variable normality and homogeneity of variances for microbial metabolic limitations, and found that neither relative C limitation nor N/P limitation conformed to normal distributions or homogeneity of variances under different vegetation types and climate zones (Supplementary Tables 2 and 3). Thus, we adopted the linear mixed-effects models to identify the effects of soil layer, forest type and climate type, as well as their interaction effects, on microbial metabolic limitations. The model was constructed using the “nlme” package in R with “soil layer”, “forest type” and “climate type” as the fixed factors and “sampling site” as a random factor50. The standardized major axis (type II) regression via the “lmodel2” package was used to identify relationships between extracellular enzymatic activities51 (Supplementary Note 1).
For examining latitudinal patterns of microbial metabolic limitations, we first identified the latitudinal breakpoints from a piecewise linear regression analysis. The regression relationships of the relative C and N/P limitations with latitude were fitted and tested with linear models using the “segmented” package52. The confidence intervals of the breakpoints were calculated via 1000 bootstrap samples using the “SiZer” package53. We then used general linear models to fit the relationships of the relative C and N/P limitations with upper and lower latitudes (before and after breakpoints), respectively. Generalized linear models also examined relationships of microbial N/P limitation with plant N/P limitation. Pearson correlations examined relationships of microbial metabolic limitations with potential predictor variables.
Variation partitioning analysis was conducted using the “vegan” package54 to identify the relative importance of climate, vegetation and soil properties to variations in microbial N/P limitation in different soil layers. We further determined the relative influence of each variable using the random-forest algorithm using the “RandomForest” package55. In these random-forest models, fourteen variables served as predictors for the microbial N/P limitation (Supplementary Table 7). To estimate the importance of these variables, we used percentage increases in the MSE (mean squared error) of variables, and higher MSE values indicate more important variables. Significance of the models and cross-validated R2 values were assessed with 1500 permutations of the response variable, by using the “A3” package. Similarly, the significance of each predictor to the response variables was assessed with the “rfPermute” package (Supplementary Note 2).
Furthermore, to explore cascading relationships between microbial N/P limitation and these key drivers, all significant variables (drivers) were divided into seven categories (hydrothermal condition (aridity index and MAT), N deposition, climate seasonality (PSEA and TSEA), vegetation, soil pH, SOC and TP). Partial least squares path models were used to identify possible pathways whereby variables control microbial N/P limitation in different soil layers. The VIF of all variables in the model requires <10 for removing the internal factor collinearity of the block (module). Meanwhile, loading in Outer Models requires >0.7 for removing those independent variables with small contributions. The models were constructed using the “plspm” package56. All statistical analyses were performed using the R software (v.3.3.2)57.
Reporting summary
Further information on research design is available in the Nature Research Reporting Summary linked to this article.
Data availability
All data related to the prediction of soil microbial metabolic limitation in this study are available at Figshare (https://doi.org/10.6084/m9.figshare.20327697.v1).
References
Bastin, J. F. et al. The global tree restoration potential. Science 364, 76–79 (2019).
Lewis, S. L., Wheeler, C. E., Mitchard, E. T. A. & Koch, A. Restoring natural forests is the best way to remove atmospheric carbon. Nature 568, 25–28 (2019).
Bonan, G. B. Forests and climate change: forcings, feedbacks, and the climate benefits of forests. Science 320, 1444–1449 (2008).
Beer, C. et al. Terrestrial gross carbon dioxide uptake: global distribution and covariation with climate. Science 329, 834–838 (2010).
Vitousek, P. M., Porder, S., Houlton, B. Z. & Chadwick, O. A. Terrestrial phosphorus limitation: mechanisms, implications, and nitrogen-phosphorus interactions. Ecol. Appl. 20, 5–15 (2010).
Camenzind, T., Httenschwiler, S., Treseder, K. K., Lehmann, A. & Rillig, M. C. Nutrient limitation of soil microbial processes in tropical forests. Ecol. Monogr. 88, 4–21 (2018).
Hou, E., Luo, Y., Kuang, Y., Chen, C. & Wen, D. Global meta-analysis shows pervasive phosphorus limitation of aboveground plant production in natural terrestrial ecosystems. Nat. Commun. 11, 1–9 (2020).
Du, E. et al. Global patterns of terrestrial nitrogen and phosphorus limitation. Nat. Geosci. 13, 221–226 (2020).
Sinsabaugh, R. L. & Follstad Shah, J. J. Ecoenzymatic stoichiometry and ecological theory. Annu. Rev. Ecol. Evol. Syst. 43, 313–343 (2012).
Houghton, R. A. Balancing the global carbon budget. Annu. Rev. Earth Planet. Sci. 35, 313–347 (2007).
Chen, J. et al. Differential responses of carbon-degrading enzyme activities to warming: implications for soil respiration. Global Change Biol. 24, 4816–4826 (2018).
Waring, B. G., Weintraub, S. R. & Sinsabaugh, R. L. Ecoenzymatic stoichiometry of microbial nutrient acquisition in tropical soils. Biogeochemistry 117, 101–113 (2014).
Mori, T., Lu, X., Aoyagi, R. & Mo, J. Reconsidering the phosphorus limitation of soil microbial activity in tropical forests. Funct. Ecol. 32, 1145–1154 (2018).
Gallardo, A. & Schlesinger, W. H. Factors limiting microbial biomass in the mineral soil and forest floor of a warm-temperate forest. Soil Biol. Biochem. 26, 1409–1415 (1994).
Feng, J. et al. Coupling and decoupling of soil carbon and nutrient cycles across an aridity gradient in the drylands of northern China: Evidence from ecoenzymatic stoichiometry. Global Biogeochem. Cycles. 33, 559–569 (2019).
Cui, Y. et al. Patterns of soil microbial nutrient limitations and their roles in the variation of soil organic carbon across a precipitation gradient in an arid and semi-arid region. Sci. Total Environ. 658, 1440–1451 (2019).
Jing, X. et al. Soil microbial carbon and nutrient constraints are driven more by climate and soil physicochemical properties than by nutrient addition in forest ecosystems. Soil Biol. Biochem. 141, 107657 (2020).
Fernández-Martínez, M. et al. Nutrient availability as the key regulator of global forest carbon balance. Nat. Clim. Change 4, 471–476 (2014).
Zhou, L. et al. Soil extracellular enzyme activity and stoichiometry in China’s forests. Funct. Ecol. 34, 1461–1471 (2020).
Fang, J., Chen, A., Peng, C., Zhao, S. & Ci, L. Changes in forest biomass carbon storage in China between 1949 and 1998. Science 292, 2320–2322 (2001).
Zhu, J. et al. Carbon stocks and changes of dead organic matter in China’s forests. Nat. Commun. 8, 1–10 (2017).
Fang, J., Yu, G., Liu, L., Hu, S. & Chapin, F. S. Climate change, human impacts, and carbon sequestration in China. Proc. Natl. Acad. Sci. USA 115, 4015–4020 (2018).
Sinsabaugh, R. L. et al. Stoichiometry of soil enzyme activity at global scale. Ecol. Lett. 11, 1252–1264 (2008).
Sinsabaugh, R. L., Hill, B. H. & Follstad Shah, J. J. Ecoenzymatic stoichiometry of microbial organic nutrient acquisition in soil and sediment. Nature 462, 795–798 (2009).
Moorhead, D. L., Sinsabaugh, R. L., Hill, B. H. & Weintraub, M. N. Vector analysis of ecoenzyme activities reveal constraints on coupled C, N and P dynamics. Soil Biol. Biochem. 93, 1–7 (2016).
Cui, Y. et al. Stoichiometric models of microbial metabolic limitation in soil systems. Global Ecol. Biogeogr. 30, 2297–2311 (2021).
Elser, J. J. et al. Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine, and terrestrial ecosystems. Ecol. Lett. 10, 1135–1142 (2007).
Schulte-Uebbing, L. & Vries, W. D. Global-scale impacts of nitrogen deposition on tree carbon sequestration in tropical, temperate, and boreal forests: a meta-analysis. Global Change Biol. 24, 416–431 (2017).
Richardson, S. J., Peltzer, D. A., Allen, R. B. & Parfitt, M. G. L. Rapid development of phosphorus limitation in temperate rainforest along the Franz josef soil chronosequence. Oecologia 139, 267–276 (2004).
Augusto, L., Achat, D. L., Jonard, M., Vidal, D. & Ringeval, B. Soil parent material-a major driver of plant nutrient limitations in terrestrial ecosystems. Global Change Biol. 23, 3808–3824 (2017).
Yao, Q. et al. Community proteogenomics reveals the systemic impact of phosphorus availability on microbial functions in tropical soil. Nat. Ecol. Evol. 2, 499–509 (2018).
Philippot, L., Raaijmakers, J. M., Lemanceau, P. & Putten, W. H. Going back to the roots: the microbial ecology of the rhizosphere. Nat. Rev. Microbiol. 11, 789–799 (2013).
Kuzyakov, Y. & Xu, X. Competition between roots and microorganisms for nitrogen: mechanisms and ecological relevance. New Phytol. 198, 656–669 (2013).
Cui, Y. et al. Ecoenzymatic stoichiometry and microbial nutrient limitation in rhizosphere soil in the arid area of the northern Loess Plateau, China. Soil Biol. Biochem. 116, 11–21 (2018).
Cui, Y. et al. Soil moisture mediates microbial carbon and phosphorus metabolism during vegetation succession in a semiarid region. Soil Biol. Biochem. 147, 107814 (2020).
Johnson, J. et al. The response of soil solution chemistry in european forests to decreasing acid deposition. Global Change Biol. 24, 3603–3619 (2018).
Janssens, I. A. et al. Reduction of forest soil respiration in response to nitrogen deposition. Nat. Geosci. 3, 315–322 (2010).
Penuelas, J. et al. Human-induced nitrogen-phosphorus imbalances alter natural and managed ecosystems across the globe. Nat. Commun. 4, 1–10 (2013).
Yu, G. et al. Stabilization of atmospheric nitrogen deposition in china over the past decade. Nat. Geosci. 12, 424–429 (2019).
Cui, Y. et al. Decreasing microbial phosphorus limitation increases soil carbon release. Geoderma 419, 115868 (2022).
Sinsabaugh, R. L., Moorhead, D. L., Xu, X. & Litvak, M. E. Plant, microbial and ecosystem carbon use efficiencies interact to stabilize microbial growth as a fraction of gross primary production. New Phytol. 214, 1518–1526 (2017).
Craig, M. E., Mayes, M. A., Sulman, B. N. & Walker, A. P. Biological mechanisms may contribute to soil carbon saturation patterns. Global Change Biol. 27, 2633–2644 (2021).
Friggens, N. L., Hester, A. J., Mitchell, R. J., Parker, T. C. & Wookey, P. A. Tree planting in organic soils does not result in net carbon sequestration on decadal timescales. Global Change Biol. 26, 5178–5188 (2020).
Jiang, M. et al. The fate of carbon in a mature forest under carbon dioxide enrichment. Nature 580, 227–231 (2020).
Rosinger, C., Rousk, J. & Sandén, H. Can enzymatic stoichiometry be used to determine growth-limiting nutrients for microorganisms?-A critical assessment in two subtropical soils. Soil Biol. Biochem. 128, 115–126 (2019).
Mori, T. Does ecoenzymatic stoichiometry really determine microbial nutrient limitations? Soil Biol. Biochem. 146, 107816 (2020).
Allison, S. D., Wallenstein, M. D. & Bradford, M. A. Soil-carbon response to warming dependent on microbial physiology. Nat. Geosci. 3, 336–340 (2010).
Saiya-Cork, K. R., Sinsabaugh, R. L. & Zak, D. R. The effects of long term nitrogen deposition on extracellular enzyme activity in an acer saccharum, forest soil. Soil Biol. Biochem. 34, 1309–1315 (2002).
German, D. P. et al. Optimization of hydrolytic and oxidative enzyme methods for ecosystem studies. Soil Biol. Biochem. 43, 1387–1397 (2011).
Lindstrom, M. J. & Bates, D. M. Newton-Raphson and EM Algorithms for Linear Mixed-Effects Models for Repeated-Measures Data. J. Am. Stat. Assoc. 83, 1014–1022 (1988).
Legendre, P. & Legendre, L. Numerical ecology, 2nd English edition. Elsevier Science BV, Amsterdam (1998).
Muggeo, V. M. R. Segmented: an R package to fit regression models with broken-line relationships. R News 8/1, 20–25 (2008).
Toms, J. D. & Lesperance, M. Piecewise regression: a tool for identifying ecological thresholds. Ecology 84, 2034–2041 (2003).
Borcard, D., Legendre, P. & Drapeau, P. Partialling out the spatial component of ecological variation. Ecology 73, 1045–1055 (1992).
Breiman, L. Random forests. Machine Learning 45, 5–32 (2001).
Sanchez, G., Trinchera, L. & Russolillo, G. plspm: Tools for Partial Least Squares Path Modeling (PLS-PM). R package version 0.4.7 edn (2016).
Development Core Team R. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2016).
Acknowledgements
Acknowledgment for the data support from “National Earth System Science Data Center, National Science & Technology Infrastructure of China” (http://www.geodata.cn). This study was financially supported by the Strategic Priority Research Program of Chinese Academy of Sciences (XDB40000000), Funds for International Cooperation and Exchange of National Natural Science Foundation of China (32061123007), National Natural Science Foundation of China (41977031), Program of State Key Laboratory of Loess and Quaternary Geology CAS (SKLLQGZR1803). Contributions from Dr. Chen are funded by H2020 Marie Skłodowska-Curie Actions (No. 839806). M.D.-B. acknowledges support from the Spanish Ministry of Science and Innovation for the I+D+i project PID2020-115813RA-I00 funded by MCIN/AEI/10.13039/501100011033. M.D.-B. is also supported by a project of the Fondo Europeo de Desarrollo Regional (FEDER) and the Consejería de Transformación Económica, Industria, Conocimiento y Universidades of the Junta de Andalucía (FEDER Andalucía 2014-2020 Objetivo temático “01–Refuerzo de la investigación, el desarrollo tecnológico y la innovación”) associated with the research project P20_00879 (ANDABIOMA).
Author information
Authors and Affiliations
Contributions
Y.C., H.B., and L.F. designed the research; Y.C., L.Y. performed analysis; Y.C. made all the figures; Y.C., J.Y., S.Z., and X.W. completed the experiment, Y.C., D.L.M., M.D.-B., L.Y., S.P., X.G., B.Z., and J.C. contributed to the interpretation of the results, and Y.C., L.M., M.D.-B., S.P., X.G., B.Z., J.C., W.T., Y.W., X.Z., and L.F. finished the editing of manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Communications Earth & Environment thanks Congde Huang, Tessa Camenzind and Zubin Xie for their contribution to the peer review of this work. Primary Handling Editors: Clare Davis.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This 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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Cui, Y., Bing, H., Moorhead, D.L. et al. Ecoenzymatic stoichiometry reveals widespread soil phosphorus limitation to microbial metabolism across Chinese forests. Commun Earth Environ 3, 184 (2022). https://doi.org/10.1038/s43247-022-00523-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s43247-022-00523-5
This article is cited by
-
Variation patterns and their driving factors in soil extracellular enzyme activities and stoichiometry along a 49-years vegetation restoration chronosequence
Plant and Soil (2024)
-
Regulation of soil C–N–P stoichiometry by intercropping mitigates microbial resource limitations and contributes to maize productivity
Plant and Soil (2023)
-
Divergent responses of soil microbial metabolic limitations to cropland revegetation at erosion and deposition topographies in the hilly-gully region of the northern Loess Plateau, China
Plant and Soil (2023)
-
Trade-off between microbial carbon use efficiency and specific nutrient-acquiring extracellular enzyme activities under reduced oxygen
Soil Ecology Letters (2023)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.