Different profiles of soil phosphorous compounds depending on tree species and availability of soil phosphorus in a tropical rainforest in French Guiana.

BACKGROUND: The availability of soil phosphorus (P) often limits the productivities of wet tropical lowland forests. Little is known, however, about the metabolomic profile of different chemical P compounds with potentially different uses and about the cycling of P and their variability across space under different tree species in highly diverse tropical rainforests. RESULTS: We hypothesised that the different strategies of the competing tree species to retranslocate, mineralise, mobilise, and take up P from the soil would promote distinct soil 31P profiles. We tested this hypothesis by performing a metabolomic analysis of the soils in two rainforests in French Guiana using 31P nuclear magnetic resonance (NMR). We analysed 31P NMR chemical shifts in soil solutions of model P compounds, including inorganic phosphates, orthophosphate mono- and diesters, phosphonates, and organic polyphosphates. The identity of the tree species (growing above the soil samples) explained > 53% of the total variance of the 31P NMR metabolomic profiles of the soils, suggesting species-specific ecological niches and/or species-specific interactions with the soil microbiome and soil trophic web structure and functionality determining the use and production of P compounds. Differences at regional and topographic levels also explained some part of the the total variance of the 31P NMR profiles, although less than the influence of the tree species. Multivariate analyses of soil 31P NMR metabolomics data indicated higher soil concentrations of P biomolecules involved in the active use of P (nucleic acids and molecules involved with energy and anabolism) in soils with lower concentrations of total soil P and higher concentrations of P-storing biomolecules in soils with higher concentrations of total P. CONCLUSIONS: The results strongly suggest "niches" of soil P profiles associated with physical gradients, mostly topographic position, and with the specific distribution of species along this gradient, which is associated with species-specific strategies of soil P mineralisation, mobilisation, use, and uptake.

Soil warming increases the number of growing bacterial taxa but not their growth rates.

Soil microorganisms control the fate of soil organic carbon. Warming may accelerate their activities putting large carbon stocks at risk of decomposition. Existing knowledge about microbial responses to warming is based on community-level measurements, leaving the underlying mechanisms unexplored and hindering predictions. In a long-term soil warming experiment in a Subarctic grassland, we investigated how active populations of bacteria and archaea responded to elevated soil temperatures (+6°C) and the influence of plant roots, by measuring taxon-specific growth rates using quantitative stable isotope probing and 18O water vapor equilibration. Contrary to prior assumptions, increased community growth was associated with a greater number of active bacterial taxa rather than generally faster-growing populations. We also found that root presence enhanced bacterial growth at ambient temperatures but not at elevated temperatures, indicating a shift in plant-microbe interactions. Our results, thus, reveal a mechanism of how soil bacteria respond to warming that cannot be inferred from community-level measurements.

Divergent rhizosphere and non-rhizosphere soil microbial structure and function in long-term warmed steppe due to altered root exudation.

While there is an extensive body of research on the influence of climate warming on total soil microbial communities, our understanding of how rhizosphere and non-rhizosphere soil microorganisms respond to warming remains limited. To address this knowledge gap, we investigated the impact of 4 years of soil warming on the diversity and composition of microbial communities in the rhizosphere and non-rhizosphere soil of a temperate steppe, focusing on changes in root exudation rates and exudate compositions. We used open top chambers to simulate warming conditions, resulting in an average soil temperature increase of 1.1°C over a span of 4 years. Our results showed that, in the non-rhizosphere soil, warming had no significant impact on dissolved organic carbon concentrations, compositions, or the abundance of soil microbial functional genes related to carbon and nitrogen cycling. Moreover, soil microbial diversity and community composition remained largely unaffected, although warming resulted in increased complexity of soil bacteria and fungi in the non-rhizosphere soil. In contrast, warming resulted in a substantial decrease in root exudate carbon (by 19%) and nitrogen (by 12%) concentrations and induced changes in root exudate compositions, primarily characterized by a reduction in the abundance in alcohols, coenzymes and vitamins, and phenylpropanoids and polyketides. These changes in root exudation rates and exudate compositions resulted in significant shifts in rhizosphere soil microbial diversity and community composition, ultimately leading to a reduction in the complexity of rhizosphere bacterial and fungal community networks. Altered root exudation and rhizosphere microbial community composition therefore decreased the expression of functional genes related to soil carbon and nitrogen cycling. Interestingly, we found that changes in soil carbon-related genes were primarily driven by the fungal communities and their responses to warming, both in the rhizosphere and non-rhizosphere soil. The study of soil microbial structure and function in rhizosphere and non-rhizosphere soil provides an ideal setting for understanding mechanisms for governing rhizosphere and non-rhizosphere soil carbon and nitrogen cycles. Our results highlight the distinctly varied responses of soil microorganisms in the rhizosphere and non-rhizosphere soil to climate warming. This suggests the need for models to address these processes individually, enabling more accurate predictions of the impacts of climate change on terrestrial carbon cycling.

Design and Construction of an Experimental Setup to Enhance Mineral Weathering through the Activity of Soil Organisms.

Enhanced weathering (EW) is an emerging carbon dioxide (CO2) removal technology that can contribute to climate change mitigation. This technology relies on accelerating the natural process of mineral weathering in soils by manipulating the abiotic variables that govern this process, in particular mineral grain size and exposure to acids dissolved in water. EW mainly aims at reducing atmospheric CO2 concentrations by enhancing inorganic carbon sequestration. Until now, knowledge of EW has been mainly gained through experiments that focused on the abiotic variables known for stimulating mineral weathering, thereby neglecting the potential influence of biotic components. While bacteria, fungi, and earthworms are known to increase mineral weathering rates, the use of soil organisms in the context of EW remains underexplored. This protocol describes the design and construction of an experimental setup developed to enhance mineral weathering rates through soil organisms while concurrently controlling abiotic conditions. The setup is designed to maximize weathering rates while maintaining soil organisms' activity. It consists of a large number of columns filled with rock powder and organic material, located in a climate chamber and with water applied via a downflow irrigation system. Columns are placed above a fridge containing jerrycans to collect the leachate. Representative results demonstrate that this setup is suitable to ensure the activity of soil organisms and quantify their effect on inorganic carbon sequestration. Challenges remain in minimizing leachate losses, ensuring homogeneous ventilation through the climate chamber, and avoiding flooding of the columns. With this setup, an innovative and promising approach is proposed to enhance mineral weathering rates through the activity of soil biota and disentangle the effect of biotic and abiotic factors as drivers of EW.

Decadal soil warming decreased vascular plant above and belowground production in a subarctic grassland by inducing nitrogen limitation.

Below and aboveground vegetation dynamics are crucial in understanding how climate warming may affect terrestrial ecosystem carbon cycling. In contrast to aboveground biomass, the response of belowground biomass to long-term warming has been poorly studied. Here, we characterized the impacts of decadal geothermal warming at two levels (on average +3.3°C and +7.9°C) on below and aboveground plant biomass stocks and production in a subarctic grassland. Soil warming did not change standing root biomass and even decreased fine root production and reduced aboveground biomass and production. Decadal soil warming also did not significantly alter the root-shoot ratio. The linear stepwise regression model suggested that following 10 yr of soil warming, temperature was no longer the direct driver of these responses, but losses of soil N were. Soil N losses, due to warming-induced decreases in organic matter and water retention capacity, were identified as key driver of the decreased above and belowground production. The reduction in fine root production was accompanied by thinner roots with increased specific root area. These results indicate that after a decade of soil warming, plant productivity in the studied subarctic grassland was affected by soil warming mainly by the reduction in soil N.

Individual and interactive effects of warming and nitrogen supply on CO2 fluxes and carbon allocation in subarctic grassland.

Climate warming has been suggested to impact high latitude grasslands severely, potentially causing considerable carbon (C) losses from soil. Warming can also stimulate nitrogen (N) turnover, but it is largely unclear whether and how altered N availability impacts belowground C dynamics. Even less is known about the individual and interactive effects of warming and N availability on the fate of recently photosynthesized C in soil. On a 10-year geothermal warming gradient in Iceland, we studied the effects of soil warming and N addition on CO2 fluxes and the fate of recently photosynthesized C through CO2 flux measurements and a 13 CO2 pulse-labeling experiment. Under warming, ecosystem respiration exceeded maximum gross primary productivity, causing increased net CO2 emissions. N addition treatments revealed that, surprisingly, the plants in the warmed soil were N limited, which constrained primary productivity and decreased recently assimilated C in shoots and roots. In soil, microbes were increasingly C limited under warming and increased microbial uptake of recent C. Soil respiration was increased by warming and was fueled by increased belowground inputs and turnover of recently photosynthesized C. Our findings suggest that a decade of warming seemed to have induced a N limitation in plants and a C limitation by soil microbes. This caused a decrease in net ecosystem CO2 uptake and accelerated the respiratory release of photosynthesized C, which decreased the C sequestration potential of the grassland. Our study highlights the importance of belowground C allocation and C-N interactions in the C dynamics of subarctic ecosystems in a warmer world.

Foliar elementome and functional traits relationships identify tree species niche in French Guiana rainforests.

Biogeochemical niche (BN) hypothesis aims to relate species/genotype elemental composition with its niche based on the fact that different elements are involved differentially in distinct plant functions. We here test the BN hypothesis through the analysis of the 10 foliar elemental concentrations and 20 functional-morphological of 60 tree species in a French Guiana tropical forest. We observed strong legacy (phylogenic + species) signals in the species-specific foliar elemental composition (elementome) and, for the first time, provide empirical evidence for a relationship between species-specific foliar elementome and functional traits. Our study thus supports the BN hypothesis and confirms the general niche segregation process through which the species-specific use of bio-elements drives the high levels of α-diversity in this tropical forest. We show that the simple analysis of foliar elementomes may be used to test for BNs of co-occurring species in highly diverse ecosystems, such as tropical rainforests. Although cause and effect mechanisms of leaf functional and morphological traits in species-specific use of bio-elements require confirmation, we posit the hypothesis that divergences in functional-morphological niches and species-specific biogeochemical use are likely to have co-evolved.

Climate change increases carbon allocation to leaves in early leaf green-up.

Global greening, characterized by an increase in leaf area index (LAI), implies an increase in foliar carbon (C). Whether this increase in foliar C under climate change is due to higher photosynthesis or to higher allocation of C to leaves remains unknown. Here, we explored the trends in foliar C accumulation and allocation during leaf green-up from 2000 to 2017 using satellite-derived LAI and solar-induced chlorophyll fluorescence (SIF) across the Northern Hemisphere. The accumulation of foliar C accelerated in the early green-up period due to both increased photosynthesis and higher foliar C allocation driven by climate change. In the late stage of green-up, however, we detected decreasing trends in foliar C accumulation and foliar C allocation. Such stage-dependent trends in the accumulation and allocation of foliar C are not represented in current terrestrial biosphere models. Our results highlight that a better representation of C allocation should be incorporated into models.

Diagnosing destabilization risk in global land carbon sinks.

Global net land carbon uptake or net biome production (NBP) has increased during recent decades1. Whether its temporal variability and autocorrelation have changed during this period, however, remains elusive, even though an increase in both could indicate an increased potential for a destabilized carbon sink2,3. Here, we investigate the trends and controls of net terrestrial carbon uptake and its temporal variability and autocorrelation from 1981 to 2018 using two atmospheric-inversion models, the amplitude of the seasonal cycle of atmospheric CO2 concentration derived from nine monitoring stations distributed across the Pacific Ocean and dynamic global vegetation models. We find that annual NBP and its interdecadal variability increased globally whereas temporal autocorrelation decreased. We observe a separation of regions characterized by increasingly variable NBP, associated with warm regions and increasingly variable temperatures, lower and weaker positive trends in NBP and regions where NBP became stronger and less variable. Plant species richness presented a concave-down parabolic spatial relationship with NBP and its variability at the global scale whereas nitrogen deposition generally increased NBP. Increasing temperature and its increasing variability appear as the most important drivers of declining and increasingly variable NBP. Our results show increasing variability of NBP regionally that can be mostly attributed to climate change and that may point to destabilization of the coupled carbon-climate system.

Radiation-constrained boundaries cause nonuniform responses of the carbon uptake phenology to climatic warming in the Northern Hemisphere.

Climatic warming has lengthened the photosynthetically active season in recent decades, thus affecting the functioning and biogeochemistry of ecosystems, the global carbon cycle and climate. Temperature response of carbon uptake phenology varies spatially and temporally, even within species, and daily total intensity of radiation may play a role. We empirically modelled the thresholds of temperature and radiation under which daily carbon uptake is constrained in the temperate and cold regions of the Northern Hemisphere, which include temperate forests, boreal forests, alpine and tundra biomes. The two-dimensionality of the temperature-radiation constraint was reduced to one single variable, θ, which represents the angle in a polar coordinate system for the temperature-radiation observations during the start and end of the growing season. We found that radiation will constrain the trend towards longer growing seasons with future warming but differently during the start and end of season and depending on the biome type and region. We revealed that radiation is a major factor limiting photosynthetic activity that constrains the phenology response to temperature during the end-of-season. In contrast, the start of the carbon uptake is overall highly sensitive to temperature but not constrained by radiation at the hemispheric scale. This study thus revealed that while at the end-of-season the phenology response to warming is constrained at the hemispheric scale, at the start-of-season the advance of spring onset may continue, even if it is at a slower pace.

Delayed use of bioenergy crops might threaten climate and food security.

The potential of mitigation actions to limit global warming within 2 °C (ref. 1) might rely on the abundant supply of biomass for large-scale bioenergy with carbon capture and storage (BECCS) that is assumed to scale up markedly in the future2-5. However, the detrimental effects of climate change on crop yields may reduce the capacity of BECCS and threaten food security6-8, thus creating an unrecognized positive feedback loop on global warming. We quantified the strength of this feedback by implementing the responses of crop yields to increases in growing-season temperature, atmospheric CO2 concentration and intensity of nitrogen (N) fertilization in a compact Earth system model9. Exceeding a threshold of climate change would cause transformative changes in social-ecological systems by jeopardizing climate stability and threatening food security. If global mitigation alongside large-scale BECCS is delayed to 2060 when global warming exceeds about 2.5 °C, then the yields of agricultural residues for BECCS would be too low to meet the Paris goal of 2 °C by 2200. This risk of failure is amplified by the sustained demand for food, leading to an expansion of cropland or intensification of N fertilization to compensate for climate-induced yield losses. Our findings thereby reinforce the urgency of early mitigation, preferably by 2040, to avoid irreversible climate change and serious food crises unless other negative-emission technologies become available in the near future to compensate for the reduced capacity of BECCS.

Convergence in phosphorus constraints to photosynthesis in forests around the world.

Tropical forests take up more carbon (C) from the atmosphere per annum by photosynthesis than any other type of vegetation. Phosphorus (P) limitations to C uptake are paramount for tropical and subtropical forests around the globe. Yet the generality of photosynthesis-P relationships underlying these limitations are in question, and hence are not represented well in terrestrial biosphere models. Here we demonstrate the dependence of photosynthesis and underlying processes on both leaf N and P concentrations. The regulation of photosynthetic capacity by P was similar across four continents. Implementing P constraints in the ORCHIDEE-CNP model, gross photosynthesis was reduced by 36% across the tropics and subtropics relative to traditional N constraints and unlimiting leaf P. Our results provide a quantitative relationship for the P dependence for photosynthesis for the front-end of global terrestrial C models that is consistent with canopy leaf measurements.

Increasing calcium scarcity along Afrotropical forest succession.

Secondary forests constitute an increasingly important component of tropical forests worldwide. Although cycling of essential nutrients affects recovery trajectories of secondary forests, the effect of nutrient limitation on forest regrowth is poorly constrained. Here we use three lines of evidence from secondary forest succession sequences in central Africa to identify potential nutrient limitation in regrowing forests. First, we show that atmospheric phosphorus supply exceeds demand along forest succession, whereas forests rely on soil stocks to meet their base cation demands. Second, soil nutrient metrics indicate that available phosphorus increases along the succession, whereas available cations decrease. Finally, fine root, foliar and litter stoichiometry show that tissue calcium concentrations decline relative to those of nitrogen and phosphorus during succession. Taken together, these observations suggest that calcium becomes an increasingly scarce resource in central African forests during secondary succession. Furthermore, ecosystem calcium storage shifts from soil to woody biomass over succession, making it a vulnerable nutrient in the wake of land-use change scenarios that involve woody biomass export. Our results thus call for a broadened focus on elements other than nitrogen and phosphorus regarding tropical forest biogeochemical cycles and identify calcium as a scarce and potentially limiting nutrient in an increasingly disturbed and dynamic tropical forest landscape.

Is the climate change mitigation effect of enhanced silicate weathering governed by biological processes?

A number of negative emission technologies (NETs) have been proposed to actively remove CO2 from the atmosphere, with enhanced silicate weathering (ESW) as a relatively new NET with considerable climate change mitigation potential. Models calibrated to ESW rates in lab experiments estimate the global potential for inorganic carbon sequestration by ESW at about 0.5-5 Gt CO2  year-1 , suggesting ESW could be an important component of the future NETs mix. In real soils, however, weathering rates may differ strongly from lab conditions. Research on natural weathering has shown that biota such as plants, microbes, and macro-invertebrates can strongly affect weathering rates, but biotic effects were excluded from most ESW lab assessments. Moreover, ESW may alter soil organic carbon sequestration and greenhouse gas emissions by influencing physicochemical and biological processes, which holds the potential to perpetuate even larger negative emissions. Here, we argue that it is likely that the climate change mitigation effect of ESW will be governed by biological processes, emphasizing the need to put these processes on the agenda of this emerging research field.

Global maps and factors driving forest foliar elemental composition: the importance of evolutionary history.

Consistent information on the current elemental composition of vegetation at global scale and the variables that determine it is lacking. To fill this gap, we gathered a total of 30 912 georeferenced records on woody plants foliar concentrations of nitrogen (N), phosphorus (P) and potassium (K) from published databases, and produced global maps of foliar N, P and K concentrations for woody plants using neural networks at a resolution of 1 km2 . We used data for climate, atmospheric deposition, soil and morphoclimatic groups to train the neural networks. Foliar N, P and K do not follow clear global latitudinal patterns but are consistent with the hypothesis of soil substrate age. We additionally built generalized linear mixed models to investigate the evolutionary history effect together with the effects of environmental effects. In this comparison, evolutionary history effects explained most of the variability in all cases (mostly > 60%). These results emphasize the determinant role of evolutionary history in foliar elemental composition, which should be incorporated in upcoming dynamic global vegetation models.

Contrasting nitrogen and phosphorus fertilization effects on soil terpene exchanges in a tropical forest.

Production, emission, and absorption of biogenic volatile organic compounds (BVOCs) in ecosystem soils and associated impacts of nutrient availability are unclear; thus, predictions of effects of global change on source-sink dynamic under increased atmospheric N deposition and nutrition imbalances are limited. Here, we report the dynamics of soil BVOCs under field conditions from two undisturbed tropical rainforests from French Guiana. We analyzed effects of experimental soil applications of nitrogen (N), phosphorus (P), and N + P on soil BVOC exchanges (in particular of total terpenes, monoterpenes, and sesquiterpenes), to determine source and sink dynamics between seasons (dry and wet) and elevations (upper and lower elevations corresponding to top of the hills (30 m high) and bottom of the valley). We identified 45 soil terpenoids compounds emitted to the atmosphere, comprising 26 monoterpenes and 19 sesquiterpenes; of these, it was possible to identify 13 and 7 compounds, respectively. Under ambient conditions, soils acted as sinks of these BVOCs, with greatest soil uptake recorded for sesquiterpenes at upper elevations during the wet season (-282 µg m-2 h-1). Fertilization shifted soils from a sink to source, with greatest levels of terpene emissions recorded at upper elevations during the wet season, following the addition of N (monoterpenes: 406 µg m-2 h-1) and P (sesquiterpenes: 210 µg m-2 h-1). Total soil terpene emission rates were negatively correlated with total atmospheric terpene concentrations. These results indicate likely shifts in tropical soils from sink to source of atmospheric terpenes under projected increases in N deposition under global change, with potential impacts on regional-scale atmospheric chemistry balance and ecosystem function.

Decay of similarity across tropical forest communities: integrating spatial distance with soil nutrients.

Understanding the mechanisms that drive the change of biotic assemblages over space and time is the main quest of community ecology. Assessing the relative importance of dispersal and environmental species selection in a range of organismic sizes and motilities has been a fruitful strategy. A consensus for whether spatial and environmental distances operate similarly across spatial scales and taxa, however, has yet to emerge. We used censuses of four major groups of organisms (soil bacteria, fungi, ground insects, and trees) at two observation scales (1-m2 sampling point vs. 2,500-m2 plots) in a topographically standardized sampling design replicated in two tropical rainforests with contrasting relationships between spatial distance and nutrient availability. We modeled the decay of assemblage similarity for each taxon set and site to assess the relative contributions of spatial distance and nutrient availability distance. Then, we evaluated the potentially structuring effect of tree composition over all other taxa. The similarity of nutrient content in the litter and topsoil had a stronger and more consistent selective effect than did dispersal limitation, particularly for bacteria, fungi, and trees at the plot level. Ground insects, the only group assessed with the capacity of active dispersal, had the highest species turnover and the flattest nonsignificant distance-decay relationship, suggesting that neither dispersal limitation nor nutrient availability were fundamental drivers of their community assembly at this scale of analysis. Only the fungal communities at one of our study sites were clearly coordinated with tree composition. The spatial distance at the smallest scale was more important than nutrient selection for the bacteria, fungi, and insects. The lower initial similarity and the moderate variation in composition identified by these distance-decay models, however, suggested that the effects of stochastic sampling were important at this smaller spatial scale. Our results highlight the importance of nutrients as one of the main environmental drivers of rainforest communities irrespective of organismic or propagule size and how the overriding effect of the analytical scale influences the interpretation, leading to the perception of greater importance of dispersal limitation and ecological drift over selection associated with environmental niches at decreasing observation scales.

Latitudinal shifts of soil microbial biomass seasonality.

Soil microbes ultimately drive the mineralization of soil organic carbon and thus ecosystem functions. We compiled a dataset of the seasonality of microbial biomass carbon (MBC) and developed a semi-mechanistic model to map monthly MBC across the globe. MBC exhibits an equatorially symmetric seasonality between the Northern and Southern Hemispheres. In the Northern Hemisphere, MBC peaks in autumn and is minimal in spring at low latitudes (<25°N), peaks in the spring and is minimal in autumn at mid-latitudes (25°N to 50°N), while peaks in autumn and is minimal in spring at high latitudes (>50°N). This latitudinal shift of MBC seasonality is attributed to an interaction of soil temperature, soil moisture, and substrate availability. The MBC seasonality is inconsistent with patterns of heterotrophic respiration, indicating that MBC as a proxy for microbial activity is inappropriate at this resolution. This study highlights the need to explicitly represent microbial physiology in microbial models. The interactive controls of environments and substrate on microbial seasonality provide insights for better representing microbial mechanisms in simulating ecosystem functions at the seasonal scale.

Tree Species and Epiphyte Taxa Determine the "Metabolomic niche" of Canopy Suspended Soils in a Species-Rich Lowland Tropical Rainforest.

Tropical forests are biodiversity hotspots, but it is not well understood how this diversity is structured and maintained. One hypothesis rests on the generation of a range of metabolic niches, with varied composition, supporting a high species diversity. Characterizing soil metabolomes can reveal fine-scale differences in composition and potentially help explain variation across these habitats. In particular, little is known about canopy soils, which are unique habitats that are likely to be sources of additional biodiversity and biogeochemical cycling in tropical forests. We studied the effects of diverse tree species and epiphytes on soil metabolomic profiles of forest floor and canopy suspended soils in a French Guianese rainforest. We found that the metabolomic profiles of canopy suspended soils were distinct from those of forest floor soils, differing between epiphyte-associated and non-epiphyte suspended soils, and the metabolomic profiles of suspended soils varied with host tree species, regardless of association with epiphyte. Thus, tree species is a key driver of rainforest suspended soil metabolomics. We found greater abundance of metabolites in suspended soils, particularly in groups associated with plants, such as phenolic compounds, and with metabolic pathways related to amino acids, nucleotides, and energy metabolism, due to the greater relative proportion of tree and epiphyte organic material derived from litter and root exudates, indicating a strong legacy of parent biological material. Our study provides evidence for the role of tree and epiphyte species in canopy soil metabolomic composition and in maintaining the high levels of soil metabolome diversity in this tropical rainforest. It is likely that a wide array of canopy microsite-level environmental conditions, which reflect interactions between trees and epiphytes, increase the microscale diversity in suspended soil metabolomes.

The three major axes of terrestrial ecosystem function.

The leaf economics spectrum1,2 and the global spectrum of plant forms and functions3 revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species2. Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities4. However, the axes of variation of ecosystem functions are largely unknown, which limits our understanding of how ecosystems respond as a whole to anthropogenic drivers, climate and environmental variability4,5. Here we derive a set of ecosystem functions6 from a dataset of surface gas exchange measurements across major terrestrial biomes. We find that most of the variability within ecosystem functions (71.8%) is captured by three key axes. The first axis reflects maximum ecosystem productivity and is mostly explained by vegetation structure. The second axis reflects ecosystem water-use strategies and is jointly explained by variation in vegetation height and climate. The third axis, which represents ecosystem carbon-use efficiency, features a gradient related to aridity, and is explained primarily by variation in vegetation structure. We show that two state-of-the-art land surface models reproduce the first and most important axis of ecosystem functions. However, the models tend to simulate more strongly correlated functions than those observed, which limits their ability to accurately predict the full range of responses to environmental changes in carbon, water and energy cycling in terrestrial ecosystems7,8.