In this study, to improve the accuracy of land surface temperature (LST) products over barren surfaces, we present an operational algorithm to retrieve the LST from Moderate-Resolution Imaging Spectroradiometer (MODIS) thermal infrared data using physically retrieved emissivity products. The LST algorithm involved two steps. First, the emissivity in the two MODIS split-window (SW) channels was estimated using the vegetation cover method, with the bare soil component emissivity derived from the ASTER global emissivity data set. Then, the LST was retrieved using a modified generalized SW algorithm. This algorithm was implemented in the MUlti-source data SYnergized Quantitative (MuSyQ) remote sensing product system. The MuSyQ MODIS LST product and the Collection 6 MODIS LST product (MxD11_L2) were compared and validated using ground measurements collected from four barren surface sites in Northwest China during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment from June 2012 to December 2015. In total, 2268 and 2715 clear-sky samples were used in the validation for Terra and Aqua, respectively. The evaluation results indicate that the MuSyQ LST products provide better accuracy than the C6 MxD11 product during both daytime and nighttime at all four sites. For the daytime results, the LST is underestimated by the C6 MxD11 products at all four sites, with a mean bias of -1.78 and -2.86 K and a mean root-mean-square error (RMSE) of 3.16 and 3.94 K for Terra and Aqua, respectively, whereas the mean biases of the MuSyQ LST products are within 1 K, with a mean bias of -0.26 and -1.03 K and a mean RMSE of 2.45 and 2.71 K for Terra and Aqua, respectively. For the nighttime results, the LST is also underestimated by the C6 MxD11 products at all four sites, with a mean bias of -1.60 and -1.26 K and a mean RMSE of 1.93 and 1.60 K for Terra and Aqua, respectively, whereas the mean biases of the MuSyQ LST products are 0.16 and 0.58 K and the mean RMSEs are 1.12 and 1.25 K for Terra and Aqua, respectively. The results indicate that the underestimation of the C6 MxD11 LST product at all four sites mainly results from the overestimation of the emissivities in MODIS bands 31 and 32. This study demonstrates that physically retrieved emissivity products are a useful source for LST retrieval over barren surfaces and can be used to improve the accuracy of global LST products.

Land surface temperature (LST) is a vital parameter for studying global ecological, climatic, and environmental changes. Although various LST retrieval algorithms have been proposed, including split-window (SW), dual-window (DW), three-channel (TC), and dual-angle (DA) algorithms, few studies have compared these algorithms using the same satellite observations. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard Sentinel-3A provides a unique opportunity to conduct this comparison owing to its dual-angle viewing capability and multiple thermal infrared (TIR) and mid-infrared (MIR) channels. Here, we implemented two SW algorithms, one DW algorithm, two TC algorithms and one DA algorithm for the SLSTR data. The LST retrievals from these six algorithms were validated, along with the SLSTR operational LST product based on an emissivity-implicit SW algorithm. Temperature-based and radiance-based validation methods were used to evaluate different LST retrievals across different land cover types. The results indicated that the proposed SW algorithm had the highest accuracy, followed by the Pérez-Planells SW and the official algorithms. The overall root-mean-square errors (RMSEs) of these three SW algorithms were 1.42 K, 1.79 K and 2.05 K, respectively. The three algorithms involving the MIR channel (one DW and two TC algorithms) were more suitable for nighttime LST retrieval and had similar performances to the three SW algorithms, with a nighttime RMSE of approximately 1.36 K. The LST retrieval accuracy of the DA algorithm had the highest uncertainty and was closely related to the angular variation in surface emissivity and brightness temperature. The findings of this study contribute to a better understanding of the different LST retrieval algorithms and facilitate potential improvements in the official LST retrieval algorithm for SLSTR.

Many physical models have been proposed to simulate the directional anisotropy in the thermal infrared (TIR) region over vegetation canopies to produce angular corrected directional brightness temperature or land surface temperature. However, too many input parameters obstruct their operational use. Semiempirical kernel-driven models are designed to be a tradeoff between physical accuracy and operationality. Recently, four kernel-driven models have been proposed: the first two are direct extensions of kernel models in the visible- and near-infrared region and the last two were directly designed for the TIR region. In this paper, 153 continuous and 153 discrete canopies with varying structures and temperature distributions were considered in order to evaluate their accuracies against two physical models (4SAIL and DART). Their error distribution, scatterplots, and directional anisotropy patterns are compared. LSF-Li model, followed by Ross-Li, Vinnikov, and RL model, gave the best fitting results for all the scenes. The R² of all four kernel models can reach up to 0.82 for discrete scenes; however, the kernel-driven models underestimate the hotspot effect from continuous scenes; therefore, further improvements are necessary for operational use with future TIR satellite missions.

This paper uses the refined Generalized Split-Window (GSW) algorithm to derive the land surface temperature (LST) from the data acquired by the Visible and Infrared Radiometer on FengYun 3B (FY-3B/VIRR). The coefficients in the GSW algorithm corresponding to a series of overlapping ranges for the mean emissivity, the atmospheric Water Vapor Content (WVC), and the LST are derived using a statistical regression method from the numerical values simulated with an accurate atmospheric radiative transfer model MODTRAN 4 over a wide range of atmospheric and surface conditions. The GSW algorithm is applied to retrieve LST from FY-3B/VIRR data in an arid area in northwestern China. Three emissivity databases are used to evaluate the accuracy of different emissivity databases for LST retrieval, including the ASTER Global Emissivity Database (ASTER_GED) at a 1-km spatial resolution (AG1km), an average of twelve ASTER emissivity data in the 2012 summer and emissivity spectra extracted from spectral libraries. The LSTs retrieved from the three emissivity databases are evaluated with ground-measured LST at four barren surface sites from June 2012 to December 2013 collected during the HiWATER field campaign. The results indicate that using emissivity extracted from ASTER_GED can achieve the highest accuracy with an average bias of 1.26 and −0.04 K and an average root mean square error (RMSE) of 2.69 and 1.38 K for the four sites during daytime and nighttime, respectively. This result indicates that ASTER_GED is a useful emissivity database for generating global LST products from different thermal infrared data and that using FY-3B/VIRR data can produce reliable LST products for other research areas.

An operational split-window (SW) algorithm was developed to retrieve high-temporal-resolution land surface temperature (LST) from global geostationary (GEO) satellite data. First, the MODTRAN 5.2 and SeeBor V5.0 atmospheric profiles were used to establish a simulation database to derive the SW algorithm coefficients for GEO satellites. Then, the dynamic land surface emissivities (LSEs) in the two SW bands were estimated using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Emissivity Dataset (GED), fractional vegetation cover (FVC), and snow cover products. Here, the proposed SW algorithm was applied to Himawari-8 Advanced Himawari Imager (AHI) observations. LST estimates were retrieved in January, April, July, and October 2016, and three validation methods were used to evaluate the LST retrievals, including the temperature-based (T-based) method, radiance-based (R-based) method, and intercomparison method. The in situ night-time observations from two Heihe Watershed Allied Telemetry Experimental Research (HiWATER) sites and four Terrestrial Ecosystem Research Network (TERN) OzFlux sites were used in the T-based validation, where a mean bias of −0.70 K and a mean root-mean-square error (RMSE) of 2.29 K were achieved. In the R-based validation, the biases were 0.14 and −0.13 K and RMSEs were 0.83 and 0.86 K for the daytime and nighttime, respectively, over four forest sites, four desert sites, and two inland water sites. Additionally, the AHI LST estimates were compared with the Collection 6 MYD11_L2 and MYD21_L2 LST products over southeastern China and the Australian continent, and the results indicated that the AHI LST was more consistent with the MYD21 LST and was generally higher than the MYD11 LST. The pronounced discrepancy between the AHI and MYD11 LST could be mainly caused by the differences in the emissivities used. We conclude that the developed SW algorithm is of high accuracy and shows promise in producing LST data with global coverage using observations from a constellation of GEO satellites.

Surface upward longwave radiation (SULR) is a critical component in the calculation of the Earth’s surface radiation budget. Multiple clear-sky SULR estimation methods have been developed for high-spatial resolution satellite observations. Here, we comprehensively evaluated six SULR estimation methods, including the temperature-emissivity physical methods with the input of the MYD11/MYD21 (TE-MYD11/TE-MYD21), the hybrid methods with top-of-atmosphere (TOA) linear/nonlinear/artificial neural network regressions (TOA-LIN/TOA-NLIN/TOA-ANN), and the hybrid method with bottom-of-atmosphere (BOA) linear regression (BOA-LIN). The recently released MYD21 product and the BOA-LIN—a newly developed method that considers the spatial heterogeneity of the atmosphere—is used initially to estimate SULR. In addition, the four hybrid methods were compared with simulated datasets. All the six methods were evaluated using the Moderate Resolution Imaging Spectroradiometer (MODIS) products and the Surface Radiation Budget Network (SURFRAD) in situ measurements. Simulation analysis shows that the BOA-LIN is the best one among four hybrid methods with accurate atmospheric profiles as input. Comparison of all the six methods shows that the TE-MYD21 performed the best, with a root mean square error (RMSE) and mean bias error (MBE) of 14.0 and −0.2 W/m², respectively. The RMSE of BOA-LIN, TOA-NLIN, TOA-LIN, TOA-ANN, and TE-MYD11 are equal to 15.2, 16.1, 17.2, 21.2, and 18.5 W/m², respectively. TE-MYD21 has a much better accuracy than the TE-MYD11 over barren surfaces (NDVI < 0.3) and a similar accuracy over non-barren surfaces (NDVI ≥ 0.3). BOA-LIN is more stable over varying water vapor conditions, compared to other hybrid methods. We conclude that this study provides a valuable reference for choosing the suitable estimation method in the SULR product generation.