The leaf chlorophyll content (Chl_leaf) is a crucial vegetation parameter in carbon cycle modelling and agricultural monitoring at local, regional and global scales. The red‐edge spectral region is sensitive to variations in Chl_leaf. An increasing number of sensors are capable of sampling red‐edge bands, providing opportunities to estimate Chl_leaf. However, the contributions of canopy/foliar/soil factors are always combined in the reflectance signal, which limits the generalizability of vegetation index (VI)‐based Chl_leaf inversions. This study aims to propose a new red‐edge chlorophyll index to decouple the effects of the canopy and soil background from the Chl_leaf estimation. The chlorophyll sensitive index (CSI) was proposed, and the regression equations between the CSI and Chl_leaf were acquired using PROSAIL (PROSPECT + SAIL) and the 4‐Scale‐PROSPECT model. Sensitivity analyses showed that the CSI is resistant to variations in the canopy structure and soil background. Validation results obtained using 308 ground‐measured samples over nine sites world‐wide revealed that CSI improves the Chl_leaf retrieval accuracy (root mean square error (RMSE = 9.39 μg cm⁻²) compared with the existing Medium Resolution Imaging Spectrometer (MERIS) terrestrial chlorophyll index (MTCI; RMSE = 13.00 μg cm⁻²). Moreover, the CSI method steadily achieves a highly accurate inversion under different LAI and Chl_leaf conditions. Based on the CSI regression method, a Chl_leaf product with a 30‐m/10‐day resolution across China was generated. The CSI is sensitive to Chl_leaf but resistant to canopy structure and soil moisture parameters, and it has the potential to explicitly retrieve leaf‐scale biochemistry in ecosystem modelling and ecological applications.

Gross primary productivity (GPP) is the most important component of terrestrial carbon flux. Red-edge (680–780 nm) reflectance is sensitive to leaf chlorophyll content, which is directly correlated with photosynthesis as the pigment pool, and it has the potential to improve GPP estimation. The European Space Agency (ESA) Sentinel-2A and B satellites provide red-edge bands at 20-m spatial resolution on a five-day revisit period, which can be used for global estimation of GPP. Previous studies focused mostly on improving cropland GPP estimation using red-edge bands. In this study, we firstly evaluated the relationship between eight vegetation indices (VIs) retrieved from Sentinel-2 imagery in association with incident photosynthetic active radiation (PAR_in) and carbon flux tower GPP (GPP_EC) across three forest and two grassland sites in Australia. We derived a time series of five red-edge VIs and three non-red-edge VIs over the CO₂ flux tower footprints at 16-day time intervals and compared both temporal and spatial variations. The results showed that the relationship between the red-edge index (CIr, $\frac{\mathsf{\rho }783}{\mathsf{\rho }705}-1$ ) multiplied by PAR_in and GPP_EC had the highest correlation (R² = 0.77, root-mean-square error (RMSE) = 0.81 gC∙m⁻²∙day⁻¹) at the two grassland sites. The CIr also showed consistency (rRMSE defined as RMSE/mean GPP, lower than 0.25) across forest and grassland sites. The high spatial resolution of the Sentinel-2 data provided more detailed information to adequately characterize the GPP variance at spatially heterogeneous areas. The high revisit period of Sentinel-2 exhibited temporal variance in GPP at the grassland sites; however, at forest sites, the flux-tower-based GPP variance could not be fully tracked by the limited satellite images. These results suggest that the high-spatial-resolution red-edge index from Sentinel-2 can improve large-scale spatio-temporal GPP assessments.

Global biophysical products at decametric resolution derived from Sentinel-2 imagery have emerged as a promising dataset for fine-scale ecosystem modeling and agricultural monitoring. Evaluating uncertainties of different Sentinel-2 biophysical products over various regions and vegetation types is pivotal in the application of land surface models. In this study, we quantified the performance of Sentinel-2-derived Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), and Fractional Vegetation Cover (FVC) estimates using global ground observations with consistent measurement criteria. Our results show that the accuracy of vegetation and non-vegetated classification based on Sentinel-2 surface reflectance products is greater than 95%, which indicates the vegetation identification is favorable for the practical application of biophysical estimates, as several LAI, FAPAR, and FVC retrievals were derived for non-vegetated pixels. The rate of best retrievals is similar between LAI and FAPAR estimates, both accounting for 87% of all vegetation pixels, while it is almost 100% for FVC estimates. Additionally, the Sentinel-2 FAPAR and FVC estimates agree well with ground-measurements-derived (GMD) reference maps, whereas a large discrepancy is observed for Sentinel-2 LAI estimates by comparing with both GMD effective LAI (LAI_e) and actual LAI (LAI) reference maps. Furthermore, the uncertainties of Sentinel-2 LAI, FAPAR and FVC estimates are 1.09 m²/m², 1.14 m²/m², 0.13 and 0.17 through comparisons to ground LAI_e, LAI, FAPAR, and FVC measurements, respectively. Given the temporal difference between Sentinel-2 observations and ground measurements, Sentinel-2 LAI estimates are more consistent with LAI_e than LAI values. The robustness of evaluation results can be further improved as long as more multi-temporal ground measurements across different regions are obtained. Overall, this study provides fundamental information about the performance of Sentinel-2 LAI, FAPAR, and FVC estimates, which imbues our confidence in the broad applications of these decametric products.

Physically-based approaches for estimating Leaf Area Index (LAI) using remote sensing data rely on radiative transfer (RT) models. Currently, many RT models are freely available, but determining the appropriate RT model for LAI retrieval is still problematic. This study aims to evaluate the necessity of RT model selection for LAI retrieval and to propose a retrieval methodology using different RT models for different vegetation types. Both actual experimental observations and RT model simulations were used to conduct the evaluation. Each of them includes needleleaf forests and croplands, which have contrasting structural attributes. The scattering from arbitrarily inclined leaves (SAIL) model and the four-scale model, which are 1D and 3D RT models, respectively, were used to simulate the synthetic test datasets. The experimental test dataset was established through two field campaigns conducted in the Heihe River Basin. The results show that the realistic representation of canopy structure in RT models is very important for LAI retrieval. If an unsuitable RT model is used, then the root mean squared error (RMSE) will increase from 0.43 to 0.60 in croplands and from 0.52 to 0.63 in forests. In addition, an RT model’s potential to retrieve LAI is limited by the availability of a priori information on RT model parameters. 3D RT models require more a priori information, which makes them have poorer generalization capability than 1D models. Therefore, physically-based retrieval algorithms should embed more than one RT model to account for the availability of a priori information and variations in structural attributes among different vegetation types.

Spatial heterogeneity is present in the land surface at every scale and is one of the key factors that introduces inherent uncertainty into simulations of land surface processes and parameter retrieval based on remotely sensed data. Because of a lack of understanding of the heterogeneous characteristics of global mixed pixels, few studies have focused on modeling and inversion algorithms in heterogeneous areas. This paper presents a parameterization scheme to describe land cover heterogeneity quantitatively by composition and boundary information based on high-resolution land cover products. Global heterogeneity features at the 1-km scale are extracted from the ‘GlobeLand30’ land cover dataset with a spatial resolution of 30 m. The composition analysis of global mixed pixels shows that only 35% of pixels over the land surface of Earth are covered by a single land cover type, namely, pure pixels, and only 25.8% are located in vegetated areas. Pixels mixed with water are more common than pixels mixed with any other non-vegetation type. The fragmentation analysis of typical biomes based on the boundary length shows that the savanna is the most heterogeneous biome, while the evergreen broadleaf forest is the least heterogeneous. Deciduous needleleaf forests are significantly affected by canopy height differences, while crop and grass biomes are less affected. Lastly, the strengths and limitations of the method and the application of the land cover heterogeneity characteristics extracted in this study are discussed.