*2.1. Study Site*

The study site was a commercial vineyard located in Ripperdan, Madera County, California, USA (Figure 2; 36◦50024.3900 N , 120◦12043.2900 W). The analysis evaluated a 31.23 hectare block, consisting of Chardonnay vines (*Vitis vinifera*) and planted in 2009 on 1103P rootstock and double vertical canopy trellising. Primary viticultural techniques included use of drip irrigation and annual cover-cropping with oats (45%), barley (45%), and mustard (10%). Site soils consisted of Dinuba-El Peco alkaline fine sandy loam for the majority of the vineyard block, and a Grangeville sandy loam distributed in what appears to be remnant stream bed coursing through the northwest corner. Data were collected on 10 site visits for the period of May–September 2018. Across the growing season, the air temperature ranged 5.2–40.1 ◦C, minimum to maximum daily measured hourly, and resulted in 2058 degree days. For the calendar year 2018, the site received 733 mm of rainfall. Site observations and agronomic records indicate uniformity of vine canopy growth, and hence soil exposure, for each date of site visitation and across the growing season.

**Figure 2.** Site map of Ripperdan 760 Block in Madera County, California, USA.

## *2.2. Analytical Methods*

#### 2.2.1. Reference Evapotranspiration

Reference evapotranspiration describes the fraction of ET from a reference crop given measured weather conditions. In the METRIC energy balance method, reference ET (*ETref*) is typically based on the ratio of irrigated alfalfa, the reference crop, to *ET<sup>o</sup>* for grass. *ETref* was collected from Spatial CIMIS [40], to obtain daily *ET<sup>o</sup>* for the study vineyard throughout the growing season. Spatial CIMIS renders *ET<sup>o</sup>* at a 2-km spatial resolution by coupling daily solar radiation from NOAA GOES satellite with the Heliosat-II method [41] and interpolating air temperature *T*, relative humidity *RH*, and wind speed *u<sup>2</sup>* from CIMIS stations. *ET<sup>o</sup>* is calculated using a modified ASCE–Penman Montieth equation for clipped grass reference from [42] (Equation (1)):

$$ET\_o = \frac{0.408\Delta (R\_{\rm nl} - G) + \gamma \frac{900}{T + 273} \mu\_2 (e\_{\rm s} - e\_{\rm a})}{\Delta + \gamma (1 + 0.34 \mu\_2)} \tag{1}$$

where *ET<sup>o</sup>* is reference crop evapotranspiration (mm day−<sup>1</sup> , *R<sup>n</sup>* is net radiation at the crop surface (MJ m−<sup>2</sup> day−<sup>1</sup> , *G* is soil heat flux density at the soil surface (MJ m−<sup>2</sup> day−<sup>1</sup> , *T* is the daily mean air temperature at 2 m (ms−<sup>1</sup> , *u*<sup>2</sup> is wind speed at 2 m, the difference *es*–*e<sup>a</sup>* is the saturation vapor pressure deficit (kPa), *γ* is the psychrometric constant (∼66 Pa K−<sup>1</sup> , and ∆ is the slope of the saturation vapor pressure-temperature curve (kPa ◦C −1 ).

#### 2.2.2. Earth Engine Flux

Earth Engine Flux (EEFlux) is an implementation of METRIC [29] in Google Earth Engine for Landsat or MODIS imagery with automatic calibration using groundlevel weather station data, land-use data, and a digital elevation model that calculates the residual latent energy used to evaporate water. EEFlux uses GridMET weather data to calculate daily *ETref* with the Standardized ASCE-PM equation. EEFlux allows users to select a region of interest, range of dates, and satellite overpasses to calculate ET using the METRIC surface energy balance method in a matter of minutes. We obtained 11 EEFlux Landsat 8 daily ET products from 16 April to 23 September for the 2018 growing season.

#### 2.2.3. OpenET

OpenET (https://etdata.org/) is a collection of Google Earth Engine-based ET models. OpenET is currently in beta with some features accessible on Github, including SSEBop and an NDVI-based ET model. Both NDVI-ET and SSEBop models are run via Google's Earth Engine Python API (https://github.com/google/earthengine-api). While NDVI can be used as a proxy for photosynthetic activity, it has also been shown to vary linearly with *ETo*. Daily *OpenETaNDV I* rasters were used to produce a dimensionless ET fraction (*ET<sup>f</sup>* ) raster by using a linear model adjusted for *ETo*. NDVI was calculated using the Red and NIR bands from each platform (see Figure 1 and Equation (2)) (S2: Bands 4 (650–680 nm) and 8 (785–899 nm); sUAS: Red (663–673 nm) and NIR (820–860 nm)). A slope coefficient of 1.25 and intercept of 0 were set to conform to OpenET NDVI-ET-Beta model parameters [31]. Daily ET rasters were created for both models (NDVI-ET and SSEBop) and used for model comparison. EOS imagery was processed through Earth Engine. sUAS imagery was processed in ArcGIS Pro 2.4 following OpenET's NDVI-ET method.

The linear transfer function for *OpenETaNDV I* is as follows:

$$OpenET\_{aNDVI} = 1.25(\frac{\sigma\_{\text{NIR}} - \sigma\_{\text{RED}}}{\sigma\_{\text{NIR}} + \sigma\_{\text{RED}}})ET\_o \tag{2}$$

where *σ*NIR is reflectance in the near infrared band and *σ*RED is reflectance in the red band (see Figure 1). *ET<sup>o</sup>* was derived from from Spatial CIMIS, such that *ET<sup>a</sup>* is actual ET (mm day−<sup>1</sup> ) and *ET<sup>o</sup>* is daily reference ET from Spatial CIMIS (mm day−<sup>1</sup> ).

#### *2.3. Data Products*

#### 2.3.1. Satellite Imagery

Sentinel 2 imagery was collected from USGS EarthExplorer for each coinciding overpass with a field survey and overpasses within 3 days of a Landsat overpass. Level 1C Top of Atmosphere products were run through the Sen2Cor plugin (ESA, v.2.5.5) in SNAP Toolbox (ESA, v.6.0). Sen2Cor was used to correct for any atmospheric effects and to produce Level 2A Bottom of Atmosphere reflectance products. Corrected Sentinel data were transformed from JPEG2000 to GeoTiff format using GDAL, and subsequently analyzed in ArcGIS and R for cross-platform comparisons. No attempt was made to harmonize Landsat and Sentinel imagery vis-à-vis Claverie et al. [43] given the limited availability of such data to grower practitioners.

#### 2.3.2. sUAS Imagery and Ancillary Data

Data Acquisition Flights: All sUAS missions were planned with Ardupilot Mission Planner to maintain a consistent flight pattern over the 10 acquisition dates between April and September 2018. Each flight was operated under FAA Part 107 flight rules by a certified Remote Pilot. The early season flights in April and May were completed using the Finwing Sabre; a foam fixed-wing sUAS, carrying a Parrot Sequioa sensor. Subsequent flights were flown with a DJI S1000 multirotor sUAS carrying a Micasense Red Edge-M sensor. Multirotor flights (DJI S1000) were flown in a simple grid pattern with

transects in a north-south orientation (90 m Altitude, 80% Sidelap, 80% Frontlap, 9 m s−<sup>1</sup> speed, 20 m flightline spacing). Most dates required four flights to cover the entire vineyard block at the desired altitude, front and side overlap. Additional flights were performed in case of sensor drop-outs or wind reducing flight endurance. Flights were completed as close to solar noon to minimize shadows.

Canopy Volume: We used a voxel approach to calculate the volume of the canopy for each of the RSE flight over the growing season. Canopy volume was calculated for each flight date with a densified point cloud from Pix4D. The densified point cloud was then classified in CloudCompare using the Cloth Simulation Filter [44]. Once ground and non-ground points were identified, a 25 cm Digital Terrain Model (DTM) was made from using only ground points and a 25 cm Digital Surface Model (DSM) was made from using all points to model the canopy as well as the bare soil elevation. The DTM was subtracted from the DSM to produce a Canopy Height Model (CHM) for each applicable date. The CHM was then multiplied by the area of each cell to calculate the canopy volume of the full field. The CHM was also used to create an above ground canopy mask used to exclude bare soil and cover crop in subsequent analyses.

Reflectance Imagery: Reflectance images from the Parrot Sequoia and Micasense Red Edge-M were georeferenced and stitched through the photogrammetry software Pix4DMapper Pix4D, version 4.3.3, 2018-04-09). RTK GPS-surveyed Ground Control Points (GCPs) were used for each flight to ensure accurate orthomosaic geolocation. Our sUAS ET model used orthomosaics of reflectance values and DSMs produced for each flight date. Model results were upscaled to 10 m spatial resolution to conform with Sentinel 2 pixel geometry. The upscaled outputs included mean, maximum, minimum, and standard deviation of the finer sUAS pixel values within each Sentinel 2 pixel.

#### **3. Results**

#### *3.1. Irrigation Delivery and Canopy Growth*

Observed irrigation practices for the growing season resulted in a average of 5.26 mm day−<sup>1</sup> across the growing season and 9.96 mm day−<sup>1</sup> for irrigation days, but varied considerably across the growing season. The longest period without irrigation extended for 19 days (13 May to 31 May), compared to a maximum application of 28.07 mm on 8-May (Figure 3). In total, 264,338 m<sup>3</sup> of water was applied to the 31.23 ha vineyard block over the course of the growing season April–September.

**Figure 3.** Daily air temperature ◦C (right axis) for the 2018 growing season is shown for mid-day (red line) as well as maxima-minima (gray band). Applied irrigation for each day in mm (left axis) is shown as blue bars.

#### *3.2. Remote Sensing of Evapotranspiration*

For this study, we focused on the comparative utility of EOS and sUAS ET estimates over a single growing season (2018). We flew sUAS on ten dates April–September (Table 1), and calculated ET with EOS imagery on 20 dates (Table 2). We observe the change in daily ET over the course of the growing season (Figure 4), peaking in mid-July. During the early portion of the growing season April–June there is a close tracking relationship between NDVI-ET and EEFlux estimates, rising from 1.4 mm day−<sup>1</sup> on 16 April to approximately 8 mm day−<sup>1</sup> 06 June. EOS observations diverge in July–September, where EEFlux estimates peak 06 August at 8.26 mm day−<sup>1</sup> and remain high while NDVI-ET values drop after harvest. The effect of canopy masking on sUAS NDVI-ET estimates resulted in a mean positive difference of 1.23 mm day−<sup>1</sup> ET, presumably due to fewer mixed pixels. We observe the smallest difference between masked and unmasked sUAS NDVI-ET on 19 September with a mean difference of 0.49 mm day−<sup>1</sup> . The greatest deviation occurred on 16 July with a mean difference of 1.75 mm day−<sup>1</sup> . Masked sUAS NDVI-ET and EOS S2 NDVI-ET estimates were often the most similar between all methods. During the season EOS ET ranged from a minimum pixel value of 0.93 mm day−<sup>1</sup> on 16 April to the maximum 9.74 mm day−<sup>1</sup> on 21 June. The minimum observed EOS mean daily ET of 1.45 mm day−<sup>1</sup> occurred on 16 April using L8 EEFlux. The maximum observed EOS mean daily ET of 8.14 mm day−<sup>1</sup> occurred on 06 August (Table 2). sUAS NDVI-ET ranged from 0.49 mm day−<sup>1</sup> on 07 May to 8.79 mm day−<sup>1</sup> on 16 July. Mean sUAS NDVI-ET was lowest on 19 September at 3.58 mm day−<sup>1</sup> and greatest on 16 July at 8.04 mm day−<sup>1</sup> .

**Table 1.** Flight dates and platforms.


**Figure 4.** Daily Evapotranspiration (ET) in mm (x-axis) ridgeline density plot for each overpass date (y-axis) comparing each method used. Landsat 8 EEFlux is colored red, Sentinel 2 NDVI-ET is shown in orange, sUAS unmasked NDVI-ET is light blue, and sUAS masked NDVI-ET is green.


**Table 2.** Summary of ET (mm) estimates from sUAS and EOS observation for 2018.

Large irrigation events occurred on 22–23 September with a cumulative applied water of approximately 52 mm for those two days (Figure 3) coinciding with the high reported ET from EEFlux on 23 September. On 19 September observations by sUAS and EOS are much lower than observed on 23 September from EEFlux. On 19 September it had been eight days since last applied water with a noon temperature of 24.7 ◦C. Similarly, 23 September had a noon temperature over 5 ◦greater at 30.5 ◦C.

Seasonal Evapotranspiration

Mean ET<sup>a</sup> from Landsat 8 EEFlux peaks on 25 August at 7.41 mm day−<sup>1</sup> (Figure 5). The other RSE methods peak much earlier in the season: Sentinel 2 NDVI-ET peaks on 27 June at 7.74 mm day−<sup>1</sup> ; sUAS Masked NDVI-ET peaks on 7 July at 7.71 mm day−<sup>1</sup> ; and sUAS Unmasked NDVI-ET peaks lower than the other methods on 3 June at 6.32 mm day−<sup>1</sup> . Landsat 8 EEFlux and sUAS Unmasked NDVI-ET follow a similar trend until the end of June where EEFlux values rise to above 7 mm day−<sup>1</sup> and remain there for the remainder of the season. The cumulative ET<sup>a</sup> over the growing season is shown in Figure 6. Cumulative ET<sup>a</sup> from EEFlux and sUAS Unmasked NDVI-ET are nearly the same until mid-July when EEFlux quickly rises and overtakes all other methods for a season total ET<sup>a</sup> of 313,634 m<sup>3</sup> (Table 3). The unmasked sUAS ET<sup>a</sup> was within 5000 m<sup>3</sup> of the volume of applied irrigation and was 95 percent confident within a range of 50,000 m<sup>3</sup> . In contrast both EOS methods had about double the range of confidence and the masked sUAS ET<sup>a</sup> was between with a range of about 80,000 m<sup>3</sup> .

**Figure 5.** Seasonal Daily ET<sup>a</sup> in mm (y-axis) interpolated from loess model fit. Landsat 8 EEFlux is colored purple, Sentinel 2 NDVI-ET is shown in green, sUAS masked NDVI-ET is blue, and sUAS unmasked NDVI-ET is red.

**Figure 6.** Cumulative ET<sup>a</sup> in m<sup>3</sup> (y-axis) interpolated from loess model fit. Landsat 8 EEFlux is colored purple, Sentinel 2 NDVI-ET is shown in green, sUAS masked NDVI-ET is blue, and sUAS unmasked NDVI-ET is red. Cumulative Applied Water from irrigation is shown in light blue columns.


**Table 3.** Cumulative ET (10<sup>3</sup> m<sup>3</sup> ) with 95% confidence interval bounds, and applied irrigation, for 2018 growing season.

#### **4. Discussion**

#### *4.1. Viticultural Considerations*

The application of RSE in viticulture is important for several reasons. Viticulture, the cultivation of wine and table grapes, is a global enterprise with increasing acreage globally [45]. It is increasingly challenged by changing hydroclimates and atmospheric warming. Impacts of warming atmospheric temperatures not only have the potential to affect ET in vineyards, but also impact fruit quality [46]. Observed air temperatures have increased markedly across major winegrowing regions globally [47], and winegrape growers have started to adopt adaptation strategies aimed at minimizing the harshest impacts of elevated temperature [48], namely dehydration and early ripening. One adaptation strategy is to reduce canopy temperatures through evaporative cooling from overhead misting, which would further increase total crop water use due to this additional application. Given water scarcity in many winegrape growing regions, minimizing the water footprint of viticulture is a top priority for the industry [49,50]. As such, grower-ready tools are necessary to better manage irrigation strategies. While subscription services to provide guidance on daily crop water demand are increasingly available, some growers may choose to use satellite or sUAS monitoring protocols as to control their own data and/or perform custom, on-demand analysis. RSE of vineyards has been the focus of several recent studies [51], in part due to the high commercial value of the crop and increasing sustainability practices [52]. In this study, we examined the temporal distribution of an energy balance EOS satellite based RSE approach to index-based RSE methods, including commercial-grade sUAS platforms that growers may want to employ. Our results indicate that there is general agreement between between EOS and sUAS NDVI-ET estimates and, moreover, that masking sUAS imagery to canopy pixels reduced the variance of RSE values (Figure 6) and generally conformed to the amount of applied water within the vineyard block (Figure 3). The processing of sUAS imagery is not trivial (see Montazar et al. [53] for another example), and smaller grower operations may be challenged to implement sUAS RSE on a routine basis.

#### *4.2. RSE Considerations*

The offset nature of EOS overpasses prevented a direct comparison between RSE products, and daily ET estimates especially given external factors such as weather and irrigation. Across most of the growing season we observe a similar trend between the EEFlux and NDVI-ET methods with low ET early in the season, peaking around July, then dropping back down again after harvest. Late in the season, however, EEFlux estimates remain consistent with early/mid summer values rather than a decrease to coincide with leaf senescence. We expect that later in the season, EEFlux is overestimating the actual ET due to recent irrigation applied on 22–23 September in conjunction with days 5 ◦C greater than 19 September during the S2 overpass and sUAS flights (Figure 3). Masking sUAS imagery yielded similar results to the EOS methods and appears to overestimate the actual amount of ET. The unmasked sUAS imagery yielded a lower estimate that was very close to the amount of applied water during the season.

As expected, soil masking increased the mean ET and decreased variance since isolating the grape canopy removed ET influence from bare soil and cover crop in the image (Table 2). We also found that the difference between masked and unmasked sUAS ET changes temporally with the largest difference occurring during peak ET in mid-July and the smallest difference after harvest. Assuming the masked sUAS ET is more representative of the real *ETa*, it implies that lower spatial resolution (i.e., mixed pixel) EOS RSE values are underestimating during the period of peak ET, and moreoever the EEFlux values are compromised during the end of season period that experience high air temperature, full field saturation from irrigation, and vine senescence. These observations are consistent with other studies, such as Gowda et al. [54] that showed high errors EOS RSE of post-harvest cotton, and more generally the theoretical basis by Li et al. [55] showing high estimates of *ET<sup>a</sup>* under saturated soil conditions and sensed high surface temperatures. We also observed large variations between the minimum and maximum RSE values (Table 2), with the lowest values attributed to the unmasked sUAS and EEFlux consistently had very high minimum values (2–3 mm day−<sup>1</sup> higher than S2 NDVI-ET estimates). The lack of variance in EEFlux estimates could have consequences for measures of centrality and confidence for growing season totals (Table 3).

#### *4.3. Land Use Considerations*

This study focused on the use of commercial grade "off-the-shelf" approaches to sUAS RSE for the purposes of improving irrigation regimes in vineyards. The study location in California's San Joaquin Valley has a number of concomitant land use stressors that require remediation. As new RSE and similar technologies become more readily available and deployable, it is highly likely that these technologies can be used to mitigate issues beyond water scarcity. Current over-reliance on groundwater reserves for crop applied water has stressed California's aquifers [56]. Furthermore, other forms of agriculture have generated excess nitrogen throughout the region [57,58], which can lead to myriad impacts on environmental [59] and agroeconomic systems [60], and ultimately human well-being. Remedies in this region include enhancing viticulture in active recharge zones around groundwater dependent communities [61], constructing distributed desalination systems for brackish irrigation return water [62], incorporating wildlife friendly farming practices [63], and integrated management of forests and watersheds for clean water supply [64]. Despite current challenges in the operational deployment of sUAS [65] for environmental management, sUAS remote sensing can provide critical data to inform such land use management decisions beyond better estimates of ET, including biomass and carbon quantification for regenerative agriculture [66]. In this sense, strategic land use planning and decision making in the agricultural sector would benefit from broader adoption of sUAS technologies.

#### *4.4. Study Limitations*

A limitation of the NDVI-ET approach used here, aside from its comparative simplicity, is that it relies on external sources for ETo. While there was a nearby weather station with a similar climate and land use, the radiometer was unreliable so reference ET from a local site was unavailable. Therefore, we relied on the Spatial CIMIS daily ET<sup>o</sup> product since the nearest CIMIS stations were over 30 km away and were potentially unrepresentative of the ET<sup>o</sup> at the study site.

EOS ET estimation is limited by infrequent return times and low spatial resolution; which can be improved through the use of sUAS. Seasonal ET estimates can be heavily biased by outlying overpass dates where conditions may not represent the majority of days between the overpasses, like we observed towards the end of the season with Landsat 8. We also see that EOS estimates overestimate the actual ET for our study site. Using fine resolution sUAS imagery up-scaled to S2 resolution (10 m) reduced the mean ET<sup>a</sup> value by an average 1.233 mm for the season which accounts for a substantial volume for a large area. The EOS ET estimates are most similar to the sUAS ET when bare soil is masked out, especially during the peak of the season when the canopy is fully grown and covers most of the field area. It seems likely that the lower spatial resolution from EOS platforms leads to an overestimation of ET for vineyards as the peak ET from the canopy gets generalized for the entire cell area, which in reality contains areas with much less ET like we can observe with the sUAS imagery. Furthermore, future studies would benefit from harmonized EOS data reconciling differences in sensor sensitivity and spectral band widths [39,43].
