Topsoil removal, among other restoration measures, has been recognized as one of the most successful methods to restore biodiversity and ecosystem functioning in European grasslands. However, knowledge about how removal as well as other restoration methods influence interactions between plant and microbial communities is very limited. The aims of the current study were to understand the impact of topsoil removal on plant-microorganism interactions and on soil nitrogen (N) mineralization, as one example of ecosystem functioning. We examined how three different grassland restoration methods, namely 'Harvest only', 'Topsoil removal' and 'Topsoil removal + Propagules (plant seed addition)', affected i) the interactions between plants and soil microorganisms, ii) soil microbial community assembly processes, and iii) soil N mineralization. We compared the outcome of these three restoration methods to initial degraded and target semi-natural grasslands in the Canton of Zurich, Switzerland. We were able to show that 'Topsoil removal' and 'Topsoil removal + Propagules', but not 'Harvest only', reduced the soil total N pool and available N concentration, but increased soil N mineralization and strengthened the plant-microorganism interactions. Microbial community assembly processes shifted towards more deterministic after both topsoil removal treatments. These shifts could be attributed to an increase in dispersal limitation and selection due to stronger interactions between plants and soil microorganisms. The negative relationship between soil N mineralization and microbial community stochasticity indicated that microbial assembly processes, to some extent, can be incorporated into model predictions of soil functions. Overall, the results suggest that topsoil removal may change the microbial assembly processes and thus the functioning of grassland ecosystems by enhancing the interaction between plants and soil microorganisms.

Aims: We introduce ReSurveyEurope - a new data source of resurveyed vegetation plots in Europe, compiled by a collaborative network of vegetation scientists. We describe the scope of this initiative, provide an overview of currently available data, governance, data contribution rules, and accessibility. In addition, we outline further steps, including potential research questions. Results: ReSurveyEurope includes resurveyed vegetation plots from all habitats. Version 1.0 of ReSurveyEurope contains 283,135 observations (i.e., individual surveys of each plot) from 79,190 plots sampled in 449 independent resurvey projects. Of these, 62,139 (78%) are permanent plots, that is, marked in situ, or located with GPS, which allow for high spatial accuracy in resurvey. The remaining 17,051 (22%) plots are from studies in which plots from the initial survey could not be exactly relocated. Four data sets, which together account for 28,470 (36%) plots, provide only presence/absence information on plant species, while the remaining 50,720 (64%) plots contain abundance information (e.g., percentage cover or cover-abundance classes such as variants of the Braun-Blanquet scale). The oldest plots were sampled in 1911 in the Swiss Alps, while most plots were sampled between 1950 and 2020. Conclusions: ReSurveyEurope is a new resource to address a wide range of research questions on fine-scale changes in European vegetation. The initiative is devoted to an inclusive and transparent governance and data usage approach, based on slightly adapted rules of the well-established European Vegetation Archive (EVA). ReSurvey:Europe data are ready for use, and proposals for analyses of the data set can be submitted at any time to the coordinators. Still, further data contributions are highly welcome.

Ecological theory posits that temporal stability patterns in plant populations are associated with differences in species' ecological strategies. However, empirical evidence is lacking about which traits, or trade-offs, underlie species stability, especially across different biomes. We compiled a worldwide collection of long-term permanent vegetation records (greater than 7000 plots from 78 datasets) from a large range of habitats which we combined with existing trait databases. We tested whether the observed inter-annual variability in species abundance (coefficient of variation) was related to multiple individual traits. We found that populations with greater leaf dry matter content and seed mass were more stable over time. Despite the variability explained by these traits being low, their effect was consistent across different datasets. Other traits played a significant, albeit weaker, role in species stability, and the inclusion of multi-variate axes or phylogeny did not substantially modify nor improve predictions. These results provide empirical evidence and highlight the relevance of specific ecological trade-offs, i.e. in different resource-use and dispersal strategies, for plant populations stability across multiple biomes. Further research is, however, necessary to integrate and evaluate the role of other specific traits, often not available in databases, and intraspecific trait variability in modulating species stability.

Aim: The microbial metabolic quotient (MMQ; mg CO2- C/mg MBC/h), defined as the amount of microbial CO2 respired (MR; mg CO2- C/kg soil/h) per unit of microbial biomass C (MBC; mg C/kg soil), is a key parameter for understanding the microbial regulation of the carbon (C) cycle, including soil C sequestration. Here, we experimentally tested hypotheses about the individual and interactive effects of multiple nutrient addition (nitrogen + phosphorus + potassium + micronutrients) and herbivore exclusion on MR, MBC and MMQ across 23 sites (five continents). Our sites encompassed a wide range of edaphoclimatic conditions; thus, we assessed which edaphoclimatic variables affected MMQ the most and how they interacted with our treatments.Location: Australia, Asia, Europe, North/South America.Time period: 2015-2016.Major taxa: Soil microbes.Methods: Soils were collected from plots with established experimental treatments. MR was assessed in a 5 week laboratory incubation without glucose addition, MBC via substrate-induced respiration. MMQ was calculated as MR/MBC and corrected for soil temperatures (MMQsoil). Using linear mixed effects models (LMMs) and structural equation models (SEMs), we analysed how edaphoclimatic characteristics and treatments interactively affected MMQsoil.Results: MMQsoil was higher in locations with higher mean annual temperature, lower water holding capacity and lower soil organic C concentration, but did not respond to our treatments across sites as neither MR nor MBC changed. We attributed this relative homeostasis to our treatments to the modulating influence of edaphoclimatic variables. For example, herbivore exclusion, regardless of fertilization, led to greater MMQsoil only at sites with lower soil organic C (< 1.7%).Main conclusions: Our results pinpoint the main variables related to MMQsoil across grasslands and emphasize the importance of the local edaphoclimatic conditions in controlling the response of the C cycle to anthropogenic stressors. By testing hypotheses about MMQsoil across global edaphoclimatic gradients, this work also helps to align the conflicting results of prior studies.

Plant damage by invertebrate herbivores and pathogens influences the dynamics of grassland ecosystems, but anthropogenic changes in nitrogen and phosphorus availability can modify these relationships. Using a globally distributed experiment, we describe leaf damage on 153 plant taxa from 27 grasslands worldwide, under ambient conditions and with experimentally elevated nitrogen and phosphorus. Invertebrate damage significantly increased with nitrogen addition, especially in grasses and non‐leguminous forbs. Pathogen damage increased with nitrogen in grasses and legumes but not forbs. Effects of phosphorus were generally weaker. Damage was higher in grasslands with more precipitation, but climatic conditions did not change effects of nutrients on leaf damage. On average, invertebrate damage was relatively higher on legumes and pathogen damage was relatively higher on grasses. Community‐weighted mean damage reflected these functional group patterns, with no effects of N on community‐weighted pathogen damage (due to opposing responses of grasses and forbs) but stronger effects of N on community‐weighted invertebrate damage (due to consistent responses of grasses and forbs). Synthesis . As human‐induced inputs of nitrogen and phosphorus continue to increase, understanding their impacts on invertebrate and pathogen damage becomes increasingly important. Our results demonstrate that eutrophication frequently increases plant damage and that damage increases with precipitation across a wide array of grasslands. Invertebrate and pathogen damage in grasslands is likely to increase in the future, with potential consequences for plant, invertebrate and pathogen communities, as well as the transfer of energy and nutrients across trophic levels.

Grassland ecosystems cover around 37% of the ice-free land surface on Earth and have critical socioeconomic importance globally. As in many terrestrial ecosystems, biological dinitrogen (N-2) fixation represents an essential natural source of nitrogen (N). The ability to fix atmospheric N-2 is limited to diazotrophs, a diverse guild of bacteria and archaea. To elucidate the abiotic (climatic, edaphic), biotic (vegetation), and spatial factors that govern diazotrophic community composition in global grassland soils, amplicon sequencing of the dinitrogenase reductase gene-nifH-was performed on samples from a replicated standardized nutrient [N, phosphorus (P)] addition experiment in 23 grassland sites spanning four continents. Sites harbored distinct and diverse diazotrophic communities, with most of reads assigned to diazotrophic taxa within the Alphaproteobacteria (e.g., Rhizobiales), Cyanobacteria (e.g., Nostocales), and Deltaproteobacteria (e.g., Desulforomonadales) groups. Likely because of the wide range of climatic and edaphic conditions and spatial distance among sampling sites, only a few of the taxa were present at all sites. The best model describing the variation among soil diazotrophic communities at the OTU level combined climate seasonality (temperature in the wettest quarter and precipitation in the warmest quarter) with edaphic (C:N ratio, soil texture) and vegetation factors (various perennial plant covers). Additionally, spatial variables (geographic distance) correlated with diazotrophic community variation, suggesting an interplay of environmental variables and spatial distance. The diazotrophic communities appeared to be resilient to elevated nutrient levels, as 2-4 years of chronic N and P additions had little effect on the community composition. However, it remains to be seen, whether changes in the community composition occur after exposure to long-term, chronic fertilization regimes.

Background and aims The amount of nitrogen (N) derived from symbiotic N 2 fixation by legumes in grasslands might be affected by anthropogenic N and phosphorus (P) inputs, but the underlying mechanisms are not known. Methods We evaluated symbiotic N 2 fixation in 17 natural and semi-natural grasslands on four continents that are subjected to the same full-factorial N and P addition experiment, using the 15 N natural abundance method. Results N as well as combined N and P (NP) addition reduced aboveground legume biomass by 65% and 45%, respectively, compared to the control, whereas P addition had no significant impact. Addition of N and/or P had no significant effect on the symbiotic N 2 fixation per unit legume biomass. In consequence, the amount of N fixed annually per grassland area was less than half in the N addition treatments compared to control and P addition, irrespective of whether the dominant legumes were annuals or perennials. Conclusion Our results reveal that N addition mainly impacts symbiotic N 2 fixation via reduced biomass of legumes rather than changes in N 2 fixation per unit legume biomass. The results show that soil N enrichment by anthropogenic activities significantly reduces N 2 fixation in grasslands, and these effects cannot be reversed by additional P amendment.

Ecological models predict that the effects of mammalian herbivore exclusion on plant diversity depend on resource availability and plant exposure to ungulate grazing over evolutionary time. Using an experiment replicated in 57 grasslands on six continents, with contrasting evolutionary history of grazing, we tested how resources (mean annual precipitation and soil nutrients) determine herbivore exclusion effects on plant diversity, richness and evenness. Here we show that at sites with a long history of ungulate grazing, herbivore exclusion reduced plant diversity by reducing both richness and evenness and the responses of richness and diversity to herbivore exclusion decreased with mean annual precipitation. At sites with a short history of grazing, the effects of herbivore exclusion were not related to precipitation but differed for native and exotic plant richness. Thus, plant species' evolutionary history of grazing continues to shape the response of the world's grasslands to changing mammalian herbivory. A NutNet experiment in 57 grasslands across six continents shows that when herbivores are excluded from grasslands with a long coevolutionary history of grazing plant diversity is reduced, while in grasslands without a long grazing history the evolutionary history of the plant species regulates the response of plant diversity.

Analysing temporal patterns in plant communities is extremely important to quantify the extent and the consequences of ecological changes, especially considering the current biodiversity crisis. Long-term data collected through the regular sampling of permanent plots represent the most accurate resource to study ecological succession, analyse the stability of a community over time and understand the mechanisms driving vegetation change. We hereby present the LOng-Term Vegetation Sampling (LOTVS) initiative, a global collection of vegetation time-series derived from the regular monitoring of vascular plants in permanent plots. With 79 datasets from five continents and 7789 vegetation time-series monitored for at least six years and mostly on an annual basis, LOTVS possibly represents the largest collection of temporally fine-grained vegetation time-series derived from permanent plots and made accessible to the research community. As such, it has an outstanding potential to support innovative research in the fields of vegetation science, plant ecology and temporal ecology.

To evaluate how increased anthropogenic nutrient inputs alter carbon cycling in grasslands, we conducted a litter decomposition study across 20 temperate grasslands on three continents within the Nutrient Network, a globally distributed nutrient enrichment experiment We determined the effects of addition of experimental nitrogen (N), phosphorus (P) and potassium plus micronutrient (K-mu) on decomposition of a common tree leaf litter in a long-term study (maximum of 7 years; exact deployment period varied across sites). The use of higher order decomposition models allowed us to distinguish between the effects of nutrients on early- versus late-stage decomposition. Across continents, the addition of N (but not other nutrients) accelerated early-stage decomposition and slowed late-stage decomposition, increasing the slowly decomposing fraction by 28% and the overall litter mean residence time by 58%. Synthesis. Using a novel, long-term cross-site experiment, we found widespread evidence that N enhances the early stages of above-ground plant litter decomposition across diverse and widespread temperate grassland sites but slows late-stage decomposition. These findings were corroborated by fitting the data to multiple decomposition models and have implications for N effects on soil organic matter formation. For example, following N enrichment, increased microbial processing of litter substrates early in decomposition could promote the production and transfer of low molecular weight compounds to soils and potentially enhance the stabilization of mineral-associated organic matter. By contrast, by slowing late-stage decomposition, N enrichment could promote particulate organic matter (POM) accumulation. Such hypotheses deserve further testing.

Forest ecosystems sequester approximately half of the world’s organic carbon (C), most of it in the soil. The amount of soil C stored depends on the input and decomposition rate of soil organic matter (OM), which is controlled by the abundance and composition of the microbial and invertebrate communities, soil physico-chemical properties, and (micro)-climatic conditions. Although many studies have assessed how these site-specific climatic and soil properties affect the decomposition of fresh OM, differences in the type and quality of the OM substrate used, make it difficult to compare and extrapolate results across larger scales. Here, we used standard wood stakes made from aspen (Populus tremuloides Michx.) and loblolly pine (Pinus taeda L.) to explore how climate and abiotic soil properties affect wood decomposition across 44 unharvested forest stands located across the northern hemisphere. Stakes were placed in three locations: (i) on top of the surface organic horizons (surface), (ii) at the interface between the surface organic horizons and mineral soil (interface), and (iii) into the mineral soil (mineral). Decomposition rates of both wood species was greatest for mineral stakes and lowest for stakes placed on the surface organic horizons, but aspen stakes decomposed faster than pine stakes. Our models explained 44 and 36% of the total variation in decomposition for aspen surface and interface stakes, but only 0.1% (surface), 12% (interface), 7% (mineral) for pine, and 7% for mineral aspen stakes. Generally, air temperature was positively, precipitation negatively related to wood stake decomposition. Climatic variables were stronger predictors of decomposition than soil properties (surface C:nitrogen ratio, mineral C concentration, and pH), regardless of stake location or wood species. However, climate-only models failed in explaining wood decomposition, pointing toward the importance of including local-site properties when predicting wood decomposition. The difficulties we had in explaining the variability in wood decomposition, especially for pine and mineral soil stakes, highlight the need to continue assessing drivers of decomposition across large global scales to better understand and estimate surface and belowground C cycling, and understand the drivers and mechanisms that affect C pools, CO2 emissions, and nutrient cycles.

1. Evaluation of restoration activities is indispensable to assess the extent to which targets have been reached. Usually, the main goal of ecological restoration is to restore biodiversity and ecosystem functioning, but validation is often based on a single indicator, which may or may not cope with whole-ecosystem dynamics. Network analyses are, however, powerful tools, allowing to examine both the recovery of various biotic and abiotic properties and the integrated response at community and ecosystem level. 2. We used restoration sites where topsoil was removed from former intensively managed grassland and seeds were added. These sites were between 3 and 32 years old. We assessed how plants, soil biota, soil properties and correlation-based interactions between biotic communities and their abiotic environment developed over time and compared the results with (i) intensively managed (not restored), and (ii) well-preserved targeted semi-natural grasslands. 3. Plant, nematode, fungal and prokaryotic diversity and community structures of the restored grasslands revealed clear successional patterns and followed similar trajectories towards targeted semi-natural grasslands. All biotic communities reached targeted diversity levels no later than 18 years post-restoration. 4. Ecological networks of intensively managed and short-term (similar to 4 years) restored grasslands were less tightly connected compared to those found in mid- and long-term (similar to 18-30 years) restored and target grasslands. Restoration specifically enhanced interactions among biotic communities, but reduced interactions between biotic communities and their abiotic environment as well as interactions among abiotic properties in the short- and mid-term. 5. Synthesis and applications: Overall, our study demonstrated that topsoil removal and seed addition were successful in restoring diverse, tightly coupled and well-connected biotic communities above- and below-ground similar to those found in the semi-natural grasslands that were restoration targets. Network analyses proved to be powerful in examining the long-term re-establishment of functionally connected biotic communities in restored ecosystems. Thus, we provide an approach to holistically assess restoration activities by notably considering the complexity of ecosystems, much in contrast to most traditional approaches.

Abstract Increasing evidence suggests that remotely sensed spectral diversity is linked to plant species richness. However, a conflicting spectral diversity–biodiversity relationship in grasslands has been found in previous studies. In particular, it remains unclear how well the spectral diversity–biodiversity relationship holds in naturally assembled species‐rich grasslands. To address the linkage between spectral diversity and plant species richness in a species‐rich alpine grassland ecosystem, we investigated (i) the trade‐off between spectral and spatial resolution in remote sensing data; (ii) the suitability of three different spectral metrics to describe spectral diversity (coefficient of variation, convex hull volume and spectral species richness) and (iii) the importance of confounding effects of live plant biomass, dead plant biomass and plant life forms on the spectral diversity–biodiversity relationship. We addressed these questions using remote sensing data collected with consumer‐grade cameras with four spectral bands and 10 cm spatial resolution on an unmanned aerial vehicle (UAV), airborne imaging spectrometer data (AVIRIS‐NG) with 372 bands and 2.5 m spatial resolution, and a fused data product of both datasets. Our findings suggest that a fused dataset can cope with the requirement of both high spatial‐ and spectral resolution to remotely measure biodiversity. However, in contrast to several previous studies, we found a negative correlation between plant species richness and spectral metrics based on the spectral information content (i.e. spectral complexity). The spectral diversity calculated based on the spectral complexity was sensitive to live and dead plant biomass. Overall, our results suggest that remote sensing of plant species diversity requires a high spatial resolution, the use of classification‐based spectral metrics, such as spectral species richness, and awareness of confounding factors (e.g. plant biomass), which may be ecosystem specific.

SummaryEcological theory posits that temporal stability patterns in plant populations are associated with differences in species’ ecological strategies. However, empirical evidence is lacking about which traits, or trade-offs, underlie species stability, specially across different ecosystems.To address this, we compiled a global collection of long-term permanent vegetation records (>7000 plots from 78 datasets) from a wide range of habitats and combined this with existing trait databases. We tested whether the observed inter-annual variability in species abundance (coefficient of variation) was related to multiple individual traits and multivariate axes of trait variations (PCoA axes).We found that species with greater leaf dry matter content and seed mass were consistently more stable over time (lower variability in species abundance) although other leaf traits played a significant role as well, albeit weaker. Using multivariate axes did not improve predictions by specific traits.Our results confirm existing theory, providing compelling empirical evidence on the importance of specific traits, which point at ecological trade-offs in different resource use and dispersal strategies, on the stability of plant populations worldwide.

Nutrients and herbivores are well-known drivers of grassland diversity and stability in local communities. However, whether they interact to impact the stability of aboveground biomass and whether these effects depend on spatial scales remain unknown. It is also unclear whether nutrients and herbivores impact stability via different facets of plant diversity including species richness, evenness, and changes in community composition through time and space. We used a replicated experiment adding nutrients and excluding herbivores for 5 years in 34 global grasslands to explore these questions. We found that both nutrient addition and herbivore exclusion alone reduced stability at the larger spatial scale (aggregated local communities; gamma stability), but through different pathways. Nutrient addition reduced gamma stability primarily by increasing changes in local community composition over time, which was mainly driven by species replacement. Herbivore exclusion reduced gamma stability primarily by decreasing asynchronous dynamics among local communities (spatial asynchrony). Their interaction weakly increased gamma stability by increasing spatial asynchrony. Our findings indicate that disentangling the processes operating at different spatial scales may improve conservation and management aiming at maintaining the ability of ecosystems to reliably provide functions and services for humanity.

Fertilisation experiments have demonstrated that nutrient availability is a key determinant of biomass production and carbon sequestration in grasslands. However, the influence of nutrients in explaining spatial variation in grassland biomass production has rarely been assessed. Using a global dataset comprising 72 sites on six continents, we investigated which of 16 soil factors that shape nutrient availability associate most strongly with variation in grassland aboveground biomass. Climate and N deposition were also considered. Based on theory-driven structural equation modelling, we found that soil micronutrients (particularly Zn and Fe) were important predictors of biomass and, together with soil physicochemical properties and C:N, they explained more unique variation (32%) than climate and N deposition (24%). However, the association between micronutrients and biomass was absent in grasslands limited by NP. These results highlight soil properties as key predictors of global grassland biomass production and point to serial co-limitation by NP and micronutrients.

Human‐caused loss of vertebrate and invertebrate animals, defaunation, is increasing, and potentially affects plant community structure of diverse grassland ecosystems world‐wide. We experimentally simulated defaunation using size‐selective fences to progressively exclude large‐, medium‐ and small‐sized mammals, and invertebrates from two subalpine vegetation types in the Swiss National Park (SNP): intensively grazed short‐grass and moderately grazed tall‐grass vegetation. We assessed plant community properties yearly from 2009 to 2013, and examined treatment effects on plant community structure in the two grassland types. In the short‐grass vegetation, the exclusion of large mammals increased total plant biomass, while the exclusion of large and medium‐sized mammals increased total, grass and forb biomass compared to when all animals had access. These increases became stronger when invertebrates were also excluded. The exclusion of all mammals and invertebrates increased biomass of grasses by 205%, forbs by 100% and total plant biomass by 118% compared to when all animals had access, hence enhancing relative biomass of grasses from 43.6% to 60%, changing plant species composition and lowering richness of forbs by 16%, the number of plant families by 13% and family‐level Shannon diversity by 23%. In contrast to these significant community‐level responses found in the short‐grass vegetation, there was no evidence that the size‐selective exclusion of animals altered the plant community structure of the tall‐grass vegetation. The contrasting results were due to the difference in plant community composition prior to our experiment, which were related to differences in quantity and quality of forage and in grazing intensities of herbivores between the two grassland types. Synthesis. Our results showed that different‐sized animals, in particular large mammals and invertebrates, contributed to maintain the plant community structure in the short‐grass vegetation, highlighting the importance of multiple, functionally different animal groups for ecosystem functioning and stability. In contrast to the short‐grass vegetation, we could not detect such a top‐down control by animals in the tall‐grass vegetation. Our results suggest that potential defaunation effects on grassland plant community structure depend on the degree of grazing pressure release and grassland vegetation type.

Interannual variability in grassland primary production is strongly driven by precipitation, nutrient availability and herbivory, but there is no general consensus on the mechanisms linking these variables. If grassland biomass is limited by the single most limiting resource at a given time, then we expect that nutrient addition will not affect biomass production at arid sites. We conducted a distributed experiment manipulating nutrients and herbivores at 44 grassland sites in 8 regions around the world, spanning a broad range in aridity. We estimated the effects of 5-11 years of nutrient addition and herbivore exclusion treatments on precipitation sensitivity of biomass (proportional change in biomass relative to proportional change in rainfall among years), and the biomass in the driest year (to measure treatment effects when water was most limiting) at each site. Grazer exclusion did not interact with nutrients to influence driest year biomass or sensitivity. Nutrient addition increased driest year biomass by 74% and sensitivity by 0.12 (proportional units), and that effect did not change across the range of aridity spanned by our sites. Grazer exclusion did not interact with nutrients to influence sensitivity or driest year biomass. At almost half of our sites, the previous year's rainfall explained as much variation in biomass as current year precipitation. Overall, our distributed fertilization experiment detected co-limitation between nutrients and water governing grasslands, with biomass sensitivity to precipitation being limited by nutrient availability irrespective of site aridity and herbivory. Our findings refute the classical ideas that grassland plant performance is limited by the single most limiting resource at a site. This suggests that nutrient eutrophication will destabilize grassland ecosystems through increased sensitivity to precipitation variation.

Anthropogenic nutrient enrichment is driving global biodiversity decline and modifying ecosystem functions. Theory suggests that plant functional types that fix atmospheric nitrogen have a competitive advantage in nitrogen-poor soils, but lose this advantage with increasing nitrogen supply. By contrast, the addition of phosphorus, potassium, and other nutrients may benefit such species in lownutrient environments by enhancing their nitrogen-fixing capacity. We present a global-scale experiment confirming these predictions for nitrogen-fixing legumes (Fabaceae) across 45 grasslands on six continents. Nitrogen addition reduced legume cover, richness, and biomass, particularly in nitrogen-poor soils, while cover of non-nitrogenfixing plants increased. The addition of phosphorous, potassium, and other nutrients enhanced legume abundance, but did not mitigate the negative effects of nitrogen addition. Increasing nitrogen supply thus has the potential to decrease the diversity and abundance of grassland legumes worldwide regardless of the availability of other nutrients, with consequences for biodiversity, food webs, ecosystem resilience, and genetic improvement of protein-rich agricultural plant species.