We present the first comparison between new fAPAR and LAI products derived from the GlobAlbedo dataset and the widely-used MODIS fAPAR and LAI products. The GlobAlbedo-derived products are produced using a 1D two-stream radiative transfer (RT) scheme designed explicitly for global parameter retrieval from albedo, with consistency between RT model assumptions and observations, as well as with typical large-scale land surface model RT schemes. The approach does not require biome-specific structural assumptions (e.g., cover, clumping, understory), unlike more detailed 3D RT model approaches. GlobAlbedo-derived values of fAPAR and LAI are compared with MODIS values over 2002–2011 at multiple flux tower sites within selected biomes, over 1200 × 1200 km regions and globally. GlobAlbedo-derived fAPAR and LAI values are temporally more stable than the MODIS values due to the smoothness of the underlying albedo, derived via optimal estimation (assimilation) using an a priori estimate of albedo derived from an albedo “climatology” (composited multi-year albedo observations). Parameters agree closely in timing but with GlobAlbedo values consistently lower than MODIS, particularly for LAI. Larger differences occur in winter (when values are lower) and in the Southern hemisphere. Globally, we find that: GlobAlbedo-derived fAPAR is ~0.9–1.01 × MODIS fAPAR with an intercept of ~0.03; GlobAlbedo-derived LAI is ~0.6 × MODIS LAI with an intercept of ~0.2. Differences arise due to the RT model assumptions underlying the products, meaning care is required in interpreting either set of values, particularly when comparing to fine-scale ground-based estimates. We present global transformations between GlobAlbedo-derived and MODIS products.

The carbon cycle of the terrestrial biosphere plays a vital role in controlling the global carbon balance and, consequently, climate change. Reliably modeled CO₂ fluxes between the terrestrial biosphere and the atmosphere are necessary in projections of policy strategies aiming at constraining carbon emissions and of future climate change. In this study, SMOS (Soil Moisture and Ocean Salinity) L3 soil moisture and JRC-TIP FAPAR (Joint Research Centre—Two-stream Inversion Package Fraction of Absorbed Photosynthetically Active Radiation) data with respective original resolutions at 10 sites were used to constrain the process-based terrestrial biosphere model, BETHY (Biosphere, Energy Transfer and Hydrology), using the carbon cycle data assimilation system (CCDAS). We find that simultaneous assimilation of these two datasets jointly at all 10 sites yields a set of model parameters that achieve the best model performance in terms of independent observations of carbon fluxes as well as soil moisture. Assimilation in a single-site mode or using only a single dataset tends to over-adjust related parameters and deteriorates the model performance of a number of processes. The optimized parameter set derived from multi-site assimilation with soil moisture and FAPAR also improves, when applied at global scale simulations, the model-data fit against atmospheric CO₂. This study demonstrates the potential of satellite-derived soil moisture and FAPAR when assimilated simultaneously in a model of the terrestrial carbon cycle to constrain terrestrial carbon fluxes. It furthermore shows that assimilation of soil moisture data helps to identity structural problems in the underlying model, i.e., missing management processes at sites covered by crops and grasslands.

Due to the substantial gross exchange fluxes with the atmosphere, the terrestrial carbon cycle plays a significant role in the global carbon budget. Drought commonly affects terrestrial carbon absorption negatively. Terrestrial biosphere models exhibit significant uncertainties in capturing the carbon flux response to drought, which have an impact on estimates of the global carbon budget. Through plant physiological processes, soil moisture tightly regulates the carbon cycle in the environment. Therefore, accurate observations of soil moisture may enhance the modeling of carbon fluxes in a model–data fusion framework. We employ the Carbon Cycle Data Assimilation System (CCDAS) to assimilate 36-year satellite-derived surface soil moisture observations in combination with flask samples of atmospheric CO₂ concentrations. We find that, compared to the default model, the performance of optimized net ecosystem productivity (NEP) and gross primary productivity (GPP) has increased with the RMSEs reduced by 1.62 gC/m²/month and 10.84 gC/m²/month, which indicates the added value of the ESA-CCI soil moisture observations as a constraint on the terrestrial carbon cycle. Additionally, the combination of soil moisture and CO₂ concentration in this study improves the representation of inter-annual variability of terrestrial carbon fluxes as well as the atmospheric CO₂ growth rate. We thereby investigate the ability of the optimized GPP in responding to drought by comparing continentally aggregated GPP with the drought index. The assimilation of surface soil moisture has been shown to efficiently capture the influences of the sub-annual (≤9 months drought durations) and large-scale (e.g., regional to continental scales) droughts on GPP. This study highlights the significant potential of satellite soil moisture for constraining inter-annual models of the terrestrial biosphere’s carbon cycle and for illustrating how GPP responds to drought at a continental scale.