Abstract:Abstract. FLUXNET assembles globally-distributed eddy covariance-based estimates of carbon fluxes between the biosphere and the atmosphere. Since eddy covariance flux towers have a relatively small footprint and are distributed unevenly across the world, upscaling the observations is necessary in order to obtain global-scale estimates of biosphere-atmosphere exchange from the flux tower network. Based on cross-consistency checks with atmospheric inversions, sun-induced fluorescence (SIF) and dynamic global vegetation models (DGVM), we provide here a systematic assessment of the latest upscaling efforts for gross primary production (GPP) and net ecosystem exchange (NEE) of the FLUXCOM initiative, where different machine learning methods, forcing datasets, and sets of predictor variables were employed. Spatial patterns of mean GPP are consistent among FLUXCOM and DGVM ensembles (R2 > 0.94 at 1° spatial resolution) while the majority of DGVMs are outside the FLUXCOM range for 70 % of the land surface. Global mean GPP magnitudes for 2008–2010 from FLUXCOM members vary within 106 and 130 PgC yr−1 with the largest uncertainty in the tropics. Seasonal variations of independent SIF estimates agree better with FLUXCOM GPP (mean global pixel-wise R2 ~ 0.75) than with GPP from DGVMs (mean global pixel wise R2 ~ 0.6). Seasonal variations of FLUXCOM NEE show good consistency with atmospheric inversion-based net land carbon fluxes, particularly for temperate and boreal regions (R2 > 0.92). Interannual variability of global NEE in FLUXCOM is underestimated compared to inversions and DGVMs. The FLUXCOM version which uses also meteorological inputs shows a strong co-variation of interannual patterns with inversions (R2 = 0.88 for 2001–2010). Mean regional NEE from FLUXCOM shows larger uptake than inversion and DGVM-based estimates, particularly in the tropics with discrepancies of up to several hundred gC m2 yr−1. These discrepancies can only partly be reconciled by carbon loss pathways that are implicit in inversions but not captured by the flux tower measurements such as carbon emissions from fires and water bodies. We hypothesize that a combination of systematic biases in the underlying eddy covariance data, in particular in tall tropical forests, and a lack of site-history effects on NEE in FLUXCOM are likely responsible for the too strong tropical carbon sink estimated by FLUXCOM. Furthermore, as FLUXCOM does not account for CO2 fertilization effects carbon flux trends are not realistic. Overall, current FLUXCOM estimates of mean annual and seasonal cycles of GPP as well as seasonal NEE variations provide useful constraints of global carbon cycling, while interannual variability patterns from FLUXCOM are valuable but require cautious interpretation. Exploring the diversity of Earth Observation data and of machine learning concepts along with improved quality and quantity of flux tower measurements will facilitate further improvements of the FLUXCOM approach overall.

Estimation of Net Ecosystem Carbon Exchange at Climate Sites by Combing Remote Sensing Data and FLUXNET2015 Data with Machine Learning Algorithms

New data‐driven method for estimation of net ecosystem carbon exchange at meteorological stations effectively increases the global carbon flux data

Estimating the Net Ecosystem Exchange at Global FLUXNET Sites Using a Random Forest Model

Monitoring of carbon-water fluxes at Eurasian meteorological stations using random forest and remote sensing

Comparative Analysis of Two Machine Learning Algorithms in Predicting Site-Level Net Ecosystem Exchange in Major Biomes

Prediction and analysis of net ecosystem carbon exchange based on gradient boosting regression and random forest

Variability and uncertainty in flux-site-scale net ecosystem exchange simulations based on machine learning and remote sensing: a systematic evaluation

Assessing Spatial Representativeness of Global Flux Tower Eddy-Covariance Measurements Using Data from FLUXNET2015

Scaling carbon fluxes from eddy covariance sites to globe: Synthesis and evaluation of the FLUXCOM approach

How Well Can Matching High Spatial Resolution Landsat Data with Flux Tower Footprints Improve Estimates of Vegetation Gross Primary Production

Eddy Covariance CO2 Flux Gap Filling for Long Data Gaps: A Novel Framework Based on Machine Learning and Time Series Decomposition

Combination of Multiple Variables and Machine Learning for Regional Cropland Water and Carbon Fluxes Estimation: A Case Study in the Haihe River Basin

Spatial Representativeness and Uncertainty of Eddy Covariance Carbon Flux Measurements for Upscaling Net Ecosystem Productivity to the Grid Scale

An Online Deep Convolutional Model of Gross Primary Productivity and Net Ecosystem Exchange Estimation for Global Forests

Evaluating Diffferent Machine Learning Methods for Upscaling Evapotranspiration from Flux Towers to the Regional Scale

Integrating Meteorological and Remote Sensing Data to Simulate Cropland Nocturnal Evapotranspiration Using Machine Learning

Gap Filling Method and Estimation of Net Ecosystem CO2 Exchange in Alpine Wetland of Qinghai–Tibet Plateau

Global Terrestrial Carbon Fluxes of 1999–2019 Estimated by Upscaling Eddy Covariance Data with a Random Forest

Evaluating Energy Fluxes of the Common Land Model Based on FLUXNET Dataset

Mapping the net ecosystem exchange of CO2 of global terrestrial systems

Upscaling Northern Peatland CO2 Fluxes Using Satellite Remote Sensing Data