Empirical vegetation indices derived from spectral reflectance data are widely used in remote sensing of the biosphere, as they represent robust proxies for canopy structure, leaf pigment content, and, subsequently, plant photosynthetic potential. Here, we generalize the broad family of commonly used vegetation indices by exploiting all higher-order relations between the spectral channels involved. This results in a higher sensitivity to vegetation biophysical and physiological parameters. The presented nonlinear generalization of the celebrated normalized difference vegetation index (NDVI) consistently improves accuracy in monitoring key parameters, such as leaf area index, gross primary productivity, and sun-induced chlorophyll fluorescence. Results suggest that the statistical approach maximally exploits the spectral information and addresses long-standing problems in satellite Earth Observation of the terrestrial biosphere. The nonlinear NDVI will allow more accurate measures of terrestrial carbon source/sink dynamics and potentials for stabilizing atmospheric CO₂ and mitigating global climate change.

The use of deep learning (DL) approaches for the analysis of remote sensing (RS) data is rapidly increasing. DL techniques have provided excellent results in applications ranging from parameter estimation to image classification and anomaly detection. Although the vast majority of studies report precision indicators, there is a lack of studies dealing with the interpretability of the predictions. This shortcoming hampers a wider adoption of DL approaches by a wider users community, as model’s decisions are not accountable. In applications that involve the management of public budgets or policy compliance, a better interpretability of predictions is strictly required. This work aims to deepen the understanding of a recurrent neural network for land use classification based on Sentinel-2 time series in the context of the European Common Agricultural Policy (CAP). This permits to address the relevance of predictors in the classification process leading to an improved understanding of the behaviour of the network. The conducted analysis demonstrates that the red and near infrared Sentinel-2 bands convey the most useful information. With respect to the temporal information, the features derived from summer acquisitions were the most influential. These results contribute to the understanding of models used for decision making in the CAP to accomplish the European Green Deal (EGD) designed in order to counteract climate change, to protect biodiversity and ecosystems, and to ensure a fair economic return for farmers.

This paper proposes a processing chain for the derivation of global Leaf Area Index (LAI), Fraction of Absorbed Photosynthetically Active Radiation (FAPAR), Fraction Vegetation Cover (FVC), and Canopy water content (CWC) maps from 15-years of MODIS data exploiting the capabilities of the Google Earth Engine (GEE) cloud platform. The retrieval chain is based on a hybrid method inverting the PROSAIL radiative transfer model (RTM) with Random forests (RF) regression. A major feature of this work is the implementation of a retrieval chain exploiting the GEE capabilities using global and climate data records (CDR) of both MODIS surface reflectance and LAI/FAPAR datasets allowing the global estimation of biophysical variables at unprecedented timeliness. We combine a massive global compilation of leaf trait measurements (TRY), which is the baseline for more realistic leaf parametrization for the considered RTM, with large amounts of remote sensing data ingested by GEE. Moreover, the proposed retrieval chain includes the estimation of both FVC and CWC, which are not operationally produced for the MODIS sensor. The derived global estimates are validated over the BELMANIP2.1 sites network by means of an inter-comparison with the MODIS LAI/FAPAR product available in GEE. Overall, the retrieval chain exhibits great consistency with the reference MODIS product (R${}^2$ = 0.87, RMSE = 0.54 m${}^2$/m${}^2$ and ME = 0.03 m${}^2$/m${}^2$ in the case of LAI, and R${}^2$ = 0.92, RMSE = 0.09 and ME = 0.05 in the case of FAPAR). The analysis of the results by land cover type shows the lowest correlations between our retrievals and the MODIS reference estimates (R${}^2$ = 0.42 and R${}^2$ = 0.41 for LAI and FAPAR, respectively) for evergreen broadleaf forests. These discrepancies could be attributed mainly to different product definitions according to the literature. The provided results proof that GEE is a suitable high performance processing tool for global biophysical variable retrieval for a wide range of applications.