sUAS Remote Sensing of Vineyard Evapotranspiration Quantifies Spatiotemporal Uncertainty in Satellite-Borne ET Estimates

Small Unmanned Aerial Systems (sUAS) show promise in being able to collect high resolution spatiotemporal data over small extents. Use of such remote sensing platforms also show promise for quantifying uncertainty in more ubiquitous Earth Observation System (EOS) data, such as evapotranspiration and consumptive use of water in agricultural systems. This study compares measurements of evapotranspiration (ET) from a commercial vineyard in California using data collected from sUAS and EOS sources for 10 events over a growing season using multiple ET estimation methods. Results indicate that sUAS ET estimates that include non-canopy pixels are generally lower on average than EOS methods by >0.5 mm day−1. sUAS ET estimates that mask out non-canopy pixels are generally higher than EOS methods by <0.5 mm day−1. Masked sUAS ET estimates are less variable than unmasked sUAS and EOS ET estimates. This study indicates that limited deployment of sUAS can provide important estimates of uncertainty in EOS ET estimations for larger areas and to also improve irrigation management at a local scale.


Introduction
Global environmental change and anthropogenic activity have long stressed Earth's hydrological cycle [1,2]. Evapotranspiration remains the least certain quantified component of the hydrological cycle [3], with implications for not only water resources planning and management, but also for human livelihoods and supporting ecosystems [4]. Globally, irrigated agriculture represents 70% of water withdrawn from surface and ground water supplies but estimates of actual consumptive loss (the largest component of the water balance in an agricultural region) through evapotranspiration (ET) remain coarse [5]. New advances in satellite remote sensing of ET via Earth Observation Systems (EOS) show promise in providing consistent and reliable quantitative estimates with global coverage and reasonable repeat cycles. However, the spatiotemporal resolution of EOS platforms and sensors may obscure finer spatial resolution phenomena or weather-related aberrations observed at finer time scales. Small Unmanned Aerial Systems (sUAS) operated under the Federal Aviation Administration (FAA) Part 107 licensing outfitted with optical sensors are increasingly used in vegetation remote sensing [6] and precision agriculture applications [7]. More recently, however, comparative studies of

Remote Sensing of Evapotranspiration
Evapotranspiration has long been the focus of hydrological and agrometeorological studies [23,24]. Zhang et al. [25] reviewed the state of science in remote sensing of evapotranspiration (RSE), described its use to estimate and map ET at regional to continental scales, and highlighted existing major RSE estimation methods. Because so many studies and reviews have been conducted on ET measurements, modeling and RSE methods of retrieval, they are not repeated here. Rather, we focus on the fact that RSE from EOS provides relatively frequent, consistent, and spatially contiguous measurements for global estimation, monitoring, and mapping of ET flux; however, due to the relatively coarse granularity of such data at the field level, we explore how high spatiotemporal resolution data from sUAS could quantify inherent variability in such estimates.
In almost all cases, RSE relies on energy balance methods. Most RSE methods focused on energy balance are rooted in Surface Energy Balance Algorithm for Land (SEBAL) [26][27][28] an image-processing model to quantify surface energy balance components, at both local and regional scales using empirical relationships and physical parameterization. SEBAL requires digital imagery data collected by any satellite sensor measuring visible, near-infrared, and thermal infrared radiation, and their derivative products including surface temperature, normalized difference vegetation index (NDVI), and albedo. SEBAL was a methodological precursor to METRIC (Mapping EvapoTranspiration at high Resolution with Internalized Calibration), developed by Allen et al. [29], which is now a standard remote sensing estimation approach using Landsat 5, 7 and 8 EOS imagery (including thermal bands), local weather station data, and calibrated from either nearby alfalfa or pasture, and bare soil land cover types. Landsat-based METRIC, and other similar approaches are described in Irmak et al. [30].
More recently, a number of open-source methods have emerged as Google Earth Engine-based collaborative [31]. One such method includes EEFlux, a version of the METRIC model as developed by Morton et al. [31] employing algorithms by Irmak et al. [30]; EEFlux has been used because of expansive access to EOS imagery and computational assets, including built-in calibration using reference ET from the North America Land Data Assimilation System [32,33] and GridMET [34] for the contiguous United States. In implementation, however, RSE methods such as these are explicitly reliant upon thermal infrared sensors (λ = 10,500-12,510 nm) onboard EOS platforms, which largely prohibits inclusion of field-deployed sUAS RSE methods of comparison because of the limited performance and calibration difficulty experienced by small thermal infrared sensors [35,36]. Full-scale incorporation of sUAS thermal-based RSE remains elusive to most practitioners due to the fact that transmissivity and atmospheric radiance vary dramatically throughout the day [37] and that most sUAS data collection missions can last several hours.
A more direct comparison of RSE between EOS and sUAS approaches, due in part to the prevalence of low-cost commercial multispectral sensors for sUAS, would provide practitioners a means by which to evaluate coarse EOS RSE estimates to finer resolution RSE estimates from local collections. For cross-platform RSE comparison, without inclusion of thermal measurements, it is necessary then to focus on methods that use photosynthetic productivity as a proxy for ET using available spectral response in the 400-900 nm domain ( Figure 1). This transfer function approach to RSE is increasingly driven by studies of robust empirical relation between image-derived estimates of vegetation greenness (i.e., NDVI) and ET [38], wherein cellular light energy in red wavelengths is absorbed by chlorophylls and near-infrared light energy is reflected by plant lignin and cellulose. It should be noted here that robust comparisons between EOS data products requires harmonization techniques [39] not employed here. For the purposes of this study, we also explore the utility of OpenET NDVI-ET and SSEBop estimates as described in more detail below.

Research Questions
Given the importance of ET in water management generally, and the potential uncertainty of EOS ET estimation for large seemingly homogeneous cropping areas with temporal variance, we explored several research questions to help parameterize the spatiotemporal uncertainty of EOS ET estimation using sUAS RSE. Answers to these research questions can provide practitioners valuable qualifiers of uncertainty as EOS ET estimates are now more routinely used in determining irrigation volumes and scheduling. Through cross-platform comparisons, we asked: 1. Do EOS and sUAS ET estimates fundamentally differ over same period of time? In other words, over the course of the growing season, does either the mean or variance of sUAS and EOS ET estimates differ? While we expect that ET estimates to similarly track the growing season (e.g., peak ET demand in mid-summer), it remains unknown if the variance in ET estimates also track either the growing season or each other, sUAS compared to EOS. By using the fixed temporal domain of the growing season and coincident measurements, we can compare upscaled sUAS ET to EOS ET estimates to determine if temporal variance is uniformly distributed or varies across time. 2. Do EOS and sUAS ET estimates fundamentally differ over the same study domain and same period of time? While we expect that the mean ET should be comparable from pixel to pixel, it remains unknown if the variance will differ in either within an EOS pixel or across the pixel domain. By using the fixed spatial domain of the vineyard block and by comparing upscaled sUAS ET to EOS ET estimates within and across pixels, we can retain the spatial variance of sUAS to determine if sUAS spatial variance is greater than EOS pixel to pixel variance. 3. Do "mixed pixels" inherent in EOS data obscure important signals in ET estimation? Over the course of a growing season, biomass and leaf area index can change, altering reflectance and ultimately energy balance models. We evaluated the canopy fraction using high spatial resolution sUAS imagery and vineyard canopy structure to determine if unobscured plant reflectance values were better proxies for ET estimation.

Materials and Methods
We evaluated cross-platform data products and methods to estimate ET in a high-yield commercial winegrape vineyard.
We used small Unmanned Aerial Systems (sUAS) to capture high-spatiotemporal, multispectral imagery to compare to well-established EOS Landsat 8 (https://landsat.gsfc.nasa.gov/landsat-8/) and Sentinel 2 (https://sentinel.esa.int/web/sentinel/ missions/sentinel-2) systems. We timed our sUAS field campaigns to coincide with EOS overpasses. We used EOS imagery to calculate ET using both a reflectance based NDVI approach and the emerging EEFlux approach. For EOS ET estimates, Sentinel 2 Bottom of Atmosphere reflectance imagery was used for the OpenET NDVI model and Landsat 8 was used for EEFlux. sUAS imagery was calibrated, stitched, and georeferenced to produce reflectance and structure-from-motion (SfM) surface model products. The OpenET NDVI model approach [31] was used to generate high resolution sUAS-based ET estimates from reflectance, and surface models were used to identify canopy and ground pixels to mask soil in subsequent analysis. As single date LiDAR collection was used to validate SfM canopy models. Each platform's raster-based ET estimates were resampled to have consistent spatial origin, extent, resolution and index to track coincident pixel value and change throughout the season and across data source. We used a combination of Ardupilot Mission Planner (v.

Study Site
The study site was a commercial vineyard located in Ripperdan, Madera County, California, USA (

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 (ET ref ) is typically based on the ratio of irrigated alfalfa, the reference crop, to ET o for grass. ET ref was collected from Spatial CIMIS [40], to obtain daily ET o for the study vineyard throughout the growing season. Spatial CIMIS renders ET o 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 2 from CIMIS stations. ET o is calculated using a modified ASCE-Penman Montieth equation for clipped grass reference from [42] (Equation (1)): where ET o is reference crop evapotranspiration (mm day −1 , R n is net radiation at the crop surface (MJ m −2 day −1 , G is soil heat flux density at the soil surface (MJ m −2 day −1 , T is the daily mean air temperature at 2 m (ms −1 , u 2 is wind speed at 2 m, the difference e s -e a is the saturation vapor pressure deficit (kPa), γ is the psychrometric constant (∼66 Pa K −1 , and ∆ is the slope of the saturation vapor pressure-temperature curve (kPa • C −1 ).

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 ET ref 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.

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 ET o . Daily OpenET aNDV I rasters were used to produce a dimensionless ET fraction (ET f ) raster by using a linear model adjusted for ET o . 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 OpenET aNDV I is as follows: where σ NIR is reflectance in the near infrared band and σ RED is reflectance in the red band (see Figure 1). ET o was derived from from Spatial CIMIS, such that ET a is actual ET (mm day −1 ) and ET o is daily reference ET from Spatial CIMIS (mm day −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.

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 −1 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.

Irrigation Delivery and Canopy Growth
Observed irrigation practices for the growing season resulted in a average of 5.26 mm day −1 across the growing season and 9.96 mm day −1 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 3 of water was applied to the 31.23 ha vineyard block over the course of the growing season April-September.

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 −1 on 16 April to approximately 8 mm day −1 06 June. EOS observations diverge in July-September, where EEFlux estimates peak 06 August at 8.26 mm day −1 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 −1 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 −1 . The greatest deviation occurred on 16 July with a mean difference of 1.75 mm day −1 . 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 −1 on 16 April to the maximum 9.74 mm day −1 on 21 June. The minimum observed EOS mean daily ET of 1.45 mm day −1 occurred on 16 April using L8 EEFlux. The maximum observed EOS mean daily ET of 8.14 mm day −1 occurred on 06 August (Table 2). sUAS NDVI-ET ranged from 0.49 mm day −1 on 07 May to 8.79 mm day −1 on 16 July. Mean sUAS NDVI-ET was lowest on 19 September at 3.58 mm day −1 and greatest on 16 July at 8.04 mm day −1 .   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 a from Landsat 8 EEFlux peaks on 25 August at 7.41 mm day −1 (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 −1 ; sUAS Masked NDVI-ET peaks on 7 July at 7.71 mm day −1 ; and sUAS Unmasked NDVI-ET peaks lower than the other methods on 3 June at 6.32 mm day −1 . 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 −1 and remain there for the remainder of the season. The cumulative ET a over the growing season is shown in Figure 6. Cumulative ET a 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 a of 313,634 m 3 ( Table 3). The unmasked sUAS ET a was within 5000 m 3 of the volume of applied irrigation and was 95 percent confident within a range of 50,000 m 3 . In contrast both EOS methods had about double the range of confidence and the masked sUAS ET a was between with a range of about 80,000 m 3 . Figure 5. Seasonal Daily ET a 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 a in m 3 (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.

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.

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 ET a , 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 a 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 −1 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).

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.

Study Limitations
A limitation of the NDVI-ET approach used here, aside from its comparative simplicity, is that it relies on external sources for ET o . 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 o product since the nearest CIMIS stations were over 30 km away and were potentially unrepresentative of the ET o 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 a 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].

Conclusions
In this study, we evaluated whether EOS and sUAS ET estimates fundamentally differ over same period of time. We showed that more traditional heat balance approaches to EOS RSE can be compromised by field conditions, such as soil saturation and excessive air temperature, and lead to overestimation of ET compared to other methods. Further, we showed that over the course of the growing season, sUAS values largely follow seasonal patterns and track EOS RSE values, but that ET variance is reduced by isolating canopy pixels. While EOS and sUAS ET estimates do not fundamentally differ over the same study domain and same period of time, comparing upscaled sUAS ET to EOS ET estimates within and across pixels, we showed that unmasked high resolution sUAS ET estimates have the highest variance but that EOS EEFlux ET values have high variance too. These "mixed pixels" inherent in EOS data obscure important signals in ET estimation through inclusion of soil fraction, and given that vineyards change over the course of a growing season in terms of biomass and leaf area index can change, caution should be used when evaluating ET whether using traditional energy balance models or RSE proxies. These findings can provide growers with guidance on their own use of new RSE products and sUAS platforms for calibration and verification.