Use of Sentinel-2 Satellite for Spatially Variable Rate Fertiliser Management in a Sicilian Vineyard

: Satellites can be used for producing maps of within-ﬁeld crop and soil parameters and, consequentially, spatially variable rate crop input application maps. The plant vegetative vigour index (i.e., Normalised Difference Vegetation Index—NDVI) and the leaf water content index (i.e., Normalised Difference Water Index—NDWI) maps were used to study—through both time and space—the phenological phases of two plots, with Syrah and Nero d’Avola grapevine varieties, in a Sicilian vineyard farm, located in Naro (Agrigento, Sicily, Italy). The aim of this work is to produce spatially variable rate nitrogen fertiliser maps to be applied in the two vineyard plots under study as well as to understand when they should be fertilised or not according to their target crop yields. The average plant vegetative vigour and leaf water content of both the plots showed a high temporal and spatial variability during all phenological phases and, according to these results, the optimal fertilisation time should have been 12 April 2021. In fact, this crop operation is aimed at supporting the vegetative activity but must be performed when the soil water and, therefore, the plant leaf water content are high. Therefore, spatially variable rate fertilisation should have been performed around 12 April 2021 in both plots, using previous NDVI maps and taking into consideration two management zones. This work demonstrates the usefulness of remote sensing data as Decision Support Systems (DSS) for nitrogen fertilisation in order to reduce the production cost, environmental impact and climate footprints per kg of produced grapes, according to the European Green Deal challenges.


Introduction
Food production is paramount for everybody, and producing food in a cost-effective way is the goal of any farmer, large-scale farm manager and regional agricultural agency. A farmer needs to be informed in order to be efficient; he must know information about his products in order to forge a viable strategy for crop operations. These tools will help him understand the soil parameters, potential yield, crop health and the extent of eventual diseases and/or stress damages. Commodity brokers are also very interested in how well farmers are producing, as crop yield (both quantity and quality) estimates all products' control price and worldwide trade.
Moreover, spatial variability in vineyards is an important issue that affects both farmers and winemakers, especially when trying to obtain the best potential from their vines. Land (topography and soil), climate and diseases affect crop development and, therefore, grape yield and quality in a vineyard. Some of these factors of spatial variation remain stable over time. The study of these factors can lead winemakers to better understand the images remotely sensed by this satellite must follow for developing a Decision Support System (DSS). The final aim is to produce spatially variable rate crop input application maps which must be adjusted in every crop season in order to become increasingly fitting to the crop needs.
Among several vegetation indexes that can be calculated, NDVI (Normalised Difference Vegetation Index) and NDWI (Normalised Difference Water Index) have some practical advantages. In fact, NDVI can be applied to evaluate the plant photosynthetic activity and, therefore, its vegetative vigour and soil coverage during growth phases, so that it is measured by optical sensors (e.g., Arvatec OptRx and Trimble GreenSeeker) for planning fertilisation [12,[32][33][34][35][36][37] and irrigation. Therefore, NDVI mapping remains the most affordable way to easily and objectively obtain vineyard spatial information [3], as this index has proved to be a useful tool for monitoring crop yield [18] and table grape quality characteristics [14]. Moreover, NDWI can be applied to evaluate plant water stress and leaf water content during growth phases, so that irrigation can be monitored in near real-time in order to improve agriculture efficiency, especially in areas where meeting the plant need for water is more challenging [38].
The aim of this work is to produce spatially variable rate nitrogen fertiliser maps to be applied in two vineyard plots under study, as well as to understand when they should be fertilised or not according to their target crop yields in order to provide the highest economic and environmental benefits.

Study Area
During the first visit to the Sicilian vineyard farm "Chilometro7", located at a Longitude of 13 • 53 34 , Latitude of 37 • 18 26 and Altitude of 634 m. a.s.l. in the territory of Naro (Agrigento, Sicily, Italy), the main vineyard (Vitis vinifera L.) characteristics (e.g., vine varieties, fruit tree forms and plant distances) and the performed crop operations (e.g., fertilisation, green pruning and harvest) of two plots, less than 100 m far from each other, were monitored and noted: (1) plot having an area of 1.75 ha ca., where the vine variety Syrah is cultivated, using the fruit tree form step-over espalier with plant distances of 2.5 × 1 m ( Figure 1); (2) plot having an area of 1 ha ca., where the vine variety Nero d'Avola is cultivated, using the fruit tree form marquee with plant distances of 2.8 × 2.8 m ( Figure 2).

Satellite Image Sensing and Statistical Data Analysis
It was possible to study and implement a protocol for converting the images remote sensed by the satellite Copernicus Sentinel-2 into maps of within-field crop and soil p rameters and then spatially variable rate crop input application maps.
It was first necessary to select the vegetation indexes to be remotely sensed by t satellite Sentinel-2, i.e., NDVI and NDWI, in order to produce maps of each vegetati index for each plot [27].
In fact, the plant vegetative vigour is expressed by NDVI, i.e., Normalised Differen Vegetation Index, which is defined as: where: -ρRED is the radiance (in reflectance units) of a red channel near 0.66 µm; and -ρNIR is the radiance (in reflectance units) of a near-IR channel around 0.86 µm [39 As NDVI relates to the reflectance in the red and near-infrared wavelengths, it is co sidered an indicator of vegetation growth and green biomass accumulation [36]. Differe studies using optical remote sensing demonstrated the relationship between canopy flectance and biomass production [37]. In fact, Johnson et al. [40] showed a significa correlation (R 2 = 0.74) between NDVI and LAI (Leaf Area Index) in vines, even if this ha pens for LAI values lower than 3.  Table 1 shows the physical and chemical soil parameters of Syrah and Nero d'Avola plots. According to the current spatially uniform rate vineyard management, fertilisation is carried out in Syrah plot during winter, e.g., between 14 and 24 February 2021, by applying the rate of 1143 kg ha −1 of Fomet powder organic fertiliser (constituted by cattle and poultry manure converted into humus) ( Table 2), followed by no irrigation. No green pruning was carried out in the Syrah plot.
Fertilisation is also carried out in the Nero d'Avola plot during winter, e.g., between 14 and 24 February 2021, by applying the rate of 1000 kg ha −1 of the above organic fertiliser, followed by no irrigation. Green pruning was carried out only in the Nero d'Avola plot in June.
The grapes of Syrah and Nero d'Avola are harvested in different periods (e.g., 9 and 21 September 2021, respectively) and sold to the wine-making industry "Cummo Vini-Casa Vinicola Sicania", located in the territory of Canicattì (Agrigento, Sicily, Italy), for the production of two different brands of wine-"Desiderio" and "Nero d'Avola Bio", respectively.

Satellite Image Sensing and Statistical Data Analysis
It was possible to study and implement a protocol for converting the images remotely sensed by the satellite Copernicus Sentinel-2 into maps of within-field crop and soil parameters and then spatially variable rate crop input application maps.
It was first necessary to select the vegetation indexes to be remotely sensed by the satellite Sentinel-2, i.e., NDVI and NDWI, in order to produce maps of each vegetation index for each plot [27].
In fact, the plant vegetative vigour is expressed by NDVI, i.e., Normalised Difference Vegetation Index, which is defined as: where: -ρRED is the radiance (in reflectance units) of a red channel near 0.66 µm; and -ρNIR is the radiance (in reflectance units) of a near-IR channel around 0.86 µm [39].
As NDVI relates to the reflectance in the red and near-infrared wavelengths, it is considered an indicator of vegetation growth and green biomass accumulation [36]. Different studies using optical remote sensing demonstrated the relationship between canopy reflectance and biomass production [37]. In fact, Johnson et al. [40] showed a significant correlation (R 2 = 0.74) between NDVI and LAI (Leaf Area Index) in vines, even if this happens for LAI values lower than 3.
In order to sense the NDVI related to vineyard farms, it is necessary to retrieve the satellite images taken during the phenological phases when the value of this vegetation index is supposed to change [36]: before green pruning (May); 2.
In fact, the plant vegetative vigour and the plant leaf water content maps of the two surveyed plots were studied during the above four phenological phases. It was possible to retrieve the satellite images including the two plots where the vine plants of Syrah cultivar and those of Nero d'Avola variety are cultivated, respectively.
Bottom-of-Atmosphere (BOA) reflectance values of Sentinel-2 data were obtained from ESA (European Space Agency) servers, taking into consideration Copernicus program. The image spatial resolution was 10 m for B4 and B8 bands, while it was 20 m for B12 band, as well as the temporal resolution was five days. The obtained images were atmospherically corrected by means of the Sen2Cor algorithm, while the ArcGIS 10.x (ESRI, Redlands, CA, USA) software was used.
Image statistical data (e.g., minimum, maximum, average, standard deviation, Coefficient of Variation of image pixels distribution) of both NDVI ("Plant") and NDWI ("Water"), from the beginning of April to the end of September 2021 were calculated. Therefore, the graphs of both indexes were plotted for the above period and each plot was divided into two management zones according to the calculated statistical data.

Results
The results of this work, expressed by graphs and maps, based on Sentinel-2 data, are structured as follows: vegetative vigour; -leaf water content; -criterion for deciding the optimal fertilisation time; -temporal variability; and spatial variability. Table 3 and Figure 3 show the meteorological data in the surveyed area during the considered period (Sicilian Agrometeorological Information System, Sicilian Region-Department of Agricultural and Food Resources-Section Infrastructural Interventions). The total rainfall from 01 April 2021 to 31 October 2021 was 254.8 mm.

Vegetative Vigour
In the graph of average NDVI (Figure 4), it is possible to observe a similar trend for both plots (i.e., vine cultivars). In fact, the plant vegetative vigour highly increases from April to June and decreases at the beginning of June, more so in the Nero d'Avola cultivar rather than the Syrah one, as a consequence of green pruning, carried out in this month only in Nero d'Avola plot-where it caused a reduction of vegetative vigour. Then, the vegetative vigour becomes almost constant between mid-June and the beginning of July and decreases in July, as the productive activity (i.e., grape ripening) is prevalent on the vegetative one and the high temperature causes plant water stress. Finally, the vigour increases in August, as a slightly lower temperature causes lower plant water stress. It suddenly decreases in September, because the rest period (i.e., plant senescence or low vegetative activity) starts before for Syrah cultivar (harvested on 9 September) and after for the Nero d'Avola variety (harvested on 21 September).

Vegetative Vigour
In the graph of average NDVI (Figure 4), it is possible to observe a similar trend f both plots (i.e., vine cultivars). In fact, the plant vegetative vigour highly increases fro April to June and decreases at the beginning of June, more so in the Nero d'Avola cultiv rather than the Syrah one, as a consequence of green pruning, carried out in this mon only in Nero d'Avola plot-where it caused a reduction of vegetative vigour. Then, t vegetative vigour becomes almost constant between mid-June and the beginning of Ju and decreases in July, as the productive activity (i.e., grape ripening) is prevalent on t vegetative one and the high temperature causes plant water stress. Finally, the vigour creases in August, as a slightly lower temperature causes lower plant water stress. It su denly decreases in September, because the rest period (i.e., plant senescence or low veg tative activity) starts before for Syrah cultivar (harvested on 9 September) and after for t Nero d'Avola variety (harvested on 21 September).
However, the Nero d'Avola plot shows slightly higher values of NDVI rather th the Syrah one, from the mid of May to the end of September. The soil is similar in the tw plots, except for organic matter content, which is 2.71 and 1.88% in the Syrah and Ne d'Avola plots, respectively, and nitrogen content, which is 1.25 and 1.10 g kg −1 in the Syr and Nero d'Avola plots, respectively (Tables 1 and 2). The difference in NDVI is due the genetic characteristics of the cultivars and the fruit tree form. In fact, Nero d'Avo cultivar has a genetic attitude of higher vegetative vigour than the Syrah variety, as w as a marquee tree form that promotes higher vegetative vigour in the former plot rath However, the Nero d'Avola plot shows slightly higher values of NDVI rather than the Syrah one, from the mid of May to the end of September. The soil is similar in the two plots, except for organic matter content, which is 2.71 and 1.88% in the Syrah and Nero d'Avola plots, respectively, and nitrogen content, which is 1.25 and 1.10 g kg −1 in the Syrah and Nero d'Avola plots, respectively (Tables 1 and 2). The difference in NDVI is due to the genetic characteristics of the cultivars and the fruit tree form. In fact, Nero d'Avola cultivar has a genetic attitude of higher vegetative vigour than the Syrah variety, as well as a marquee tree form that promotes higher vegetative vigour in the former plot rather than step-over espalier form in the latter one.
As the main objective of the considered vineyard farm is to maximise the crop yield (which, in 2021, resulted in 7.4 t ha −1 ca. in the Syrah plot and 11 t ha −1 ca. in the Nero d'Avola plot), it is necessary to increase the crop yield of Syrah cultivar to levels similar to those of Nero d'Avola by implementing fertilisation in the right time and with the right rate; i.e., the rate required for satisfying the need of the species for nutrients such as nitrogen. In fact, if fertilisation is performed at the right time-not too late, when the soil water content is low-it is possible to prevent root burning and minimise the emissions of greenhouse gases (GHGs), e.g., nitrous oxides (NO x ) from nitrate (NO 3 ). Moreover, if fertilisation is carried out at the right rate, it is possible to provide the vine plants with the nutrient amounts needed for vegetative activity well-balanced with productive activity. Thus, green pruning, which is aimed at reducing vegetative activity and, therefore, promoting productive activity, can be avoided in both plots in order to reduce production costs. tion is carried out at the right rate, it is possible to provide the vine plants with the nutrient amounts needed for vegetative activity well-balanced with productive activity. Thus, green pruning, which is aimed at reducing vegetative activity and, therefore, promoting productive activity, can be avoided in both plots in order to reduce production costs.

Leaf Water Content
In the graph of average NDWI ( Figure 5), it is possible to observe an almost identical trend for both plots (i.e., cultivars), even if Nero d'Avola cultivar has a higher plant leaf water content rather than Syrah during most of the surveyed period, meaning that the former has lower water stress rather than the latter. In fact, the plant leaf water content increases in the first half of April, due to the last rains of the season, and decreases in the second half of this month. Then, the leaf water content increases, decreases and increases in May-as a result of the rain distribution in this month-becomes almost constant in June, slightly decreases in July and increases in August. Finally, the leaf water content decreases in the first half of September and then suddenly increases because, even if the rest period starts, vine plants keep their leaves during winter but low temperature causes low transpiration.

Leaf Water Content
In the graph of average NDWI ( Figure 5), it is possible to observe an almost identical trend for both plots (i.e., cultivars), even if Nero d'Avola cultivar has a higher plant leaf water content rather than Syrah during most of the surveyed period, meaning that the former has lower water stress rather than the latter. In fact, the plant leaf water content increases in the first half of April, due to the last rains of the season, and decreases in the second half of this month. Then, the leaf water content increases, decreases and increases in May-as a result of the rain distribution in this month-becomes almost constant in June, slightly decreases in July and increases in August. Finally, the leaf water content decreases in the first half of September and then suddenly increases because, even if the rest period starts, vine plants keep their leaves during winter but low temperature causes low transpiration.

Criterion for Deciding the Optimal Fertilisation Time
In the graphs of average NDVI and NDWI for the Syrah and Nero d'Avola plots ( Figure 6), it is possible to observe that vegetative vigour suddenly increases from 12 May 2021, while the leaf water content suddenly increases until 12 April 2021. Thus, as fertilisation is aimed at supporting vegetative activity and must be performed when the soil water content and consequently plant leaf water content is high, the optimal time for this crop operation should have been 12 April 2021, above all because no irrigation was performed. In fact, in this case, a low soil water content would cause root burning and emissions of GHGs (i.e., nitrous oxides from nitrate). Therefore, the criterion for deciding the optimal fertilisation time is based on the comparison between NDVI and NDWI graphs for each plot. In fact, the optimal fertilisation time can be different if a farmer wants to promote vegetative vigour or fruit quality. However, in the case study, the objective is to produce the right amount of vineyard Leaf Area Content (LAC) that, for example, can maximise the crop yield. Even if different vine cultivars can have different needs for nutrients, they generally need higher amounts of available nutrients and water when plant LAC increases. In this period, if the plants cannot find the right amounts of nutrients and water in the soil to adsorb for growing, this is a factor limiting both vegetative and productive activities. Thus, plotting the curves of NDVI and NDWI for both cultivars ( Figure 6) allows one to define the optimal fertilisation time, which is normally before the highest plant growth rate is achieved and is fundamental for providing the plants with the best nutritional soil conditions, i.e., amount of both nutrients and water to be adsorbed by the plants.

Criterion for Deciding the Optimal Fertilisation Time
In the graphs of average NDVI and NDWI for the Syrah and Nero d'Avola plots (Figure 6), it is possible to observe that vegetative vigour suddenly increases from 12 May 2021, while the leaf water content suddenly increases until 12 April 2021. Thus, as fertilisation is aimed at supporting vegetative activity and must be performed when the soil water content and consequently plant leaf water content is high, the optimal time for this crop operation should have been 12 April 2021, above all because no irrigation was performed. In fact, in this case, a low soil water content would cause root burning and emissions of GHGs (i.e., nitrous oxides from nitrate). Therefore, the criterion for deciding the optimal fertilisation time is based on the comparison between NDVI and NDWI graphs for each plot. In fact, the optimal fertilisation time can be different if a farmer wants to promote vegetative vigour or fruit quality. However, in the case study, the objective is to produce the right amount of vineyard Leaf Area Content (LAC) that, for example, can maximise the crop yield. Even if different vine cultivars can have different needs for nutrients, they generally need higher amounts of available nutrients and water when plant

Temporal Variability
The temporal variability of the within-field crop and soil parameters can be assessed by means of the Coefficient of Variation (CV), which is a measure of relative variability, i.e., a standardised measure of the dispersion of a probability distribution or frequency distribution. It is defined as the ratio of the standard deviation (σ) to the average (µ), often expressed as a percentage: The temporal variability of vegetative vigour and leaf water content in both Syrah and Nero d'Avola plots was assessed based on the CV of NDVI and NDWI, respectively, from the beginning of April to the end of October 2021.
The CV of NDVI and NDWI of both the Syrah (Figure 7) and Nero D'Avola plots ( Figure 8) show a high temporal variability during the four surveyed phenological phases (before green pruning, after green pruning, before harvest, after harvest).

Temporal Variability
The temporal variability of the within-field crop and soil parameters can be assessed by means of the Coefficient of Variation (CV), which is a measure of relative variability, i.e., a standardised measure of the dispersion of a probability distribution or frequency distribution. It is defined as the ratio of the standard deviation (σ) to the average (µ), often expressed as a percentage: The temporal variability of vegetative vigour and leaf water content in both Syrah and Nero d'Avola plots was assessed based on the CV of NDVI and NDWI, respectively, from the beginning of April to the end of October 2021.
The CV of NDVI and NDWI of both the Syrah (Figure 7) and Nero D'Avola plots ( Figure 8) show a high temporal variability during the four surveyed phenological phases (before green pruning, after green pruning, before harvest, after harvest).
Moreover, the Syrah plot shows a higher CV, above-average NDWI, than the Nero d'Avola plot, meaning that the former has a higher variation of water stress than the latter.     Moreover, the Syrah plot shows a higher CV, above-average NDWI, than the Nero d'Avola plot, meaning that the former has a higher variation of water stress than the latter.

Spatial Variability
The spatial variation of the NDVI in a plot can be measured using the Coefficient of Variation (CV) [41]. The threshold is set as the median value of the coefficients of all tested vineyard fields [42] at fertilisation time, i.e., 12 April 2021. Thus, from 06 February 2021 to 07 April 2021, the CV of NDVI ranges from a minimum value of 4.949% to a maximum of 27.859% for the Syrah plot, while it ranges from a minimum value of 3.976% to a maximum one of 17.395% for the Nero d'Avola plot. In this case, the median value is 19.28% for the Syrah plot, while it is 10.322% for Nero d'Avola one.
Therefore, a spatially variable rate fertiliser application should have been planned in both the plots, when the CV of NDVI was higher than 19.28% for the Syrah plot, i.e., on 2 April 2021, and higher than 10.322% for the Nero d'Avola plot, i.e., on 2 April 2021 or 22 April 2021. Thus, according to the principles of precision agriculture, a higher rate should have been applied where the plant vegetative vigour was lower, while a lower rate should have been applied where the vigour was higher in order to have an uniform within-field vigour and, therefore, crop yield.
As fertilisation should have been carried out on 12 April 2021 and in a geo-referenced way, for both the plots, the map of the average NDVI of 7 April 2021-immediately before 12 April 2021-should have been used to produce the spatially variable rate fertiliser map. Indeed, if an organic fertiliser should have been used rather than a chemical one, as was done in 2021, the map of 12 April 2021 should have been used so that over five days, the nutrients contained in this crop input could have been converted in soluble forms to be adsorbed by the plant roots.
Then, based on the maps of the average vegetative vigour and leaf water content of 12 April 2021, the spatially variable rate fertilisation maps, including two management zones, were produced for the Syrah (Figure 9) and Nero d'Avola ( Figure 10) plots, by also considering the spatial and temporal image statistical data for each plot. In fact, some scientists, e.g., Bonilla et al. [3], discovered significant differences for the most part of the analysed traits of grape chemical composition between two areas; thus, they suggested two zones for practical differential management. Moreover, Allegro et al. [34] discovered significant differences in vegetative vigour between two areas; thus, they applied a fertiliser rate of 60 kg ha −1 in the low-vigour management zone and another of 30 kg ha −1 in the high-vigour zone.
way, for both the plots, the map of the average NDVI of 07 April 2021-immediately before 12 April 2021-should have been used to produce the spatially variable rate fertiliser map. Indeed, if an organic fertiliser should have been used rather than a chemical one, as was done in 2021, the map of 12 April 2021 should have been used so that over five days, the nutrients contained in this crop input could have been converted in soluble forms to be adsorbed by the plant roots.
Then, based on the maps of the average vegetative vigour and leaf water content of 12 April 2021, the spatially variable rate fertilisation maps, including two management zones, were produced for the Syrah (Figure 9) and Nero d'Avola ( Figure 10) plots, by also considering the spatial and temporal image statistical data for each plot. In fact, some scientists, e.g., Bonilla et al. [3], discovered significant differences for the most part of the analysed traits of grape chemical composition between two areas; thus, they suggested two zones for practical differential management. Moreover, Allegro et al. [34] discovered significant differences in vegetative vigour between two areas; thus, they applied a fertiliser rate of 60 kg ha −1 in the low-vigour management zone and another of 30 kg ha −1 in the high-vigour zone.  Moreover, the above management zones could be used for soil sampling in order to determine the soil parameters, both chemical (e.g., water and nutrients contents) and physical (e.g., texture), in each zone. Furthermore, the above vineyard management zones could allow to plan one of two harvest strategies, i.e., temporally variable harvest in the zones having different ripening periods, for the production of uniform quality grapes and therefore wine, or spatially variable harvest (by means of a grape harvester equipped with two hoppers) for the production of different quality grapes and therefore wines [1,2,43]. In fact, under a climate change scenario, it may be interesting to consider the NDVI zones for a selective harvest based on grape juice pH [3].
As another option, fertiliser rate can be associated with each range of the average NDVI shown in the map of 7 April 2021 (or 12 April 2021) in order to produce a spatially variable rate fertiliser map for each of the two surveyed plots by taking into account the best plant vegetative vigour/crop yield ratio. Moreover, the above management zones could be used for soil sampling in order determine the soil parameters, both chemical (e.g., water and nutrients contents) a physical (e.g., texture), in each zone. Furthermore, the above vineyard management zones could allow to plan one of tw harvest strategies, i.e., temporally variable harvest in the zones having different ripeni periods, for the production of uniform quality grapes and therefore wine, or spatially v iable harvest (by means of a grape harvester equipped with two hoppers) for the produ tion of different quality grapes and therefore wines [1,2,43]. In fact, under a climate chan scenario, it may be interesting to consider the NDVI zones for a selective harvest based grape juice pH [3].
As another option, fertiliser rate can be associated with each range of the avera NDVI shown in the map of 07 April 2021 (or 12 April 2021) in order to produce a spatia variable rate fertiliser map for each of the two surveyed plots by taking into account t best plant vegetative vigour/crop yield ratio.

Discussion
The Syrah variety has a much higher plant density (area of 1.75 ha, plant distances 2.5 × 1 m, plant density of 7000 plants ha −1 ) when compared with Nero d'Avola (area o ha, plant distances of 2.8 × 2.8 m, plant density of 1276 plants ha −1 ). By taking into accou this difference in plant density, it would be normal that the Syrah variety had a high NDVI during the whole year, as a consequence of the much higher plant density. Instea the opposite is shown in Figure 4. Thus, in a dry (not irrigated) vineyard, a higher pla density determines a higher competition for resources (i.e., water and nutrients) amo plants that, therefore, are less vigorous. However, a higher plant density (in the Syr

Discussion
The Syrah variety has a much higher plant density (area of 1.75 ha, plant distances of 2.5 × 1 m, plant density of 7000 plants ha −1 ) when compared with Nero d'Avola (area of 1 ha, plant distances of 2.8 × 2.8 m, plant density of 1276 plants ha −1 ). By taking into account this difference in plant density, it would be normal that the Syrah variety had a higher NDVI during the whole year, as a consequence of the much higher plant density. Instead, the opposite is shown in Figure 4. Thus, in a dry (not irrigated) vineyard, a higher plant density determines a higher competition for resources (i.e., water and nutrients) among plants that, therefore, are less vigorous. However, a higher plant density (in the Syrah plot) determines a lower grape production per plant (1.06 vs. 8.62 kg plant −1 ) and a lower crop yield (7.4 vs. 11 t ha −1 ), on average. Under the surveyed soil parameters and climate conditions, even if a higher plant density (in Syrah plot) decreases both the production per plant (kg plant −1 ) and crop yield (t ha −1 ), it reduces the vigour (expressed by a lower NDVI in Figure 4), so that no green pruning is needed and, finally, the production cost is reduced. Furthermore, as a lower vegetative vigour (in the Syrah plot) promotes the productive activity and both the production per plant and crop yield are lower, a higher quality of grapes and, therefore, wine, is expected.
As the Nero d'Avola cultivar has a higher plant leaf water content during most of the surveyed period ( Figure 5) and therefore, more water stress rather than the Syrah cultivar, higher plant density increases plant competition for resources, e.g., soil water. In fact, water stress plays a key role in vineyard management, above all in dry (not irrigated) plots such as those under study. Thus, less water per plant (in the Syrah plot, having 7000 plants ha −1 ) determines a lower vegetative vigour and root volume and, therefore, a lower soil nutrient adsorption, resulting in a lower NDVI.
The NDVI temporal evolution for the two surveyed cultivars differs from that found by other scientists. In fact, Vélez et al. [18], for the Verdejo cultivar, found a lower peak in July and a higher one in September, after harvest (instead of the highest peak found in June in this study). Moreover, Heemann Junges et al. [36] found, on average, a peak when the berries begin to swell in the second half of October (which, for climate conditions, is very similar to June in the boreal hemisphere), for the Chardonnay variety, and a peak when the berries touch each other in the first half of December (which, for climate conditions, is very similar to July-August in the boreal hemisphere), for the Cabernet Sauvignon variety. These differences express different benchmark curves for both cultivars, depending on genetic characteristics, and different climate conditions, depending on the year.
Very few scientists have also studied NDWI temporal evolution, such as Borgogno-Mondino et al. [25], who performed preliminary tests on table grape vineyards covered with polyethylene sheets, using Sentinel-2 free optical data, and compared their spectral response with that of similar uncovered vineyards in order to assess if a reliable spectral signal is detectable through the plastic cover. The maximum NDVI was obtained in the first half of September or in the first half of April, while the maximum NDWI was recorded in the first half of October or in the first half of April, according to the surveyed covered vineyard. Moreover, their results demonstrated that no significant limitation (for both NDVI and NDWI) was introduced by plastic sheets while monitoring the spectral behaviour of covered vineyards.
Based on the previous studies and this work, NDVI and NDWI can be monitoring tools of plant behaviour and health; they can be used for planning fertilisation, irrigation, water stress management, pest and disease control.
However, the interpretation of the maps remotely sensed from the satellite Sentinel-2 for DSS is the most difficult step of any precision viticulture platform for spatially variable rate vineyard management. The maps remotely sensed from the satellite Sentinel-2 must be compared to each other in order to find eventual correlations among different crops (e.g., yield, plant vegetative vigour, expressed by NDVI, and plant leaf water content, expressed by NDWI) and soil parameters. The final aim is the spatially variable rate application of crop inputs (e.g., fertiliser), in order to achieve more sustainable management, in terms of crop inputs and production cost optimisation, as well as in terms of environmental and climate footprints according to the European Green Deal objectives.
In fact, climate change has also affected vine physiology [19] and contributed to increasing the berry ripeness process. Therefore, wines present higher alcohol content and tend to have higher pH. Generally, for most of the winemaking regions, the warming trend from the last decades to date has led to higher quality wines (with higher concentrations in anthocyanins and polyphenols), especially in the Northern regions. However, the average increase of +2 • C-forecasted by 2049-should affect vine phenology, leading to high-temperature ripening processes [3].

Conclusions
This work demonstrates the usefulness of remote sensing as DSS for precision viticulture and traditional viticulture, i.e., for deciding if, when and how spatially variable rate fertiliser application must be performed in a field, according to the principles of precision agriculture. In fact, remote sensing can be used for: deciding if spatially variable rate fertiliser management must be performed in a field, based on the eventual within-field spatial variability of plant vegetative vigour; -identifying the optimal fertilisation time, based on the comparison between the graph of plant vegetative vigour and that of plant leaf water content (also for spatially uniform rate fertilisation, within traditional viticulture); and determining the spatially variable fertiliser rate to be applied in each management zone, based on the within-field spatial variability of plant vegetative vigour.
Practically, the method of remote sensing tested in this work can be used to identify two or three levels of vegetative vigour inside each plot, e.g., high, medium and low. Then, the grapes harvested in the two or three different vigour level areas can be processed for producing different wines having different quality.
However, the correlation between plant vegetative vigour and leaf water content, as well as that between crop yield and other crop and soil parameters, is still one of the challenges to be overcome in precision viticulture in order to produce spatially variable rate crop input application maps that can increasingly satisfy crop needs.