One of the primary applications of satellite Land Surface Temperature (LST) observations lies in their utilization for modeling of actual evapotranspiration (ET) in agricultural crops, with the primary goals of monitoring and enhancing irrigation practices and improving crop water use productivity, as stipulated by Sustainable Development Goal (SDG) indicator 6.4.1. Evapotranspiration is a complex and dynamic process, both temporally and spatially, necessitating LST observations with high spatio-temporal resolution. Presently, none of the existing spaceborne thermal sensors can provide quasi-daily field-scale LST observations, prompting the development of methods for data fusion (thermal sharpening) of observations from various shortwave and thermal sensors to meet this spatio-temporal requirement. Previous research has demonstrated the effectiveness of combining shortwave-multispectral Sentinel-2 observations with thermal-infrared Sentinel-3 observations to derive daily, field-scale LST and ET estimates. However, these studies also highlighted limitations in capturing the distinct thermal contrast between cooler LST in irrigated agricultural areas and the hotter, adjacent dry regions. In this study, we aim to address this limitation by incorporating information on thermal spatial variability observed by Landsat satellites into the data fusion process, without being constrained by infrequent or cloudy Landsat thermal observations and while retaining the longwave radiance emission captured by the Sentinel-3 thermal sensor at its native resolution. Two approaches are evaluated, both individually and as a complementary combination, and validated against in situ LST measurements. The best performing approach, which leads to reduction in root mean square error of up to 1.5 K when compared to previous research, is subsequently used to estimate parcel-level actual evapotranspiration. The ET modeling process has also undergone various improvements regarding the gap-filling