Laboratory for Earth Observation, Image Processing Laboratory - Scientific Park, Department of Earth Physics and Thermodynamics, Faculty of Physics, University of Valencia, C/ Dr. Moliner, 50, 46100 Burjassot, Valencia, Spain
Prof. Dr. Jose Moreno
Laboratory for Earth Observation, Image Processing Laboratory - Scientific Park, Department of Earth Physics and Thermodynamics, Faculty of Physics, University of Valencia, C/ Dr. Moliner, 50, 46100 Burjassot, Valencia, Spain
ESA’s upcoming satellite Sentinel-2 will provide Earth images of high spatial, spectral and temporal resolution and aims to ensure continuity for Landsat and SPOT observations. In comparison to the latter sensors, Sentinel-2 incorporates three new spectral bands in the red-edge region, which are centered at 705, 740 and 783 nm. This study addresses the importance of these new bands for the retrieval and monitoring of two important biophysical parameters: green leaf area index (LAI) and chlorophyll content (Ch). With data from several ESA field campaigns over agricultural sites (SPARC, AgriSAR, CEFLES2) we have evaluated the efficacy of two empirical methods that specifically make use of the new Sentinel-2 bands. First, it was shown that LAI can be derived from a generic normalized difference index (NDI) using hyperspectral data, with 674 nm with 712 nm as best performing bands. These bands are positioned closely to the Sentinel-2 B4 (665 nm) and the new red-edge B5 (705 nm) band. The method has been applied to simulated Sentinel-2 data. The resulting green LAI map was validated against field data of various crop types, thereby spanning a LAI between 0 and 6, and yielded a RMSE of 0.6. Second, the recently developed “Normalized Area Over reflectance Curve” (NAOC), an index that derives Ch from hyperspectral data, was studied on its compatibility with simulated Sentinel-2 data. This index integrates the reflectance curve between 643 and 795 nm, thereby including the new Sentinel-2 bands in the red-edge region. We found that these new bands significantly improve the accuracy of Ch estimation. Both methods emphasize the importance of red-edge bands for operational estimation of biophysical parameters from Sentinel-2.
In this paper we compare two different methodologies for Fractional Vegetation Cover (FVC) retrieval from Compact High Resolution Imaging Spectrometer (CHRIS) data onboard the European Space Agency (ESA) Project for On-Board Autonomy (PROBA) platform. The first methodology is based on empirical approaches using Vegetation Indices (VIs), in particular the Normalized Difference Vegetation Index (NDVI) and the Variable Atmospherically Resistant Index (VARI). The second methodology is based on the Spectral Mixture Analysis (SMA) technique, in which a Linear Spectral Unmixing model has been considered in order to retrieve the abundance of the different constituent materials within pixel elements, called Endmembers (EMs). These EMs were extracted from the image using three different methods: i) manual extraction using a land cover map, ii) Pixel Purity Index (PPI) and iii) Automated Morphological Endmember Extraction (AMEE). The different methodologies for FVC retrieval were applied to one PROBA/CHRIS image acquired over an agricultural area in Spain, and they were calibrated and tested against in situ measurements of FVC estimated with hemispherical photographs. The results obtained from VIs show that VARI correlates better with FVC than NDVI does, with standard errors of estimation of less than 8% in the case of VARI and less than 13% in the case of NDVI when calibrated using the in situ measurements. The results obtained from the SMA-LSU technique show Root Mean Square Errors (RMSE) below 12% when EMs are extracted from the AMEE method and around 9% when extracted from the PPI method. A RMSE value below 9% was obtained for manual extraction of EMs using a land cover use map.
Lookup-table (LUT)-based radiative transfer model inversion is considered a physically-sound and robust method to retrieve biophysical parameters from Earth observation data but regularization strategies are needed to mitigate the drawback of ill-posedness. We systematically evaluated various regularization options to improve leaf chlorophyll content (LCC) and leaf area index (LAI) retrievals over agricultural lands, including the role of (1) cost functions (CFs); (2) added noise; and (3) multiple solutions in LUT-based inversion. Three families of CFs were compared: information measures, M-estimates and minimum contrast methods. We have only selected CFs without additional parameters to be tuned, and thus they can be immediately implemented in processing chains. The coupled leaf/canopy model PROSAIL was inverted against simulated Sentinel-2 imagery at 20 m spatial resolution (8 bands) and validated against field data from the ESA-led SPARC (Barrax, Spain) campaign. For all 18 considered CFs with noise introduction and opting for the mean of multiple best solutions considerably improved retrievals; relative errors can be twice reduced as opposed to those without these regularization options. M-estimates were found most successful, but also data normalization influences the accuracy of the retrievals. Here, best LCC retrievals were obtained using a normalized “L1 -estimate” function with a relative error of 17.6% (r2 : 0.73), while best LAI retrievals were obtained through non-normalized “least-squares estimator” (LSE) with a relative error of 15.3% (r2 : 0.74).
The spatial quantification of green leaf area index (LAIgreen), the total green photosynthetically active leaf area per ground area, is a crucial biophysical variable for agroecosystem monitoring. The Sentinel-2 mission is with (1) a temporal resolution lower than a week, (2) a spatial resolution of up to 10 m, and (3) narrow bands in the red and red-edge region, a highly promising mission for agricultural monitoring. The aim of this work is to define an easy implementable LAIgreen index for the Sentinel-2 mission. Two large and independent multi-crop datasets of in situ collected LAIgreen measurements were used. Commonly used LAIgreen indices applied on the Sentinel-2 10 m × 10 m pixel resulted in a validation R2 lower than 0.6. By calculating all Sentinel-2 band combinations to identify high correlation and physical basis with LAIgreen, the new Sentinel-2 LAIgreen Index (SeLI) was defined. SeLI is a normalized index that uses the 705 nm and 865 nm centered bands, exploiting the red-edge region for low-saturating absorption sensitivity to photosynthetic vegetation. A R2 of 0.708 (root mean squared error (RMSE) = 0.67) and a R2 of 0.732 (RMSE = 0.69) were obtained with a linear fitting for the calibration and validation datasets, respectively, outperforming established indices. Sentinel-2 LAIgreen maps are presented.
Leaf Area Index (LAI) is a key variable that bridges remote sensing observations to the quantification of agroecosystem processes. In this study, we assessed the universality of the relationships between crop LAI and remotely sensed Vegetation Indices (VIs). We first compiled a global dataset of 1459 in situ quality-controlled crop LAI measurements and collected Landsat satellite images to derive five different VIs including Simple Ratio (SR), Normalized Difference Vegetation Index (NDVI), two versions of the Enhanced Vegetation Index (EVI and EVI2), and Green Chlorophyll Index (CIGreen). Based on this dataset, we developed global LAI-VI relationships for each crop type and VI using symbolic regression and Theil-Sen (TS) robust estimator. Results suggest that the global LAI-VI relationships are statistically significant, crop-specific, and mostly non-linear. These relationships explain more than half of the total variance in ground LAI observations (R2 > 0.5), and provide LAI estimates with RMSE below 1.2 m2/m2. Among the five VIs, EVI/EVI2 are the most effective, and the crop-specific LAI-EVI and LAI-EVI2 relationships constructed by TS, are robust when tested by three independent validation datasets of varied spatial scales. While the heterogeneity of agricultural landscapes leads to a diverse set of local LAI-VI relationships, the relationships provided here represent global universality on an average basis, allowing the generation of large-scale spatial-explicit LAI maps. This study contributes to the operationalization of large-area crop modeling and, by extension, has relevance to both fundamental and applied agroecosystem research.