Search

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 informationThe online version contains supplementary material available at 10.1007/s00271-022-00808-9.

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 informationThe online version contains supplementary material available at 10.1007/s00271-022-00787-x.

Our ability to measure and map plant function at multiple ecological scales is critical for understanding current and future changes in Earth's ecosystems and the global carbon budget. Conventional plant functional types (cPFTs) based on a few productivity-related traits have been previously used to simplify and represent major differences in global plant functions, but more recent research has directly focused on the use of functional trait information. Still, sampling limitations have constrained efforts to truly understand the variance and covariance of functional traits globally. Reflectance spectra offer a fast, repeatable, simultaneous measurement of a wide variety of leaf functional traits and could be used to optically define leaf functional types. To evaluate this concept, we measured leaf reflectance from a wide range of species in a diverse set of ecosystems across central and northern California, including observations from multiple individuals, sites, and seasons. Using principal components analysis, we analyzed spectral variation in relation to categorical attributes such as species and cPFTs, as well as to a set of functional trait metrics calculated from the spectra. We found the first three principal components (PCs) to be weakly related to categorical attributes and more strongly related to spectrally-derived functional metrics. Each PC was more strongly associated with different portions of the spectrum and contained different functional information. We applied a hybrid clustering algorithm to the PC coordinates of the observations to define potential optical leaf functional types. Twelve spectral clusters were identified, and these did not correspond directly to either single cPFTs or species. However, each cluster had a unique functional metric profile. Clusters represented both inter- and intra-species and cPFT functional differences driven by taxonomy, trait evolution and environmental responses, demonstrating their value as optical leaf functional types and the value of the clustering approach used here for defining optical types from leaf spectra. Our findings support the notion that cPFTs do not adequately capture differences in leaf function. They demonstrate that spectral measurements can be used to improve both the definition of PFTs as well as our knowledge regarding the covariance of functional traits within these classes.

The role of root systems in drought tolerance is a subject of very limited information compared with above-ground responses. Adjustments to the ability of roots to supply water relative to shoot transpiration demand is proposed as a major means for woody perennial plants to tolerate drought, and is often expressed as changes in the ratios of leaf to root area (A(L):A(R)). Seasonal root proliferation in a directed manner could increase the water supply function of roots independent of total root area (A(R)) and represents a mechanism whereby water supply to demand could be increased. To address this issue, seasonal root proliferation, stomatal conductance (g(s)) and whole root system hydraulic conductance (k(r)) were investigated for a drought-tolerant grape root system (Vitis berlandieri×V. rupestris cv. 1103P) and a non-drought-tolerant root system (Vitis riparia×V. rupestris cv. 101-14Mgt), upon which had been grafted the same drought-sensitive clone of Vitis vinifera cv. Merlot. Leaf water potentials (ψ(L)) for Merlot grafted onto the 1103P root system (-0.91±0.02 MPa) were +0.15 MPa higher than Merlot on 101-14Mgt (-1.06±0.03 MPa) during spring, but dropped by approximately -0.4 MPa from spring to autumn, and were significantly lower by -0.15 MPa (-1.43±0.02 MPa) than for Merlot on 101-14Mgt (at -1.28±0.02 MPa). Surprisingly, g(s) of Merlot on the drought-tolerant root system (1103P) was less down-regulated and canopies maintained evaporative fluxes ranging from 35-20 mmol vine(-1) s(-1) during the diurnal peak from spring to autumn, respectively, three times greater than those measured for Merlot on the drought-sensitive rootstock 101-14Mgt. The drought-tolerant root system grew more roots at depth during the warm summer dry period, and the whole root system conductance (k(r)) increased from 0.004 to 0.009 kg MPa(-1) s(-1) during that same time period. The changes in k(r) could not be explained by xylem anatomy or conductivity changes of individual root segments. Thus, the manner in which drought tolerance was conveyed to the drought-sensitive clone appeared to arise from deep root proliferation during the hottest and driest part of the season, rather than through changes in xylem structure, xylem density or stomatal regulation. This information can be useful to growers on a site-specific basis in selecting rootstocks for grape clonal material (scions) grafted to them.

The magnitude of chlorophyll fluorescence emission represents both chlorophyll content and energy quenching processes enabling its application to serve as a proxy of photosynthetic activity. Thus, there is interest in advancing methods for canopy‐scale monitoring of chlorophyll fluorescence. Remotely sensed solar‐induced fluorescence (SIF) retrievals offer daytime monitoring of chlorophyll fluorescence, which can serve as an indicator of photosynthesis. However, it represents an instantaneous measurement during the day, which is strongly influenced by incoming radiation, solar angle, and sun/shade fraction—making it difficult to tease out baseline information on plant health and potential photosynthetic capacity—which could be tracked by changes in fluorescence yield (independent of sunlight). Recent advances have demonstrated the potential for inducing nighttime chlorophyll fluorescence via LED light sources at the canopy‐scale, which can be retrieved as LED‐induced chlorophyll fluorescence (LEDIF), potentially serving as a baseline indicator of plant health and photochemical capacity, independent of daytime conditions. In this study, we explored two methods of LEDIF retrievals: (a) hyperspectral sensor (1.33 nm full‐width half max) and (b) low‐cost Red‐Far‐Red photodiode sensor. LEDIF retrieved by the hyperspectral sensor demonstrated strong correlations with daytime SIF and gross primary productivity during mid to end of season phenology (R2 > 0.70). In contrast, phenological dynamics of LEDIF retrieved by the photodiode sensor was more subtle, likely due to weaker signal‐to‐noise ratio, but still demonstrated some potential. Overall, LEDIF offers a technique to monitor nighttime chlorophyll fluorescence emissions (and changes in its spectral shape with a hyperspectral sensor) to assess canopy‐scale phenology of photosynthetic potential. Plain Language Summary: Monitoring the baseline health status of plants is difficult because of their dynamic responses to the environment over the course of a day to sunlight, temperature, and water availability. Measuring plants at night offers an alternative method to assess a plant's baseline health status. Here we use blue LED lights at night to induce the reemissions of light from leaves via a red‐far‐red chlorophyll fluorescence signal detected by a low‐cost system. We found that this signal closely tracks the seasonal variation of photosynthesis. Because this signal was measured at night, unconfounded by daytime variation of sunlight, temperature, and water availability, our low‐cost system offers an approach to assess a plant's photosynthetic capacity and baseline health. Key Points: Low‐cost hyperspectral and red‐far‐red photodiode sensors with blue LED lights offer nighttime retrieval of chlorophyll fluorescence (LEDIF)LEDIF can track photosynthetic phenology similar to tower‐based solar‐induced fluorescence and gross primary productivityLEDIF may serve as a more direct assessment of baseline plant status unaffected by incoming radiation [ABSTRACT FROM AUTHOR]

Accurate characterization of evapotranspiration (ET) is imperative in water-limited cropping systems such as California vineyards and almond orchards. Satellite-based ET modeling techniques, including the atmosphere–land exchange inverse model (ALEXI) and associated flux disaggregation technique (DisALEXI), have proven reliable in determining field scale ET. However, validation efforts typically focus on ET and omit an evaluation of partitioned evaporation (E) and transpiration (T). ALEXI/DisALEXI is based on the two-source energy balance (TSEB) model, making it uniquely qualified to derive E and T individually. The current study evaluated E and T estimates derived using two formulations of DisALEXI; one based on Priestley-Taylor (DisALEXI-PT) and the other on Penman-Monteith (DisALEXI-PM). The modeled values were validated against partitioned fluxes derived from the conditional eddy covariance (CEC) approach using EC flux towers in three wine grape vineyards and three almond orchards for the year 2021. Modeled estimates were derived using Landsat 8 Collection 2 thermal infrared and surface reflectance imagery as well as Harmonized Landsat and Sentinel-2 surface reflectance datasets as input into DisALEXI. The results indicated that the modeled total ET fluxes were similar between the two methods, but the partitioned values diverged, with DisALEXI-PT overestimating E and slightly underestimating T when compared to CEC estimates. Conversely, DisALEXI-PM agreed better with CEC-derived E and overestimated T estimates under non-advective conditions. Compared to one another, DisALEXI-PM estimated canopy temperatures ~5 °C cooler and soil temperatures ~5 °C warmer than DisALEXI-PT, causing differences in E and T of −2.6 mm day−1 and +2.6 mm day−1, respectively. The evaluation of the iterative process required for DisALEXI indicates DisALEXI-PM ET values converge on ALEXI ET with proportionate adjustments to E and T, while DisALEXI-PT convergence is driven by adjustments to E. The analysis presented here can potentially drive improvements in the modeling framework to provide specific soil and canopy consumptive water use information in unique canopy structures, allowing for improved irrigation and water use efficiencies in these water-limited systems. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In viticulture, deficit irrigation strategies are often implemented to control vine canopy growth and to impose stress at critical stages of vine growth to improve wine grape quality. To support deficit irrigation scheduling, remote sensing technologies can be employed in the mapping of evapotranspiration (ET) at the field to sub-field scales, quantifying time-varying vineyard water requirements and actual water use. In the current study, we investigate the utility of ET maps derived from thermal infrared satellite imagery over a vineyard in the Central Valley of California equipped with a variable rate drip irrigation (VRDI) system which enables differential water applications at the 30 × 30 m scale. To support irrigation management at that scale, we utilized a thermal-based multi-sensor data fusion approach to generate weekly total actual ET (ETa) estimates at 30 m spatial resolution, coinciding with the resolution of the Landsat reflectance bands. Crop water requirements (ETc) were defined with a vegetative index (VI)-based approach. To test capacity to capture stress signals, the vineyard was sub-divided into four blocks with different irrigation management strategies and goals, inducing varying degrees of stress during the growing season. Results indicate derived weekly total ET from the thermal-based data fusion approach match well with observations. The thermal-based method was also able to capture the spatial heterogeneity in ET over the vineyard due to a water stress event imposed on two of the four vineyard blocks. This transient stress event was not reflected in the VI-based ETc estimate, highlighting the value of thermal band imaging. While the data fusion system provided valuable information, latency in current satellite data availability, particularly from Landsat, impacts operational applications over the course of a growing season. [ABSTRACT FROM AUTHOR]

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.

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.

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 (T s and T c ) to estimate surface energy fluxes including R n , H, LE, and G. The estimation of T s and T c 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 T s and T c 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 NDVI s and NDVI c . The observed relative error of LE corresponding to NDVI s uncertainty was between -1% and 2%, while for NDVI c uncertainty, the relative error was between -2.2% and 1.2%. However, when the combined NDVI s and NDVI c uncertainties were used simultaneously, the domain of the observed relative error corresponding to the absolute values of |ΔLE| was between 0% and 4%., Competing Interests: Conflicts of Interest: The authors declare no conflict of interest.

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 (u star ), 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., Competing Interests: Conflicts of Interest: The authors declare no conflict of interest.