Satellite-based land surface temperature (LST) products are frequently affected by cloud cover and atmospheric conditions, resulting in missing data that significantly limits the continuous monitoring of the thermal environment in complex terrains, such as the Tibetan Plateau. Existing spatiotemporal interpolation methods face clear accuracy limitations when addressing extensive data gaps, while physical models often struggle due to insufficient meteorological inputs in complex landscapes. Moreover, conventional data-driven approaches usually overlook local spatial variations, resulting in smoothed thermal patterns and systematic errors. To overcome these issues, we propose a Physically Constrained Spatial Residual Learning framework. In this framework, we use the Enhanced Annual Temperature Cycle (EATC) model to capture the temporal baseline of LST first. Then, we integrate multi-source auxiliary data into the Geographical-XGBoost (G-XGBoost) algorithm to model spatial nonlinear residuals. Using simulated cloud masks on the 2017 MODIS LST dataset from the Qinghai–Tibet Engineering Corridor, we show that the hybrid model outperforms both individual physical models and global machine learning models in accuracy and spatial detail recovery. Validation results yield an R2 of 0.88, an RMSE of 1.92 K, and a mean bias of 0.07 K. Seasonal evaluations indicate best performance in winter (RMSE = 1.19 K) with robust performance in summer. Furthermore, the framework reduces boundary artifacts and accurately reproduces thermal spatial patterns in complex terrain through adaptive local bandwidth and weight adjustments. This approach provides a reliable method for high-precision LST reconstruction over heterogeneous alpine surfaces.
Xu et al. (Tue,) studied this question.