Feasibility of Using the Two-Source Energy Balance Model (TSEB) with Sentinel-2 and Sentinel-3 Images to Analyze the Spatio-Temporal Variability of Vine Water Status in a Vineyard

In viticulture, detailed spatial information about actual evapotranspiration (ETa) and vine water status within a vineyard may be of particular utility when applying site-specific, precision irrigation management. Over recent decades, extensive research has been carried out in the use of remote sensing energy balance models to estimate and monitor ETa at the field level. However, one of the major limitations remains the coarse spatial resolution in the thermal infrared (TIR) domain. In this context, the recent advent of the Sentinel missions of the European Space Agency (ESA) has greatly improved the possibility of monitoring crop parameters and estimating ETa at higher temporal and spatial resolutions. In order to bridge the gap between the coarse-resolution Sentinel-3 thermal and the fine-resolution Sentinel-2 shortwave data, sharpening techniques have been used to downscale the Sentinel-3 land surface temperature (LST) from 1 km to 20 m. However, the accurate estimates of high-resolution LST through sharpening techniques are still unclear, particularly when intended to be used for detecting crop water stress. The goal of this study was to assess the feasibility of the two-source energy balance model (TSEB) using sharpened LST images from Sentinel-2 and Sentinel-3 (TSEB-PTS2+3) to estimate the spatio-temporal variability of actual transpiration (T) and water stress in a vineyard. T and crop water stress index (CWSI) estimates were evaluated against a vine water consumption model and regressed with in situ stem water potential (Ψstem). Two different TSEB approaches, using very high-resolution airborne thermal imagery, were also included in the analysis as benchmarks for TSEB-PTS2+3. One of them uses aggregated TIR data at the vine+inter-row level (TSEB-PTairb), while the other is based on a contextual method that directly, although separately, retrieves soil and canopy temperatures (TSEB-2T). The results obtained demonstrated that when comparing airborne Trad and sharpened S2+3 LST, the latter tend to be underestimated. This complicates the use of TSEB-PTS2+3 to detect crop water stress. TSEB-2T appeared to outperform all the other methods. This was shown by a higher R2 and slightly lower RMSD when compared with modelled T. In addition, regressions between T and CWSI-2T with Ψstem also produced the highest R2.


Introduction
In a scenario of climate change and water scarcity, viticulturists will have to rely more on efficient irrigation management and adopting regulated deficit irrigation (RDI) strategies for successful grape production and wine quality [1][2][3][4]. To achieve these targets, accurate estimates of actual evapotranspiration (ET a ) and water status at the sub-field scale are necessary, particularly for specialty crops. Water requirements within a vineyard are variable due to differences in soil components [5][6][7] and administering uniform irrigation may result in either the over-or under-watering of some areas within the vineyard. Obtaining detailed spatial information about ET a and vine water status within a given vineyard may therefore be of particular interest, particularly for applying site-specific, precision irrigation management.
Over recent decades, remote sensing techniques have been widely used for mapping and monitoring ET a at the field scale, using satellite and airborne imagery. Several satellites have been used to estimate ET a based on thermal infrared (TIR) data, whose spatial and temporal resolutions may vary [8]. Satellites such as Landsat (60-120 m), MODerate resolution Imaging Spectrometer (MODIS) (1 km), ASTER (90 m) and GOES (4 km) have been widely applied for mapping ET a in irrigated agricultural fields [9][10][11][12]. Although extensive research has been carried out on ET a estimates made with satellites, one of the major limitations of this approach is the coarse spatial resolution of satellite TIR images. In this context, the recent advent of Sentinel-2 and Sentinel-3 missions from the European Space Agency (ESA), as part of the Copernicus program (http://www.copernicus.eu/), has greatly improved the possibilities of monitoring crop parameters and estimating ET a at higher temporal and spatial resolutions. Although Sentinel-3 carries a TIR sensor with a coarse spatial resolution of 1000 m at nadir, several other techniques could be used to sharpen the lower spatial resolution TIR observations with higher spatial resolution visible and near-infrared (VISNIR) observations [13][14][15][16]. These techniques involve empirically relating TIR to VISNIR spectral signals within the same scene, at the coarse pixel resolution of the thermal band, and then applying this relationship to the fine pixel resolution VISNIR bands in order to produce sharpened thermal band imagery at the same, with high resolution [16]. Guzinski et al. [17] evaluated the feasibility of sharpening daily Sentinel-3 (SLSTR sensor~1000 m TIR data) satellite data with Sentinel-2 MIS (20 m VISNIR data). They obtained an acceptable level of accuracy and improvements in the estimates of fluxes, with respect to low-resolution satellite data. However, since the sharpened high-resolution land surface temperature (LST) is estimated from statistical relations using optical data (e.g., the correlation between LST and the NIR band), low LST used to correspond with pixels with high biomass and vice versa. One inherent limitation is the need to detect crop water stress in cases in which crops, and particularly those with heterogeneous or complex canopies, are stressed but without provoking a reduction in their leaf biomass.
On the other hand, very high-resolution thermal imagery, usually acquired from unmanned aerial vehicles (UAVs) or manned aircrafts, can potentially provide the spatial resolutions required for precision irrigation applications. Most current applications mainly focus on detecting crop water status in various woody trees [18][19][20][21] or on scheduling irrigation based on estimates of leaf water potential (Ψ L ) [22]. This last study estimated Ψ L from a crop water stress index (CWSI) using an empirical approach which consisted of developing a baseline by normalizing differences between the canopy (T c ) and air temperatures (T a ) based on evaporative demand (by means of the vapor pressure deficit) [23,24]. Other studies have also used very high-resolution thermal infrared and multispectral data from airborne platforms to map ET a over a vineyard [25][26][27]. One of the main shortcomings associated with using very high-resolution imagery would, however, be the need to provide this information on a routine basis over relatively large regions, given the need to calibrate this by crop, phenology and weather conditions [28].
Many different surface energy balance models have been widely used to estimate consumptive water use from ET a (SEBAL, [29][30][31]; METRIC, [32]; TSEB, [33]; ALEXI, [34]; DISALEXI, [35]; DTD, [36]; ETEML, [37]). However, precision irrigation should focus on crop water fluxes and hence models are required that have the ability to separate ET a into its component parts: soil evaporation (E) and canopy Remote Sens. 2020, 12,2299 3 of 25 transpiration (T). One such model is the two-source energy balance (TSEB) land surface scheme [33], which has proven be fairly robust for a wide range of landscape and weather conditions [38][39][40]. The performance of the TSEB with satellite data relies on accurate observations of directional radiometric land surface temperature (LST). Canopy (T c ) and soil (T s ) temperatures are then estimated in an iterative process in which it is first assumed that the green canopy (expressed as the fraction of LAI that is green, f g ) transpires at a potential rate, which is based on the Priestley-Taylor formulation [41] (TSEB-PT), and that after that, a step-wise reduction of canopy transpiration can be performed until realistic daytime fluxes for both soil and canopy values are retrieved.
While waiting for the launch of a high-resolution TIR mission with a high revisit time, and given the open data access policies of the Copernicus satellites, it would seem reasonable to think-at least a priori-that the use of the TSEB-PT approach with Sentinel data and the adoption of disaggregation methods for downscaling spatial resolution could offer the most feasible option for a routine way of assessing water status at serviceable spatial and temporal scales and be used for precision irrigation purposes on an operational basis.
This paper evaluates the use of TSEB-PT with sharpened LST images from Sentinel-2 and Sentinel-3 for supporting operational irrigation decision-making in a vineyard throughout a full growing season. It also evaluates the accuracy of ET a and T estimates with respect to very high-resolution airborne imagery and in situ vine physiological measurements. To the best of our knowledge, no other study has benchmarked the feasibility of using this methodology in an operational setting or quantified the feasibility of using it to detect vines under water stress. Based on the previously mentioned studies, which pointed to a good level of agreement between CWSI and Ψ L using very high-resolution airborne imagery and demonstrated its workability for scheduling irrigation, one hypothesis of this study focuses on whether the theoretically based CWSI-which is defined as one minus the ratio of T, over potential transpiration (T 0 ) [24]-calculated with Sentinel data could provide an appropriate approach for detecting differences in vine water status.

Retrieval of TSEB Approaches
This study first evaluated estimates of actual transpiration (T) obtained with the TSEB model. We did this using three different methodologies which all make use of land surface temperature (LST): (i) TSEB-PT S2+3 (PT for Priestley-Taylor), which uses sharpened LST from Sentinel-2 and Sentinel-3; (ii) TSEB-PT airb , which is similar to TSEB-PT S2+3 , but which obtains LST from very high-resolution airborne thermal imagery; and (iii) TSEB-2T, which is a contextual approach that directly uses T c and T s , obtained from very high-resolution thermal imagery. Due to the difficulty involved in validating ET a and its separate components, owing to the lack of eddy covariance or sap flow sensors, a spatio-temporal assessment of both actual and potential transpiration (T and T 0 respectively) was conducted through a vine water consumption model. Secondly, we analyzed different methodologies used to obtain the crop water stress index (CWSI) and related the results to stem water potential (Ψ stem ) measurements. We also used maps showing the spatial distribution of the ET components to assess the accumulated water consumption of the whole vineyard.

Two-Source Energy Balance (TSEB) Model
The TSEB land surface energy balance scheme has the advantage that it partitions estimates of sensible heat fluxes into canopy and soil layers. These depend on soil and canopy temperatures (T s and T c , respectively). The approach is therefore able to estimate canopy transpiration and soil evaporation separately. TSEB is based on the resistance network of the Shuttleworth and Wallace (1985) energy combination model, which has also been used in this study for computing potential transpiration (T 0 ). For more information, the full Python code is available online (https://github.com/hectornieto/pyTSEB,

Priestley-Taylor Iterative Retrieval, TSEB-PT
The PT approach is used when TSEB is run at a coarse resolution and it is not possible to separate temperatures from pure canopy (T c ) and bare soil (T s ) pixels. In such cases, the observed directional radiometric temperature, T RAD (θ), is partitioned into soil and canopy temperature (T s and T c , respectively), based on the vegetation cover fraction at the thermal sensor view angle, f c (θ): T c and T s are used to separately compute the surface energy budgets for the soil and canopy components where R NS is the net radiation at the soil surface; R NC is the net radiation in the canopy layer; LE S is the soil evaporation rate; LE C is the canopy transpiration rate; G is the soil heat flux; and H S and H C are the soil and canopy sensible heat fluxes (Equations (4) and (5)), which are directly related to T s and T c , respectively. ρ is the air density; c p is the heat capacity of the air; T ac is the aerodynamic temperature, or the air temperature at the sink-source height; and r s and r x are, respectively, the bulk canopy and soil resistances to heat transport. In order to derive T c and T s from Equation (1), TSEB uses a first approximation of canopy latent heat flux, LE C . This includes a first guess of canopy transpiration at the potential rates retrieved from the Priestley and Taylor model [33] (Equations (4) and (5)): where α PT is the Priestley-Taylor coefficient, which was initially set to 1.26 [43] but subsequently automatically reduced for stress conditions; ∆ is the slope of the saturation vapor pressure vs. temperature; and γ is the psychrometric constant. With this initial value of LEc, the other soil and canopy flux components are retrieved from Equations (2)-(5). However, given that plants may be water-limited, TSEB-PT iteratively reduces the canopy transpiration value until realistic fluxes are obtained (i.e., allowing for non-negative latent heat fluxes during daytime, [33]) by looping between Equations (1)-(5).

Data Sharpening Scheme
In this study, the PT approach was compared using LST obtained by sharpening Sentinel-2 and Sentinel-3 images (TSEB-PT S2+3 ) and from very high-resolution airborne thermal images (TSEB-PT airb ). The data sharpening method used in this study was based on the data mining sharpener (DMS) introduced by Gao et al. [16] and used by Guzinski et al. [13,17] with Sentinel 2 and Sentinel 3 data. This approach was also used by the European Space Agency (ESA) to develop the evapotranspiration plugin for the Sentinel Application Platform (SNAP) (www.esa-sen4et.org). The methodology consists of improving the coarse spatial resolution in LST retrievals from Sentinel-3 by sharpening the thermal band imagery using shortwave multispectral data. In addition, a regression tree approach was used to relate LST to a suite of shortwave spectral reflectance and ancillary data (e.g., elevation, illumination), at a coarse resolution, and then to apply it to a fine (shortwave band) pixel resolution to produce sharpened thermal band imagery. In this study, Sentinel-3 LST was sharpened to 20 m using nine bands of Sentinel-2. More information about the approach and code used is available online (https://github.com/radosuav/pyDMS, last accessed: 22 September 2019).

Contextual TSEB (TSEB-2T)
If high-resolution thermal imagery is available, it is possible to obtain T s and T c directly, without the need to compute an initial canopy transpiration [27]. In this study, we obtained the average canopy and soil temperature for each 20 × 20 m grid (coinciding with pixels of Sentinel-2) from very high-resolution airborne thermal images acquired on three different dates. The selection of pure canopy vegetation and soil pixels was obtained using a supervised classification method that took advantage of the semi-automatic classification plugin (SCP) of QGIS (QGIS 2.18.22) [44]. More details of the TSEB-2T model can be found in [27].

Biophysical Parameters of the Vegetation and Ancillary Data
The biophysical parameters of the vegetation, such as its leaf area index (LAI) and the fraction of vegetation cover (f c ) used for running the TSEB and Shuttleworth and Wallace (S-W) models, were retrieved from atmospherically corrected Sentinel-2 reflectances with a spatial resolution of 20 m, and by using the biophysical processor available in the SNAP software v.6.0.1 [45]. This method creates a synthetic dataset of biophysical parameters covering the whole possible range of expected values; the reflectances of these pixels are then simulated using radiative transfer models (Prospect+SAIL, [46]) and convolved to the Sentinel-2 bands. Once these two datasets (of biophysical parameters and reflectances) have been generated, a neuronal network is built in order to predict the biophysical parameters from the Sentinel-2 band reflectances. The resulting estimates were used with the three methodologies tested (TSEB-PT S2+3 , TSEB-PT airb and TSEB-2T). Other structural parameters of the vegetation, including its canopy height (h c ) and leaf width or canopy width to height ratio (w/h c ), were measured for the 32 vines that were distributed in the vineyard. Meteorological inputs were obtained from the local weather station, which was located 1 km from the study site [47].

Study Site
The study was conducted during the 2018 growing season, at a 6.6-ha vineyard of the Tempranillo grape variety, located at Raimat (Lleida, Spain) (41 • 41 19.40" N; 0 • 29 37.16" E, 280 m elevation) ( Figure 1). The climate in the area is Mediterranean, and its average accumulated rainfall and reference evapotranspiration (ET ref ) during the growing season (1 April to 30 September) were 214 and 793 mm, respectively. The vines were 21 years old and were planted with a spacing of 1.6 × 2.5 m and a NE-SW row orientation. The vines were cordon-trained through a vertical shoot positioning (VSP) trellis system, with a bilateral, spur-pruned cordon located at a height of 0.9 m above the ground. The soil had a loamy texture, with an effective depth ranging from 0.6 to 1.2 m. The vines were drip-irrigated using drippers spaced 0.60 m apart along a single drip line per vine row. These lines discharged 4 l h −1 , with an irrigation frequency which varied from 2 to 3 days per week. The vineyard had a single irrigation sector and its irrigation scheduling was managed by the Raimat winery. The vine water requirements were calculated based on a water balance method to replace crop evapotranspiration (ET c ), in which ET c was calculated from the ET ref of Penman-Monteith [48] and crop coefficients derived from previous experiments [49,50], as follows: K c1 = 0.2 (bud break on 5 April), K c2 = 0.5 (bloom on 25 May to fruit-set on 5 June), K c3 = 0.8 (mid-season, from veraison on 20 July until harvest) and K c4 = 0.3 (at leaf fall at the end of October). During post-veraison, a regulated deficit irrigation (RDI) strategy was applied that consisted of reducing the amount of water administered by up to 50%. A water meter (CZ2000-3M Contazara, Zaragoza, Spain) was installed on the dripline to measure the volume of irrigation water applied. The total amount of water applied in the vineyard throughout the growing season was 280 mm. Disease control and vine nutrition management were

Satellite Data and Outputs of the TSEB model
A total of 28 Sentinel-2 and 44 Sentinel-3 cloud-free scenes were downloaded for the period 25 th March to 15 th October 2018 (Table 1). Sentinel-2 images were downloaded at level-1 (top of atmosphere reflectance) and atmospherically corrected using the Sen2Cor radiometric corrector [51]. These images were then resampled at 20 m in order to obtain the biophysical parameters of the vegetation. The land surface temperature (LST) was first obtained at the Sentinel-3 resolution, following the split-window approach proposed by Sobrino et al. [52]. Land surface energy fluxes and LST sharpening were derived for those days with available Sentinel-3 LST images that were within a 4-day difference of an S2 scene. The evapotranspiration components obtained were: actual evapotranspiration (ETa); actual transpiration (T); soil evaporation (E); crop water stress index (CWSI), computed as one minus the ratio between the actual and potential transpiration; and the partitioning of the evapotranspiration expressed as the T/ETa ratio. All the fluxes and metrics were integrated into daily values. To achieve this, the instantaneous latent flux estimates, expressed in W m -2 , were upscaled to daily water fluxes, expressed in units of mm/day, by multiplying the instantaneous ratio between the latent heat flux and solar irradiance by the average daily solar irradiance [53].

Satellite Data and Outputs of the TSEB model
A total of 28 Sentinel-2 and 44 Sentinel-3 cloud-free scenes were downloaded for the period 25 March to 15 October 2018 (Table 1). Sentinel-2 images were downloaded at level-1 (top of atmosphere reflectance) and atmospherically corrected using the Sen2Cor radiometric corrector [51]. These images were then resampled at 20 m in order to obtain the biophysical parameters of the vegetation. The land surface temperature (LST) was first obtained at the Sentinel-3 resolution, following the split-window approach proposed by Sobrino et al. [52]. Land surface energy fluxes and LST sharpening were derived for those days with available Sentinel-3 LST images that were within a 4-day difference of an S2 scene. The evapotranspiration components obtained were: actual evapotranspiration (ET a ); actual transpiration (T); soil evaporation (E); crop water stress index (CWSI), computed as one minus the ratio between the actual and potential transpiration; and the partitioning of the evapotranspiration expressed as the T/ET a ratio. All the fluxes and metrics were integrated into daily values. To achieve this, the instantaneous latent flux estimates, expressed in W m −2 , were upscaled to daily water fluxes, expressed in units of mm/day, by multiplying the instantaneous ratio between the latent heat flux and solar irradiance by the average daily solar irradiance [53]. Table 1. Acquisition dates of the Sentinel 2 and Sentinel-3 cloud-free images used in the analysis (Sentinel 2: tile 31TBG, Sentinel 3: 6059) and airborne thermal imagery. For each Sentinel 2 acquisition, there is generally one or more Sentinel 3 acquisitions, with a maximum delay of 4 days relative to the Sentinel-2 acquisition date.

Shuttleworth-Wallace (S-W) Model and Crop Water Stress
In order to be consistent with TSEB, the two-layer Shuttleworth-Wallace (S-W) model [54] was used to estimate the potential evapotranspiration (ET 0 ) and to separate its component parts. The theoretical base of the S-W model is provided by the Penman-Monteith energy combination equation, which has two parts: one for the soil surface, and the other for the plant surface. The potential evapotranspiration and transpiration were computed with the S-W model by setting a minimum stomatal resistance value of 100 sm −1 . These values were then used as the base for estimating the theoretical CWSI.
In this study, CWSI was also computed using a baseline empirical approach (CWSI e ) [23]. This consisted of normalizing the difference between the air and canopy temperatures (T c − T a ) using the vapor pressure deficit (VPD). In this study, the following baselines, developed by Bellvert et al. [55] for "Tempranillo" grapevines, were used to calculate the CWSI e : Lower limit: y = −1.780x + 1.253 Upper limit: y = 0.466x + 5.317 (8)  where y corresponds to T c − T a , and x represents VPD. A seasonal spatio-temporal analysis was needed to better understand the feasibility of using this methodology throughout a whole growing season. However, due to the lack of available S2 and S3 data for certain days in the growing season (due to cloud cover), gaps had to be filled using daily continuous estimates of ET, which were obtained by calculating a crop water stress coefficient (K s ). K s was calculated as the ratio between ET a and ET 0 , or as T/T 0 , for those days with available satellite data. It was then extrapolated for those days within the same week with gaps in satellite data. This was achieved by multiplying the previously obtained value for K s by the ET 0 or T 0 corresponding to dates with cloud cover.

Airborne Campaign
The airborne campaign consisted of three flights conducted on 18 (DOY 199) and 31 July (DOY 212) and 22 August (DOY 234) 2018. The vapor pressure deficits (VPD) and air temperatures (T a ) at the moment of image acquisition were 2.49 kPa and 29.9 ºC for DOY 199, 2.18 kPa and 31.2 ºC for DOY 212 and 2.51 kPa and 31.4 ºC for DOY 234. It is important to highlight that during the period from 12 to 17 August, there was a series of rainfall events at the study site which accounted for a total of 20 mm of precipitation. The aeroplane used was equipped with a thermal camera (FLIR SC655, FLIR Systems, Wilsonville, OR, USA) with a resolution of 640 x 480 pixels and a 13.1 mm focal length, with an FOV angle of 45 • . The spectral response was in the range of 7.5-13 µm. The times of acquisition were between approximately 10:30 and 11:00 UTC (coordinated universal time), coinciding with the S3 satellite passing overhead. The plane flew over at a height of 180 m agl, capturing images with a resolution of 0.25 m pixels. At the beginning of each growing season, the thermal camera was radiometrically calibrated by the manufacturer (FLIR Systems s.r.l) using a very stable blackbody radiator. During the flight, in situ temperatures were acquired over homogeneous cold and hot targets, using a portable IR-gun (Fluke 572 infrared thermometer, Everett, WA, USA) in order to correct the atmospheric contribution to the signal. Image orthorectification and mosaicking were performed using off-the-shelf photogrammetric software. This was completed with the aid of ground control point (GCP) targets, after their geolocations on the ground had been collected with a handheld GPS (Global Positioning System) (Geo7x, Trimble GeoExplorer series, Sunnyvale, CA, USA).

Field Measurements
Physiological and structural measurements of the vines were conducted at a total of 32 points which were uniformly distributed across the vineyard. Each point related to three vines. Midday stem water potential (Ψ stem ) was measured for each vine (two measurements per point), using a pressure chamber (Soil Moisture plant water status console 3005 Corp. Sta. Barbara, CA, USA). To ensure a balance between the leaf and the stem attached to it, shaded leaves located near the main trunk were wrapped in plastic bags covered with aluminum foil for 1 h prior to taking the Ψ stem measurements. All the measurements were taken in less than 1 hour, with their timing coinciding with that of the image acquisition on the airborne campaign days. The fraction of PAR light intercepted by the crop (fiPAR) and the LAI were determined as indicators of vegetative development, using a portable ceptometer (Accupar, Decagon Devices Inc., Pullman, WA, USA). Measurements were taken on each vine, using a 12-point grid that was measured at ground level and by covering the vine spacing distance. Incident radiation readings were measured above the vines. There was a total of seven measuring days, three of which coincided with the image acquisition dates. Daily fiPAR (fiPAR d ) was calculated using the hourly model for light interception [56] and by adjusting the porosity parameter in Oyarzun's model so that the noon fiPAR value of the hourly simulation matched that measured in the field concomitant to where the image was acquired. Structural parameters such as vine height, and vine width, both perpendicular to and along rows, were also measured, using a ruler. Effective LAI was obtained from fiPAR, using the approach presented by Norman et al. [57] and setting a leaf absorptivity for light of 0.85.

Spatio-Temporal Validations with the Vine Water Consumption Model
The estimates of actual and potential transpiration obtained with the TSEB and S-W models throughout the growing season were compared with those generated by a vine water consumption crop model. This model has been developed for grapevines and successfully used for conducting precision irrigation scheduling at a commercial vineyard over two consecutive growing seasons. A full description of the model has been provided by Bellvert et al. [58]. Although this model is able to simulate evapotranspiration, in this study, we only used it for estimates of transpiration. The reason for this is that soil evaporation is related to the total size of the wet soil surface, and this depends on the timing of applying irrigation. Since water meter readings were taken manually, on a weekly basis, we did not have all the information required to accurately estimate the soil evaporation component. In situ measurements of the fraction of intercepted radiation (fiPAR) and the midday stem water potential (Ψ stem ) are required in order to obtain modelled values for T and T 0 .
The spatio-temporal assessments of crop stress using the different methodologies (TSEB-PT S2+3 , TSEB-PT airb and TSEB-2T) were conducted for 32 vines and on three different image acquisition dates (DOY 199, 212 and 234). In addition, the temporal evaluation of TSEB-PT S2+3 was carried out at two representative extreme zones within the vineyard. This was done throughout the growing season and covered well-watered and stressed grapevines.

Spatial Assessment of Transpiration
Spatially distributed, remotely sensed estimates of actual and potential transpiration were obtained using three different methodologies and compared with those modelled using a vine water consumption model. Table 2 presents a summary of the statistics used to evaluate the performance of the different methods. It shows that estimates of T 0 obtained with the S-W model yielded results that closely agreed with the modelled T 0 values for all three dates and that when the data were analyzed together, R 2 ranged from 0.71 to 0.78. TSEB-2T also appeared to outperform all the other models, yielding closer agreement with modelled T. This was shown by the higher R 2 value and slightly lower root mean square deviation (RMSD) for the three measured dates and also by encompassing all the data together. This suggested that whenever it is possible to retrieve T c and T s directly from very high-resolution thermal imagery, the accuracy of the estimates of T could significantly improve in comparison with other tested methodologies. Although the transpiration derived from TSEB-PT airb did not provide the same accuracy as fluxes derived from TSEB-2T, it provided improved values compared with TSEB-PT S2+3 . In addition to lower R 2 values, the latter also exhibited slopes that significantly differed from 1 in all cases, suggesting that this approach tends to overestimate low transpiration rates. We also noted that all three methodologies (TSEB-PT S2+3 , TSEB-PT airb and TSEB-2T) tended to give slightly lower T values at rates higher than approximately 2.5 mm day −1 , with this bias being more pronounced for the two methods which used airborne thermal imagery. This could probably be explained as follows: (i) in TSEB-2T, the average T c for all the vines within the 20 × 20 m pixel (~20 vines) could have been higher than it should have been for the three measured vines due to either the selection of mixed pixels combining vegetation and soil or because of within-pixel spatial variability in water status; (ii) as TSEB-PT airb uses the Priestley-Taylor approach and this is based upon an initial guess that assumes that the vines transpire at their potential rate, it is essential to have a robust characterization of the canopy potential transpiration; and (iii) bias is more pronounced in methods which use airborne thermal imagery (TSEB-2T and TSEB-PT airb ) than in TSEB-PT S2+3 because the sharpened LST of the latter tends to be overestimated ( Figure 2) and, as a consequence, the estimated T rates tend to be higher. Table 2. Accuracy statistics for regressions between modelled transpiration, obtained from the vine water consumption model, and estimates of potential and actual transpiration. The latter were respectively estimated using: the Shuttleworth and Wallace model (T 0 S-W); two-source energy balance (TSEB) with (i) Priestley-Taylor with sharpened LST from Sentinel-2 and Sentinel-3 (T TSEB-PT S2+3 ), (ii) PT with LST obtained from very high-resolution airborne thermal imagery resampled at 20 m (T TSEB-PT airb ) and (iii) the separation of canopy (T c ) and soil (T s ) temperatures obtained from very high-resolution airborne thermal imagery (T TSEB-2T), and correlations with the crop water stress index (CWSI), calculated as 1-T/T 0 . CWSI e represents the CWSI obtained with the empirical baseline approach for cv.Tempranillo, using the following equations: y = −1.780x + 1.253 and y = 0.466x + 5.317 (Bellvert et al. 2015), to establish the lower and upper limits. Root mean square deviation (RMSD) and bias (modelled without the estimated values) are expressed in mm day −1 ; the coefficient of variation (C V -RMSE divided by the mean of the modelled values), equation (estimated = a modelled + b) and coefficient of determination (R 2 ) are unitless. Only statistically significant regressions between the remote sensing results and water consumption model are shown.   As with transpiration, the modelled CWSI also agreed with the values obtained using the contextual (CWSI-2T) method, although the closest agreement was obtained using the baseline empirical approach (CWSI e ) ( Table 2). In both cases, the results were robust for DOY 199 and 212 with the RMSD and R 2 respectively ranging from 0.11 to 0.22 and from 0.57 to 0.71. When the data for all three dates were encompassed, CWSI-2T was also the approach with the highest R 2 and lowest RMSD values. However, DOY 234 showed only a low level of agreement for all the methodologies except CWSI-PT S2+3 . This was possibly due to the effect of the 20 mm of rainfall received in the week before the image was acquired. The effects of convective rainfall occurring a few days before image acquisition on the CWSI have been previously described by Bellvert et al. [22], who observed that rainfall may modify the VPD and air temperature within the vineyard as the amount of water evaporated from the soil and vegetation surfaces is high. In consequence, the evaporative demand within the vineyard may be lower than at the location where the meteorological data were gathered and this may affect the accuracy of estimates of CWSI. Another hypothesis was that T s was clearly lower than for previous days due to the soil being wetter: this was detected by very high-resolution imagery. Partitioning readings into E and T would have probably tended to overestimate E and to underestimate T, thereby resulting in an overestimation of the CWSI values. The regressions obtained using the CWSI-PT S2+3 method were not significant, with the partial exception of DOY 234 and when all three dates were plotted together. It seems that although the overall RMSD values were low, they tended to increase as the CWSI increased. This suggests that this approach was unable to detect severely water-stressed vines. On the other hand, the result for DOY 234 was significant, probably because of there being a lower degree of water stress, as apparently reflected by the smaller range in CWSI values associated with the rainfall event. In fact, Figure 2 shows the regression between LST with sharpened S2 and S3 and airborne LST, resampled at 20 m. This shows that estimates of high-resolution satellite LST obtained using the DMS approach tended to be underestimated, which tended to increase with higher LST values. The average RMSE values were 3.8 ºC, 3.1 ºC and 3.7 ºC, for DOY 199, 212 and 234, respectively. Figure 2b also shows that there was quite a strong correlation between RMSE and LST. This tended to increase as the LST values increased, with similar slopes for all three dates. In agreement with this, it is interesting to note that the LST values for DOY 234 were significantly lower than for the other two days. The good correlation between the sharpened and actual LST values shown in Figure 2b highlights the fact that sharpening the LST from 1 km to 20 m was able to resolve some of the spatial variability in surface temperature at the farm scale. However, its major underestimation shows that this sharpening was unable to capture the full range of temperatures registered under water-stressed conditions. This was probably due to the fact that the actual dynamic range of LST at 1 km was smaller than at 20 m and, in particular, when different rainfed and irrigated areas were present within the coarse-resolution pixel. The decision tree regression model was therefore trained using a set of temperatures with a narrower range than should normally be used to predict a larger range of temperatures. Figure 3 illustrates the temporal evolution of T (TSEB-PT S2+3 ) and T 0 (S-W) for well-watered ( Figure 3a) and stressed vines (Figure 3b) in the vineyard. It also shows comparisons with values modelled based on the vine water consumption model. The figure also includes the seasonal evolution of LAI and FAPAR estimated from Sentinel-2 using the biophysical processor and the modelled values for daily fiPAR (fiPAR d ) and the measured value for Ψ stem . Ψ fc and Ψ wilt respectively correspond to the stem water potential at field capacity and the lowest plant water potential at wilting point, when transpiration is null. Both these parameters are obtained from the vine water consumption model [58]. show seasonal measurements and simulations of stem and daily fiPAR, and estimates of leaf area index (LAI) and FAPAR based on the biophysical processor of S2. RMSD represents the root mean square deviation (mm season -1 ). ETref is the reference evapotranspiration obtained from the weather station using the Penman-Monteith equation. fc and wilt respectively correspond to the soil water potential at field capacity and to the lowest plant water potential at wilting point, when transpiration is null.

Comparison Between Methodologies
Estimates of actual evapotranspiration (ETa) and transpiration (T) obtained with TSEB-PTS2+3 showed a significant level of correlation with other methodologies. R 2 values of 0.83 and 0.65 (ETa) and of 0.78 and 0.70 (T) were obtained when they were compared with TSEB-PTairb and TSEB-2T, respectively ( Figure 4). Although the regression between TSEB-PTS2+3 and TSEB PTairb indicated a low RMSD, it tended to increase at lower ETa values (Figure 4a). This mostly occurred for DOY 199, which had the lowest evapotranspiration rates. In line with the previously mentioned underestimation of LST at high temperatures when using the TSEB-PTS2+3 method, this probably also caused an overestimation of the ETa values, particularly when the vines were subject to water stress. In fact, the sharpened S2+3 LST value was slightly lower for DOY 199 than DOY 212 (Figure 2). Estimates of T obtained using other approaches significantly declined in comparison with ETa. For instance, the averaged estimates of T for DOY 234 obtained using the TSEB-PTS2+3 approach accounted for 54% of The well-watered vines must have been transpiring at near to full potential rates. This was shown by the small differences between T 0 and T throughout the growing season and by the fact that the Ψ stem values were close to Ψ fc (Figure 3a). The seasonal average RMSDs for T and T 0 were 0.46 and 0.64 mm day −1 , respectively. Transpiration tended to increase sharply from DOY 100 (beginning of April) to 160 (mid-June), when maximum values of 4.2 mm day −1 for T 0 and 3.8 mm day −1 for T were registered. After DOY 230 (mid-August), transpiration began to moderately decline until senescence. When we compared the modelled and estimated values of T and T 0 between well-watered and stressed vines, as expected, we found significantly lower rates for the latter throughout the growing season. Although the RMSD values for T and T 0 were similar, the seasonal averaged R 2 values of the stressed vines were significantly lower (Figure 3b). The T values of stressed vines were slightly overestimated throughout the growing season, although the largest differences with respect to the modelled T values were detected during the period of highest water stress (from DOY 180 to 240). Figure 3b.2 shows how the Ψ stem values were lower during that particular period. This confirmed the hypothesis that the underestimation of sharpened satellite S2+3 LST for high LST values (Figure 2) would be reflected in an overestimation of T. Under both water regimes, the highest RMSD values were concentrated at the beginning and at the end of the growing season. This was probably because during these two periods, the cover cropping between rows tended to be higher and this was reflected in the higher estimates of satellite LAI. This was probably also linked to the fact that TSEB is only able to account for a single, actively transpiring canopy layer [59]. A similar consideration was reported by Semmens et al. [10]. Despite this, if we assume that most of the intercepted PAR was actually absorbed during the growing season, the R 2 between the estimates of FAPAR made with Sentinel-2 and the modelled fiPAR d would have been 0.92. Few studies have validated Sentinel-2 FAPAR products in heterogeneous crops such as grapevines. Although the correlation obtained probably depended on the training system, spacing distance and/or row orientation, and further research will be needed to adapt FAPAR estimates to these conditions, the preliminary results seem promising.

Comparison Between Methodologies
Estimates of actual evapotranspiration (ET a ) and transpiration (T) obtained with TSEB-PT S2+3 showed a significant level of correlation with other methodologies. R 2 values of 0.83 and 0.65 (ET a ) and of 0.78 and 0.70 (T) were obtained when they were compared with TSEB-PT airb and TSEB-2T, respectively ( Figure 4). Although the regression between TSEB-PT S2+3 and TSEB PT airb indicated a low RMSD, it tended to increase at lower ET a values (Figure 4a). This mostly occurred for DOY 199, which had the lowest evapotranspiration rates. In line with the previously mentioned underestimation of LST at high temperatures when using the TSEB-PT S2+3 method, this probably also caused an overestimation of the ET a values, particularly when the vines were subject to water stress. In fact, the sharpened S2+3 LST value was slightly lower for DOY 199 than DOY 212 (Figure 2). Estimates of T obtained using other approaches significantly declined in comparison with ET a . For instance, the averaged estimates of T for DOY 234 obtained using the TSEB-PT S2+3 approach accounted for 54% of the total ET a and for 68% and 62%, respectively, for DOY 199 and DOY 212 (Figure 4c). The T TSEB-PT S2+3 values were slightly overestimated with respect to TSEB-PT airb and TSEB-2T over their full range (Figure 4c,d). Furthermore, the T estimates obtained with TSEB-2T tended to be somewhat lower than those obtained with TSEB-PT airb (RMSD: 0.27 mm day −1 ). This overestimation of T could possibly be explained by the underestimation of LST which, in the case of the TSEB-PT S2+3 model, caused the iterative process to stop earlier at higher transpiration rates than it should have: TSEB-PT S2+3 barely reduced the initial α PT value from its original a priori value of 1.26. Figure 4 also shows that estimates of ET a and T made with TSEB-2T were lower than those obtained with the Priestley-Taylor (PT) approach over the full range of values, independently of whether S2+3 (Figure 4b) or airborne imagery (TSEB-PT airb , RMSD = 1.05 mm day −1 ) were used. In addition to the underestimation of LST with TSEB-PT S2+3 , another potential reason for this could have been the fact that when T c and T s are obtained from very high-resolution thermal images using a contextual approach, we are actually separating low-temperature pixels, corresponding to T c , from high-temperature pixels, corresponding to bare soil exposed to sunlight in the inter-row space of the vineyard. It is likely that in such a case, the estimated T s value would be higher than it should be, since the actual T s value should also include information relating to the shadowed soil and wet bulb surface because of the drip irrigation system below the vines. As TSEB-2T may have been overestimating T s in comparison with TSEB-PT S2+S3 and TSEB-PT airb (which assumed an effective T s value for the entire pixel, considering both the shaded and sunlit areas), the soil sensible heat flux (H s ) obtained was also higher and, as a consequence, both the E and ET a values would be lower than TSEB-PT S2+3 and TSEB-PT airb . This has previously been confirmed by the results obtained by Nieto et al. [27] for a Californian vineyard: they showed that TSEB-2T produced greater T/ET a and T s values than TSEB-PT when no active cover crop was present in the inter-row space.
vineyard. It is likely that in such a case, the estimated Ts value would be higher than it should be, since the actual Ts value should also include information relating to the shadowed soil and wet bulb surface because of the drip irrigation system below the vines. As TSEB-2T may have been overestimating Ts in comparison with TSEB-PTS2+S3 and TSEB-PTairb (which assumed an effective Ts value for the entire pixel, considering both the shaded and sunlit areas), the soil sensible heat flux (Hs) obtained was also higher and, as a consequence, both the E and ETa values would be lower than TSEB-PTS2+3 and TSEB-PTairb. This has previously been confirmed by the results obtained by Nieto et al. [27] for a Californian vineyard: they showed that TSEB-2T produced greater T/ETa and Ts values than TSEB-PT when no active cover crop was present in the inter-row space.

Regressions of ET a with Stem Water Potential
Stem water potential (Ψ stem ) is the most widely used plant water status indicator proposed for irrigation scheduling involving fruit trees and grapevines [60][61][62]. In this study, midday Ψ stem ranged from −0.5 to −1.3 MPa throughout the vineyard, indicating a high spatial variability in water status. However, the midday Ψ stem values for DOY 234 were not as low as for the other two days. Midday Ψ stem was significantly linearly related to ET a and T for all three dates ( Figure 5). The regressions indicated a clear tendency for both ET a and T to decrease with lower Ψ stem values. The maximum ET a values for vines transpiring at their potential rates (Ψ stem = −0.4 MPa) were found in TSEB-PT S2+3 , with values ranging from~3.9 to 4.4 mm day −1 in vines with a fiPAR of 0.4 and days with VPDs of 2.18 (DOY 212) and 2.49 kPa (DOY 199). Two other methods (TSEB-2T and TSEB-PT airb ) produced lower ET a and T estimates over the full range: this was because the sharpened S2+3 LST was also slightly underestimated at low values ( Figure 2). The slopes of the ET a TSEB-PT S2+3 vs. Ψ stem regressions were gentler than those obtained using other methods and revealed that the highest minimum ET a values (~2.6 mm day −1 ) corresponded to a Ψ stem of −1.3 MPa. On the other hand, the lowest ET a , corresponding to a Ψ stem = −1.3 MPa, was estimated using the contextual approach (~1.2 mm day −1 ). These differences in evapotranspiration rates, of about 1.4 mm day −1 for water-stressed vines, called into question the feasibility of using TSEB-PT S2+3 to detect vines under severe water-stressed conditions. Although the regressions were significant for all three dates, the y-intercept varied, which was probably due to differences in soil evaporation or atmospheric water demand. In all cases, DOY 234 produced the highest ET a values (probably due to the rainfall event) and DOY 199 produced the lowest when PT was used with either satellite or airborne imagery (TSEB-PT S2+3 and TSEB-PT airb ). obtained using other methods and revealed that the highest minimum ETa values (~2.6 mm day -1 ) corresponded to a stem of -1.3 MPa. On the other hand, the lowest ETa, corresponding to a stem=-1.3 MPa, was estimated using the contextual approach (~1.2 mm day -1 ). These differences in evapotranspiration rates, of about 1.4 mm day -1 for water-stressed vines, called into question the feasibility of using TSEB-PTS2+3 to detect vines under severe water-stressed conditions. Although the regressions were significant for all three dates, the y-intercept varied, which was probably due to differences in soil evaporation or atmospheric water demand. In all cases, DOY 234 produced the highest ETa values (probably due to the rainfall event) and DOY 199 produced the lowest when PT was used with either satellite or airborne imagery (TSEB-PTS2+3 and TSEB-PTairb). As previously explained, the ET a values estimated using TSEB-2T were clearly lower than those obtained using other methods. In addition, this method was only able to distinguish the higher ET a rates corresponding to DOY 234 due to the soil wetness and lower T s , but was not able to note the differences between days with lower soil evaporation rates (DOY 199 and 212). This confirms the previous hypothesis that the selection of bare soil exposed to sunlit pixels with high T s values, tends to overestimate the actual soil sensible heat fluxes (H s ) and that it is therefore difficult to detect differences in soil evaporation (E) between specific days. On the other hand, the regressions between T and Ψ stem showed similar patterns and non-significant differences were detected between days (Figure 5b,c,f). With the partial exception of DOY 234, the differences between ET a and T in TSEB-2T were very few: this confirmed the hypothesis that transpiration is a better indicator of crop water stress than bulk evapotranspiration, but it also suggested a significant relationship between the stem water potential and stomatal conductance in Tempranillo grapevines.

The Feasibility of Using the CWSI
Adopting regulated deficit irrigation (RDI) strategies is a common practice in grapevines in order to obtain high-quality berries. It is therefore necessary to have reliable information not only about T, but also about vine water status, in order to schedule irrigation. In this respect, the crop water stress index (CWSI) has been identified as a good crop water status indicator and proposed as a possible substitute for Ψ stem [22] due to the latter being both time-consuming and costly, particularly when used for scheduling irrigation in large vineyards. Figure 6 shows the regressions between Ψ stem and CWSI obtained with different methodologies. The results obtained produced significant correlations in all cases, with R 2 values of 0.29, 0.32 and 0.39, respectively, for CWSI-PT S2+3 , CWSI-PT airb and CWSI-2T (Figure 6a-c). As expected, the CWSI-PT S2+3 method had the lowest R 2 and produced a poor and non-linear response to the separation between well-watered and stressed vines (Figure 6a). The probable reason for the low CWSI corresponding to the vines in the vineyard which were water-stressed is that the sharpened S2+3 LST values tended to be underestimated, provoking the overestimation of T. Furthermore, shallower soil zones within the vineyard had much smaller soil water reservoirs and their vines were therefore more prone to suffer sustained water deficit. Such vines are likely to have a lower leaf area index (LAI) and this can be observed in Figure 7a. Under these conditions, it is also probable that vines will gradually modify both their hydraulic system and stomatal conductance to improve their acclimation to water deficits [63], with the net effect of modifying their LAI (i.e., their evaporative surface). Such adjustments are also driven to enable vines to maintain higher gas exchanges under conditions of low stem water potential [64]. These acclimation modifications could, at least in part, explain why the differences between T and T 0 in those zones where vines had low LAI were minimal and also why the CWSI values for low Ψ stem were not as high as expected. This could also partially explain the differences detected between the modelled and estimated T values in stressed vines (Figure 3b). Despite the low measured daily Ψ stem , vines growing under these circumstances (low LAI due to prolonged water stress conditions) probably still have high transpiration rates. As a result, in such cases, the modelled T values could also be a potential source of error.
When very high-resolution airborne imagery was used to calculate the averaged CWSI (CWSI-PT airb and CWSI-2T) of each 20 m pixel, it seems that there was a slight increase in the R 2 and CWSI values (Figure 6b,c) and, in particular, a more linear response of CWSI to Ψ stem in the case of CWSI-2T. In this case, the increase in CWSI was proportional for all levels of Ψ stem ; however, as a result, the estimated CWSI values for high Ψ stem were probably too high. This could probably be explained by the fact that the estimates of CWSI were obtained from a 20 by 20 m pixel (~20 vines), while those for Ψ stem were measured from three vines. The natural variability of the vines captured within a given pixel is likely to have a greater range of Ψ stem . This would tend to reduce the average Ψ stem and improve the regression with CWSI. In fact, the correlation between Ψ stem and the CWSI obtained using the empirical baseline approach (CWSI e ), which only used T c obtained from measured vines, produced the best R 2 value (R 2 = 0.55) and a steeper slope (Figure 6d).
Although the CWSI-Ψ stem regression using the Priestley-Taylor approach (CWSI-PT S2+3 and CWSI-PT airb ) did not produce the most desirable results in terms of absolute values, its regression was significant. This opens the door to the possibility of analyzing data in relative terms and calculating a relative CWSI (CWSI-PT S2+3,rel ), using either the highest T 0 value for the whole vineyard or irrigation management zone for each pixel. Although this methodology calculates the CWSI using the T 0 values of vines with different LAIs, it can also be used to obtain relative maps of CWSI using Sentinel imagery and distinguish zones with different levels of water status. In all our cases, the CWSI rel significantly improved the R 2 value, which increased from 0.29 to 0.61 for CWSI-PT S2+3 , from 0.32 to 0.56 for CWSI-PT airb and from 0.39 to 0.57 for CWSI-2T. In fact, Figure 7b compares the spatial distribution of the estimates of CWSI using two different approaches: (i) CWSI e resampled at 20 m, and (ii) CWSI-PT S2+3 , rel . Our results indicated a similar spatio-temporal pattern between the two methods, with averaged RMSE values of 0.07, 0.12 and 0.12 for DOY 199, 212 and 234, respectively. These results support the argument that CWSI-PT S2+3,rel maps could be used as tools for mapping spatial variability in water status when no other thermal infrared data sources, such as airborne or Landsat imagery, are available.
Remote Sens. 2020, 12, 2299 18 of 25 improve the regression with CWSI. In fact, the correlation between stem and the CWSI obtained using the empirical baseline approach (CWSIe), which only used Tc obtained from measured vines, produced the best R 2 value (R 2 =0.55) and a steeper slope (Figure 6d). Although the CWSI-stem regression using the Priestley-Taylor approach (CWSI-PTS2+3 and CWSI-PTairb) did not produce the most desirable results in terms of absolute values, its regression was significant. This opens the door to the possibility of analyzing data in relative terms and calculating a relative CWSI (CWSI-PTS2+3,rel), using either the highest T0 value for the whole vineyard or irrigation management zone for each pixel. Although this methodology calculates the CWSI using the T0 values of vines with different LAIs, it can also be used to obtain relative maps of CWSI using Sentinel imagery and distinguish zones with different levels of water status. In all our cases, the CWSIrel distribution of the estimates of CWSI using two different approaches: i) CWSIe resampled at 20 m, and ii) CWSI-PTS2+3,rel. Our results indicated a similar spatio-temporal pattern between the two methods, with averaged RMSE values of 0.07, 0.12 and 0.12 for DOY 199, 212 and 234, respectively. These results support the argument that CWSI-PTS2+3,rel maps could be used as tools for mapping spatial variability in water status when no other thermal infrared data sources, such as airborne or Landsat imagery, are available.

The Spatial Distribution of Biophysical and ET Parameters Within the Vineyard
Maps of the spatial distribution of the biophysical parameters of the vegetation and ET components in the vineyard are shown in Figure 8. Overall, there was a positive spatial agreement between the LAI or FAPAR and the ET components. One of the main advantages of TSEB is its ability to obtain separate soil evaporation and canopy transpiration values. As a result, this vineyard had a seasonal accumulated ETa of 470 ± 26 mm, 32% of which was lost through evaporation (150 ± 26 mm) and 68% due to transpiration (320 ± 33 mm). These values agreed with those presented in other studies conducted in drip-irrigated vineyards which reported E/ETa ratios of ~0.3±0.12 [26]. The

The Spatial Distribution of Biophysical and ET Parameters Within the Vineyard
Maps of the spatial distribution of the biophysical parameters of the vegetation and ET components in the vineyard are shown in Figure 8. Overall, there was a positive spatial agreement between the LAI or FAPAR and the ET components. One of the main advantages of TSEB is its ability to obtain separate soil evaporation and canopy transpiration values. As a result, this vineyard had a seasonal accumulated ET a of 470 ± 26 mm, 32% of which was lost through evaporation (150 ± 26 mm) and 68% due to transpiration (320 ± 33 mm). These values agreed with those presented in other studies conducted in drip-irrigated vineyards which reported E/ET a ratios of~0.3 ± 0.12 [26]. The partitioning of E/T ranged from 0.64 (in areas of the vineyard with low LAI and high evaporation) to 0.33 (in areas with high LAI and low evaporation). The total amount of water applied to the vineyard through irrigation (280 mm) and rainfall (214 mm) accounted for a total of 494 mm over the whole growing season (1 April to 30 September). Considering that the differences in annual ET a across the vineyard varied from 397 to 576 mm, it is therefore clear that some areas of the vineyard were over-irrigated, with water being lost through evaporation or runoff, while others were subject to water stress during certain periods of the growing season. Finally, the seasonal averaged CWSI-PT S2+3,rel showed that Remote Sens. 2020, 12, 2299 20 of 25 the maximum values for the stressed vines was 0.43 and that these were located in the areas with the lowest LAI or FAPAR values, which also coincided with the lowest Ψ stem .
through irrigation (280 mm) and rainfall (214 mm) accounted for a total of 494 mm over the whole growing season (1 st April to 30 th September). Considering that the differences in annual ETa across the vineyard varied from 397 to 576 mm, it is therefore clear that some areas of the vineyard were over-irrigated, with water being lost through evaporation or runoff, while others were subject to water stress during certain periods of the growing season. Finally, the seasonal averaged CWSI-PTS2+3,rel showed that the maximum values for the stressed vines was 0.43 and that these were located in the areas with the lowest LAI or FAPAR values, which also coincided with the lowest stem.

Conclusions and Future Perspectives
This study demonstrates that the use of TSEB with sharpened LST obtained from the Sentinel-3 and Sentinel-2 satellites and the S-W model was adequate for estimating both actual and potential transpiration from grapevines throughout a growing season. However, the T values estimated with TSEB-PTS2+3 tended to be overestimated, particularly at low transpiration rates. In contrast, the CWSI-PTS2+3 values obtained for vines under severe water stress were excessively low. This study also explains that the main reason for this was that the data mining sharpening (DMS) approach used to estimate high-resolution LST in TSEB exhibited certain limitations in its ability to capture the full range of temperatures and therefore to identify crops suffering water stress, particularly where water stress did not reduce the amount of biomass. A regression between the airborne LST and sharpened S2+3 LST values demonstrated that the latter tended to be underestimated. In this context, it is possible that the ET products derived by SEN-ET [65] (TSEB-PTS2+3) are still not valid for use in applications in which regulated deficit irrigation (RDI) strategies have to be adopted in almost realtime. This could probably, however, be a good option when these data are assimilated into crop models which either prescribe irrigation recommendations or evaluate water consumption at the regional level over medium-or long-term periods. The lack of high-resolution TIR remains an important issue at the operational setting.
Of the three methodologies tested (TSEB-PTS2+3, TSEB-PTairb and TSEB-2T), TSEB-2T was the one that appeared to perform the best. This was shown by the higher R 2 and slightly lower RMSD when

Conclusions and Future Perspectives
This study demonstrates that the use of TSEB with sharpened LST obtained from the Sentinel-3 and Sentinel-2 satellites and the S-W model was adequate for estimating both actual and potential transpiration from grapevines throughout a growing season. However, the T values estimated with TSEB-PT S2+3 tended to be overestimated, particularly at low transpiration rates. In contrast, the CWSI-PT S2+3 values obtained for vines under severe water stress were excessively low. This study also explains that the main reason for this was that the data mining sharpening (DMS) approach used to estimate high-resolution LST in TSEB exhibited certain limitations in its ability to capture the full range of temperatures and therefore to identify crops suffering water stress, particularly where water stress did not reduce the amount of biomass. A regression between the airborne LST and sharpened S2+3 LST values demonstrated that the latter tended to be underestimated. In this context, it is possible that the ET products derived by SEN-ET [65] (TSEB-PT S2+3 ) are still not valid for use in applications in which regulated deficit irrigation (RDI) strategies have to be adopted in almost real-time. This could probably, however, be a good option when these data are assimilated into crop models which either prescribe irrigation recommendations or evaluate water consumption at the regional level over medium-or long-term periods. The lack of high-resolution TIR remains an important issue at the operational setting.
Of the three methodologies tested (TSEB-PT S2+3 , TSEB-PT airb and TSEB-2T), TSEB-2T was the one that appeared to perform the best. This was shown by the higher R 2 and slightly lower RMSD when it was evaluated against modelled T. The regressions between T and CWSI-2T with stem water potential (Ψ stem ) also produced the highest R 2 values. The only constraint detected in TSEB-2T was that the selection of pixels corresponding to T s was related to bare soil that was well-exposed to sunlight in the inter-row space of the vineyard. This resulted in an overestimation of the soil sensible heat flux (H S ) and therefore a reduction in estimated soil evaporation in comparison with other methods which used the Priestley-Taylor approach, either with satellite or airborne LST. Further studies should focus on the evaluation of TSEB where vegetation and the soil temperature of shaded and unshaded areas and wet bulb surfaces created by drip irrigation are taken into account.
With regard to the use of CWSI as a plant water status indicator, the baseline empirical approach seemed to be the method which provided the best results: this was indicated by the highest R 2 value obtained in the regression with Ψ stem . When very high-resolution thermal imagery is available, this study suggests that CWSI e followed by CWSI-2T are probably the best approaches to use in order to detect differences in water status affecting grapevines. Even so, as the use of airborne thermal imagery is still far from becoming an operational reality, immediate improvements should focus on using LST satellite imagery.
This study demonstrates that to relativize CWSI-PT S2+3 , a feasible option could be to emphasize the significant spatial differences in water status that it can detect. However, it is also true that in some cases, this approach could exacerbate water stress in vines which are relatively well-watered. As this study demonstrates that TSEB-PT airb significatively improved estimates of T and CWSI compared with those made with TSEB-PT S2+3 , it is realistic to suggest that the Priestley-Taylor (PT) approach, combined with certain improvements to its parametrization of potential transpiration estimates, could potentially provide a good alternative. This could perhaps be tested in some other on-going or future satellite missions incorporating higher spatio-temporal resolution TIR imaging, such as Landsat 9, ESA-LSTM, ECOSTRESS [66] or TRISHINA [67].
There is interest in testing this methodology with other crops and environments in order to assess whether the results found in this study are repeated in other conditions. Several current works are also focused on modifying the decision tree analysis used in the DMS approach by giving greater strength to the SWIR bands in order to capture canopy water content and by merging optical and microwave SAR data. Furthermore, as one of the main requirements for obtaining accurate estimates of T c and T s is to have a good parametrization of transpiration under potential conditions, research is currently being undertaken to improve estimates of T 0 , considering different phenological stomatal responses to vapor pressure. This could also help to promote a better characterization of crop water needs. This would be achieved through the direct incorporation of satellite-derived biophysical products (such as LAI or albedo) into energy combination equations, such as those of Penman and Monteith and Shuttleworth and Wallace, and by avoiding the use of more controversial empirical crop coefficients to define crop potential ET. Funding: This study was partially supported by the Catalan Agency for Business Competitiveness (ACCIO), using FEDER funds to carry out the RIS3CAT project on Low Input Sustainable Agriculture (LISA) (COMRDI16-1-00031-05), and by the Horizon 2020 Programme for Research and Innovation (H2020) of the European Commission, in the context of the Marie Sklodowska-Curie Research and Innovation Staff Exchange (RISE) action and ACCWA project: grant agreement no.: 823965. This study has also benefitted from interaction with the Sentinels for Evapotranspiration project, which is funded by the European Space Agency (ESA), in which several evapotranspiration products have been developed using Sentinel imagery.