# Towards Vine Water Status Monitoring on a Large Scale
## Using Sentinel-2 Images

Eve Laroche-Pinel, Sylvie Duthoit, Mohanad Albughdadi, Anne D Costard, Jacques Rousseau, Véronique Cheret, Harold Clenet

► To cite this version:

Eve Laroche-Pinel, Sylvie Duthoit, Mohanad Albughdadi, Anne D Costard, Jacques Rousseau, et al. "Towards Vine Water Status Monitoring on a Large Scale Using Sentinel-2 Images." Remote Sensing, 2021, 13 (9), 21p. 10.3390/rs13091837. hal-03247343

HAL Id: hal-03247343

https://hal.inrae.fr/hal-03247343

Submitted on 3 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Towards Vine Water Status Monitoring on a Large Scale Using Sentinel-2 Images

Eve Laroche-Pinel <fnref n="1" />,<fnref n="2" />,<fnref n="3" />,* , Sylvie Duthoit <fnref n="1" />, Mohanad Albughdadi <fnref n="1" />, Anne D. Costard <fnref n="1" />, Jacques Rousseau <fnref n="4" />, Véronique Chéret <fnref n="2" />,<fnref n="3" /> and Harold Clenet <fnref n="2" />,<fnref n="3" /> ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Laroche-Pinel, E.; Duthoit S.; Albughdadi, M.; Costard A.D.; Rousseau, R.; Chéret, V.; Clenet, H. Towards Vine Water Status Monitoring on a Large Scale Using Sentinel-2 Images. *Remote Sens.* **2021**, *13*, 1837. https://doi.org/10.3390/rs13091837

Academic Editor: Frédéric Cointault

Received: 15 April 2021

Accepted: 4 May 2021

Published: 9 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Wine growing needs to adapt to confront climate change. In fact, the lack of water becomes more and more important in many regions. Whereas vineyards have been located in dry areas for decades, so they need special resilient varieties and/or a sufficient water supply at key development stages in case of severe drought. With climate change and the decrease of water availability, some vineyard regions face difficulties because of unsuitable variety, wrong vine management or due to the limited water access. Decision support tools are therefore required to optimize water use or to adapt agronomic practices. This study aimed at monitoring vine water status at a large scale with Sentinel-2 images. The goal was to provide a solution that would give spatialized and temporal information throughout the season on the water status of the vines. For this purpose, thirty six plots were monitored in total over three years (2018, 2019 and 2020). Vine water status was measured with stem water potential in field measurements from pea size to ripening stage. Simultaneously, Sentinel-2 images were downloaded and processed to extract band reflectance values and compute vegetation indices. In our study, we tested five supervised regression machine learning algorithms to find possible relationships between stem water potential and data acquired from Sentinel-2 images (bands reflectance values and vegetation indices). Regression model using Red, NIR, Red-Edge and SWIR bands gave promising result to predict stem water potential ($R^2$ = 0.40, RMSE = 0.26).

**Keywords:** vineyard; water status; Sentinel-2; precision viticulture

## 1. Introduction

Water management is an important factor for wine growing especially in the Mediterranean regions. Vineyards, in these regions, require a sufficient amount of water in key development stages such as flowering, veraison and berry growth [1,2]. However, episodes of drought have significantly impacted Mediterranean vineyards leading to yield loss and irreversible damages over the last few years [3,4]. These constraints increase the considerations on the development of irrigation systems [5]. However, in numerous regions water restrictions are set up due to the lack of water or to AOC (Protected Designation of Origin) specifications. It is therefore essential to manage efficiently water supply in order to stick to the recommendations to both preserve water resource as well as yield and quality targets on vineyards. In addition, considerable thoughts are ongoing in regions that are increasingly suffering from drought such as those mentioned in the review [6] which describes solutions to improve water use efficiency. Proposed solutions include changes in agronomic practice (increasing soil water storage capacity, reducing direct soil water loss, using cover crops) and genetic improvements. Adaptive strategies [7] are also used to cope with environmental stress such as precision viticulture (use of decision models

<footnote n="1">¹ TerraNIS, 12 Avenue de l'Europe, F-31520 Ramonville Saint-Agne, France; sylvie.duthoit@terranis.fr (S.D.); mohanad.albughdadi@terranis.fr (M.A.); anne.costard@terranis.fr (A.D.C.)</footnote>

<footnote n="2">² Ecole d'Ingénieurs de Purpan, 75 voie du TOEC, F-31076 Toulouse, France; veronique.cheret@purpan.fr (V.C.); harold.clenet@purpan.fr (H.C.)</footnote>

<footnote n="3">³ UMR DYNAFOR, INRAE, Université de Toulouse, F-31326 Castanet-Tolosan, France</footnote>

<footnote n="4">⁴ Groupe—Institut Coopératif du vin, La Jasse de Maurin, F-34970 Montpellier, France; jrousseau@icv.fr</footnote>

<footnote>* Correspondence: eve.laroche-pinel@terranis.fr; Tel.: +33-067-903-2721</footnote>
or field captors) or in field technique (application of protective compounds). Multi-level strategies are also proposed to promote sustainable water use from the plant management to the consumer [8].

In this context, decision-support tools play a very important role in water status monitoring for vineyards. These tools allow following the evolution of the vine water status through the season and focus water supply on some areas or plots which suffer the most of water stress. They also could legitimize authorization requests for water supply but also provide proof of reasonable use of irrigation. Furthermore, a spatialized information could be used to adapt some cultural practices and to map the evolution between different dates.

In order to do that, remote sensing could be advantageous in particular with the availability of Sentinel-2 (S2) constellation data [9]. Its free availability combined with its revisit frequency of 5 days are well suited for monitoring vineyard during all development stages, making it a powerful tool for a low cost service [10]. The other advantage of the S2 sensor is its spectral resolution with 3 bands in the visible domains (Blue, Green and Red), 3 bands in the Red-Edge (RE), 2 in the Near-Infrared (NIR) and 2 in the Shortwave Infrared (SWIR). This is a benefit compared to other satellites as those domains have the potential to provide critical insights on vegetation behaviour [10].

Various studies have demonstrated the capacity of S2 images to analyze vineyard vegetation. According to [11], S2 can be powerful to monitor vine development with a high temporal resolution and the utility of S2 spatial resolution have been proved by [12]. Moreover, an estimation of crop coefficient (Kc) was done with S2 time-series by [13] and another study already proved the feasibility of qualifying the impact of heatwaves on irrigated vineyard using S2 images ([14]). To our knowledge, only one recent study has evaluated the potential of S2 images to directly estimate vineyard water status [15]. This study highlighted good relationship between stem water potential (SWP) and vegetation indices (VI) using NIR and SWIR regions in 82 vineyard in 2017.

The SWIR domains are known to be related to water concentration in plants [16–19] as they correspond to the vibrations of the O–H stretch in water molecules [17]. The NIR domains also react to water status as they can be affected by leaf morphology and structure [20]. Another spectral domain often related to water content is the Red-Edge as it is related to leaf chlorophyll concentration [21–23], which can in turn be linked to the plant water concentration [24,25]. A previous study using hyperspectral field measurements has already validated these three spectral domains to be linked with vineyard stem water potential [26].

In order to monitor vegetation, VI are commonly used since decades starting with the Normalized Difference Vegetation Index (NDVI) [27]. To predict water status, a previous study relied on a relationship between NDVI and crop coefficients [28] but others studies highlighted VI with RE and SWIR positions to be more related to vine water status [15,29]. Besides, as it becomes easier to implement Artificial Intelligence algorithms, more and more studies focus on the use of machine learning that uses all available spectral bands. This allows to take into account all the information included in the spectra [30–32]. There are nevertheless still obstacles to their operational use, especially considering their integration into automatic processing chains already in place in companies for example.

This paper investigates the potential of using S2 images to predict and monitor vine water status throughout the growing season on a large scale with a high temporal resolution. The database used in this study is composed of 36 vine plots in Mediterranean regions in the south of France with several grape varieties and grass management over three years (2018, 2019, 2020). Stem water potentials allowing to estimate vine water status were measured in the field from pea-size to harvest and S2 L2A images were downloaded at the closest dates. The objective is to find the best machine learning algorithm and the best feature to use Vegetation Index or reflectance bands values to predict stem water potential values of Mediterranean vineyard monitored during 3 years.
## 2. Materials and Methods

Figure 1 presents the materials and methodology used in this study. The first step is the data acquisition including S2 images acquisition and vine water status measurements on vine subplots (Section 2.1). Then the data preprocessing (Section 2.2) consists of resampling and merging S2 images as well as average the SWP values at the subplot-level.

In order to predict the water status with S2 images, two feature types were extracted from the images. The first one is the reflectance values of the bands. The second feature type is a set of selected VIs that are known to be correlated with water status (see Section 2.3 for more details).

Finally, analyses were performed in two steps. The first step evaluates the performance of multiple supervised regression algorithms to predict the SWP using S2 bands or VIs (Section 2.4.3). The aim is to determine the best feature type (S2 bands values or VI) and the best algorithm to use in order to predict SWP values in terms of accuracy but also versatility for the whole dataset (the one that works best in most cases). The second step explores the best models (Section 2.4.4) in terms of features importance (which S2 bands are needed) and analyzes data distribution.

<figure><img src="image_5.png" /><figcaption>Figure 1. Flowchart of the method from the features extraction to the analysis step.</figcaption></figure>
## 2.1. Data Acquisition

### 2.1.1. Study Sites

The study was conducted over three years with thirty six vine plots in the South of France (Figure 2). Most of the data were acquired on Syrah variety, one of the most common red grape varieties in the Mediterranean region. However, in order to complete the database, several other varieties were added during the study. In 2018, the measurements were focused on eleven vine plots with Syrah variety in the Minervois AOC area near Carcassonne. In 2019, five vine plots have been chosen with Syrah, Grenache and Chardonnay varieties near Carcassonne, Beziers and Montpellier. In 2020, twenty vine plots with six vine varieties were monitored near these previous city but also near Marseille.

<figure><img src="image_6.png" /></figure>

Sources: Esri, Airbus DS, USGS, NGA, NASA, CGIAR, N Robinson, NCEAS, NLS, OS, NMA, Geodatastyrelsen, Rijkswaterstaat, GSA, Geoland, FEMA, Intermap and the GIS user community

<caption>**Figure 2.** Location of the study vine plots for 2018, 2019 and 2020 and position of the corresponding S2 tiles.</caption>

Plots were all about 1 ha in area and in majority ungrassed (67%) and are all managed in bilateral cordon with an inter-row distance of 2 or 2.5 m and an inter-vine distance of 1 m. An example of three inter-row management is given in Figure 3. Particular attention has been paid to choose plots with no evident sign of diseases such as esca, Black dead arm or chlorosis.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Inter-row management example. (a) Plot ungrassed. (b) Plot with grass one row on two. (c) Plot totally grassed.</figcaption></figure>

Our study spans over three years with different types of vintage being represented. The year 2018 was a particularly exceptional year in terms of precipitation until June (2nd wettest year since 1975 with a cumulative precipitation of 971 mm) but with a heatwave event in the beginning of August [33]. The year 2019 was a very dry year (4th driest year since 1950) with a cumulative rainfall of 400 mm, i.e., 2 times less than in 2018. Additionally
average temperatures from July to October were above normal with a difference of up to 2 °C [33]. The year 2020 was similar to 2019, i.e., quite dry, but with a more humid spring [34]. Table 1 gives details for the months related to the vine development. Particularities are highlighted in red for a lack of water or warm temperatures, blue for too much water and cool temperatures, green for normal data and brown for contrasting rainfall or temperatures.

<table><caption>**Table 1.** Summary of the weather from May to September for 2018, 2019 and 2020 adapted from technical reports of Departmental Council of the Hérault region. Particularities are highlighted in red for a lack of water or warm temperatures, blue for too much water and cool temperatures, green for normal data and brown for contrasting rainfall or temperatures.</caption><thead><tr><th>Year</th><th>Month</th><th>Rainfall</th><th>Temperatures</th><th>Events</th></tr></thead><tbody><tr><td rowspan="5">2018</td><td>May</td><td>Very excessive</td><td>Alternating hot and very cold</td><td>Thunderstorms, hail and snow</td></tr><tr><td>June</td><td>Very excessive</td><td>In accordance with season</td><td></td></tr><tr><td>July</td><td>No rainfall</td><td>Mild to warm</td><td></td></tr><tr><td>August</td><td>Normal</td><td>Very hot</td><td>Storms, heatwave the first week</td></tr><tr><td>September</td><td>Very low</td><td>Very hot</td><td></td></tr><tr><td rowspan="5">2019</td><td>May</td><td>Deficit</td><td>Cool to very cool</td><td>Heatwave from the 26th to the 28th</td></tr><tr><td>June</td><td>Deficit</td><td>Average to cool</td><td>Heatwave, hail and fire</td></tr><tr><td>July</td><td>Deficit</td><td>Very hot</td><td></td></tr><tr><td>August</td><td>Deficit</td><td>Hot to very hot</td><td></td></tr><tr><td>September</td><td>Contrasting</td><td>Warm to very warm</td><td></td></tr><tr><td rowspan="5">2020</td><td>May</td><td>Heterogeneous</td><td>Mild to hot</td><td></td></tr><tr><td>June</td><td>Heterogeneous</td><td>Cool, warm the last week</td><td></td></tr><tr><td>July</td><td>Poor</td><td>Warm to hot</td><td>Hail and fires</td></tr><tr><td>August</td><td>Heterogeneous</td><td>Warm</td><td></td></tr><tr><td>September</td><td>Heterogeneous</td><td>Warm</td><td></td></tr></tbody></table>

### 2.1.2. Experimental Design

In each vine plot, from one to six subplots of 20 × 20 m were defined for data acquisition. The number and location of the subplots were chosen according to three main criteria: (1) the S2 20 m pixel grid in order to have a good representation of the variability in the entire pixel, (2) plot soil and/or vegetation variability as reported by agronomic expert or vine-growers and (3) the time needed to perform the SWP measurements in the field (example in Figure 4). Vis for SWP measurements were selected in order to have a good representation of the entire pixel and according to the rows orientation inside the pixel.

Table 2 gives details about the number of plots and subplots for each year, by inter-row management and by variety. As expected, the Syrah variety is the most represented as it is the most common in the region but other grape varieties were added during the investigation in order to increase the variability within the database.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Example of 2 subplots in a plot and location of 10 SWP measurements in a subplot. The number and location of the subplots are chosen according to (1) the S2 20 m pixel grid in order to have a good representation of the variability in the entire pixel, (2) plot soil and/or vegetation variability as reported by agronomic expert or vine-growers and (3) the time needed to perform the SWP measurements in the field. Vis for SWP measurements were selected in order to have a good representation of the entire pixel and according to the rows orientation inside the pixel.</figcaption></figure>

<table><caption>**Table 2.** Number of grass, ungrass and total plot and subplot for each year and their grape variety.</caption><thead><tr><th rowspan="2"></th><th colspan="2">2018</th><th colspan="2">2019</th><th colspan="2">2020</th><th colspan="2">Total</th></tr><tr><th>Plot</th><th>Subplot</th><th>Plot</th><th>Subplot</th><th>Plot</th><th>Subplot</th><th>Plot</th><th>Subplot</th></tr></thead><tbody><tr><td>Total</td><td>11</td><td>18</td><td>5</td><td>32</td><td>20</td><td>53</td><td>36</td><td>103</td></tr><tr><td colspan="9">Inter-row management</td></tr><tr><td>Grass</td><td>4</td><td>7</td><td>0</td><td>0</td><td>8</td><td>20</td><td>12</td><td>27</td></tr><tr><td>No grass</td><td>7</td><td>11</td><td>5</td><td>32</td><td>12</td><td>33</td><td>24</td><td>76</td></tr><tr><td colspan="9">Grape Variety</td></tr><tr><td>Syrah</td><td>11</td><td>18</td><td>3</td><td>19</td><td>7</td><td>18</td><td>21</td><td>55</td></tr><tr><td>Grenache</td><td>0</td><td>0</td><td>1</td><td>6</td><td>1</td><td>2</td><td>2</td><td>8</td></tr><tr><td>Chardonnay</td><td>0</td><td>0</td><td>1</td><td>7</td><td>0</td><td>0</td><td>1</td><td>7</td></tr><tr><td>Merlot</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>3</td><td>1</td><td>3</td></tr><tr><td>Cabernet</td><td>0</td><td>0</td><td>0</td><td>0</td><td>2</td><td>5</td><td>2</td><td>5</td></tr><tr><td>Sauvignon</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>6</td><td>1</td><td>6</td></tr><tr><td>Caladoc</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>Mourvèdre</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td></tr></tbody></table>

### 2.1.3. Water Status

In order to assess the water status of the vine subplot, the stem water potential was measured on several vis (rSWP) as explained in [26]. The steps of the measurement are described below:

* Bagging mature leaves at 10 a.m. in order to close stomata and balance the sap pressure between plant and leaf;

* Removing leave with stem 3 or 4 h after, and set up quickly in the pressure chamber;

* Recording the pressure in MegaPascal (MPa) required to squeeze the first drop of sap out of the stem.

Between five and ten stems were measured (one stem by root) in each subplot as shown in Figure 4 which represents more than 3200 rSWP measurements in total for the three years. Measurements were carried out in maximum in 3 rows to optimize the time in
the field while trying to take into account the variability inside the subplot. Values range from 0 for water comfort to -2 MPa for extreme water stress.

### 2.1.4. Sentinel-2 Images

The S2 images were acquired over the studied vine plots and processed in L2A level [35,36] from the THEIA platform (theia.cnes.fr/, online accessed on the 6 January 2021). This platform was created at the end of 2012 by 9 French public institutions involved in earth observation and environmental sciences. This scientific and technical structure aims to facilitate the use of satellite images by making them available to the scientific community. It provides Sentinel-2 surface reflectances for each available date, corrected for atmospheric effects using the MAJA processing chain [37], and complemented with a cloud mask [38]. The images were chosen according to two criteria: (1) no cloud covering over the subplots and (2) a maximum of 5 days between the image date and the SWP measurement date. 80% of the samples had no more than 2 days apart between the field measurements and the date of the image. These multispectral optical images include 13 spectral bands in the visible, the near infra-red (NIR) and the shortwave-infrared (SWIR) spectral region [39]. 10 bands were used in this study with a spatial resolution of 10m and 20 m. The three other bands were not adapted for this study considering their spatial (60 m) and spectral resolution (in the atmospheric absorption domains for example). Details about the spectral and spatial resolutions of the bands used for the analysis are described in Table 3.

<table><caption>**Table 3.** Sentinel-2A spectral and spatial resolution of bands used in this study.</caption><thead><tr><td>Band</td><td>Spectral Region</td><td>Wavelength Range (nm)</td><td>Spatial Resolution</td></tr></thead><tbody><tr><td>B2</td><td>Blue</td><td>458–523</td><td>10 m</td></tr><tr><td>B3</td><td>Green</td><td>543–578</td><td>10 m</td></tr><tr><td>B4</td><td>Red</td><td>650–680</td><td>10 m</td></tr><tr><td>B5</td><td>Red-Edge 1</td><td>698–713</td><td>20 m</td></tr><tr><td>B6</td><td>Red-Edge 2</td><td>733–748</td><td>20 m</td></tr><tr><td>B7</td><td>Near Infrared</td><td>779–793</td><td>20 m</td></tr><tr><td>B8</td><td>Near Infrared</td><td>785–899</td><td>10 m</td></tr><tr><td>B8a</td><td>Near Infrared</td><td>855–875</td><td>20 m</td></tr><tr><td>B11</td><td>Shortwave Infrared</td><td>1565–1655</td><td>20 m</td></tr><tr><td>B12</td><td>Shortwave Infrared</td><td>2100–2280</td><td>20 m</td></tr></tbody></table>

In total, 306 S2 images were available from June to September on the 4 concerned tiles (Figure 2) for 2018, 2019 and 2020. Taking into account our 2 criteria (no cloud and close to field measurements dates), 38 images were kept. Details of the dates of these images for each year and each tile are given in Table 4.

<table><caption>**Table 4.** Dates of the S2 images used for 2018, 2019 and 2020 for each S2 tile. Images were chosen among all the images available as there is no cloud over the subplots and the date is close to the in field measurements (Maximum difference of 5 days).</caption><thead><tr><th>S2 Tile</th><th>Year of Study</th><th>June</th><th>July</th><th>August</th><th>September</th></tr></thead><tbody><tr><td rowspan="3">T31TDH/TDJ</td><td>2018</td><td>27th and 30th</td><td>30th</td><td>1st</td><td>8th and 10th</td></tr><tr><td>2019</td><td>17th</td><td>5th, 17th and 25th</td><td>14th</td><td></td></tr><tr><td>2020</td><td>19th and 21st</td><td>1st, 4th, 16th and 29th</td><td>10th, 15th and 25th</td><td></td></tr><tr><td rowspan="2">T31TEJ</td><td>2019</td><td>17th and 22nd</td><td>2nd and 17th</td><td>16th and 26th</td><td></td></tr><tr><td>2020</td><td>16th and 26th</td><td>6th, 11th, 16th, 26th and 31st</td><td>10th, 15th and 25th</td><td></td></tr><tr><td>T31TGJ</td><td>2020</td><td>28th</td><td>2nd</td><td></td><td></td></tr></tbody></table>
## 2.2. Sentinel-2 Images Preprocessing

The 20 m bands have been resampled to 10 m using the 8th band at 10 m resolution and then the ten bands have been merged in only one multispectral image. Figure 5 illustrates the data processing procedure.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Flowchart of the resampling and the merge of S2 bands used in this study. Example for an image with only one vine plot with three subplots. This method as been used for the 38 images and the 103 subplot.</figcaption></figure>

## 2.3. Features Extraction

### 2.3.1. Sentinel-2 Reflectance Bands Values

Once the ten S2 bands are merged together in one multispectral image, the values of each band are extracted for each subplot by averaging the value of the four 10 m pixels contained in the subplot. Pixels values of each band from THEIA platform have been divided by 10,000 in order to convert them to reflectance and to compute VI.

### 2.3.2. Vegetation Indices

Several VI have been chosen according to their ability to establish a link with the vine water status following a previous study [26]. The NDVI was also evaluated as it is a common VI usually used to monitor the vegetation. The details of the seven VI and their calculations formula are given in Table 5.

<table><caption>**Table 5.** Multispectral vegetation indices.</caption><thead><tr><td>Index (Abbreviation)</td><td>Main Use</td><td>Formula Used with S2 Bands</td><td>Reference</td></tr></thead><tbody><tr><td>Normalized Difference Vegetation Index (NDVI)</td><td>Vigor</td><td><sup>B<sub>8</sub>-B<sub>4</sub></sup>/<sub>B<sub>8</sub>+B<sub>4</sub></sub></td><td>[27]</td></tr><tr><td>Normalized Difference Red-Edge (NDRE 1/2)</td><td>Chlorophyll, Water</td><td><sup>B<sub>8</sub>-B<sub>5/6</sub></sup>/<sub>B<sub>8</sub>+B<sub>5/6</sub></sub></td><td>[40]</td></tr><tr><td>Inverted Red-Edge Chlorophyll Index (IRECI)</td><td>Chlorophyll</td><td><sup>B<sub>8</sub>-B<sub>4</sub></sup>/<sub>B<sub>5</sub></sub></td><td>[41]</td></tr><tr><td>Red-Edge Chlorophyll Absorption Index (RECAI)</td><td>Chlorophyll</td><td><sup>B<sub>8</sub>-B<sub>6</sub></sup>/<sub>B<sub>3</sub></sub> * <sup>B<sub>6</sub></sup>/<sub>B<sub>3</sub></sub></td><td>[42]</td></tr><tr><td>Normalized Difference Infrared Index (NDII)</td><td>Chlorophyll, Water</td><td><sup>B<sub>8</sub>-B<sub>11</sub></sup>/<sub>B<sub>8</sub>+B<sub>11</sub></sub></td><td>[43]</td></tr><tr><td>Red-Edge Position (REP)</td><td>Chlorophyll</td><td><sup>(<sup>B<sub>4</sub>+B<sub>8</sub></sup>/<sub>2</sub>) - B<sub>5</sub></sup>/<sub>B<sub>6</sub>-B<sub>5</sub></sub></td><td>[21]</td></tr><tr><td>Moisture Stress Index (MSI)</td><td>Water</td><td><sup>B<sub>11</sub></sup>/<sub>B<sub>8</sub></sub></td><td>[18]</td></tr></tbody></table>
## 2.4. Analysis

### 2.4.1. SWP Mean by Subplot

The five to ten rSWP measurements are averaged by subplot in order to reflect the pixel representation. In total, the database is composed of 349 observations, 50 for 2018, 146 for 2019 and 153 for 2020. An observation represents a subplot at a given date with its average value of SWP measurements.

### 2.4.2. Algorithms Description

Multiple supervised regression algorithms were used in the analysis (Table 6). These algorithms aim at mapping a set of observations (S2 reflectance values/VIs) with a continuous variable (SWP values). To be more specific, the K-nearest K-neighbors, ExtraTrees, Support Vector, linear and Bayesian Ridge model were evaluated. Supervised algorithms based on neural networks were avoided due to the limited number of training samples.

All the analysis have been performed using ScikitLearn Python package [44].

<table><caption>**Table 6.** Details of the algorithms used.</caption><thead><tr><td>Name of the Algorithms</td><td>Type of the Algorithm</td><td>Reference</td></tr></thead><tbody><tr><td>Kneighbors regressor</td><td>k-nearest neighbors</td><td>[45]</td></tr><tr><td>ExtraTrees regressor</td><td>Decision Tree</td><td>[46]</td></tr><tr><td>Support Vector Regressor</td><td>Support Vector Machine</td><td>[47]</td></tr><tr><td>Linear regression and BayesianRidge model</td><td>Linear Model</td><td>[48]</td></tr></tbody></table>

Each algorithm was tested (1) for all the database with all samples (2018, 2019, 2020, grass and ungrass), (2) for each year with grassed and ungrassed plots and (3) for each year with only ungrassed plots. Figure 6.

<figure><img src="image_10.png" /><figcaption>**Figure 6.** Illustration of the database sub-sample decomposition.</figcaption></figure>

### 2.4.3. Robustness and Precision Evaluation of Algorithms

Figure 7 gives the analysis details. In a first part, a test of the robustness and precision of each algorithms is done with each feature (S2 bands values and VI) and each sub-sample (Years and inter-row management). A stratified suffle split is used to split the dataset into 10 parts (k = 10) in order to test the stability of the model. Into each of the 10 parts, attention has been paid to keep samples from all the studied development stages (a). For each 10 parts, 75% of the sample are used for the training and 25% for the test.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Flowchart details of analysis (n=349 for the entire dataset, 1: pea-size, 2: Pre-veraison, 3: veraison, 4: ripening).</figcaption></figure>

* Regression score

For each split, the determination coefficient (R²) and the Root Mean Square Error (RMSE) are computed for each algorithm according to the following formulas:

$$R^2 = 1 - \frac{\sum_{i=1}^{n} (y_i - \hat{y}_i)^2}{\sum_{i=1}^{n} (y_i - \bar{y})^2} \quad (1)$$

$$RMSE = \sqrt{\frac{\sum_{i=1}^{n} (\hat{y}_i - y_i)^2}{n}} \quad (2)$$

where $y_i$ and $\hat{y}_i$ are the ith measurement and corresponding predicted value for i = {1, ..., n}, n is the number of measurements. and $\bar{y}$ is the mean of all the values. Note that the R² and RMSE are then averaged over the 10 splits.

In order to determine the best algorithm to use and the best feature type (S2 bands or VI), in addition to the model score, the best model is determined according to one more criterion: the versatility for all sub-samples (Years and inter-row management).

The results are shown for an RMSE > 0.3 (which corresponds to an average error of 30% for values ranging from 0 to -2 MPa).

### 2.4.4. Best Model Exploration

* Bands importance

It is important to highlight which bands of S2 are the most used in the algorithms to know which ones are the most interesting to predict vine water status.

The significance of the bands used to build the model is checked with the p-value and their coefficient in the model. Then, we recursively remove non-significant bands or bands with the smallest coefficients until scores (R² and RMSE) decline.

* Impact of experimental conditions on the result

The distribution of the values is verified with a single split of the data (75% for the training, 25% for the test) and the impact of inter-row management, years or study and variety are then examined.
# 3. Results

## 3.1. SWP Values

A box plot of the SWP values for each development stage is given in Figure 8. As expected, the mean values decrease during the season but differences are also observable for a same stage (from -0.4 MPa to -1.4 MPa at veraison for example).

<figure><img src="image_12.png" /><figcaption>**Figure 8.** Box plot of SWP values (in MegaPascal) measured in field for each development stage covered by our study. The size of the colored rectangles indicates the percentage of data per year.</figcaption></figure>

The detail for all statistics is given in Table 7 for each development stage of each year. The number of observations by development stage ranges from 15 to 49 according to the number of subplots, and the number of people available to do the measurements. The standard deviation ranges from 0.09 to 0.28 which indicates some significant variability in the measurements made in the same year at the same stage.

<table><caption>**Table 7.** Statistic description of SWP values for each development stage of each year (STD: standard deviation).</caption><thead><tr><th rowspan="2">Years</th><th colspan="4">2018</th><th colspan="4">2019</th><th colspan="4">2020</th></tr><tr><th>1</th><th>3</th><th>4</th><th>1</th><th>2</th><th>3</th><th>4</th><th>1</th><th>2</th><th>3</th><th>4</th></tr></thead><tbody><tr><td>Development Stage</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Nb_obstruction</td><td>15</td><td>17</td><td>18</td><td>26</td><td>39</td><td>32</td><td>49</td><td>27</td><td>48</td><td>37</td><td>41</td></tr><tr><td>Min</td><td>-0.68</td><td>-1.27</td><td>-1.74</td><td>-0.93</td><td>-1.17</td><td>-1.15</td><td>-1.66</td><td>-0.91</td><td>-1.34</td><td>-1.43</td><td>-1.66</td></tr><tr><td>Max</td><td>-0.33</td><td>-0.50</td><td>-0.60</td><td>-0.28</td><td>-0.53</td><td>-0.54</td><td>-0.62</td><td>-0.42</td><td>-0.49</td><td>-0.72</td><td>-0.80</td></tr><tr><td>Median</td><td>-0.46</td><td>-0.91</td><td>-1.20</td><td>-0.52</td><td>-0.92</td><td>-0.92</td><td>-1.41</td><td>-0.63</td><td>-0.81</td><td>-1.09</td><td>-1.32</td></tr><tr><td>Mean</td><td>-0.47</td><td>-0.87</td><td>-1.14</td><td>-0.55</td><td>-0.91</td><td>-0.89</td><td>-1.34</td><td>-0.63</td><td>-0.81</td><td>-1.06</td><td>-1.31</td></tr><tr><td>STD</td><td>0.09</td><td>0.21</td><td>0.28</td><td>0.18</td><td>0.16</td><td>0.16</td><td>0.25</td><td>0.13</td><td>0.15</td><td>0.18</td><td>0.25</td></tr></tbody></table>

## 3.2. Evaluation of the Robustness and Precision

The mean $R^2$ and RMSE of the 10 splits for each algorithm are computed for each feature type (VIs and bands values) for each sub-sample. An example of results is shown in Table 8 for the whole dataset (all years and all inter-row management). In this example, the best result is obtained using all S2 bands values with Bayesian Ridge and Linear Regression ($R^2$ = 0.40 and RMSE = 0.26).
<table><caption>**Table 8.** Results of scores exploration for all data (all years, all inter-row management). Best scores are highlighted in bold.</caption><thead><tr><th rowspan="2">Features</th><th rowspan="2">Scores</th><th colspan="5">Algorithms</th></tr><tr><th>K-Nearest K-neighbors</th><th>ExtraTrees</th><th>Support Vector</th><th>Linear Model</th><th>Bayesian Model</th></tr></thead><tbody><tr><td rowspan="2">All S2 bands</td><td>R2</td><td>0.29</td><td>0.39</td><td>0.37</td><td><strong>0.40</strong></td><td><strong>0.40</strong></td></tr><tr><td>RMSE</td><td>0.28</td><td>0.26</td><td>0.26</td><td><strong>0.26</strong></td><td><strong>0.26</strong></td></tr><tr><td rowspan="2">NDVI</td><td>R2</td><td>0.03</td><td>0.03</td><td>0.03</td><td>0.04</td><td>0.04</td></tr><tr><td>RMSE</td><td>0.34</td><td>0.39</td><td>0.33</td><td>0.33</td><td>0.33</td></tr><tr><td rowspan="2">NDRE1</td><td>R2</td><td>0.07</td><td>0.05</td><td>0.03</td><td>0.004</td><td>0.003</td></tr><tr><td>RMSE</td><td>0.34</td><td>0.41</td><td>0.33</td><td>0.33</td><td>0.33</td></tr><tr><td rowspan="2">NDRE2</td><td>R2</td><td>0.14</td><td>0.06</td><td>0.02</td><td>0.04</td><td>0.04</td></tr><tr><td>RMSE</td><td>0.35</td><td>0.42</td><td>0.33</td><td>0.33</td><td>0.33</td></tr><tr><td rowspan="2">IRECI</td><td>R2</td><td>0.02</td><td>0.04</td><td>0.11</td><td>0.11</td><td>0.11</td></tr><tr><td>RMSE</td><td>0.34</td><td>0.39</td><td>0.31</td><td>0.31</td><td>0.31</td></tr><tr><td rowspan="2">RECAI</td><td>R2</td><td>0.02</td><td>0.04</td><td>0.02</td><td>0.02</td><td>0.01</td></tr><tr><td>RMSE</td><td>0.36</td><td>0.33</td><td>0.34</td><td>0.34</td><td>0.33</td></tr><tr><td rowspan="2">NDII</td><td>R2</td><td>0.04</td><td>0.05</td><td>0.07</td><td>0.09</td><td>0.09</td></tr><tr><td>RMSE</td><td>0.34</td><td>0.41</td><td>0.32</td><td>0.32</td><td>0.32</td></tr><tr><td rowspan="2">REP</td><td>R2</td><td>0.02</td><td>0.03</td><td>0.15</td><td>0.19</td><td>0.17</td></tr><tr><td>RMSE</td><td>0.33</td><td>0.39</td><td>0.31</td><td>0.30</td><td>0.30</td></tr><tr><td rowspan="2">MSI</td><td>R2</td><td>0.09</td><td>0.06</td><td>0.07</td><td>0.10</td><td>0.10</td></tr><tr><td>RMSE</td><td>0.34</td><td>0.41</td><td>0.31</td><td>0.31</td><td>0.31</td></tr></tbody></table>

Table 9 summarizes all the best algorithms (best $R^2$ and RMSE) for each sub-sample (years and inter-row management). Only results with an $R^2 \ge 0.25$ and RMSE $\le 0.3$ have been reported, the only vegetation index meeting these conditions is REP (with Red, NIR and Red-Edge bands) and only for 2019. Using all S2 bands reflectance always gives a better result than this VI. The linear and Bayesian model appear to be the most accurate models in most cases (4 times on 7). The results seem better with only the plot without grass cover ($R^2$ up to 0.48 for all plots in 2020 and to 0.58 for plots without grass in 2019). The result per year is often higher for $R^2$ than with all years aggregated but the RMSE does not decrease so much and never falls below 0.2.

The best algorithm that will be explored next is the linear regression since it appeared to be the best in all cases with the whole dataset (all years with or without grass).

<table><caption>**Table 9.** Summary of the best algorithms for every subplot with $R^2 \ge 0.25$ and RMSE $\le 0.3$.</caption><thead><tr><th>Year</th><th>Inter-Row Management</th><th>REP</th><th>All S2 Bands</th></tr></thead><tbody><tr><td rowspan="2">All</td><td>All</td><td>$R^2 < 0.25$</td><td>$R^2 = 0.40$<br>RMSE = 0.26<br>Bayesian Ridge model and Linear model</td></tr><tr><td>No grass</td><td>$R^2 < 0.25$</td><td>$R^2 = 0.48$<br>RMSE = 0.24<br>Linear model</td></tr><tr><td rowspan="2">2018</td><td>All</td><td>$R^2 < 0.25$</td><td>$R^2 < 0.25$<br>$R^2 = 0.52$</td></tr><tr><td>No grass</td><td>$R^2 < 0.25$</td><td>RMSE = 0.2<br>Extra Tree regressor</td></tr></tbody></table>
<table><caption>**Table 9. Cont.**</caption><tr><td>2019</td><td>No grass</td><td>R<sup>2</sup> = 0.27<br>RMSE = 0.29<br>Bayesian Ridge model</td><td>R<sup>2</sup> = 0.58<br>RMSE = 0.22<br>Linear model</td></tr><tr><td>2020</td><td>All</td><td>R<sup>2</sup> &lt; 0.25</td><td>R<sup>2</sup> = 0.48<br>RMSE = 0.21<br>Bayesian Ridge model and Linear model</td></tr><tr><td></td><td>No grass</td><td>R<sup>2</sup> &lt; 0.25</td><td>R<sup>2</sup> = 0.56<br>RMSE = 0.21<br>Bayesian Ridge Model</td></tr></table>

For the linear model with the complete dataset, the values of R² and RMSE seem steady over the 10 splits as R² ranges from 0.43 to 0.55 and RMSE ranges from 0.22 to 0.28 (Figure 9).

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Boxplot of R² and RMSE over the 10 split for Linear model with the whole dataset (all years and all inter-row management).</figcaption></figure>

### 3.3. Best Models Exploration

#### 3.3.1. Bands Importance

According to their p-values (<0.05), all S2 bands were significant. In order to test if all the bands are really necessary to give a good result, we progressively removed the band with the smallest coefficient until an impact on the R² and RMSE scores is observed (Section 2.4.4). The result shown in Figure 10 highlights that only 4 bands are needed to perform as well as with 12 bands. Those 4 required bands are B4, B6, B8a and B12 in the Red, Red-Edge, NIR and SWIR regions, respectively.
<figure><img src="image_14.png" /><figcaption>Figure 10. Best band to use for regression model.</figcaption></figure>

### 3.3.2. Data Distribution and Impact of Experimental Conditions

This section aims at exploring the data distribution and the impact on the predictions when considering three experimental conditions: 1- Grass cover management, 2- Year of study, 3- Grapes variety. Comparison between observed and predicted SWP were plotted to test for discrepancies between reality and model's result obtained using the Linear Model based on the 4 bands of S2 highlighted in Section 3.3.1 (B4, B6, B8a and B12).

As a primary analysis, it can be noticed in Figure 11 that predicted SWP are a little overestimated until -1.4 MPa and a threshold effect seems to appear after.

* Inter-row management

Figure 11 shows the distribution of the predicted and observed SWP depending on the presence of grass in the inter-row. The scatterplot with and without grass superimpose pretty well.

<figure><img src="image_15.png" /><figcaption>Figure 11. Distribution of predicted and observed data according to inter-row management.</figcaption></figure>

* Grape variety

Figure 12 gives the data distribution according to grapes varieties. Even if the number of observations per grape variety is not the same Table 2, particular attention was paid
during the protocol so that measurements cover the whole range of potential water status (i.e., with just 4 points for the Mourvèdre, SWP values range from -0.7 MPa to -1 MPa). One can also mention that the scatterplot of all varieties superimpose pretty well.

<figure><img src="image_16.png" /><figcaption>Figure 12. Data distribution according to grapes variety.</figcaption></figure>

* Development stage

Figure 13 shows the distribution of predicted and observed SWP for each development stage. As expected, the relationship between development stage and SWP (i.e., increase of water stress (low SWP) as the season progresses from pea size to ripening) appears clearly. The second point, is the very large point cloud at the ripening stage with observed SWP spanning from -0.6 to -1.8. The prediction seems to be better for pre-veraison and veraison stage (data are closer to the straight line).

<figure><img src="image_17.png" /><figcaption>Figure 13. Data distribution according to development stage.</figcaption></figure>

* Year of study

The observation of Figure 14 highlights the good distribution of the data for the 3 years with SWP values ranging from water comfort (-0.4 MPa) to water stress (-1.8 MPa) for each year.
<figure><img src="image_18.png" /><figcaption>Figure 14. Data distribution according to study year.</figcaption></figure>

3.4. *Comparison between NDVI, Best VI (REP) and Best Model with the Four S2 Bands*

In this study, both bands values from S2 images and VI were tested to predict vine
SWP. According to the result, the commonly used NDVI (with B4 and B8) is not relevant at
all (Figure 15a) and the only significant VI is the REP (with B4, B8, B5 and B6) but with a
low R² of 0.19 for the whole dataset (Figure 15b). This VI gives less effective results than
the model using the following S2 bands B4, B6, B8a and B12 (Figure 15c).

<figure><img src="image_19.png" /><figcaption>Figure 15. *Cont.*</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>**Figure 15.** Distribution of the observed SWP according to (a) NDVI values, (b) REP values, (c) values of the model developed in this study.</figcaption></figure>

# 4. Discussion

## 4.1. *Significant Features*

As previously reported in the literature, the most discriminant bands, highlighted in Section 3.3.1, are located in Red-Edge (B6) [49], NIR (B8a) [50] and SWIR (B12) [51]. The Red domain (B4) is also highlighted in this study as in [26]. This domain is related to pigments which also react to water stress (xanthophyll and chlorophyll) according to [52]. Except the REP, no other vegetation index appears to be interesting enough in our study to predict the vine water status with the algorithms we tested. Therefore the direct use of S2 bands seems to be more relevant.

## 4.2. *Robustness of the Model*

A first observation to highlight is that a threshold effect seems to appear for values to predict lower than -1.4 MPa. This can be explained by the fact that -1.4 MPa already reflects a very high water stress [53].

In the following paragraphs we discuss the impact of the inter-row grass management, as the diversity in grape varieties, the development stage and vintage, on the robustness of the model. Overall results show that in terms of prediction there is clearly a trend in the data distribution despite varying conditions. Specific effects are discussed in the 3 next paragraphs.

### 4.2.1. Impact of Grass Cover

As shown in Table 9, the model performs better when only plots without grass in the inter-row are considered ($R^2$ = 0.48 and RMSE = 0.24 vs. $R^2$ = 0.40 and RMSE = 0.26). However, as can be seen in Figure 11, the results when also considering plots with grass show that only few points deviate from the trend. It can be explained by the fact that in Mediterranean region the vegetation in the inter-row is not always well developed especially in the middle of summer as shown by Figure 3. It would be interesting to extend the study with further investigation in the Bordelais or Burgundy regions where grass in the inter-row is more developed.

### 4.2.2. Impact of Grape Variety

Most of the measurements were made on Syrah grape variety (Table 2) because it is the most common one in the study region. Nevertheless additional varieties were added in 2019 and 2020 to expend the scope of our initial survey. This way we were able to include both isohydric (Cabernet-Sauvignon, Grenache, Caladoc) and anisohydric grape varieties (Syrah, Chardonnay, Merlot, Mourvèdre) which behave differently to water stress.
Anisohydric varieties maintain their stomata opened to ensure photosynthesis and can be damaged in case of strong water stress while isohydric varieties close their stomata when the first signs of drought appeared and thus limit their photosynthesis [54]. Despite these distinct behaviors, we do not see any real differences in the predictions as shown in Figure 12 and the model seems to be able to predict correctly the SWP both for isohydric and anisohydric varieties. To extend the study it would be interesting to focus on more isohydric varieties and to add other varieties common in other regions (e.g., Merlot in the Bordelais region).

### 4.2.3. Impact of Development Stage

Following the trends observed in Figure 13, the model appears to be less efficient when considering data at ripening stage. This may be due to other physiological mechanisms that affect the plant when under high water stress, such as sugar accumulation in berries and modification in leaf area to fruit ratios [55], that could also affect the optical properties of the leaves. The objective of the study was to develop a model that could be used for images acquired at any stage of vine development. Nevertheless, tests were made without data from ripening stage to evaluate the impact on predictions. The models results were ultimately not better (max R² = 0.30), indicating that general model is more relevant, even if it is more accurate only when considering pre-veraison and veraison stage. This is rather positive those stages are the ones where winegrowers have to be careful on water stress as it can delay the good development of the berries.

### 4.2.4. Impact of Years

Despite differences in precipitations and temperatures (detailed at the end of Section 2.1.1), one can see in Figure 14 that no clear difference can be observed between the predictions for the three years covered by our study.

## 4.3. *Considerations for a Future Operational Service*

The results of this study brings promising perspectives to develop a future operational service to monitor vineyard water status with Sentinel-2. However, when dealing with S2 images the question of spatial resolution is crucial and special attention will be needed to different points in order to be able to implement its use for a potential future service:

* Avoid the pixels at the edge which can cover both the plot and a path or a forest bordering it for example. In order to avoid its so-called "mixed" pixels, an advice would be to apply a buffer around the edge of the plot of at least 5 m to keep and interpret only the pixels fully included in the plot,

* Be careful with the size and shape of the plot in order to have a consistent number of pixels to interpret inside the plot,

* Consider the inter-row management and the soil management and/or the soil composition which can impact the observed signal since vineyard is conducted in row and the inter-row is also visible with 10 or 20 m pixel. Maybe think to use two models according to the grass cover management (one for ungrassed field and one for grassed) by improving the number of grassed vine field in the database.

In order to bring a new operational service it will be necessary to take into account these considerations but also to be able to add data over the years to improve the robustness of the model.

A next step is in progress with winegrowers in order to evaluate the spatial and temporal coherence of water status maps obtained with the model using the four S2 bands.

# 5. Conclusions

This study investigates the feasibility of using S2 images to monitor vine water status. The dataset is composed of 349 samples acquired during 3 years in the Mediterranean region. Water status in field measurements (SWP) and S2 images were acquired all along the vegetation cycle (from pea size to ripening) in 2018, 2019 and 2020. Our analyses allow
to highlight the best features to use from S2 images to predict SWP values: the four bands values in Red, Red-Edge, NIR and SWIR spectral domains. The regression linear model gives promising results to predict vine water status from different years, grape varieties and grass management. The model found will be more relevant for veraison stage than ripening. The characteristics of S2 make it possible to use this model to monitor vine water status at a large spatial scale (in a whole region for example) with 10m resolution and a little time of revisit (5 days). However, it is necessary to take into account the particularity of vine management in row and not to forget the potential impact of grass cover or soil management/composition to the result interpretation. The first winegrowers feedbacks are encouraging and further investigations will be done to increase the dataset and improve the model with other grape varieties and grass cover-management for example.

**Author Contributions:** Conceptualization, E.L.-P., H.C. and S.D.; methodology, E.L.-P., H.C., S.D. and J.R.; validation, H.C., S.D. and V.C.; formal analysis, E.L.-P. and M.A.; investigation, E.L.-P., S.D. and J.R.; resources, E.L.-P., H.C., S.D. and J.R.; data curation, E.L.-P., A.D.C. and M.A.; writing—original draft preparation, E.L.-P., H.C. and S.D.; writing—review and editing, E.L.-P., H.C., S.D., M.A. and A.D.C.; visualization, E.L.-P.; supervision, H.Clenet, S.D. and V.C.; project administration, E.L.-P.; funding acquisition, S.D. and H.C. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by ANRT, TerraNIS, ICV and EI Purpan.

**Acknowledgments:** The authors thank all the trainees of ICV and TerraNIS involved in field measurements acquisition. They also thank winegrowers for allowing the realization of measurements in their fields and the Minervois Union in particular Marie Vidal-Vigneron.

**Conflicts of Interest:** The authors declare no conflict of interest.

## Abbreviations

The following abbreviations are used in this manuscript:

MDPI Multidisciplinary Digital Publishing Institute

DOAJ Directory of open access journals

NIR Near Infrared

NDVI Normalized Difference Vegetation Index

S2 Sentinel-2

SWP Stem Water Potential

SWIR Short Wave Infra-Red

VI Vegetation indices

REP Red-Edge Position

RMSE Root Mean Square Error

## References

1. Ramos, M.C.; Pérez-Álvarez, E.P.; Peregrina, F.; Martínez de Toda, F. Relationships between grape composition of Tempranillo variety and available soil water and water stress under different weather conditions. *Sci. Hortic.* **2020**, *262*, 109063. [CrossRef]

2. van Leeuwen, C.; Darriet, P. The Impact of Climate Change on Viticulture and Wine Quality. *J. Wine Econ.* **2016**, *11*, 150–167. [CrossRef]

3. Chaves, M.; Santos, T.; Souza, C.; Ortuño, M.; Rodrigues, M.; Lopes, C.; Maroco, J.; Pereira, J. Deficit irrigation in grapevine improves water-use efficiency while controlling vigour and production quality. *Ann. Appl. Biol.* **2007**, *150*, 237–252. [CrossRef]

4. Costa, J.M.; Ortuño, M.F.; Chaves, M.M. Deficit Irrigation as a Strategy to Save Water: Physiology and Potential Application to Horticulture. *J. Integr. Plant Biol.* **2007**, *49*, 1421–1434. [CrossRef]

5. Gorguner, M.; Kavvas, M.L. Modeling impacts of future climate change on reservoir storages and irrigation water demands in a Mediterranean basin. *Sci. Total Environ.* **2020**, p. 141246. [CrossRef] [PubMed]

6. Medrano, H.; Tomás, M.; Martorell, S.; Escalona, J.M.; Pou, A.; Fuentes, S.; Flexas, J.; Bota, J. Improving water use efficiency of vineyards in semi-arid regions. A review. *Agron. Sustain. Dev.* **2015**, *35*, 499–517. [CrossRef]

7. Bernardo, S.; Dinis, L.T.; Machado, N.; Moutinho-Pereira, J. Grapevine abiotic stress assessment and search for sustainable adaptation strategies in Mediterranean-like climates. A review. *Agron. Sustain. Dev.* **2018**, *38*. [CrossRef]

8. Costa, J.M.; Vaz, M.; Escalona, J.; Egipto, R.; Lopes, C.; Medrano, H.; Chaves, M.M. Modern viticulture in southern Europe: Vulnerabilities and strategies for adaptation to water scarcity. *Agric. Water Manag.* **2016**. [CrossRef]
9. Addabbo, P.; Focareta, M.; Marcuccio, S.; Votto, C.; Ullo, S.L. Contribution of Sentinel-2 data for applications in vegetation monitoring. *Acta Imeko* **2016**, *5*, 44. [CrossRef]

10. Segarra, J.; Buchaillot, M.L.; Araus, J.L.; Kefauver, S.C. Remote Sensing for Precision Agriculture: Sentinel-2 Improved Features and Applications. *Agronomy* **2020**, *10*, 641. [CrossRef]

11. Devaux, N.; Crestey, T.; Leroux, C.; Tisseyre, B. Potential of Sentinel-2 satellite images to monitor vine fields grown at a territorial scale. *OENO One* **2019**, *53*, 51–58. [CrossRef]

12. Di Gennaro, S.F.; Dainelli, R.; Palliotti, A.; Toscano, P.; Matese, A. Sentinel-2 validation for spatial variability assessment in overhead trellis system viticulture versus UAV and agronomic data. *Remote Sens.* **2019**, *11*, 2573. [CrossRef]

13. Rozenstein, O.; Haymann, N.; Kaplan, G.; Tanny, J. Estimating cotton water consumption using a time series of Sentinel-2 imagery. *Agric. Water Manag.* **2018**, *207*, 44–52. [CrossRef]

14. Cogato, A.; Pagay, V.; Marinello, F.; Meggio, F.; Grace, P.; De, M. Assessing the feasibility of using medium-resolution imagery information to quantify the impact of the heatwaves on irrigated vineyards. *Remote Sens.* **2019**, *11*, 2869. [CrossRef]

15. Cohen, Y.; Gogumalla, P.; Bahat, I.; Netzer, Y.; Ben-Gal, A.; Lenski, I.; Michael, Y.; Helman, D. Can time series of multispectral satellite images be used to estimate stem water potential in vineyards? In *Precision Agriculture '19; Wageningen Academic: Wageningen, The Netherlands, 2019; pp. 445–451*. [CrossRef]

16. Seelig, H.D.; Hoehn, A.; Stodieck, L.S.; Klaus, D.M.; Adams, W.W.; Emery, W.J. The assessment of leaf water content using leaf reflectance ratios in the visible, near-, and short-wave-infrared. *Int. J. Remote Sens.* **2008**, *29*, 3701–3713. [CrossRef]

17. Zovko, M.; Zibrat, U.; Knapic, M.; Bubalo Kovacic, M.; Romic, D. Hyperspectral remote sensing of grapevine drought stress. *Precis. Agric.* **2019**, *20*, 335–347. [CrossRef]

18. Ceccato, P.; Flasse, S.; Tarantola, S.; Jacquemoud, S.; Grégoire, J.M. Detecting vegetation leaf water content using reflectance in the optical domain. *Remote Sens. Environ.* **2001**, *77*, 22–33. [CrossRef]

19. Das, B.; Sahoo, R.N.; Pargal, S.; Krishna, G.; Verma, R.; Viswanathan, C.; Sehgal, V.K.; Gupta, V.K. Evaluation of different water absorption bands, indices and multivariate models for water-deficit stress monitoring in rice using visible-near infrared spectroscopy. *Spectrochim. Acta Part A Mol. Biomol. Spectrosc.* **2020**, *119104*. [CrossRef]

20. Kim, D.; Zhang, H.; Zhou, H.; Du, T.; Wu, Q.; Mockler, T.; Berezin, M. Highly sensitive image-derived indices of water-stressed plants using hyperspectral imaging in SWIR and histogram analysis. *Sci. Rep.* **2015**, *5*, 15919. [CrossRef] [PubMed]

21. Main, R.; Cho, M.A.; Mathieu, R.; O'Kennedy, M.M.; Ramoelo, A.; Koch, S. An investigation into robust spectral indices for leaf chlorophyll estimation. *ISPRS J. Photogramm. Remote Sens.* **2011**, *66*, 751–761. [CrossRef]

22. Behmann, J.; Steinrücken, J.; Plümer, L. Detection of early plant stress responses in hyperspectral images. *ISPRS J. Photogramm. Remote Sens.* **2014**, *93*, 98–111. [CrossRef]

23. Matese, A.; Di Gennaro, S.F. Technology in precision viticulture: A state of the art review. *Int. J. Wine Res.* **2015**. [CrossRef]

24. Tilling, A.K.; O'Leary, G.J.; Ferwerda, J.G.; Jones, S.D.; Fitzgerald, G.J.; Rodriguez, D.; Belford, R. Remote sensing of nitrogen and water stress in wheat. *Field Crops Res.* **2007**, *104*, 77–85. [CrossRef]

25. Li, M.; Chu, R.; Yu, Q.; Islam, A.R.M.T.; Chou, S.; Shen, S. Evaluating Structural, Chlorophyll-Based and Photochemical Indices to Detect Summer Maize Responses to Continuous Water Stress. *Water* **2018**, *10*, 500. [CrossRef]

26. Laroche-Pinel, E.; Albughdadi, M.; Duthoit, S.; Chéret, V.; Rousseau, J.; Clenet, H. Understanding Vine Hyperspectral Signature through Different Irrigation Plans: A First Step to Monitor Vineyard Water Status. *Remote Sens.* **2021**, *13*, 536. [CrossRef]

27. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring vegetation systems in the Great Plains with ERTS. In *Third ERTS Symposium, NASA SP-351*; NASA: Greenbelt, MD, USA, 1973; pp. 309–317.

28. Rienth, M.; Scholasch, T. State-of-the-art of tools and methods to assess vine water status. *OENO One* **2019**, *53*, 619–637. [CrossRef]

29. Romero, M.; Luo, Y.; Su, B.; Fuentes, S. Vineyard water status estimation using multispectral imagery from an UAV platform and machine learning algorithms for irrigation scheduling management. *Comput. Electron. Agric.* **2018**, *147*, 109–117. [CrossRef]

30. Ferrant, S.; Selles, A.; Le Page, M.; Herrault, P.A.; Pelletier, C.; Al-Bitar, A.; Mermoz, S.; Gascoin, S.; Bouvet, A.; Saqalli, M.; et al. Detection of irrigated crops from Sentinel-1 and Sentinel-2 data to estimate seasonal groundwater use in South India. *Remote Sens.* **2017**, *9*, 1119. [CrossRef]

31. Loggenberg, K.; Strever, A.; Greyling, B.; Poona, N. Modelling water stress in a Shiraz vineyard using hyperspectral imaging and machine learning. *Remote Sens.* **2018**, *10*, 202. [CrossRef]

32. Schmitter, P.; Steinrücken, J.; Römer, C.; Ballvora, A.; Léon, J.; Rascher, U.; Plümer, L. Unsupervised domain adaptation for early detection of drought stress in hyperspectral images. *ISPRS J. Photogramm. Remote Sens.* **2017**, *131*, 65–76. [CrossRef]

33. Conseil départemental de l'Hérault. *Annales Climatologiques et Hydrologiques 2019*; Technical report; DREAL Occitanie; DREAL: Marseille, France, 2019.

34. Conseil départemental de l'Hérault. *Info Clim, Synthèse 2020*; Technical Report 252; Conseil départemental de l'Hérault; DREAL: Marseille, France, 2020.
35. Thales Alenia Space. *Sentinel-2 Products Specification Document*; Technical Report S2-PDGS-TAS-DI-PSD; ESA: Paris, France, 2016. Available online: https://sentinels.copernicus.eu/web/sentinel/document-library/latest-documents/-/asset publisher/EgUy8 pfXboLO/content/sentinel-2-level-1-to-level-1c-product-specifications;jsessionid=8BE6EE17FECEE9CDECD948BD1F6A852 2.jvm2?redirect=https%3A%2F%2Fsentinels.copernicus.eu%2Fweb%2Fsentinel%2Fdocument-library%2Flatest-documents% 3Bjsessionid%3D8BE6EE17FECEE9CDECD948BD1F6A8522.jvm2%3Fp_p_id%3D101_INSTANCE_EgUy8pfXboLO%26p_p_ lifecycle%3D0%26p_p_state%3Dnormal%26p_p_mode%3Dview%26p_p_col_id%3Dcolumn-1%26p_p_col_pos%3D1%26p_p_ col_count%3D2 (accessed on 10 April 2021).

36. Donadieu, J.; L'Helguen, C. *Sentinel-2A L2A Products Description*; Technical Report PSC-NT-411-0362-CNES; CNES, Toulouse, France, 2016. Available online: https://laboobs-mip.fr/wp-content-labo/uploads/sites/19/2016/09/PSC-NT-411-0362-CNES_ 01_00_SENTINEL-2A_L2A_Products Description.pdf (accessed on 10 April 2021).

37. Hagolle, O.; Huc, M.; Pascual, D.V.; Dedieu, G. A multi-temporal and multi-spectral method to estimate aerosol optical thickness over land, for the atmospheric correction of FormoSat-2, LandSat, VENUS and Sentinel-2 images. *Remote Sens.* **2015**, *7*, 2668-2691. [CrossRef]

38. Baetens, L.; Desjardins, C.; Hagolle, O. Validation of copernicus Sentinel-2 cloud masks obtained from MAJA, Sen2Cor, and FMask processors using reference cloud masks generated with a supervised active learning procedure. *Remote Sens.* **2019**, *11*, 433. [CrossRef]

39. ESA. *Sentinel-2 Spectral Response Functions (S2-SRF)*; Technical Report COPE-GSEG-EOPG-TN-15-0007; ESA: Paris, France, 2017. Available online: https://sentinel.esa.int/web/sentinel/user-guides/sentinel-2-msi/document-library/-/asset publisher/Wk0 TKajISaR/content/sentinel-2a-spectral-responses (accessed on 10 April 2021).

40. Gitelson, A.A.; Viña, A.; Ciganda, V.; Rundquist, D.C.; Arkebauer, T.J. Remote estimation of canopy chlorophyll content in crops. *Geophys. Res. Lett.* **2005**, *32*, 1-4. [CrossRef]

41. Frampton, W.J.; Dash, J.; Watmough, G.; Milton, E.J. Evaluating the capabilities of Sentinel-2 for quantitative estimation of biophysical variables in vegetation. *ISPRS J. Photogramm. Remote Sens.* **2013**, *82*, 83-92. [CrossRef]

42. Cui, B.; Zhao, Q.; Huang, W.; Song, X.; Ye, H.; Zhou, X. A New Integrated Vegetation Index for the Estimation of Winter Wheat Leaf Chlorophyll Content. *Remote Sens.* **2019**, *11*, 974. [CrossRef]

43. Lemaire, G.; Jeuffroy, M.H.; Gastal, F. Diagnosis tool for plant and crop N status in vegetative stage. Theory and practices for crop N management. *Eur. J. Agron.* **2008**, *28*, 614-624. [CrossRef]

44. Pedregosa, F.; Varoquaux, G.; Gramfort, A.; Michel, V.; Thirion, B.; Grisel, O.; Blondel, M.; Prettenhofer, P.; Weiss, R.; Dubourg, V.; et al. Scikit-learn: Machine Learning in Python. *J. Mach. Learn. Res.* **2011**, *12*, 2825-2830.

45. Kramer, O. K-Nearest Neighbors. In *Dimensionality Reduction with Unsupervised Nearest Neighbors*; Springer: Berlin/Heidelberg, Germany, 2013; pp. 13-23. [CrossRef]

46. Kazemitabar, J.; Amini, A.; Bloniarz, A.; Talwalkar, A.S. Variable Importance Using Decision Trees. In *Advances in Neural Information Processing Systems 30*; Guyon, I., Luxburg, U.V., Bengio, S., Wallach, H., Fergus, R., Vishwanathan, S., Garnett, R., Eds.; Curran Associates Inc.: Greenbelt, MD, USA, 2017; pp. 426-435.

47. Awad, M.; Khanna, R. Support Vector Regression. In *Efficient Learning Machines: Theories, Concepts, and Applications for Engineers and System Designers*; Apress: Berkeley, CA, USA, 2015; pp. 67-80. [CrossRef]

48. Pilz, J. Bayesian estimation and experimental design in linear regression models. *J. Am. Stat. Assoc.* **1992**, *87*, 1250.

49. Ballester, C.; Zarco-Tejada, P.J.; Nicolás, E.; Alarcón, J.J.; Fereres, E.; Intrigliolo, D.S.; Gonzalez-Dugo, V. Evaluating the performance of xanthophyll, chlorophyll and structure-sensitive spectral indices to detect water stress in five fruit tree species. *Precis. Agric.* **2017**, 1-16. [CrossRef]

50. Maimaitiyiming, M.; Ghulam, A.; Bozzolo, A.; Wilkins, J.L.; Kwasniewski, M.T. Early detection of plant physiological responses to different levels of water stress using reflectance spectroscopy. *Remote Sens.* **2017**, *9*, 1-23. [CrossRef]

51. Das, B.; Mahajan, G.R.; Singh, R. Hyperspectral Remote Sensing: Use in Detecting Abiotic Stresses in Agriculture. *Adv. Crop. Environ. Interact.* **2018**, 317-335. [CrossRef]

52. Rapaport, T.; Hochberg, U.; Shoshany, M.; Karnieli, A.; Rachmilevitch, S. Combining leaf physiology, hyperspectral imaging and partial least squares-regression (PLS-R) for grapevine water status assessment. *ISPRS J. Photogramm. Remote Sens.* **2015**, *109*, 88-97. [CrossRef]

53. Ojeda, H.; Saurin, N. L'irrigation de précision de la vigne : méthodes, outils et stratégies pour maximiser la qualité et les rendements de la vendange en économisant de l'eau. *Innov. Agron.* **2014**, *38*, 97-108.

54. Simonneau, T.; Ollat, N.; Pellegrino, A.; Lebon, E. Contrôle de l'état hydrique dans la plante et réponses physiologiques de la vigne à la contrainte hydrique. *Innov. Agron.* **2014**, *38*, 13-82.

55. Suter, B.; Destrac Irvine, A.; Gowdy, M.; Dai, Z.; van Leeuwen, C. Adapting Wine Grape Ripening to Global Change Requires a Multi-Trait Approach. *Front. Plant Sci.* **2021**, *12*, 1-17. [CrossRef] [PubMed]