Fertilized graminoids intensify negative drought effects on grassland productivity

Assessing vulnerability to embolism and hydraulic safety margins in reed-like Restionaceae

The African Restionaceae (Poales), the dominant graminoid layer in the megadiverse Cape Floristic Region of South Africa, are distributed across a wide range of moisture availability, yet currently there is very little known about the underlying hydraulics of this group. We tested two methods for measuring culm vulnerability to embolism, the optical and pneumatic methods, in three species of Cannomois ranging in habitat from semi-riparian (Cannomois virgata) to dryland (Cannomois parviflora and C. congesta). Estimates of culm xylem vulnerability were coupled with measures of turgor loss point (ΨTLP) and minimum field water potential (ΨMD) to assess hydraulic safety margins. The optical and pneumatic methods produced similar estimates of P50, but differed for P12 and P88. All three species were quite vulnerable to embolism, with P50 of -1.9 MPa (C. virgata), -2.3 MPa (C. congesta), and -2.4 MPa (C. parviflora). Estimates of P50, ΨTLP and ΨMD aligned with habitat moisture stress, with highest values found in the semi-riparian C. virgata. Consistent differences in P50, ΨMD and ΨTLP between species resulted in consistent hydraulic safety margins across species of 0.96 ± 0.1 MPa between ΨMD and P50, with onset of embolism occurring 0.43 ± 0.04 MPa after ΨTLP for all three species. Our study demonstrates that restio occupancy of dry environments involves more than the evolution of highly resistant xylem, suggesting that other aspects of water relations are key to understanding trait-environment relationships in this group.

Effects of Short-Term Nitrogen Additions on Biomass and Soil Phytochemical Cycling in Alpine Grasslands of Tianshan, China

The nitrogen deposition process, as an important phenomenon of global climate change and an important link in the nitrogen cycle, has had serious and far-reaching impacts on grassland ecosystems. This study aimed to investigate the survival adaptation strategies of plants of different functional groups under nitrogen deposition, and the study identified the following outcomes of differences in biomass changes by conducting in situ simulated nitrogen deposition experiments while integrating plant nutrient contents and soil physicochemical properties: (1) nitrogen addition enhanced the aboveground biomass of grassland communities, in which Poaceae were significantly affected by nitrogen addition. Additionally, nitrogen addition significantly influenced plant total nitrogen and total phosphorus; (2) nitrogen addition improved the plant growth environment, alleviated plant nitrogen limitation, and promoted plant phosphorus uptake; and (3) there was variability in the biomass responses of different functional groups to nitrogen addition. The level of nitrogen addition was the primary factor affecting differences in biomass changes, while nitrogen addition frequency was an important factor affecting changes in plant community structure.

Can permanent grassland soils with elevated organic carbon buffer negative effects of more persistent precipitation regimes on forage grass performance?

Agricultural practices enhancing soil organic carbon (SOC) show potential to buffer negative effects of climate change on forage grass performance. We tested this by subjecting five forage grass varieties differing in fodder quality and drought/flooding resistance to increased persistence in summer precipitation regimes (PR) across sandy and sandy-loam soils from either permanent (high SOC) or temporary grasslands (low SOC) in adjacent parcels. Over the course of two consecutive summers, monoculture mesocosms were subjected to rainy/dry weather alternation either every 3 days or every 30 days, whilst keeping total precipitation equal. Increased PR persistence induced species-specific drought damage and productivity declines. Soils from permanent grasslands with elevated SOC buffered plant quality, but buffering effects of SOC on drought damage, nutrient availability and yield differed between texture classes. In the more persistent PR, Festuca arundinacea FERMINA was the most productive species but had the lowest quality under both ample water supply and mild soil drought, whilst under the most intense soil droughts, Festulolium FESTILO maintained the highest yields. The hybrid Lolium × boucheanum kunth MELCOMBI had intermediate productivity and both Lolium perenne varieties showed the lowest yields under soil drought, but the highest forage quality (especially the tetraploid variety MELFORCE). Performance varied with plant maturity stage and across seasons/years and was driven by altered water and nutrient availability and related nitrogen nutrition among species during drought and upon rewetting. Moreover, whilst permanent grassland soils showed the most consistent positive effects on plant performance, their available water capacity also declined under increased PR persistence. We conclude that permanent grassland soils with historically elevated SOC likely buffer negative effects of increasing summer weather persistence on forage grass performance, but may also be more sensitive to degradation under climate change.

Cumulative nitrogen enrichment alters the drivers of grassland overyielding

Effects of plant diversity on grassland productivity, or overyielding, are found to be robust to nutrient enrichment. However, the impact of cumulative nitrogen (N) addition (total N added over time) on overyielding and its drivers are underexplored. Synthesizing data from 15 multi-year grassland biodiversity experiments with N addition, we found that N addition decreases complementarity effects and increases selection effects proportionately, resulting in no overall change in overyielding regardless of N addition rate. However, we observed a convex relationship between overyielding and cumulative N addition, driven by a shift from complementarity to selection effects. This shift suggests diminishing positive interactions and an increasing contribution of a few dominant species with increasing N accumulation. Recognizing the importance of cumulative N addition is vital for understanding its impacts on grassland overyielding, contributing essential insights for biodiversity conservation and ecosystem resilience in the face of increasing N deposition.

Plant diversity and community age stabilize ecosystem multifunctionality

It is well known that biodiversity positively affects ecosystem functioning, leading to enhanced ecosystem stability. However, this knowledge is mainly based on analyses using single ecosystem functions, while studies focusing on the stability of ecosystem multifunctionality (EMF) are rare. Taking advantage of a long-term grassland biodiversity experiment, we studied the effect of plant diversity (1-60 species) on EMF over 5 years, its temporal stability, as well as multifunctional resistance and resilience to a 2-year drought event. Using split-plot treatments, we further tested whether a shared history of plants and soil influences the studied relationships. We calculated EMF based on functions related to plants and higher-trophic levels. Plant diversity enhanced EMF in all studied years, and this effect strengthened over the study period. Moreover, plant diversity increased the temporal stability of EMF and fostered resistance to reoccurring drought events. Old plant communities with shared plant and soil history showed a stronger plant diversity-multifunctionality relationship and higher temporal stability of EMF than younger communities without shared histories. Our results highlight the importance of old and biodiverse plant communities for EMF and its stability to extreme climate events in a world increasingly threatened by global change.

Herb hydraulics: Variation and correlation for traits governing drought tolerance and efficiency of water transport

Hydraulic traits dictate plant response to drought, thus enabling better understanding of community dynamics under global climate change. Despite being intensively documented in woody species, herbaceous species (graminoids and forbs) are largely understudied, hence the distribution and correlation of hydraulic traits in herbaceous species remains unclear. Here, we collected key hydraulic traits for 436 herbaceous species from published literature, including leaf hydraulic conductivity (Kleaf), water potential inducing 50 % loss of hydraulic conductivity (P50), stomatal closure (Pclose) and turgor loss (Ptlp). Trait variation of herbs was analyzed and contrasted with angiosperm woody species within the existing global hydraulic traits database, as well as between different growth forms within herbs. Furthermore, hydraulic traits coordination was also assessed for herbaceous species. We found that herbs showed overall more negative Pclose but less negative Ptlp compared with angiosperm woody species, while P50 did not differ between functional types, regardless of the organ (leaf and stem). In addition, correlations were found between Kleaf and P50 of leaf (P50leaf), as well as between Pclose, P50leaf and Kleaf. Within herbs, graminoids generally exhibited more negative P50 and Ptlp, but lower Kleaf, relative to forbs. Within herbs, no clear pattern regarding hydraulic traits-climate relationship was found. Our analysis provided insights into herb hydraulic, and highlighted the knowledge gaps need to be filled regarding the response of herbs to drought.

Effects of Spring Drought and Nitrogen Addition on Productivity and Community Composition of Degraded Grasslands

To explore whether there were differences among the patterns of response of grasslands with different levels of degradation to extreme drought events and nitrogen addition, three grasslands along a degradation gradient (extremely, moderately, and lightly degraded) were selected in the Bashang area of northern China using the human disturbance index (HDI). A field experiment with simulated extreme spring drought, nitrogen addition, and their interaction was conducted during the growing seasons of 2020 and 2021. The soil moisture, aboveground biomass, and composition of the plant community were measured. The primary results were as follows. (1) Drought treatment caused soil drought stress, with moderately degraded grassland being the most affected, which resulted in an 80% decrease in soil moisture and a 78% decrease in aboveground biomass. The addition of nitrogen did not mitigate the impact of drought. Moreover, the aboveground net primary production (ANPP) in 2021 was less sensitive to spring drought than in 2020. (2) The community composition changed after 2 years of drought treatment, particularly for the moderately degraded grasslands with annual forbs, such as Salsola collina, increasing significantly in biomass proportion, which led to a trend of exacerbated degradation (higher HDI). This degradation trend decreased under the addition of nitrogen. (3) The variation in drought sensitivities of the ANPP was primarily determined by the proportion of plants based on the classification of degradation indicators in the community, with higher proportions of intermediate degradation indicator species exhibiting more sensitivity to spring drought. These findings can help to provide scientific evidence for the governance and restoration of regional degraded grassland under frequent extreme weather conditions.

Nutrient effects on drought responses vary across common temperate grassland species

Drought and nutrient input are two main global change drivers that threaten ecosystem function and services. Resolving the interactive effects of human-induced stressors on individual species is necessary to improve our understanding of community and ecosystem responses. This study comparatively assessed how different nutrient conditions affect whole-plant drought responses across 13 common temperate grassland species. We conducted a fully factorial drought-fertilization experiment to examine the effect of nutrient addition [nitrogen (N), phosphorus (P), and combined NP] on species' drought survival, and on drought resistance of growth as well as drought legacy effects. Drought had an overall negative effect on survival and growth, and the adverse drought effects extended into the next growing season. Neither drought resistance nor legacy effects exhibited an overall effect of nutrients. Instead, both the size and the direction of the effects differed strongly among species and between nutrient conditions. Consistently, species performance ranking under drought changed with nitrogen availability. The idiosyncratic responses of species to drought under different nutrient conditions may underlie the seemingly contradicting effects of drought in studies on grassland composition and productivity along nutrient and land-use gradients-ranging from amplifying to dampening. Differential species' responses to combinations of nutrients and drought, as observed in our study, complicate predictions of community and ecosystem responses to climate and land-use changes. Moreover, they highlight the urgent need for an improved understanding of the mechanisms that render species more or less vulnerable to drought under different nutrients.

When things get MESI: The Manipulation Experiments Synthesis Initiative-A coordinated effort to synthesize terrestrial global change experiments

Responses of the terrestrial biosphere to rapidly changing environmental conditions are a major source of uncertainty in climate projections. In an effort to reduce this uncertainty, a wide range of global change experiments have been conducted that mimic future conditions in terrestrial ecosystems, manipulating CO₂ , temperature, and nutrient and water availability. Syntheses of results across experiments provide a more general sense of ecosystem responses to global change, and help to discern the influence of background conditions such as climate and vegetation type in determining global change responses. Several independent syntheses of published data have yielded distinct databases for specific objectives. Such parallel, uncoordinated initiatives carry the risk of producing redundant data collection efforts and have led to contrasting outcomes without clarifying the underlying reason for divergence. These problems could be avoided by creating a publicly available, updatable, curated database. Here, we report on a global effort to collect and curate 57,089 treatment responses across 3644 manipulation experiments at 1145 sites, simulating elevated CO₂ , warming, nutrient addition, and precipitation changes. In the resulting Manipulation Experiments Synthesis Initiative (MESI) database, effects of experimental global change drivers on carbon and nutrient cycles are included, as well as ancillary data such as background climate, vegetation type, treatment magnitude, duration, and, unique to our database, measured soil properties. Our analysis of the database indicates that most experiments are short term (one or few growing seasons), conducted in the USA, Europe, or China, and that the most abundantly reported variable is aboveground biomass. We provide the most comprehensive multifactor global change database to date, enabling the research community to tackle open research questions, vital to global policymaking. The MESI database, freely accessible at doi.org/10.5281/zenodo.7153253, opens new avenues for model evaluation and synthesis-based understanding of how global change affects terrestrial biomes. We welcome contributions to the database on GitHub.

Assessing multiple limiting factors of seasonal biomass production and N content in a grassland with a year-round production

There is little evidence on the extent that multiple factors simultaneously limit ecosystem function of grasslands with year-round production. Here we test if multiple factors simultaneously limit (i.e., more than one factor at a time) grassland functioning in different seasons and how they interacted with N availability. In a Flooding Pampa grassland, we ran a separate factorial experiment in spring, summer, and winter with several treatments: control, mowing, shading, P addition, watering (only in summer), and warming (only in winter), each of them crossed with two nitrogen treatments: control and N addition. Grassland functioning was assessed by aboveground net primary productivity (ANPP), green and standing dead biomass, and N content at the species group level. Out of 24 potential cases (three seasons by eight response variables), 13 corresponded to just one limiting factor, 4 to multiple limiting factors, and the other 7 to no evidence of limitation. In conclusion, grassland functioning in each season was most often limited by just one factor, while multiple limiting factors were rarer. Nitrogen was the prevailing limiting factor. Our study expands our knowledge of limitations imposed by factors associated with disturbance and stress, such as mowing, shading, water availability, and warming in grasslands with year-round production.

Shifting biomass allocation and light limitation co-regulate the temporal stability of an alpine meadow under eutrophication

Eutrophication generally promotes but destabilizes grassland productivity. Under eutrophication, plants tend to decrease biomass allocation to roots but increase aboveground allocation and light limitation, likely affecting community stability. However, it remains unclear to understand how shifting plant biomass allocation and light limitation regulate grassland stability in response to eutrophication. Here, using a 5-yr multiple nutrient addition experiment in an alpine meadow, we explored the role of changes in plant biomass allocation and light limitation on its community stability under eutrophication as well as traditionally established mechanisms (i.e., plant Shannon diversity, species asynchrony and grass subcommunity stability). Our results showed that nitrogen (N) addition, rather than phosphorus (P) or potassium (K) addition, significantly reduced the temporal stability of the alpine meadow. In accordance with previous studies, we found that N addition decreased plant Shannon diversity, species asynchrony and grass subcommunity stability, further destabilizing meadow community productivity. In addition, we also found the decrease in biomass allocation to belowground by N addition, further weakening its community stability. Moreover, this shifts in plant biomass allocation from below- to aboveground, intensifying plant light limitation. Further, the light limitation reduced plant species asynchrony, which finally weakened its community stability. Overall, in addition to traditionally established mechanisms, this study highlights the role of plant biomass allocation shifting from belowground to aboveground in determining grassland community stability. These "unseen" mechanisms might improve our understanding of grassland stability in the context of ongoing eutrophication.

Low precipitation due to climate change consistently reduces multifunctionality of urban grasslands in mesocosms

Urban grasslands are crucial for biodiversity and ecosystem services in cities, while little is known about their multifunctionality under climate change. Thus, we investigated the effects of simulated climate change, i.e., increased [CO2] and temperature, and reduced precipitation, on individual functions and overall multifunctionality in mesocosm grasslands sown with forbs and grasses in four different proportions aiming at mimicking road verge grassland patches. Climate change scenarios RCP2.6 (control) and RCP8.5 (worst-case) were simulated in walk-in climate chambers of an ecotron facility, and watering was manipulated for normal vs. reduced precipitation. We measured eight indicator variables of ecosystem functions based on below- and aboveground characteristics. The young grassland communities responded to higher [CO2] and warmer conditions with increased vegetation cover, height, flower production, and soil respiration. Lower precipitation affected carbon cycling in the ecosystem by reducing biomass production and soil respiration. In turn, the water regulation capacity of the grasslands depended on precipitation interacting with climate change scenario, given the enhanced water efficiency resulting from increased [CO2] under RCP8.5. Multifunctionality was negatively affected by reduced precipitation, especially under RCP2.6. Trade-offs arose among single functions that performed best in either grass- or forb-dominated grasslands. Grasslands with an even ratio of plant functional types coped better with climate change and thus are good options for increasing the benefits of urban green infrastructure. Overall, the study provides experimental evidence of the effects of climate change on the functionality of urban ecosystems. Designing the composition of urban grasslands based on ecological theory may increase their resilience to global change.

Effects of multi-resource addition on grassland plant productivity and biodiversity along a resource gradient

The change of plant biodiversity caused by resource-enhancing global changes has greatly affected grassland productivity. However, it remains unclear how multi-resource enrichment induces the effects of multifaceted biodiversity on grassland productivity under different site resource constraints. We conducted a multiple resource addition (MRA) experiment of water and nutrients at three sites located along a resource gradient in northern China. This allowed us to assess the response of aboveground net primary productivity (ANPP), species (species richness and plant density), functional (functional richness and community-weighted mean of traits) and phylogenetic (phylogenetic richness) diversity to increasing number of MRA. We used structural equation model (SEM) to examine the direct and indirect effects of MRA and multifaceted biodiversity on ANPP. The combined addition of the four resources increased ANPP at all three sites. But with increasing number of MRA, biodiversity varied at the three sites. At the high resource constraint site, species richness, plant density and leaf nitrogen concentration (LNC) increased. At the medium resource constraint site, plant height and LNC increased, leaf dry matter content (LDMC) decreased. At the low resource constraint site, species, functional and phylogenetic richness decreased, and height increased. The SEM showed that MRA increased ANPP directly at all three sites, and indirectly by increasing plant density at the high constraint site and height at the medium constraint site. Independent of MRA, ANPP was affected by height at the high resource constraint site and LNC at the low resource constraint site. Our results illustrate that multi-resource addition positively affects productivity, while affects biodiversity depending on site resource constraint. The study highlights that site resource constraint conditions need to be taken into consideration to better predict grassland structure and function, particularly under the future multifaceted global change scenarios.

Does previous exposure to extreme precipitation regimes result in acclimated grassland communities?

Climate change will likely increase weather persistence in the mid-latitudes, resulting in precipitation regimes (PR) with longer dry and wet periods compared to historic averages. This could affect terrestrial ecosystems substantially through the increased occurrence of repeated, prolonged drought and water logging conditions. Climate history is an important determinant of ecosystem responses to consecutive environmental extremes, through direct damage, community restructuring as well as morphological and physiological acclimation in species or individuals. However, it is unclear how community restructuring and individual metabolic acclimation effects interact to determine ecosystem responses to subsequent climate extremes. Here, we investigated, if and how, differences in exposure to extreme or historically normal PR induced long-lasting (i.e. legacy) effects at the level of community (e.g., species composition), plant (e.g., biomass), and molecular composition (e.g., sugars, lipids, stress markers). Experimental grassland communities were exposed to long (extreme) or short (historically normal) dry/wet cycles in year 1 (Y1), followed by exposure to an identical PR or the opposite PR in year 2 (Y2). Results indicate that exposure to extreme PR in Y1, reduced diversity but induced apparent acclimation effects in all climate scenarios, stimulating biomass (higher productivity and structural sugar content) in Y2. In contrast, plants pre-exposed to normal PR, showed more activated stress responses (higher proline and antioxidants) under extreme PR in Y2. Overall, Y1 acclimation effects were strongest in the dominant grasses, indicating comparatively high phenotypical plasticity. However, Y2 drought intensity also correlated with grass productivity and structural sugar findings, suggesting that responses to short-term soil water deficits contributed to the observed patterns. Interactions between different legacy effects are discussed. We conclude that more extreme PR will likely alter diversity in the short-to midterm and select for acclimated grassland communities with increased productivity and attenuated molecular stress responses under future climate regimes.

Nitrogen enrichment alters climate sensitivity of biodiversity and productivity differentially and reverses the relationship between them in an alpine meadow

Biodiversity and productivity that highly determine ecosystem services are varying largely under global change. However, the climate sensitivity of them and their relationship are not well understood, especially in the context of increasing nitrogen (N) deposition. Here, based on a six-year N manipulation experiment in an alpine meadow, we quantified interannual climate sensitivity of species richness (SR) and above-ground net primary productivity (ANPP) as well as SR-ANPP relationship as affected by six N addition rate (N_rate) gradients. We found that interannual variations in ANPP and SR were mainly driven by temperature instead of precipitation. In the plots without N addition, higher temperature substantially increased ANPP but reduced SR across years, thus resulting in a negative SR-ANPP relationship. However, the negative and positive responses of SR and ANPP to temperature increased and declined significantly with increasing N_rate, respectively, leading to a shift of the negative relationship between SR and ANPP into a positive one under high N_rate. Moreover, the adverse influence of drought on SR and ANPP would be aggravated by N fertilization, as indicated by the increased positive effect of precipitation on them under N enrichment. Our findings indicate that climate sensitivity of productivity and biodiversity may be misestimated if the impact of N deposition is not considered, and the importance of biodiversity to maintain productivity would enhance as N deposition increases. This study provides a new insight to explain variation of biodiversity-productivity relationship along with environmental changes.

Fallow Land Enhances Carbon Sequestration in Glomalin and Soil Aggregates Through Regulating Diversity and Network Complexity of Arbuscular Mycorrhizal Fungi Under Climate Change in Relatively High-Latitude Regions

Soil aggregation and aggregate-associated carbon (C) play an essential function in soil health and C sequestration. Arbuscular mycorrhizal fungi (AMF) are considered to be primary soil aggregators due to the combined effect of extraradical hyphae and glomalin-related soil proteins (GRSPs). However, the effects of diversity and network complexity of AMF community on stability of soil aggregates and their associated C under long-term climate change (CC) and land-use conversion (LUC) in relatively high-latitude regions are largely unexplored. Therefore, an 8-year soil plot (with a 30-year cropping history) transplantation experiment was conducted to simulate CC and LUC from cropland to fallow land. The results showed that Glomus, Paraglomus, and Archaeospora were the most abundant genera. The diversity of AMF community in fallow land was higher than cropland and increased with increasing of mean annual temperature (MAT) and mean annual precipitation (MAP). Fallow land enhanced the network complexity of AMF community. The abundance families (Glomeraceae and Paraglomeraceae) exhibited higher values of topological features and were more often located in central ecological positions. Long-term fallow land had a significantly higher hyphal length density, GRSP, mean weight diameter (MWD), geometric mean diameter (GMD), and C concentration of GRSP (C-GRSP) than the cropland. The soil aggregate associated soil organic carbon (SOC) was 16.8, 18.6, and 13.8% higher under fallow land compared to that under cropland at HLJ, JL, and LN study sites, respectively. The structural equation model and random forest regression revealed that AMF diversity, network complexity, and their secreted GRSP mediate the effects of CC and LUC on C-GRSP and aggregate-associated SOC. This study elucidates the climate sensitivity of C within GRSP and soil aggregates which response symmetry to LUC and highlights the potential importance of AMF in C sequestration and climate change mitigation.

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.

Interspecific trait variability and local soil conditions modulate grassland model community responses to climate

Medium‐to‐high elevation grasslands provide critical services in agriculture and ecosystem stabilization, through high biodiversity and providing food for wildlife. However, these ecosystems face elevated risks of disruption due to predicted soil and climate changes. Separating the effects of soil and climate, however, is difficult in situ, with previous experiments focusing largely on monocultures instead of natural grassland communities. We experimentally exposed model grassland communities, comprised of three species grown on either local or reference soil, to varied climatic environments along an elevational gradient in the European Alps, measuring the effects on species and community traits. Although species‐specific biomass varied across soil and climate, species' proportional contributions to community‐level biomass production remained consistent. Where species experienced low survivorship, species‐level biomass production was maintained through increased productivity of surviving individuals; however, maximum species‐level biomass was obtained under high survivorship. Species responded directionally to climatic variation, spatially separating differentially by plant traits (including height, reproduction, biomass, survival, leaf dry weight, and leaf area) consistently across all climates. Local soil variation drove stochastic trait responses across all species, with high levels of interactions occurring between site and species. This soil variability obscured climate‐driven responses: we recorded no directional trait responses for soil‐corrected traits like observed for climate‐corrected traits. Our species‐based approach contributes to our understanding of grassland community stabilization and suggests that these communities show some stability under climatic variation.

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.