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.

The operational Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation absorbed by vegetation (FPAR) algorithm has been successfully implemented for Visible Infrared Imager Radiometer Suite (VIIRS) observations by optimizing a small set of configurable parameters in Look-Up-Tables (LUTs). Our preliminary evaluation showed reasonable agreement between VIIRS and MODIS LAI/FPAR retrievals. However, there is a need for a more comprehensive investigation to assure continuity of multi-sensor global LAI/FPAR time series, as the preliminary evaluation was spatiotemporally limited. In this study, we use a multi-year (2012–2016) global LAI/FPAR product generated from VIIRS and MODIS to evaluate for spatiotemporal consistency. We also quantify uncertainty of the product by utilizing available ground measurements. For both consistency and uncertainty evaluation, we account for variations in biome type and temporal resolution. Our results indicate that the LAI/FPAR retrievals from VIIRS and MODIS are consistent at different spatial (i.e., global and site) and temporal (i.e., 8-day, seasonal and annual) scales. The estimate of mean discrepancy (−0.006 ± 0.013 for LAI and −0.002 ± 0.002 for FPAR) meets the stability requirement for long-term LAI/FPAR Earth System Data Records (ESDRs) from multi-sensors as suggested by the Global Climate Observing System (GCOS). It is noteworthy that the rate of retrievals from the radiative transfer-based main algorithm is also comparable between two sensors. However, a relatively larger discrepancy over tropical forests was observed due to reflectance saturation and an unexpected interannual variation of main algorithm success was noticed due to instability in input surface reflectances. The uncertainties/relative uncertainties of VIIRS and MODIS LAI (FPAR) products assessed through comparisons to ground measurements are estimated to be 0.60/42.2% (0.10/24.4%) and 0.55/39.3% (0.11/26%), respectively. Note that the validated LAI were only distributed in low domains (~2.5), resulting in large relative uncertainty. Therefore, more ground measurements are needed to achieve a more comprehensive evaluation result of product uncertainty. The results presented here generally imbue confidence in the consistency between VIIRS and MODIS LAI/FPAR products and the feasibility of generating long-term multi-sensor LAI/FPAR ESDRs time series.

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.

Remote-sensing phenology detection can compensate for deficiencies in field observations and has the advantage of capturing the continuous expression of phenology on a large scale. However, there is some variability in the results of remote-sensing phenology detection derived from different vegetation parameters in satellite time-series data. Since the enhanced vegetation index (EVI) and the leaf area index (LAI) are the most widely used vegetation parameters for remote-sensing phenology extraction, this paper aims to assess the differences in phenological information extracted from EVI and LAI time series and to explore whether either index performs well for all vegetation types on a large scale. To this end, a GLASS (Global Land Surface Satellite Product)-LAI-based phenology product (GLP) was generated using the same algorithm as the MODIS (Moderate Resolution Imaging Spectroradiometer)-EVI phenology product (MLCD) over China from 2001 to 2012. The two phenology products were compared in China for different vegetation types and evaluated using ground observations. The results show that the ratio of missing data is 8.3% for the GLP, which is less than the 22.8% for the MLCD. The differences between the GLP and the MLCD become stronger as the latitude decreases, which also vary among different vegetation types. The start of the growing season (SOS) of the GLP is earlier than that of the MLCD in most vegetation types, and the end of the growing season (EOS) of the GLP is generally later than that of the MLCD. Based on ground observations, it can be suggested that the GLP performs better than the MLCD in evergreen needleleaved forests and croplands, while the MLCD performs better than the GLP in shrublands and grasslands.

The normalization of topographic effects on vegetation indices (VIs) is a prerequisite for their proper use in mountainous areas. We assessed the topographic effects on the normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the soil adjusted vegetation index (SAVI), and the near-infrared reflectance of terrestrial vegetation (NIRv) calculated from Sentinel-2. The evaluation was based on two criteria: the correlation with local illumination condition and the dependence on aspect. Results show that topographic effects can be neglected for the NDVI, while they heavily influence the SAVI, EVI, and NIRv: the local illumination condition explains 19.85%, 25.37%, and 26.69% of the variation of the SAVI, EVI, and NIRv, respectively, and the coefficients of variation across different aspects are, respectively, 8.13%, 10.46%, and 14.07%. We demonstrated the applicability of existing correction methods, including statistical-empirical (SE), sun-canopy-sensor with C-correction (SCS + C), and path length correction (PLC), dedicatedly designed for reflectance, to normalize topographic effects on VIs. Our study will benefit vegetation monitoring with VIs over mountainous areas.

This paper presents a simple radiative transfer model based on spectral invariant properties (SIP). The canopy structure parameters, including the leaf angle distribution and multi-angular clumping index, are explicitly described in the SIP model. The SIP model has been evaluated on its bidirectional reflectance factor (BRF) in the angular space at the radiation transfer model intercomparison platform, and in the spectrum space by the PROSPECT+SAIL (PROSAIL) model. The simulations of BRF by SIP agreed well with the reference values in both the angular space and spectrum space, with a root-mean-square-error (RMSE) of 0.006. When compared with the widely-used Soil-Canopy Observation of Photochemistry and Energy fluxes (SCOPE) model on fPAR, the RMSE was 0.006 and the R² was 0.99, which shows a high accuracy. This study also suggests the newly proposed vegetation index, the near-infrared (NIR) reflectance of vegetation (NIRv), was a good linear approximation of the canopy structure parameter, the directional area scattering factor (DASF), with an R² of 0.99. NIRv was not influenced much by the soil background contribution, but was sensitive to the leaf inclination angle. The sensitivity of NIRv to canopy structure and the robustness of NIRv to the soil background suggest NIRv is a promising index in future biophysical variable estimations with the support of the SIP model, especially for the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC) observations near the hot spot directions.