Abstract:The ongoing monitoring of terrestrial carbon fluxes (TCF) goes hand in hand with progress in technical capacities, such as the next-generation Earth observation missions of the Copernicus initiative and advanced machine learning algorithms. Proceeding along this line, we present a physically-based data-driven workflow for quantifying gross primary productivity (GPP) and net primary productivity (NPP) at a global scale from the synergy of Copernicus' Sentinel-3 (S3) Ocean and Land Color Instrument (OLCI) and the TROPOspheric Monitoring Instrument (TROPOMI) onboard Sentinel-5 Precursor (S5P), along with meteorological variables from Copernicus ERA5-Land. Specifically, we created generic hybrid Gaussian process regression (GPR) retrieval models combining S3-OLCI-derived vegetation products with the TROPOMI solar-induced fluorescence (SIF) product to capture global GPP and NPP. First, the GPR algorithms were trained on theoretical simulations through the Soil-Canopy-Observation of Photosynthesis and Energy fluxes (SCOPE) model, with the final retrieval models termed SCOPE-GPR-TCF. Second, the SCOPE-GPR-TCF models were integrated in Google Earth Engine (GEE) and fed with satellite data and products (coming from Sentinel 3 & 5P and ERA5-Land), producing global and regional (Iberian Peninsula) maps at spatial resolutions of 5 km and 300 m during the year 2019. Moderate relative uncertainties in the range between 10%–40% of the GPP and NPP estimates were achieved by the SCOPE-GPR-TCF models. Analysis of the driving variables revealed that the S3-OLCI vegetation products, i.e., leaf area index (LAI), the fraction of absorbed photosynthetically active radiation (FAPAR), and SIF provided the highest prediction strengths. Validation of GPP temporal estimates from GPR against partitioned GPP estimates at 113 flux towers located in America and Europe highlighted a good overall consistency at the local scale, with performances varying depending on the site and vegetation type. The highest scores emerged for stations located in croplands, grasslands, deciduous broad-leaf and evergreen needle-leaf forests with top R2 and rmse values above 0.8 and below 2 μmolm−2s−1 respectively. Further, benchmarking spatiotemporal analysis revealed a strong intra-annual global correlation against reference products for the same year 2019: (i) Cross-comparison against LPJ-GUESS resulted in modal values of R = 0.8 and rmse = 1.93 μmolm−2s−1 for GPP. (ii) MOD17A2H GPP and NPP estimations cross-correlated with R modal values of 0.94 and 0.92 and rmse of 1.26 and 1.05 μmolm−2s−1 , respectively. We conclude that the hybrid models integrated into the GEE cloud-computing platform facilitate streamlining the global mapping of TCF products at efficient processing costs. This is particularly promising in preparation for the upcoming Fluorescence Explorer (FLEX) mission, where the SCOPE-GPR-TCF models are foreseen to be customized to 300 m resolution FLEX SIF data streams for high-resolution global productivity monitoring.

Estimating Crop Primary Productivity with Sentinel-2 and Landsat 8 using Machine Learning Methods Trained with Radiative Transfer Simulations

Remote Estimation of Crop Gross Primary Production with Landsat Data

Estimating Crop Yield Using a Satellite-Based Light Use Efficiency Model

Remote Estimation of Gross Primary Productivity in Crops Using MODIS 250m Data

Inferring global terrestrial carbon fluxes from the synergy of Sentinel 3 & 5P with Gaussian process hybrid models

Crop Biomass Mapping Based on Ecosystem Modeling at Regional Scale Using High Resolution Sentinel-2 Data.

Integration of Sentinel-3 OLCI Land Products and MERRA2 Meteorology Data into Light Use Efficiency and Vegetation Index-Driven Models for Modeling Gross Primary Production

Synergistic Integration of Optical and Microwave Satellite Data for Crop Yield Estimation

Novel approach for computing photosynthetically active radiation for productivity modeling using remotely sensed images in the Great Plains, United States

Exploring Simple Algorithms for Estimating Gross Primary Production in Forested Areas from Satellite Data

Predicting gross primary production from the enhanced vegetation index and photosynthetically active radiation: Evaluation and calibration

Global-scale Improvement of the Estimation of Terrestrial Gross Primary Productivity by Integrating Optical and Microwave Remote Sensing with Meteorological Data

Top-of-Atmosphere Retrieval of Multiple Crop Traits Using Variational Heteroscedastic Gaussian Processes within a Hybrid Workflow

Modeling Gross Primary Productivity for Winter Wheat–maize Double Cropping System Using MODIS Time Series and CO2 Eddy Flux Tower Data

Integrating the Sentinel-1, Sentinel-2 and Topographic data into soybean yield modelling using Machine Learning

Gaussian processes retrieval of crop traits in Google Earth Engine based on Sentinel-2 top-of-atmosphere data

Modeling Gross Primary Production by Integrating Satellite Data and Coordinated Flux Measurements in Arid and Semi-Arid China

Remote sensing of terrestrial gross primary productivity: a review of advances in theoretical foundation, key parameters and methods

A Production Efficiency Model-Based Method for Satellite Estimates of Corn and Soybean Yields in the Midwestern US

Crop productivity estimation by integrating multisensor satellite, in situ, and eddy covariance data into efficiency-based model