The leaf economics spectrum^1,2 and the global spectrum of plant forms and functions³ revealed fundamental axes of variation in plant traits, which represent different ecological strategies that are shaped by the evolutionary development of plant species². Ecosystem functions depend on environmental conditions and the traits of species that comprise the ecological communities⁴. 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 variability^4,5. Here we derive a set of ecosystem functions⁶ 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 ecosystems^7,8.

Evaporation (E) and transpiration (T) respond differently to ongoing changes in climate, atmospheric composition, and land use. It is difficult to partition ecosystem-scale evapotranspiration (ET) measurements into E and T, which makes it difficult to validate satellite data and land surface models. Here, we review current progress in partitioning E and T and provide a prospectus for how to improve theory and observations going forward. Recent advancements in analytical techniques create new opportunities for partitioning E and T at the ecosystem scale, but their assumptions have yet to be fully tested. For example, many approaches to partition E and T rely on the notion that plant canopy conductance and ecosystem water use efficiency exhibit optimal responses to atmospheric vapor pressure deficit (D). We use observations from 240 eddy covariance flux towers to demonstrate that optimal ecosystem response to D is a reasonable assumption, in agreement with recent studies, but more analysis is necessary to determine the conditions for which this assumption holds. Another critical assumption for many partitioning approaches is that ET can be approximated as T during ideal transpiring conditions, which has been challenged by observational studies. We demonstrate that T can exceed 95 % of ET from certain ecosystems, but other ecosystems do not appear to reach this value, which suggests that this assumption is ecosystem-dependent with implications for partitioning. It is important to further improve approaches for partitioning E and T, yet few multi-method comparisons have been undertaken to date. Advances in our understanding of carbon–water coupling at the stomatal, leaf, and canopy level open new perspectives on how to quantify T via its strong coupling with photosynthesis. Photosynthesis can be constrained at the ecosystem and global scales with emerging data sources including solar-induced fluorescence, carbonyl sulfide flux measurements, thermography, and more. Such comparisons would improve our mechanistic understanding of ecosystem water fluxes and provide the observations necessary to validate remote sensing algorithms and land surface models to understand the changing global water cycle.

Mapping in situ eddy covariance measurements of terrestrial land–atmosphere fluxes to the globe is a key method for diagnosing the Earth system from a data-driven perspective. We describe the first global products (called X-BASE) from a newly implemented upscaling framework, FLUXCOM-X, representing an advancement from the previous generation of FLUXCOM products in terms of flexibility and technical capabilities. The X-BASE products are comprised of estimates of CO₂ net ecosystem exchange (NEE), gross primary productivity (GPP), evapotranspiration (ET), and for the first time a novel, fully data-driven global transpiration product (ET_T), at high spatial (0.05°) and temporal (hourly) resolution. X-BASE estimates the global NEE at −5.75 ± 0.33 Pg C yr⁻¹ for the period 2001–2020, showing a much higher consistency with independent atmospheric carbon cycle constraints compared to the previous versions of FLUXCOM. The improvement of global NEE was likely only possible thanks to the international effort to increase the precision and consistency of eddy covariance collection and processing pipelines, as well as to the extension of the measurements to more site years resulting in a wider coverage of bioclimatic conditions. However, X-BASE global net ecosystem exchange shows a very low interannual variability, which is common to state-of-the-art data-driven flux products and remains a scientific challenge. With 125 ± 2.1 Pg C yr⁻¹ for the same period, X-BASE GPP is slightly higher than previous FLUXCOM estimates, mostly in temperate and boreal areas. X-BASE evapotranspiration amounts to 74.7×10³ ± 0.9×10³ km³ globally for the years 2001–2020 but exceeds precipitation in many dry areas, likely indicating overestimation in these regions. On average 57 % of evapotranspiration is estimated to be transpiration, in good agreement with isotope-based approaches, but higher than estimates from many land surface models. Despite considerable improvements to the previous upscaling products, many further opportunities for development exist. Pathways of exploration include methodological choices in the selection and processing of eddy covariance and satellite observations, their ingestion into the framework, and the configuration of machine learning methods. For this, the new FLUXCOM-X framework was specifically designed to have the necessary flexibility to experiment, diagnose, and converge to more accurate global flux estimates.

Summary: Sun‐induced fluorescence (SIF) in the far‐red region provides a new noninvasive measurement approach that has the potential to quantify dynamic changes in light‐use efficiency and gross primary production (GPP). However, the mechanistic link between GPP and SIF is not completely understood. We analyzed the structural and functional factors controlling the emission of SIF at 760 nm (F₇₆₀) in a Mediterranean grassland manipulated with nutrient addition of nitrogen (N), phosphorous (P) or nitrogen–phosphorous (NP). Using the soil–canopy observation of photosynthesis and energy (SCOPE) model, we investigated how nutrient‐induced changes in canopy structure (i.e. changes in plant forms abundance that influence leaf inclination distribution function, LIDF) and functional traits (e.g. N content in dry mass of leaves, N%, Chlorophyll a+b concentration (Cab) and maximum carboxylation capacity (V_cmax)) affected the observed linear relationship between F₇₆₀ and GPP. We conclude that the addition of nutrients imposed a change in the abundance of different plant forms and biochemistry of the canopy that controls F₇₆₀. Changes in canopy structure mainly control the GPP–F₇₆₀ relationship, with a secondary effect of Cab and V_cmax. In order to exploit F₇₆₀ data to model GPP at the global/regional scale, canopy structural variability, biodiversity and functional traits are important factors that have to be considered.

Plant transpiration links physiological responses of vegetation to water supply and demand with hydrological, energy, and carbon budgets at the land–atmosphere interface. However, despite being the main land evaporative flux at the global scale, transpiration and its response to environmental drivers are currently not well constrained by observations. Here we introduce the first global compilation of whole-plant transpiration data from sap flow measurements (SAPFLUXNET, https://sapfluxnet.creaf.cat/, last access: 8 June 2021). We harmonized and quality-controlled individual datasets supplied by contributors worldwide in a semi-automatic data workflow implemented in the R programming language. Datasets include sub-daily time series of sap flow and hydrometeorological drivers for one or more growing seasons, as well as metadata on the stand characteristics, plant attributes, and technical details of the measurements. SAPFLUXNET contains 202 globally distributed datasets with sap flow time series for 2714 plants, mostly trees, of 174 species. SAPFLUXNET has a broad bioclimatic coverage, with woodland/shrubland and temperate forest biomes especially well represented (80 % of the datasets). The measurements cover a wide variety of stand structural characteristics and plant sizes. The datasets encompass the period between 1995 and 2018, with 50 % of the datasets being at least 3 years long. Accompanying radiation and vapour pressure deficit data are available for most of the datasets, while on-site soil water content is available for 56 % of the datasets. Many datasets contain data for species that make up 90 % or more of the total stand basal area, allowing the estimation of stand transpiration in diverse ecological settings. SAPFLUXNET adds to existing plant trait datasets, ecosystem flux networks, and remote sensing products to help increase our understanding of plant water use, plant responses to drought, and ecohydrological processes. SAPFLUXNET version 0.1.5 is freely available from the Zenodo repository (https://doi.org/10.5281/zenodo.3971689; Poyatos et al., 2020a). The “sapfluxnetr” R package – designed to access, visualize, and process SAPFLUXNET data – is available from CRAN.

Sun-induced fluorescence in the far-red region (SIF) is increasingly used as a remote and proximal-sensing tool capable of tracking vegetation gross primary production (GPP). However, the use of SIF to probe changes in GPP is challenged during extreme climatic events, such as heatwaves. Here, we examined how the 2018 European heatwave (HW) affected the GPP-SIF relationship in evergreen broadleaved trees with a relatively invariant canopy structure. To do so, we combined canopy scale SIF measurements, GPP estimated from an eddy covariance tower, and active PAM fluorescence. The HW caused an inversion of the photosynthesis-fluorescence relationship at both the canopy and leaf scales. The highly nonlinear relationship was strongly shaped by nonphotochemical quenching (NPQ), i.e., a dissipation mechanism to protect from the adverse effects of high light intensity. During the extreme heat stress, plants experienced a saturation of NPQ, causing a change in the allocation of energy dissipation pathways towards SIF. Our results show the complex modulation of the NPQ-SIF-GPP relationship at an extreme level of heat stress, which is not completely represented in state-of-the-art coupled radiative transfer and photosynthesis models.