Influence of Model Grid Size on the Estimation of Surface Fluxes Using the Two Source Energy Balance Model and sUAS Imagery in Vineyards

Building a UAV Based System to Acquire High Spatial Resolution Thermal Imagery for Energy Balance Modelling

High spatial resolution and geolocation accuracy canopy evapotranspiration (ET) maps are well suited tools for evaluation of small plot field trials. While creating such a map by use of an energy balance model is routinely performed, the acquisition of the necessary imagery at a suitable quality is still challenging. An UAV based thermal/RGB integrated imaging system was built using the RaspberryPi (RPi) microcomputer as a central unit. The imagery served as input to the two-source energy balance model pyTSEB to derive the ET map. The setup's flexibility and modularity are based on the multiple interfaces provided by the RPi and the software development kit (SDK) provided for the thermal camera. The SDK was installed on the RPi and used to trigger cameras, retrieve and store images and geolocation information from an onboard GNSS rover for PPK processing. The system allows acquisition of 8 cm spatial resolution thermal imagery from a 60 m height of flight and less than 7 cm geolocation accuracy of the mosaicked RGB imagery. Modelled latent heat flux data have been validated against latent heat fluxes measured by eddy covariance stations at two locations with RMSE of 75 W/m² over a two-year study period.

How High to Fly? Mapping Evapotranspiration from Remotely Piloted Aircrafts at Different Elevations

Recent advancements in remotely piloted aircrafts (RPAs) have made frequent, low-flying imagery collection more economical and feasible than ever before. The goal of this work was to create, compare, and quantify uncertainty associated with evapotranspiration (ET) maps generated from different conditions and image capture elevations. We collected optical and thermal data from a commercially irrigated potato (Solanum tuberosum) field in the Wisconsin Central Sands using a quadcopter RPA system and combined multispectral/thermal camera. We conducted eight mission sets (24 total missions) during the 2019 growing season. Each mission set included flights at 90, 60, and 30 m above ground level. Ground reference measurements of surface temperature and soil moisture were collected throughout the domain within 15 min of each RPA mission set. Evapotranspiration values were modeled from the flight data using the High-Resolution Mapping of Evapotranspiration (HRMET) model. We compared HRMET-derived ET estimates to an Eddy Covariance system within the flight domain. Additionally, we assessed uncertainty for each flight using a Monte Carlo approach. Results indicate that the primary source of uncertainty in ET estimates was the optical and thermal data. Despite some additional detectable features at low elevation, we conclude that the tradeoff in resources and computation does not currently justify low elevation flights for annual vegetable crop management in the Midwest USA.

Using Remote Sensing to Estimate Scales of Spatial Heterogeneity to Analyze Evapotranspiration Modeling in a Natural Ecosystem

Understanding the spatial variability in highly heterogeneous natural environments such as savannas and river corridors is an important issue in characterizing and modeling energy fluxes, particularly for evapotranspiration (ET) estimates. Currently, remote-sensing-based surface energy balance (SEB) models are applied widely and routinely in agricultural settings to obtain ET information on an operational basis for use in water resources management. However, the application of these models in natural environments is challenging due to spatial heterogeneity in vegetation cover and complexity in the number of vegetation species existing within a biome. In this research effort, small unmanned aerial systems (sUAS) data were used to study the influence of land surface spatial heterogeneity on the modeling of ET using the Two-Source Energy Balance (TSEB) model. The study area is the San Rafael River corridor in Utah, which is a part of the Upper Colorado River Basin that is characterized by arid conditions and variations in soil moisture status and the type and height of vegetation. First, a spatial variability analysis was performed using a discrete wavelet transform (DWT) to identify a representative spatial resolution/model grid size for adequately solving energy balance components to derive ET. The results indicated a maximum wavelet energy between 6.4 m and 12.8 m for the river corridor area, while the non-river corridor area, which is characterized by different surface types and random vegetation, does not show a peak value. Next, to evaluate the effect of spatial resolution on latent heat flux (LE) estimation using the TSEB model, spatial scales of 6 m and 15 m instead of 6.4 m and 12.8 m, respectively, were used to simplify the derivation of model inputs. The results indicated small differences in the LE values between 6 m and 15 m resolutions, with a slight decrease in detail at 15 m due to losses in spatial variability. Lastly, the instantaneous (hourly) LE was extrapolated/upscaled to daily ET values using the incoming solar radiation (Rs) method. The results indicated that willow and cottonwood have the highest ET rates, followed by grass/shrubs and treated tamarisk. Although most of the treated tamarisk vegetation is in dead/dry condition, the green vegetation growing underneath resulted in a magnitude value of ET.

Evapotranspiration partitioning assessment using a machine-learning-based leaf area index and the two-source energy balance model with sUAV information

Accurate quantification of the partitioning of evapotranspiration (ET) into transpiration and evaporation fluxes is necessary to understanding ecosystem interactions among carbon, water, and energy flux components. ET partitioning can also support the description of atmosphere and land interactions and provide unique insights into vegetation water status. Previous studies have identified leaf area index (LAI) estimation as a key descriptor of biomass conditions needed for the estimation of transpiration and evaporation. LAI estimation in clumped vegetation systems, such as vineyards and orchards, has proven challenging and is strongly related to crop phenological status and canopy management. In this study, a feature extraction model based on previous research was built to generate a total of 202 preliminary variables at a 3.6-by-3.6-meter-grid scale based on submeter-resolution information from a small Unmanned Aerial Vehicle (sUAV) in four commercial vineyards across California. Using these variables, a machine learning model called eXtreme Gradient Boosting (XGBoost) was successfully built for LAI estimation. The XGBoost built-in function requires only six variables relating to vegetation indices and temperature to produce high-accuracy LAI estimation for the vineyard. Using the six-variable XGBoost-based LAI map, two versions of the Two-Source Energy Balance (TSEB) model, TSEB-PT and TSEB-2T were used for energy balance and ET partitioning. Comparing these results with the Eddy-Covariance (EC) tower data, showed that TSEB-PT outperforms TSEB-2T on the estimation of sensible heat flux (within 13% relative error) and surface heat flux (within 34% relative error), while TSEB-2T outperforms TSEB-PT on the estimation of net radiation (within 14% relative error) and latent heat flux (within 2% relative error). For the mature vineyard (north block), TSEB-2T performs better than TSEB-PT in partitioning the canopy latent heat flux with 6.8% relative error and soil latent heat flux with 21.7% relative error; however, for the younger vineyard (south block), TSEB-PT performs better than TSEB-2T in partitioning the canopy latent heat flux with 11.7% relative error and soil latent heat flux with 39.3% relative error.

Assessing Daily Evapotranspiration Methodologies from One-Time-of-Day sUAS and EC Information in the GRAPEX Project

Daily evapotranspiration (ET_d) plays a key role in irrigation water management and is particularly important in drought-stricken areas, such as California and high-value crops. Remote sensing allows for the cost-effective estimation of spatial evapotranspiration (ET), and the advent of small unmanned aerial systems (sUAS) technology has made it possible to estimate instantaneous high-resolution ET at the plant, row, and subfield scales. sUAS estimates ET using “instantaneous” remote sensing measurements with half-hourly/hourly forcing micrometeorological data, yielding hourly fluxes in W/m² that are then translated to a daily scale (mm/day) under two assumptions: (a) relative rates, such as the ratios of ET-to-net radiation (R_n) or ET-to-solar radiation (R_s), are assumed to be constant rather than absolute, and (b) nighttime evaporation (E) and transpiration (T) contributions are negligible. While assumption (a) may be reasonable for unstressed, full cover crops (no exposed soil), the E and T rates may significantly vary over the course of the day for partially vegetated cover conditions due to diurnal variations of soil and crop temperatures and interactions between soil and vegetation elements in agricultural environments, such as vineyards and orchards. In this study, five existing extrapolation approaches that compute the daily ET from the “instantaneous” remotely sensed sUAS ET estimates and the eddy covariance (EC) flux tower measurements were evaluated under different weather, grapevine variety, and trellis designs. Per assumption (b), the nighttime ET contribution was ignored. Each extrapolation technique (evaporative fraction (EF), solar radiation (R_s), net radiation-to-solar radiation (R_n/R_s) ratio, Gaussian (GA), and Sine) makes use of clear skies and quasi-sinusoidal diurnal variations of hourly ET and other meteorological parameters. The sUAS ET estimates and EC ET measurements were collected over multiple years and times from different vineyard sites in California as part of the USDA Agricultural Research Service Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX). Optical and thermal sUAS imagery data at 10 cm and 60 cm, respectively, were collected by the Utah State University AggieAir sUAS Program and used in the Two-Source Energy Balance (TSEB) model to estimate the instantaneous or hourly sUAS ET at overpass time. The hourly ET from the EC measurements was also used to validate the extrapolation techniques. Overall, the analysis using EC measurements indicates that the R_s, EF, and GA approaches presented the best goodness-of-fit statistics for a window of time between 1030 and 1330 PST (Pacific Standard Time), with the R_s approach yielding better agreement with the EC measurements. Similar results were found using TSEB and sUAS data. The 1030–1330 time window also provided the greatest agreement between the actual daily EC ET and the extrapolated TSEB daily ET, with the R_s approach again yielding better agreement with the ground measurements. The expected accuracy of the upscaled TSEB daily ET estimates across all vineyard sites in California is below 0.5 mm/day, (EC extrapolation accuracy was found to be 0.34 mm/day), making the daily scale results from TSEB reliable and suitable for day-to-day water management applications.

Actual Evapotranspiration from UAV Images: A Multi-Sensor Data Fusion Approach

Multispectral imaging using Unmanned Aerial Vehicles (UAVs) has changed the pace of precision agriculture. Actual evapotranspiration (ETa) from the very high spatial resolution of UAV images over agricultural fields can help farmers increase their production at the lowest possible cost. ETa estimation using UAVs requires a full package of sensors capturing the visible/infrared and thermal portions of the spectrum. Therefore, this study focused on a multi-sensor data fusion approach for ETa estimation (MSDF-ET) independent of thermal sensors. The method was based on sharpening the Landsat 8 pixels to UAV spatial resolution by considering the relationship between reference ETa fraction (ETrf) and a Vegetation Index (VI). Four Landsat 8 images were processed to calculate ETa of three UAV images over three almond fields. Two flights coincided with the overpasses and one was in between two consecutive Landsat 8 images. ETrf was chosen instead of ETa to interpolate the Landsat 8-derived ETrf images to obtain an ETrf image on the UAV flight. ETrf was defined as the ratio of ETa to grass reference evapotranspiration (ETr), and the VIs tested in this study included the Normalized Difference Vegetation Index (NDVI), Soil Adjusted Vegetation Index (SAVI), Enhanced Vegetation Index (EVI), Normalized Difference Water Index (NDWI), and Land Surface Water Index (LSWI). NDVI performed better under the study conditions. The MSDF-ET-derived ETa showed strong correlations against measured ETa, UAV- and Landsat 8-based METRIC ETa. Also, visual comparison of the MSDF-ET ETa maps was indicative of a promising performance of the method. In sum, the resulting ETa had a higher spatial resolution compared with thermal-based ETa without the need for the Albedo and hot/cold pixels selection procedure. However, wet soils were poorly detected, and in cases of continuous cloudy Landsat pixels the long interval between the images may cause biases in ETa estimation from the MSDF-ET method. Generally, the MSDF-ET method reduces the need for very high resolution thermal information from the ground, and the calculations can be conducted on a moderate-performance computer system because the main image processing is applied on Landsat images with coarser spatial resolutions.

UAVs for Vegetation Monitoring: Overview and Recent Scientific Contributions

This paper reviewed a set of twenty-one original and innovative papers included in a special issue on UAVs for vegetation monitoring, which proposed new methods and techniques applied to diverse agricultural and forestry scenarios. Three general categories were considered: (1) sensors and vegetation indices used, (2) technological goals pursued, and (3) agroforestry applications. Some investigations focused on issues related to UAV flight operations, spatial resolution requirements, and computation and data analytics, while others studied the ability of UAVs for characterizing relevant vegetation features (mainly canopy cover and crop height) or for detecting different plant/crop stressors, such as nutrient content/deficiencies, water needs, weeds, and diseases. The general goal was proposing UAV-based technological solutions for a better use of agricultural and forestry resources and more efficient production with relevant economic and environmental benefits.

Development of High Performance Computing Tools for Estimation of High-Resolution Surface Energy Balance Products Using sUAS Information

sUAS (small-Unmanned Aircraft System) and advanced surface energy balance models allow detailed assessment and monitoring (at plant scale) of different (agricultural, urban, and natural) environments. Significant progress has been made in the understanding and modeling of atmosphere-plant-soil interactions and numerical quantification of the internal processes at plant scale. Similarly, progress has been made in ground truth information comparison and validation models. An example of this progress is the application of sUAS information using the Two-Source Surface Energy Balance (TSEB) model in commercial vineyards by the Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment - GRAPEX Project in California. With advances in frequent sUAS data collection for larger areas, sUAS information processing becomes computationally expensive on local computers. Additionally, fragmentation of different models and tools necessary to process the data and validate the results is a limiting factor. For example, in the referred GRAPEX project, commercial software (ArcGIS and MS Excel) and Python and Matlab code are needed to complete the analysis. There is a need to assess and integrate research conducted with sUAS and surface energy balance models in a sharing platform to be easily migrated to high performance computing (HPC) resources. This research, sponsored by the National Science Foundation FAIR Cyber Training Fellowships, is integrating disparate software and code under a unified language (Python). The Python code for estimating the surface energy fluxes using TSEB2T model as well as the EC footprint analysis code for ground truth information comparison were hosted in myGeoHub site https://mygeohub.org/ to be reproducible and replicable.

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.

Implications of Soil and Canopy Temperature Uncertainty in the Estimation of Surface Energy Fluxes Using TSEB2T and High-resolution Imagery in Commercial Vineyards

Estimation of surface energy fluxes using thermal remote sensing-based energy balance models (e.g., TSEB2T) involves the use of local micrometeorological input data of air temperature, wind speed, and incoming solar radiation, as well as vegetation cover and accurate land surface temperature (LST). The physically based Two-source Energy Balance with a Dual Temperature (TSEB2T) model separates soil and canopy temperature (Ts and Tc) to estimate surface energy fluxes including Rn, H, LE, and G. The estimation of Ts and Tc components for the TSEB2T model relies on the linear relationship between the composite land surface temperature and a vegetation index, namely NDVI. While canopy and soil temperatures are controlling variables in the TSEB2T model, they are influenced by the NDVI threshold values, where the uncertainties in their estimation can degrade the accuracy of surface energy flux estimation. Therefore, in this research effort, the effect of uncertainty in Ts and Tc estimation on surface energy fluxes will be examined by applying a Monte Carlo simulation on NDVI thresholds used to define canopy and soil temperatures. The spatial information used is available from multispectral imagery acquired by the AggieAir sUAS Program at Utah State University over vineyards near Lodi, California as part of the ARS-USDA Agricultural Research Service's Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX) project. The results indicate that LE is slightly sensitive to the uncertainty of NDVIs and NDVIc. The observed relative error of LE corresponding to NDVIs uncertainty was between -1% and 2%, while for NDVIc uncertainty, the relative error was between -2.2% and 1.2%. However, when the combined NDVIs and NDVIc uncertainties were used simultaneously, the domain of the observed relative error corresponding to the absolute values of |ΔLE| was between 0% and 4%.

To What Extend Does the Eddy Covariance Footprint Cutoff Influence the Estimation of Surface Energy Fluxes Using Two Source Energy Balance Model and High-Resolution Imagery in Commercial Vineyards?

Validation of surface energy fluxes from remote sensing sources is performed using instantaneous field measurements obtained from eddy covariance (EC) instrumentation. An eddy covariance measurement is characterized by a footprint function / weighted area function that describes the mathematical relationship between the spatial distribution of surface flux sources and their corresponding magnitude. The orientation and size of each flux footprint / source area depends on the micro-meteorological conditions at the site as measured by the EC towers, including turbulence fluxes, friction velocity (ustar), and wind speed, all of which influence the dimensions and orientation of the footprint. The total statistical weight of the footprint is equal to unity. However, due to the large size of the source area / footprint, a statistical weight cutoff of less than one is considered, ranging between 0.85 and 0.95, to ensure that the footprint model is located inside the study area. This results in a degree of uncertainty when comparing the modeled fluxes from remote sensing energy models (i.e., TSEB2T) against the EC field measurements. In this research effort, the sensitivity of instantaneous and daily surface energy flux estimates to footprint weight cutoffs are evaluated using energy balance fluxes estimated with multispectral imagery acquired by AggieAir sUAS (small Unmanned Aerial Vehicle) over commercial vineyards near Lodi, California, as part of the ARS-USDA Agricultural Research Service's Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX) project. The instantaneous fluxes from the eddy covariance tower will be compared against instantaneous fluxes obtained from different TSEB2T aggregated footprint weights (cutoffs). The results indicate that the size, shape, and weight of pixels inside the footprint source area are strongly influenced by the cutoff values. Small cutoff values, such as 0.3 and 0.35, yielded high weights for pixels located within the footprint domain, while large cutoffs, such as 0.9 and 0.95, result in low weights. The results also indicate that the distribution of modelled LE values within the footprint source area are influenced by the cutoff values. A wide variation in LE was observed at high cutoffs, such as 0.90 and 0.95, while a low variation was observed at small cutoff values, such as 0.3. This happens due to the large number of pixel units involved inside the footprint domain when using high cutoff values, whereas a limited number of pixels are obtained at lower cutoff values.