Satellite land surface temperature (LST) is vital for climatological and environmental studies. However, LST datasets are not continuous in time and space mainly due to cloud cover. Here we combine LST with Climate Forecast System Version 2 (CFSv2) modeled temperatures to derive a continuous gap filled global LST dataset at a spatial resolution of 1 km. Temporal Fourier analysis is used to derive the seasonality (climatology) on a pixel-by-pixel basis, for LST and CFSv2 temperatures. Gaps are filled by adding the CFSv2 temperature anomaly to climatological LST. The accuracy is evaluated in nine regions across the globe using cloud-free LST (mean values: R2 = 0.93, Root Mean Square Error (RMSE) = 2.7 °C, Mean Absolute Error (MAE) = 2.1 °C). The provided dataset contains day, night, and daily mean LST for the Eastern Mediterranean. We provide a Google Earth Engine code and a web app that generates gap filled LST in any part of the world, alongside a pixel-based evaluation of the data in terms of MAE, RMSE and Pearson’s r.