Microbial competition for phosphorus limits the CO2 response of a mature forest.

The capacity for terrestrial ecosystems to sequester additional carbon (C) with rising CO2 concentrations depends on soil nutrient availability1,2. Previous evidence suggested that mature forests growing on phosphorus (P)-deprived soils had limited capacity to sequester extra biomass under elevated CO2 (refs. 3-6), but uncertainty about ecosystem P cycling and its CO2 response represents a crucial bottleneck for mechanistic prediction of the land C sink under climate change7. Here, by compiling the first comprehensive P budget for a P-limited mature forest exposed to elevated CO2, we show a high likelihood that P captured by soil microorganisms constrains ecosystem P recycling and availability for plant uptake. Trees used P efficiently, but microbial pre-emption of mineralized soil P seemed to limit the capacity of trees for increased P uptake and assimilation under elevated CO2 and, therefore, their capacity to sequester extra C. Plant strategies to stimulate microbial P cycling and plant P uptake, such as increasing rhizosphere C release to soil, will probably be necessary for P-limited forests to increase C capture into new biomass. Our results identify the key mechanisms by which P availability limits CO2 fertilization of tree growth and will guide the development of Earth system models to predict future long-term C storage.

The fate of carbon in a mature forest under carbon dioxide enrichment.

Atmospheric carbon dioxide enrichment (eCO2) can enhance plant carbon uptake and growth1-5, thereby providing an important negative feedback to climate change by slowing the rate of increase of the atmospheric CO2 concentration6. Although evidence gathered from young aggrading forests has generally indicated a strong CO2 fertilization effect on biomass growth3-5, it is unclear whether mature forests respond to eCO2 in a similar way. In mature trees and forest stands7-10, photosynthetic uptake has been found to increase under eCO2 without any apparent accompanying growth response, leaving the fate of additional carbon fixed under eCO2 unclear4,5,7-11. Here using data from the first ecosystem-scale Free-Air CO2 Enrichment (FACE) experiment in a mature forest, we constructed a comprehensive ecosystem carbon budget to track the fate of carbon as the forest responded to four years of eCO2 exposure. We show that, although the eCO2 treatment of +150 parts per million (+38 per cent) above ambient levels induced a 12 per cent (+247 grams of carbon per square metre per year) increase in carbon uptake through gross primary production, this additional carbon uptake did not lead to increased carbon sequestration at the ecosystem level. Instead, the majority of the extra carbon was emitted back into the atmosphere via several respiratory fluxes, with increased soil respiration alone accounting for half of the total uptake surplus. Our results call into question the predominant thinking that the capacity of forests to act as carbon sinks will be generally enhanced under eCO2, and challenge the efficacy of climate mitigation strategies that rely on ubiquitous CO2 fertilization as a driver of increased carbon sinks in global forests.

Decoupled responses of above- and below-ground beta-diversity to nitrogen enrichment in a typical steppe.

Increased atmospheric nitrogen (N) deposition affects biodiversity in terrestrial ecosystems. However, we do not know whether the effects of N on above-ground plant ß-diversity are coupled with changes occurring in the soil seed bank. We conducted a long-term N-addition experiment in a typical steppe and found that above-ground ß-diversity increased and then decreased with increasing N addition, whereas below-ground ß-diversity decreased linearly. This suggests decoupled dynamics of plant communities and their soil seed bank under N enrichment. Species substitution determined above- and below-ground ß-diversity change via an increasing role of deterministic processes with N addition. These effects were mostly driven by differential responses of the above-ground vegetation and the soil seed bank ß-diversities to N-induced changes in environmental heterogeneity, increased soil inorganic N concentrations and soil acidification. Our findings highlight the importance of considering above- and below-ground processes simultaneously for effectively conserving grassland ecosystems under N enrichment.

Depth-dependent responses of soil organic carbon under nitrogen deposition.

Emerging evidence points out that the responses of soil organic carbon (SOC) to nitrogen (N) addition differ along the soil profile, highlighting the importance of synthesizing results from different soil layers. Here, using a global meta-analysis, we found that N addition significantly enhanced topsoil (0-30 cm) SOC by 3.7% (±1.4%) in forests and grasslands. In contrast, SOC in the subsoil (30-100 cm) initially increased with N addition but decreased over time. The model selection analysis revealed that experimental duration and vegetation type are among the most important predictors across a wide range of climatic, environmental, and edaphic variables. The contrasting responses of SOC to N addition indicate the importance of considering deep soil layers, particularly for long-term continuous N deposition. Finally, the lack of depth-dependent SOC responses to N addition in experimental and modeling frameworks has likely resulted in the overestimation of changes in SOC storage under enhanced N deposition.

Deciphering microbiomes dozens of meters under our feet and their edaphoclimatic and spatial drivers.

Microbes inhabiting deep soil layers are known to be different from their counterpart in topsoil yet remain under investigation in terms of their structure, function, and how their diversity is shaped. The microbiome of deep soils (>1 m) is expected to be relatively stable and highly independent from climatic conditions. Much less is known, however, on how these microbial communities vary along climate gradients. Here, we used amplicon sequencing to investigate bacteria, archaea, and fungi along fifteen 18-m depth profiles at 20-50-cm intervals across contrasting aridity conditions in semi-arid forest ecosystems of China's Loess Plateau. Our results showed that bacterial and fungal α diversity and bacterial and archaeal community similarity declined dramatically in topsoil and remained relatively stable in deep soil. Nevertheless, deep soil microbiome still showed the functional potential of N cycling, plant-derived organic matter degradation, resource exchange, and water coordination. The deep soil microbiome had closer taxa-taxa and bacteria-fungi associations and more influence of dispersal limitation than topsoil microbiome. Geographic distance was more influential in deep soil bacteria and archaea than in topsoil. We further showed that aridity was negatively correlated with deep-soil archaeal and fungal richness, archaeal community similarity, relative abundance of plant saprotroph, and bacteria-fungi associations, but increased the relative abundance of aerobic ammonia oxidation, manganese oxidation, and arbuscular mycorrhizal in the deep soils. Root depth, complexity, soil volumetric moisture, and clay play bridging roles in the indirect effects of aridity on microbes in deep soils. Our work indicates that, even microbial communities and nutrient cycling in deep soil are susceptible to changes in water availability, with consequences for understanding the sustainability of dryland ecosystems and the whole-soil in response to aridification. Moreover, we propose that neglecting soil depth may underestimate the role of soil moisture in dryland ecosystems under future climate scenarios.

Unearthing the soil-borne microbiome of land plants.

Plant-soil biodiversity interactions are fundamental for the functioning of terrestrial ecosystems. Yet, the existence of a set of globally distributed topsoil microbial and small invertebrate organisms consistently associated with land plants (i.e., their consistent soil-borne microbiome), together with the environmental preferences and functional capabilities of these organisms, remains unknown. We conducted a standardized field survey under 150 species of land plants, including 58 species of bryophytes and 92 of vascular plants, across 124 locations from all continents. We found that, despite the immense biodiversity of soil organisms, the land plants evaluated only shared a small fraction (less than 1%) of all microbial and invertebrate taxa that were present across contrasting climatic and soil conditions and vegetation types. These consistent taxa were dominated by generalist decomposers and phagotrophs and their presence was positively correlated with the abundance of functional genes linked to mineralization. Finally, we showed that crossing environmental thresholds in aridity (aridity index of 0.65, i.e., the transition from mesic to dry ecosystems), soil pH (5.5; i.e., the transition from acidic to strongly acidic soils), and carbon (less than 2%, the lower limit of fertile soils) can result in drastic disruptions in the associations between land plants and soil organisms, with potential implications for the delivery of soil ecosystem processes under ongoing global environmental change.

Soils in warmer and less developed countries have less micronutrients globally.

Soil micronutrients are capital for the delivery of ecosystem functioning and food provision worldwide. Yet, despite their importance, the global biogeography and ecological drivers of soil micronutrients remain virtually unknown, limiting our capacity to anticipate abrupt unexpected changes in soil micronutrients in the face of climate change. Here, we analyzed >1300 topsoil samples to examine the global distribution of six metallic micronutrients (Cu, Fe, Mn, Zn, Co and Ni) across all continents, climates and vegetation types. We found that warmer arid and tropical ecosystems, present in the least developed countries, sustain the lowest contents of multiple soil micronutrients. We further provide evidence that temperature increases may potentially result in abrupt and simultaneous reductions in the content of multiple soil micronutrients when a temperature threshold of 12-14°C is crossed, which may be occurring on 3% of the planet over the next century. Altogether, our findings provide fundamental understanding of the global distribution of soil micronutrients, with direct implications for the maintenance of ecosystem functioning, rangeland management and food production in the warmest and poorest regions of the planet.

Litter and soil biodiversity jointly drive ecosystem functions.

The decomposition of litter and the supply of nutrients into and from the soil are two fundamental processes through which the above- and belowground world interact. Microbial biodiversity, and especially that of decomposers, plays a key role in these processes by helping litter decomposition. Yet the relative contribution of litter diversity and soil biodiversity in supporting multiple ecosystem services remains virtually unknown. Here we conducted a mesocosm experiment where leaf litter and soil biodiversity were manipulated to investigate their influence on plant productivity, litter decomposition, soil respiration, and enzymatic activity in the littersphere. We showed that both leaf litter diversity and soil microbial diversity (richness and community composition) independently contributed to explain multiple ecosystem functions. Fungal saprobes community composition was especially important for supporting ecosystem multifunctionality (EMF), plant production, litter decomposition, and activity of soil phosphatase when compared with bacteria or other fungal functional groups and litter species richness. Moreover, leaf litter diversity and soil microbial diversity exerted previously undescribed and significantly interactive effects on EMF and multiple individual ecosystem functions, such as litter decomposition and plant production. Together, our work provides experimental evidence supporting the independent and interactive roles of litter and belowground soil biodiversity to maintain ecosystem functions and multiple services.

Mowing enhances the positive effects of nitrogen addition on ecosystem carbon fluxes and water use efficiency in a semi-arid meadow steppe.

Grasslands are now facing a continuously increasing supply of nitrogen (N) fertilizers, resulting in alterations in ecosystem functioning, including changes in carbon (C) and water cycling. Mowing, one of the most widely used grassland management techniques, has been shown to mitigate the negative impacts of increased N availability on species richness. However, knowledge of how N addition and mowing, alone and/or in combination, affect ecosystem-level C fluxes and water use efficiency (WN) is still limited. We experimentally manipulated N fertilization (0 and 10 g N m-2 yr-1) and mowing (once per year at the end of the growing season) following a randomized block design in a meadow steppe characterized by salinization and alkalinization in northeastern China. We found that, compared to the control plots, N addition, mowing, and their interaction increased net ecosystem CO2 exchange by 65.1%, 14.7%, and 133%, and WN by 40.7%, 18.5%, and 96.1%, respectively. Nitrogen enrichment also decreased soil pH, which resulted in greater aboveground biomass (AGB). Moreover, N addition indirectly increased AGB by inducing changes in species richness. Our results indicate that mowing enhances the positive effects of N addition on ecosystem C fluxes and WN. Therefore, appropriate grassland management practices are essential to improve ecosystem C sequestration, WN, and mitigate future species diversity declines due to ecosystem eutrophication.

Soil element coupling is driven by ecological context and atomic mass.

The biogeochemical cycling of multiple soil elements is fundamental for life on Earth. Here, we conducted a global field survey across 16 chronosequences from contrasting biomes with soil ages ranging from centuries to millions of years. For this, we collected and analysed 435 topsoil samples (0-10 cm) from 87 locations. We showed that high levels of topsoil element coupling, defined as the average correlation among nineteen soil elements, are maintained over geological timescales globally. Cross-biome changes in plant biodiversity, soil microbial structure, weathering, soil pH and texture, and mineral-free unprotected organic matter content largely controlled multi-element coupling. Moreover, elements with heavier atomic mass were naturally more decoupled and unpredictable in space than those with lighter mass. Only the coupling of carbon, nitrogen and phosphorus, which are essential to life on Earth, deviated from this predictable pattern, suggesting that this anomaly may be an undeniable fingerprint of life in terrestrial soils.

Biogeography of global drylands.

Despite their extent and socio-ecological importance, a comprehensive biogeographical synthesis of drylands is lacking. Here we synthesize the biogeography of key organisms (vascular and nonvascular vegetation and soil microorganisms), attributes (functional traits, spatial patterns, plant-plant and plant-soil interactions) and processes (productivity and land cover) across global drylands. These areas have a long evolutionary history, are centers of diversification for many plant lineages and include important plant diversity hotspots. This diversity captures a strikingly high portion of the variation in leaf functional diversity observed globally. Part of this functional diversity is associated with the large variation in response and effect traits in the shrubs encroaching dryland grasslands. Aridity and its interplay with the traits of interacting plant species largely shape biogeographical patterns in plant-plant and plant-soil interactions, and in plant spatial patterns. Aridity also drives the composition of biocrust communities and vegetation productivity, which shows large geographical variation. We finish our review by discussing major research gaps, which include: studying regular vegetation spatial patterns; establishing large-scale plant and biocrust field surveys assessing individual-level trait measurements; knowing whether the impacts of plant-plant and plant-soil interactions on biodiversity are predictable; and assessing how elevated CO2 modulates future aridity conditions and plant productivity.

Global impacts of fertilization and herbivore removal on soil net nitrogen mineralization are modulated by local climate and soil properties.

Soil nitrogen (N) availability is critical for grassland functioning. However, human activities have increased the supply of biologically limiting nutrients, and changed the density and identity of mammalian herbivores. These anthropogenic changes may alter net soil N mineralization (soil net Nmin ), that is, the net balance between N mineralization and immobilization, which could severely impact grassland structure and functioning. Yet, to date, little is known about how fertilization and herbivore removal individually, or jointly, affect soil net Nmin across a wide range of grasslands that vary in soil and climatic properties. Here we collected data from 22 grasslands on five continents, all part of a globally replicated experiment, to assess how fertilization and herbivore removal affected potential (laboratory-based) and realized (field-based) soil net Nmin . Herbivore removal in the absence of fertilization did not alter potential and realized soil net Nmin . However, fertilization alone and in combination with herbivore removal consistently increased potential soil net Nmin. Realized soil net Nmin , in contrast, significantly decreased in fertilized plots where herbivores were removed. Treatment effects on potential and realized soil net Nmin were contingent on site-specific soil and climatic properties. Fertilization effects on potential soil net Nmin were larger at sites with higher mean annual precipitation (MAP) and temperature of the wettest quarter (T.q.wet). Reciprocally, realized soil net Nmin declined most strongly with fertilization and herbivore removal at sites with lower MAP and higher T.q.wet. In summary, our findings show that anthropogenic nutrient enrichment, herbivore exclusion and alterations in future climatic conditions can negatively impact soil net Nmin across global grasslands under realistic field conditions. This is an important context-dependent knowledge for grassland management worldwide.

Microbial processing of plant remains is co-limited by multiple nutrients in global grasslands.

Microbial processing of aggregate-unprotected organic matter inputs is key for soil fertility, long-term ecosystem carbon and nutrient sequestration and sustainable agriculture. We investigated the effects of adding multiple nutrients (nitrogen, phosphorus and potassium plus nine essential macro- and micro-nutrients) on decomposition and biochemical transformation of standard plant materials buried in 21 grasslands from four continents. Addition of multiple nutrients weakly but consistently increased decomposition and biochemical transformation of plant remains during the peak-season, concurrent with changes in microbial exoenzymatic activity. Higher mean annual precipitation and lower mean annual temperature were the main climatic drivers of higher decomposition rates, while biochemical transformation of plant remains was negatively related to temperature of the wettest quarter. Nutrients enhanced decomposition most at cool, high rainfall sites, indicating that in a warmer and drier future fertilized grassland soils will have an even more limited potential for microbial processing of plant remains.

Drought consistently alters the composition of soil fungal and bacterial communities in grasslands from two continents.

The effects of short-term drought on soil microbial communities remain largely unexplored, particularly at large scales and under field conditions. We used seven experimental sites from two continents (North America and Australia) to evaluate the impacts of imposed extreme drought on the abundance, community composition, richness, and function of soil bacterial and fungal communities. The sites encompassed different grassland ecosystems spanning a wide range of climatic and soil properties. Drought significantly altered the community composition of soil bacteria and, to a lesser extent, fungi in grasslands from two continents. The magnitude of the fungal community change was directly proportional to the precipitation gradient. This greater fungal sensitivity to drought at more mesic sites contrasts with the generally observed pattern of greater drought sensitivity of plant communities in more arid grasslands, suggesting that plant and microbial communities may respond differently along precipitation gradients. Actinobateria, and Chloroflexi, bacterial phyla typically dominant in dry environments, increased their relative abundance in response to drought, whereas Glomeromycetes, a fungal class regarded as widely symbiotic, decreased in relative abundance. The response of Chlamydiae and Tenericutes, two phyla of mostly pathogenic species, decreased and increased along the precipitation gradient, respectively. Soil enzyme activity consistently increased under drought, a response that was attributed to drought-induced changes in microbial community structure rather than to changes in abundance and diversity. Our results provide evidence that drought has a widespread effect on the assembly of microbial communities, one of the major drivers of soil function in terrestrial ecosystems. Such responses may have important implications for the provision of key ecosystem services, including nutrient cycling, and may result in the weakening of plant-microbial interactions and a greater incidence of certain soil-borne diseases.

Nonlinear disruption of ecological interactions in response to nitrogen deposition.

Global environmental change (GEC) is affecting species interactions and causing a rapid decline in biodiversity. In this study, I present a new Ecosystem Disruption Index to quantify the impacts of simulated nitrogen (N) deposition (0, 10, 20, and 50 kg N·ha-1 ·yr-1  + 6-7 kg N·ha-1 ·yr-1 background) on abiotic and biotic ecological interactions. This comparative index is based on pairwise linear and quadratic regression matrices. These matrices, calculated at the N treatment level, were constructed using a range of abiotic and biotic ecosystem constituents: soil pH, shrub cover, and the first component of several separate principal component analyses using soil fertility data (total carbon and N) and community data (annual plants, microorganisms, biocrusts, edaphic fauna) for a total of seven ecosystem constituents. Four years of N fertilization in a semiarid shrubland completely disrupted the network of ecological interactions, with a greater proportional increase in ecosystem disruption at low N addition levels. Biotic interactions, particularly those involving microbes, shrubs, and edaphic fauna, were more prone to be lost in response to N, whereas interactions involving soil properties were more resilient. In contrast, edaphic fauna was the only group directly affected by N addition, with mites and collembolans increasing their abundance with up to 20 kg N·ha-1 ·yr-1 and then decreasing, which supports the idea of higher-trophic-level organisms being more sensitive to disturbance due to more complex links with other ecosystem constituents. Future experimental studies evaluating the impacts of N deposition, and possibly other GEC drivers, on biodiversity and biotic and abiotic interactions may be able to explain results more effectively in the context of ecological networks as a key feature of ecosystem sensitivity.

Comparison of trends in habitat and resource selection by the Spanish Festoon, Zerynthia rumina, and the whole butterfly community in a semiarid Mediterranean ecosystem.

Butterfly community and single species based approaches were taken to establish conservation priorities within a nature reserve in Central Spain. In this study, patch type (sclerophyllous, halophilous, or disturbed), potential herbaceous nectar availability, potential woody plant nectar availability, total nectar availability, and two approximations to plant diversity (herbaceous and woody plant diversity) were evaluated as variables that account for adult butterfly density. Butterfly communities in the reserve, which consist mostly of generalist species, were denser in relatively wet areas dominated by halophilous vegetation. Diversity did not significantly vary between ecologically different transects. Total nectar availability correlated with higher butterfly densities within both undisturbed and disturbed areas, which could be primarily explained by the lack of water typical of semiarid Mediterranean climates, where fresh, nectariferous vegetation is scarce. Woody plants were also found to be important sources of nectar and shelter. In the dryer sclerophyllous sites, adult butterfly density was best explained by herbaceous plant diversity, suggesting better quality of available resources. The endangered specialist Zerynthia rumina (L.) (Lepidoptera: Papilionidae) was only present at the sclerophyllous sites. Its density was very low in all sampled transects, excluding one relatively isolated transect with high larval hostplant density. In contrast to the community-based approach, density of Z. rumina adults is better explained by the density of its larval hostplant than by nectar availability, a trend previously described for other sedentary species. Management strategies for protecting insect-rich areas should consider the specific ecological requirements of endangered species.

Biogeochemical indicators of elevated nitrogen deposition in semiarid Mediterranean ecosystems.

Nitrogen (N) deposition has doubled the natural N inputs received by ecosystems through biological N fixation and is currently a global problem that is affecting the Mediterranean regions. We evaluated the existing relationships between increased atmospheric N deposition and biogeochemical indicators related to soil chemical factors and cryptogam species across semiarid central, southern, and eastern Spain. The cryptogam species studied were the biocrust-forming species Pleurochaete squarrosa (moss) and Cladonia foliacea (lichen). Sampling sites were chosen in Quercus coccifera (kermes oak) shrublands and Pinus halepensis (Aleppo pine) forests to cover a range of inorganic N deposition representative of the levels found in the Iberian Peninsula (between 4.4 and 8.1 kg N ha(-1) year(-1)). We extended the ambient N deposition gradient by including experimental plots to which N had been added for 3 years at rates of 10, 20, and 50 kg N ha(-1) year(-1). Overall, N deposition (extant plus simulated) increased soil inorganic N availability and caused soil acidification. Nitrogen deposition increased phosphomonoesterase (PME) enzyme activity and PME/nitrate reductase (NR) ratio in both species, whereas the NR activity was reduced only in the moss. Responses of PME and NR activities were attributed to an induced N to phosphorus imbalance and to N saturation, respectively. When only considering the ambient N deposition, soil organic C and N contents were positively related to N deposition, a response driven by pine forests. The PME/NR ratios of the moss were better predictors of N deposition rates than PME or NR activities alone in shrublands, whereas no correlation between N deposition and the lichen physiology was observed. We conclude that integrative physiological measurements, such as PME/NR ratios, measured on sensitive species such as P. squarrosa, can provide useful data for national-scale biomonitoring programs, whereas soil acidification and soil C and N storage could be useful as additional corroborating ecosystem indicators of chronic N pollution.

Carbon-phosphorus cycle models overestimate CO2 enrichment response in a mature Eucalyptus forest.

The importance of phosphorus (P) in regulating ecosystem responses to climate change has fostered P-cycle implementation in land surface models, but their CO2 effects predictions have not been evaluated against measurements. Here, we perform a data-driven model evaluation where simulations of eight widely used P-enabled models were confronted with observations from a long-term free-air CO2 enrichment experiment in a mature, P-limited Eucalyptus forest. We show that most models predicted the correct sign and magnitude of the CO2 effect on ecosystem carbon (C) sequestration, but they generally overestimated the effects on plant C uptake and growth. We identify leaf-to-canopy scaling of photosynthesis, plant tissue stoichiometry, plant belowground C allocation, and the subsequent consequences for plant-microbial interaction as key areas in which models of ecosystem C-P interaction can be improved. Together, this data-model intercomparison reveals data-driven insights into the performance and functionality of P-enabled models and adds to the existing evidence that the global CO2-driven carbon sink is overestimated by models.

Protist Diversity Responses to Experimental N Deposition in Biological Crusts of a Semiarid Mediterranean Ecosystem.

Biological soil crusts (BSC) are associations of different macro and microorganisms and aggregated soil particles located on the surface of soils in many different habitats. BSC harbour a diverse and complex community of ciliates and testate amoebae. These phagotrophic protists play an important role in C and N recycling in soil ecosystems but have not been frequently studied in BSC. In this context, the effects of three increasing N inputs on ciliates and testate amoebae in crusts from a semi-arid Mediterranean ecosystem were evaluated. A field experiment with artificial N-deposition was designed to mimic the effects caused by anthropogenic N depositions. The results have shown that the protist populations of these semi-arid Mediterranean environments have lower species richness than other soil environments. The increase in N produces a net loss of diversity in the populations studied and shifts in the community structure. It has also been shown that some ciliates and testate amoebae, due to their population responses to increased N concentrations, could potentially be used as bio-indicators of N contamination in these BSCs.

Ecosystem consequences of invertebrate decline.

Human activities cause substantial changes in biodiversity.1,2 Despite ongoing concern about the implications of invertebrate decline,3,4,5,6,7 few empirical studies have examined the ecosystem consequences of invertebrate biomass loss. Here, we test the responses of six ecosystem services informed by 30 above- and belowground ecosystem variables to three levels of aboveground (i.e., vegetation associated) invertebrate community biomass (100%, 36%, and 0% of ambient biomass) in experimental grassland mesocosms in a controlled Ecotron facility. In line with recent reports on invertebrate biomass loss over the last decade, our 36% biomass treatment also represented a decrease in invertebrate abundance (-70%) and richness (-44%). Moreover, we simulated the pronounced change in invertebrate biomass and turnover in community composition across the season. We found that the loss of invertebrate biomass decreases ecosystem multifunctionality, including two critical ecosystem services, aboveground pest control and belowground decomposition, while harvested plant biomass increases, likely because less energy was channeled up the food chain. Moreover, communities and ecosystem functions become decoupled with a lower biomass of invertebrates. Our study shows that invertebrate loss threatens the integrity of grasslands by decoupling ecosystem processes and decreasing ecosystem-service supply.