Abstract
Background and aims:
Nitrogen (N) enrichment usually weakens the stabilizing effect of biodiversity on productivity. However, previous studies focused on plant species richness and thus largely ignored the potential contributions of plant functional traits to stability, even though evidence is increasing that functional traits are stronger predictors than species richness of ecosystem functions.
Methods:
We conducted a common garden experiment manipulating plant species richness and N addition levels to quantify effects of N addition on relations between species richness and functional trait identity and diversity underpinning the ‘fast–slow’ economics spectrum and community stability.
Results:
Nitrogen addition had a minor effect on community stability but increased the positive effects of species richness on community stability. Increasing community stability was found in the species-rich communities dominated by fast species due to substantially increasing temporal mean productivity relative to its standard deviation. Furthermore, enhancement in ‘fast–slow’ functional diversity in species-rich communities dominated by fast species under N addition increased species asynchrony, resulting in a robust biodiversity–stability relationship under N addition the artificial grassland communities.
Conclusion:
The findings demonstrate mechanistic links between plant species richness, ‘fast–slow’ functional traits, and community stability under N addition, suggesting that dynamics of biodiversity–stability relations under global changes are the results of species-specific responses of ‘fast–slow’ traits on the plant economics spectrum.
1 Introduction
Nitrogen (N) deposition has increased considerably because of increases in anthropogenic activities, including increased combustion of fossil fuels and direct agricultural inputs (; ; ). The resulting eutrophication can weaken the stabilizing effect of biodiversity on primary productivity (; Zhang et al., 2016; ; ). Stability of primary productivity is usually calculated as the ratio of mean productivity to its temporal standard deviation (; ; Yan et al., 2021). The stabilizing effect of biodiversity is the result of two nonexclusive mechanisms: high species stability, characterized by stability of each species, and high species asynchrony, characterized by variations in interspecific responses to environmental fluctuations, with high species richness (; ; ).
Although most studies focus on plant species richness as a measure of plant diversity, plant functional trait identity and diversity along the ‘fast–slow’ economics spectrum are also important (; ; ). Variation in leaf traits occurs along a single plant economics axis that is divided into acquisitive species with relatively fast growth and resource uptake (hereafter ‘fast’ species) and conservative species with relatively slow growth and resource uptake (hereafter ‘slow’ species) (Wright et al., 2004; ). Fast species with high specific leaf area, leaf N concentration, and nitrogen/phosphorus (N/P) ratio and low leaf dry matter content could enhance community stability by increasing community temporal mean productivity, whereas slow species have opposite characteristics could enhance community stability by increasing community persistence (; ; van der Plas et al., 2020). Plant functional traits stabilize community productivity because growth-related ‘fast–slow’ functional traits strongly influence plant productivity and persistence in grasslands (; ). Therefore, ecosystem functions may depend more on ‘fast–slow’ functional traits of species than on plant species richness (; ; ).
Previous studies assessing N addition effects on community stability in grassland focused on plant species richness. In those studies, nitrogen addition usually reduces community stability by decreasing species richness-mediated species asynchrony and species stability (; Zhang et al., 2016; ; ). However, previous studies ignored potential contributions of plant functional traits to community stability (). There are multiple mechanisms by which functional traits can affect community stability. First, the species-rich communities dominated by species with greater acquisitive trait can maintain higher community productivity than those dominated by species with higher conservative traits (; ; ), which can result in increasing community stability because of increased temporal mean productivity. However, the species-rich communities dominated by species with high slow trait values (i.e., high leaf dry matter content) are generally more stable over time (; ; ; ), with relatively high community stability resulting from decreased variation of productivity from year to year. Moreover, high ‘fast-slow’ functional diversity (fast-slow FD), characterized by interspecific differences in resource utilization, can be correlated with high community stability because of overyielding (; Yan et al., 2021). Similarly, high fast-slow FD also indicates that species can fluctuate independently from one another in fluctuating environments in the species-rich communities (i.e., relatively high species asynchrony) (Vergnon et al., 2009; ; ; ).
Communities composed of species with different fast–slow traits respond divergently to N addition (; ; ). The mechanisms by which N addition affects community stability can be explained by differences in functional traits among plant species. Plants with higher acquisitive trait values responded sensitively to fluctuating environments (), and plants in communities may shift from species with slow leaf economics to species with fast traits under N addition (Zhou et al., 2018; Zhang et al., 2019; Zheng et al., 2022), resulting in decreased species stability. In addition, the effects of N addition on community stability are correlated with either increased or decreased species asynchrony depending on whether N addition increased or decreased fast–slow FD in species-rich communities (). Nitrogen addition can decrease fast–slow FD by reducing species richness, which leads to a decline in species asynchrony (). However, nitrogen addition can also increase fast–slow FD because of divergent responses of different fast–slow functional traits when N addition does not change species richness (), which leads to an increase in species asynchrony.
In this study, we used a common garden experiment to examine the effects of N addition on relations between plant species richness, ‘fast–slow’ functional trait identity and diversity, and community stability of productivity in assembled grassland communities. The common garden experiment manipulated plant species richness at five levels (1, 2, 4, 6, and 8 species) and N addition at two levels (0 and 6 g N m−2 year−1). The experiment was designed to answer the following two questions: (1) Does N addition affect relations between plant species richness, functional traits (‘fast–slow’ trait identity and diversity), and community stability (species stability and species asynchrony)? (2) What are the mechanisms by which N addition affects those relations (see Table 1)?
Table 1
| Pathway | Hypotheses and mechanisms | References |
|---|---|---|
| Fast–slow trait identity→species stability | 1. Species with high fast trait values increase species stability by increasing productivity of individual species over time; 2. Species with high slow trait values increase the stability of species productivity by decreasing variation in the productivity of individual species over time. | ; ; |
| Fast–slow trait diversity→species asynchrony | The possibility for asynchronous responses to fluctuating environments increases when the diversity of fast–slow traits is high in species-rich communities. | ; |
| Species stability→community stability | High species stability increases community stability because of high fast trait-induced high productivity or high slow trait-induced low variation in productivity over time. | |
| Species asynchrony→community stability | High species asynchrony in response to fluctuating environments increases community stability because decreases in fast–slow trait values of some species are compensated for by increases in those of others, thus buffering temporal variations in the community-wide metric of fast–slow traits. |
Underlying mechanisms for relations between plant species richness, functional traits, and community stability.
2 Materials and methods
2.1 Study site and experimental design
The study site was the National Grasses Variety Regional Experimental Station (40°10’45’’N, 116°26’13’’E, 50 m above sea level) on Xiaotang Mountain in Beijing, China. Mean annual temperature and precipitation is 11.8°C and 526 mm (2000–2018), respectively. In 2019, a common garden experiment was established within an 8 m × 12 m area in a split-plot experiment design with N addition and plant species richness as the main treatment factors. The experiment was composed of the following: two N levels (0 and 6 g N m−2 year−1) × four replicate blocks × 5 species richness levels per N treatment/replicate (four combinations per level of diversity with one, two, four, and six species and one with eight species, in total 17 combinations) = 136 pots (Figure S1). Each pot (30-cm diameter, 50-cm height) was filled with uniformly mixed soil (natural native grassland soil with stones and roots removed by sieving) and sand (3:1 soil: sand ratio), three holes (3-cm diameter) were set in the bottom to allow the pots to drain freely, and then buried in the ground without any shelter.
Eight plant species were used to build the artificial grassland communities (Wang et al., 2023). Seeds for each species were purchased commercially and germination rates were determined in a greenhouse (Figure S1A). To maintain the consistency of plant communities at the beginning of the experiment, sowing rate was determined based on the germination rate of each species in order to manipulate the species to have the same original proportion in the species-rich communities (Figure S1A). In each block, well-mixed seeds of the different combinations at different levels of species richness were distributed randomly among the 17 pots under ambient condition and N addition treatment, respectively. Community abundance was maintained at approximately 60 individuals in each pot by manual removal based on native grassland community abundance. Plants in pots were not re-sown after 2019 because the species were perennial. Weeds in pots were removed by hand once a month. Nitrogen was added as urea (12.86 g m−2 year-1, twice the amount of N deposition in Beijing) dissolved in water and applied by spraying on 1 May each year; the ambient conditions received an equal amount of N-free water.
2.2 Productivity and ‘fast–slow’ functional trait measurements
First, the abundance of each plant species in each pot was measured in 2019–2021. Second, plant height and leaf traits of three healthy individuals of each plant species in each pot were measured. Leaf traits of two leaves on individual plants were measured once at the peak of the growing season. Leaf traits were measured using standard methods according to . Spread leaves were scanned, and leaf mean width (LW) and area (LA) were determined using a Li−Cor 310 (Li−Cor Inc. USA). Leaves were weighed fresh, dried for 48 h at 65°C, and then weighed again to calculate leaf dry matter content (LDMC, g dry mass g−1 fresh mass). Specific leaf area (SLA, cm2 dry mass g−1) was calculated as the ratio of leaf area to its dry mass. Third, oven-dried leaf samples were ground, and leaf carbon (LCC, %) and N (LNC, %) concentration were determined with an elemental analyzer (Vario MACRO Cube. Germany). Leaf phosphorus concentration (LPC, %) was determined by inductively coupled plasma atomic emission spectroscopy (AutoAnalyzer 3. Germany). Finally, the leaf C/N, C/P, and N/P ratios were calculated. All the traits were measured in 2019 and 2020.
To measure community and individual productivity, the aboveground part of each pot was clipped in early September (peak of the growing season) from 2019 to 2021. All plants were clipped in each pot, sorted by species, oven−dried at 65°C to constant weight, and then community productivity was calculated. Average productivity from 2019 to 2021 was defined as temporal mean productivity. Because plants were not re-sown after the experiment began, clipping might have affected aboveground biomass in the following year and thus affected community stability (; Wang et al., 2021). Fortunately, the plants selected in our experiment were all perennial species, aboveground biomass was clipped in the same way in all pots, and plant species richness did not change significantly across the three years (Figure S2). Hence, calculating community stability and its abiotic factors based on community productivity was reliable in our analysis.
2.3 Quantifying functional trait identity and diversity
Eight traits associated with the ‘fast-slow’ continuum of plant resource captures and usage (; ) were selected to calculate the community-weighted mean (CWM) for each pot using the mean abundance-weighted values for each species. Then, the trait values were converted to standardized units (Z-scores). First, principal component analysis (PCA) was used to distinguish between fast and slow traits based on community-weighted mean traits using the “vegan” package in R software (Figure S3). The first PCA axis of CWM of functional traits (CWM fast-slow) captured 57.0% of the variation in the 8 traits and represented the fast–slow plant economics spectrum of communities, from those dominated by fast species with high plant height, SLA, leaf N, and N/P ratio and low LDMC and LA to those dominated by slow species with low plant height, SLA, leaf N, and N/P ratio and high LDMC and LA (Figure S3). Principal component analysis was also used to separate plant species into those with fast or with slow growth based on species-specific traits (Figure S4), it separated plant species into two slow species (Chrysanthemum maximum and Allium schoenoprasum) and six fast species (Poa annua, Carex breviculmis, Medicago sativa, Astragalus adsurgens, Penstemon campanulatus, Dianthus barbatus) based on relations with ‘fast–slow’ traits (Figures S1, S4), which indicated that the assembled grassland communities were fast growth and sensitive to N addition. Last, fast-slow FD was calculated as mean abundance-weighted distance of species to represent the range and distribution of traits in multidimensional niche space using the ‘FD’ package ().
2.4 Quantifying stability
Community stability of productivity was quantified as the ratio of community temporal mean productivity (μT) to its standard deviation (σT) in each pot over three years (2019–2021) as follows: Scom= μT/σT. The ratio was widely used in many other studies (; ). According to recent theory, community stability can be mathematically partitioned into species stability and species asynchrony (Wang et al., 2019). Species stability (Ssp) was calculated as the reciprocal of mean species variability (Wang et al., 2019): , where σi represents the temporal standard deviation of species i. Species asynchrony () was defined as the reciprocal of a community-wide metric of species synchrony (Wang et al., 2019): .
2.5 Statistical analyses
Linear mixed-effects models using the “lme” function (package “nlme”) () were applied to test the effects of N addition and plant species richness on community productivity, CWM fast-slow, and fast-slow CWM. Nitrogen addition and plant species richness were fixed effects, and year, block, and plant combinations were random effects. Linear mixed-effects models were also used to examine the effects of N addition and plant species richness on community stability, species stability, and species asynchrony. Nitrogen addition and plant species richness were set as fixed effects, and block and plant combinations were set as random effects. Linear mixed-effects models were also used to test the relationship between community stability and abiotic factors. Nitrogen addition were set as fixed effects, and block and plant combinations were set as random effects. In order to test the interaction of N addition and plant species richness, analysis of covariance (ANCOVA) was used to test the differences in slopes of linear regressions between the ambient conditions and N addition treatments in the ‘car’ package in R software (). Moreover, we compared the plant functional traits between the ambient conditions and N addition treatment in different species by paired-samples t test.
A multi-group structural equation model (SEMm) was used to examine direct and indirect relations between plant species richness, CWM fast-slow and fast-slow FD (temporal mean across two years), and community stability under different levels of N addition. The ‘piecewise-SEM’ package in R software was used to construct the SEMm (). The SEMm was also constructed to determine whether the stabilizing effects on productivity operated via the two components of the invariability measure, temporal mean and standard deviation of productivity under ambient condition and N addition treatment. In the second model, direct and indirect paths from species richness, CWM fast–slow, fast–slow FD, and species asynchrony were added to the temporal mean and standard deviation of productivity. The models were assessed by Fisher’s C statistic, Akaike information criterion (AIC), p value, and degrees of freedom (df) values. Last, models with the lowest Fisher’s C and AIC values and higher p value (p≥0.05) were selected. In addition, pathways showed one arrow in the final SEMm when the relationship showed no significant difference between the ambient condition and N addition treatment, but it divided into two arrows when the relationship under N addition treatment significantly different from that in ambient condition.
All statistical analyses and graph preparations were performed using R 3.2.2 (). Effects were considered significant at p ≤ 0.05 or marginally significant at 0.05< p ≤ 0.10.
3 Results
3.1 Productivity and stability response to nitrogen addition and plant species richness
Over the 3-year experimental period, nitrogen addition significantly increased community productivity, and species richness was positively associated with temporal mean community productivity, although the effects were largely independent of one another (Figures 1A, B; Tables S1, S2). Neither N addition nor species richness significantly affected standard deviation of productivity (Figure 1C; Table S1). Although N addition had minor effects on community stability (Table S1), species richness was positively associated with community stability (Figure 1D; Table S4). Nitrogen addition strengthened the relation between species richness and community stability and species asynchrony, resulting in a significant N addition × species richness interaction term (Figures 1D, 2A; Tables S1, S2). In addition, nitrogen addition and species richness did not affect species stability (Figures 2A, B; Table S1). However, nitrogen addition marginally strengthened the relation between species richness and species asynchrony, resulting in a significant N addition × species richness interaction term (Figure 2A; Tables S1, S2), indicating the positive effect of species richness on species asynchrony enhanced under N addition in the species-rich communities.
Figure 1
Figure 2
3.2 Functional traits response to nitrogen addition and plant species richness
Nitrogen addition had a minor effect on fast–slow FD, but species richness was positively associated with fast–slow FD (Figure 2C; Tables S1, S5). Nitrogen addition strengthened the relation between species richness and fast–slow FD, resulting in a significant N addition × species richness interaction term (Figure 2C; Tables S1, S2), indicating the positive effect of species richness on fast–slow FD increased under N addition in species-rich communities. Nitrogen addition did not affect the CWM fast–slow (Figure 2D; Table S1), species richness was positively associated with CWM fast–slow (Figure 2D; Tables S1, S5), but the effects of N addition were not interacted with plant species richness (Figure 2D; Table S1). Moreover, nitrogen addition did not affect most CWM trait values, but CWM of LA and LW increased and CWM of SLA decreased with N addition (Table S3). Species richness was positively associated with most CWM trait values (Table S3).
3.3 Ecological factors influencing community stability
Community stability was positively associated with species stability and CWM fast–slow (higher CWM fast–slow represents greater acquisitive trait) under ambient condition (Figures 3B, D; Table S6). Under N addition, community stability was positively associated with species asynchrony (Figure 3A), species stability (Figure 3B), fast–slow FD (Figure 3C), and CWM fast–slow (Figure 3D; Table S6). Species asynchrony was positively associated with the fast–slow FD under N addition (Figure 4A; Table S7), and species stability was positively associated with the CWM fast–slow traits in both control and N addition (Figure 4B; Table S7). In addition, temporal mean productivity was positively associated with the CWM fast–slow and fast–slow FD in both ambient condition and N addition treatment (Figures S4A, C). It was positively associated with CWM of plant height, SLA, and leaf N content but negatively associated with CWM of LDMC (Figure S5). Moreover, the standard deviation of productivity exhibited no significant relationships with the CWM of plant height, LA, LDMC, SLA, or leaf N content or the N/P ratio (Figure S6).
Figure 3
Figure 4
The SEMm indicated that community stability–species stability, community stability–species asynchrony, and fast–slow FD–species richness relations were significantly affected by N addition (Figures 5, 6; Table S9). Specifically, both species stability and species asynchrony were key factors mediating the relation between species richness and community stability, with species stability driven by the CWM fast–slow and species asynchrony driven by fast–slow FD under N addition (Figure 5). Moreover, species richness increased temporal mean productivity by increasing the CWM fast–slow and decreased productivity variability, which resulted in a positive relation between biodiversity and stability under ambient condition (Figure 6; Table S9). However, species asynchrony–community stability, standard deviation of temporal mean productivity–community stability, and species richness–fast-slow FD relations were significantly affected by N addition (Figure S7B; Table S9). Species richness reduced the standard deviation of temporal mean productivity by increasing fast–slow FD, which resulted in a strengthened biodiversity–stability relation under N addition (Figure 6).
Figure 5
Figure 6
4 Discussion
The intrinsic links among species richness, functional traits (i.e., CWM fast–slow and fast–slow FD), and community stability (i.e., species stability and species asynchrony) were identified, and mechanisms underlying effects of N addition on the biodiversity–stability relation were clarified. In particular, species richness increased community stability mainly by increasing the CWM fast–slow under ambient condition. Although N addition did not significantly affect fast–slow FD and community stability, there was a positive interaction with species richness. Importantly, nitrogen addition strengthened the relation between species richness and community stability by increasing fast–slow FD and species asynchrony in the species-rich communities dominated by fast species. The results provide a novel perspective to understand potential mechanisms linking functional traits and community stability under global change.
4.1 Biodiversity increased community stability by increasing fast–slow trait identity
Species richness increased community stability by increasing species stability under ambient condition (Figure 5). The result is similar to theoretical predictions and experimental investigations that show a weak relation between biodiversity and stability when community stability is driven by dominant species (Tilman, 1996; ; Yang et al., 2012; Xu et al., 2015; Wang Y. et al., 2020). In addition, species stability was positively associated with CWM of plant height, leaf N content, and leaf N/P ratio (Figure S9). Thus, the species-rich communities dominated by species with relatively high fast trait values maintained relatively high species stability. Our results showed that species richness increased species stability by increasing the CWM fast–slow (Figure 5). The result contrasts with previous studies who showed that communities dominated by slow species exist stronger resistance response to fluctuating environments (; ; ). In further analysis, CWM of plant height, SLA, and leaf N content were positively associated but LDMC was negatively associated with temporal mean productivity rather than its standard deviation over time (Figures S5, S6). Therefore, the species-rich communities dominated by species with high acquisitive traits increased community stability by increasing temporal mean productivity (Figure S9). As with the effects of species richness on temporal mean productivity over time, those effects might be explained by selection effects in the species-rich communities (), and thus, the positive association between community stability and CWM fast–slow might be because more fast species were selected in the species-rich communities.
4.2 Nitrogen addition strengthened the biodiversity−stability relations
One key finding in our study was that the interaction between N addition and species richness positively affected community stability, with N addition increasing community stability in the species-rich communities dominated by fast species (Figures 1D, 5). Nitrogen addition strengthened the biodiversity−stability relation via a combination of two processes (Figure 5; Table S8). First, positive relations between species stability and community stability were not significantly affected by N addition. Second, the weak relation between community stability and species asynchrony under ambient condition became a significant positive relation after N addition. The results are not consistent with those of ; who found that N addition weakened the contribution of species richness to community stability. proposed that a relatively weak positive relation between species richness and species asynchrony was the result of decreasing species richness under N addition; whereas a relatively strong negative relation between species richness and species stability under N addition may be caused by shifts in functional traits in communities from species with conservative traits to those with acquisition traits (; ). There are two primary differences between our experiment and previous studies. First, the assembled communities in our experiment were dominated by fast growth species, nitrogen addition diversified species acquisitive traits and increased functional diversity in species-rich communities, which led to greater asynchronous dynamics to fluctuating environments and lower productivity variation over time in species-rich communities, resulting in a stronger relation between community stability and species asynchrony under N addition. Second, nitrogen addition did not significantly affect species richness and CWM fast–slow in our study (Table S1), resulting in no shifts in functional traits in species-rich communities, and thus and not changes in the relation between species richness and species stability (Figure 2A).
Asynchrony between species is generally attributed to divergent responses of species with different functional traits to environmental fluctuations (; ; ). Hence, species asynchrony could have increased because of divergent responses of species-specific fast–slow traits to N addition in our common garden experiment (Figure S10; Table S10) (; ). It has been proved that species with various fast–slow traits can stabilize the productivity of mixed-species communities by divergent responses to N addition (; ; ). For example, the legumes (Medicago sativa and Astragalus adsurgens) in the species-rich communities can perform symbiotic N fixation (). Thus, because growth is not limited by N, legume fast–slow traits might not be influenced by N addition (Figure S10). However, most grassland ecosystems are limited by N (; ; ), and growth of most species in grasslands (i.e., grasses and forbs) is closely associated with soil available nutrients (). To maximize nutrient use for plant growth (), acquisitive trait values (i.e., higher specific leaf area) of grasses, sedges, and forbs in the species-rich communities enhanced under N addition (Figure S10), possibly because of increases in soil available N and P (; Wang C. et al., 2020; ). Thus, the divergent responses of functional traits under N addition led to the increase in fast–slow FD in the species-rich communities.
Moreover, higher fast-slow FD in the species-rich communities under N addition indicated that an increase in dissimilarity in functional traits among plant species, higher differentiation in functional traits in species-rich communities exhibited greater asynchronous dynamics in fluctuating environments (; ; ). Hence, greater species asynchrony stabilized productivity because of its stronger negative effects on the standard deviation of productivity over time than its positive effects on temporal mean productivity under N addition (Figure 6; Table S9). Overall, the strengthened biodiversity–stability relation under N addition was caused by reduction in the standard deviation of productivity in fluctuating environments that was associated with the increase in fast-slow FD in the species-rich communities dominated by fast species.
4.3 Limitations
Furthermore, our study provides novel perspectives from which to understand potential links among biodiversity and community stability under N addition, but the conclusion in our common garden experiment differs from that in most previous studies was partly attributed to experimental duration, which could affect the response of community stability to N addition in two aspects. First, the artificial grassland communities were built through random combinations of different species, but some species in the mixed community were not native species (i.e., Chrysanthemum maximum). Hence, the artificial communities did not develop to s stable stage in 3 years, there were significant differences between the artificial and natural grasslands, which might have led to differences in results with previous studies. Hence, our results provide a new sight to investigate the relationship between plant functional traits and community temporal stability, but simply applying our conclusions to natural grasslands requires caution. Second, dynamic of climate conditions (i.e., mean annual temperature and precipitation) during 3 period years may not cover the climate variation in this site, which could also lead to the results in our experiment different from the others. Fortunately, each treatment in our experiment was conducted under the same climate condition, and our analysis focused on the comparison in the effects of biodiversity under ambient condition and nitrogen addition along the plant species richness gradients, which lead our conclusion different from the conclusion from other grassland communities but convinced.
5 Conclusions
The common garden experiment was an important step forward in understanding the biodiversity–stability relation because it demonstrated mechanistic links between species richness, functional traits (i.e., functional trait identity and diversity), and community stability (i.e., species stability and species asynchrony) under different levels of N addition. Specifically, a positive relation was identified between species stability and fast–slow trait identity, highlighting that species richness increased community stability because of high temporal mean productivity in the species-rich communities dominated by fast species. In particular, nitrogen addition strengthened the biodiversity–stability relation by increasing fast–slow functional diversity and species asynchrony in the species-rich communities dominated by fast species in the artificial grassland communities. The links suggest that N addition may increase grassland ecosystem stability when it does not alter species richness, and thus, appropriate nutrient addition may be an effective method to restore the degraded grasslands by using legume species.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
JS: Data curation, Writing – review & editing, Investigation, Resources. XLu: Formal Analysis, Writing – original draft, Data curation, Investigation. XLi: Writing – review & editing. YH: Writing – review & editing. CW: Formal Analysis, Writing – original draft, Methodology, Visualization, Conceptualization.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by Natural Science Foundation of Beijing Municipality (Grant No. 5232006), National Natural Science Foundation of China (Grant No. 31901173 and 32160285), and Beijing Academy of Agriculture and Forestry Sciences Special Project for Innovation (Grant No. KJCX20230111).
Acknowledgments
We thank the staff at the Xiaotangshan station of Institute of Grassland, Flowers and Ecology for providing logistic support in the field. This study was supported by Natural Science Foundation of Beijing Municipality (Grant No. 5232006), National Natural Science Foundation of China (Grant No. 31901173 and 32160285), and Beijing Academy of Agriculture and Forestry Sciences Special Project for Innovation (Grant No. KJCX20230111).
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.2023.1301461/full#supplementary-material
References
1
AdlerP. B.FajardoA.KleinhesselinkA. R.KraftN. J. (2013). Trait-based tests of coexistence mechanisms. Ecol. Lett.16, 1294–1306. doi: 10.1111/ele.12157
2
AvolioM. L.KoernerS. E.La PierreK. J.WilcoxK. R.WilsonG. W.SmithM. D.et al. (2014). Changes in plant community composition, not diversity, during a decade of nitrogen and phosphorus additions drive above-ground productivity in a tallgrass prairie. J. Ecol.102, 1649–1660. doi: 10.1111/1365-2745.12312
3
BongersF. J.SchmidB.BruelheideH.BongersF.LiS.von OheimbG.et al. (2021). Functional diversity effects on productivity increase with age in a forest biodiversity experiment. Nat. Ecol. Evol.5, 1594–1603. doi: 10.1038/s41559-021-01564-3
4
BordeleauL.PrévostD. (1994). Nodulation and nitrogen fixation in extreme environments. Plant Soil161, 115–125. doi: 10.1007/BF02183092
5
CravenD.EisenhauerN.PearseW. D.HautierY.IsbellF.RoscherC.et al. (2018). Multiple facets of biodiversity drive the diversity–stability relationship. Nat. Ecol. Evol.2, 1579–1587. doi: 10.1038/s41559-018-0647-7
6
de BelloF.LavorelS.DíazS.HarringtonR.CornelissenJ. H.BardgettR. D.et al. (2010). Towards an assessment of multiple ecosystem processes and services via functional traits. Biodivers. Conserv.19, 2873–2893. doi: 10.1007/s10531-010-9850-9
7
de BelloF. D.LavorelS.HallettL. M.ValenciaE.GarnierE.RoscherC.et al. (2021). Functional trait effects on ecosystem stability: assembling the jigsaw puzzle. Trends Ecol. Evol.36, 822–836. doi: 10.1016/j.tree.2021.05.001
8
DiazS.CabidoM. (2001). Vive la difference: plant functional diversity matters to ecosystem processes. Trends Ecol. Evol.16, 646–655. doi: 10.1016/S0169-5347(01)02283-2
9
DíazS.KattgeJ.CornelissenJ. H.WrightI. J.LavorelS.DrayS.et al. (2016). The global spectrum of plant form and function. Nature529, 167–171. doi: 10.1038/nature16489
10
ElserJ. J.BrackenM. E.ClelandE. E.GrunerD. S.HarpoleW. S.HillebrandH.et al. (2007). Global analysis of nitrogen and phosphorus limitation of primary producers in freshwater, marine and terrestrial ecosystems. Ecol. Lett.10, 1135–1142. doi: 10.1111/j.1461-0248.2007.01113.x
11
FayP. A.ProberS. M.HarpoleW. S.KnopsJ. M.BakkerJ. D.BorerE. T.et al. (2015). Grassland productivity limited by multiple nutrients. Nat. Plants1, 15080. doi: 10.1038/nplants.2015.80
12
FoxJ.WeisbergS. (2019). An R companion to applied regression, Third Edition, Sage Publications, Thousand Oaks, CA.
13
GallowayJ. N.TownsendA. R.ErismanJ. W.BekundaM.CaiZ.FreneyJ. R.et al. (2008). Transformation of the nitrogen cycle: recent trends, questions, and potential solutions. Science320I, 889–892. doi: 10.1126/science.1136674
14
GarnierE.NavasM. L.GrigulisK. (2016). Plant Functional Diversity. Organism Traits, Community Structure, and Ecosystem Properties (Oxford: Oxford University Press).
15
GrigulisK.LavorelS.KrainerU.LegayN.BaxendaleC.DumontM.et al. (2013). Relative contributions of plant traits and soil microbial properties to mountain grassland ecosystem services. J. Ecol.101, 47–57. doi: 10.1111/1365-2745.12014
16
GrimeJ. P. (1998). Benefits of plant diversity to ecosystems: immediate, filter and founder effects. J. Ecol.86, 902–910. doi: 10.1046/j.1365-2745.1998.00306.x
17
GrmanE.LauJ. A.SchoolmasterD. R.Jr.GrossK. L. (2010). Mechanisms contributing to stability in ecosystem function depend on the environmental context. Ecol. Lett.13, 1400–1410. doi: 10.1111/j.1461-0248.2010.01533.x
18
HarpoleW. S.NgaiJ. T.ClelandE. E.SeabloomE. W.BorerE. T.BrackenM. E.et al. (2011). Nutrient co-limitation of primary producer communities. Ecol. Lett.14, 852–862. doi: 10.1111/j.1461-0248.2011.01651.x
19
HautierY.SeabloomE. W.BorerE. T.AdlerP. B.HarpoleW. S.HillebrandH.et al. (2014). Eutrophication weakens stabilizing effects of diversity in natural grasslands. Nature508, 521–525. doi: 10.1038/nature13014
20
HautierY.TilmanD.IsbellF.SeabloomE. W.BorerE. T.ReichP. B. (2015). Anthropogenic environmental changes affect ecosystem stability via biodiversity. Science348, 336–340. doi: 10.1126/science.aaa1788
21
HautierY.ZhangP.LoreauM.WilcoxK. R.SeabloomE. W.BorerE. T.et al. (2020). General destabilizing effects of eutrophication on grassland productivity at multiple spatial scales. Nat. Commun.11, 5375. doi: 10.1038/s41467-021-20997-9
22
HeF.TongZ.HannawayD. B.LiX. (2021). Erratic precipitation and clipping frequency reshape the community structure and species stability of leymus chinensis steppe. Ecol. Indic.133, 108432. doi: 10.1016/j.ecolind.2021.108432
23
HectorA.HautierY.SanerP.WackerL.BagchiR.JoshiJ.et al. (2010). General stabilizing effects of plant diversity on grassland productivity through population asynchrony and overyielding. Ecology91, 2213–2220. doi: 10.1890/09-1162
24
HooperD. U.ChapinF. S.EwelJ. J.HectorA.InchaustiP.LavorelS.et al. (2005). Effects of biodiversity on ecosystem functioning: a consensus of current knowledge. Ecol. Monogr.75, 3–35. doi: 10.1016/j.drugalcdep.2010.01.006
25
IsbellF.CravenD.ConnollyJ.LoreauM.SchmidB.BeierkuhnleinC.et al. (2015). Biodiversity increases the resistance of ecosystem productivity to climate extremes. Nature526, 574–577. doi: 10.1038/nature15374
26
IvesA. R.GrossK.KlugJ. L. (1999). Stability and variability in competitive communities. Science286, 542–544. doi: 10.1126/science.286.5439.542
27
JiaX.TaoD.KeY.LiW.YangT.YangY.et al. (2022). Dominant species control effects of nitrogen addition on ecosystem stability. Sci. Total Environ.838, 156060. doi: 10.1016/j.scitotenv.2022.156060
28
LaliberteE.LegendreP. (2010). A distance-based framework for measuring functional diversity from multiple traits. Ecology91, 299–305. doi: 10.1890/08-2244.1
29
LavorelS.GrigulisK. (2012). How fundamental plant functional trait relationships scale-up to trade-offs and synergies in ecosystem services. J. Ecol.100, 128–140. doi: 10.1111/j.1365-2745.2011.01914.x
30
LeBauerD.TresederK. (2008). Nitrogen limitation of net primary productivity in terrestrial ecosystems is globally distributed. Ecology89, 371–379. doi: 10.1890/06-2057.1
31
LefcheckJ. S. (2016). piecewiseSEM: Piecewise structural equation modelling in R for ecology, evolution, and systematics. Methods Ecol. Evol.7, 573–579.12512. doi: 10.1111/2041-210X.12512
32
LepšJ.MájekováM.VítováA.DoležalJ.de BelloF. (2018). Stabilizing effects in temporal fluctuations: management, traits, and species richness in high-diversity communities. Ecology99, 360–371. doi: 10.1002/ecy.2065
33
LepšJ.Osbornová-KosinováJ.RejmánekM. (1982). Community stability, complexity and species life history strategies. Vegetatio50, 53–63. doi: 10.1007/BF00120678
34
LiY.ReichP. B.SchmidB.ShresthaN.FengX.LyuT.et al. (2020). Leaf size of woody dicots predicts ecosystem primary productivity. Ecol. Lett.23, 1003–1013. doi: 10.1111/ele.13503
35
LiuJ.LiX.MaQ.ZhangX.ChenY.IsbellF.et al. (2019). Nitrogen addition reduced ecosystem stability regardless of its impacts on plant diversity. J. Ecol.107, 2427–2435. doi: 10.1111/1365-2745.13187
36
LiuX.ZhangY.HanW.TangA.ShenJ.CuiZ.et al. (2013). Enhanced nitrogen deposition over China. Nature494, 459–462. doi: 10.1038/nature11917
37
LoreauM.de MazancourtC. (2008). Species synchrony and its drivers: Neutral and nonneutral community dynamics in fluctuating environments. Am. Nat.172, 48–66. doi: 10.1086/589746
38
LoreauM.de MazancourtC. (2013). Biodiversity and ecosystem stability: A synthesis of underlying mechanisms. Ecol. Lett.16, 106–115. doi: 10.1111/ele.12073
39
LoreauM.HectorA. (2001). Partitioning selection and complementarity in biodiversity experiments. Nature412, 72–76. doi: 10.1038/35083573
40
MaZ.LiuH.MiZ.ZhangZ.WangY.XuW.et al. (2017). Climate warming reduces the temporal stability of plant community biomass production. Nat. Commun.8, 15378. doi: 10.1038/ncomms15378
41
MájekováM.de BelloF.DolezalJ.LepšJ. (2014). Plant functional traits as determinants of population stability. Ecology95, 2369–2374. doi: 10.1890/13-1880.1
42
PalmrothS.Holm BachL.NordinA.PalmqvistK. (2014). Nitrogen-addition effects on leaf traits and photosynthetic carbon gain of boreal forest understory shrubs. Oecologia175, 457–470. doi: 10.1007/s00442-014-2923-9
43
PenningsS. C.ClarkC. M.ClelandE. E.CollinsS. L.GoughL.GrossK. L.et al. (2005). Do individual plant species show predictable responses to nitrogen addition across multiple experiments? Oikos110, 547–555. doi: 10.1111/j.0030-1299.2005.13792.x
44
PenuelasJ.PoulterB.SardansJ.CiaisP.Van Der VeldeM.BoppL.et al. (2013). Human-induced nitrogen-phosphorus imbalances alter natural and managed ecosystems across the globe. Nat. Commun.4, 1–10. doi: 10.1038/ncomms3934
45
Pérez-HarguindeguyN.DíazS.GarnierE.LavorelS.PoorterH.JaureguiberryP.et al. (2013). New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot.64, 715. doi: 10.1071/BT12225
46
PinheiroJ.BatesD.DebRoyS.SarkarD.R Core Team. (2007). Linear and nonlinear mixed effects models. R Package Version.3, 57.
47
PolleyH. W.YangC.WilseyB. J.FayP. A. (2020). Temporal stability of grassland metacommunities is regulated more by community functional traits than species diversity. Ecosphere11, e03178. doi: 10.1002/ecs2.3178
48
RadujkovićD.VerbruggenE.SeabloomE. W.BahnM.BiedermanL. A.BorerE. 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
49
R Development Core Team. (2018). R: A language and environment for statistical computing (Vienna, Austria: R Foundation for Statistical Computing).
50
ReichP. B.CornelissenH. (2014). The world-wide 'fast-slow' plant economics spectrum: A traits manifesto. J. Ecol.102, 275–301. doi: 10.1111/1365-2745.12211
51
RenF.SongW.ChenL.MiZ.ZhangZ.ZhuW.et al. (2016). Phosphorus does not alleviate the negative effect of nitrogen enrichment on legume performance in an alpine grassland. J. Plant Ecol.6, 1–9. doi: 10.1093/jpe/rtw089
52
RochaM. R.GaedkeU.VasseurD. A. (2011). Functionally similar species have similar dynamics. J. Ecol.99, 1453–1459. doi: 10.1111/j.1365-2745.2011.01893.x
53
SchnabelF.LiuX.KunzM.BarryK. E.BongersF. J.BruelheideH.et al. (2021). Species richness stabilizes productivity via asynchrony and drought-tolerance diversity in a large-scale tree biodiversity experiment. Sci. Adv.7, eabk1643. doi: 10.1126/sciadv.abk1643
54
SongM.YuF. (2015). Reduced compensatory effects explain the nitrogen-mediated reduction in stability of an alpine meadow on the tibetan plateau. New Phytol.207, 70–77. doi: 10.1111/nph.13329
55
SudingK. N.CollinsS. L.GoughL.ClarkC.ClelandE. E.GrossK. L.et al. (2005). Functional- and abundance-based mechanisms explain diversity loss due to N fertilization. Proc. Natl. Acad. Sci.102, 4387–4392. doi: 10.1073/pnas.0408648102
56
SunS.JinD.LiR. (2006). Leaf emergence in relation to leaf traits of woody species in eastern Chinese Quercus fabri forests. Acta Oecol.30, 212–222. doi: 10.1016/j.actao.2006.04.001
57
ThibautL. M.ConnollyS. R. (2013). Understanding diversity-stability relationships: Towards a unified model of portfolio effects. Ecol. Lett.16, 140–150. doi: 10.1111/ele.12019
58
TilmanD. (1996). Biodiversity population versus ecosystem stability. Ecology77, 358–363. doi: 10.2307/2265614
59
van der PlasF.Schröder-GeorgiT.WeigeltA.BarryK.MeyerS.AlzateA.et al. (2020). Plant traits alone are poor predictors of ecosystem properties and long-term ecosystem functioning. Nat. Ecol. Evol.4, 1–14. doi: 10.1038/s41559-020-01316-9
60
VergnonR.DulvyN. K.FreckletonR. P. (2009). Niches versus neutrality: uncovering the drivers of diversity in a species-rich community. Ecol. Lett.12, 1079–1090. doi: 10.1111/j.1461-0248.2009.01364.x
61
WangT.GuoC.SangS.LiuY.LiuG.QiD.et al. (2021). Temporal stability and maintenance mechanisms of alpine meadow communities under clipping and fertilization. Ecol. Evol.11, 15545–15555. doi: 10.1002/ece3.8128
62
WangS.LamyT.HallettL. M.LoreauM. (2019). Stability and synchrony across ecological hierarchies in heterogeneous metacommunities: Linking theory to data. Ecography42, 1–12. doi: 10.1111/ecog.04290
63
WangC.LiX.HuY.ZhengR.HouY. (2023). Nitrogen addition weakens the biodiversity–multifunctionality relationships across soil profiles in a grassland assemblage. Agr. Ecosyst. Environ.342, 108241. doi: 10.1016/j.agee.2022.108241
64
WangY.NiuX.ZhaoL.LiangC.MiaoB.ZhangQ.et al. (2020). Biotic stability mechanisms in Inner Mongolian grassland. P Roy Soc. B Biol. Sci.287, 20200675. doi: 10.1098/rspb.2020.0675
65
WangC.RenF.ZhouX.MaW.LiangC.WangJ.et al. (2020). Variations in the nitrogen saturation threshold of soil respiration in grassland ecosystems. Biogeochemistry148, 1–14. doi: 10.1007/s10533-020-00661-y
66
WrightI. J.ReichP. B.WestobyM.AckerlyD. D.BaruchZ.BongersF.et al. (2004). The worldwide leaf economics spectrum. Nature428, 821–827. doi: 10.1038/nature02403
67
XuZ.RenH.LiM. H.van RuijvenJ.HanX.WanS.et al. (2015). Environmental changes drive the temporal stability of semi-arid natural grasslands through altering species asynchrony. J. Ecol.103, 1308–1316. doi: 10.1111/1365-2745.12441
68
YanY.ConnollyJ.LiangM.JiangL.WangS. (2021). Mechanistic links between biodiversity effects on ecosystem functioning and stability in a multi-site grassland experiment. J. Ecol.109, 3370–3378. doi: 10.1111/1365-2745.13725
69
YangH.JiangL.LiL.LiA.WuM.WanS. (2012). Diversity-dependent stability under mowing and nutrient addition: evidence from a 7-year grassland experiment. Ecol. Lett.15, 619–626. doi: 10.1111/j.1461-0248.2012.01778.x
70
ZhangY.FengJ.LoreauM.HeN.HanX.JiangL. (2019). Nitrogen addition does not reduce the role of spatial asynchrony in stabilising grassland communities. Ecol. Lett.22, 563–571. doi: 10.1111/ele.13212
71
ZhangY.LoreauM.LuX.HeN.ZhangG.HanX. (2016). Nitrogen enrichment weakens ecosystem stability through decreased species asynchrony and population stability in a temperate grassland. Glob. Change Biol.22, 1445–1455. doi: 10.1111/gcb.13140
72
ZhengS.ChiY.YangX.LiW.LanZ.BaiY. (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
73
ZhouX.GuoZ.ZhangP.DuG. (2018). Shift in community functional composition following nitrogen fertilization in an alpine meadow through intraspecific trait variation and community composition change. Plant Soil431, 289–302. doi: 10.1007/s11104-018-3771-x
Summary
Keywords
community stability, fast–slow traits, functional diversity, nitrogen addition, species stability, species asynchrony
Citation
Suonan J, Lu X, Li X, Hautier Y and Wang C (2023) Nitrogen addition strengthens the stabilizing effect of biodiversity on productivity by increasing plant trait diversity and species asynchrony in the artificial grassland communities. Front. Plant Sci. 14:1301461. doi: 10.3389/fpls.2023.1301461
Received
25 September 2023
Accepted
03 November 2023
Published
20 November 2023
Volume
14 - 2023
Edited by
Yu-Long Feng, Shenyang Agricultural University, China
Reviewed by
Huiying Liu, East China Normal University, China; Zhenhua Zhang, Chinese Academy of Sciences (CAS), China
Updates
Copyright
© 2023 Suonan, Lu, Li, Hautier and Wang.
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: Yann Hautier, yannhautier@gmail.com; Chao Wang, wangchao@grass-env.com
†These authors have contributed equally to this work
Disclaimer
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.