Next Article in Journal
Improved Bound Fit Algorithm for Fine Delay Scheduling in a Multi-Group Scan of Ultrasonic Phased Arrays
Next Article in Special Issue
Application of a Waveguide-Mode Sensor to Blood Testing for Hepatitis B Virus, Hepatitis C Virus, Human Immunodeficiency Virus and Treponema pallidum Infection
Previous Article in Journal
A Novel Method for Polar Form of Any Degree of Multivariate Polynomials with Applications in IoT
Previous Article in Special Issue
Monte Carlo Analysis of Optical Interactions in Reflectance and Transmittance Finger Photoplethysmography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI)

1
Image Processing Laboratory (IPL), University of Valencia, 46980 Valencia, Spain
2
Council for Agricultural Research and Economics—Research Centre for Cereal and Industrial Crops, S.S. 673 km 25, 200, 71122 Foggia, Italy
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(4), 904; https://doi.org/10.3390/s19040904
Submission received: 14 January 2019 / Revised: 5 February 2019 / Accepted: 18 February 2019 / Published: 21 February 2019
(This article belongs to the Special Issue Optical Bio Sensing)

Abstract

:
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.

1. Introduction

Leaf area index (LAI), or the total one-sided leaf area per unit of ground area (m2 leaf per m2 surface or dimensionless), can be distinguished in two types. On the one hand, there is the green leaf area index (LAIgreen), representing the leaves which are photosynthetically active, being the most common type of LAI [1], and, on the other hand, there is the brown leaf area index (LAIbrown), representing the leaf area normalized which is senescent and losing photosynthetic function [2]. The Sentinel-2 mission from the European Space Agency (ESA) has, with the improved optical sensor bands in the red-edge, an increased sensitivity towards LAIgreen [2], while the shortwave infrared bands are sensitive to cellulose and lignin (dry matter) absorption [2]. Such improved capabilities to obtain more accurate quantifications of LAIgreen over large areas provides an important aspect in climatic [3], ecological [4] and biogeochemical [5] cycles models, as well as for estimating crop vegetation status [6], developing soil maps [7] and estimating light-use efficiency [8]. Its determination is crucial for the understanding of biophysical processes of crop canopies, being the main morphological parameter used for determining crop growth through the correlation with crop productivity [1,9,10]. In the context of agricultural monitoring, there is a strong interest in estimate LAIgreen parameter. Near real-time LAIgreen estimates provides the tool for farmers to obtain the crop health and growth status, further improving the effective technical support in farming practices such as fertilizer application and water management. In this way, increased crop yields and reduced costs and input resources for the agricultural sector are envisaged [11,12]. Remote sensing from satellite, aerial and unmanned aerial vehicle platforms has become a popular technique in monitoring crop LAIgreen because of its ability to acquire synoptic information at different times and spatial scales [13,14,15]. For agricultural monitoring by remote sensing, the spatial resolution should be at least 20 m and, preferably, 10 m in order to make site-specific management possible [16]. A temporal resolution of less than a week would be required to follow-up acute changes in crop condition and provide timely response in management practices. These requirements are fulfilled by the ESA’s Sentinel-2 mission, providing 10 m pixel size products with a 10-day temporal resolution. Sentinel-2 is a polar-orbiting, superspectral high-resolution imaging mission with twin polar-orbiting satellites, Sentinel-2A and 2B. The mission’s main objective is providing quality information for agricultural and forestry practices and, hence, helping manage food security [17]. With Sentinel-2A in orbit (launched 23 June 2015), the temporal resolution was not yet sufficient for real applications at the individual farmer’s level. But with the additional availability of data from Sentinel-2B (launched 7 March 2017) the revisit period goes down to five days under cloud-free conditions.
LAIgreen is functionally linked to the canopy spectral reflectance, so its retrieval from optical remote-sensing data has prompted many studies using various techniques [18,19]. Essentially, these retrieval techniques can be classified into two groups, i.e., (1) empirical retrieval methods, which typically consist of relating the biophysical parameter of interest against spectral data through linear (e.g., vegetation indices) or nonlinear (e.g., machine learning approaches) regression techniques [20,21,22,23] and (2) physically-based retrieval methods, which refers to inversion of radiative transfer models (RTMs) against remote sensing observations [24,25,26]. Concerning physical models, experimental studies using RTMs have shown great flexibility in retrieving plant cover variables, because of being able to parameterize these models to a wide range of land cover situations and sensor configurations [27,28]. However, two main drawbacks limit the use of the inversion of RTMs for operational applications. First, RTM approaches typically require some ancillary information to enable the parameterization of the physical model, which may not always be available [13,29]. An additional problem hereby is that if uncertainties are introduced the likelihood increases that the model inversion will not lead to a unique solution and extra steps are required to overcome the ill-posed problem [30]. Second, regardless of the availability of auxiliary data, there is the intrinsic risk of oversimplifying the architecture of canopy for those RTMs fast enough for operational applications. The difficulty in describing canopy structure increases in heterogeneous scenes, such as mosaics of crops at different phenological stages or complex mixtures of woodlands and/or grasslands [2,31,32]. Non-linear regression techniques are standardly used for operational LAI products. For Sentinel-2 an operational LAI product, associated with a quality indicator, is provided through the SNAP (Sentinel Application Platform) toolbox and produced through a neural network which has been trained by simulated spectra generated from well-known RTMs [33]. The algorithm is trained with simulated LAIgreen values generated from the SAIL radiative transfer model [34], which describes the canopy as a homogenous and horizontal turbid-medium, and the PROSPECT radiative transfer model [35], which considers the leaf as a succession of absorption layers. However, the accuracy of this product is shown improvable [36]. Other machine learning algorithms than neural networks have been proposed to study the retrieval opportunities of LAI from Sentinel-2 and -3 [37], solving the black box problem. However, although machine learning approaches can be fast and can capture the non-linear relationship between different parameters, they are time variant and location dependent [38].
Alternatively, linear empirical models, i.e., vegetation indices (VIs), are one of the most straightforward implementable method in an operational data processing chain. These indices relate a few spectral bands with the biophysical parameter of interest [39] in a way that enhances the spectral characteristics of a given vegetation property while minimizing the soil, atmospheric, and sun-target-sensor geometry effects [22]. Despite the positive aspects of VIs developed for LAI retrieval, their major weakness is the lack of a generally applicable index for multiple vegetation types. The best way to find efficient and robust indices is to use large and diverse field datasets, with a large variety of canopy structures [22,40]. Early studies identified the red and near-infrared (NIR) regions as sensitive to LAIgreen, resulting in the common use of the reflectance broad-bands in these regions through simple ratios [41] or normalized difference ratios [42,43]. It should be mentioned that while these indices were found to be sensitive to low LAIgreen values, they usually lose sensitivity as LAIgreen increases (typically above 2–3 according to Haboudane et al. [19]). This saturation of the reflectance at moderate to high LAIgreen values in the red range (600–700 nm) is due to the high chlorophyll absorption in this spectral range [9]. The wavelength region located in the visible–near infrared (VIS-NIR) transition, i.e., between 690 and 750 nm, generally referred as the “red-edge”, is the region between maximum chlorophyll absorption in the red, and maximum reflection (high scattering) in the NIR caused by leaf cellular structure abundance, i.e., LAI [44,45,46]. Reflectance in the red-edge transition region is much higher than in the visible range especially for these moderate to high LAIgreen values, where the upwelling radiance in the red-edge range provides a higher and less noisy signal compared to the low values in the red region. It has been specifically demonstrated, through real [47,48] and simulated spectral data [42,46], that the shape of the red-edge region and mainly the slope is strongly influenced by chlorophyll density and, hence, by LAIgreen. Despite this well-known sensitivity, practically no established indices use the red-edge region for the LAIgreen retrieval as until now no free operational satellites had narrow-bands in this region. With the Sentinel-2 satellites (13 spectral bands) not only optimal and temporal resolution for crop monitoring is guaranteed, but, moreover, also spectral configuration in the red-edge is improved, with narrow-bands centred at 705 nm (B5) and 740 nm (B6). Recent studies have explored the potential of Sentinel-2 for the LAIgreen retrieval based on simulated datasets [49,50]. But at this moment, few studies have used real Sentinel-2 images in combination with in situ datasets for agricultural applications. Moreover, these studies using the red-edge Sentinel-2 bands for LAIgreen retrieval, calibrated and validated their products for only a few crop types [51,52], leaving the robustness of a generic retrieval application still an open issue.
In this respect, we aim to develop a simple, accurate empirical algorithm for deriving LAIgreen from Sentinel-2 real data of multi-crop agricultural fields, using two large in situ field datasets. The first objective is to determine if the commonly used VIs for estimating LAIgreen may be applicable for a variety of crop types. Secondly, we want to identify the Sentinel-2 spectral bands that present the highest correlation for the estimation of a wide variation in crop LAIgreen. Based on this analysis and on a parallel study of the importance of the new Sentinel-2 red-edge bands, a new robust LAIgreen index is defined. The performance of the new index and established VIs indices are validated and applied over two distinct agricultural test sites.

2. Materials and Methods

2.1. Study Sites

Field data was collected at two study sites in the Mediterranean region (Figure 1). The first site located in Valencia (Spain), named Huerta of Valencia, is an area in an alluvial plain between the Turia river, the Mediterranean Sea and Albufera lake, with an approximate area of 12,000 hectares. The climate is typically Mediterranean with mild, wet winters and hot dry summers, and a yearly average temperature around 18 °C. Seasonal rainfall is minimal in summer and maximal in autumn and spring, with an average annual value of 454 mm [53]. A complex historical irrigation system based on irrigation ditches brings water to this fertile soil in which cereals, vineyards and olive trees were originally the main crops and nowadays accompanied by rice, tigernut and new species of vegetables and citric orchards. All these crop types are currently cultivated in small plots of size 40–100 m. For this study, crop fields were measured in the Burjassot, Moncada and Alboraya municipalities (study site central coordinates 39°31’11.73’’ N, 0°23’20.48’’ W, 18 m above sea level, Datum WGS84).
The second test site is an Italian agricultural area located near Foggia (study site central coordinates 41°27’36.76’’ N, 15°32’45.33’’ E, 70 m a.s.l., Datum WGS84). The climate is Mediterranean, but with a marked continental influence being 30 km distant from the coast. This promotes abrupt seasonal and daily temperature changes, sometimes as high as 20 °C. The average annual temperature is 15.8 °C but summers can be very hot and dry, with temperatures easily exceeding 35 °C, and winters with temperatures close to 0 °C. Rainfall is usually between 350 mm and 700 mm (average 469 mm), occurring mainly during autumn and winter [54]. The agricultural sector is the mainstay of Foggia’s economy, where grapefruit, olives, durum wheat and tomato are the majority crops for centuries. The industries present are mostly devoted to food processing, with tomato processing the major industry branch.

2.2. Green Leaf Area Index (LAIgreen) Datasets

At each study site, a large LAIgreen dataset was collected with the LAI-2200 Plant Canopy Analyzer [55], which uses a fish-eyes lens with a hemispheric field of view (±45°). The detector is composed of five concentric rings (sensitive to radiation below 490 nm). Each ring responds over a different range of zenith angles and radiation is, thus, azimuthally integrated. The measurements were collected in one sensor mode using a 180° view cap, in clear sky condition, to avoid interferences from a user’s shadow, following the Land European Remote-Sensing Instruments (VALERI) field protocol (http://w3.avignon.inra.fr/valeri/). The VALERI protocol is a sampling strategy corresponding to high spatial-resolution satellite imagery, choosing elementary sampling units (ESUs) of 20 m × 20 m for each measuring plot. Each ESU was chosen in the middle of the crop field which had minimum dimensions of 40 m × 40 m in Valencia and 100 m × 100 m in Foggia. A minimum distance of 20 m from the edges of the field was kept. To account for the spatial LAIgreen variability within each ESU, measuring points were sampled following a square spatial sampling with 5 measurements at each point (A, B, C, D and E), providing a statistically mean LAIgreen estimate per ESU (Figure 2). The centre of the ESU (sampling point A) was geo-located using a GPS providing an accuracy of less than 5 m for later matching the mean LAIgreen estimate with the corresponding Sentinel-2 reflectance data. The field protocol was identic for both study sites.
The field dataset collected in Valencia (Figure 1a), hereafter called “VLC17_ES”, consists of 79 average LAIgreen values for the respective ESUs sampled on 22 and 23 May; 18 and 19 July; and 8 and 9 November 2017. Several dates throughout the season were selected to cover a wider variety of growth stages. In total, LAIgreen data from 79 ESUs were taken containing orange tree (Citrus x sinensis), collard (Brassica oleracea), tigernut (Cyperus esculentus), potato (Solanum tuberosum), artichoke (Cynara scolymus), squash (Cucurbita pepo), alfalfa (Medicago sativa), broad bean (Vicia faba), watermelon (Citrullus lanatus), pumpkin (Cucurbita maxima), onion (Allium cepa), celery (Apium graveolens) and lettuce (Lactuca sativa). The number of ESUs, classified by LAIgreen value range, is given in the corresponding histogram of the VLC17_ES dataset (Figure 3a). Furthermore, 13 bare soil ESUs were included (LAIgreen = 0), with the aim to create a more robust and general method.
The field dataset collected in Foggia (Figure 1b), hereafter called “FOG17_IT”, was taken in the framework of the H2020 SENSAGRI (Sentinels Synergy for Agriculture, http://sensagri.eu/) project. This dataset consists of 99 average LAIgreen values collected on the 16, 21, 22 and 29 March; 5 and 13 April; 11, 17 and 30 May; 12, 15 and 21 June 2017. Mean LAIgreen data of durum wheat (Triticum durum), tomato (Solanum lycopersicum) and horse bean (Vicia faba) ESUs were measured, with the number of ESU’s specified in the histogram of Foggia test site (Figure 3b). In addition, 10 bare soil ESUs were included.
The two standard and independently collected field datasets VLC17_ES (n = 79) and FOG17_IT (n = 99) were, respectively, used as testing and validation dataset for the LAIgreen index testing and development. Both datasets are covering a wide range of crop LAI values, i.e., from 0 to 4.5, providing an optimal experimental dataset for the definition of a new general methodology.

2.3. Sentinel-2 Imagery and Sentinel Application Platform (SNAP) LAI Product

All field campaigns were carried out on days close, with a maximum of five days’ difference, to overpass dates of Sentinel-2 over the study area. Each Sentinel-2 satellite carries a Multi-Spectral Imager (MSI) instrument/sensor with a swath of 290 km on board. The MSI provides a versatile set of 13 spectral bands spanning from the visible and NIR to the shortwave infrared (SWIR) (443–2190 nm), featuring four bands at 10 m (VIS and NIR bands), six bands at 20 m (red-edge and SWIR) and three bands at 60 m spatial resolution for atmospheric correction (Table 1).
The images were downloaded directly and free of charge from the ESA server (https://scihub.copernicus.eu/). ESA provides Level-1C images, being geometrically corrected [56], with top-of-atmosphere (TOA) reflectance; and Level-2A images, being geometrically and atmospherically corrected, with top-of-canopy (TOC) reflectance. We downloaded 11 available cloud-free Level-1C acquisitions of Sentinel-2 over Valencia and Foggia study areas (Table 2). In addition, we used the SNAP toolbox to process these Level-1C images into Level-2A data, with the retrieval of the LAIgreen product accompanied by a product quality indicator [33].
For each ESU the TOA reflectance spectrum was obtained from the central pixel of the corresponding plot of the Sentinel-2 image. These images were atmospherically corrected using the Sen2Cor procedure available in the Sentinel-2 SNAP toolbox, converting TOA reflectance into TOC reflectance [57]. The Sentinel images were resampled to 10 m pixel size with all selected pixels falling entirely inside the corresponding ESU. Subsetting was done to reduce the image size and the processing time, and to cover only the study areas.

2.4. Established Vegetation Indices Analysis

The performance of commonly used LAIgreen indices was tested (Table 3, indices shaded), with the specific bands as defined by the original authors. The accuracies of each index was specifically analyzed with linear (f(x) = ax + b), polynomial of second order (f(x) = ax2 + bx + c), exponential (f(x) = a × exp(bx)) and exponential of second order fitting (f(x) = a × exp(bx) + c × exp(dx)). Additionally, we introduced several generic index formulations, i.e., with undefined specific bands, into a Matlab-based graphical user interface toolbox called ARTMO (Automated Radiative Transfer Models Operator) [58]. ARTMO consists of multiple RTMs and several retrieval toolboxes that enable the development and optimization of retrieval algorithms to convert optical images into maps of vegetation properties. The indices formulation introduced in ARTMO was based on commonly used LAIgreen indices, among other VIs typically used to estimate various biophysical variables such as chlorophyll (Table 3, non-shaded indices). The spectral indices assessment toolbox [59] was used to calibrate and validate generic indices, i.e., looking for those wavelengths that provide the best correlation with LAIgreen, using the VLC17_ES testing dataset. A cross-validation method with the k-fold technique was used to ensure more robust results [60]. This method divides the available data into k subsets. From these k sub-datasets, k-1 sub-datasets are selected as a calibration dataset and a single k sub-dataset is used for model validation. The cross-validation process is then repeated k times, with each of the k sub-datasets used as validation dataset. Thus, all VLC17_ES field data are used for both calibration and validation. Here, we used a 4-fold (k = 4) cross-validation procedure.
The selection of the index and, accordingly, the best performing bands based on the testing dataset will rely on a series of criteria. First, it should be a simple index, preferably using only two bands to minimize processing time and improve its operational use. Secondly, the index must have physical sense, that is, it should be based on areas of the spectrum with high influence of the LAIgreen parameter. As mentioned earlier, the red, red-edge and NIR regions are the areas with the greatest influence of this parameter, so the Sentinel-2 bands used by the proposed index will be from these spectral regions. Third, it must present good statistics when applied to the completely independent validation dataset (FOG17_IT). The coefficient of determination (R2) and root mean squared error (RMSE) will be selected as indicators of the accuracy of the statistical estimation models.

3. Results

In this section we first present the quality results of Sentinel-2 LAI product, comparing them with both the Valencia and Foggia LAIgreen in situ datasets. In the same way, we show the performance of the common indices used by the bibliography, presenting graphs with the best settings for each index. Next, the best band combination for each commonly used index, obtained with ARTMO toolbox, are shown, defining finally the new SeLI index and showing its validation assessment with the Foggia dataset.

3.1. Performance of the Sentinel-2 Level-2A LAIgreen Product

The LAIgreen products of both Valencia and Foggia Sentinel-2 images were obtained and compared with the in situ LAIgreen of the corresponding pixel. The analysis was carried out with only 31 LAIgreen VLC17_ES values and 99 FOG17_ES values, because only those corresponding pixels indicated with a good product quality flag were selected. Figure 4 shows the Sentinel-2 LAIgreen product values and the in situ LAIgreen values, with the 1:1 line (black), its statistics and the linear fitting (red). A clear underestimation of the product values is shown in Figure 4.

3.2. Performance of Common LAIgreen Indices for a Multi-Crop Dataset

Commonly used LAIgreen indices (Table 3, indices shaded) were evaluated with their default bands using the multi-crop VLC17_ES dataset. The R2 obtained with different types of fitting functions ranged between 0.234–0.663 when applying the respective indices on the multi-crop dataset (Table 4). Hence, the accuracies of each index obtained with linear, polynomial of second order, exponential and exponential of second order fitting, were rather low. The p-value is <0.001 in all cases except for the ratio index (RI).
Each established index was represented as a function of the VLC17_ES LAIgreen values, to show its predicting performance (Figure 5). As can be seen, they generally present a scattered performance and all indices present a saturation problem with high and/or low LAIgreen values. Concretely, the RI index (Figure 5a) overestimates to an extreme extent at low and high LAIgreen values. The TRBI (Figure 5d) shows a saturation process with LAIgreen values close to 3. And both normalized indices (Figure 5b,c) already present their greater value (the unit), with also LAIgreen values close to 3. The normalized index defined by Delegido et al. (2011) uses a red-edge band (705 nm), but as it is used in combination with the red band (665 nm) [42], the saturation problem at high LAIgreen values appears. So, the effectiveness of LAIgreen indices depends entirely on the combination of bands used.
Another demonstration of the saturation produced by the band located at the red region (B4: 665 nm) is Figure 6. In this Figure, real Sentinel-2 spectra of artichoke (LAI = 2.8 in blue, 3.4 in orange) and orange tree crops (LAI = 2.3 in blue, 3.9 in orange) with moderate-high LAIgreen values of the VLC17_ES dataset are represented. Values higher than 2 have been chosen because it is usually the limit at which saturation process starts with common LAIgreen indices (Figure 5). As can be seen, they have the same reflectance value in the 665 nm band, that is, at moderate-high LAIgreen values, chlorophyll maximally absorbs in this region, producing saturation. At the same time, the red-edge band close to the red (B5: 705 nm) does not show such entire absorption saturation, which brings the advantage of using red-edge Sentinel-2 bands for LAIgreen estimation.

3.3. Sensitivity of Spectral Bands against LAIgreen Parameter

In order to investigate more generic index options, the following step involved systematically calculating all band combinations with the ARTMO toolbox. Table 5 lists the 10 best statistical results obtained for each generic index and the corresponding bands with a linear fit, using the VLC17_ES dataset and a cross-validation process. The analysis was performed with linear fitting, because the aim of this study is to define a simple relationship between the LAIgreen parameter and Sentinel-2 bands, analysing whether there is a linear relationship between the LAIgreen in situ values and the estimated values. Furthermore, as seen in Table 4, linear fitting produces one of the best statistical results.
Comparing to the tested established indices, these results already show a more promising correlation with a R2 ranging between 0.701 and 0.737. However, questions arose when evaluating the obtained wavelengths of the resulting best-performing bands from a physical point of view. In the majority of cases, the selected bands were physically not only influenced by chlorophyll absorption, but mainly by other leaf constituents such as lignin, cellulose and water (e.g., 1610 nm, 2190 nm) affecting the scattering properties in the NIR and SWIR [69,70]. The only case where physical chlorophyll-related bands were chosen, was the NDGI, using one of the new red-edge bands (705 nm, in the tail of the chlorophyll absorption peak) and the other in the NIR region (865 nm). The red band (B4) appeared not to be chosen at all.

3.4. Optimized Simple Index for LAIgreen Retrieval from Sentinel-2 Data: SeLI

According to the previous outcome of the NDGI, we used this index structure and analysed all band combinations in red, red-edge and NIR regions, i.e., the bands 4, 5, 6, 7, 8 and 8a (Table 1). Table 6 summarizes the normalized difference ratio index band combinations with their corresponding statistics for the independent testing and validation datasets. The statistical results are ranked according their performance for the testing dataset. All indices perform better compared to the common used indices, apart from those listed at the bottom. The top four performing combinations do not use the red band, but instead, all use the red-edge band at 705 nm in combination with a far-red or NIR band.
As expected, the best result for the testing dataset is again the 705 (B5)–865 (B8a) nm combination (R2 = 0.708), and, moreover, confirmed by the validation dataset (R2 = 0.732). These both independent datasets on real in situ data give us the strong experimental proof that the LAIgreen parameter presents a linear behavior up to values of five, proposing therefore a physiologically-based LAIgreen index, hereafter called the Sentinel-2 LAIgreen Index (SeLI) (Equation (1)). For the in situ VLC17_ES LAIgreen dataset (Figure 7a), it is observed that SeLI values vary from 0.03, corresponding to bare soils, to a maximum of 0.76, corresponding to LAIgreen values up to 4.5, with potato and alfalfa crops showing the highest LAIgreen values. We do not observe any saturation at these high LAIgreen values, while previously shown indices did. Hence, SeLI showed the potential to be used in a unified algorithm for LAIgreen estimation in different crop types.
SeLI   = R 865     R 705 R 865   +   R 705
From Figure 7a, the LAIgreen estimation equation through SeLI index is extracted, with a linear fit (Equation (2)). We have tested, besides the linear fitting (R2 = 0.708, RMSE = 0.67), also an exponential fitting (R2 = 0.603, RMSE = 0.78) and a polynomial of second order (R2 = 0.727, RMSE = 0.65) to fit the index with in situ LAI values. The polynomial fitting presents slightly higher statistics, but the adjustment is negative, not presenting much physical sense.
LAIgreen = 5.405 × SeLI − 0.114
Figure 7b shows the LAIgreen values estimated with Equation (2) together to the validation dataset taken in Foggia, obtaining a correlation R2 of 0.732 and RMSE of 0.69 (1:1 line), showing strong validation statistics.
Finally, we applied the SeLI to both field sites, using the 26 May 2017 and 8 April 2017 Sentinel-2 images for the Huerta of Valencia and Foggia site, respectively. The resulting maps are shown in Figure 8, demonstrating the applicability of SeLI at high spatial resolution for two distinct agricultural areas. In brown, the pixels with the lowest value of LAIgreen are shown, mainly corresponding to bare soils, and in green colour, the pixels with the highest LAIgreen value, corresponding to potato and alfalfa crops in the case of Valencia, and wheat in the case of Foggia. A zoom of each map is shown to demonstrate the high spatial resolution of Sentinel-2, being able to observe individual agricultural plots of 40 m × 40 m (Valencia) and 100 m × 100 m (Foggia) based on LAI, while also observing slight LAI variability within these plots.

4. Discussion

With the availability of a narrow band in the red-edge region by Sentinel-2, an improved and simple estimation of LAIgreen based on a simple index becomes possible at high spatial resolution. The proposed SeLI index shows a significant improvement towards indices using the saturating bands in the red (B4 in Sentinel-2). Moreover, no saturation appeared in the obtained LAIgreen product based on the red-edge bands (B5: 705 nm and B6: 740 nm). The B4 (665 nm) saturation at high LAIgreen values is clearly shown with real Sentinel-2 TOC reflectance spectra for different LAIgreen values, higher than 2 (Figure 5 and Figure 6). The red-edge bands (B5: 705 nm and B6: 740 nm) in contrast are both affected by higher scattering, whereby the B5 band is still driven by chlorophyll absorption. This agrees with numerous authors who emphasize the importance of the red-edge bands for the estimation of biophysical parameters, mainly the LAIgreen and chlorophyll estimation [42,46,71]. The proposed Sentinel-2 LAI Index (SeLI) exploits the B5 red-edge band, which has been widely demonstrated that is highly influenced by the LAIgreen parameter [46,47,48], and the B8a NIR band, which is driven by the scattering changes in moderate-to-high LAI values in crops [72]. Very few previous indices have used bands in the red-edge region because no free operational previous sensors had narrow bands in this spectral area [42]. Both linear and non-linear empirical regression techniques have been tested for the LAI retrieval on simulated spectrally resampled airborne data [50,73] and recently on real Sentinel-2 data [10,52]. The band selection obtained from these methodologies appeared to favor (1) green and SWIR bands in the case of linear regression by VIs, and (2) red, NIR and SWIR bands in the case of non-linear regression by machine learning approaches [74]. These bands were also chosen by several of our tested VIs (Table 5), with the difference that the NDGI formulation indicated the use of red-edge (705 nm) and a NIR band (865 nm) as best band selection.
The robustness and generality of the SeLI index is demonstrated by applying it to an independent in situ field dataset from a distinct geographical location with crop types different from those included in the testing dataset, obtaining equally good statistics (R2 of 0.732, RMSE of 0.69). Specifically, SeLI does not present problems of saturation when it is applied to a multi-crop Valencia in situ dataset composed of 13 different crop types and LAIgreen values that go up to 4.5, obtaining R2 of 0.708 and RMSE of 0.67. Furthermore, when SeLI is applied to the Foggia region, characterized by high LAIgreen values, the limits of the different crop fields and LAIgreen variability within the crop field appears even for the high value ranges, indicating variable growing conditions. Such a clear distinction in LAIgreen variability allows evaluation of management practices at the field level. Hence, it is shown that the SeLI index generally can be applied for LAIgreen retrieval of different crop types and distinct areas. A limitation of the index is that it has been calibrated and validated with LAI data up to 5, so it is only applicable to agricultural areas with this range of values, although it is the common range of in situ LAI measured values in a lot of studies with a great variety of crop types; such as wheat [75,76], corn [19], potato [10] and sugar beet [2,47]. Currently, ongoing scientific debate is taking place on the discussion if there is a linear relationship or not between in situ LAI values and estimated values [77]. Our study shows that in situ data are linearly related to SeLI, in the value range of 0 to 5. This result is in accordance with other results in which also LAI of agricultural areas is estimated through indices and linear models in similar ranges [10,75,78,79,80]. Generally, no values higher than 5 are used in these studies constraining the model applicability to this range. As SeLI has a physiological foundation, the index will be applicable to a higher value range, but the SeLI-LAI fitting relationship might change depending on the dataset range. To verify this, a LAI in situ data range >5 would be required. However, one must also consider the instrumental limitations for in situ measurements. The LAI-2200 Plant Canopy Analyzer [55] calculates LAI by comparing differential light measurements above and below canopy. The maximum measurable LAI is generally lower for these devices measuring gap fraction with LAI reaching an asymptotic saturation level at a value of about 5, compared to that assessed via destructive methods. The cause for this is gap fraction saturation as LAI approaches five or six [81,82,83].
In this work, the Sentinel-2 LAIgreen product obtained from the SNAP toolbox was also tested with the multi-crop dataset from Valencia. The results show underestimated LAI estimations (R2 = 0.475, RMSE = 0.91). There are some studies, which have also compared this Sentinel-2 LAIgreen product with in situ LAIgreen crop data, that obtain better R2 [84,85], but they only tested the product with few crop types. When the product is analysed in different areas and plant species, the results can be improved [36]. This finding could be explained by the fact that the Sentinel-2 algorithm used for land surface parameters, including LAIgreen product, ingests almost all spectral bands and applies a nonlinear regression to estimate each parameter [33], in addition to the fact that it has been proven that there is a substantial sensitivity of Sentinel-2 biophysical products to the implemented rugged terrain corrections [36].
The other main challenge in the retrieval of biophysical parameters with vegetation indices is the difficulty of finding a simple index with such a general character that it can estimate the parameter of a wide variety of crop types. In this work it has also been shown that the established indices do not present this general character. This may be because they were developed and calibrated based on limited experimental data in terms of species, presenting improvable statistics (R2 between 0.234–0.663) when applied to multi-crop datasets. In an attempt to improve estimations over this multi-crop dataset, all band combinations were systematically calculated for each index in order to achieve the highest possible correlation for the estimation of LAIgreen. More promising results were obtained, with a R2 between 0.701 and 0.737. However, when inspecting these sensitive bands whether they are physically meaningful, i.e., if the selected bands are actually influenced only or mostly by LAIgreen, then these indices turned out to be questionable. In the majority of cases, the selected bands were influenced by leaf constituents such as lignin, cellulose and water (e.g., 1610 nm, 2190 nm) affecting the scattering properties in the NIR and SWIR [69,70], and being less related to photosynthetically based LAIgreen. At the field or landscape scale, canopy reflectance patterns represent the integrated effects of all biophysical parameters. Co-variation mechanisms of leaf constituents is typically causing the selection of bands related to other covarying biochemicals such as pigments or lignin due to their high effect on spectral variability [86]. Similarly, it was earlier observed that due to the covariation between water content and chlorophyll content (related with LAI parameter), typically bands in the water content absorption region are selected as most sensitive [87]. To improve the estimation of LAIgreen (aside from LAIbrown), bands only affected by structural leaf components should be omitted. Structurally-related NIR and SWIR bands may improve the LAIgreen retrieval when the model is trained on healthy vegetation [74] but may be less generally applicable for scenarios with different structural types or stress conditions. With the band selection B5 and B8a, SeLI is functionally related to green LAI, avoiding absorption saturation in the red region.
It should be mentioned that this is the first time that this kind of LAIgreen retrieval can be carried out for agricultural areas with plots sizes of only 40–100 m, such as Huerta of Valencia, due to the lack of an operational satellite with the required spatial and temporal resolution. ESA’s satellite Sentinel-2 aims to replace and improve the older generation of satellite sensors such as Landsat and SPOT, with improved spectral and spatial capabilities. Therefore, the Sentinel-2 satellites provide a great opportunity for global vegetation monitoring, and specifically crop field monitoring, due to its enhanced spatial, spectral and temporal characteristics [42].
Finally, further validation is required with other field campaigns and synthetic Sentinel-2 data to reinforce findings. Considering appropriate instrumental tools, the index behavior for LAIgreen values higher than 5 should be tested, as well as the fitting behavior of these further ranges.

5. Conclusions

Numerous VIs have been proposed for the estimation of green leaf area index (LAIgreen) over various crop types, but the general problem appears when they are applied to multi-crop datasets, obtaining low estimation accuracies and additionally at LAIgreen values higher than 2, saturation problems appear. Based on the availability of Sentinel-2 narrow red-edge bands, we explored new index possibilities for accurate LAIgreen estimations for heterogeneous agricultural areas, based on spectral areas influenced mainly by photosynthesis-related absorption regions. The proposed Sentinel-2 LAI Index (SeLI) is a normalized index that uses the new Sentinel-2 narrow B5-band located at the beginning of the red-edge region (705 nm), a spectral area which balances the influence of strong chlorophyll absorption and minimal scattering at moderate-high LAIgreen values, and a NIR band (865 nm) influenced by leaf scattering, as a reference band. SeLI was calibrated with a multi-crop Valencia dataset composed of LAIgreen values of 13 different crop types, obtaining a R2 of 0.708 and RMSE of 0.67. The validation of SeLI proved good, with a R2 of 0.732 and RMSE of 0.69. This work demonstrated the existence of a linear relationship between in situ LAIgreen and the spectral information of Sentinel-2, in the range of 0 to 5. Sentinel-2 generated maps over the Valencia and Foggia test sites convincingly illustrate the great potential of high spatial resolution LAIgreen monitoring at the single agricultural plot level, even for small- and medium-scale agricultural activities.

Author Contributions

N.P., J.D. and M.R. designed the field campaigns and collected in situ LAI multi-crop data; N.P. and J.D. analyzed the data and obtained the results; N.P. and S.V.W. wrote the paper and J.M. designed the experiments and guided and supervised the work.

Funding

This research was funded by SENSAGRI (Sentinels Synergy for Agriculture) project (H2020-EO-2016-730074).

Acknowledgments

This work was supported by the pre-doctoral scholarship of the Generalitat Valenciana Vali+d (ACIF/2016/378) of the first author. S.V.W. acknowledges the European Union’s Marie Sklodowska-Curie fellowship (MSCA-IF-2015-701815) and the Generalitat Valenciana Vali+d (APOSTD/2018/162). We would also like to thank the SENSAGRI project (H2020-EO-2016-730074) and especially to the EOLAB spin-off from the University of Valencia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Daughtry, C.S.T.; Gallo, K.P.; Goward, S.N.; Prince, S.D.; Kustas, W.P. Spectral estimates of absorbed radiation and phytomass production in corn and soybean canopies. Remote Sens. Environ. 1992, 39, 141–152. [Google Scholar] [CrossRef]
  2. Delegido, J.; Verrelst, J.; Rivera, J.P.; Ruiz-Verdú, A.; Moreno, J. Brown and green LAI mapping through spectral indices. Int. J. Appl. Earth Obs. Geoinf. 2015, 35, 350–358. [Google Scholar] [CrossRef]
  3. Zaroug, M.A.H.; Sylla, M.B.; Giorgi, F.; Eltahir, E.A.B.; Aggarwal, P.K. A sensitivity study on the role of the swamps of southern Sudan in the summer climate of North Africa using a regional climate model. Theor. Appl. Climatol. 2013, 113, 63–81. [Google Scholar] [CrossRef]
  4. Richardson, A.D.; Dail, D.B.; Hollinger, D.Y. Leaf area index uncertainty estimates for model-data fusion applications. Agric. For. Meteorol. 2011, 151, 1287–1292. [Google Scholar] [CrossRef]
  5. Bréda, N.J.J. Ground-based measurements of leaf area index: A review of methods, instruments and current controversies. J. Exp. Bot. 2003, 54, 2403–2417. [Google Scholar] [CrossRef] [PubMed]
  6. Sakamoto, T.; Gitelson, A.; Nguy-Robertson, A.; Arkebauer, T.; Wardlow, B.; Suyker, A.; Verma, S.; Shibayama, M. An alternative method using digital cameras for continuous monitoring of crop status. Agric. For. Meteorol. 2012, 154, 113–126. [Google Scholar] [CrossRef]
  7. Coops, N.C.; Waring, R.H.; Hilker, T. Prediction of soil properties using a process-based forest growth model to match satellite-derived estimates of leaf area index. Remote Sens. Environ. 2012, 126, 160–173. [Google Scholar] [CrossRef]
  8. Boegh, E.; Soegaard, H.; Broge, N.; Hasager, C.B.; Jensen, N.O.; Schelde, K.; Thomsen, A. Airborne multispectral data for quantifying leaf area index, nitrogen concentration, and photosynthetic efficiency in agriculture. Remote Sens. Environ. 2002, 81, 179–193. [Google Scholar] [CrossRef]
  9. Kira, O.; Nguy-Robertson, A.L.; Arkebauer, T.J.; Linker, R.; Gitelson, A.A. Informative spectral bands for remote green LAI estimation in C3 and C4 crops. Agric. For. Meteorol. 2016, 218–219, 243–249. [Google Scholar] [CrossRef]
  10. Clevers, J.G.P.W.; Kooistra, L.; van den Brande, M.M.M. Using Sentinel-2 data for retrieving LAI and leaf and canopy chlorophyll content of a potato crop. Remote Sens. 2017, 9, 405. [Google Scholar] [CrossRef]
  11. Brisco, B.; Brown, R.J.; Hirose, T.; McNairn, H.; Staenz, K. Precision agriculture and the role of remote sensing: A review. Can. J. Remote Sens. 1998, 24, 315–327. [Google Scholar] [CrossRef]
  12. Houlès, V.; Guérif, M.; Mary, B. Elaboration of a nitrogen nutrition indicator for winter wheat based on leaf area index and chlorophyll content for making nitrogen recommendations. Eur. J. Agron. 2007, 27, 1–11. [Google Scholar] [CrossRef]
  13. Bsaibes, A.; Courault, D.; Baret, F.; Weiss, M.; Olioso, A.; Jacob, F.; Hagolle, O.; Marloie, O.; Bertrand, N.; Desfond, V.; et al. Albedo and LAI estimates from FORMOSAT-2 data for crop monitoring. Remote Sens. Environ. 2009, 113, 716–729. [Google Scholar] [CrossRef]
  14. Yao, X.; Wang, N.; Liu, Y.; Cheng, T.; Tian, Y.; Chen, Q.; Zhu, Y. Estimation of wheat LAI at middle to high levels using unmanned aerial vehicle narrowband multispectral imagery. Remote Sens. 2017, 9, 1304. [Google Scholar] [CrossRef]
  15. Campos-Taberner, M.; García-Haro, F.J.; Camps-Valls, G.; Grau-Muedra, G.; Nutini, F.; Crema, A.; Boschetti, M. Multitemporal and multiresolution leaf area index retrieval for operational local rice crop monitoring. Remote Sens. Environ. 2016, 187, 102–118. [Google Scholar] [CrossRef]
  16. Mulla, D.J. Twenty five years of remote sensing in precision agriculture: Key advances and remaining knowledge gaps. Biosyst. Eng. 2013, 114, 358–371. [Google Scholar] [CrossRef]
  17. Drusch, M.; Del Bello, U.; Carlier, S.; Colin, O.; Fernandez, V.; Gascon, F.; Hoersch, B.; Isola, C.; Laberinti, P.; Martimort, P.; et al. Sentinel-2: ESA’s optical high-resolution mission for GMES operational services. Remote Sens. Environ. 2012, 120, 25–36. [Google Scholar] [CrossRef]
  18. Baret, F.; Guyot, G. Potential and limitations of vegetation indices for LAI and APAR assessment. Remote Sens. Environ. 1991, 35, 161–173. [Google Scholar] [CrossRef]
  19. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  20. Verrelst, J.; Romijn, E.; Kooistra, L. Mapping vegetation density in a heterogeneous river floodplain ecosystem using pointable CHRIS/PROBA data. Remote Sens. 2012, 4, 2866–2889. [Google Scholar] [CrossRef]
  21. Broge, N.H.; Mortensen, J.V. Deriving green crop area index and canopy chlorophyll density of winter wheat from spectral reflectance data. Remote Sens. Environ. 2002, 81, 45–57. [Google Scholar] [CrossRef]
  22. Glenn, E.P.; Huete, A.R.; Nagler, P.L.; Nelson, S.G. Relationship between remotely-sensed vegetation indices, canopy attributes and plant physiological processes: What vegetation indices can and cannot tell us about the landscape. Sensors 2008, 8, 2136–2160. [Google Scholar] [CrossRef] [PubMed]
  23. Myneni, R.B.; Maggion, S.; Iaquinta, J.; Privette, J.L.; Gobron, N.; Pinty, B.; Kimes, D.S.; Verstraete, M.M.; William, D.L. Optical remote sensing of vegetation: Modeling, caveates, and algorithms. Remote Sens. Environ. 1995, 51, 169–188. [Google Scholar] [CrossRef]
  24. Gobron, N.; Pinty, B.; Verstraete, M.M.; Widlowski, J.L. Advanced vegetation indices optimized for up-coming sensors: Design, performance, and applications. IEEE Trans. Geosci. Remote Sens. 2000, 38, 2489–2505. [Google Scholar] [CrossRef]
  25. Houborg, R.; Boegh, E. Mapping leaf chlorophyll and leaf area index using inverse and forward canopy reflectance modeling and SPOT reflectance data. Remote Sens. Environ. 2008, 112, 186–202. [Google Scholar] [CrossRef]
  26. Jacquemoud, S.; Baret, F.; Andrieu, B.; Danson, M.; Jaggard, K. Extraction of vegetation biophysical parameters by inversion of the PROSPECT+SAIL model on sugar beet canopy reflectance data: Application to TM and AVIRIS sensors. Remote Sens. Environ. 1995, 52, 163–172. [Google Scholar] [CrossRef]
  27. Zheng, G.; Moskal, L.M. Retrieving leaf area index (LAI) using remote sensing: Theories, methods and sensors. Sensors 2009, 9, 2719–2745. [Google Scholar] [CrossRef]
  28. González-Sanpedro, M.C.; Le Toan, T.; Moreno, J.; Kergoat, L.; Rubio, E. Seasonal variations of leaf area index of agricultural fields retrieved from Landsat data. Remote Sens. Environ. 2008, 112, 810–824. [Google Scholar] [CrossRef]
  29. Combal, B.; Baret, F.; Weiss, M.; Trubuil, A.; Macé, D.; Pragnère, A.; Myneni, R.; Knyazikhin, Y.; Wang, L. Retrieval of canopy biophysical variables from bidirectional reflectance: Using prior information to solve the ill-posed inverse problem. Remote Sens. Environ. 2002, 84, 1–15. [Google Scholar] [CrossRef]
  30. Atzberger, C. Object-based retrieval of biophysical canopy variables using artificial neural nets and radiative transfer models. Remote Sens. Environ. 2004, 93, 53–67. [Google Scholar] [CrossRef]
  31. España, M.L.; Baret, F.; Aries, F.; Chelle, M.; Andrieu, B.; Prévot, L. Modeling maize canopy 3D architecture: Application to reflectance simulation. Ecol. Model. 1999, 122, 25–43. [Google Scholar] [CrossRef]
  32. Casa, R.; Baret, F.; Buis, S.; Lopez-Lozano, R.; Pascucci, S.; Palombo, A.; Jones, H.G. Estimation of maize canopy properties from remote sensing by inversion of 1-D and 4-D models. Precis. Agric. 2010, 11, 319–334. [Google Scholar] [CrossRef] [Green Version]
  33. Weiss, M.; Baret, F. S2ToolBox Level 2 products: LAI, FAPAR, FCOVER. Available online: http://step.esa.int/docs/extra/ATBD_S2ToolBox_L2B_V1.1.pdf (accessed on 9 November 2018).
  34. Verhoef, W. Light scattering by leaf layers with application to canopy reflectance modeling: The SAIL model. Remote Sens. Environ. 1984, 16, 125–141. [Google Scholar] [CrossRef] [Green Version]
  35. Jacquemoud, S.; Baret, F. PROSPECT: A model of leaf optical properties spectra. Remote Sens. Environ. 1990, 34, 75–91. [Google Scholar] [CrossRef]
  36. Djamai, N.; Fernandes, R. Comparison of SNAP-derived Sentinel-2A L2A product to ESA product over Europe. Remote Sens. 2018, 10, 926. [Google Scholar] [CrossRef]
  37. Verrelst, J.; Muñoz, J.; Alonso, L.; Delegido, J.; Rivera, J.P.; Camps-Valls, G.; Moreno, J. Machine learning regression algorithms for biophysical parameter retrieval: Opportunities for Sentinel-2 and -3. Remote Sens. Environ. 2012, 118, 127–139. [Google Scholar] [CrossRef]
  38. Cui, Z.; Kerekes, J.P. Potential of red edge spectral bands in future landsat satellites on agroecosystem canopy green leaf area index retrieval. Remote Sens. 2018, 10, 1458. [Google Scholar] [CrossRef]
  39. Verrelst, J.; Camps-Valls, G.; Muñoz-Marí, J.; Rivera, J.P.; Veroustraete, F.; Clevers, J.G.P.W.; Moreno, J. Optical remote sensing and the retrieval of terrestrial vegetation bio-geophysical properties—A review. ISPRS J. Photogramm. Remote Sens. 2015, 108, 273–290. [Google Scholar] [CrossRef]
  40. le Maire, G.; François, C.; Soudani, K.; Berveiller, D.; Pontailler, J.Y.; Bréda, N.; Genet, H.; Davi, H.; Dufrêne, E. Calibration and validation of hyperspectral indices for the estimation of broadleaved forest leaf chlorophyll content, leaf mass per area, leaf area index and leaf canopy biomass. Remote Sens. Environ. 2008, 112, 3846–3864. [Google Scholar] [CrossRef]
  41. Jordan, C.F. Derivation of leaf area index from quality of light on the forest floor. Ecology 1969, 50, 663–666. [Google Scholar] [CrossRef]
  42. Delegido, J.; Verrelst, J.; Alonso, L.; Moreno, J. Evaluation of sentinel-2 red-edge bands for empirical estimation of green LAI and chlorophyll content. Sensors 2011, 11, 7063–7081. [Google Scholar] [CrossRef] [PubMed]
  43. Rouse, J.W.; Hass, R.H.; Schell, J.A.; Deering, D.W. Monitoring vegetation systems in the great plains with ERTS. In Third Earth Resources Technology Satellite-1 Symposium; NASA: Washington, DC, USA, 1973; Volume 1, pp. 309–317. [Google Scholar]
  44. Guyot, G.; Baret, F.; Jacquemond, S. Imaging spectroscopy for vegetation studies. Imaging Spectrosc. Fundam. Prospect. Appl. 1992, 2, 145–165. [Google Scholar]
  45. Baret, F.; Jacquemoud, S.; Guyot, G.; Leprieur, C. Modeled analysis of the biophysical nature of spectral shifts and comparison with information content of broad bands. Remote Sens. Environ. 1992, 41, 133–142. [Google Scholar] [CrossRef]
  46. Liu, J.; Miller, J.R.; Haboudane, D.; Pattey, E. Exploring the relationship between red edge parameters and crop variables for precision agriculture. In Proceedings of the IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium, Anchorage, AK, USA, 20–24 September 2004; Volume 2, pp. 1276–1279. [Google Scholar] [CrossRef]
  47. Delegido, J.; Fernandez, G.; Gandia, S.; Moreno, J. Retrieval of chlorophyll content and LAI of crops using hyperspectral techniques: Application to PROBA/CHRIS data. Int. J. Remote Sens. 2008, 29, 7107–7127. [Google Scholar] [CrossRef]
  48. Herrmann, I.; Pimstein, A.; Karnieli, A.; Cohen, Y.; Alchanatis, V.; Bonfil, D.J. LAI assessment of wheat and potato crops by VENμS and Sentinel-2 bands. Remote Sens. Environ. 2011, 115, 2141–2151. [Google Scholar] [CrossRef]
  49. Campos-Taberner, M.; García-Haro, F.J.; Camps-Valls, G.; Grau-Muedra, G.; Nutini, F.; Busetto, L.; Katsantonis, D.; Stavrakoudis, D.; Minakou, C.; Gatti, L.; et al. Exploitation of SAR and optical sentinel data to detect rice crop and estimate seasonal dynamics of leaf area index. Remote Sens. 2017, 9, 248. [Google Scholar] [CrossRef]
  50. Kira, O.; Nguy-Robertson, A.L.; Arkebauer, T.J.; Linker, R.; Gitelson, A.A. Toward generic models for green LAI estimation in maize and soybean: Satellite observations. Remote Sens. 2017, 9, 318. [Google Scholar] [CrossRef]
  51. Nguy-Robertson, A.L.; Peng, Y.; Gitelson, A.A.; Arkebauer, T.J.; Pimstein, A.; Herrmann, I.; Karnieli, A.; Rundquist, D.C.; Bonfil, D.J. Estimating green LAI in four crops: Potential of determining optimal spectral bands for a universal algorithm. Agric. For. Meteorol. 2014, 192–193, 140–148. [Google Scholar] [CrossRef]
  52. Xie, Q.; Dash, J.; Huang, W.; Peng, D.; Qin, Q.; Mortimer, H.; Casa, R.; Pignatti, S.; Laneve, G.; Pascucci, S.; et al. Vegetation indices combining the red and red-edge spectral information for Leaf Area Index retrieval. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 1482–1493. [Google Scholar] [CrossRef]
  53. AEMET. Agencia Estatal de Meteorología. Available online: http://www.aemet.es/es/portada (accessed on 20 November 2018).
  54. METEOBLUE Weather. Climate Foggia. Available online: https://www.meteoblue.com/en/weather/forecast/modelclimate/foggia_italy_3176885 (accessed on 20 November 2018).
  55. LI-COR. Licor 2200 Instruction Manual; LI-COR: Lincoln, NE, USA, 2012; ISBN 4024672819. [Google Scholar]
  56. Pandžić, M.; Mihajlović, D.; Pandžić, J.; Pfeifer, N. Assessment of the geometric quality of sentinel-2 data. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. ISPRS Arch. 2016, 41, 489–494. [Google Scholar] [CrossRef]
  57. Louis, J.; Debaecker, V.; Pflug, B.; Main-Knorn, M.; Bieniarz, J.; Mueller-Wilm, U.; Cadau, E.; Gascon, F. Sentinel-2 SEN2COR: L2A processor for users. In Proceedings of the Living Planet Symposium, Prague, Czech Republic, 9–13 May 2016; pp. 1–8. [Google Scholar]
  58. Verrelst, J.; Rivera, J.P.; Alonso, L.; Moreno, J. ARTMO: An Automated Radiative Transfer Models Operator toolbox for automated retrieval of biophysical parameters through model inversion. In Proceedings of the EARSeL 7th SIG-Imaging Spectroscopy Workshop, Edinburgh, UK, 11–13 April 2011; pp. 11–13. [Google Scholar]
  59. Rivera, J.; Verrelst, J.; Delegido, J.; Veroustraete, F.; Moreno, J. On the semi-automatic retrieval of biophysical parameters based on spectral index optimization. Remote Sens. 2014, 6, 4927–4951. [Google Scholar] [CrossRef]
  60. Snee, R. Validation and regression models: Methods and examples. Technometrics 1977, 19, 415–428. [Google Scholar] [CrossRef]
  61. Vincini, M.; Frazzi, E.; D’Alessio, P. Comparison of narrow-band and broad-band vegetation indices for canopy chlorophyll density estimation in sugar beet. In Proceedings of the 6th European Conference on Precision Agriculture, Skiathos, Greece, 3–6 June 2007; pp. 189–196. [Google Scholar]
  62. Merzlyak, M.N.; Gitelson, A.A.; Chivkunova, O.B.; Rakitin, V.Y. Non-destructive optical detection of leaf senescence and fruit ripening. Physiol. Plant. 1999, 106, 135–141. [Google Scholar] [CrossRef]
  63. Dash, J.; Curran, P.J. The MERIS terrestrial chlorophyll index. Int. J. Remote Sens. 2004, 2523, 5403–5413. [Google Scholar] [CrossRef]
  64. Sims, D.A.; Gamon, J.A. Relationships between leaf pigment content and spectral reflectance across a wide range of species, leaf structures and developmental stages. Remote Sens. Environ. 2002, 81, 337–354. [Google Scholar] [CrossRef]
  65. Gitelson, A.A.; Kaufman, Y.J.; Stark, R.; Rundquist, D. Novel algorithms for remote estimation of vegetation fraction. Remote Sens. Environ. 2002, 80, 76–87. [Google Scholar] [CrossRef] [Green Version]
  66. Gower, J.F.R. Observations of in situ fluorescence of chlorophyll-a in Saanich Inlet. Bound. Layer Meteorol. 1980, 18, 235–245. [Google Scholar] [CrossRef]
  67. Broge, N.H.; Leblanc, E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sens. Environ. 2000, 76, 156–172. [Google Scholar] [CrossRef]
  68. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; De Colstoun, E.B.; McMurtrey, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  69. Pasqualotto, N.; Delegido, J.; Van Wittenberghe, S.; Verrelst, J.; Rivera, J.P.; Moreno, J. Retrieval of canopy water content of different crop types with two new hyperspectral indices: Water Absorption Area Index and Depth Water Index. Int. J. Appl. Earth Obs. Geoinf. 2018, 67, 69–78. [Google Scholar] [CrossRef]
  70. Dawson, T.P.; Curran, P.J.; North, P.R.J.; Plummer, S.E. The propagation of foliar biochemical absorption features in forest canopy reflectance: A theoretical analysis. Remote Sens. Environ. 1999, 67, 147–159. [Google Scholar] [CrossRef]
  71. Clevers, J.G.P.W.; Gitelson, A.A. Remote estimation of crop and grass chlorophyll and nitrogen content using red-edge bands on Sentinel-2 and -3. Int. J. Appl. Earth Obs. Geoinf. 2013, 23, 344–351. [Google Scholar] [CrossRef]
  72. Gitelson, A.A. Wide dynamic range vegetation index for remote quantification of biophysical characteristics of vegetation. J. Plant Physiol. 2004, 161, 165–173. [Google Scholar] [CrossRef] [PubMed]
  73. Liang, L.; Di, L.; Zhang, L.; Deng, M.; Qin, Z.; Zhao, S.; Lin, H. Estimation of crop LAI using hyperspectral vegetation indices and a hybrid inversion method. Remote Sens. Environ. 2015, 165, 123–134. [Google Scholar] [CrossRef]
  74. Verrelst, J.; Rivera, J.P.; Veroustraete, F.; Muñoz-Marí, J.; Clevers, J.G.P.W.; Camps-Valls, G.; Moreno, J. Experimental Sentinel-2 LAI estimation using parametric, non-parametric and physical retrieval methods—A comparison. ISPRS J. Photogramm. Remote Sens. 2015, 108, 260–272. [Google Scholar] [CrossRef]
  75. Tanaka, S.; Kawamura, K.; Maki, M.; Muramoto, Y.; Yoshida, K.; Akiyama, T. Spectral index for quantifying leaf area index of winter wheat by field hyperspectral measurements: A case study in Gifu Prefecture, Central Japan. Remote Sens. 2015, 7, 5329–5346. [Google Scholar] [CrossRef]
  76. Duveiller, G.; Baret, F.; Defourny, P. Remotely sensed green area index for winter wheat crop monitoring: 10-Year assessment at regional scale over a fragmented landscape. Agric. For. Meteorol. 2012, 166–167, 156–168. [Google Scholar] [CrossRef]
  77. Fernandes, R.; Omari, K.; Canisius, F.; Maloley, M.; Rochdi, N. A theoretical basis for the observed linear relationship between LAI and normalized differences of red edge reflectance with implications for LAI retrieval from Sentinel-2. In Proceedings of the Recent Advances in Quantitative Remote Sensing RAQRS, Torrent, Spain, 18–22 September 2017. [Google Scholar]
  78. Darvishzadeh, R.; Skidmore, A.; Atzberger, C.; van Wieren, S. Estimation of vegetation LAI from hyperspectral reflectance data: Effects of soil type and plant architecture. Int. J. Appl. Earth Obs. Geoinf. 2008, 10, 358–373. [Google Scholar] [CrossRef]
  79. Wu, C.; Han, X.; Niu, Z.; Dong, J. An evaluation of EO-1 hyperspectral Hyperion data for chlorophyll content and leaf area index estimation. Int. J. Remote Sens. 2010, 31, 1079–1086. [Google Scholar] [CrossRef]
  80. Fernandes, R.; Omari, K.; Canisius, F.; Rochdi, N.; Baret, F. Robust leaf area index retrieval using Sentinel-2 red edge bands. In Proceedings of the First Sentinel-2 Preparatory Symposium, Frascati, Italy, 23–27 April 2012. [Google Scholar]
  81. Gower, S.T.; Kucharik, C.J.; Norman, J.M. Direct and indirect estimation of leaf area index, f(APAR), and net primary production of terrestrial ecosystems. Remote Sens. Environ. 1999, 70, 29–51. [Google Scholar] [CrossRef]
  82. Küßner, R.; Mosandl, R. Comparison of direct and indirect estimation of leaf area index in mature Norway spruce stands of eastern Germany. Can. J. For. Res. 2000, 30, 440–447. [Google Scholar] [CrossRef]
  83. Leblanc, S.G.; Chen, J.M. A practical scheme for correcting multiple scattering effects on optical LAI measurements. Agric. For. Meteorol. 2001, 110, 125–139. [Google Scholar] [CrossRef] [Green Version]
  84. Vuolo, F.; Zóltak, M.; Pipitone, C.; Zappa, L.; Wenng, H.; Immitzer, M.; Weiss, M.; Baret, F.; Atzberger, C. Data service platform for Sentinel-2 surface reflectance and value-added products: System use and examples. Remote Sens. 2016, 8, 938. [Google Scholar] [CrossRef]
  85. Vanino, S.; Nino, P.; De Michele, C.; Falanga, S.; Urso, G.D.; Di, C.; Pennelli, B.; Vuolo, F.; Farina, R.; Pulighe, G. Capability of Sentinel-2 data for estimating maximum evapotranspiration and irrigation requirements for tomato crop in Central Italy. Remote Sens. Environ. 2018, 215, 452–470. [Google Scholar] [CrossRef]
  86. Ollinger, S.V. Sources of variability in canopy reflectance and the convergent properties of plants. New Phytol. 2011, 189, 375–394. [Google Scholar] [CrossRef] [PubMed]
  87. Van Wittenberghe, S.; Verrelst, J.; Rivera, J.P.; Alonso, L.; Moreno, J.; Samson, R. Gaussian processes retrieval of leaf parameters from a multi-species reflectance, absorbance and fluorescence dataset. J. Photochem. Photobiol. B Biol. 2014, 134, 37–48. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Test site locations. (a) Valencia (Spain) from the Sentinel-2 image of 6 April 2017, UTM-WGS84 Zone 30S (b) Foggia (Italy) from the Sentinel-2 image of 8 April 2017, UTM-WGS84 Zone 33T.
Figure 1. Test site locations. (a) Valencia (Spain) from the Sentinel-2 image of 6 April 2017, UTM-WGS84 Zone 30S (b) Foggia (Italy) from the Sentinel-2 image of 8 April 2017, UTM-WGS84 Zone 33T.
Sensors 19 00904 g001
Figure 2. Sampling approach for each elementary sampling unit (ESU).
Figure 2. Sampling approach for each elementary sampling unit (ESU).
Sensors 19 00904 g002
Figure 3. ESUs frequency histogram classified for LAIgreen value range, for (a) the VLC17_ES testing dataset and (b) the FOG17_IT validation dataset. Bare soil ESUs are not included.
Figure 3. ESUs frequency histogram classified for LAIgreen value range, for (a) the VLC17_ES testing dataset and (b) the FOG17_IT validation dataset. Bare soil ESUs are not included.
Sensors 19 00904 g003
Figure 4. Green leaf area index (LAIgreen) in situ values in comparison to the Sentinel-2 LAIgreen product, with 1:1 line (black), its statistics and linear fitting (red).
Figure 4. Green leaf area index (LAIgreen) in situ values in comparison to the Sentinel-2 LAIgreen product, with 1:1 line (black), its statistics and linear fitting (red).
Sensors 19 00904 g004
Figure 5. Graphic representation of the LAIgreen of VLC17_ES dataset as a function of the established LAIgreen indices, with the best fitting function and 95% of confidence interval. (a) RI [41]; (b) NDGI [43]; (c) NDGI [42]; (d) TRBI [61].
Figure 5. Graphic representation of the LAIgreen of VLC17_ES dataset as a function of the established LAIgreen indices, with the best fitting function and 95% of confidence interval. (a) RI [41]; (b) NDGI [43]; (c) NDGI [42]; (d) TRBI [61].
Sensors 19 00904 g005aSensors 19 00904 g005b
Figure 6. Representation of real Sentinel-2 top-of-canopy (TOC) reflectance spectra of (a) artichoke crops (20/07/2017 Sentinel-2 image used) and (b) orange tree crops (07/11/2017 Sentinel-2 image used). The red and red-edge Sentinel-2 bands (B4, B5 and B6) are represented with their corresponding bandwidth.
Figure 6. Representation of real Sentinel-2 top-of-canopy (TOC) reflectance spectra of (a) artichoke crops (20/07/2017 Sentinel-2 image used) and (b) orange tree crops (07/11/2017 Sentinel-2 image used). The red and red-edge Sentinel-2 bands (B4, B5 and B6) are represented with their corresponding bandwidth.
Sensors 19 00904 g006
Figure 7. (a) Graphic representation of the LAIgreen of the testing VLC17_ES dataset (n = 92) as a function of SeLI, the index proposed, with linear fit and 95% of confidence interval, (b) LAIgreen in situ validation data from FOG17_IT (n = 109) as a function of LAIgreen estimated with SeLI, with 1:1 line and its statistics.
Figure 7. (a) Graphic representation of the LAIgreen of the testing VLC17_ES dataset (n = 92) as a function of SeLI, the index proposed, with linear fit and 95% of confidence interval, (b) LAIgreen in situ validation data from FOG17_IT (n = 109) as a function of LAIgreen estimated with SeLI, with 1:1 line and its statistics.
Sensors 19 00904 g007
Figure 8. LAIgreen maps estimated with SeLI proposed index in Valencia test site using the Sentinel-2 image of 26 May 2017 (left) and Foggia test site using the Sentinel-2 image of 8 April 2017 (right). Note the different scale of both study sites and their respective agricultural areas.
Figure 8. LAIgreen maps estimated with SeLI proposed index in Valencia test site using the Sentinel-2 image of 26 May 2017 (left) and Foggia test site using the Sentinel-2 image of 8 April 2017 (right). Note the different scale of both study sites and their respective agricultural areas.
Sensors 19 00904 g008
Table 1. Sentinel-2 bands setting.
Table 1. Sentinel-2 bands setting.
Band NumberFunctionCentral Wavelength (nm)Bandwidth (nm)Spatial Resolution (m)
1Coastal aerosol4432760
2Blue4909810
3Green5604510
4Red6653810
5Vegetation red-edge7051920
6Vegetation red-edge7401820
7Vegetation red-edge7832820
8Near infrared (NIR)84214510
8aVegetation red-edge8653320
9Water vapour9452660
10Shortwave infrared (SWIR)-cirrus13807560
11SWIR161014320
12SWIR219024220
Table 2. Sentinel-2 images used in each field campaign.
Table 2. Sentinel-2 images used in each field campaign.
LocationField Work Date (2017 year)Sentinel-2 Image Code
Valencia22 and 23 May S2A_MSIL1C_20170526T105031_N0205_R051_T30SYJ_20170526T105518
18 and 19 JulyS2A_MSIL1C_20170720T105029_N0205_R051_T30SYJ_20170720T105641
8 and 9 NovemberS2A_MSIL1C_20171107T105229_N0206_R051_T30SYJ_20171107T131035
Foggia16 MarchS2A_MSIL1C_20170319T095021_N0204_R079_T33TWF_20170319T095021
21 and 22 MarchS2A_MSIL1C_20170319T095021_N0204_R079_T33TWG_20170319T095021
29 MarchS2A_MSIL1C_20170329T095021_N0204_R079_T33TWF_20170329T095024
5 and 13 AprilS2A_MSIL1C_20170408T095031_N0204_R079_T33TWF_20170408T095711
11 and 17 MayS2A_MSIL1C_20170518T095031_N0205_R079_T33TWF_20170518T095716
3 MayS2A_MSIL1C_20170528T095031_N0205_R079_T33TWF_20170528T095531
12 JuneS2A_MSIL1C_20170607T095031_N0205_R079_T33TWF_20170607T095031
15 and 21 JuneS2A_MSIL1C_20170617T095031_N0205_R079_T33TWF_20170617T095546
Table 3. Generic vegetation indices introduced in Automated Radiative Transfer Models Operator (ARTMO), where indices based on commonly LAIgreen indices are shown shaded and typically indices used to estimate other biophysical parameters are shown non-shaded. Rλ represents reflectance at the wavelength λ. The generic name of each index has been established in this study.
Table 3. Generic vegetation indices introduced in Automated Radiative Transfer Models Operator (ARTMO), where indices based on commonly LAIgreen indices are shown shaded and typically indices used to estimate other biophysical parameters are shown non-shaded. Rλ represents reflectance at the wavelength λ. The generic name of each index has been established in this study.
Based ReferenceFormulaGeneric Name AbbreviationGeneric Formula
[41] R 800 R 675 Ratio IndexRI R 1 R 2
[43] R 800 R 670 R 800 + R 670 Normalized Difference Generic IndexNDGI R 1 R 2 R 1 +   R 2
[42] R 704 R 665 R 704 + R 665
[61] R 520 600   + R 630 690 R 760 900 Three Ratio Band IndexTRBI R 1 + R 2 R 3
[62] R 680   R 500 R 750 Three Difference Band IndexTDBI R 1 R 2 R 3
[63] R 750 R 710 R 710 R 680 MERIS Terrestrial Generic IndexMTGI R 1 R 2 R 2 R 3
Normalized Difference 3 bandND3b R 1 R 2 R 2 + R 3
[64] R 750 R 705 R 750 + R 705 2 R 445 Multi-band Normalized IndexMNI R 1 R 2 R 1 + R 2 R 3
[65] R 550 R 670 R 550 + R 670 R 480
[66] R 676 0.5 ( R 746 + R 665 ) Generic Line HeightGLH R 1 0.5 ( R 2 + R 3 )
[67] 0.5 [ 120 ( R 750 R 550 ) 200 ( R 670 R 550 ) ]Triangular Generic IndexTGI 0.5 [ 120 ( R 1 R 2 ) 200 ( R 3 R 2 ) ]
[68] [ ( R 700 R 670 ) 0.2 ( R 700 R 550 ) ] ( R 700 / R 670 ) Modified Chlorophyll Generic IndexMCGI [ ( R 1 R 2 ) 0.2 ( R 1 R 3 ) ] ( R 1 / R 2 )
Table 4. Statistics obtained with a linear, polynomial of second order, exponential and exponential of second order fitting for each index. The best fitting is boldfaced.
Table 4. Statistics obtained with a linear, polynomial of second order, exponential and exponential of second order fitting for each index. The best fitting is boldfaced.
IndexReferencesLinear FittingPolynomial Fitting, Second OrderExponential FittingExponential Fitting, Second Order
R2RMSER2RMSER2RMSER2RMSE
RI [41]0.3550.910.3140.930.2341.130.6630.75
y = 0.15x + 1.11y = −0.01x2 + 0.33x + 0.55y = 1.51exp(0.04x)y = 2.59exp(0.01x) − 6.19exp(-0.71x)
NDGI [43]0.6590.720.6120.740.5710.850.6290.79
y = 3.93x − 0.18y = −1.98x2 + 5.93x − 0.55y = 0.68exp(1.78x)y = − 1547exp(3.96x) + 1548exp(3.96x)
[42]0.4020.860.3890.890.3101.070.5490.88
y = 3.58x + 0.91y = − 6.51x2 + 9.33x + 0.16y = 1.33exp(1.18x)y = −3.61exp(−0.08x) − 3.84exp(−4.21x)
TRBI [61]0.6630.750.6390.750.6250.790.6590.76
y = −2.55x + 3.62y = 0.44x2 − 3.27x + 3.84y = 4.39exp(−1.42x)y = −7.13exp(−3.34x) − exp(−3.34x)
Table 5. Best band combination for each generic vegetation index introduced in ARTMO from cross-validation of the testing dataset (VLC17_ES), ordered from highest to lowest R2, with a linear fit.
Table 5. Best band combination for each generic vegetation index introduced in ARTMO from cross-validation of the testing dataset (VLC17_ES), ordered from highest to lowest R2, with a linear fit.
Index BandsR2RMSENRMSE (%)p-Value
TRBI2190;740;8650.7370.6314<0.001
TDBI2190;865;7400.7320.6415<0.001
ND3b2190;865;7400.7310.6415<0.001
MNI2190;865;16100.7310.6415<0.001
RI2190;8650.7280.6415<0.001
MCGI1610;865;7400.7170.6515<0.001
TGI1610;842;21900.7130.6615<0.001
GLH2190;1610;8420.7080.6715<0.001
NDGI865;7050.7080.6715<0.001
MTGI2190;865;4900.7010.6815<0.001
Table 6. Linear cross-validation fitting results of the normalized difference ratio index with different band combinations and LAIgreen values from the testing and validation datasets, ordered from highest to lowest R2 according to the testing dataset.
Table 6. Linear cross-validation fitting results of the normalized difference ratio index with different band combinations and LAIgreen values from the testing and validation datasets, ordered from highest to lowest R2 according to the testing dataset.
BandsTesting (VLC17_ES)Validation (FOG17_IT)
R2RMSER2RMSE
865;7050.7080.670.7320.69
783;7050.7020.680.7110.71
842;7050.6880.690.7170.71
740;7050.6850.710.6860.74
783;6650.6750.710.6780.75
842;6650.6650.720.6840.74
783;7400.5310.850.6740.76

Share and Cite

MDPI and ACS Style

Pasqualotto, N.; Delegido, J.; Van Wittenberghe, S.; Rinaldi, M.; Moreno, J. Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI). Sensors 2019, 19, 904. https://doi.org/10.3390/s19040904

AMA Style

Pasqualotto N, Delegido J, Van Wittenberghe S, Rinaldi M, Moreno J. Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI). Sensors. 2019; 19(4):904. https://doi.org/10.3390/s19040904

Chicago/Turabian Style

Pasqualotto, Nieves, Jesús Delegido, Shari Van Wittenberghe, Michele Rinaldi, and José Moreno. 2019. "Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI)" Sensors 19, no. 4: 904. https://doi.org/10.3390/s19040904

APA Style

Pasqualotto, N., Delegido, J., Van Wittenberghe, S., Rinaldi, M., & Moreno, J. (2019). Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI). Sensors, 19(4), 904. https://doi.org/10.3390/s19040904

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop