The demand for freshwater is increasing with population growth and rapid socio-economic development. It is more and more important for refined irrigation water management to conduct research on crop evapotranspiration (ET) data with a high spatiotemporal resolution in agricultural regions. We propose the unmixing–weight ET image fusion model (UWET), which integrates the advantages of the unmixing method in spatial downscaling and the weight-based method in temporal prediction to produce daily ET maps with a high spatial resolution. The Landsat-ET and MODIS-ET datasets for the UWET fusion data are retrieved from Landsat and MODIS images based on the surface energy balance model. The UWET model considers the effects of crop phenology, precipitation, and land cover in the process of the ET image fusion. The precision evaluation is conducted on the UWET results, and the measured ET values are monitored by eddy covariance at the Luancheng station, with average MAE values of 0.57 mm/day. The image results of UWET show fine spatial details and capture the dynamic ET changes. The seasonal ET values of winter wheat from the ET map mainly range from 350 to 660 mm in 2019–2020 and from 300 to 620 mm in 2020–2021. The average seasonal ET in 2019–2020 is 499.89 mm, and in 2020–2021, it is 459.44 mm. The performance of UWET is compared with two other fusion models: the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) and the Spatial and Temporal Reflectance Unmixing Model (STRUM). UWET performs better in the spatial details than the STARFM and is better in the temporal characteristics than the STRUM. The results indicate that UWET is suitable for generating ET products with a high spatial–temporal resolution in agricultural regions.
High-spatiotemporal-resolution evapotranspiration (ET) is critical for precision irrigation management and water resource regulation. Regarding the existing spatiotemporal fusion methods suffering from sparse high-resolution observations and coarse land surface temperature (LST), this study took winter wheat in Luancheng District, Hebei Province, as the research object, and proposed a remote sensing ET retrieval method based on the LST sharpening model. The Data Mining Sharpener (DMS) algorithm combined with Sentinel-2 multispectral data was used to downscale MODIS LST from 1000 m to 10 m, with auxiliary variables (DEM, albedo, NDVI, land cover) integrated into the Cubist regression tree to improve the physical rationality and spatial details of MODIS LST. The 10 m resolution ET was estimated from 10 m sharpened LST and Sentinel-2 multispectral data using the surface energy balance model, and the unmixing–weight ET image fusion model (UWET) was adopted to fuse the 10 m resolution ET with MODIS low-resolution ET to generate a daily 10 m ET dataset covering the entire winter wheat growing season. Validation with eddy covariance flux measurements showed that the correlation coefficient R = 0.921, RMSE = 0.779 mm/day during 2019–2020, and R = 0.900, RMSE = 0.831 mm/day during 2020–2021. The results demonstrate that auxiliary variables significantly enhance the spatial reality of LST, LST sharpening effectively improves the spatial heterogeneity of ET, and Sentinel-2 data compensates for the temporal deficiency of Landsat, thereby greatly promoting the accuracy of spatiotemporal fusion. This method can provide reliable high-spatiotemporal-resolution data support for refined farmland irrigation management and water resources regulation.