Lian-Lian Fan1,2,3
Okhonniyozov Mekrovar1,2,3
Yao-Ming Li1,2,3
Kai-Hui Li1,2,3Xue-Xi Ma1,2
Jie-Fei Mao1,2,3*
- 1State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi, China
- 2Research Center for Ecology and Environment of Central Asia, Chinese Academy of Sciences, Urumqi, China
- 3University of Chinese Academy of Sciences, Beijing, China
Grasslands provide key resource for the millions of people who are dependent on livestock in Tajikistan. Productivity and species richness (SR) are important characteristics of grassland ecosystems and are greatly affected by nutrient inputs. The effect that climate change might have on these characteristics remains unclear. Here, an in situ nitrogen (N) and phosphorus (P) fertilization experiment was conducted at four sites along with an elevational gradient (650, 1,100, 1,250, and 2,000 m) in western Tajikistan over 2 years (2018 and 2019) to examine the influences of nutrient availability and climate change on aboveground biomass (AGB) and SR; precipitation and temperature were also considered to analyze the responses. It demonstrated that enrichment with N, P, and their combinations significantly increased AGB along with an elevational gradient (p < 0.05). AGB increased as the concentrations of nutrients added increased. The maximum AGB, which was 2-fold higher compared with control, was observed when 90 kg N ha–1year–1 and 30 kg P ha–1year–1 were added. In addition, nitrogen addition alone stimulated greater AGB than P addition, although no significant difference was observed between these two treatments. Enrichment with N, P, and their combination had no significant effect on SR; however, SR significantly changed at different elevation. Elevation had direct effect on precipitation and temperature, which, in turn, resulted in variation in AGB and SR. Moreover, both nutrient and elevation had significant effect on AGB and SR, but there was no interaction effect of them. AGB and SR interacted with significant negative correlation. In the high-elevation area, plants grew better in the warmer year (2018); this indicates that grasslands in high mountain areas in Tajikistan might have higher productivity as the climate warms, which will positively affect the economic development of the country.
Introduction
Grasslands are important terrestrial ecosystems, accounting for almost 30% of the land area of earth (Bi et al., 2020). They are important providers of various ecosystem services including the production of livestock forage and biodiversity conservation. Currently, grassland ecosystems have experienced severe changes, for example, the decline of grassland productivity and the loss of species diversity, as both climate change and anthropogenic activities have contributed to dynamic variations in these ecosystems (Jiang et al., 2015; Su et al., 2021). Plant productivity and species richness (SR) are both fundamental properties of grassland and those two properties can reflect the structure and function of grassland ecosystems (Tilman et al., 2006; Wu et al., 2020). However, the degree to which these two properties might change under climate change remains unclear.
Nitrogen (N) and phosphorus (P) play vital roles in plant development and are the most important limiting nutrients in various grasslands (Fay et al., 2015). Climate change and human activities have greatly exacerbated grassland degradation. Fertilization is the primary means for improving the productivity of grassland systems (Ineichen et al., 2020). The effect of nitrogen and phosphorus application on the productivity and diversity of grassland ecosystems has been examined by previous studies and the results of these studies are variable (Alhamad et al., 2012; Harpole et al., 2017; Soons et al., 2017). Increase in nitrogen is known to significantly increase the above-ground biomass (AGB) and reduce the species diversity (Zhao et al., 2019; Gao et al., 2021). However, some studies have shown that although nitrogen addition reduced the species diversity, it increased the functional richness (Niu et al., 2014; Mao et al., 2017). Enrichment with phosphorus fertilizer has been shown to have either increased or no effect on ecosystem productivity and species diversity (Menge and Field, 2007; Gao et al., 2021). The effect of phosphorus fertilizer on ecosystems can also be regulated by nitrogen fertilizer (Bracken et al., 2015).
Tajikistan is a mountainous country located in the region of Central Asia, where the grasslands are widely distributed. The main grassland distribution ranges about from 300 to 4,000 m, resulting in different microclimate and grassland types (Zhang et al., 2014). Grasslands are an important resource supporting the economy of Tajikistan. However, the rapid increase in the human population during the last three decades has greatly increased the demand for dairy products and meat; consequently, the overgrazing and overuse of grasslands have resulted in serious grassland degradation (Han et al., 2016). As fertilization is commonly used to improve grassland quality (Ineichen et al., 2020), understanding how fertilization affects grassland productivity and diversity is essential for managing grassland and ensuring the sustainable development and utilization of grassland resources. Meanwhile, the vegetation in this arid area is sensitive to variation in climate, as changes in vegetation in this region are correlated with annual precipitation and temperature fluctuation (Yin et al., 2016). Therefore, in this study, we explored the impacts of nutrient (N and P) addition along with an elevational gradient in western Tajikistan on aboveground biomass and SR in 2018 and 2019. Furthermore, precipitation and temperature were also considered to analyze the responses.
Materials and Methods
Study Area Description
This study was conducted in Tajikistan (total area: 1.431 × 103 km2), which covers a latitudinal range of 36°40′ to 41°04′N and a longitudinal range of 67°22′ to 75°13′E (Figure 1). Tajikistan is known as “the country of mountains” and its mountainous area accounts for about 93% of the total land area (Hafiz, 2008). The area of natural grassland in Tajikistan is 3.6895 million hm2 and the distribution of plant communities shows obvious zonal characteristics (Zhang et al., 2014). Depending on vegetation composition, altitude, and utilization, Tajikistan has more than 10 pasture types, which are important resources for the economic development (Ahmadov, 2012). This study was conducted at four sites varying in elevation. The study area experiences a continental climate. The patterns of temperature change at the four sites are similar (Figure 2). The annual average temperature increases from January, peaks in July, and then decreases. The distribution of precipitation is similar among the four experimental sites (Figures 2E–H). Most precipitation is concentrated in winter and spring; summer and autumn are the dry seasons. In Table 1, the annual total precipitation increased with elevation: 172.6 mm in Tabakchi (650 m) in 2018 compared with 1,290.00 and 692.9 mm in Luchob (1,250 m) and Ziddi (2,000 m), respectively.
Figure 1. Map of the study area and the experimental design. N0P0–the control; N3—addition of 30 kg N ha−1year−1, P3—addition of 30 kg P ha−1year−1, N3P3—addition of 30 kg N ha−1year−1 and 30 kg P ha−1year−1, N9P3—addition of 90 kg N ha−1year−1 and 30 kg P ha−1year−1. The fertilizer was added via a single application. The Shuttle Radar Topographic Mission (SRTM) 30 m digital elevation model (DEM) dataset was obtained from USGS (https://glovis.usgs.gov/).
Figure 2. The average monthly temperature (A–G) and monthly precipitation (E–H) in the four experimental sites.
Table 1. The location of the four experimental sites and the total precipitation and annual average temperature among the four sites.
Experimental Design
Our in situ field experiments were conducted in West Tajikistan from 2018 to 2019. Experiments were conducted at four sites (Figure 1) along with an elevational gradient: Tabakchi (650 m), Blakhi (1,110 m), Luchob (1,250 m), and Ziddi site (2,000 m) to examine the effect of nutrient (N and P) addition on the productivity and the SR of pasture ecosystems. At each site, four blocks (20 m × 3 m) were setup at each site with a 1-m wide buffer zone between the blocks. Each block was further divided into 20 plots for different nutrition treatment (Figure 1). The size of each plot was 3 m × 3 m with 1 m wide buffer zones between the plots. Five nutrition treatments were conducted in each block and the four replicates of each five treatments were randomly assigned: CK—control; N3—addition 30 kg Nha–1 year–1, P3—addition 30 kg Pha–1 year–1, N3P3—addition 30 kg Nha–1 year–1 and 30 kg Pha–1 year–1, and N9P3—addition 90 kg Nha–1 year–1 and 30 kg Pha–1 year–1. Nutrient addition was added via a single application after snow melting. Stainless steel sheets were used to prevent water infiltration between treatments, 40 cm of the sheets were inserted into the ground and 10 cm of the sheets were left above ground. The meteorological data were obtained from the metrological station near the experimental sites.
Plant Sampling
In 2018 and 2019, the name and number of plant species were recorded in each 1 m × 1 m quadrats under different nutrient treatments along with an elevational gradient each year during the peak of the growing season. Above-ground plants were harvested at the end of the growth period in 1 m × 1 m quadrats. AGB was determined by oven drying at 68°C until a constant weight was achieved. Total SR was recorded by noting all the plant species present within 1 m × 1 m quadrats.
Statistical Analysis
A linear mixed-effects models were used to analyze the effect of different nutrition addition treatments on plant AGB and SR along with an elevational gradient, with nutrient and elevation as independent variable and years as a random factor. Structural equation modeling (SEM) was used to evaluate and quantify the effects of elevation, precipitation, temperature, and nutrient addition on plant AGB. The SEM was fit using the “lavaan” package in R version 3.3.2 (R Core Team, 2016). The model was considered to be a good fit based on the following criteria: an insignificant (p > 0.05) chi-square test statistic, chi-square/df < 3, root mean square error of approximation (RMSEA) < 0.05, goodness of fit index (GFI) > 0.95, and standardized root mean square residual (SRMR) < 0.09 (Liu et al., 2020). The Tukey’s post hoc test following the one-way ANOVA was used to compare the differences of mean values at p < 0.05. All the statistical analyses were conducted in R version 3.3.2 (R Core Team, 2016). Origin 2015 software (Origin Lab, Northampton, MA, United States) was used to prepare figures.
Results
Above-Ground Biomass and Species Richness Under Different Nutrition Addition Treatment Along With an Elevational Gradient
In both 2018 and 2019, the addition of nutrients greatly increased AGB as the elevation increased (Figures 3A,B and Table 2; p < 0.05). AGB was higher under N addition than under P addition; however, no significant differences were observed between these treatments. AGB was higher under nutrient addition compared with the control and the maximum AGB was observed under N9P3. Nutrient addition significantly increased AGB 2-fold at Tabakchi and Ziddi relative to the control; AGB at these sites ranged from 37.65 to 99.50 g/m2 and from 451.30 to 918.86 g/m2, respectively. The effect of nutrient addition on species diversity did not significantly differ at these two sites. However, elevation had a significant effect on species diversity (Figures 3C,D and Table 2). SR was lowest at Luchob (6 species/m2) (Figures 3C,D). There was no interaction effect between nutrients and elevations and they independently affected AGB and SR (Table 2).
Figure 3. Changes in the above-ground biomass (A,B) and species richness (C,D) changing under different nutrient treatments at different elevations in 2018 (left) and 2019 (right). Within each panel, groups with the same capital letters above them are not significantly different at p < 0.05 (Duncan’s test). Within each group of five bars, bars sharing the different letters are significantly different from each other at p < 0.05. Letters are not presented when none of the five bars differs significantly from any of the others (Duncan’s test), ns means that differences are not significant. The data were showed as mean ± SE.
Table 2. Summary of the linear mixed-effects models relating fixed factor (N, E) and random factor (Y) for the grassland above-ground biomass and species richness.
Influence of Elevation P and T on Aboveground Biomass and Species Richness
The SEM analysis showed that elevation had significant (p < 0.01) positive direct and indirect effects on plant AGB and SR. Elevation had direct effects on annual total precipitation and annual average temperature, but positive direct effect on precipitation and negative direct effect on temperature. Annual total precipitation and annual average temperature were significant negatively correlated (p < 0.01). Nutrient addition had positive direct effect on plant AGB, but nutrient addition had no significant direct effect on SR (p > 0.05). AGB and SR had negative correlation (Figure 4). AGB was observed higher in 2018 (more precipitation) at lower and middle elevation areas (Figures 5A,B). Contrary to expectations, it showed the opposite result at the high elevation (Ziddi). In the warm year (2018), AGB at higher elevation was more than in the cold year (2019) though with less precipitation (Figures 5A,B). SR had no significant between years at Tabakchi and Blakhi. However, SR had significant difference between years at Luchob and Ziddi with higher SR in 2018 (Figure 5C).
Figure 4. Structural equation model for the effect of environmental factors and nutrient inputs on the plant above-ground biomass. Red and black arrows indicate positive and negative relationships, respectively. Single-headed arrows represent causal relationships and two-headed arrows indicate correlation. Dotted arrows represent non-significant paths (p > 0.05). Values adjacent to arrows indicate standardized path coefficients. The proportion of variance explained (r2) appears alongside response variable in the model. Significance levels are denoted with **p < 0.01.
Figure 5. Annual total precipitation and annual average temperature at different elevations (A), the average above-ground biomass (B), and species richness (C) in 2018 and 2019. Data of the control treatment were used (mean ± SE). Different capital letters indicate significant difference between years at the same elevation (p < 0.05), ns, no significant difference.
Discussion
Effect of Nutrients on the Grass Aboveground Biomass and Species Richness
Nitrogen and phosphorus are two limiting nutrients that regulate the development of grassland ecosystems (Fay et al., 2015). Nitrogen and phosphorus addition significantly increased the AGB as the elevation increased and AGB increased as the concentrations of nutrients applied to the soil increased. These results are consistent with previous studies in grassland ecosystems showing that nitrogen and phosphorus additions stimulated primary production (Li et al., 2015; Zhao et al., 2019). Moreover, our findings indicated that nutrient addition was more effective at promoting plant growth compared to the control. The effect was enhanced by the increased precipitation at the middle to lower elevations (Tabqichi, Blakhi, and Luchob). The AGB was higher in 2019 than in 2018 in the nutrient addition treatments because precipitation was higher in 2019 compared with 2018 (Figures 3A,B and Table 1). Other studies have shown that increase in water promotes nutrient availability and increases primary productivity (Yang et al., 2011; Fan et al., 2013). A global scale experiment of nutrient addition in grasslands has shown that N addition increases biomass to a greater degree than P addition. The addition of both N and P was more beneficial for biomass increasing than adding only one of these two elements (Fay et al., 2015). It is in an agreement with our 2-year nutrient fertilizer experiment. The amount of the AGB with adding only nitrogen was higher than that adding only phosphorus. Combined NP enrichment was more beneficial for the plant growth than either along nitrogen or phosphorus addition. Furthermore, the AGB showed a positive response to nutrient concentration, especially to the nitrogen. These results indicated that nitrogen is more likely to be the limiting nutrient for biomass production in Tajikistan grassland as the elevation increases and the effect of phosphorus addition on the plant growth was regulated by nitrogen nutrient.
Species diversity plays an important role in maintaining the structure of ecosystems (Hautier et al., 2020). Increase in nitrogen is thought to significantly increase the AGB and reduces the species diversity (Humbert et al., 2016). Several meta-analyses have shown that nitrogen enrichment generally reduced plant SR in terrestrial ecosystems, especially in grassland ecosystems (Soons et al., 2017; Yue et al., 2020). However, the results of this study were not fully consistent with the results of these previous studies. Nitrogen, phosphorus, or their combinations addition did not increase or decrease SR, but SR exhibited significant responses within elevations. Among with an elevational gradient, SR showed significant different with Lucholb that had the lowest SR. The reason for that might be explained by short-term nutrient addition. Long-term experiments may be necessary for characterizing the effect of nutrient addition on species diversity. The previous study found that the greater increase in plant productivity with long-term nitrogen addition, the greater decline in plant diversity, and vice versa (Yue et al., 2020). This study conducted nutrient addition for 2 years; the SR to nutrient availability may be not so sensitive, which needs to take longer time fertilizer experiment to analyze the response. Overall, nutrient limitation often increases with increasing elevation (Johnson et al., 2000), but some recent observational studies suggested that the relative importance of different limiting nutrients changes with elevation (Fisher et al., 2013; Sundqvist et al., 2014). In this study, though both nutrient and elevation had significant effect on the AGB or SR, there was no interaction effect of them on the AGB or SR and they independently affected AGB and SR (Table 2).
Changes in Biomass and Diversity Driven by Environmental Variables
Precipitation and temperature are considered to be the major limiting factors in arid and semiarid region and high latitude ecosystems and these systems are, therefore, predicted to be particularly sensitive to climate change (IPCC, 2014). Most previous studies have found productivity to be linearly correlated with mean annual precipitation (Hu et al., 2007; Mowll et al., 2015), as the amount of precipitation largely determines the productivity of the grassland ecosystem (Esteban and Osvaldo, 2000; Bai et al., 2004). Other studies indicated that declining temperatures associated with increasing elevation can directly influence ecosystem properties (Bahram et al., 2012; Sundqvist et al., 2014). Our study area spans different elevations, so there are differences in rainfall and temperature. Through the SEM analysis, we can know that annual total precipitation and annual average temperature variation were directly significantly affected via elevational gradient, with positive on precipitation and negative on temperature. Meantime, precipitation and temperature were negatively correlated (Figure 4). Therefore, consistent with previous studies, the AGB and SR were affected not only by elevation, but also by annual total precipitation and annual average temperature (Figure 5). Different effects of precipitation and temperature can be observed at different elevations correlated with the vegetation. In this study, vegetation changed from desert grassland to mountain grassland from the south to the north experimental site and correspondingly, the AGB increased. Specifically, the AGB in desert grassland was more sensitive to variation of precipitation, which can be more than quadrupled in 2019 (the year with more precipitation) compared with 2018. Besides precipitation, annual average temperature also had a significant positive relationship with the grassland AGB. Generally, warming had a negative indirect effect on plant productivity because of increased soil drought associated with higher evapotranspiration and temperature (Dulamsuren et al., 2013). Nevertheless, in our high-elevation experimental site Ziddi, the year (2018) with higher temperature had the higher AGB compared with the year with relative lower temperature (2019). This result indicated that warming might be beneficial for the plant growth. Since global warming is an indisputable fact (IPCC, 2014), the grassland in high mountain areas in Tajikistan may have higher productivity in the future based on current warming trends and this would have a positive effect on the economic development of Tajikistan.
Productivity and diversity are two major properties of ecosystems and environmental factors can affect the productivity; the variation will lead to the change of diversity simultaneously and vice versa (Gross and Cardinale, 2007; Wang et al., 2019). In arid areas, water largely determines the plant growth, thus competition among plants is mainly access to water resources (Mu et al., 2021). When precipitation increases, the rapid growth of some species will further lead to the competition for light, which may lead to the decrease in diversity. Therefore, the relationship between growth ability and species diversity is mainly caused by the competition of resources (Miki, 2009). The competition for water and light, thus, provides the best explanation of our finding.
Conclusion
The AGB and SR in this mountainous arid community were largely affected by nutrient availability and elevation. Nutrient enrichment was beneficial for the plant growth. The AGB increased with the more nutrient concentration along with an elevational gradients and the maximum increment of the AGB can be quadrupled compared with CK. Nitrogen adding alone was more beneficial for accumulation of the AGB compared with phosphorus adding. Enrichment with N, P, and their combination had no significant effect on SR; however, SR significantly changed at different elevation. Moreover, both nutrient and elevation had significant effect on the AGB and SR, but there was no interaction effect of them. The AGB and SR interacted with significant negative correlation. Elevation had direct effect on precipitation and temperature, which, in turn, caused variation in the AGB and SR.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author Contributions
Y-ML and J-FM conceived the study. L-LF conducted the experiments, analyzed the data, and drafted the manuscript. OM, K-HL, and X-XM performed the experiments. All authors showed in the article contributed to study conceptualization, read, and approved the final version of this manuscript.
Funding
This study was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDA20020101), the “West Light” Talent Program of the Chinese Academy of Sciences (2019-XBQNXZ-B-005), and National Natural Science Foundation of China (42077327). Strategic Priority Research Program of the Chinese Academy of Sciences (XDA2005010304).
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.
Acknowledgments
This study was conducted in Tajikistan and we thank all the staffs at the Research Center for Ecology and Environment of Central Asia (Dushanbe subcenter) for their great and effective assistance. A special thanks to A. A. Madaminov for his field work and comments. We also thank TopEdit (www.topeditsci.com) for linguistic assistance during the preparation of this manuscript.
References
Ahmadov, H. M. (2012). Influence of over-grazing on biodiversity and productivity of the pasture on high mountain zone Tajikistan. J. Acad. Soil Sci. 2, 33–42, (In Russia).
Alhamad, N. M., Alrababah, A. M., and Gharaibeh, A. M. (2012). Impact of burning and fertilization on dry Mediterranean grassland productivity and diversity. Acta Oecol. 40, 19–26. doi: 10.1016/j.actao.2012.02.005
Bahram, M., Põlme, S., Kõljalg, U., Zarre, S., and Tedersoo, L. (2012). Regional and local patterns of ectomycorrhizal fungal diversity and community structure along an altitudinal gradient in the Hyrcanian forests of northern Iran. New Phytol. 193, 465–473. doi: 10.1111/j.1469-8137.2011.03927.x
Bai, Y. F., Han, X., Wu, J., Chen, Z., and Li, L. (2004). Ecosystem stability and compensatory effects in the Inner Mongolia grassland. Nature 431, 181–184. doi: 10.1038/nature02850
Bi, X., Li, B., Zhang, L. X., and Nan, B. (2020). Response of grassland productivity to climate change and anthropogenic activities in arid regions of Central Asia. PeerJ 8:e9797. doi: 10.7717/peerj.9797
Bracken, M. E. S., Hillebrand, H., Borer, E. T., Seabloom, E. W., Cebrian, J., Cleland, E. E., et al. (2015). Signatures of nutrient limitation and co-limitation: responses of autotroph internal nutrient concentrations to nitrogen and phosphorus additions. Oikos 124, 113–121. doi: 10.1111/oik.01215
Dulamsuren, C., Wommelsdorf, T., Zhao, F. J., Xue, Y. Q., Zhumadilov, B. Z., Leuschner, C., et al. (2013). Increased summer temperatures reduce the growth and regeneration of Larix sibirica in southern boreal forests of Eastern Kazakhstan. Ecosystems 16, 1536–1549. doi: 10.1007/s10021-013-9700-1
Esteban, G. J., and Osvaldo, E. S. (2000). Controls of Grass and shrub aboveground production in the patagonian steppe. Ecol. Appl. 10, 541–549. doi: 10.2307/2641113
Fan, L. L., Li, Y., Tang, L. S., and Ma, J. (2013). Combined effects of snow depth and nitrogen addition on ephemeral growth at the southern edge of the Gurbantunggut Desert, China. J. Arid Land. 5, 500–510. doi: 10.1007/s40333-013-0185-8
Fay, P. A., Prober, S., Harpole, W. S., Knops, J. M. H., Bakker, J. D., Borer, E. T., et al. (2015). Grassland productivity limited by multiple nutrients. Nat. Plants 1:15080. doi: 10.1038/nplants.2015.80
Fisher, J. B., Malhi, Y., Torres, I. C., Metcalfe, D. B., van de Weg, M. J., Meir, P., et al. (2013). Nutrient limitation in rainforests and cloud forests along a 3,000-m elevation gradient in the Peruvian Andes. Oecologia 172, 889–902. doi: 10.1007/s00442-012-2522-6
Gao, L. M., Su, J., Chen, C., Tian, Q., and Shen, Y. X. (2021). Increases in forage legume biomass as a response to nitrogen input depend on temperature, soil characters and planting system: a meta-analysis. Grass For. Sci. 76, 309–319. doi: 10.1111/gfs.12511
Gross, K., and Cardinale, B. (2007). Does species richness drive community production or vice versa? reconciling historical and contemporary paradigms in competitive communities. Am. Natur. 170, 207–220. doi: 10.1086/518950
Hafiz, M. (2008). State of plant genetic resources for food and agriculture (pgrfa) in the Republic of Tajikistan. Count. Rep. 8, 4–10.
Han, Q. F., Luo, G. P., Li, C. F., Shakir, A., Wu, M., and Saidov, A. (2016). Simulated grazing effects on carbon emission in Central Asia. Agricult. For. Meteorol. 216, 203–214. doi: 10.1016/j.agrformet.2015.10.007
Harpole, W. S., Sullivan, L. L., Lind, E. M., Firn, J., Adler, P. B., Borer, E. T., et al. (2017). Out of the shadows: multiple nutrient limitations drive relationships among biomass, light and plant diversity. Funct. Ecol. 31, 1839–1846. doi: 10.1111/1365-2435.12967
Hautier, Y., Zhang, P., Loreau, M., Wilcox, K. R., Seabloom, E. W., Borer, E. T., et al. (2020). General destabilizing effects of eutrophication on grassland productivity at multiple spatial scales. Nat. Commun. 11:1234567890. doi: 10.1038/s41467-020-19252-4
Hu, Z., Fan, J., Zhong, H., and Yu, G. (2007). Spatiotemporal dynamics of aboveground primary productivity along a precipitation gradient in Chinese temperate grassland. Sci. China Ser. D Earth Sci. 50, 754–764. doi: 10.1007/s11430-007-0010-3
Humbert, J., Dwyer, A., Andrey, A., and Arlettaz, A. (2016). Impacts of nitrogen addition on plant biodiversity in mountain grasslands depend on dose, application duration and climate: a systematic review. Glob. Change Biol. 22, 110–120. doi: 10.1111/gcb.12986
Ineichen, S., Marquardt, S., Kreuzer, M., and Reidy, B. (2020). Forage quality of species-rich mountain grasslands subjected to zero, PK and NPK mineral fertilization for decades. Grass Forage Sci. 75, 385–397. doi: 10.1111/gfs.12488
IPCC (2014). Climate Change 2014: Synthesis Report: Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Geneva: IPCC.
Jiang, W. G., Yuan, L. H., Wang, W. J., Cao, R., Zhang, Y. F., and Shen, W. M. (2015). Spatio-temporal analysis of vegetation variation in the yellow river basin. Ecol. Indicat. 51, 117–126. doi: 10.1016/j.ecolind.2014.07.031
Johnson, C. E., Driscoll, C. T., Siccama, T. G., and Likens, G. E. (2000). Element fluxes and landscape position in a northern hardwood forest watershed ecosystem. Ecosystems 3, 159–184. doi: 10.1007/s100210000017
Li, K. H., Liu, X. J., Song, L., Gong, Y. M., Lu, C., Yue, P., et al. (2015). Response of alpine grassland to elevated nitrogen deposition and water supply in China. Oecologia 177, 65–72. doi: 10.1007/s00442-014-3122-4
Liu, X., Lamb, E. G., and Zhang, S. (2020). Nitrogen addition impacts on soil microbial stoichiometry are driven by changes in plant resource stoichiometry not by the composition of main microbial groups in an alpine meadow. Biol. Fertil. Soils 56:1423. doi: 10.1007/s00374-019-01423-1
Mao, W., Felton, A. J., and Zhang, T. H. (2017). Linking changes to intraspecific trait diversity to community functional diversity and biomass in response to snow and nitrogen addition within an inner mongolian grassland. Front. Plant Sci. 8:339. doi: 10.3389/fpls.2017.00339
Menge, D. N. L., and Field, C. B. (2007). Simulated global changes alter phosphorus demand in annual grassland. Glob. Change Biol. 13, 2582–2591. doi: 10.1111/j.1365-2486.2007.01456.x
Miki, T. (2009). A new graphical model for untangling complex relationships among environment, biodiversity, and ecosystem functioning. Ecol. Res. 24, 937–941. doi: 10.1007/s11284-008-0552-7
Mowll, W., Blumenthal, D. M., Cherwin, K., Smith, A., Symstad, A. J., Vermeire, L. T., et al. (2015). Climatic controls of aboveground net primary production in semi-arid grasslands along a latitudinal gradient portend low sensitivity to warming. Oecologia 177, 959–969. doi: 10.1007/s00442-015-3232-7
Mu, X. H., Huang, G., Li, Y., Zheng, X. J., Xu, G. Q., Wu, X., et al. (2021). Population dynamics and life history response to precipitation changes for a desert ephemeral plant with biseasonal germination. Front. Plant Sci. 12:625475. doi: 10.3389/fpls.2021.625475
Niu, K. C., Choler, P., Bello, F. D., Mirotchnickg, N., Du, J., and Sun, S. (2014). Fertilization decreases species diversity but increases functional diversity: a three-year experiment in a Tibetan alpine meadow. Agricult. Ecosyst. Environ. 182, 106–122. doi: 10.1016/j.agee.2013.07.015
R Core Team (2016). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing. Available online at: https://www.r-project.org/
Soons, M. B., Hefting, M. M., Dorland, E., Lamers, L. P. M., Versteeg, C., and Bobbink, R. (2017). Nitrogen effects on plant species richness in herbaceous communities are more widespread and stronger than those of phosphorus. Biol. Conserv. 212(Pt B), 390–397. doi: 10.1016/j.biocon.2016.12.006
Su, Y., Ma, X. F., Le, J. J., Li, K. H., Han, W. X., and Liu, X. J. (2021). Decoupling of nitrogen and phosphorus in dominant grass species in response to long-term nitrogen addition in an Alpine Grassland in Central Asia. Plant Ecol. 9, 1–14. doi: 10.1007/s11258-020-01103-3
Sundqvist, M. K., Liu, Z. F., Giesler, R., and Wardle, D. (2014). Plant and microbial responses to nitrogen and phosphorus addition across an elevational gradient in subarctic tundra. Ecology 95, 1819–1835. doi: 10.1890/13-0869.1
Tilman, D., Reich, P. B., and Knops, J. M. H. (2006). Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature 441, 629–632. doi: 10.1038/nature04742
Wang, Z., Chiarucci, A., and Arratia, J. F. (2019). Integrative models explain the relationships between species richness and productivity in plant communities. Sci. Rep. 9:13730. doi: 10.1038/s41598-019-50016-3
Wu, Q., Ren, H. Y., Wang, Z. W., Li, Z. G., Liu, Y. H., Wang, Z., et al. (2020). Additive negative effects of decadal warming and nitrogen addition on grassland community stability. J. Ecol. 108:13363. doi: 10.1111/1365-2745.13363
Yang, H., Li, Y., Wu, M., and Zhang, Z. (2011). Plant community responses to nitrogen addition and increased precipitation: the importance of water availability and species traits. Glob. Change Biol. 17, 2936–2944. doi: 10.1111/j.1365-2486.2011.02423.x
Yin, G., Hu, Z. Y., Chen, X., and Tashpolat, Y. (2016). Vegetation dynamics and its response to climate change in Central Asia. J. Arid Land 8, 375–388. doi: 10.1007/s40333-016-0043-6
Yue, K., Jarvie, S., Senior, A. M., Meerbeek, K. V., Peng, Y., Ni, X., et al. (2020). Changes in plant diversity and its relationship with productivity in response to nitrogen addition, warming and increased rainfall. Oikos 129, 939–952. doi: 10.1111/oik.07006
Zhang, L. P., Li, X. S., and Ding, A. Y. (2014). The status and conservation of grassland resources in Tajikistan. Grass Feed. Livestock 1, 26–28, (In Chinese).
Keywords: climate change, meteorological factor, nutrient addition, above-ground biomass (AGB), Central Asia
Citation: Fan L-L, Mekrovar O, Li Y-M, Li K-H, Ma X-X and Mao J-F (2021) Effect of Nutrient Addition on the Productivity and Species Richness of Grassland Along With an Elevational Gradient in Tajikistan. Front. Plant Sci. 12:765077. doi: 10.3389/fpls.2021.765077
Received: 26 August 2021; Accepted: 11 November 2021;
Published: 10 December 2021.
Edited by:
Huakun Zhou, Key Laboratory of Restoration Ecology in Cold Regions, Northwest Institute of Plateau Biology, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Jianshuang Wu, Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of Agricultural Sciences (CAAS), ChinaHui Fu, Hunan Agricultural University, China
Copyright © 2021 Fan, Mekrovar, Li, Li, Ma and Mao. 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-Fei Mao, mjf@ms.xjb.ac.cn