Negative effects of nitrogen override positive effects of phosphorus on grassland legumes worldwide.

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 low-nutrient 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-nitrogen-fixing 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.

Grassland ecosystem type drives AM fungal diversity and functional guild distribution in North American grasslands.

Nutrient exchange forms the basis of the ancient symbiotic relationship that occurs between most land plants and arbuscular mycorrhizal (AM) fungi. Plants provide carbon (C) to AM fungi and fungi provide the plant with nutrients such as nitrogen (N) and phosphorous (P). Nutrient addition can alter this symbiotic coupling in key ways, such as reducing AM fungal root colonization and changing the AM fungal community composition. However, environmental parameters that differentiate ecosystems and drive plant distribution patterns (e.g., pH, moisture), are also known to impact AM fungal communities. Identifying the relative contribution of environmental factors impacting AM fungal distribution patterns is important for predicting biogeochemical cycling patterns and plant-microbe relationships across ecosystems. To evaluate the relative impacts of local environmental conditions and long-term nutrient addition on AM fungal abundance and composition across grasslands, we studied experimental plots amended for 10 years with N, P, or N and P fertilizer in different grassland ecosystem types, including tallgrass prairie, montane, shortgrass prairie, and desert grasslands. Contrary to our hypothesis, we found ecosystem type, not nutrient treatment, was the main driver of AM fungal root colonization, diversity, and community composition, even when accounting for site-specific nutrient limitations. We identified several important environmental drivers of grassland ecosystem AM fungal distribution patterns, including aridity, mean annual temperature, root moisture, and soil pH. This work provides empirical evidence for niche partitioning strategies of AM fungal functional guilds and emphasizes the importance of long-term, large scale research projects to provide ecologically relevant context to nutrient addition studies.

Nutrient identity modifies the destabilising effects of eutrophication in grasslands.

Nutrient enrichment can simultaneously increase and destabilise plant biomass production, with co-limitation by multiple nutrients potentially intensifying these effects. Here, we test how factorial additions of nitrogen (N), phosphorus (P) and potassium with essential nutrients (K+) affect the stability (mean/standard deviation) of aboveground biomass in 34 grasslands over 7 years. Destabilisation with fertilisation was prevalent but was driven by single nutrients, not synergistic nutrient interactions. On average, N-based treatments increased mean biomass production by 21-51% but increased its standard deviation by 40-68% and so consistently reduced stability. Adding P increased interannual variability and reduced stability without altering mean biomass, while K+ had no general effects. Declines in stability were largest in the most nutrient-limited grasslands, or where nutrients reduced species richness or intensified species synchrony. We show that nutrients can differentially impact the stability of biomass production, with N and P in particular disproportionately increasing its interannual variability.

Linking changes in species composition and biomass in a globally distributed grassland experiment.

Global change drivers, such as anthropogenic nutrient inputs, are increasing globally. Nutrient deposition simultaneously alters plant biodiversity, species composition and ecosystem processes like aboveground biomass production. These changes are underpinned by species extinction, colonisation and shifting relative abundance. Here, we use the Price equation to quantify and link the contributions of species that are lost, gained or that persist to change in aboveground biomass in 59 experimental grassland sites. Under ambient (control) conditions, compositional and biomass turnover was high, and losses (i.e. local extinctions) were balanced by gains (i.e. colonisation). Under fertilisation, the decline in species richness resulted from increased species loss and decreases in species gained. Biomass increase under fertilisation resulted mostly from species that persist and to a lesser extent from species gained. Drivers of ecological change can interact relatively independently with diversity, composition and ecosystem processes and functions such as aboveground biomass due to the individual contributions of species lost, gained or persisting.

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

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.

Soil carbon stocks in temperate grasslands differ strongly across sites but are insensitive to decade-long fertilization.

Enhancing soil carbon (C) storage has the potential to offset human-caused increases in atmospheric CO2 . Rising CO2 has occurred concurrently with increasing supply rates of biologically limiting nutrients such as nitrogen (N) and phosphorus (P). However, it is unclear how increased supplies of N and P will alter soil C sequestration, particularly in grasslands, which make up nearly a third of non-agricultural land worldwide. Here, we leverage a globally distributed nutrient addition experiment (the Nutrient Network) to examine how a decade of N and P fertilization (alone and in combination) influenced soil C and N stocks at nine grassland sites spanning the continental United States. We measured changes in bulk soil C and N stocks and in three soil C fractions (light and heavy particulate organic matter, and mineral-associated organic matter fractions). Nutrient amendment had variable effects on soil C and N pools that ranged from strongly positive to strongly negative, while soil C and N pool sizes varied by more than an order of magnitude across sites. Piecewise SEM clarified that small increases in plant C inputs with fertilization did not translate to greater soil C storage. Nevertheless, peak season aboveground plant biomass (but not root biomass or production) was strongly positively related to soil C storage at seven of the nine sites, and across all nine sites, soil C covaried with moisture index and soil mineralogy, regardless of fertilization. Overall, we show that site factors such as moisture index, plant productivity, soil texture, and mineralogy were key predictors of cross-site soil C, while nutrient amendment had weaker and site-specific effects on C sequestration. This suggests that prioritizing the protection of highly productive temperate grasslands is critical for reducing future greenhouse gas losses arising from land use change.

Soil properties as key predictors of global grassland production: Have we overlooked micronutrients?

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.

Fertilized graminoids intensify negative drought effects on grassland productivity.

Droughts can strongly affect grassland productivity and biodiversity, but responses differ widely. Nutrient availability may be a critical factor explaining this variation, but is often ignored in analyses of drought responses. Here, we used a standardized nutrient addition experiment covering 10 European grasslands to test if full-factorial nitrogen, phosphorus, and potassium addition affected plant community responses to inter-annual variation in drought stress and to the extreme summer drought of 2018 in Europe. We found that nutrient addition amplified detrimental drought effects on community aboveground biomass production. Drought effects also differed between functional groups, with a negative effect on graminoid but not forb biomass production. Our results imply that eutrophication in grasslands, which promotes dominance of drought-sensitive graminoids over forbs, amplifies detrimental drought effects. In terms of climate change adaptation, agricultural management would benefit from taking into account differential drought impacts on fertilized versus unfertilized grasslands, which differ in ecosystem services they provide to society.

Nutrient availability controls the impact of mammalian herbivores on soil carbon and nitrogen pools in grasslands.

Grasslands are subject to considerable alteration due to human activities globally, including widespread changes in populations and composition of large mammalian herbivores and elevated supply of nutrients. Grassland soils remain important reservoirs of carbon (C) and nitrogen (N). Herbivores may affect both C and N pools and these changes likely interact with increases in soil nutrient availability. Given the scale of grassland soil fluxes, such changes can have striking consequences for atmospheric C concentrations and the climate. Here, we use the Nutrient Network experiment to examine the responses of soil C and N pools to mammalian herbivore exclusion across 22 grasslands, under ambient and elevated nutrient availabilities (fertilized with NPK + micronutrients). We show that the impact of herbivore exclusion on soil C and N pools depends on fertilization. Under ambient nutrient conditions, we observed no effect of herbivore exclusion, but under elevated nutrient supply, pools are smaller upon herbivore exclusion. The highest mean soil C and N pools were found in grazed and fertilized plots. The decrease in soil C and N upon herbivore exclusion in combination with fertilization correlated with a decrease in aboveground plant biomass and microbial activity, indicating a reduced storage of organic matter and microbial residues as soil C and N. The response of soil C and N pools to herbivore exclusion was contingent on temperature - herbivores likely cause losses of C and N in colder sites and increases in warmer sites. Additionally, grasslands that contain mammalian herbivores have the potential to sequester more N under increased temperature variability and nutrient enrichment than ungrazed grasslands. Our study highlights the importance of conserving mammalian herbivore populations in grasslands worldwide. We need to incorporate local-scale herbivory, and its interaction with nutrient enrichment and climate, within global-scale models to better predict land-atmosphere interactions under future climate change.

Sensitivity of global soil carbon stocks to combined nutrient enrichment.

Soil stores approximately twice as much carbon as the atmosphere and fluctuations in the size of the soil carbon pool directly influence climate conditions. We used the Nutrient Network global change experiment to examine how anthropogenic nutrient enrichment might influence grassland soil carbon storage at a global scale. In isolation, enrichment of nitrogen and phosphorous had minimal impacts on soil carbon storage. However, when these nutrients were added in combination with potassium and micronutrients, soil carbon stocks changed considerably, with an average increase of 0.04 KgCm-2  year-1 (standard deviation 0.18 KgCm-2  year-1 ). These effects did not correlate with changes in primary productivity, suggesting that soil carbon decomposition may have been restricted. Although nutrient enrichment caused soil carbon gains most dry, sandy regions, considerable absolute losses of soil carbon may occur in high-latitude regions that store the majority of the world's soil carbon. These mechanistic insights into the sensitivity of grassland carbon stocks to nutrient enrichment can facilitate biochemical modelling efforts to project carbon cycling under future climate scenarios.

More salt, please: global patterns, responses and impacts of foliar sodium in grasslands.

Sodium is unique among abundant elemental nutrients, because most plant species do not require it for growth or development, whereas animals physiologically require sodium. Foliar sodium influences consumption rates by animals and can structure herbivores across landscapes. We quantified foliar sodium in 201 locally abundant, herbaceous species representing 32 families and, at 26 sites on four continents, experimentally manipulated vertebrate herbivores and elemental nutrients to determine their effect on foliar sodium. Foliar sodium varied taxonomically and geographically, spanning five orders of magnitude. Site-level foliar sodium increased most strongly with site aridity and soil sodium; nutrient addition weakened the relationship between aridity and mean foliar sodium. Within sites, high sodium plants declined in abundance with fertilisation, whereas low sodium plants increased. Herbivory provided an explanation: herbivores selectively reduced high nutrient, high sodium plants. Thus, interactions among climate, nutrients and the resulting nutritional value for herbivores determine foliar sodium biogeography in herbaceous-dominated systems.

Spatial heterogeneity in species composition constrains plant community responses to herbivory and fertilisation.

Environmental change can result in substantial shifts in community composition. The associated immigration and extinction events are likely constrained by the spatial distribution of species. Still, studies on environmental change typically quantify biotic responses at single spatial (time series within a single plot) or temporal (spatial beta diversity at single time points) scales, ignoring their potential interdependence. Here, we use data from a global network of grassland experiments to determine how turnover responses to two major forms of environmental change - fertilisation and herbivore loss - are affected by species pool size and spatial compositional heterogeneity. Fertilisation led to higher rates of local extinction, whereas turnover in herbivore exclusion plots was driven by species replacement. Overall, sites with more spatially heterogeneous composition showed significantly higher rates of annual turnover, independent of species pool size and treatment. Taking into account spatial biodiversity aspects will therefore improve our understanding of consequences of global and anthropogenic change on community dynamics.

Herbivory and eutrophication mediate grassland plant nutrient responses across a global climatic gradient.

Plant stoichiometry, the relative concentration of elements, is a key regulator of ecosystem functioning and is also being altered by human activities. In this paper we sought to understand the global drivers of plant stoichiometry and compare the relative contribution of climatic vs. anthropogenic effects. We addressed this goal by measuring plant elemental (C, N, P and K) responses to eutrophication and vertebrate herbivore exclusion at eighteen sites on six continents. Across sites, climate and atmospheric N deposition emerged as strong predictors of plot-level tissue nutrients, mediated by biomass and plant chemistry. Within sites, fertilization increased total plant nutrient pools, but results were contingent on soil fertility and the proportion of grass biomass relative to other functional types. Total plant nutrient pools diverged strongly in response to herbivore exclusion when fertilized; responses were largest in ungrazed plots at low rainfall, whereas herbivore grazing dampened the plant community nutrient responses to fertilization. Our study highlights (1) the importance of climate in determining plant nutrient concentrations mediated through effects on plant biomass, (2) that eutrophication affects grassland nutrient pools via both soil and atmospheric pathways and (3) that interactions among soils, herbivores and eutrophication drive plant nutrient responses at small scales, especially at water-limited sites.

Increased grassland arthropod production with mammalian herbivory and eutrophication: a test of mediation pathways.

Increases in nutrient availability and alterations to mammalian herbivore communities are a hallmark of the Anthropocene, with consequences for the primary producer communities in many ecosystems. While progress has advanced understanding of plant community responses to these perturbations, the consequences for energy flow to higher trophic levels in the form of secondary production are less well understood. We quantified arthropod biomass after manipulating soil nutrient availability and wild mammalian herbivory, using identical methods across 13 temperate grasslands. Of experimental increases in nitrogen, phosphorus, and potassium, only treatments including nitrogen resulted in significantly increased arthropod biomass. Wild mammalian herbivore removal had a marginal, negative effect on arthropod biomass, with no interaction with nutrient availability. Path analysis including all sites implicated nutrient content of the primary producers as a driver of increased arthropod mean size, which we confirmed using 10 sites for which we had foliar nutrient data. Plant biomass and physical structure mediated the increase in arthropod abundance, while the nitrogen treatments accounted for additional variation not explained by our measured plant variables. The mean size of arthropod individuals was 2.5 times more influential on the plot-level total arthropod biomass than was the number of individuals. The eutrophication of grasslands through human activity, especially nitrogen deposition, thus may contribute to higher production of arthropod consumers through increases in nutrient availability across trophic levels.

Joint effects of nutrient addition and enemy exclusion on exotic plant success.

Worldwide, ecosystems are increasingly dominated by exotic plant species, a shift hypothesized to result from numerous ecological factors. Two of these, increased resource availability and enemy release, may act in concert to increase exotic success in plant communities (Resource-Enemy Release Hypothesis, R-ERH). To test this, we manipulated the availability of soil nutrients and access of vertebrate herbivores, insect herbivores, and fungal pathogens to intact grassland communities containing both native and exotic species. Our results supported both conditions necessary for R-ERH. First, exotics were less damaged than natives, experiencing less foliar damage (insect herbivory and fungal disease) than native species, particularly in communities where soil nutrients were added. Second, fertilization increased foliar damage on native species, but not exotic species. As well as fulfilling both conditions for R-ERH, these results demonstrate the importance of considering the effects of resource availability when testing for enemy release. When both conditions are fulfilled, R-ERH predicts that increasing resource availability will increase exotic abundance only in the presence of enemies. Our results fully supported this prediction for vertebrate herbivores: fertilization increased exotic cover only in communities exposed to vertebrate herbivores. Additionally, the prediction was partially supported for insect herbivores and fungal pathogens, excluding these enemies reduced exotic cover as predicted, but inconsistent with R-ERH, this effect occurred only in unfertilized communities. These results highlight the need to consider the influence of multiple enemy guilds on community processes like exotic plant invasions. Moreover, this study experimentally demonstrates that resource availability and natural enemies can jointly influence exotic success in plant communities.

Stronger fertilization effects on aboveground versus belowground plant properties across nine U.S. grasslands.

Increased nutrient inputs due to anthropogenic activity are expected to increase primary productivity across terrestrial ecosystems, but changes in allocation aboveground versus belowground with nutrient addition have different implications for soil carbon (C) storage. Thus, given that roots are major contributors to soil C storage, understanding belowground net primary productivity (BNPP) and biomass responses to changes in nutrient availability is essential to predicting carbon-climate feedbacks in the context of interacting global environmental changes. To address this knowledge gap, we tested whether a decade of nitrogen (N) and phosphorus (P) fertilization consistently influenced aboveground and belowground biomass and productivity at nine grassland sites spanning a wide range of climatic and edaphic conditions in the continental United States. Fertilization effects were strong aboveground, with both N and P addition stimulating aboveground biomass at nearly all sites (by 30% and 36%, respectively, on average). P addition consistently increased root production (by 15% on average), whereas other belowground responses to fertilization were more variable, ranging from positive to negative across sites. Site-specific responses to P were not predicted by the measured covariates. Atmospheric N deposition mediated the effect of N fertilization on root biomass and turnover. Specifically, atmospheric N deposition was positively correlated with root turnover rates, and this relationship was amplified with N addition. Nitrogen addition increased root biomass at sites with low N deposition but decreased it at sites with high N deposition. Overall, these results suggest that the effects of nutrient supply on belowground plant properties are context dependent, particularly with regard to background N supply rates, demonstrating that site conditions must be considered when predicting how grassland ecosystems will respond to increased nutrient loading from anthropogenic activity.

Nitrogen deposition and climate: an integrated synthesis.

Human activities have more than doubled reactive nitrogen (N) deposited in ecosystems, perturbing the N cycle and considerably impacting plant, animal, and microbial communities. However, biotic responses to N deposition can vary widely depending on factors including local climate and soils, limiting our ability to predict ecosystem responses. Here, we synthesize reported impacts of elevated N on grasslands and draw upon evidence from the globally distributed Nutrient Network experiment (NutNet) to provide insight into causes of variation and their relative importance across scales. This synthesis highlights that climate and elevated N frequently interact, modifying biotic responses to N. It also demonstrates the importance of edaphic context and widespread interactions with other limiting nutrients in controlling biotic responses to N deposition.

Global Grassland Diazotrophic Communities Are Structured by Combined Abiotic, Biotic, and Spatial Distance Factors but Resilient to Fertilization.

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 (N2) fixation represents an essential natural source of nitrogen (N). The ability to fix atmospheric N2 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.

Increasing effects of chronic nutrient enrichment on plant diversity loss and ecosystem productivity over time.

Human activities are enriching many of Earth's ecosystems with biologically limiting mineral nutrients such as nitrogen (N) and phosphorus (P). In grasslands, this enrichment generally reduces plant diversity and increases productivity. The widely demonstrated positive effect of diversity on productivity suggests a potential negative feedback, whereby nutrient-induced declines in diversity reduce the initial gains in productivity arising from nutrient enrichment. In addition, plant productivity and diversity can be inhibited by accumulations of dead biomass, which may be altered by nutrient enrichment. Over longer time frames, nutrient addition may increase soil fertility by increasing soil organic matter and nutrient pools. We examined the effects of 5-11 yr of nutrient addition at 47 grasslands in 12 countries. Nutrient enrichment increased aboveground live biomass and reduced plant diversity at nearly all sites, and these effects became stronger over time. We did not find evidence that nutrient-induced losses of diversity reduced the positive effects of nutrients on biomass; however, nutrient effects on live biomass increased more slowly at sites where litter was also increasing, regardless of plant diversity. This work suggests that short-term experiments may underestimate the long-term nutrient enrichment effects on global grassland ecosystems.

Opposing community assembly patterns for dominant and nondominant plant species in herbaceous ecosystems globally.

Biotic and abiotic factors interact with dominant plants-the locally most frequent or with the largest coverage-and nondominant plants differently, partially because dominant plants modify the environment where nondominant plants grow. For instance, if dominant plants compete strongly, they will deplete most resources, forcing nondominant plants into a narrower niche space. Conversely, if dominant plants are constrained by the environment, they might not exhaust available resources but instead may ameliorate environmental stressors that usually limit nondominants. Hence, the nature of interactions among nondominant species could be modified by dominant species. Furthermore, these differences could translate into a disparity in the phylogenetic relatedness among dominants compared to the relatedness among nondominants. By estimating phylogenetic dispersion in 78 grasslands across five continents, we found that dominant species were clustered (e.g., co-dominant grasses), suggesting dominant species are likely organized by environmental filtering, and that nondominant species were either randomly assembled or overdispersed. Traits showed similar trends for those sites (<50%) with sufficient trait data. Furthermore, several lineages scattered in the phylogeny had more nondominant species than expected at random, suggesting that traits common in nondominants are phylogenetically conserved and have evolved multiple times. We also explored environmental drivers of the dominant/nondominant disparity. We found different assembly patterns for dominants and nondominants, consistent with asymmetries in assembly mechanisms. Among the different postulated mechanisms, our results suggest two complementary hypotheses seldom explored: (1) Nondominant species include lineages adapted to thrive in the environment generated by dominant species. (2) Even when dominant species reduce resources to nondominant ones, dominant species could have a stronger positive effect on some nondominants by ameliorating environmental stressors affecting them, than by depleting resources and increasing the environmental stress to those nondominants. These results show that the dominant/nondominant asymmetry has ecological and evolutionary consequences fundamental to understand plant communities.