Effects of meteorological and land surface modeling uncertainty on errors in winegrape ET calculated with SIMS.

Characterization of model errors is important when applying satellite-driven evapotranspiration (ET) models to water resource management problems. This study examines how uncertainty in meteorological forcing data and land surface modeling propagate through to errors in final ET data calculated using the Satellite Irrigation Management Support (SIMS) model, a computationally efficient ET model driven with satellite surface reflectance values. The model is applied to three instrumented winegrape vineyards over the 2017-2020 time period and the spatial and temporal variation in errors are analyzed. We illustrate how meteorological data inputs can introduce biases that vary in space and at seasonal timescales, but that can persist from year to year. We also observe that errors in SIMS estimates of land surface conductance can have a particularly strong dependence on time of year. Overall, meteorological inputs introduced RMSE of 0.33-0.65 mm/day (7-27%) across sites, while SIMS introduced RMSE of 0.55-0.83 mm/day (19-24%). The relative error contribution from meteorological inputs versus SIMS varied across sites; errors from SIMS were larger at one site, errors from meteorological inputs were larger at a second site, and the error contributions were of equal magnitude at the third site. The similar magnitude of error contributions is significant given that many satellite-driven ET models differ in their approaches to estimating land surface conductance, but often rely on similar or identical meteorological forcing data. The finding is particularly notable given that SIMS makes assumptions about the land surface (no soil evaporation or plant water stress) that do not always hold in practice. The results of this study show that improving SIMS by eliminating these assumptions would result in meteorological inputs dominating the error budget of the model on the whole. This finding underscores the need for further work on characterizing spatial uncertainty in the meteorological forcing of ET. Supplementary Information: The online version contains supplementary material available at 10.1007/s00271-022-00808-9.

Improving the spatiotemporal resolution of remotely sensed ET information for water management through Landsat, Sentinel-2, ECOSTRESS and VIIRS data fusion.

Robust information on consumptive water use (evapotranspiration, ET) derived from remote sensing can significantly benefit water decision-making in agriculture, informing irrigation schedules and water management plans over extended regions. To be of optimal utility for operational usage, these remote sensing ET data should be generated at the sub-field spatial resolution and daily-to-weekly timesteps commensurate with the scales of water management activities. However, current methods for field-scale ET retrieval based on thermal infrared (TIR) imaging, a valuable diagnostic of canopy stress and surface moisture status, are limited by the temporal revisit of available medium-resolution (100 m or finer) thermal satellite sensors. This study investigates the efficacy of a data fusion method for combining information from multiple medium-resolution sensors toward generating high spatiotemporal resolution ET products for water management. TIR data from Landsat and ECOSTRESS (both at ~ 100-m native resolution), and VIIRS (375-m native) are sharpened to a common 30-m grid using surface reflectance data from the Harmonized Landsat-Sentinel dataset. Periodic 30-m ET retrievals from these combined thermal data sources are fused with daily retrievals from unsharpened VIIRS to generate daily, 30-m ET image timeseries. The accuracy of this mapping method is tested over several irrigated cropping systems in the Central Valley of California in comparison with flux tower observations, including measurements over irrigated vineyards collected in the GRAPEX campaign. Results demonstrate the operational value added by the augmented TIR sensor suite compared to Landsat alone, in terms of capturing daily ET variability and reduced latency for real-time applications. The method also provides means for incorporating new sources of imaging from future planned thermal missions, further improving our ability to map rapid changes in crop water use at field scales.

Application of a remote-sensing three-source energy balance model to improve evapotranspiration partitioning in vineyards.

Improved accuracy of evapotranspiration (ET) estimation, including its partitioning between transpiration (T) and surface evaporation (E), is key to monitor agricultural water use in vineyards, especially to enhance water use efficiency in semi-arid regions such as California, USA. Remote-sensing methods have shown great utility in retrieving ET from surface energy balance models based on thermal infrared data. Notably, the two-source energy balance (TSEB) has been widely and robustly applied in numerous landscapes, including vineyards. However, vineyards add an additional complexity where the landscape is essentially made up of two distinct zones: the grapevine and the interrow, which is often seasonally covered by an herbaceous cover crop. Therefore, it becomes more complex to disentangle the various contributions of the different vegetation elements to total ET, especially through TSEB, which assumes a single vegetation source over a soil layer. As such, a remote-sensing-based three-source energy balance (3SEB) model, which essentially adds a vegetation source to TSEB, was applied in an experimental vineyard located in California's Central Valley to investigate whether it improves the depiction of the grapevine-interrow system. The model was applied in four different blocks in 2019 and 2020, where each block had an eddy-covariance (EC) tower collecting continuous flux, radiometric, and meteorological measurements. 3SEB's latent and sensible heat flux retrievals were accurate with an overall RMSD ~ 50 W/m2 compared to EC measurements. 3SEB improved upon TSEB simulations, with the largest differences being concentrated in the spring season, when there is greater mixing between grapevine foliage and the cover crop. Additionally, 3SEB's modeled ET partitioning (T/ET) compared well against an EC T/ET retrieval method, being only slightly underestimated. Overall, these promising results indicate 3SEB can be of great utility to vineyard irrigation management, especially to improve T/ET estimations and to quantify the contribution of the cover crop to ET. Improved knowledge of T/ET can enhance grapevine water stress detection to support irrigation and water resource management. Supplementary Information: The online version contains supplementary material available at 10.1007/s00271-022-00787-x.

Intermittency of water vapor fluxes from vineyards during light wind and convective conditions.

Vineyards in many semi-arid regions globally face limited water resources. Monitoring évapotranspiration (ET) of vineyards is critical for water resource management, but remains difficult due to the complex biophysics of the surfaces. Both measurement and modeling approaches for estimating turbulent water vapor transport rely on implicit assumptions that exchanges occur in a reasonably regular fashion over the time scales generally used for averaging. However, heterogeneous vegetation in semi-arid climates, such as many vineyards, presents inherent factors, including canopy row/row space structure and frequent periods of light wind, unstable conditions, that can create episodic transport characteristics. Eddy covariance data were collected above and within the canopy of two vineyards in the Central Valley of California during the Grape Remote sensing Atmospheric Profile & Evapotranspiration experiment (GRAPEX). The goal was to document and quantify the existence of intermittent turbulence transport of water vapor, and associated episodic canopy venting. These effects were found to correlate with periods light winds and highly unstable/convective conditions. Power and cross-spectra for intermittent periods documented enhancement of low-frequency water vapor exchange events compared to more steady periods, and diminished time scale correlation between humidity within the canopy and above the canopy. Analyses show that intermittent cases can necessitate longer flux-averaging periods (up to 2 h) than more steady conditions. Episodic exchange events were isolated and summed to determine their relative contribution to the overall water vapor flux. Since light wind, unstable conditions are relatively common in many arid vineyard regions, these findings have implications for mechanistic ET models that rely on time-averaged vertical gradients, which implies reasonably steady transport.

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.

Particulate emissions calculations from fall tillage operations using point and remote sensors.

Soil preparation for agricultural crops produces aerosols that may significantly contribute to seasonal atmospheric particulate matter (PM). Efforts to reduce PM emissions from tillage through a variety of conservation management practices (CMPs) have been made, but the reductions from many of these practices have not been measured in the field. A study was conducted in California's San Joaquin Valley to quantify emissions reductions from fall tillage CMP. Emissions were measured from conventional tillage methods and from a "combined operations" CMP, which combines several implements to reduce tractor passes. Measurements were made of soil moisture, bulk density, meteorological profiles, filter-based total suspended PM (TSP), concentrations of PM with an equivalent aerodynamic diameter ≤10 µm (PM) and PM with an equivalent aerodynamic diameter ≤2.5 µm (PM), and aerosol size distribution. A mass-calibrated, scanning, three-wavelength light detection and ranging (LIDAR) procedure estimated PM through a series of algorithms. Emissions were calculated via inverse modeling with mass concentration measurements and applying a mass balance to LIDAR data. Inverse modeling emission estimates were higher, often with statistically significant differences. Derived PM emissions for conventional operations generally agree with literature values. Sampling irregularities with a few filter-based samples prevented calculation of a complete set of emissions through inverse modeling; however, the LIDAR-based emissions dataset was complete. The CMP control effectiveness was calculated based on LIDAR-derived emissions to be 29 ± 2%, 60 ± 1%, and 25 ± 1% for PM, PM, and TSP size fractions, respectively. Implementation of this CMP provides an effective method for the reduction of PM emissions.

Particulate emissions from a beef cattle feedlot using the flux-gradient technique.

Data on air emissions from open-lot beef cattle () feedlots are limited. This research was conducted to determine fluxes of particulate matter with an aerodynamic diameter ≤10 µm (PM) from a commercial beef cattle feedlot in Kansas using the flux-gradient technique, a widely used micrometeorological method for air emissions from open sources. Vertical PM concentration profiles and micrometeorological parameters were measured at the feedlot using tapered element oscillating microbalance PM samplers and eddy covariance instrumentations (i.e., sonic anemometer and infrared hygrometer), respectively, from May 2010 through September 2011, representing feedlot conditions with air temperatures ranging from -24 to 39°C. Calculated hourly PM fluxes varied diurnally and seasonally, ranging up to 272 mg m h, with an overall median of 36 mg m h. For warm conditions (air temperature of 21 ± 10°C), the highest hourly PM fluxes (range 116-146 mg m h) were observed during the early evening period, from 2000 to 2100 h. For cold conditions (air temperature of -2 ± 8°C), the highest PM fluxes (range 14-27 mg m h) were observed in the afternoon, from 1100 to 1500 h. Changes in the hourly trend of PM fluxes coincided with changes in friction velocity, air temperature, sensible heat flux, and surface roughness. The PM emission was also affected by the pen surface water content, where a water content of at least 20% (wet basis) would be sufficient to effectively reduce PM emissions from pens by as much as 60%.

Comparison of AERMOD and WindTrax dispersion models in determining PM10 emission rates from a beef cattle feedlot.

UNLABELLED: Reverse dispersion modeling has been used to determine air emission fluxes from ground-level area sources, including open-lot beef cattle feedlots. This research compared Gaussian-based AERMOD, the preferred regulatory dispersion model of the US. Environmental Protection Agency (EPA), and WindTrax, a backward Lagrangian stochastic-based dispersion model, in determining PM10 emission rates for a large beef cattle feedlot in Kansas. The effect of the type of meteorological data was also evaluated. Meteorological conditions and PM10 concentrations at the feedlot were measured with micrometeorological/eddy covariance instrumentation and tapered element oscillating microbalance (TEOM) PM10 monitors, respectively, from May 2010 through September 2011. Using the measured meteorological conditions and assuming a unit emission flux (i.e., 1 microg/m2-sec), each model was used to calculate PM10 concentrations (referred to as unit-flux concentrations). PM10 emission fluxes were then back-calculated using the measured and calculated unit-flux PM10 concentrations. For AERMOD, results showed that the PM10 emission fluxes determined using the two different meteorological data sets evaluated (eddy covariance-derived and AERMET-generated) were basically the same. For WindTrax, the two meteorological data sets (sonic anemometer data set, a three-variable data set composed of wind parameters, surface roughness, and atmospheric stability) also produced basically the same PM10 emission fluxes. Back-calculated emission fluxes from AERMOD were 32 to 69% higher than those from WindTrax. IMPLICATIONS: This work compared the PM10 emission rates determined from a large commercial cattle feedlot in Kansas by reverse dispersion modeling using AERMOD and WindTrax. Emission fluxes derived from AERMOD were greater than those from WindTrax by mean factors of 1.3 to 1.6. Based on the high linearity observed between the two models, emission fluxes derived from one dispersion model for the purpose of simulating dispersion could be applied to the other model using appropriate conversion factors.

Comparison of field-scale herbicide runoff and volatilization losses: an eight-year field investigation.

An 8-yr study was conducted to better understand factors influencing year-to-year variability in field-scale herbicide volatilization and surface runoff losses. The 21-ha research site is located at the USDA-ARS Beltsville Agricultural Research Center in Beltsville, MD. Site location, herbicide formulations, and agricultural management practices remained unchanged throughout the duration of the study. Metolachlor [2-chloro--(2-ethyl-6-methylphenyl)--(2-methoxy-1-methylethyl) acetamide] and atrazine [6-chloro--ethyl--(1-methylethyl)-1,3,5-triazine-2,4-diamine] were coapplied as a surface broadcast spray. Herbicide runoff was monitored from a month before application through harvest. A flux gradient technique was used to compute volatilization fluxes for the first 5 d after application using herbicide concentration profiles and turbulent fluxes of heat and water vapor as determined from eddy covariance measurements. Results demonstrated that volatilization losses for these two herbicides were significantly greater than runoff losses ( < 0.007), even though both have relatively low vapor pressures. The largest annual runoff loss for metolachlor never exceeded 2.5%, whereas atrazine runoff never exceeded 3% of that applied. On the other hand, herbicide cumulative volatilization losses after 5 d ranged from about 5 to 63% of that applied for metolachlor and about 2 to 12% of that applied for atrazine. Additionally, daytime herbicide volatilization losses were significantly greater than nighttime vapor losses ( < 0.05). This research confirmed that vapor losses for some commonly used herbicides frequently exceeds runoff losses and herbicide vapor losses on the same site and with the same management practices can vary significantly year to year depending on local environmental conditions.

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.

Estimation of Evapotranspiration and Energy Fluxes using a Deep-Learning based High-Resolution Emissivity Model and the Two-Source Energy Balance Model with sUAS information.

Surface temperature is necessary for the estimation of energy fluxes and evapotranspiration from satellites and airborne data sources. For example, the Two-Source Energy Balance (TSEB) model uses thermal information to quantify canopy and soil temperatures as well as their respective energy balance components. While surface (also called kinematic) temperature is desirable for energy balance analysis, obtaining this temperature is not straightforward due to a lack of spatially estimated narrowband (sensor-specific) and broadband emissivities of vegetation and soil, further complicated by spectral characteristics of the UAV thermal camera. This study presents an effort to spatially model narrowband and broadband emissivities for a microbolometer thermal camera at UAV information resolution (~0.15 m) based on Landsat and NASA HyTES information using a deep learning (DL) model. The DL model is calibrated using equivalent optical Landsat / UAV spectral information to spatially estimate narrowband emissivity values of vegetation and soil in the 7-14-nm range at UAV resolution. The resulting DL narrowband emissivity values were then used to estimate broadband emissivity based on a developed narrowband-broadband emissivity relationship using the MODIS UCSB Emissivity Library database. The narrowband and broadband emissivities were incorporated into the TSEB model to determine their impact on the estimation of instantaneous energy balance components against ground measurements. The proposed effort was applied to information collected by the Utah State University AggieAir small Unmanned Aerial Systems (sUAS) Program as part of the ARS-USDA GRAPEX Project (Grape Remote sensing Atmospheric Profile and Evapotranspiration eXperiment) over a vineyard located in Lodi, California. A comparison of resulting energy balance component estimates, with and without the inclusion of high-resolution narrowband and broadband emissivities, against eddy covariance (EC) measurements under different scenarios are presented and discussed.

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.

Soil moisture and metolachlor volatilization observations over three years.

A 3-yr study was conducted to focus on the impact of surface soil water content on metolachlor (2-chloro-N-(2-ethyl-6-methylphenyl)-N-(2-methoxy-1-methylethyl) acetamide) volatilization from a field with different surface soil water regimes created by subsurface water flow paths. Metolachlor vapor fluxes were measured at two locations within the field where local meteorological and soil conditions were relatively constant, except for surface soil water content, which differed significantly. Surface soil water content at the two sites differed in response to the presence of subsurface flow pathways. Detailed soil moisture observations over the duration of the study showed that for the first 2 yr (2004 and 2005), surface soil water contents at the dry location (V1) were nearly half those at the wetter location (V2). Cumulative metolachlor vapor fluxes during 2004 and 2005 at V1 were also about half that at V2. In the third year (2006), early-season drought conditions rendered the soil water content at the two locations to be nearly identical, resulting in similar metolachlor volatilization losses. Analysis of infrared soil surface temperatures suggests a correlation between surface soil temperatures and metolachlor volatilization when soils are wet (2004 and 2005) but not when the soils are dry (2006). Field-averaged metolachlor volatilization losses were highly correlated with increasing surface soil water contents (r(2) = 0.995).

Solar radiation, relative humidity, and soil water effects on metolachlor volatilization.

Pesticide volatilization is a significant loss pathway that may have unintended consequences in nontarget environments. Field-scale pesticide volatilization involves the interaction of a number of complex variables. There is a need to acquire pesticide volatilization fluxes from a location where several of these variables can be held constant. Accordingly, soil properties, tillage practices, surface residue management, and pesticide formulations were held constant while fundamental information regarding metolachlor volatilization (a pre-emergent pesticide) was monitored over a five-year period as influenced by meteorological variables and soil water content. Metolachlor vapor concentrations were measured continuously for 120 h after each application using polyurethane foam plugs in a logarithmic profile above the soil surface. A flux gradient technique was used to compute volatilization fluxes from metolachlor concentration profiles and turbulent fluxes of heat and water vapor (as determined from eddy covariance measurements). Differences in meteorological conditions and surface soil water contents resulted in variability of the volatilization losses over the years studied. The peak volatilization losses for each year occurred during the first 24 h after application with a maximum flux rate in 2001 (1500 ng m(-2) s(-1)) associated with wet surface soil conditions combined with warm temperatures. The cumulative volatilization losses for the 120-hour period following metolachlor application varied over the years from 5 to 25% of the applied active ingredient, with approximately 87% of the losses occurring during the first 72 h. In all of the years studied, volatilization occurred diurnally and accounted for between 43 and 86% during the day and 14 and 57% during the night of the total measured loss. The results suggest that metolachlor volatilization is influenced by multiple factors involving meteorological, surface soil, and chemical factors.