Nutrients and herbivores impact grassland stability across spatial scales through different pathways

Local nutrient addition drives plant diversity losses but not biotic homogenization in global grasslands

Nutrient enrichment typically causes local plant diversity declines. A common but untested expectation is that nutrient enrichment also reduces variation in nutrient conditions among localities and selects for a smaller pool of species, causing greater diversity declines at larger than local scales and thus biotic homogenization. Here we apply a framework that links changes in species richness across scales to changes in the numbers of spatially restricted and widespread species for a standardized nutrient addition experiment across 72 grasslands on six continents. Overall, we find proportionally similar species loss at local and larger scales, suggesting similar declines of spatially restricted and widespread species, and no biotic homogenization after 4 years and up to 14 years of treatment. These patterns of diversity changes are generally consistent across species groups. Thus, nutrient enrichment poses threats to plant diversity, including for widespread species that are often critical for ecosystem functions.

The spatiotemporal stability of plant diversity is disconnected from biomass stability in response to human activities in a South American temperate grassland

Human activities alter biomass, nutrient availability, and species dominance in grasslands, impacting their richness, composition, and biomass production. Stability (invariability in time or space) can inform the predictability of plant communities in response to human activities. However, this measure has been simplistically analyzed for temporal (interannual) changes in live biomass, disregarding their spatial stability and the temporal stability of other plant community attributes. Moreover, the simultaneous analysis of temporal and spatial stabilities of plant communities has been scarcely assessed. Here, we test how biomass removal and nutrient addition simultaneously modify the temporal and spatial stabilities of plant richness (α diversity), composition dissimilarity (β diversity), aboveground live biomass, and the role of plant species dominance in the stability responses. We conducted a factorial experiment of biomass removal (grazing, mowing, or intact -no removal-) and nutrient addition (unfertilized or fertilized with nitrogen, phosphorus, and potassium) in a temperate grassland of Argentina, South America. We replicated the experiment in 6 blocks over 10 years to estimate the temporal and spatial stabilities of the plant community. The spatiotemporal stability of plant richness and composition dissimilarity decreased in the intact grassland, while the temporal stability of live biomass increased, compared to the grazed and mowed grasslands. Nutrient addition reduced the spatiotemporal stability of live biomass and the spatial stability of plant richness. The stabilities of species richness as well as that of composition dissimilarity were negatively associated with plant dominance, while the live biomass stability was not. Our results suggest that simplifying the effect of biomass removal and nutrient addition on grassland stability is not feasible, as plant diversity stability responses are not surrogates for biomass stability. The contrasting spatiotemporal stability responses of plant diversity and biomass represent a step forward in predicting human activities' impact over time and across space in temperate grasslands.

Atmospheric nitrogen deposition is related to plant biodiversity loss at multiple spatial scales

Due to various human activities, including intensive agriculture, traffic, and the burning of fossil fuels, in many parts of the world, current levels of reactive nitrogen emissions strongly exceed pre-industrial levels. Previous studies have shown that the atmospheric deposition of these excess nitrogen compounds onto semi-natural terrestrial environments has negative consequences for plant diversity. However, these previous studies mostly investigated biodiversity loss at local spatial scales, that is, at the scales of plots of typically a few square meters. Whether increased atmospheric nitrogen deposition also affects plant diversity at larger spatial scales remains unknown. Here, using grassland plant community data collected in 765 plots, across 153 different sites and 9 countries in northwestern Europe, we investigate whether relationships between atmospheric nitrogen deposition and plant biodiversity are scale-dependent. We found that high levels of atmospheric nitrogen deposition were associated with low levels of plant species richness at the plot scale but also at the scale of sites and regions. The presence of 39% of plant species was negatively associated with increasing levels of nitrogen deposition at large (site) scales, while only 1.5% of the species became more common with increasing nitrogen deposition, indicating that large-scale biodiversity changes were mostly driven by "loser" species, while "winner" species profiting from high N deposition were rare. Some of the "loser" species whose site presence was negatively associated with atmospheric nitrogen deposition are listed as "threatened" in at least some EU member states, suggesting that nitrogen deposition may be a key contributor to their threat status. Hence, reductions in reactive nitrogen emissions will likely benefit plant diversity not only at local but also at larger spatial scales.

Herbivore exclusion stabilizes alpine grassland biomass production across spatial scales

There is growing evidence that land-use management practices such as livestock grazing can strongly impact the local diversity, functioning, and stability of grassland communities. However, whether these impacts depend on environmental condition and propagate to larger spatial scales remains unclear. Using an 8-year grassland exclosure experiment conducted at nine sites in the Tibetan Plateau covering a large precipitation gradient, we found that herbivore exclusion increased the temporal stability of alpine grassland biomass production at both the local and larger (site) spatial scales. Higher local community stability was attributed to greater stability of dominant species, whereas higher stability at the larger scale was linked to higher spatial asynchrony of productivity among local communities. Additionally, sites with higher mean annual precipitation had lower dominant species stability and lower grassland stability at both the spatial scales considered. Our study provides novel evidence that livestock grazing can impair grassland stability across spatial scales and climatic gradients.

Disease decreases variation in host community structure in an old-field grassland

Disease may drive variation in host community structure by modifying the interplay of deterministic and stochastic processes that shape communities. For instance, deterministic processes like ecological selection can benefit species less impacted by disease. When communities have higher levels of disease and disease consistently selects for certain host species, this can reduce variation in host community composition. On the other hand, when host communities are less impacted by disease and selection is weaker, stochastic processes (e.g., drift, dispersal) may play a bigger role in host community structure, which can increase variation among communities. While effects of disease on host community structure have been quantified in field experiments, few have addressed the role of disease in modulating variation in structure among host communities. To address this, we conducted a field experiment spanning three years, using a tractable system: foliar fungal pathogens in an old-field grassland community dominated by the grass Lolium arundinaceum, tall fescue. We reduced foliar fungal disease burden in replicate host communities (experimental plots in intact vegetation) in three fungicide regimens that varied in the seasonal duration of fungicide treatment and included a fungicide-free control. We measured host diversity, biomass, and variation in community structure among replicate communities. Disease reduction generally decreased plant richness and increased aboveground biomass relative to communities experiencing ambient levels of disease. These changes in richness and aboveground biomass were consistent across years despite changes in structure of the plant communities over the experiment's three years. Importantly, disease reduction amplified host community variation, suggesting that disease diminished the degree to which host communities were structured by stochastic processes. These results of experimental disease reduction both highlight the potential importance of stochastic processes in plant communities and reveal the potential for disease to regulate variation in host community structure.

Multidimensional responses of grassland stability to eutrophication

Eutrophication usually impacts grassland biodiversity, community composition, and biomass production, but its impact on the stability of these community aspects is unclear. One challenge is that stability has many facets that can be tightly correlated (low dimensionality) or highly disparate (high dimensionality). Using standardized experiments in 55 grassland sites from a globally distributed experiment (NutNet), we quantify the effects of nutrient addition on five facets of stability (temporal invariability, resistance during dry and wet growing seasons, recovery after dry and wet growing seasons), measured on three community aspects (aboveground biomass, community composition, and species richness). Nutrient addition reduces the temporal invariability and resistance of species richness and community composition during dry and wet growing seasons, but does not affect those of biomass. Different stability measures are largely uncorrelated under both ambient and eutrophic conditions, indicating consistently high dimensionality. Harnessing the dimensionality of ecological stability provides insights for predicting grassland responses to global environmental change.

Diversity-stability relationships across organism groups and ecosystem types become decoupled across spatial scales

The relationship between biodiversity and stability, or its inverse, temporal variability, is multidimensional and complex. Temporal variability in aggregate properties, like total biomass or abundance, is typically lower in communities with higher species diversity (i.e., the diversity-stability relationship [DSR]). At broader spatial extents, regional-scale aggregate variability is also lower with higher regional diversity (in plant systems) and with lower spatial synchrony. However, focusing exclusively on aggregate properties of communities may overlook potentially destabilizing compositional shifts. It is not yet clear how diversity is related to different components of variability across spatial scales, nor whether regional DSRs emerge across a broad range of organisms and ecosystem types. To test these questions, we compiled a large collection of long-term metacommunity data spanning a wide range of taxonomic groups (e.g., birds, fish, plants, invertebrates) and ecosystem types (e.g., deserts, forests, oceans). We applied a newly developed quantitative framework for jointly analyzing aggregate and compositional variability across scales. We quantified DSRs for composition and aggregate variability in local communities and metacommunities. At the local scale, more diverse communities were less variable, but this effect was stronger for aggregate than compositional properties. We found no stabilizing effect of γ-diversity on metacommunity variability, but β-diversity played a strong role in reducing compositional spatial synchrony, which reduced regional variability. Spatial synchrony differed among taxa, suggesting differences in stabilization by spatial processes. However, metacommunity variability was more strongly driven by local variability than by spatial synchrony. Across a broader range of taxa, our results suggest that high γ-diversity does not consistently stabilize aggregate properties at regional scales without sufficient spatial β-diversity to reduce spatial synchrony.

Negative effects of phosphorus addition outweigh effects of arbuscular mycorrhizal fungi and nitrogen addition on grassland temporal stability in the eastern Eurasian desert steppe

The temporal stability of grassland plant communities is substantially affected by soil nutrient enrichment. However, the potential main and interactive effects of arbuscular mycorrhizal fungi (AMF) and soil nitrogen (N) and phosphorus (P) enrichment on the stability of plant productivity have not yet been clarified. We combined a three-year in situ field experiment to assess the impacts of soil fertilization and AMF on the stability of plant productivity. P addition decreased the stability of plant productivity by increasing the standard deviation relative to the mean of plant productivity. However, compared to species richness, the stability of C3 grasses and other functional groups asynchrony were the most important drivers changing the stability of plant productivity. The negative impacts of P addition overrode the impacts of AMF on the stability of plant productivity. Overall, our study suggests the importance of soil nutrient availability over AMF in terms of shaping the stability of plant productivity. Our results also suggest that three-year anthropogenic soil nutrient enrichment could reduce the stability of plant communities in grassland regardless of AMF in the P-limited grassland ecosystem.

Universal beta-diversity-functioning relationships are neither observed nor expected

Widespread evidence shows that local species richness (α-diversity) loss hampers the biomass production and stability of ecosystems. β-Diversity, namely the variation of species compositions among different ecological communities, represents another important biodiversity component, but studies on how it drives ecosystem functioning show mixed results. We argue that to better understand the importance of β-diversity we need to consider it across contexts. We focus on three scenarios that cause gradients in β-diversity: changes in (i) abiotic heterogeneity, (ii) habitat isolation, and (iii) species pool richness. We show that across these scenarios we should not expect universally positive relationships between β-diversity, production, and ecosystem stability. Nevertheless, predictable relationships between β-diversity and ecosystem functioning do exist in specific contexts, and can reconcile seemingly contrasting empirical relationships.

Increasing connections among temporal invariability, resistance and resilience of alpine grasslands on the Tibetan Plateau

Ecological stability contains multiple components, such as temporal invariability, resistance and resilience. Understanding the response of stability components to perturbations is beneficial for optimizing the management of biodiversity and ecosystem functioning. Although previous studies have investigated the effects of multiple perturbations on each stability component, few studies simultaneously measure the multiple stability components and their relationships. Alpine grasslands on the Tibetan Plateau are exposed to co-occurring perturbations, including climate change and human activities. Here, we quantified three stability components (temporal invariability, resistance, and resilience) of alpine grasslands on the Tibetan Plateau during periods of high (2000-2008) and low (2009-2017) human activity intensity, respectively. We focused on the effects of climate variables (temperature, precipitation, radiation) and human activities (grazing intensity) on covariation among stability components. The results show that (1) for periods of high and low human activity, temporal invariability was positively correlated with resistance and resilience, while resistance was independent of resilience; (2) the dimensionality of alpine grasslands decreased by almost 10%, from 0.61 in the first period to 0.55 in the second period, suggesting the increasing connections among temporal invariability, resistance and resilience of alpine grasslands; and (3) temperature but not grazing intensity dominated the changes in the dimensionality of stability. These findings improve our understanding of multi-dimensional stability and highlight the importance of climate variability on alpine grassland stability on the Tibetan Plateau.

Loss of grazing by large mammalian herbivores can destabilize the soil carbon pool

Grazing by mammalian herbivores can be a climate mitigation strategy as it influences the size and stability of a large soil carbon (soil-C) pool (more than 500 Pg C in the world's grasslands, steppes, and savannas). With continuing declines in the numbers of large mammalian herbivores, the resultant loss in grazer functions can be consequential for this soil-C pool and ultimately for the global carbon cycle. While herbivore effects on the size of the soil-C pool and the conditions under which they lead to gain or loss in soil-C are becoming increasingly clear, their effect on the equally important aspect of stability of soil-C remains unknown. We used a replicated long-term field experiment in the Trans-Himalayan grazing ecosystem to evaluate the consequences of herbivore exclusion on interannual fluctuations in soil-C (2006 to 2021). Interannual fluctuations in soil-C and soil-N were 30 to 40% higher after herbivore exclusion than under grazing. Structural equation modeling suggested that grazing appears to mediate the stabilizing versus destabilizing influences of nitrogen (N) on soil-C. This may explain why N addition stimulates soil-C loss in the absence of herbivores around the world. Herbivore loss, and the consequent decline in grazer functions, can therefore undermine the stability of soil-C. Soil-C is not inert but a very dynamic pool. It can provide nature-based climate solutions by conserving and restoring a functional role of large mammalian herbivores that extends to the stoichiometric coupling between soil-C and soil-N.