Accurate and frequent monitoring of evapotranspiration (ET) at sub-field scales can provide valuable information for agricultural water management, quantifying crop water use and stress toward the goal of increasing crop water use efficiency and production. Using land-surface temperature (LST) data retrieved from Landsat thermal infrared (TIR) imagery, along with surface reflectance data describing albedo and vegetation cover fraction, surface energy balance models can generate ET maps down to a 30 m spatial resolution. However, the temporal sampling by such maps can be limited by the relatively infrequent revisit period of Landsat data (8 days for combined Landsats 7 and 8), especially in cloudy areas experiencing rapid changes in moisture status. The Sentinel-2 (S2) satellites, as a good complement to the Landsat system, provide surface reflectance data at 10–20 m spatial resolution and 5 day revisit period but do not have a thermal sensor. On the other hand, the Visible Infrared Imaging Radiometer Suite (VIIRS) provides TIR data on a near-daily basis with 375 m resolution, which can be refined through thermal sharpening using S2 reflectances. This study assesses the utility of augmenting the Harmonized Landsat and Sentinel-2 (HLS) dataset with S2-sharpened VIIRS as a thermal proxy source on S2 overpass days, enabling 30 m ET mapping at a potential combined frequency of 2–3 days (including Landsat). The value added by including VIIRS-S2 is assessed both retrospectively and operationally in comparison with flux tower observations collected from several U.S. agricultural sites covering a range of crop types. In particular, we evaluate the performance of VIIRS-S2 ET estimates as a function of VIIRS view angle and cloud masking approach. VIIRS-S2 ET retrievals (MAE of 0.49 mm d⁻¹ against observations) generally show comparable accuracy to Landsat ET (0.45 mm d⁻¹) on days of commensurate overpass, but with decreasing performance at large VIIRS view angles. Low-quality VIIRS-S2 ET retrievals linked to imperfect VIIRS/S2 cloud masking are also discussed, and caution is required when applying such data for generating ET timeseries. Fused daily ET time series benefited during the peak growing season from the improved multi-source temporal sampling afforded by VIIRS-S2, particularly in cloudy regions and over surfaces with rapidly changing vegetation conditions, and value added for real-time monitoring applications is discussed. This work demonstrates the utility and feasibility of augmenting the HLS dataset with sharpened VIIRS TIR imagery on S2 overpass dates for generating high spatiotemporal resolution ET products.

Evapotranspiration (ET) represents crop water use and is a key indicator of crop health. Accurate estimation of ET is critical for agricultural irrigation and water resource management. ET retrieval using energy balance methods with remotely sensed thermal infrared data as the key input has been widely applied for irrigation scheduling, yield prediction, drought monitoring and so on. However, limitations on the spatial and temporal resolution of available thermal satellite data combined with the effects of cloud contamination constrain the amount of detail that a single satellite can provide. Fusing satellite data from different satellites with varying spatial and temporal resolutions can provide a more continuous estimation of daily ET at field scale. In this study, we applied an ET fusion modeling system, which uses a surface energy balance model to retrieve ET using both Landsat and Moderate Resolution Imaging Spectroradiometer (MODIS) data and then fuses the Landsat and MODIS ET retrieval timeseries using the Spatial-Temporal Adaptive Reflectance Fusion Model (STARFM). In this paper, we compared different STARFM ET fusion implementation strategies over various crop lands in the central California. In particular, the use of single versus two Landsat-MODIS pair images to constrain the fusion is explored in cases of rapidly changing crop conditions, as in frequently harvested alfalfa fields, as well as an improved dual-pair method. The daily 30 m ET retrievals are evaluated with flux tower observations and analyzed based on land cover type. This study demonstrates improvement using the new dual-pair STARFM method compared with the standard one-pair STARFM method in estimating daily field scale ET for all the major crop types in the study area.