Evaluation of TSEB turbulent fluxes using different methods for the retrieval of soil and canopy component temperatures from UAV thermal and multispectral imagery.

The thermal-based Two-Source Energy Balance (TSEB) model partitions the evapotranspiration (ET) and energy fluxes from vegetation and soil components providing the capability for estimating soil evaporation (E) and canopy transpiration (T). However, it is crucial for ET partitioning to retrieve reliable estimates of canopy and soil temperatures and net radiation, as the latter determines the available energy for water and heat exchange from soil and canopy sources. These two factors become especially relevant in row crops with wide spacing and strongly clumped vegetation such as vineyards and orchards. To better understand these effects, very high spatial resolution remote-sensing data from an unmanned aerial vehicle were collected over vineyards in California, as part of the Grape Remote sensing and Atmospheric Profile and Evapotranspiration eXperiment and used in four different TSEB approaches to estimate the component soil and canopy temperatures, and ET partitioning between soil and canopy. Two approaches rely on the use of composite T rad, and assume initially that the canopy transpires at the Priestley-Taylor potential rate. The other two algorithms are based on the contextual relationship between optical and thermal imagery partition T rad into soil and canopy component temperatures, which are then used to drive the TSEB without requiring a priori assumptions regarding initial canopy transpiration rate. The results showed that a simple contextual algorithm based on the inverse relationship of a vegetation index and T rad to derive soil and canopy temperatures yielded the closest agreement with flux tower measurements. The utility in very high-resolution remote-sensing data for estimating ET and E and T partitioning at the canopy level is also discussed.

Assessment of different methods for shadow detection in high-resolution optical imagery and evaluation of shadow impact on calculation of NDVI, and evapotranspiration.

Significant efforts have been made recently in the application of high-resolution remote sensing imagery (i.e., sub-meter) captured by unmanned aerial vehicles (UAVs) for precision agricultural applications for high-value crops such as wine grapes. However, at such high resolution, shadows will appear in the optical imagery effectively reducing the reflectance and emission signal received by imaging sensors. To date, research that evaluates procedures to identify the occurrence of shadows in imagery produced by UAVs is limited. In this study, the performance of four different shadow detection methods used in satellite imagery was evaluated for high-resolution UAV imagery collected over a California vineyard during the Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX) field campaigns. The performance of the shadow detection methods was compared and impacts of shadowed areas on the normalized difference vegetation index (NDVI) and estimated evapotranspiration (ET) using the Two-Source Energy Balance (TSEB) model are presented. The results indicated that two of the shadow detection methods, the supervised classification and index-based methods, had better performance than two other methods. Furthermore, assessment of shadowed pixels in the vine canopy led to significant differences in the calculated NDVI and ET in areas affected by shadows in the high-resolution imagery. Shadows are shown to have the greatest impact on modeled soil heat flux, while net radiation and sensible heat flux are less affected. Shadows also have an impact on the modeled Bowen ratio (ratio of sensible to latent heat) which can be used as an indicator of vine stress level.

Spatial Scale Gap Filling Using an Unmanned Aerial System: A Statistical Downscaling Method for Applications in Precision Agriculture.

Applications of satellite-borne observations in precision agriculture (PA) are often limited due to the coarse spatial resolution of satellite imagery. This paper uses high-resolution airborne observations to increase the spatial resolution of satellite data for related applications in PA. A new variational downscaling scheme is presented that uses coincident aerial imagery products from "AggieAir", an unmanned aerial system, to increase the spatial resolution of Landsat satellite data. This approach is primarily tested for downscaling individual band Landsat images that can be used to derive normalized difference vegetation index (NDVI) and surface soil moisture (SSM). Quantitative and qualitative results demonstrate promising capabilities of the downscaling approach enabling effective increase of the spatial resolution of Landsat imageries by orders of 2 to 4. Specifically, the downscaling scheme retrieved the missing high-resolution feature of the imageries and reduced the root mean squared error by 15, 11, and 10 percent in visual, near infrared, and thermal infrared bands, respectively. This metric is reduced by 9% in the derived NDVI and remains negligibly for the soil moisture products.

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.

Incorporation of Unmanned Aerial Vehicle (UAV) Point Cloud Products into Remote Sensing Evapotranspiration Models.

In recent years, the deployment of satellites and unmanned aerial vehicles (UAVs) has led to production of enormous amounts of data and to novel data processing and analysis techniques for monitoring crop conditions. One overlooked data source amid these efforts, however, is incorporation of 3D information derived from multi-spectral imagery and photogrammetry algorithms into crop monitoring algorithms. Few studies and algorithms have taken advantage of 3D UAV information in monitoring and assessment of plant conditions. In this study, different aspects of UAV point cloud information for enhancing remote sensing evapotranspiration (ET) models, particularly the Two-Source Energy Balance Model (TSEB), over a commercial vineyard located in California are presented. Toward this end, an innovative algorithm called Vegetation Structural-Spectral Information eXtraction Algorithm (VSSIXA) has been developed. This algorithm is able to accurately estimate height, volume, surface area, and projected surface area of the plant canopy solely based on point cloud information. In addition to biomass information, it can add multi-spectral UAV information to point clouds and provide spectral-structural canopy properties. The biomass information is used to assess its relationship with in situ Leaf Area Index (LAI), which is a crucial input for ET models. In addition, instead of using nominal field values of plant parameters, spatial information of fractional cover, canopy height, and canopy width are input to the TSEB model. Therefore, the two main objectives for incorporating point cloud information into remote sensing ET models for this study are to (1) evaluate the possible improvement in the estimation of LAI and biomass parameters from point cloud information in order to create robust LAI maps at the model resolution and (2) assess the sensitivity of the TSEB model to using average/nominal values versus spatially-distributed canopy fractional cover, height, and width information derived from point cloud data. The proposed algorithm is tested on imagery from the Utah State University AggieAir sUAS Program as part of the ARS-USDA GRAPEX Project (Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment) collected since 2014 over multiple vineyards located in California. The results indicate a robust relationship between in situ LAI measurements and estimated biomass parameters from the point cloud data, and improvement in the agreement between TSEB model output of ET with tower measurements when employing LAI and spatially-distributed canopy structure parameters derived from the point cloud data.

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

Evapotranspiration (ET) is a key variable for hydrology and irrigation water management, with significant importance in drought-stricken regions of the western US. This is particularly true for California, which grows much of the high-value perennial crops in the US. The advent of small Unmanned Aerial System (sUAS) with sensor technology similar to satellite platforms allows for the estimation of high-resolution ET at plant spacing scale for individual fields. However, while multiple efforts have been made to estimate ET from sUAS products, the sensitivity of ET models to different model grid size/resolution in complex canopies, such as vineyards, is still unknown. The variability of row spacing, canopy structure, and distance between fields makes this information necessary because additional complexity processing individual fields. Therefore, processing the entire image at a fixed resolution that is potentially larger than the plant-row separation is more efficient. From a computational perspective, there would be an advantage to running models at much coarser resolutions than the very fine native pixel size from sUAS imagery for operational applications. In this study, the Two-Source Energy Balance with a dual temperature (TSEB2T) model, which uses remotely sensed soil/substrate and canopy temperature from sUAS imagery, was used to estimate ET and identify the impact of spatial domain scale under different vine phenological conditions. The analysis relies upon high-resolution imagery collected during multiple years and times by the Utah State University AggieAir™ sUAS program over a commercial vineyard located near Lodi, California. This project is part of the USDA-Agricultural Research Service Grape Remote Sensing Atmospheric Profile and Evapotranspiration eXperiment (GRAPEX). Original spectral and thermal imagery data from sUAS were at 10 cm and 60 cm per pixel, respectively, and multiple spatial domain scales (3.6, 7.2, 14.4, and 30 m) were evaluated and compared against eddy covariance (EC) measurements. Results indicated that the TSEB2T model is only slightly affected in the estimation of the net radiation (R n ) and the soil heat flux (G) at different spatial resolutions, while the sensible and latent heat fluxes (H and LE, respectively) are significantly affected by coarse grid sizes. The results indicated overestimation of H and underestimation of LE values, particularly at Landsat scale (30 m). This refers to the non-linear relationship between the land surface temperature (LST) and the normalized difference vegetation index (NDVI) at coarse model resolution. Another predominant reason for LE reduction in TSEB2T was the decrease in the aerodynamic resistance (R a ), which is a function of the friction velocity F*) that varies with mean canopy height and roughness length. While a small increase in grid size can be implemented, this increase should be limited to less than twice the smallest row spacing present in the sUAS imagery. The results also indicated that the mean LE at field scale is reduced by 10% to 20% at coarser resolutions, while the with-in field variability in LE values decreased significantly at the larger grid sizes and ranged between approximately 15% and 45%. This implies that, while the field-scale values of LE are fairly reliable at larger grid sizes, the with-in field variability limits its use for precision agriculture applications.

The impact of shadows on partitioning of radiometric temperature to canopy and soil temperature based on the contextual two-source energy balance model (TSEB-2T).

Tests of the most recent version of the two-source energy balance model have demonstrated that canopy and soil temperatures can be retrieved from high-resolution thermal imagery captured by an unmanned aerial vehicle (UAV). This work has assumed a linear relationship between vegetation indices (VIs) and radiometric temperature in a square grid (i.e., 3.6 m × 3.6 m) that is coarser than the resolution of the imagery acquired by the UAV. In this method, with visible, near infrared (VNIR), and thermal bands available at the same high-resolution, a linear fit can be obtained over the pixels located in a grid, where the x-axis is a vegetation index (VI) and the y-axis is radiometric temperature. Next, with an accurate VI threshold that separates soil and vegetation pixels from one another, the corresponding soil and vegetation temperatures can be extracted from the linear equation. Although this method is simpler than other approaches, such as TSEB with Priestly-Taylor (TSEB-PT), it could be sensitive to VIs and the parameters that affect VIs, such as shadows. Recent studies have revealed that, on average, the values of VIs, such as normalized difference vegetation index (NDVI) and leaf area index (LAI), that are located in sunlit areas are greater than those in shaded areas. This means that involving or compensating for shadows will affect the linear relationship parameters (slope and bias) between radiometric temperature and VI, as well as thresholds that separate soil and vegetation pixels. This study evaluates the impact of shadows on the retrieval of canopy and soil temperature data from four UAV images before and after applying shadow compensation techniques. The retrieved temperatures, using the TSEB-2T approach, both before and after shadow correction, are compared to the average temperature values for both soil and canopy in each grid. The imagery was acquired by the Utah State University AggieAir UAV system over a commercial vineyard located in California as part of the USDA Agricultural Research Service Grape Remote sensing Atmospheric Profile and Evapotranspiration Experiment (GRAPEX) Program during 2014 to 2016. The results of this study show when it is necessary to employ shadow compensation methods to retrieve vegetation and soil temperature directly.

Validation of digital surface models (DSMs) retrieved from unmanned aerial vehicle (UAV) point clouds using geometrical information from shadows.

Theoretically, the appearance of shadows in aerial imagery is not desirable for researchers because it leads to errors in object classification and bias in the calculation of indices. In contrast, shadows contain useful geometrical information about the objects blocking the light. Several studies have focused on estimation of building heights in urban areas using the length of shadows. This type of information can be used to predict the population of a region, water demand, etc., in urban areas. With the emergence of unmanned aerial vehicles (UAVs) and the availability of high- to super-high-resolution imagery, the important questions relating to shadows have received more attention. Three-dimensional imagery generated using UAV-based photogrammetric techniques can be very useful, particularly in agricultural applications such as in the development of an empirical equation between biomass or yield and the geometrical information of canopies or crops. However, evaluating the accuracy of the canopy or crop height requires labor-intensive efforts. In contrast, the geometrical relationship between the length of the shadows and the crop or canopy height can be inversely solved using the shadow length measured. In this study, object heights retrieved from UAV point clouds are validated using the geometrical shadow information retrieved from three sets of high-resolution imagery captured by Utah State University's AggieAir UAV system. These flights were conducted in 2014 and 2015 over a commercial vineyard located in California for the USDA Agricultural Research Service Grape Remote sensing Atmospheric Profile and Evapotranspiration Experiment (GRAPEX) Program. The results showed that, although this approach could be computationally expensive, it is faster than fieldwork and does not require an expensive and accurate instrument such as a real-time kinematic (RTK) GPS.

Estimation of surface thermal emissivity in a vineyard for UAV microbolometer thermal cameras using NASA HyTES hyperspectral thermal, Landsat and AggieAir optical data.

Microbolometer thermal cameras in UAVs and manned aircraft allow for the acquisition of high-resolution temperature data, which, along with optical reflectance, contributes to monitoring and modeling of agricultural and natural environments. Furthermore, these temperature measurements have facilitated the development of advanced models of crop water stress and evapotranspiration in precision agriculture and heat fluxes exchanges in small river streams and corridors. Microbolometer cameras capture thermal information at blackbody or radiometric settings (narrowband emissivity equates to unity). While it is customary that the modeler uses assumed emissivity values (e.g. 0.99-0.96 for agricultural and environmental settings); some applications (e.g. Vegetation Health Index), and complex models such as energy balance-based models (e.g. evapotranspiration) could benefit from spatial estimates of surface emissivity for true or kinetic temperature mapping. In that regard, this work presents an analysis of the spectral characteristics of a microbolometer camera with regard to emissivity, along with a methodology to infer thermal emissivity spatially based on the spectral characteristics of the microbolometer camera. For this work, the MODIS UCBS Emissivity Library, NASA HyTES hyperspectral emissivity, Landsat, and Utah State University AggieAir UAV surface reflectance products are employed. The methodology is applied to a commercial vineyard agricultural setting located in Lodi, California, where HyTES, Landsat, and AggieAir UAV spatial data were collected in the 2014 growing season. Assessment of the microbolometer spectral response with regards to emissivity and emissivity modeling performance for the area of study are presented and discussed.

Implications of sensor inconsistencies and remote sensing error in the use of small unmanned aerial systems for generation of information products for agricultural management.

Small, unmanned aerial systems (sUAS) for remote sensing represent a relatively new and growing technology to support decisions for agricultural operations. The size and power limitations of these systems present challenges for the weight, size, and capability of the sensors that can be carried, as well as the geographical coverage that is possible. These factors, together with a lack of standards for sensor technology, its deployment, and data analysis, lead to uncertainties in data quality that can be difficult to detect or characterize. These, in turn, limit comparability between data from different sources and, more importantly, imply limits on the analyses that can be accomplished with the data that are acquired with sUAS. This paper offers a simple statistical examination of the implications toward information products of an array of sensor data uncertainty issues. The analysis relies upon high-resolution data collected in 2016 over a commercial vineyard, located near Lodi, California, for the USD A Agricultural Research Service Grape Remote sensing Atmospheric Profile and Evapotranspiration experiment (GRAPEX) Program. A Monte Carlo analysis is offered of how uncertainty in sensor spectral response and/or orthorectification accuracy can affect the estimation of information products of potential interest to growers, as illustrated in the form of common vegetation indices.

Behavior of vegetation/soil indices in shaded and sunlit pixels and evaluation of different shadow compensation methods using UAV high-resolution imagery over vineyards.

In high-resolution imagery, shadows may cause problems in object segmentation and recognition due to their low reflectance. For instance, the spectral reflectance of shadows and water are similar, particularly in the visible band. In precision agriculture, the vegetation condition in terms of plant water use, plant water stress, and chlorophyll content can be estimated using vegetation indices. Normalized Difference Vegetation Index (NDVI), Leaf Area Index (LAI) and Enhanced Vegetation Index (EVI) are widely used vegetation indices for characterizing the condition of the vegetation. In addition, many soil indices have been developed for describing soil characteristics, such as Soil-Adjusted Vegetation Index (SAVI). However, shadows can have an influence on the performance of these vegetation and soil indices.Moreover, enhancing spatial resolution heightens the impact of shadows in the imagery. In this study, the behavior of vegetation and soil indices are evaluated using four sets of high-resolution imagery captured by the Utah State University AggieAir unmanned aerial vehicle (UAV) system. These indices were obtained from flights conducted in 2014, 2015, and 2016 over a commercial vineyard located in California for the USDA Agricultural Research Service Grape Remote sensing Atmospheric Profile and Evapotranspiration Experiment (GRAPEX) Program. Different shadow restoration methods are used to alleviate the impact of shadows in information products that might be developed from the high-resolution imagery. The histogram pattern of vegetation and soil indices before and after shadow compensation, are compared using sanalysi of variance (ANOVA). The results of this study indicate how shadows can affect the vegetation/soil indices and whether shadow compensation methods are able to remove the statistical difference between sunlit and shadowed vegetation/soil indices.

Inter-comparison of thermal measurements using ground-based sensors, UAV thermal cameras, and eddy covariance radiometers.

With the increasing availability of thermal proximity sensors, UAV-borne cameras, and eddy covariance radiometers there may be an assumption that information produced by these sensors is interchangeable or compatible. This assumption is often held for estimation of agricultural parameters such as canopy and soil temperature, energy balance components, and evapotranspiration. Nevertheless, environmental conditions, calibration, and ground settings may affect the relationship between measurements from each of these thermal sensors. This work presents a comparison between proximity infrared radiometer (IRT) sensors, microbolometer thermal cameras used in UAVs, and thermal radiometers used in eddy covariance towers in an agricultural setting. The information was collected in the 2015 and 2016 irrigation seasons at a commercial vineyard located in California for the USDA Agricultural Research Service Grape Remote Sensing Atmospheric Profile and Evapotranspiration Experiment (GRAPEX) Program. Information was captured at different times during diurnal cycles, and IRT and radiometer footprint areas were calculated for comparison with UAV thermal raster information. Issues such as sensor accuracy, the location of IRT sensors, diurnal temperature changes, and surface characterizations are presented.

Support vector machines (SVMs) for monitoring network design.

In this paper we present a hydrologic application of a new statistical learning methodology called support vector machines (SVMs). SVMs are based on minimization of a bound on the generalized error (risk) model, rather than just the mean square error over a training set. Due to Mercer's conditions on the kernels, the corresponding optimization problems are convex and hence have no local minima. In this paper, SVMs are illustratively used to reproduce the behavior of Monte Carlo-based flow and transport models that are in turn used in the design of a ground water contamination detection monitoring system. The traditional approach, which is based on solving transient transport equations for each new configuration of a conductivity field, is too time consuming in practical applications. Thus, there is a need to capture the behavior of the transport phenomenon in random media in a relatively simple manner. The objective of the exercise is to maximize the probability of detecting contaminants that exceed some regulatory standard before they reach a compliance boundary, while minimizing cost (i.e., number of monitoring wells). Application of the method at a generic site showed a rather promising performance, which leads us to believe that SVMs could be successfully employed in other areas of hydrology. The SVM was trained using 510 monitoring configuration samples generated from 200 Monte Carlo flow and transport realizations. The best configurations of well networks selected by the SVM were identical with the ones obtained from the physical model, but the reliabilities provided by the respective networks differ slightly.