The methods for accurately fusing medium- and high-spatial-resolution satellite reflectance are vital for monitoring vegetation biomass, agricultural irrigation, ecological processes and climate change. However, the currently existing fusion methods cannot accurately capture the temporal variation in reflectance for heterogeneous landscapes. In this study, we proposed a new method, the spatial and temporal reflectance fusion method based on the unmixing theory and a fuzzy C-clustering model (FCMSTRFM), to generate Landsat-like time-series surface reflectance. Unlike other data fusion models, the FCMSTRFM improved the similarity of pixels grouped together by combining land cover maps and time-series data cluster algorithms to define endmembers. The proposed method was tested over a 2000 km² study area in Heilongjiang Provence, China, in 2017 and 2018 using ten images. The results show that the accuracy of the FCMSTRFM is better than that of the popular enhanced spatial and temporal adaptive reflectance fusion model (ESTARFM) (correlation coefficient (R): 0.8413 vs. 0.7589; root mean square error (RMSE): 0.0267 vs. 0.0401) and the spatial-temporal data fusion approach (STDFA) (R: 0.8413 vs. 0.7666; RMSE: 0.0267 vs. 0.0307). Importantly, the FCMSTRFM was able to maintain the details of temporal variations in complicated landscapes. The proposed method provides an alternative method to monitor the dynamics of land surface variables over complicated heterogeneous regions.

An accurate and spatially continuous estimation of terrestrial latent heat flux (LE) is fundamental and crucial for the rational utilization of water resources in the Haihe River Basin (HRB). However, the sparsity of flux observation sites hinders the accurate characterization of spatiotemporal LE patterns over the HRB. In this study, we estimated the daily LE across the HRB using the gradient boosting regression tree (GBRT) from global land surface satellite NDVI data, reanalysis data and eddy covariance data. Compared with the random forests (RF) and extra tree regressor (ETR) methods, the GBRT obtains the best results, with R² = 0.86 and root mean square error (RMSE = 18.1 W/m². Then, we applied the GBRT algorithm to map the average annual terrestrial LE of the HRB from 2016 to 2018 with a spatial resolution of 0.05°. When compared with the Global Land Surface Satellite (GLASS) and Moderate Resolution Imaging Spectroradiometer (MODIS) LE products, the difference between the terrestrial LE estimated by the GBRT algorithm and the GLASS and MODIS products was less than 20 W/m² in most areas; thus, the GBRT algorithm was reliable and reasonable for estimating the long-term LE estimation over the HRB.