Reconstruction of land surface temperature under cloudy conditions from Landsat 8 data using annual temperature cycle model

Land surface temperature (LST) is an important parameter in the processes of energy exchange and water cycle between the land surface and the atmosphere. The impact of cloud cover leads to spatially incomplete of thermal infrared (TIR)-based LST products, which seriously hinders the applications of LST products in various fields. Several methods have been developed to reconstruct LST under cloudy conditions in previous studies, but there is a lack of an effective method for the reconstruction of cloudy LST at the spatial resolution of Landsat pixel (30 m). In this study, a novel method was proposed to reconstruct LST under cloudy conditions from Landsat 8 data. The LST reconstruction method includes four main steps: (1) identification of cloud-free, cloud-shadow, cloud-obscured, and cloud-covered pixels by integrating the Fmask method with a cloud-shape matching method; (2) calculation of annual temperature cycle (ATC)-based reference LST by fitting an ATC model to all available Landsat 8 LST product during 2013-2020; (3) estimation of LST residual from spatially adjacent similar pixels; and (4) estimation of reconstructed LST in terms of the sum of ATC-based reference LST and LST residual. The performance of the LST reconstruction method was evaluated using Landsat 8 LST images under clear-sky conditions as reference data. The root mean squared error (RMSE) between reconstructed LST and Landsat 8 reference LST ranges from 0.9 K to 2.5 K. The LST reconstruction method was further applied to reconstruct actual Landsat 8 LST images under cloudy conditions. Compared with original Landsat 8 LST images, the spatial distribution of reconstructed LST images is more complete. The pattern of reconstructed LST images reflects the spatial variability of LST well. The accuracy of the LST reconstruction method was validated against in situ LST measurements at six SURFRAD (Surface Radiation Budget Network) sites. The overall bias and RMSE between reconstructed LST and in situ LST at all sites are approximately −0.3 K and 3.5 K, respectively. The LST reconstruction method has great potentials to improve the applications of Landsat LST product in urban thermal environment monitoring and crop water stress monitoring.