- 1College of Life Sciences, Xinjiang Normal University, Urumqi, China
- 2Institute of Ecology and Key Laboratory of Earth Surface Processes of Ministry of Education, College of Urban and Environmental Sciences, Peking University, Beijing, China
Changes in net primary productivity (NPP) to global change have been studied, yet the relative impacts of global change on grassland productivity at large scales remain poorly understood. Using 182 grassland samples established in 17 alpine meadows (AM) and 21 desert steppes (DS) in China, we show that NPP of AM was significantly higher than that of DS. NPP increased significantly with increasing leaf nitrogen content (LN) and leaf phosphorus content (LP) but decreased significantly with increasing leaf dry matter content (LDMC). Among all abiotic factors, soil nutrient factor was the dominant factor affecting the variation of NPP of AM, while the NPP of DS was mainly influenced by the changing of precipitation. All abiotic factors accounted for 62.4% of the spatial variation in the NPP of AM, which was higher than the ability to explain the spatial variation in the NPP of DS (43.5%). Leaf traits together with soil nutrients and climatic factors determined the changes of the grassland productivity, but the relative contributions varied somewhat among different grassland types. We quantified the effects of biotic and abiotic factors on grassland NPP, and provided theoretical guidance for predicting the impacts of global change on the NPP of grasslands.
Introduction
Grassland accounts for about a quarter of the global land area and is an important component of the terrestrial ecosystem (Pecina et al., 2019). Grasslands can not only effectively contributes to the global carbon cycle, but also plays a key role in regulating climate and soil conservation (Wang et al., 2019; Zhang et al., 2020). Grassland ecosystems are monospecific and structurally simple, making them more climate-sensitive than forest ecosystems (Liang and Gornish, 2019). Different types of grasslands (such as meadow steppe, typical steppe, desert steppe, and alpine steppe) also have obvious differences in ecosystem functions. Therefore, it is of great significance for us to investigate how environmental factors affect the functions of various grasslands in the context of global change (Anderegg et al., 2015; Seddon et al., 2016).
Net primary productivity (NPP) reflects the amount of carbon fixed by photosynthesis in a certain ecosystem per unit of time and space (Gang et al., 2018). NPP represents the energy that plants can use for growth, development, and reproduction, and is also the material basis for the survival and reproduction of biological groups (Sun and Du, 2017). Numerous studies have found that biological and abiotic factors jointly determine the temporal and spatial changes of NPP in grassland ecosystems (Shi et al., 2022). For example, Henry et al. (2018) found that precipitation controlled the spatio-temporal variation of grassland NPP, while Gong H. et al. (2020) believed that the availability of soil N and P content were the key factors that shaping the distribution pattern of grassland NPP. In recent years, with the development of functional trait research, leaf traits were widely used to predict the changes of NPP (Forrestel et al., 2017; Li et al., 2020; Schrader et al., 2021; Jiao et al., 2022). There remains no consensus on the dominant factors of spatial-temporal variation in grassland NPP (Henry et al., 2018; Wei et al., 2020).
Alpine meadows (AM) are mainly distributed in the alpine plateau region, which located at higher altitude regions. While desert steppes (DS) are mainly located in temperate regions with little mean annual precipitation (MAP) and are mainly distributed in northwest China (Zhou et al., 2020). Environmental factors driving the NPP of various grasslands vary considerably due to the difference in hydrothermal conditions and the altitude at which they are located (Sun et al., 2021). Drought stress significantly affects the photosynthesis and metabolism of the plants in DS (Grilli et al., 2017), therefore, precipitation is a key limiting factor of the NPP in DS. With the decline in MAP, the NPP of the DS decreased significantly (Chen et al., 2022). However, many studies have found that rainfall does not directly affect plant metabolism, but further influences plant uptake of soil nutrients by acting on soil microbial respiration and accelerating soil nutrient leaching and transformation (Gao et al., 2021). Soil is the direct living environment for plants. Plants absorb water and nutrients from the soil for photosynthesis. Soil nutrients control the photosynthesis of grass plants and affect the NPP of grasslands (Wieder et al., 2015). Soil N content can mitigate the effects on desert grassland NPP due to reduced precipitation, while soil P content is an essential element for plant energy transfer (Baldarelli et al., 2021), and soil pH represents soil fertility, which is a key soil factor to maintain NPP in desert grassland ecosystem (Zhang et al., 2021). Plants adapt to climate change by altering their traits (Wang et al., 2017b). Plants adapt to drought stress by reducing the leaf dry matter content (LDMC) and increased the leaf nitrogen content (LNC; Mahaut et al., 2020). Plants with lower specific leaf area (SLA) had a stronger ability to retain water and generally have a longer life span (Gong and Gao, 2019; Luo et al., 2020).
Compared with the DS, the structure and function of the AM are more sensitive to climate change (Pries et al., 2017). Zhang et al. (2019) found that light intensity significantly affects the NPP of the AM through the brightening experiment. NPP of the AM is also limited by low temperature and has no significant relationship with precipitation (Pan et al., 2021). Leaf traits represent the adaptability of alpine grassland function to low temperature, and plants with smaller leaf area (LA) are better able to avoid frost and heat damage (Dong et al., 2017). LN and LP are essential for plant growth and development (Heineman et al., 2016), and grassland NPP largely depends on plant availability to nitrogen and phosphorus (Wieder et al., 2015). LDMC is also closely related to grassland NPP, and Xu et al. (2018) found that grassland NPP is positively correlated with LDMC. The spatio-temporal variation of NPP of different grassland types is mainly influenced by a combination of environmental and biological factors, but there remains no consensus on the relative contribution of these factors (Gong Y. H. et al., 2020; Jochum et al., 2020; Shi et al., 2020).
Based on the field data from 182 grasslands in 17 AM and 21 DS in China, we sought to identify the main drivers of NPP in different grassland ecosystems. To answer this question, we proposed the following hypotheses: (1) The NPP of AM is significantly higher than DS. (2) Precipitation is the driving factor of NPP of DS, while NPP of AM is the combined effect of precipitation and temperature. (3) Climate factors are the key environmental factors dominating NPP of AM and that of DS, however, soil nutrient factors also play a non-negligible role.
Materials and methods
Study area and plot description
China has a vast territory and a diverse climate, with possesses many types of grasslands. These grasslands are mainly distributed in arid and semi-arid areas in northern China and the Qinghai-Tibet Plateau (Yu et al., 2021; Shen et al., 2022). The large spatial scale and the richness of climate types provide conditions for exploring the spatio-temporal variation of NPP of grasslands (Huang et al., 2021). The study used data from 182 grasslands in 17 AM in the southwest region and 21 DS in the north central region in China (Figure 1A). Our study sites range in latitude from 19.1° to 53.5° north and longitude from 79.72° to 121.1° east, with elevations ranging from 13 to 5,000 m.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996313/fpls-13-996313-HTML-r1/image_m/fpls-13-996313-g001.jpg)
Figure 1. Comparison of the variability of different grassland NPP. (A) Geographical location of the grassland plots in this study; (B) variability of the NPP of AM and DS. AM represents alpine meadows and DS represents desert steppes. ****Represents P < 0.0001.
Community leaf traits data
Leaf traits are widely used to detect the adaptability of plants to the living environment (Hikosaka et al., 2021). The traits we selected included LA, SLA, LDMC, LN, LP, and leaf N/P ratio. These traits are closely related to light interception, resource utilization, growth, and development of plants (Brun et al., 2022). Previous studies at leaf traits usually focused on the species level (Butler et al., 2017). However, the effects of competitive asymmetry (i.e., the magnitude and consequences of competition affecting each party) related to interspecies functional characteristics has been neglected (Vázquez et al., 2007). We used community weighted average trait (CWMi) to represent average trait values in grassland.
where CWMi represents the weighted trait value of the community, and Di represents the abundance of the target species.
Environmental factors and net primary productivity data
Numerous studies have found that grassland NPP is influenced by climate factors, especially mean annual temperature (MAT) and MAP (Sloat et al., 2018). MAT, mean coldest monthly temperature (MCMT), mean warmest monthly temperature (MWMT), and MAP were extracted from the WorldClim global climate layer at a spatial resolution of 1 km (Conradi et al., 2020). Light is a key climatic factor affecting photosynthesis (Kramer, 1981). We hypothesized that annual sunshine hours (ASD) may be a key predictor of NPP. ASD and mean annual evaporation (MAE) were also obtained from the Meteorological Data Center of the China Meteorological Administration.1
We extracted from 250 m resolution of the grid in the top 30 cm soil layer soil pH, soil nitrogen,2 and the content of soil effective phosphorus.3
NPP data from the national aeronautics and space administration (NASA4), the site offers from 2000 to 2015, a resolution of 250 m by 250 m of NPP. This dataset was derived from the widely used Medium Resolution Imaging Spectroradiometer (MOD13Q1) product, calculated using the C5 MOD17 algorithm, and verified by data from flflux towers (Li et al., 2020).
Data analysis
We used a significant difference test at the 0.05 significance level to test whether there was a significant difference between NPP in AM and DS (Figure 1). We also tested CWMi differences between different grassland types at α = 0.05 (Chen et al., 2020; Supplementary Figure 1). Significant difference tests were performed using the R package agricolae (version 4.1.0, R Core Team, 2020).
We tested the effects of environmental factors and community functional traits on NPP using general linearity. R2 represents the goodness of fit of the model. Linear regression was performed using the R package lme4 (R Core Team, 2020).
Our preliminary results indicate that most functional traits have no significant linear relationship with NPP (Figure 2). Therefore, in subsequent studies, we mainly explore the impact of environmental factors on NPP.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996313/fpls-13-996313-HTML-r1/image_m/fpls-13-996313-g002.jpg)
Figure 2. General linear correlations between community leaf traits and the NPP of AM and DS. (A) General linear relationship between LA and the NPP of AM and DS; (B) general linear relationship between SLA and the NPP of AM and DS; (C) general linear relationship between LDMC and the NPP of AM and DS; (D) general linear relationship between LN and the NPP of AM and DS; (E) general linear relationship between LP and the NPP of AM and DS; (F) general linear relationship between N/P and the NPP of AM and DS. LA represents leaf area; SLA represents specific leaf area; LDMC represents leaf dry matter content; LN represents leaf nitrogen content; LP represents leaf phosphorus content; N/P represents leaf nitrogen and phosphorus ratio. The red line represents the NPP of AM and the green line represents the NPP of DS; R2 represents the fit of the model, and P represents the correlation.
We used a generalized additive model (GAM) to test the effects of climate and soil factors on NPP. This model is composed of parametric variables and non-parametric variables (Zou et al., 2016). Non-metric multidimensional scaling analysis (NMDS) was used to reflect the fit of GAM (Zou et al., 2016).
where g(●) denotes the connection function, the form of which depends on the specific form and can be interpreted as the variable Y distribution. E is a random error term that can be interpreted as a normally distributed function named constant variable connection, and the connection function takes the form g(u) = u, u = E(Y | X), E(E| X) = 0. Xi is the explanatory variable that strictly follows the parametric form of the explanatory variables, βi is the corresponding parameter, and fj(●) is the smoothing function corresponding to the explanatory variable Xj that follows a non-parametric form. In our study, the spline smoothing function S(●) is chosen to fit, thin-slab spline smoothing is chosen to fit the function between different nodes, and least squares is used to estimate each smoothing function S(●).
Results
Effects of leaf traits on net primary productivity of grasslands
The NPP of AM was significantly higher than that of DS (P < 0.001), and most functional traits of AM and DS were also significantly different (Figure 1B). Except that the N/P of AM was lower than that of DS, other leaf traits were significantly higher than that of DS (Supplementary Figure 1). The NPP of different grassland types increased significantly with the increase of LDMC (Figure 2C), and decreased significantly with the increase of LN and LP (Figures 2D,E). NPP of AM and DS showed opposite trends with the increase of LA, SLA, and leaf N/P ratio (Supplementary Figures 1A,B). With the increase of the above leaf traits, NPP increased significantly of AM, while NPP decreased significantly of DS. Among all leaf traits, LDMC had better predictive power for NPP of AM (R2 = 0.35, P < 0.001; Figure 2C) and LA had better predictive power for NPP of DS (R2 = 0.13, P < 0.001; Figure 2A).
Effects of climate and soil nutrients on net primary productivity of grasslands
The NPP of AM was significantly positively correlated with MAT and MCMT and negatively correlated with MAE. The NPP of DS was significantly positively correlated with both MWMT and MAP (Figure 3). MAT (R2 = 0.28, P < 0.001) was the best climate factor in predicting the change of NPP of AM. MAP (R2 = 0.22, P < 0.001) was the best climate factor in predicting the change of NPP of DS.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996313/fpls-13-996313-HTML-r1/image_m/fpls-13-996313-g003.jpg)
Figure 3. General linear regression analysis between climate factors and the NPP of AM and DS. (A) General linear relationship between MAT and the NPP of AM and DS; (B) general linear relationship between MCMT and the NPP of AM and DS; (C) general linear relationship between MWMT and the NPP of AM and DS; (D) general linear relationship between MAP and the NPP of AM and DS; (E) general linear relationship between ASD and the NPP of AM and DS; (F) general linear relationship between MAE and the NPP of AM and DS. MAT represents mean annual temperature; MCMT represents mean coldest monthly temperature; MWMT represents mean warmest monthly temperature; MAP represents mean annual precipitation; ASD represents annual sunshine hours; MAE represents mean annual evaporation. The shaded area shows a 95% confidence interval. The red line represents the NPP of AM, and the green line represents the NPP of DS; R2 represents the model fit, and P represents the correlation.
We further analyzed the effects of soil nutrient factors on the NPP of different types of grasslands (Figure 4). The prediction ability of soil N, soil P, and soil pH on NPP of AM were 0.30, 0.30, and 0.40, respectively. Among all soil nutrient factors, except the positive correlation between soil N and the NPP of DS, the other factors have no significant effect (P > 0.05).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996313/fpls-13-996313-HTML-r1/image_m/fpls-13-996313-g004.jpg)
Figure 4. General linear regression analysis between soil nutrient factors and the NPP (soil N, soil P and soil PH) of AM and DS. (A) General linear relationship between soil N and the NPP of AM and DS; (B) general linear relationship between soil P and the NPP of AM and DS; (C) general linear relationship between soil pH and the NPP of AM and DS. The shaded area shows a 95% confidence interval. The red line represents the NPP of AM and the green line represents the NPP of DS; R2 represents the model fit and P represents the correlation.
Soil factors dominate the variation of grassland net primary productivity
We analyzed the effects of abiotic factors on the NPP of AM and DS based on a generalized additive model with non-metric multidimensional scaling (NMDS) ordering. All environmental factors jointly explained 62.4% of NPP in AM and 43.5% of NPP in DS (Figures 5C,F). Among them, the soil factor explained the variation of NPP in different types of grasslands to a greater extent than the climate factor (Figures 5A,B,D,E). Soil factors are the dominant factors driving the spatial and temporal variation of NPP in grasslands, and climate factors also play an important role.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/996313/fpls-13-996313-HTML-r1/image_m/fpls-13-996313-g005.jpg)
Figure 5. The results of the generalized additive model (GAM) fitted the effects of environmental factors and climate factors on grassland NPP. (A) NMDS ranking of soil factors and the NPP of AM; (B) NMDS ranking of climate factors and the NPP of AM; (C) NMDS ranking of environmental factors and the NPP of AM; (D) NMDS ranking of soil factors and the NPP of DS; (E) NMDS ranking of climate factors and the NPP of DS; (F) NMDS ranking of environmental factors with the NPP of DS. The color spline shows the fit of trait data from high (green) to low (yellow) values during the ranking. Trait overlays indicate that abiotic factors represented by points on the NMD are associated with higher or lower trait values, consistent with the colored trait gradient. Note that if the relationship between NPP and abiotic factors is linear, the gradient spline lines will be parallel. non-linear relationships between NPP and abiotic factors are represented by curve splines; de indicates the deviation explained by the corresponding model.
Discussion
There were significant spatial differences in NPP of different grassland types, and the NPP of AM was significantly higher than that of DS. Significant differences in NPP were closely related to habitat heterogeneity (Hikosaka et al., 2021). Due to the differences in biological or abiotic factors in different regions, there are significant differences in NPP of grassland in different regions (Li et al., 2019). AM are mainly distributed in the cold and high-altitude area of southwest China, with low temperature, strong radiation, thin soil layer, long soil freezing period, and neutral soil reaction. And DS are mainly distributed in the arid and semi-arid areas in the middle and north of China, which is characterized by lack of rainfall, strong evaporation, loose soil, lack of water, and low organic matter content. Therefore, differences in local hydrothermal conditions and community structure may be important reasons for the significant differences in NPP between AM and DS.
A large number of studies have found that the community leaf traits are closely related to the NPP of grasslands (Wang et al., 2017b). Our results found that the NPP of different grasslands increased significantly with the decreased of LN and LP, and increased significantly with the increase of LDMC. Nitrogen is the main component of chlorophyll, which is used to accumulate organic matter through light reactions. Phosphorus promotes the accumulation of glucose phosphate and organic matter through phosphorylation. Heineman et al. (2016) found that LN and LP limit the accumulation of plant organic matter. The NPP of the grasslands is largely determined by LN and LP contents (Oldroyd and Leyser, 2020). Studies have shown that LN is closely related to plant growth and defense, and that the relationship between grassland NPP and LN and LP is affected by the relative distribution of these nutrients among plant tissues (Zheng et al., 2022). The higher resource utilization and competitiveness of plants are beneficial to improve NPP, however, plants in arid, water-deficient, and alpine environments usually have lower LN and higher LDMC (Yang et al., 2022). Therefore, even with an increase in LN and LP, NPP did not increase with it (Xu et al., 2018), which explains our results well. We also found that NPP of grassland ecosystems was positively correlated with LDMC (Forrestel et al., 2017). As LDMC increases, free water required for plant growth and development is gradually converted to bound water, which is involved in plant metabolic processes (e.g., photosynthesis). The NPP of DS decreased significantly with the increase of LA and SLA. Plants with small and thick leaves (smaller SLA) have better drought tolerance and were more conducive to cell material transport rate and photosynthetic efficiency (Poorter and Markesteijn, 2008; Osnas et al., 2018). These traits reflect the optimal utilization of limited resources by plants (Gupta et al., 2020). Although MAT of AM was significantly lower than DS (Supplementary Figure 2A), MAP (Supplementary Figure 2D), soil nitrogen content (Supplementary Figure 3A) and soil phosphorus content (Supplementary Figure 3B) were significantly higher than DS. This reflects that relatively abundant rainwater and soil nutrients can alleviate the impact of low-temperature stress on NPP (Berdugo et al., 2020).
Climatic factors are the main driving factors of spatial variation of grassland NPP (Wang et al., 2017a). We found that NPP of AM was significantly positively correlated with MAT and MCMT, and significantly negatively correlated with MAE. Numerous studies have shown that changes in NPP of AM are strongly correlated with MAT and MCMT (Keenan and Riley, 2018; Sloat et al., 2018; Shi et al., 2022). Temperature can not only directly affect the enzyme activities related to photosynthesis and respiration of plants, but also affects the production of organic compounds such as starch, protein, and lipid. Niu et al. (2019) found that global warming can improve the productivity of grassland in alpine regions for a long time. For NPP of AM, the influence of solar radiation is no less than that of temperature (Zheng et al., 2020), and the intensity of solar radiation is also one of the important factors affecting plant photosynthesis, which is the energy source of the grassland ecosystem (Brun et al., 2022). However, excessive evaporation triggered by too high radiation zone intensity can weaken the positive response of plants to high temperatures, causing plant water deficit and thus reducing NPP. Therefore, with the increase of MAE, the NPP of AM decreased significantly, and precipitation was the main climate factor leading to the temporal and spatial variation of the NPP of DS. Our conclusion also has been verified by a large number of previous studies (Liang et al., 2015; Erb et al., 2018; Guo et al., 2022). As the MAP increases, the NPP of AM increases significantly. Water restriction will seriously reduce the material transport efficiency and photosynthetic productivity efficiency of plants in arid areas, and reduce the organic matter yield of leaves and significantly affecting the productivity of the whole region (Sun et al., 2021). Significant difference analysis of climate factors of different grassland types found that the differences in climate factors (e.g., MAT, MWMT, MAP, and MAE; Supplementary Figure 2) and the differences in NPP (Figure 1B) showed surprising consistency, which also indirectly indicated that climate factors were the key environmental factors affecting grassland NPP (Yang et al., 2016).
Soil is a natural “storage box,” providing most of the nutrients needed for plants’ growth and reproduction. Soil nutrient elements are closely related to plant leaf nutrient elements (Gao et al., 2019). Soil nutrient factors were significantly positively correlated with community leaf traits such as LA, SLA, and LDMC (Huang et al., 2021). Climate affects the function of grassland ecosystems by changing the physicochemical properties of the soil. These results indicate that soil factors interact with other environmental factors to jointly affect grasslands NPP (Gong and Gao, 2019). With the increase of soil N and P, the NPP of AM decreased significantly. Soil N and P are essential nutrients for plant growth and development, which limit plant productivity (Gong H. et al., 2020). The alpine environment limits the decomposition capacity of soil microorganisms, and relatively low soil mineralization rate and soil nutrient availability inhibits the NPP of AM (Wang et al., 2017c). Nitrogen addition has been found to inhibit microbial respiration and thus affect the rate of microbial decomposition (Wang et al., 2018). Soil nutrient availability, the main soil factor, leading to change the NPP of AM (Blowes et al., 2019). Higher pH reduces soil nutrient availability by affecting soil organic matter storage (Chen et al., 2018). The availability of soil water and soil nutrients limit the NPP of DS. Precipitation can promote mineralization and improve soil nutrients, thereby increasing plant NPP, which is the limiting factor for plants to absorb soil nutrients. Jiang et al. (2017) also showed that NPP in arid areas is positively correlated with precipitation, and soil water can promote the decomposition of soil organic matter and thus improve the availability of soil nutrients. These findings explain the weak predictive power of soil N, soil P and soil pH for NPP of DS. Water shortage and strong evaporation lead to slow soil development and high soil pH, resulting in insufficient soil fertility and water (Zhao et al., 2022).
The NPP of AM (de = 62.4%) has stronger environment plasticity than DS (de = 43.5%). Regardless of AM or DS, the contribution of soil nutrients in shaping the NPP pattern was slightly higher than that of climate factors (Figure 5). Climatic factors, especially MAP and MAT, have long been considered key factors in grassland NPP (Erb et al., 2018). However, more and more studies have proved that soil properties play a dominant role in grassland NPP (Radujković et al., 2021). Temperature and precipitation control plant growth by changing microbial activity and soil physicochemical properties, thus affecting the NPP of grassland ecosystems (Newsham et al., 2016). Drought and frozen soil will lead to the decrease of soil mineralization rate and soil fertility (Nogueira et al., 2018). Soil is the habitat for plants, where providing water, inorganic matter, and organic matter for plant growth and development, and is one of the core elements linking the entire ecosystem. Climate change directly affects the soil environment and indirectly affects grassland NPP. Therefore, soil factors dominated the temporal and spatial changes of NPP in different grassland types on a large scale.
The wide study area and large amount of data in this paper ensure the credibility of our experimental results. Our study quantifies the effects of plant functional traits and environmental factors on NPP of AM and DS, providing important theoretical guidance for addressing global climate change and better management of grasslands. In future research, we will further study and explore more soil factors, such as soil microorganisms and soil N and P effectiveness on grassland NPP.
Conclusion
We examined the relative roles of biological and abiotic factors in shaping NPP patterns using the data from 182 grassland plots in 17 DS and 21 AM in China. Our results confirmed that there are significant differences in productivity among different grasslands. Leaf traits, as well as climate and soil nutrient factors, jointly affect the changes of the grassland NPP. Among them, soil nutrients play the most critical role. Soil nutrient factors explained 61.8% of the spatial variation in NPP of AM and 44.4% of the spatial variation in NPP of DS. Quantifying the main driving factors of NPP of different grasslands is crucial for predicting future grassland dynamics and proposing reasonable management strategies.
Data availability statement
The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
JG: conceptualization, methodology, and investigation. XW, RW, and JG: formal analysis. XW: writing – original draft. All authors have read and agreed to the published version of the manuscript.
Acknowledgments
We thank reviewers for their valuable comments on this article. We also thank Xinjiang Normal University for providing scientific research funds for this research.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.996313/full#supplementary-material
Supplementary Figure 1 | Comparison of functional traits among different grassland types. (A) Variability between LA and the NPP of AM and DS; (B) variability between SLA and the NPP of AM and DS; (C) variability between LDMC and the NPP of AM and DS; (D) variability between LN and the NPP of AM and DS; (E) variability between LP and the NPP of AM and DS; (F) variability between N/P and the NPP of AM and DS. AM represents alpine meadows; DS represents deserts steppe; ns represents the insignificant difference; and *represents 0.01 < P < 0.05. SLA represents specific leaf area; LDMC represents leaf dry matter content; LN represents leaf nitrogen content; LP represents leaf phosphorus content; N/P represents leaf nitrogen and phosphorus ratio.
Supplementary Figure 2 | Comparison of AM and DS with different climate factors. (A) Variability between MAT and the NPP of AM and DS; (B) variability between MCMT and the NPP of AM and DS; (C) variability between MWMT and the NPP of AM and DS; (D) variability between MAP and the NPP of AM and DS; (E) variability between ASD and the NPP of AM and DS; (F) variability between MAE and the NPP of AM and DS. AM represents alpine meadows; DS represents desert steppes; ns represents no significant difference between them; and ****represents P < 0.0001. MAT represents annual mean temperature; MCMT represents mean coldest monthly temperature; MWMT represents mean warmest monthly temperature; MAP represents mean annual precipitation; ASD represents annual sunshine hours; MAE represents average annual evaporation.
Supplementary Figure 3 | Comparison of different soil nutrient factors (soil N, soil P, and soil PH) between AM and DS. (A) Variability between soil N and the NPP of AM and DS; (B) variability between soil P and the NPP of AM and DS; (C) variability between soil pH and the NPP of AM and DS. AM represents alpine meadows; DS represents desert steppes; ns represents the insignificant difference between the two; and ****represents P < 0.0001.
Footnotes
- ^ https://data.cma.cn/
- ^ http://www.csdn.store
- ^ https://www.osgeo.cn/data/wc137
- ^ https://search.earthdata.NASA.gov/search
References
Anderegg, W. R. L., Schwalm, C., Biondi, F., Camarero, J. J., Koch, G., and Litvak, M. (2015). Pervasive drought legacies in forest ecosystems and their implications for carbon cycle models. Science 349, 528–532. doi: 10.1126/science.aab1833
Baldarelli, L. M., Throop, H. L., Collins, S. L., and Ward, D. (2021). Nutrient additions have direct and indirect effects on biocrust biomass in a long-term Chihuahuan Desert grassland experiment. J. Arid Environ. 184:104317. doi: 10.1016/j.jaridenv.2020.104317
Berdugo, M., Delgado-Baquerizo, M., Soliveres, S., Hernández-Clemente, R., Zhao, Y. C., and Gaitán, J. J. (2020). Global ecosystem thresholds driven by aridity. Science 367, 787–790. doi: 10.1126/science.aay5958
Blowes, S. A., Supp, S. R., Antão, L. H., Bates, A., Bruelheide, H., and Chase, J. M. (2019). The geography of biodiversity change in marine and terrestrial assemblages. Science 366, 339–345. doi: 10.1126/science.aaw1620
Brun, P., Violle, C., Mouillot, D., Mouquet, N., Enquist, B. J., and Munoz, F. (2022). Plant community impact on productivity: Trait diversity or key (stone) species effects? Ecol. Lett. 25, 913–925. doi: 10.1111/ele.13968
Butler, E. E., Datta, A., Flores-Moreno, H., and Reich, P. B. (2017). Mapping local and global variability in plant trait distributions. Proc. Natl. Acad. Sci. U.S.A. 114:E10937–E10946. doi: 10.1073/pnas.1708984114
Chen, L., Hänninen, H., Rossi, S., Smith, N. G., Pau, S., and Liu, Z. Y. (2020). Leaf senescence exhibits stronger climatic responses during warm than during cold autumns. Nat. Clim. Chang. 10, 777–780. doi: 10.1038/s41558-020-0820-2
Chen, S. P., Wang, W. T., Xu, W. T., Wang, Y., Wan, H. W., and Bai, Y. F. (2018). Plant diversity enhances productivity and soil carbon storage. Proc. Natl. Acad. Sci. U.S.A. 115, 4027–4032. doi: 10.1073/pnas.1700298114
Chen, X. Y., Liu, P., Zhao, B., Zhang, J. W., Ren, B. Z., and Li, Z. (2022). Root physiological adaptations that enhance the grain yield and nutrient use efficiency of maize (Zea mays L) and their dependency on phosphorus placement depth. Field Crops Res. 276:108378. doi: 10.1016/j.fcr.2021.108378
Conradi, T., Meerbeek, K. V., Ordonez, A., and Svenning, J. C. (2020). Biogeographic historical legacies in the net primary productivity of Northern Hemisphere forests. Ecol. Lett. 23, 800–810. doi: 10.1111/ele.13481
Dong, N., Prentice, I. C., Harrison, S. P., Song, Q. H., and Zhang, Y. P. (2017). Biophysical homoeostasis of leaf temperature: A neglected process for vegetation and land-surface modelling. Glob. Ecol. Biogeogr. 26, 998–1007. doi: 10.1111/geb.12614
Erb, K. H., Kastner, T., Plutzar, C., Bais, A. L. S., Carvalhais, N., and Fetzel, T. (2018). Unexpectedly large impact of forest management and grazing on global vegetation biomass. Nature 553, 73–76. doi: 10.1038/nature25138
Forrestel, E. J., Donoghue, M. J., Edwards, E. J., Jetz, W., and Smith, M. D. (2017). Different clades and traits yield similar grassland functional responses. Proc. Natl. Acad. Sci. U.S.A. 114, 705–710. doi: 10.1073/pnas.1612909114
Gang, C. C., Zhao, W., Zhao, T., Zhang, Y., Gao, X. R., and Wen, Z. M. (2018). The impacts of land conversion and management measures on the grassland net primary productivity over the Loess Plateau. Northern China. Sci. Total Environ. 645, 827–836. doi: 10.1016/j.scitotenv.2018.07.161
Gao, J., Song, Z. P., and Liu, Y. H. (2019). Response mechanisms of leaf nutrients of endangered plant (Acer catalpifolium) to environmental factors varied at different growth stages. Glob. Ecol. Conserv. 17:e00521. doi: 10.1016/j.gecco.2019.e00521
Gao, X. L., Li, X. G., Zhao, L., and Kuzyakov, Y. (2021). Shrubs magnify soil phosphorus depletion in Tibetan meadows: Conclusions from C: N: P stoichiometry and deep soil profiles. Sci. Total Environ. 785:147320. doi: 10.1016/j.scitotenv.2021.147320
Gong, H., and Gao, J. (2019). Soil and climatic drivers of plant SLA (specific leaf area). Glob. Ecol. Conserv. 20:e00696. doi: 10.1016/j.gecco.2019.e00696
Gong, H., Li, Y. Y., Yu, T., Zhang, S. H., Gao, J., and Zhang, S. K. (2020). Soil and climate effects on leaf nitrogen and phosphorus stoichiometry along elevational gradients. Glob. Ecol. Conserv. 23:e01138. doi: 10.1016/j.gecco.2020.e01138
Gong, Y. H., Zhao, D. M., Ke, W. B., Fang, C., Pei, J. Y., and Ye, J. S. (2020). Legacy effects of precipitation amount and frequency on the aboveground plant biomass of a semi-arid grassland. Sci. Total Environ. 705:135899. doi: 10.1016/j.scitotenv.2019.135899
Grilli, J., Barabas, G., Michalska-Smith, M. J., and Allesina, S. (2017). Higher-order interactions stabilize dynamics in competitive network models. Nature 548, 210–213. doi: 10.1038/nature23273
Guo, X. X., Zuo, X. A., Yue, P., Li, X. Y., Hu, Y., and Chen, M. (2022). Direct and indirect effects of precipitation change and nutrients addition on desert steppe productivity in Inner Mongolia, Northern China. Plant Soil 471, 527–540. doi: 10.1007/s11104-021-05206-2
Gupta, A., Rico-Medina, A., and Caño-Delgado, A. I. (2020). The physiology of plant responses to drought. Science 368, 266–269. doi: 10.1126/science.aaz7614
Heineman, K. D., Turner, B. L., and Dalling, J. W. (2016). Variation in wood nutrients along a tropical soil fertility gradient. New Phytol. 211, 440–454. doi: 10.1111/nph.13904
Henry, H. A. L., Abedi, M., Alados, C. L., Beard, K. H., Fraser, L. H., and Jentsch, A. (2018). Increased soil frost versus summer drought as drivers of plant biomass responses to reduced precipitation: Results from a globally coordinated field experiment. Ecosystems 21, 1432–1444. doi: 10.1007/s10021-018-0231-7
Hikosaka, K., Kurokawa, H., Arai, T., Takayanagi, S., Tanaka, H. O., and Nagano, S. (2021). Intraspecific variations in leaf traits, productivity and resource use efficiencies in the dominant species of subalpine evergreen coniferous and deciduous broad-leaved forests along the altitudinal gradient. J. Ecol. 109, 1804–1818. doi: 10.1111/1365-2745.13603
Huang, W. J., Wang, W., Cao, M., Fu, G., Xia, J. Y., and Wang, Z. X. (2021). Local climate and biodiversity affect the stability of China’s grasslands in response to drought. Sci. Total Environ. 768:145482. doi: 10.1016/j.scitotenv.2021.145482
Jiang, L. L., Jiapaer, G., Bao, A. M., Guo, H., and Ndayisaba, F. (2017). Vegetation dynamics and responses to climate change and human activities in Central Asia. Sci. Total Environ. 599, 967–980. doi: 10.1016/j.scitotenv.2017.05.012
Jiao, Y., Niklas, K. J., Wang, L., Yu, K., Li, Y., and Shi, P. (2022). Influence of leaf age on the scaling relationships of lamina mass vs. area. Front. Plant Sci. 13:860206. doi: 10.3389/fpls.2022.860206
Jochum, M., Fischer, M., Isbell, F., Roscher, C., Plas, F. V. D., and Boch, S. (2020). The results of biodiversity-ecosystem functioning experiments are realistic. Nat. Ecol. Evol. 4, 1485–1494. doi: 10.1038/s41559-020-1280-9
Keenan, T. F., and Riley, W. J. (2018). Greening of the land surface in the world’s cold regions consistent with recent warming. Nat. Clim. Chang. 8, 825–828. doi: 10.1038/s41558-018-0258-y
Kramer, P. J. (1981). Carbon dioxide concentration, photosynthesis, and dry matter production. Bioscience 31, 29–33. doi: 10.2307/1308175
Li, F. R., Liu, L. L., Liu, J. L., and Yang, K. (2019). Abiotic and biotic controls on dynamics of labile phosphorus fractions in calcareous soils under agricultural cultivation. Sci. Total Environ. 681, 163–174. doi: 10.1016/j.scitotenv.2019.05.091
Li, Y. Q., Reich, P. B., Schmid, B., Shrestha, N., Feng, X., and Wang, Z. H. (2020). Leaf size of woody dicots predicts ecosystem primary productivity. Ecol. Lett. 23, 1003–1013. doi: 10.1111/ele.13503
Liang, M. W., and Gornish, E. S. (2019). Rainfall regulation of grazed grasslands. Proc. Natl. Acad. Sci. U.S.A. 116, 23887–23888. doi: 10.1073/pnas.1915399116
Liang, W., Yang, Y. T., Fan, D. M., Guan, H. D., Zhang, T., and Bai, D. (2015). Analysis of spatial and temporal patterns of net primary production and their climate controls in China from 1982 to 2010. Agric. For. Meteorol. 204, 22–36. doi: 10.1016/j.agrformet.2015.01.015
Luo, W. T., Zuo, X. A., Griffin-Nolan, R. J., Xu, C., Sardans, J., and Yu, Q. (2020). Chronic and intense droughts differentially influence grassland carbon-nutrient dynamics along a natural aridity gradient. Plant Soil. 473, 37–148. doi: 10.1007/s11104-020-04571-8
Mahaut, L., Fort, F., Violle, C., and Freschet, G. T. (2020). Multiple facets of diversity effects on plant productivity: Species richness, functional diversity, species identity and intraspecific competition. Funct. Ecol. 34, 287–298. doi: 10.1111/1365-2435.13473
Newsham, K. K., Hopkins, D. W., Carvalhais, L. C., Fretwell, P. T., Rushton, S. P. and O’Donnell, A. G. (2016). Relationship between soil fungal diversity and temperature in the maritime Antarctic. Nat. Clim. Change 6, 182–186. doi: 10.1038/nclimate2806
Niu, Y. J., Yang, S. W., Zhou, J. W., Chu, B., Ma, S. J., and Hua, L. M. (2019). Vegetation distribution along mountain environmental gradient predicts shifts in plant community response to climate change in alpine meadow on the Tibetan Plateau. Sci. Total Environ. 650, 505–514. doi: 10.1016/j.scitotenv.2018.08.390
Nogueira, C., Nunes, A., Bugalho, M. N., Branquinho, C., McCulley, R. L., and Caldeira, M. C. (2018). Nutrient addition and drought interact to change the structure and decrease the functional diversity of a Mediterranean grassland. Front. Ecol. Evol. 6:00155. doi: 10.3389/fevo.2018.00155
Oldroyd, G. E. D., and Leyser, O. (2020). A plant’s diet, surviving in a variable nutrient environment. Science 368:eaba0196. doi: 10.1126/science.aba0196
Osnas, J. L., Katabuchi, M., Kitajima, K., Wright, S. J., Reich, P. B., and Lichstein, J. W. (2018). Divergent drivers of leaf trait variation within species, among species, and among functional groups. Proc. Natl. Acad. Sci. U.S.A. 115, 5480–5485. doi: 10.1073/pnas.1803989115
Pan, Q. M., Symstad, A. J., Bai, Y. F., Huang, J. H., Wu, J. G., and Han, X. G. (2021). Biodiversity–productivity relationships in a natural grassland community vary under diversity loss scenarios. J. Ecol. 110, 210–220. doi: 10.1111/1365-2745.13797
Pecina, M. V., Ward, R. D., Bunce, R. G. H., Sepp, K., Kuusemets, V., and Luuk, O. (2019). Country-scale mapping of ecosystem services provided by semi-natural grasslands. Sci. Total Environ. 661, 212–225. doi: 10.1016/j.scitotenv.2019.01.174
Poorter, L., and Markesteijn, L. (2008). Seedling traits determine drought tolerance of tropical tree species. Biotropica 40, 321–331. doi: 10.1111/j.1744-7429.2007.00380.x
Pries, C. H., Castanha, C., Porras, R. C., and Torn, M. S. (2017). The whole-soil carbon flux in response to warming. Science 355, 1420–1423. doi: 10.1126/science.aal1319
R Core Team. (2020). A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Radujković, D., Verbruggen, E., Seabloom, E. W., Bahn, M., Biederman, L. A., and Borer, E. T. (2021). Soil properties as key predictors of global grassland production: Have we overlooked micronutrients? Ecol. Lett. 24, 2713–2725. doi: 10.1111/ele.13894
Schrader, J., Shi, P. J., Royer, D. L., Peppe, D. J., Gallagher, R. V., and Li, Y. R. (2021). Leaf size estimation based on leaf length, width and shape. Ann. Bot. 128, 395–406. doi: 10.1093/aob/mcab078
Seddon, A. W., Macias-Fauria, M., Long, P. R., Benz, D., and Willis, K. J. (2016). Sensitivity of global terrestrial ecosystems to climate variability. Nature 531, 229–232. doi: 10.1038/nature16986
Shen, X. J., Liu, B. H., Henderson, M., Wang, L., Jing, M., and Lu, X. G. (2022). Vegetation greening, extended growing seasons, and temperature feedbacks in warming temperate grasslands of China. J. Clim. 35, 5103–5117. doi: 10.1175/JCLI-D-21-0325.1
Shi, F., Liu, M. X., Qiu, J., Zhang, Y. L., Su, H. Y., and Mao, X. P. (2022). Assessing land cover and ecological quality changes in the Forest-Steppe Ecotone of the Greater Khingan Mountains, Northeast China, from Landsat and MODIS observations from 2000 to 2018. Remote Sens. 14, 725–725. doi: 10.3390/rs14030725
Shi, P. J., Preisler, H. K., Quinn, B. K., Zhao, J., Huang, W., and Röll, A. (2020). Precipitation is the most crucial factor determining the distribution of moso bamboo in Mainland China. Glob. Ecol. Conserv. 22:e00924. doi: 10.1016/j.gecco.2020.e00924
Sloat, L. L., Gerber, J. S., Samberg, L. H., Smith, W. K., Herrero, M., and Ferreira, L. G. (2018). Increasing importance of precipitation variability on global livestock grazing lands. Nat. Clim. Chang. 8, 214–218. doi: 10.1038/s41558-018-0081-5
Sun, J., and Du, W. P. (2017). Effects of precipitation and temperature on net primary productivity and precipitation use efficiency across China’s grasslands. GIsci Remote Sens. 54, 881–897. doi: 10.1080/15481603.2017.1351147
Sun, Y. F., Yang, Y. H., Zhao, X., Tang, Z. Y., Wang, S. P., and Fang, J. Y. (2021). Global patterns and climatic drivers of above-and belowground net primary productivity in grasslands. Sci. China Life Sci. 64, 739–751. doi: 10.1007/s11427-020-1837-9
Vázquez, D. P., Melián, C. J., Williams, N. M., Blüthgen, N., Krasnov, B. R., and Poulin, R. (2007). Species abundance and asymmetric interaction strength in ecological networks. Oikos 116, 1120–1127. doi: 10.1111/j.0030-1299.2007.15828.x
Wang, C., Liu, D. W., and Bai, E. (2018). Decreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition. Soil Biol. Biochem. 120, 126–133. doi: 10.1016/j.soilbio.2018.02.003
Wang, R. L., Wang, Q. F., Zhao, N., Yu, G. R., and He, N. P. (2017b). Complex trait relationships between leaves and absorptive roots: Coordination in tissue N concentration but divergence in morphology. Ecol. Evol. 7, 2697–2705. doi: 10.1002/ece3.2895
Wang, J., Dong, J., Yi, Y., Lu, G., Oyler, J., and Smith, W. K. (2017a). Decreasing net primary production due to drought and slight decreases in solar radiation in China from 2000 to 2012. J. Geophys. Res. Biogeosci. 122, 261–278. doi: 10.1002/2016JG003417
Wang, X., Xu, Z. E., Lü, X. T., Wang, R. Z., Cai, J. P., and Jiang, Y. (2017c). Responses of litter decomposition and nutrient release rate to water and nitrogen addition differed among three plant species dominated in a semi-arid grassland. Plant Soil. 418, 241–253. doi: 10.1007/s11104-017-3288-8
Wang, N., Quesada, B., Xia, L. L., Butterbach-Bahl, K., Goodale, C. L., and Kiese, R. (2019). Effects of climate warming on carbon fluxes in grasslands—A global meta-analysis. Glob. Chang. Biol. 25, 1839–1851. doi: 10.1111/gcb.14603
Wei, D., Zhao, H., Zhang, J. X., Qi, Y. H., and Wang, X. D. (2020). Human activities alter response of alpine grasslands on Tibetan Plateau to climate change. J. Environ. Manag. 262:110335. doi: 10.1016/j.jenvman.2020.110335
Wieder, W. R., Cleveland, C. C., Smith, W. K., and Todd-Brown, K. (2015). Future productivity and carbon storage limited by terrestrial nutrient availability. Nat. Geosci. 8, 441–444. doi: 10.1038/ngeo2413
Xu, Z. W., Li, M. H., Zimmermann, N. E., Li, S. P., Li, H., and Ren, H. Y. (2018). Plant functional diversity modulates global environmental change effects on grassland productivity. J. Ecol. 106, 1941–1951. doi: 10.1111/1365-2745.12951
Yang, Y. T., Chen, Z. F., Xu, B. C., Wei, J. Q., Zhu, X. X., and Yao, H. B. (2022). Using Trait-Based Methods to Study the Response of Grassland to Fertilization in the Grassland in Semiarid Areas in the Loess Plateau of China. Plants 11:2045. doi: 10.3390/plants11152045
Yang, Y., Wang, Z. Q., Li, J. L., Gang, C. C., Zhang, Y. Z., and Zhang, Y. (2016). Comparative assessment of grassland degradation dynamics in response to climate variation and human activities in China, Mongolia, Pakistan and Uzbekistan from 2000 to 2013. J. Arid Environ. 135, 164–172. doi: 10.1016/j.jaridenv.2016.09.004
Yu, H. L., Ding, Q. N., Meng, B. P., Lv, Y. Y., Liu, C., and Zhang, X. Y. (2021). The Relative Contributions of Climate and Grazing on the Dynamics of Grassland NPP and PUE on the Qinghai-Tibet Plateau. Remote Sens. 13:3424. doi: 10.3390/rs13173424
Zhang, D. Y., Peng, Y. F., Li, F., Yang, G. B., Wang, J., and Yu, J. C. (2019). Trait identity and functional diversity co-drive response of ecosystem productivity to nitrogen enrichment. J. Ecol. 107, 2402–2414. doi: 10.1111/1365-2745.13184
Zhang, Q. Y., Wang, Q. F., Zhu, J. X., Xu, L., Li, M. X., and Rengel, Z. (2021). Higher soil acidification risk in southeastern Tibetan Plateau. Sci. Total Environ. 755:143372. doi: 10.1016/j.scitotenv.2020.143372
Zhang, Z. C., Sun, J., Liu, M., Xu, M., Wang, Y., and Wei, T. X. (2020). Don’t judge toxic weeds on whether they are native but on their ecological effects. Ecol. Evol. 10, 9014–9025. doi: 10.1002/ece3.6609
Zhao, Y. F., Wang, X., Jiang, S. L., Xiao, J. J., Li, J., and Zhou, X. H. (2022). Soil development mediates precipitation control on plant productivity and diversity in alpine grasslands. Geoderma 412, 115721. doi: 10.1016/j.geoderma.2022.115721
Zheng, S. X., Chi, Y. G., Yang, X. J., Li, W. H., Lan, Z. C., and Bai, Y. F. (2022). Direct and indirect effects of nitrogen enrichment and grazing on grassland productivity through intraspecific trait variability. J. Appl. Ecol. 59, 598–610. doi: 10.1111/1365-2664.14078
Zheng, Z. T., Zhu, W. Q., and Zhang, Y. J. (2020). Seasonally and spatially varied controls of climatic factors on net primary productivity in alpine grasslands on the Tibetan Plateau. Glob. Ecol. Conserv. 21:e00814. doi: 10.1016/j.gecco.2019.e00814
Zhou, W., Yang, H., Huang, L., Chen, C., Lin, X. S., and Hu, Z. J. (2020). Grassland degradation remote sensing monitoring and driving factors quantitative assessment in China from 1982 to 2010. Ecol. Indic. 83, 303–313. doi: 10.1016/j.ecolind.2017.08.019
Keywords: grassland productivity, leaf traits, climate change, soil nutrients, desert steppe
Citation: Wang X, Wang R and Gao J (2022) Precipitation and soil nutrients determine the spatial variability of grassland productivity at large scales in China. Front. Plant Sci. 13:996313. doi: 10.3389/fpls.2022.996313
Received: 17 July 2022; Accepted: 17 August 2022;
Published: 09 September 2022.
Edited by:
Peijian Shi, Nanjing Forestry University, ChinaReviewed by:
Xiangjin Shen, Northeast Institute of Geography and Agroecology (CAS), ChinaRomà Ogaya, Ecological and Forestry Applications Research Center (CREAF), Spain
Zhuoya Zhang, Southwest Forestry University, China
Yuyang Zhang, Tarim University, China
Copyright © 2022 Wang, Wang and Gao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jie Gao, jiegao72@gmail.com