Modeling lateral carbon fluxes for agroecosystems in the Mid-Atlantic region: Control factors and importance for carbon budget.

Estimating lateral carbon fluxes in agroecosystems presents challenges due to intricate anthropogenic and biophysical interactions. We used a modeling technique to enhance our comprehension of the determinants influencing lateral carbon fluxes and their significance in agroecosystem carbon budgets. The SWAT-C model was refined by incorporating a dynamic dissolved inorganic carbon (DIC) module, enhancing our ability to accurately quantify total lateral carbon fluxes. This improved model was calibrated using observed data on riverine particulate organic carbon (POC) and dissolved organic carbon (DOC) fluxes, as well as net ecosystem exchange (NEE) data monitored by a flux tower situated in a representative agricultural watershed, the Tuckahoe Watershed (TW) of the Chesapeake Bay's coastal plain. We assessed the losses of POC, DOC, and DIC across five primary rotation types: C (continuous carbon), CS (corn-soybean), CSS (corn-soybean-soybean), CWS (corn-wheat-soybean), and CWSCS (corn-wheat-soybean-corn-soybean). Our study revealed notable variations in the average annual fluxes of POC (ranging between 152 and 198 kg ha-1), DOC (74-85 kg ha-1), and DIC (93-156 kg ha-1) across the five rotation types. The primary influencing factor for annual POC fluxes was identified as sediment yield. While both annual DOC and DIC fluxes displayed a marked correlation with surface runoff across all crop rotation schemes, soil respiration also significantly influenced annual DIC fluxes. Total lateral carbon fluxes (POC + DOC+DIC) constituted roughly 11 % of both net ecosystem production (NEP) and NEE, yet they represented a striking 95 % of net biome production (NBP) in the TW's agroecosystem. Grain yield carbon accounted for 80 % of both NEP and NEE and was nearly seven times that of NBP. Our findings suggest that introducing soybeans into cornfields tends to reduce NEP, NEE, and also NBP. Conversely, integrating winter wheat into the corn-soybean rotation significantly boosted NEP, NEE, and NBP values, with NBP even surpassing the levels in continuous corn cultivation.

Effectiveness and diurnal variations of vegetative environmental buffers (VEBs) for mitigating NH3 and PM emissions from poultry houses.

Air pollutants from poultry production, such as ammonia (NH3) and particulate matter (PM), have raised concerns due to their potential negative impacts on human health and the environment. Vegetative environmental buffers (VEBs), consisting of trees and/or grasses planted around poultry houses, have been investigated as a mitigation strategy for these emissions. Although previous research demonstrated that VEBs can reduce NH3 and PM emissions, these studies used a limited number of samplers and did not examine concentration profiles. Moreover, the differences between daytime and nighttime emissions have not been investigated. In this study, we characterized emission profiles from a commercial poultry house using an array with multiple sampling heights and explored the differences between daytime and nighttime NH3 and PM profiles. We conducted three sampling campaigns, each with ten sampling events (five daytime and five nighttime), at a VEB-equipped poultry production facility. NH3 and PM samples were collected downwind from the ventilation tunnel fans before, within, and after the VEB. Results showed that ground-level concentrations beyond the VEB decreased to 8.0% ± 2.7% for NH3, 13% ± 4% for TSP, 13% ± 4% for PM10, and 2.4% ± 2.8% for PM2.5 of the original concentrations from the exhaust tunnel fan, with greater reduction efficiency during daytime than nighttime. Furthermore, pollutant concentrations were positively intercorrelated. These findings will be valuable for developing more effective pollutant remediation strategies in poultry house emissions.

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.

Evaluation of satellite Leaf Area Index in California vineyards for improving water use estimation.

Remote sensing estimation of evapotranspiration (ET) directly quantifies plant water consumption and provides essential information for irrigation scheduling, which is a pressing need for California vineyards as extreme droughts become more frequent. Many ET models take satellite-derived Leaf Area Index (LAI) as a major input, but how uncertainties of LAI estimations propagate to ET and the partitioning between evaporation and transpiration is poorly understood. Here we assessed six satellite-based LAI estimation approaches using Landsat and Sentinel-2 images against ground measurements from four vineyards in California and evaluated ET sensitivity to LAI in the thermal-based two-source energy balance (TSEB) model. We found that radiative transfer modeling-based approaches predicted low to medium LAI well, but they significantly underestimated high LAI in highly clumped vine canopies (RMSE ~ 0.97 to 1.27). Cubist regression models trained with ground LAI measurements from all vineyards achieved high accuracy (RMSE ~ 0.3 to 0.48), but these empirical models did not generalize well between sites. Red edge bands and the related vegetation index (VI) from the Sentinel-2 satellite contain complementary information of LAI to VIs based on near-infrared and red bands. TSEB ET was more sensitive to positive LAI biases than negative ones. Positive LAI errors of 50% resulted in up to 50% changes in ET, while negative biases of 50% in LAI caused less than 10% deviations in ET. However, even when ET changes were minimal, negative LAI errors of 50% led to up to a 40% reduction in modeled transpiration, as soil evaporation and plant transpiration responded to LAI change divergently. These findings call for careful consideration of satellite LAI uncertainties for ET modeling, especially for the partitioning of water loss between vine and soil or cover crop for effective vineyard irrigation management.

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.

Scaling photosynthetic function and CO2 dynamics from leaf to canopy level for maize - dataset combining diurnal and seasonal measurements of vegetation fluorescence, reflectance and vegetation indices with canopy gross ecosystem productivity.

Recent advances in leaf fluorescence measurements and canopy proximal remote sensing currently enable the non-destructive collection of rich diurnal and seasonal time series, which are required for monitoring vegetation function at the temporal and spatial scales relevant to the natural dynamics of photosynthesis. Remote sensing assessments of vegetation function have traditionally used actively excited foliar chlorophyll fluorescence measurements, canopy optical reflectance data and vegetation indices (VIs), and only recently passive solar induced chlorophyll fluorescence (SIF) measurements. In general, reflectance data are more sensitive to the seasonal variations in canopy chlorophyll content and foliar biomass, while fluorescence observations more closely relate to the dynamic changes in plant photosynthetic function. With this dataset we link leaf level actively excited chlorophyll fluorescence, canopy proximal reflectance and SIF, with eddy covariance measurements of gross ecosystem productivity (GEP). The dataset was collected during the 2017 growing season on maize, using three automated systems (i.e., Monitoring Pulse-Amplitude-Modulation fluorimeter, Moni-PAM; Fluorescence Box, FloX; and from eddy covariance tower). The data were quality checked, filtered and collated to a common 30 minutes timestep. We derived vegetation indices related to canopy functioning (e.g., Photochemical Reflectance Index, PRI; Normalized Difference Vegetation Index, NDVI; Chlorophyll Red-edge, Clre) to investigate how SIF and VIs can be coupled for monitoring vegetation photosynthesis. The raw datasets and the filtered and collated data are provided to enable new processing and analyses.

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.

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.

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/m2 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.

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.

Data Assimilation of High-Resolution Thermal and Radar Remote Sensing Retrievals for Soil Moisture Monitoring in a Drip-Irrigated Vineyard.

Efficient water use assessment and irrigation management is critical for the sustainability of irrigated agriculture, especially under changing climate conditions. Due to the impracticality of maintaining ground instrumentation over wide geographic areas, remote sensing and numerical model-based fine-scale mapping of soil water conditions have been applied for water resource applications at a range of spatial scales. Here, we present a prototype framework for integrating high-resolution thermal infrared (TIR) and synthetic aperture radar (SAR) remote sensing data into a soil-vegetation-atmosphere-transfer (SVAT) model with the aim of providing improved estimates of surface- and root-zone soil moisture that can support optimized irrigation management strategies. Specifically, remotely-sensed estimates of water stress (from TIR) and surface soil moisture retrievals (from SAR) are assimilated into a 30-m resolution SVAT model over a vineyard site in the Central Valley of California, U.S. The efficacy of our data assimilation algorithm is investigated via both the synthetic and real data experiments. Results demonstrate that a particle filtering approach is superior to an ensemble Kalman filter for handling the nonlinear relationship between model states and observations. In addition, biophysical conditions such as leaf area index are shown to impact the relationship between observations and states and must therefore be represented accurately in the assimilation model. Overall, both surface and root-zone soil moisture predicted via the SVAT model are enhanced through the assimilation of thermal and radar-based retrievals, suggesting the potential for improving irrigation management at the agricultural sub-field scale using a data assimilation strategy.

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.

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.

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.