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.

Spatio-temporal fields of land–atmosphere fluxes derived from data-driven models can complement simulations by process-based land surface models. While a number of strategies for empirical models with eddy-covariance flux data have been applied, a systematic intercomparison of these methods has been missing so far. In this study, we performed a cross-validation experiment for predicting carbon dioxide, latent heat, sensible heat and net radiation fluxes across different ecosystem types with 11 machine learning (ML) methods from four different classes (kernel methods, neural networks, tree methods, and regression splines). We applied two complementary setups: (1) 8-day average fluxes based on remotely sensed data and (2) daily mean fluxes based on meteorological data and a mean seasonal cycle of remotely sensed variables. The patterns of predictions from different ML and experimental setups were highly consistent. There were systematic differences in performance among the fluxes, with the following ascending order: net ecosystem exchange (R² < 0.5), ecosystem respiration (R² > 0.6), gross primary production (R²> 0.7), latent heat (R² > 0.7), sensible heat (R² > 0.7), and net radiation (R² > 0.8). The ML methods predicted the across-site variability and the mean seasonal cycle of the observed fluxes very well (R² > 0.7), while the 8-day deviations from the mean seasonal cycle were not well predicted (R² < 0.5). Fluxes were better predicted at forested and temperate climate sites than at sites in extreme climates or less represented by training data (e.g., the tropics). The evaluated large ensemble of ML-based models will be the basis of new global flux products.

Significance: Carbon uptake by terrestrial ecosystems mitigates the impact of anthropogenic fossil fuel emissions on atmospheric CO ₂ concentrations, but the strength of this carbon sink is highly sensitive to large-scale extreme climate events. In 2012, the United States experienced the most severe drought since the Dust Bowl period, along with the warmest spring on record. Here, we quantify the impact of this climate anomaly on the carbon cycle. Our results show that warming-induced earlier vegetation activity increased spring carbon uptake, and thus compensated for reduced carbon uptake during the summer drought in 2012. This compensation, however, came at the cost of soil moisture depletion from increased spring evapotranspiration that likely enhanced summer heating through land-atmosphere coupling.

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.

Low soil moisture and high vapour pressure deficit (VPD) cause plant water stress and lead to a variety of drought responses, including a reduction in transpiration and photosynthesis^1,2. When soils dry below critical soil moisture thresholds, ecosystems transition from energy to water limitation as stomata close to alleviate water stress^3,4. However, the mechanisms behind these thresholds remain poorly defined at the ecosystem scale. Here, by analysing observations of critical soil moisture thresholds globally, we show the prominent role of soil texture in modulating the onset of ecosystem water limitation through the soil hydraulic conductivity curve, whose steepness increases with sand fraction. This clarifies how ecosystem sensitivity to VPD versus soil moisture is shaped by soil texture, with ecosystems in sandy soils being relatively more sensitive to soil drying, whereas ecosystems in clayey soils are relatively more sensitive to VPD. For the same reason, plants in sandy soils have limited potential to adjust to water limitations, which has an impact on how climate change affects terrestrial ecosystems. In summary, although vegetation–atmosphere exchanges are driven by atmospheric conditions and mediated by plant adjustments, their fate is ultimately dependent on the soil.