Influence of Different Satellite Imagery on the Analysis of Riparian Leaf Density in a Mountain Stream

: In recent decades, technological advancements in sensors have generated increasing interest in remote sensing data for the study of vegetation features. Image pixel resolution can a ﬀ ect data analysis and results. This study evaluated the potential of three satellite images of di ﬀ ering resolution (Landsat 8, 30 m; Sentinel-2, 10 m; and Pleiades 1A, 2 m) in assessing the Leaf Area Index (LAI) of riparian vegetation in two Mediterranean streams, and in both a winter wheat ﬁeld and a deciduous forest used to compare the accuracy of the results. In this study, three di ﬀ erent retrieval methods—the Caraux-Garson, the Lambert-Beer, and the Campbell and Norman equations—are used to estimate LAI from the Normalized Di ﬀ erence Vegetation Index (NDVI). To validate sensor data, LAI values were measured in the ﬁeld using the LAI 2200 Plant Canopy Analyzer. The statistical indices showed a better performance for Pleiades 1A and Landsat 8 images, the former particularly in sites characterized by high canopy closure, such as deciduous forests, or in areas with stable riparian vegetation, the latter where stable reaches of riparian vegetation cover are almost absent or very homogenous, as in winter wheat ﬁelds. Sentinel-2 images provided more accurate results in terms of the range of LAI values. Considering the di ﬀ erent types of satellite imagery, the Lambert-Beer equation generally performed best in estimating LAI from the NDVI, especially in areas that are geomorphologically stable or have a denser vegetation cover, such as deciduous forests.


Introduction
Remote sensing is a technique for acquiring information on the Earth using different types of sensors located at variable distances [1]. Ever since the first remote sensing imagery became available, it has been used for environmental applications, such as land use, hydrology, geology, weather, climate, and vegetation studies [2,3]. Satellites and aircrafts equipped with optical, radar, or laser sensors (i.e., satellite platforms) are the most common means for collecting data on and studying Earth's surface [4]. Along with these traditional remote sensing technologies, drones have recently been used to conduct environmental studies [5][6][7][8]. Satellite remote sensing, however, is still used to monitor large areas at regular intervals of time [9,10]. The main advantages of satellite imagery are that some platforms offer free access to data and the availability of different spatial and temporal resolutions [11]. The use of satellite imagery is, therefore, essential in contexts where the size of the area makes in-field monitoring expensive and time-consuming, as in the case of watersheds [12].
There are currently numerous satellite platforms operating; these carry different types of sensors with spatial and temporal resolutions that range from sub-metric to kilometric and from 30 min to weeks or months, respectively [13]. AVHRR, Sentinel, MODIS, Landsat, SPOT, Pleiades, and Ikonos are among the satellite platforms most commonly used to monitor vegetation features [13][14][15][16]. Given the wide variety of satellite imagery, it is important to choose the one most suited to the objectives of Downy oak and hardwood (Quercus pubescens Wild and Quercus. Cerris L., respectively), along with pasture, meadow, and coniferous species, cover the higher slopes. Wheat cultivation, and to a lesser extent, agricultural crops and olive groves are the predominant cultivation systems in the plain and low hilly areas [48]. LAI values were measured in-situ at the Cammarota (CR) and Vallone della Madonna (R5) streams, both tributaries of the Carapelle Torrent, and in two surrounding areas, one covered by winter wheat (WW) fields and the other by a deciduous forest (DF). The watercourses, with narrow cross-sections and stratified vegetation, differ from one another in geomorphological stability, and consequently, in vegetation density. The CR is characterized by different vegetation patterns and degrees of erosion and by the presence of numerous check dams, which are severely damaged or destroyed. Along the monitored transect of the CR, there are 28 check dams, of which 15 are Downy oak and hardwood (Quercus pubescens Wild and Quercus. Cerris L., respectively), along with pasture, meadow, and coniferous species, cover the higher slopes. Wheat cultivation, and to a Remote Sens. 2020, 12, 3376 4 of 16 lesser extent, agricultural crops and olive groves are the predominant cultivation systems in the plain and low hilly areas [48].
LAI values were measured in-situ at the Cammarota (CR) and Vallone della Madonna (R5) streams, both tributaries of the Carapelle Torrent, and in two surrounding areas, one covered by winter wheat (WW) fields and the other by a deciduous forest (DF). The watercourses, with narrow cross-sections and stratified vegetation, differ from one another in geomorphological stability, and consequently, in vegetation density. The CR is characterized by different vegetation patterns and degrees of erosion and by the presence of numerous check dams, which are severely damaged or destroyed. Along the monitored transect of the CR, there are 28 check dams, of which 15 are destroyed, 3 are damaged, and 10 are intact. The destroyed check dams are the main cause of active erosion processes along the stream [27]. Based on the severity of the erosion and on the conservation state of the check dams, the stream was divided into four reaches (R1-R4, Figure 1): Their features are described in Table 1. Along with R5, instead, check dams are intact, and vegetation is denser and well stratified, with no signs of riverbed erosion.

Field Data and Satellite Imagery
LAI values were measured in the field from May to July 2015, during the growing season and before crop harvest. A large LAI dataset was collected at each study site using the LAI 2200 Plant Canopy Analyzer [44], a sensor composed of five concentric rings. Readings were collected along several cross-sections of the monitored watercourses, for a width ranging from 10 m to 20 m, in one sensor mode for above and below readings, using a 90 • view cap to avoid interference caused by the user's shadow. One above-canopy reading was detected for every four below-canopy readings. Moreover, to take into account the slope of the study area the "measuring on slope" protocol was applied [44]. Other optical corrections were not applied, since the effect on the measurement is negligible for areas with a slope less than 30 • [49]. The transects were conducted from one bank to the other at the site of a check dam and between two consecutive check dams. In the wheat fields and the forest area, readings were taken along a single direction, to obtain the most suitable overlap between field measurements and the satellite imagery acquired by the different sensors. A total of 82 field datasets were collected: 53 in CR, 11 in R5, 8 in DF, and 10 in WW. The acronym RV (Riparian Vegetation) is used henceforth to refer to CR and R5 together. A Garmin GPSMAP 76CSx was used to register the location of transects and check dam positions.
Multispectral data for this study is from medium-to high-resolution imagery provided by Landsat 8 (L8), Sentinel-2 (S-2), and Pleiades 1A (P1A) satellites. Considering the aims of this study, images with coarser resolutions were not used because larger pixel sizes could not detect riparian vegetation features adequately. The adopted L8 (30 m) and S-2 (10 m) medium resolution platforms are open access and guarantee excellent spatial and temporal coverage (Landsat from 1972, Sentinel from 2015). For these reasons, this satellite imagery is commonly used to carry out Vis studies and is of particular interest to both technicians and scientists [50]. In addition, P1A (2 m) high-resolution imagery was also used because it yields accurate Vis in forest environments [51][52][53].
The reduction of between scenes variability, due to different sensors acquisitions, was achieved by converting the imagery spectral digital number to top of atmosphere (TOA) reflectance. The advantages of using TOA reflectance arise from the necessity of better comparing remote sensing data acquired by different sensors at different times [48]. Also, TOA reflectance compensates different radiance values caused by different spectral bands, and lastly, corrects the variation in the Earth-Sun distance over the acquisition dates [54]. As suitable atmospheric aerosol data were not available for every image acquisition, all the sets of image data were not atmospherically corrected to top of canopy reflectance values.
Landsat data (https://landsat.usgs.gov/landsat-8) has been frequently used for land cover studies [55], but rarely for the analysis of riparian vegetation in streams with narrow cross-sections [27]. The cloud-free imagery used in this work, acquired by Landsat OLI sensors on 1 June and 3 July 2015, was downloaded using the Earth Explorer Get Data tool of the U.S. Geological Survey Landsat web page (http://earthexplorer.usgs.gov). L8 imagery data were corrected to TOA spectral radiance, using the radiance rescaling factors in the MTL file, and then converted to TOA reflectance, using the rescaling coefficients and the correction for the sun angle (https: //www.usgs.gov/land-resources/nli/landsat/using-usgs-landsat-level-1-data-product).
The Sentinel satellite is an orbital wide-swath (290 km) imaging mission carrying an innovative Multispectral Imager (MSI), which provides a versatile set of 13 spectral bands from the visible to the shortwave infrared at 10 to 60 m spatial resolution. For the present work, a cloud-free Level 1C (Top-of-Atmosphere reflectance) image acquired on 28 July 2015 was downloaded from the ESA Sentinel-2 Pre-Operations Hub (https://scihub.copernicus.eu). The Level 1C of S2 imagery provided data already in TOA reflectance along with the parameters to eventually transform them into radiances (https://sentinel.esa.int/web/sentinel/user-guides/sentinel-2-msi/product-types/level-1c). The SNAP toolbox was used to process the visible, the near-infrared, and the shortwave infrared bands (10 m, 10 m, and 20 m spatial resolution, respectively) to produce a subset covering 3.7 km × 1.94 km centered at 41 • 12 35 N, 15 • 23 05 E and to resample all the selected bands to a 10 m pixel size.
The Pleiades constellation provides multispectral data at 2.8 m pixel resolution for the visible and the near-infrared bands, at 0.7 m pixel resolution for the panchromatic band. A satellite operator resampled the images to 2 m and 0.5 m, respectively [56]. The satellite imagery used in this work was acquired by the Pleiades devices on 20 July 2015.
P1-A imagery data were first transformed to exoatmospheric (TOA) radiance using the calibration coefficients provided in the image metadata. TOA reflectance was then computed in each band considering the sun-earth distance in astronomical units, the extraterrestrial solar irradiance, and the solar zenith angle at the center of the image [56].
Considering the smooth slope of the study area and that NDVI is less sensitive to the topographical features [57,58], the topographic correction was not applied to all of the three satellite imagery.

LAI Calculation
LAI values were calculated from NDVI data retrieved from the three selected satellite images using the following Equation (1) [59]: In the case of L8 imagery, the corrected Normalized Difference Vegetation Index (NDVIc) was calculated to improve the correlation accuracy between the NDVI and LAI. The NDVIc uses the mid-infrared band as a filter to eliminate or reduce the contribution of ground vegetation and yields more accurate canopy values [60]. Although the correction was also tested for S-2, it was not applied because it yielded very low NDVI, probably due to differences in band wavelengths with respect to L8 [61]. Moreover, the correction was not performed for P1A data, due to the lack of information in the shortwave infrared band [62].

Equation Algorithm Source
Caraux-Garson (CG) Lacaze et al. [64] (b) Campbell and Norman (CN) Walthall et al. [65] (c) NDVI can (canopy) and NDVI back (background) are the asymptotic value of NDVI corresponding, respectively, to the higher LAI values and to bare soil; k is the extinction coefficient. In this work NDVI can = 0.859, NDVI back = 0.224, k = 0.213. fc represents the fraction crop endmember.

Results
The observed LAIs (LAI obs ) differed significantly at the monitored sites, both in range and median value ( Figure 3). Overall, LAI obs values along RV transects ranged from 0.26 to 5.93. In particular, CR showed a variation of 5.45 (0.26-5.71 range) and a median value of 2.97. The lower LAI values were mostly found in R2 and R3, the higher ones in R1, R2, and R4. In contrast, more uniform LAI values were detected in R5, where the variation was only 1.85 (4.08-5.93 range), and the median was 4.46.
Remote Sens. 2020, 12, x FOR PEER REVIEW 7 of 16 . Although this, the median, was in absolute the highest detected (4.77), which was very similar to the median obtained in R5. Moreover, the higher DF values were comparable to the higher LAIobs in R1, R2, and R3.
The PBIAS (percent bias), R 2 (coefficient of determination), and NSE (Nash-Sutcliffe efficiency index) statistical indices were used to assess the accuracy of LAI values calculated based on NDVI values derived from the different remote sensing images (Landsat 8, Sentinel-2, and Pleiades 1A) using Equations (a)-(c). This analysis was carried out considering the different types of vegetation monitored (riparian vegetation, deciduous forest, and winter wheat fields). In particular, to determine which combination of image and equation is most suited to the study of variability in riparian vegetation, the performance was evaluated for every single reach (R1-R5) of RV. The results are reported in Table 3.  The higher LAI values detected in R5 were comparable to those obtained in R1, R2, and R3 of the CR ( Figure 3).
As expected, WW showed the narrowest range of LAI obs (1.11-2.50) and the lowest median value (1.81). However, values detected in WW are similar to the lower values measured in R2 and R3. The variability of LAI obs values in DF was 4.26 (1.18-5.44 range). Although this, the median, was in absolute the highest detected (4.77), which was very similar to the median obtained in R5. Moreover, the higher DF values were comparable to the higher LAI obs in R1, R2, and R3.
The PBIAS (percent bias), R 2 (coefficient of determination), and NSE (Nash-Sutcliffe efficiency index) statistical indices were used to assess the accuracy of LAI values calculated based on NDVI values derived from the different remote sensing images (Landsat 8, Sentinel-2, and Pleiades 1A) using Equations (a)-(c). This analysis was carried out considering the different types of vegetation monitored (riparian vegetation, deciduous forest, and winter wheat fields). In particular, to determine which combination of image and equation is most suited to the study of variability in riparian vegetation, the performance was evaluated for every single reach (R1-R5) of RV. The results are reported in Table 3.  On the whole, all three satellite images underestimated LAI obs (PBIAS > 0) in R1. The equation which yielded the best accuracy was LB for L8 and P1A and CG for S-2. In particular, L8 and P1A performed comparably well, with a slightly finer accuracy for P1A (lower RMSE) and better than S-2. The extreme variability of LAIs measured in R2 (the incipient erosion reach) was not detected correctly by any imagery and was generally underestimated. In any case, LB is the equation that performed best for all the satellite platforms. L8, followed by P1A and S-2, when processed with equation CN, yielded the best results in R3. In particular, LAI obs values were slightly overestimated by L8 (PBIAS < 0) and underestimated by P1A and S-2. L8 and equation LB, gave the best results in R4. In this case, the calculated values slightly underestimated the LAI obs . In R5, the last reach of RV, P1A performed best, with slight overestimation. S-2 and L8, instead, tended to underestimate the LAI obs . Here, equation LB yielded the best results for all the satellite platforms. Considering the other two vegetation types, P1A showed the best statistical indices, followed by L8 and S-2 in DF, when all processed with LB. In general, all the images overestimated LAI obs values. In WW, both P1A and S-2 performed poorly, greatly underestimating the LAI obs . L8 performed better statistically, although LAI values were slightly overestimated. Equation CG worked best. Figure 4 represents the scatterplots between LAI obs and LAI retrieved values for each satellite imagery, which give a visual assessment of the validation results reported in Table 3. The plots confirm a good prediction ability for the LB equation, except for WW in P1A and S-2, while showing a tendency to underestimate LAI values for the CN and CG equations. performed poorly, greatly underestimating the LAIobs. L8 performed better statistically, although LAI values were slightly overestimated. Equation CG worked best. Figure 4 represents the scatterplots between LAIobs and LAI retrieved values for each satellite imagery, which give a visual assessment of the validation results reported in Table 3. The plots confirm a good prediction ability for the LB equation, except for WW in P1A and S-2, while showing a tendency to underestimate LAI values for the CN and CG equations.  For comparison. The ranges and median values of both LAI obs and LAI derived from all three satellite images (Landsat 8, Sentinel-2, and Pleiades 1A) were plotted using the equation that best suited each type of vegetation (Table 3), (Figures 5 and 6).
Remote Sens. 2020, 12, x FOR PEER REVIEW 9 of 16 For comparison. The ranges and median values of both LAIobs and LAI derived from all three satellite images (Landsat 8, Sentinel-2, and Pleiades 1A) were plotted using the equation that best suited each type of vegetation (Table 3), (Figures 5 and 6).     LAI obs values generally showed wider ranges than LAI retrieved from satellite imagery, especially in the reaches of the CR. None of the three satellite images detected the high variability of the LAI obs in R1, R2, and R3. In R4, S-2 yielded a value range very similar to that of LAI obs; in particular, the highest LAI obs and calculated LAI values were 5.60 and 5.68, respectively, whereas the lowest were 2.29 and 2.70, respectively ( Figure 5). S-2 also detected ranges well in R5, although with a tendency to underestimate both higher and lower LAI values ( Figure 6). In this reach, P1A yielded a more accurate minimum (4.08, LAI obs and 4.27, LAI retrieved) and maximum (5.93, LAI obs and 5.68, LAI retrieved) values ( Figure 6). P1A also accurately detected LAI variability in DF. The ranges were quite similar (4.26 for LAI obs and 3.75 for retrieved LAI), although retrieved LAI values slightly overestimated both the lower and higher values ( Figure 6). WW was the only land use category in which L8 performed best in terms of both lower and higher values ( Figure 6). However, the maximum value was underestimated (2.50 for LAI obs and 2.37 for retrieved LAI). P1A and L8 median values (3.54 and 3.39, respectively) were close to the observed one (3.73) in R1 ( Figure 5). In R2, instead, none of the satellite platforms considered yielded median LAI values in good agreement with the median LAI obs ( Figure 5). In R3, both L8 and P1A yielded median LAI values (1.41 and 1.43, respectively) similar to the LAI obs (1.32) ( Figure 5). S-2 also yielded accurate results, although the median LAI retrieved was underestimated (1.06). In R4, the median LAI obs was 3.13. L8 yielded the best results (3.10), whereas P1A and S-2 tended to overestimate actual values (3.37 and 3.50, respectively) ( Figure 5). In R5, where the median LAI obs was 4.46, P1A yielded the most accurate results (4.76), followed by L8 (4.12) and S-2 (4.86). As for the two other types of vegetation cover, DF and WW, only L8 yielded a median value comparable to the observed value. In DF, LAI obs was 4.77, and the retrieved LAI was 5.14 ( Figure 6); whereas, in WW, the LAI obs was 1.81, and the retrieved LAI was 1.97 ( Figure 6).

Discussion
The LAI of RV was measured in-situ (LAI 2200 Plant Canopy Analyzer) and retrieved from satellite imagery to test the performance of three satellite images with different spatial resolutions (Pleiades 1-A, Sentinel-2, and Landsat 8). DF and WW land cover types were also analyzed to compare results with those obtained for RV. As the NDVI is one of the Vis most commonly used to detect canopy cover variability [31], the LAI-NDVI approach was used to derive LAI values using three different equations (Lambert-Beer, Caraux-Garson, and Campbell-Norman).
The wide variability of in-field LAI obs measurements for RV reflects the composition and structure of such complex environments. In particular, in R2 and in R3, characterized, respectively, by partially or totally destroyed check dams, LAI obs values were lower, highlighting a correlation between the vegetation status and geomorphological channel stability [12,38,69]. R1 and R4, where there is almost no erosion and vegetation is more stratified, showed a higher maximum LAI obs comparable to the values measured in the DF. R5, characterized by intact check dams and well-structured riparian vegetation, has median LAI obs values similar to those detected in R4. Higher values are detected in DF because the vegetation here is mostly single story, without the presence of ground vegetation. As expected, LAI obs values for the WW field are low and show less variability, due to the unstratified cover and the simple leaf structure.
As for the performance of the three satellite images in the study area, L8 and P1A showed similar results in R1, R3, R4, R5, and DF. Where erosion processes influence the presence of vegetation (R3), P1A and L8 yielded very similar ranges and median values. In R1 and R4, instead, although statistical values were similar between L8 and P1A, there was a better match between value ranges for P1A.
In R5 and DF, the better statistical performance of P1A was due to the high-resolution, which limited the negative influence of "mixed pixels" [70,71], which instead affected L8 [27] (Figure 6), as well as S-2, whose image resolution (10 m) was too coarse to avoid the "mixed pixels" effect. In R4 and WW, instead, L8 performed better than P1A. In Figures 5 and 6, LAI obs and retrieved LAI show a better match in terms of median values for both R4 and WW, and the best match of value ranges for WW. In both cases, P1A underestimated LAI, especially the lower range values.
The reason why P1A show no noticeable improvements with respect to L8, especially considering R4 and WW, may be due to the different satellite data acquisition period (1 June and 3 July for L8, 20 July for P1A). Satellite imagery are sensible to biological factors [72], such as the phenological phases of vegetation [73], and imagery acquired at the right time for the purposes of a study can yield better results, even with images at a coarser scale [74]. At the end of July, the herbaceous plants characterizing R4 riverbed vegetation were in a water stress phase, and all crops had been harvested in the WW fields. Moreover, although the higher spatial resolution, P1A is characterized by a lower spectral resolution with respect to L8 and S-2 (4 bands in P1A) [53]. This means that the P1A red and NIR bandwidth are greater respect to the same bands of L8 and S-2. Teillet et al. [75] demonstrated that an increase in the bandwidth is directly related to a decrease in the values of the NDVI.
In general, the statistical performances of S-2 were worse than those of the other sensors, except for R5, where it performed comparably to L8 and P1A, while it was the only image to capture the extreme variability of LAI obs values in CR and R5 (Figures 5 and 6). The worse statistical performances of S-2 could depend on several factors. First, as mentioned above, the acquisition date (28 July) certainly influences the accuracy of results [74]. Moreover, several authors [76][77][78] highlight that although NDVI is commonly used to retrieve LAI values, it is not always the VI with the highest correlation. In particular, for S-2, the newest of the three satellites, alternative Vis may be considered to improve results [32,[79][80][81].
None of the sensors was a good descriptor of the high variability of LAI obs values in R2, resulting in a general underestimation. In this reach, where the vegetation was sparse or totally absent [27], L8 and S-2 were affected by the "mixed pixel" effect [71]. P1A, which is characterized by a lower spectral resolution, and consequently, by higher bandwidth, suffered a spectral interference of the soil background, due to the surface heterogeneity [82,83].
As regards the use of different equations to derive LAI from the NDVI, results indicate that, in almost all cases, the Lambert-Beer (LB) equation was the most suitable. In particular, when the vegetation was more homogenous (R1, R4, R5, and DF), LB worked best in processing both the L8 and the P1A imagery. Moreover, LB (coupled with P1A) improved the results in R2, although with unsatisfactory statistics, confirming that the more complex features in R2 are not easily detected using satellite imagery (Romano et al., 2020). However, it should be considered that at each specific site, there may be variations related to canopy characteristics (i.e., density, leaf angle distribution) and soil optical properties [84,85]. For this reason, the match between the LAI obs and retrieved LAI could still be improved by modifying the De Jong coefficients [86], generic for Mediterranean areas, used in the LB equation. In R3, confirming findings by Ricci et al. [27] for L8, the equation CN performed well also in processing S-2 and P1A. This equation requires identifying two NDVI values representing the value of bare soil and the maximum value for vegetation [65]. In a complex environment, such as a Mediterranean watershed that includes different types of vegetation, the first value is easy to determine, whereas the second is difficult to determine because the maximum value of the vegetation depends, within the same image, on changes in vegetation type [87,88]. For this reason, CN yielded acceptable results in R3 only, where only LAI obs values similar to those of the bare soil were estimated correctly ( Figure 5). Linear equation CG yielded the best accuracy for WW because, where the observed LAI values were lower than 3, results were not affected by the effect of saturation on greenness, which for higher values determines an asymptotic trend in retrieved LAIs [31,85]. Generally, differences in results could also be ascribed to the fact that the adopted equations were developed before the launch of both the P1A and the S-2 missions, so that the correlation coefficients may not be optimal. In addition, since LAI obs measurements are required to validate retrieved data, the sampling design is an important aspect to consider [89]. In this work, transects were selected to obtain the most suitable overlap between field measurements and the spatial resolutions of satellite images. For this reason, P1A results could be improved using a different sampling design for medium resolution imagery (L8 and S-2) and high-resolution imagery (P1A). The obtained results, confirm what was observed by Wang et al. [53], or that the spatial resolution is an important factor to consider when choosing a satellite product, but to the differences in the satellite sensors bands (spectral resolution) could also impact the results.

Conclusions
This study assessed the potential of three satellite images with different resolutions (Landsat 8, Sentinel-2, and Pleiades 1A) in studying vegetation features in a mountain riparian environment of Southern Italy. Pleiades 1A and Landsat 8 imagery yielded values closer to the observed ones. In particular, the higher spatial resolution of the former limited the influence of the mixed pixel over the study area and improved the results both in terms of statistics and retrieved LAI value ranges. However, the lower spectral resolution of the Pleiades 1A leads to spectral interference, due to the surface heterogeneity, and to an underestimation of the results, especially in the transects where the vegetation is variable. Sentinel-2 imagery showed poor statistical performances attributable, mainly, to both the acquisition date of the image and to the vegetation index considered, which may not be the best for this type of sensor. The results obtained suggest that, apart from the spatial resolution, other factors have an impact on the LAI retrieval. The spectral resolution could negatively influence the results, especially when the environment is more complex. The time of acquisition should be as close as possible to the in-field measurement period, since the satellite multispectral images are sensitive to the phenology of vegetation. The sampling design could improve the accuracy of the results if it is specifically designed for the satellite sensor. A cost analysis is another key issue, as some remote sensing images can be downloaded free of charge (i.e., Landsat 8 and Sentinel 2), while others cannot (i.e., Pleiades 1A). Future studies could aim to improve the coefficient of the Lambert-Beer equation for the Mediterranean areas, and/or test the performance of different vegetation indices-especially for the Sentinel-2 imagery. Funding: The Pleiades acquisition was founded by Fondazione Cassa di Risparmio di Puglia, Project "Impiego delle immagini da satellite nella gestione della risorsa idrica a scala di bacino", Scientific Coordinator prof. Francesco Gentile, University of Bari.