More than half of the solar energy absorbed by land surfaces is currently used to evaporate water¹. Climate change is expected to intensify the hydrological cycle² and to alter evapotranspiration, with implications for ecosystem services and feedback to regional and global climate. Evapotranspiration changes may already be under way, but direct observational constraints are lacking at the global scale. Until such evidence is available, changes in the water cycle on land—a key diagnostic criterion of the effects of climate change and variability—remain uncertain. Here we provide a data-driven estimate of global land evapotranspiration from 1982 to 2008, compiled using a global monitoring network³, meteorological and remote-sensing observations, and a machine-learning algorithm⁴. In addition, we have assessed evapotranspiration variations over the same time period using an ensemble of process-based land-surface models. Our results suggest that global annual evapotranspiration increased on average by 7.1 ± 1.0 millimetres per year per decade from 1982 to 1997. After that, coincident with the last major El Niño event in 1998, the global evapotranspiration increase seems to have ceased until 2008. This change was driven primarily by moisture limitation in the Southern Hemisphere, particularly Africa and Australia. In these regions, microwave satellite observations indicate that soil moisture decreased from 1998 to 2008. Hence, increasing soil-moisture limitations on evapotranspiration largely explain the recent decline of the global land-evapotranspiration trend. Whether the changing behaviour of evapotranspiration is representative of natural climate variability or reflects a more permanent reorganization of the land water cycle is a key question for earth system science.

The FLUXNET2015 dataset provides ecosystem-scale data on CO₂, water, and energy exchange between the biosphere and the atmosphere, and other meteorological and biological measurements, from 212 sites around the globe (over 1500 site-years, up to and including year 2014). These sites, independently managed and operated, voluntarily contributed their data to create global datasets. Data were quality controlled and processed using uniform methods, to improve consistency and intercomparability across sites. The dataset is already being used in a number of applications, including ecophysiology studies, remote sensing studies, and development of ecosystem and Earth system models. FLUXNET2015 includes derived-data products, such as gap-filled time series, ecosystem respiration and photosynthetic uptake estimates, estimation of uncertainties, and metadata about the measurements, presented for the first time in this paper. In addition, 206 of these sites are for the first time distributed under a Creative Commons (CC-BY 4.0) license. This paper details this enhanced dataset and the processing methods, now made available as open-source codes, making the dataset more accessible, transparent, and reproducible.

[1] We upscaled FLUXNET observations of carbon dioxide, water, and energy fluxes to the global scale using the machine learning technique, model tree ensembles (MTE). We trained MTE to predict site‐level gross primary productivity (GPP), terrestrial ecosystem respiration (TER), net ecosystem exchange (NEE), latent energy (LE), and sensible heat (H) based on remote sensing indices, climate and meteorological data, and information on land use. We applied the trained MTEs to generate global flux fields at a 0.5° × 0.5° spatial resolution and a monthly temporal resolution from 1982 to 2008. Cross‐validation analyses revealed good performance of MTE in predicting among‐site flux variability with modeling efficiencies (MEf) between 0.64 and 0.84, except for NEE (MEf = 0.32). Performance was also good for predicting seasonal patterns (MEf between 0.84 and 0.89, except for NEE (0.64)). By comparison, predictions of monthly anomalies were not as strong (MEf between 0.29 and 0.52). Improved accounting of disturbance and lagged environmental effects, along with improved characterization of errors in the training data set, would contribute most to further reducing uncertainties. Our global estimates of LE (158 ± 7 J × 10¹⁸ yr⁻¹), H (164 ± 15 J × 10¹⁸ yr⁻¹), and GPP (119 ± 6 Pg C yr⁻¹) were similar to independent estimates. Our global TER estimate (96 ± 6 Pg C yr⁻¹) was likely underestimated by 5–10%. Hot spot regions of interannual variability in carbon fluxes occurred in semiarid to semihumid regions and were controlled by moisture supply. Overall, GPP was more important to interannual variability in NEE than TER. Our empirically derived fluxes may be used for calibration and evaluation of land surface process models and for exploratory and diagnostic assessments of the biosphere.

Terrestrial ecosystems sequester 2.1 Pg of atmospheric carbon annually. A large amount of the terrestrial sink is realized by forests. However, considerable uncertainties remain regarding the fate of this carbon over both short and long timescales. Relevant data to address these uncertainties are being collected at many sites around the world, but syntheses of these data are still sparse. To facilitate future synthesis activities, we have assembled a comprehensive global database for forest ecosystems, which includes carbon budget variables (fluxes and stocks), ecosystem traits (e.g. leaf area index, age), as well as ancillary site information such as management regime, climate, and soil characteristics. This publicly available database can be used to quantify global, regional or biome‐specific carbon budgets; to re‐examine established relationships; to test emerging hypotheses about ecosystem functioning [e.g. a constant net ecosystem production (NEP) to gross primary production (GPP) ratio]; and as benchmarks for model evaluations. In this paper, we present the first analysis of this database. We discuss the climatic influences on GPP, net primary production (NPP) and NEP and present the CO₂ balances for boreal, temperate, and tropical forest biomes based on micrometeorological, ecophysiological, and biometric flux and inventory estimates. Globally, GPP of forests benefited from higher temperatures and precipitation whereas NPP saturated above either a threshold of 1500 mm precipitation or a mean annual temperature of 10 °C. The global pattern in NEP was insensitive to climate and is hypothesized to be mainly determined by nonclimatic conditions such as successional stage, management, site history, and site disturbance. In all biomes, closing the CO₂ balance required the introduction of substantial biome‐specific closure terms. Nonclosure was taken as an indication that respiratory processes, advection, and non‐CO₂ carbon fluxes are not presently being adequately accounted for.

We use eddy covariance measurements of net ecosystem productivity (NEP) from 21 FLUXNET sites (153 site-years of data) to investigate relationships between phenology and productivity (in terms of both NEP and gross ecosystem photosynthesis, GEP) in temperate and boreal forests. Results are used to evaluate the plausibility of four different conceptual models. Phenological indicators were derived from the eddy covariance time series, and from remote sensing and models. We examine spatial patterns (across sites) and temporal patterns (across years); an important conclusion is that it is likely that neither of these accurately represents how productivity will respond to future phenological shifts resulting from ongoing climate change. In spring and autumn, increased GEP resulting from an ‘extra’ day tends to be offset by concurrent, but smaller, increases in ecosystem respiration, and thus the effect on NEP is still positive. Spring productivity anomalies appear to have carry-over effects that translate to productivity anomalies in the following autumn, but it is not clear that these result directly from phenological anomalies. Finally, the productivity of evergreen needleleaf forests is less sensitive to phenology than is productivity of deciduous broadleaf forests. This has implications for how climate change may drive shifts in competition within mixed-species stands.