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.
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.
18 September 2023 in IEEE Transactions on Geoscience and Remote Sensing
Land surface temperature (LST) is an important parameter that characterizes the energy balance of the land surface, and it is widely used in various research fields. This paper proposes an operational split-window (SW) algorithm for use with the Chinese Fengyun-3 (FY-3) series satellite data, with the purpose of generating long-term global LST products. The algorithm primarily involves three steps. First, the brightness temperatures of the FY-3 Visible and Infra-Red Radiometer (VIRR) were recalibrated using historical recalibration coefficients to improve the accuracy of the absolute radiometric calibration. Second, daily dynamic emissivity maps were estimated using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) global emissivity dataset (GED) and vegetation/snow cover products based on the vegetation cover method. Finally, the coefficients of the SW algorithm were simulated using MODTRAN 5 combined with the SeeBor V5.0 atmospheric profile library and ASTER spectral library, and then the coefficients were stratified by the view zenith angle and atmospheric water vapor content to improve the fitting accuracy. The proposed SW algorithm was integrated into the MUlti-source data SYnergized Quantitative (MUSYQ) remote sensing production system to then generate FY-3 VIRR LST products. Ten land surface sites from the HiWATER and SURFRAD networks and nine water surface sites from the National Data Buoy Center (NDBC) were used to evaluate the accuracy of the FY-3 VIRR LST products. The results demonstrated that the accuracy of the historical recalibration coefficients of the FY-3A/B VIRR is higher than that of the operational calibration coefficients for LST retrieval. The evaluation results revealed that the FY-3A VIRR LST products (2009-2013) had a bias of 0.13 K and an RMSE of 2.77 K, and the FY-3B VIRR LST products (2011-2020) had a bias of -0.07 K and an RMSE of 2.83 K. These results demonstrate that the proposed operational SW algorithm has reasonable accuracy and can be used to produce global LST products from the FY-3 VIRR data.
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.
1 July 2021 in IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing
Land surface temperature (LST) is a key parameter in numerous environmental studies. However, currently, there is no satellite sensor that can completely provide LST data with both high spatial and high temporal resolutions simultaneously. LST downscaling is regarded as an effective remedy for improving the temporal and spatial resolutions of LST data. In this article, a geographically and temporally weighted autoregressive (GTWAR) model of LST downscaling is that comprehensively considers the spatial heterogeneity, spatial autoregression, and temporality of LST is newly proposed. The normalized difference water index, the normalized difference built-up index, and the normalized difference vegetation index were selected as explanatory variables to downscale the moderate resolution imaging spectroradiometer (MODIS) LST from 1000 to 100 m, while the Landsat 8 LST was selected as the reference data. Compared with the thermal data sharpening (TsHARP), the geographically weighted regression (GWR), the geographically weighted autoregressive (GWAR) and the geographically and temporally weighted regression (GTWR) downscaling methods, the proposed method was superior based on quantitative indices, with the lowest root mean square error (Zhangye: 1.57 °C, Beijing: 1.22 °C) and mean absolute error (Zhangye: 1.06 °C, Beijing: 0.85 °C). The downscaling model of GTWAR will facilitate improvements in the accuracy of downscaling for temporal series of LST data.
28 December 2022 in IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing
The Thermal Infrared Sensor-2 (TIRS-2) carried on Landsat 9 is the newest thermal infrared (TIR) sensor for the Landsat project and provides two adjacent TIR bands, which greatly benefits the land surface temperature (LST) retrieval at high spatial resolution. In this study, a radiance based split window (RBSW) algorithm for retrieving LST from Landsat 9 TIRS-2 data was proposed. In addition, the split-window covariance-variance ratio (SWCVR) algorithm was improved and applied to Landsat 9 TIRS-2 data for estimating atmospheric water vapor (AWV) that is required for accurate LST retrieval. The performance of the proposed method was assessed using the simulation data and satellite observations. Results reveal that the retrieved LST using the RBSW algorithm has a bias of 0.06 K and root-mean-square error (RMSE) of 0.51 K based on validation with the simulation data. The sensitivity analysis exhibited a LST error of <1.75 K using the RBSW algorithm when the uncertainties in input parameters (i.e., AWV, emissivity and at-sensor radiance) were considered. There is a marginal discrepancy between LST retrievals using the estimated AWV and MODIS AWV, with a difference in bias of -0.14 K and RMSE of 0.22 K, which indicates that the improved SWCVR method can provide an optional means to obtain AWV inputted in LST retrieval. With regard to the validation using the in-situ measurements, the retrieved LST from Landsat 9 TIRS-2 data exhibits a bias of 0.44 K and RMSE of 1.98 K, respectively, showing a higher accuracy than the USGS Landsat LST product with bias of 1.21 and RMSE of 2.56 K. In conclusion, the proposed algorithm combining RBSW algorithm and improved SWCVR algorithm for LST retrieval from Landsat 9 has a good accuracy without dependence on external atmospheric data, and it is expected to be a reliable method for LST generation from Landsat 9 TIRS-2 data.