Abstract Land surface temperature (LST) is an important factor in land monitoring studies, but due to the presence of clouds, dust and sensor issues, there are missing values. The aims of this research are to determine the optimal parameters for the reconstruction of Landsat-LST images, required in many applications, by the harmonic analysis of time series algorithm (HANTS) and to investigate the possibility of improving LST reconstruction accuracy using Landsat 8 and 9 images simultaneously. For these aims, 91 Landsat 8 and 9 images with 100 m spatial resolution in 2022 and 2023 are employed, covering Yazd-Ardakan plain in Iran. Three methods are used for evaluation. In method one, a part of LST image is considered as a gap and is compared with the initial value after reconstruction. In method two, on a cloudy day and a cloudless day, surface temperature values are measured using thermometers at fifty points in plain lands, and the difference between gap-filled satellite measurements and ground measurements is calculated. In method three, all the reconstructed LST images are compared with the original images. In method one, the root mean square error (RMSE) of reconstructed LST reduces by 1.3°C when using the combined Landsat 8 and 9 images. In method two, RMSEs of reconstructed LST images are 6.1°C when using Landsat 8 and 5.4°C when using the combined Landsat 8 and 9. Method three shows that 41% of the study region has RMSE of less than 2°C when using only Landsat 8, while this value becomes 72% when combining Landsat 8 and 9. In general, the combined use of Landsat 8 and 9 LST images improves the accuracy of reconstruction using HANTS. The findings of this research are crucial for regional applications and remote monitoring of surface temperature in areas with limited weather stations.