Cloudy-sky land surface temperature from VIIRS and MODIS satellite data using a surface energy balance-based method

Land surface temperature (LST) has been effectively retrieved from thermal infrared (TIR) satellite measurements under clear-sky conditions. However, TIR satellite data are often severely contaminated by clouds, which cause spatiotemporal discontinuities and low retrieval accuracy in the LST products. Several solutions have been proposed to fill the "gaps"; however, a majority of these possess constraints. For example, fusion methods with microwave data suffer from coarse spatial resolution and diverse land cover types while spatial-temporal interpolation methods neglect cloudy cooling effects. We developed a novel method to estimate cloudy-sky LST from polar-orbiting satellite data based on the surface energy balance (SEB) principle. First, the hypothetical clear-sky LST of missing or likely cloud-contaminated pixels was reconstructed by assimilating high-quality satellite retrievals into a time-evolving model built from reanalysis data using a Kalman filter data assimilation algorithm. Second, clear-sky LST was hypothetically corrected by accounting for cloud cooling based on SEB theory. The proposed method was applied to Visible Infrared Imaging Radiometer Suite (VIIRS) and Moderate Resolution Imaging Spectroradiometer (MODIS) data, and further validated using ground measurements of fourteen sites from SURFRAD, BSRN, and AmeriFlux in 2013. VIIRS LST recovered from cloud gaps exhibited a root mean square error (RMSE) of 3.54 K, a bias of −0.36 K, R2 of 0.94, and sample size (N) of 2411, comparable to the accuracy of clear-sky LST products and cloudy-sky LST estimation from MODIS (RMSE of 3.69 K, bias of −0.45 K, R2 of 0.93, and N of 2398). Thus, the proposed method performs well across different sensors, seasons, and land cover types. The abnormal retrieval values caused by cloud contamination were also corrected in the proposed method. The overall accuracy was better than the downscaled cloudy-sky LST retrieved from passive microwave (PMW) observations and former SEB-based cloudy-sky LST estimation methods. Validation using time-series measurements showed that the all-sky LST time series, including both clear- and cloudy-sky retrievals, can capture realistic variability without sudden abruptions or discontinuities. RMSE values for the all-sky LST varied from 2.54 to 4.15 K at the fourteen sites. Spatially continuous LST maps over the Contiguous United States were compared with corresponding maps from PMW data in the winter and summer of 2018, exhibiting similar spatial patterns but with additional spatial details. Moreover, sensitivity analysis suggested that the reconstruction of clear-sky LST dominantly impacts the accuracy of cloudy-sky LST estimation. The proposed method can be potentially implemented in similar satellite sensors for global real-time production.