In this study, two collection 6 (C6) Moderate Resolution Imaging Spectroradiometer (MODIS) level-2 land surface temperature (LST) products (MYD11_L2 and MYD21_L2) from the Aqua satellite were evaluated using temperature-based (T-based) and radiance-based (R-based) validation methods over barren surfaces in Northwestern China. The ground measurements collected at four barren surface sites from June 2012 to September 2018 during the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment were used to perform the T-based evaluation. Ten sand dune sites were selected in six large deserts in Northwestern China to carry out an R-based validation from 2012 to 2018. The T-based validation results indicate that the C6 MYD21 LST product has a better accuracy than the C6 MYD11 product during both daytime and nighttime. The LST is underestimated by the C6 MYD11 products at the four T-based sites during the daytime, with a mean bias of -2.82 K and a mean RMSE of 3.82 K, whereas the MYD21 LST product has a mean bias and RMSE of -0.51 and 2.53 K, respectively. The LST is also underestimated at night by the C6 MYD11 products at the four T-based sites, with a mean bias of -1.40 K and a mean RMSE of 1.72 K, whereas the MYD21 LST product has a mean bias and RMSE of 0.23 and 1.01 K, respectively. For the R-based validation, the MYD11 results are associated with large negative biases during both daytime and nighttime at three sand dune sites and biases within 1 K at the other seven sites, whereas the MYD21 results are more consistent at all ten sand dune sites, with a mean bias of 0.45 and 0.70 K for daytime and nighttime, respectively. The emissivities for these two products in MODIS bands 31 and 32 were compared with each other and then compared with the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) emissivity and laboratory emissivity. The results indicate that the emissivities in MODIS bands 31 and 32 of MYD11 at the four T-based and three of the R-based validation sites are overestimated and result in LST underestimation, whereas the emissivities of MYD21 are more consistent with the laboratory emissivity. Besides, an experiment was carried out to demonstrate that the physically retrieved dynamic emissivity of the MYD21 product can be utilized to improve the accuracy of the split-window (SW) algorithm for barren surfaces, making it a valuable data source for retrieving LST from different remote sensing data.

A new directional canopy emissivity model (CE-P) based on spectral invariants is proposed in this paper. First, we prove the existence of the spectral invariant properties in the thermal infrared (TIR) band using a Monte Carlo model. Based on it, the equation of the new model is derived from the perspective of absorption. In this expression, single-scattering and multiscattering effects are separated analytically in the TIR band. We find that the overall contribution of multiple scatterings is less than 0.005 when the component emissivities are over 0.90, and the overall contribution decreases with increasing leaf or soil emissivity. Furthermore, the new model can avoid the logical difficulty encountered when using the traditional cavity effect factor to simulate the emissivity of a sparse vegetation canopy. The results of 4SAIL and Discrete Anisotropic Radiative Transfer (DART) are selected to do cross validation. The CE-P can achieve a high accuracy compared with 4SAIL and DART, with an absolute bias less than 0.002 when the leaf (soil) emissivity is equal to 0.98 (0.94). Four widely used analytical models are selected for comparison. The resulting accuracies of these models are ordered from CE-P to REN15, FR97, FR02, and VALOR96 with the most serious error up to 0.002, 0.002, 0.007, 0.013, and 0.014, respectively. Three main conclusions are obtained through the sensitivity analysis: the multiscattering between vegetation and the background can be ignored when the leaf (soil) emissivity is no less than 0.94 (0.90), the second and higher order scattering within the vegetation can also be ignored when the leaf (soil) emissivity is no less than 0.94 (0.90), and the single-scattering effect within the canopy should be considered which can be calculated using three view factors.

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.