Next Article in Journal
Robust Underwater Direction-of-Arrival Estimation Method Using Acoustic Sensor Array under Unknown Swing Deviation Elements
Next Article in Special Issue
Melt Pond Evolution along the MOSAiC Drift: Insights from Remote Sensing and Modeling
Previous Article in Journal
A New Robust Lunar Landing Selection Method Using the Bayesian Optimization of Extreme Gradient Boosting Model (BO-XGBoost)
Previous Article in Special Issue
Stripe Noise Elimination with a Novel Trend Repair Method for Push-Broom Thermal Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation and Validation of Split-Window Algorithms for Estimating Land Surface Temperature from Landsat 9 TIRS-2 Data

by
Qinghua Su
1,
Xiangchen Meng
1,* and
Lin Sun
2
1
School of Geography and Tourism, Qufu Normal University, Rizhao 276826, China
2
College of Geodesy and Geomatics, Shandong University of Science and Technology, Qingdao 266590, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(19), 3633; https://doi.org/10.3390/rs16193633
Submission received: 28 August 2024 / Revised: 23 September 2024 / Accepted: 25 September 2024 / Published: 29 September 2024
(This article belongs to the Special Issue Surface Radiative Transfer: Modeling, Inversion, and Applications)

Abstract

:
Land surface temperature (LST) is important in a variety of applications, such as urban thermal environment monitoring and water resource management. In this paper, eleven candidate split-window (SW) algorithms were adapted to Thermal Infrared Sensor-2 (TIRS-2) data of the Landsat 9 satellite for estimating the LST. The simulated dataset produced by extensive radiative transfer modeling and five global atmospheric profile databases was used to determine the SW algorithm coefficients. Ground measurements gathered at Surface Radiation Budget Network sites were used to confirm the efficiency of the SW algorithms after their performance was initially examined using the independent simulation dataset. Five atmospheric profile databases perform similarly in training accuracy under various subranges of total water vapor. The candidate SW algorithms demonstrate superior performance compared to the radiative transfer equation algorithm, exhibiting a reduction in overall bias and RMSE by 1.30 K and 1.0 K, respectively. It is expected to provide guidance for the generation of the Landsat 9 LST using the SW algorithms.

Graphical Abstract

1. Introduction

Land surface temperature and emissivity (LST&E) are the key parameters to control the water and energy budget balance and have been identified as important Earth System Data Records (ESDRs) by the National Aeronautics and Space Administration (NASA) and many other international organizations [1,2,3,4,5]. The LST product generated from thermal infrared (TIR) remote sensing provides a unique opportunity to obtain the LST at regional and global scales and has been widely used in hydrology, meteorology, and surface energy balance [6,7,8].
With the launch of Landsat 9 on 27 September 2021, the Landsat series of satellites have been providing TIR observations for 40 years, thus constituting a highly valuable data resource for the generation of LST records [9,10,11]. The integration of the Landsat 9 and Landsat 8 satellites will facilitate cooperative observations with an eight-day revisit period. Landsat 9 is equipped with the TIR Sensor-2 (TIRS-2), which features two TIR channels similar to the TIRS on Landsat 8. However, the TIRS-2 has significantly reduced stray light, thereby enabling more accurate LST [12,13]. Researchers have studied the accuracy of Landsat 9 LST products produced by the radiative transfer equation (RTE) algorithm since the United States Geological Survey (USGS) made these products available. In comparison to ground measurements, the average bias and root mean square error (RMSE) were found to be 0.24 K and 3.42 K, respectively. Additionally, the Landsat 9 LST product showed good agreement with the Landsat 7/8 LSTs [10]. The ground campaign conducted by the Rochester Institute of Technology (RIT) also indicated that Landsat 9 is exhibiting consistent behavior with Landsat 8, the ground reference, and unmanned aircraft system (UAS) measurements [9]. It should be noted that the LST differences between Landsat 8 and Landsat 9 are typically more pronounced in barren desert regions and areas with a high concentration of built-up land [14].
Since 1 March 2023, the Landsat 9 Collection 2 data from the first year of operation have been reprocessed to incorporate calibration updates identified by the USGS/NASA calibration and validation team. Using accurate ground measurements as a reference, Niclòs et al. performed vicarious calibrations on the originally calibrated and reprocessed Landsat 9 TIRS-2 data and assessed the impact of the reprocessing on the LST retrieval [15]. No significant calibration errors were observed, and both single-channel and split-window (SW) algorithms can provide LST with uncertainties around 1 K. Theoretically, the split-window algorithm performs better than the single-channel algorithm over global conditions without needing high-precision atmospheric profiles [16,17,18]. The researchers also developed SW and ensemble learning algorithms for Landsat 9 LST retrieval [11,19,20]. For example, a radiance-based split-window (RBSW) algorithm was proposed to retrieve Landsat 9 LST and the validation result indicated that a reduction of about 0.8 K in RMSE was achieved compared to the official LST product [11]. Taking full advantage of the ensemble learning method, a new algorithm is proposed to retrieve LST directly from the Landsat 9 top-of-atmosphere brightness temperature without requiring external parameters [20]. Compared to ground measurements, the algorithm achieves a bias and RMSE of −0.465 K and 2.508 K, respectively.
Over the past 30 years, various SW algorithms have been developed to perform LST retrieval using different TIR sensors. Among the published SW algorithms, which one is the most reliable and accurate is unclear because every sensor has a different filter response and central wavelength. Extensive investigations are needed to validate the accuracy of Landsat 9 LSTs derived from various SW algorithms. In this paper, we adapted eleven published SW algorithms to TIRS-2 data and evaluated their performance using ground measurements. This paper has the following structure. In Section 2, we describe the SW algorithms, the land surface emissivity estimation, the simulation dataset, and the ground LST estimation. The results and analysis are presented in Section 3. The discussion and conclusions are given in Section 4 and Section 5.

2. Datasets and Methods

2.1. Split-Window LST Algorithms

After conducting a comprehensive review of LST algorithms, we selected eleven SW algorithms from the literature as candidate algorithms for the TIRS-2, which were also summarized in previous research [21,22,23]. Table 1 presents the formulas for the eleven candidate LST algorithms, which utilize the land surface emissivity and the brightness temperatures of two TIRS-2 bands centered at 11 and 12 μm, respectively. Although the total water vapor (TWV) can be explicitly considered in the SW algorithm, it is not considered here. On the one hand, the uncertainty in the TWV will affect the accuracy of the retrieved LST. On the other hand, previous research suggested that this led to only minor enhancements in the accuracy of SW algorithms [15,24]. Additionally, the TIRS-2 has a near-zero view zenith angle, so the path length correction term is not necessary for the SW algorithms [25].
In this study, we utilized both Landsat 9 Collection 2 (C2) Level-1 (L1) and Level-2 (L2) products obtained from https://earthexplorer.usgs.gov (accessed on 1 June 2024) for LST retrieval. The brightness temperatures of channels 10 and 11 were calculated from digital numbers provided in the Landsat 9 C2 L1 product. The Normalized Difference Vegetation Index (NDVI) was calculated using the land surface reflectance (LSR) product included in the Landsat 9 C2 L2 product. To eliminate the impact of clouds and cloud shadows, the quality assessment (QA) band was employed.

2.2. Land Surface Emissivity Estimation

Taking into account the emissivity increment from the cavity effect caused by the multiple scattering in the pixel, the Landsat 9 LSE can be calculated from the improved NDVI threshold method, which can be expressed by the following equation [37,38,39]:
ε i = a 1 i + j = 2 7 a j i ρ j N D V I < N D V I s ε v i P v + ε s i 1 P v + 4 d ε P v 1 P v N D V I s N D V I  
P v = N D V I N D V I s N D V I v N D V I s 2
d ε = ε v i ( 0.435 ε s i + 0.4343 ) / 0.985
where ε i is the Landsat LSE of channel 10 or 11,   ε v i and ε s i are the vegetation and soil component emissivities, and ρ j is the LSR of the Operational Land Imager 2(OLI-2) channel j. According to the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) spectral library,   ε v i and ε s i are set to 0.9847 and 0.9706, 0.9854 and 0.9769 for channels 10 and 11, respectively. The coefficients a j i used to calculate the soil emissivity for Landsat 9 channels 10 and 11 are shown in Table 2, with corresponding root mean square errors (RMSEs) of 0.0043 and 0.0029 for channels 10 and 11, respectively.

2.3. Simulation Dataset

We utilized the moderate spectral resolution atmospheric transmittance model version 5.2 (MODTRAN5.2), global atmospheric profiles, and the ASTER/MODIS spectral library to build a comprehensive simulation database for deriving coefficients essential for the mentioned SW algorithms. The specific process was as follows:
1:
Global atmospheric profiles and the spectral response functions of TIRS-2 channels were input into the MODTRAN5.2 to obtain three atmospheric parameters (atmospheric upward radiance L i , downward radiance L i , and transmittance τ i ) for TIRS-2 channels [40,41]. In previous research, the global atmospheric profile databases such as Thermodynamic Initial Guess Retrieval (TIGR) [42], Global Atmospheric Profiles from Reanalysis Information (GAPRI) [43], SeeBor V5.0 training database of global profiles (SeeBor for simplicity as follows) [44], and Cloudless Land Atmosphere Radiosounding (CLAR) [45] were frequently utilized to compute the SW algorithm coefficients for thermal sensors onboard different satellites [46,47,48,49,50,51,52]. Recently, the comprehensive clear-sky database based on the European Centre for Medium-Range Forecast (ECMWF) version-5 reanalysis (ERA5 for simplicity as follows) [53] was also used to calculate atmospheric parameters [54]. This research utilized the five atmospheric profile databases mentioned above to calculate the SW algorithm coefficients. Additionally, the radiosounding profiles collected from global radiosonde stations [40,55] were used independently for validation purposes.
2:
Emissivity spectra selected from the ASTER/MODIS spectral library were convolved with the spectral response function of TIRS-2 channels to obtain the channel-integrated LSE ( ε i ). In total, 110 emissivity spectra were selected for simulation, including 19 vegetation types, 39 soil types, 5 snow/water types, 35 rock types, and 12 manmade material types [24].
ε i = λ 1 λ 2 ε ( λ ) f i ( λ ) d λ λ 1 λ 2 f i ( λ ) d λ
where ε i is the channel-integrated LSE of channel i, ε ( λ ) is the measured emissivity spectra of various types, and f i λ is the spectral response function of TIRS-2 channel i.
3:
The top-of-atmosphere radiation observed by TIRS-2 channels ( L i T O A ) can be calculated from the following radiative transfer equation:
L i T O A = τ i ε i B i ( T s ) + ( 1 ε i ) L i + L i
where L i , L i , a n d   τ i are the atmospheric upward radiance, downward radiance, and transmittance of channel i, respectively. B i is the channel-averaged Planck’s radiance function of channel i.  T s is the LST used for the simulation. We followed the method described in [16] and set the T s as T 0 − 5, T 0 , T 0 + 5, T 0 + 10, T 0 + 15, and T 0 + 20, with T 0 representing the air temperature at the first layer of the atmospheric profiles [56].
4:
The brightness temperatures of TIRS-2 channels ( T i ) can be calculated using the inversion of the Plank function [57]:
T i = h c k λ i 1 ln ( 2 h c 2 λ i 5 / L i T O A + 1 )
where h , c ,   a n d   k are Plank’s constant, the light speed, and Boltzmann’s constant, respectively. λ i is the central wavelength of channel i. For the TIRS-2 channels, the values are 10.8 and 12.0 μ m , respectively [58].
5:
The SW algorithm coefficients in Table 1 are determined by ordinary least squares linear regression. The candidate SW algorithm coefficients obtained from the SeeBor databases are given in Appendix A. Distribution of the air temperature at the first layer and the total water vapor of the six atmospheric profile databases are given in Appendix B.

2.4. Ground Measurements

The estimated Landsat LST was validated using ground measurements from the Surface Radiation Budget Network (SURFRAD), which is extensively used in evaluating LST and surface longwave radiation products [59,60,61,62]. Specific information on each in situ site is shown in Table 3. With an uncertainty of about 5 W/m2 and an effective footprint of about 70 × 70 m2, the pyrgeometer measures upward and downward longwave radiances at 3-minute intervals before January 2009 and every minute after that [63,64]. The in situ LSTs were calculated from the site-measured upward and downward longwave radiation using the following equation [65]:
T s = F ( 1 ε b ) F ε b σ 1 4
where T s is the surface LST, F is the measured upward longwave radiation, F is the measured downward longwave radiation, σ is Stefan–Boltzmann’s constant (5.67 × 10−8 W/m2/K4), and ε b is the BBE, which can be calculated from the in situ measured emissivity or the five ASTER spectral emissivities using the following equation [66]:
ε = 0.197 + 0.025 ε 10 + 0.057 ε 11 + 0.237 ε 12 + 0.333 ε 13 + 0.146 ε 14

3. Results and Analyses

3.1. SW Algorithm Training

Figure 1 illustrates the performance of the eleven candidate SW algorithms under different TWV subranges. As illustrated in Figure 1, when the TWV is below 1.5 g/cm2, the SW algorithms show similar performance, with RMSEs ranging between 0.24 and 0.44 K, except for the SW5 and SW8 algorithms. The SW algorithms trained using five atmospheric profile databases also show consistent performance in uncertainty. When the TWV is between 1.5 and 3.0 g/cm2, the SW algorithms trained using the CLAR database exhibit the best results, with RMSEs of about 0.80 K for the SW5 and SW8 algorithms and less than 0.50 K for the other algorithms. The GAPRI, SeeBor, and TIGR databases demonstrate similar RMSEs in various SW algorithms, while the ERA5 database shows slightly inferior results with RMSEs greater than 0.52 K for all the algorithms. When the TWV is between 3.0 and 4.5 g/cm2, the SW algorithms trained using the CLAR database exhibit the best results, followed by the SeeBor, GAPRI, and ERA5 databases. The TIGR database produces the least favorable results, with RMSEs ranging from 0.70 to 0.94 K. When the TWV exceeds 4.5 g/cm2, the GAPRI database comes from behind and performs best, with RMSEs of the SW algorithms between 0.65 and 0.84 K. The RMSE values for the eleven SW algorithms trained using the TIGR database are all above 1.26 K. With the increase in the TWV, SW5 and SW8 demonstrate comparable performance to other SW algorithms, with an RMSE difference of less than 0.1 K. Without classification of the TWV, the Seebor database has the highest simulation accuracy, followed by the GAPRI and TIGR databases. The SW2 and SW11 algorithms, which utilize coefficients simulated from the SeeBor database, demonstrate superior performance compared to other SW algorithms that employ coefficients simulated from alternative databases.
The above-mentioned results show that most SW algorithms demonstrate slight deviations in training accuracy when operating under specific TWV subranges. The SW1, SW2, SW3, and SW4 algorithms have been found to demonstrate slightly superior training accuracy across all four TWV subranges. In contrast, the SW5 and SW8 algorithms have been identified as exhibiting less optimal training accuracy. In the event that the TWV is less than 4.5 g/cm2, any of the five atmospheric profile databases may be employed for simulation purposes, as they yield comparable results. Conversely, the GAPRI database is recommended for simulation when the TWV is above 4.5 g/cm2.

3.2. Sensitivity Analysis

The sensor noise, LSE, and TWV represent three principal sources of uncertainty that contribute to errors in LST retrieval [7,24]. To address this, we conducted an assessment of the sensitivity of the candidate SW algorithms to these three variables. The eleven candidate SW algorithm coefficients, calculated from the CLAR database, were tested with an independent dataset simulated using radiosounding profiles. Normal distributed random noise with a mean of 0 and a standard deviation of 0.07, 0.15, and 0.40 K was incorporated into the simulated TIRS-2 brightness temperature [58]. Furthermore, uncertainties of 0.005, 0.01, and 0.015 (unitless) and 0.5, 0.8, and 1.0 g/cm2 were added to the LSE and TWV datasets, respectively [19]. Figure 2 shows the performance of the eleven candidate SW algorithms before and after the addition of random noise under different TWV subranges.
As illustrated in Figure 2a,d, the SW5 and SW8 algorithms exhibit the least sensitivity to emissivity uncertainty, while the SW10 and SW11 algorithms demonstrate the second-least sensitivity. In contrast, the SW1 and SW2 algorithms exhibit the greatest sensitivity to emissivity uncertainty. The LST uncertainty increases as the emissivity noise increases. This increase is more pronounced for lower TWV than for higher TWV. For example, when the emissivity uncertainty is 0.01, the RMSE of the SW2 algorithm increases from 0.27 K to 2.22 K and from 1.04 K to 1.41 K when the TWV is less than 1.5 g/cm2 and greater than 4.5 g/cm2, respectively. This suggests that the uncertainty in the emissivity has a more significant impact on the retrieval of the LST under dry atmospheric conditions, such as in northwest China and western Australia.
As shown in Figure 2b,e, when the uncertainty in the brightness temperature is equal to 0.15 K, the increase ratio of the RMSE is less than 0.5% for the majority of algorithms, with values above 1.0% for the SW1 and SW2 algorithms under a lower TWV subrange. Furthermore, as the TWV increases, the ratio of the increase declines. Based on the performance of on-orbit calibrations, the noise level in both channels of TIRS-2 is below 0.1 K [15,67]. Therefore, the uncertainty in the brightness temperature has a slight impact on the retrieval of the LST, both in dry and humid atmospheric conditions.
As indicated in Figure 2c,f, when the uncertainty in the TWV is less than 1.0 g/cm2, the increase ratio of the RMSE is less than 0.6% for the eleven candidate algorithms. In other words, the uncertainty in the TWV has a negligible effect on the retrieval of the LST. But as indicated in previous research [24,68], the incorrect selection of the TWV subrange led to a significant increase in LST retrieval. For example, when the TWV were below 1.5 g/cm2, the bias and RMSE for the SW2 algorithm were 0.0 K and 0.41 K, respectively. However, when a subrange of coefficients spanning 1.5–3.0 g/cm2 was employed, the values shifted to −0.66 K and 0.79 K, respectively.

3.3. T-Based Validation Results

Figure 3 presents the bias and RMSE between the retrieved Landsat 9 LST and the in situ LST for all the matchups at the SURFRAD sites. After removing invalid matchups affected by clouds and cloud shadows based on the Landsat level-2 pixel quality channel, 234 matchups were obtained at the SURFRAD sites. In addition, the official Landsat 9 LST that was retrieved by the RTE algorithm was also validated. As illustrated in Figure 3, the Landsat 9 LSTs retrieved using the SW algorithm achieved a lower RMSE than those retrieved using the RTE algorithm when compared to the ground measurements, except for the PSU site. The GWN and SXF sites had a root mean square error (RMSE) of approximately 1.3 and 1.6 K, followed by the BND, PSU, and DRA sites with an RMSE ranging between 2.0 and 2.7 K. The FPK and TBL sites achieved an RMSE above 5.0 K for most of the SW algorithms. The RTE algorithm performed slightly better than the SW algorithm regarding the RMSE at the PSU station and provided consistently worse results at the other stations. In terms of the average difference (bias) between the retrieved Landsat 9 LST and the in situ LST, the Landsat 9 LST retrieved from the SW algorithm were underestimated at the PSU, GWN, and DRA sites. Conversely, the corresponding LSTs were overestimated at the BND, SXF, TBL, and FPK sites. However, the RTE algorithm underestimated the LST at the PSU site and overestimated it elsewhere.
The validation results above indicate that the candidate SW algorithm is suitable for estimating the LST from the Landsat 9 TIR data. There was no significant difference between the eleven candidate SW algorithms, with overall biases and RMSEs ranging from 1.44 to 1.61 K and 3.50 to 3.73 K, whereas these values were 2.82 K and 4.69 K for the RTE algorithm.

4. Discussion

Since March 2023, Landsat 9 Collection 2 data have been reprocessed based on the initial calibration results, which may affect the accuracy of the retrieved LST. Based on the research of Niclòs et al. [15], the official LST product retrieved with the original and reprocessed TIRS-2 data have similar accuracy in Valencia, but a systematic difference of about −0.6 K was observed for the SW LST. We also calculated the statistics between the official Landsat 9 LST product and the in situ LST obtained at the SURFRAD sites. The bias (RMSE) for the original and reprocessed TIRS-2 data were 1.44 (2.94) K and 1.42 (2.93) K. The validation result in Section 3 indicates that the SW algorithm can achieve a more accurate LST than the RTE algorithm, which is consistent with previous conclusions [11,19]. The uncertainty of the LST algorithm at the DRA, FPK, and TBL sites was greater than 2.5 K, which can be explained by the spatial heterogeneity in the cover types [64]. After excluding three sites with high spatial heterogeneity, the overall biases (RMSEs) of the Landsat 9 LSTs were between −0.18 (1.65) and 0.07 (1.71) K for the SW algorithms, while the bias and RMSE were 1.06 K and 2.38 K, respectively, for the RTE algorithm.
To assess the accuracy of the retrieved Landsat 9 LSTs at the FPK and TBL sites, we used the MODIS LSTs produced by the Temperature-Emissivity Separation (TES) algorithms as a reference for evaluating the performance of the retrieved LSTs. Figure 4 presents the cross-validation result for Landsat 9 LST at the FPK and TBL sites. It can be seen that the retrieved Landsat 9 LSTs have higher uncertainties during the summer months. The RMSEs range from 3.00 to 3.30 K at the FPK site and from 2.36 to 2.71 K at the TBL site for the SW algorithms, while the values are 3.96 K and 3.00 K for the RTE algorithm. As for the bias, the values are between −0.72 K and 0.37 K for the SW algorithms, but the biases are 3.08 K and 2.34 K for the RTE algorithm. It seems that the imprecise LSE could be the reason for the poor performance of the official Landsat 9 LST. Because the official LSEs at the FPK and TBL sites were invariant with the change in the NDVIs, they were unable to accurately reflect the changes in the vegetation cover over time [69,70]. In contrast to the summer, the validation results for the different algorithms show similar accuracy in the autumn, with an RMSE of about 1.20 K. Figure 5 shows the spatial distribution of the difference between the retrieved Landsat 9 LST and the MODIS LST at the FPK site (Figure 5a–d) and the TBL site (Figure 5e–h). As shown in Figure 5, the discrepancies between the MODIS LST and the Landsat 9 LST retrieved by the SW1 and SW3 algorithms exhibit a comparable spatial distribution. Conversely, the differences between the MODIS LST and the Landsat 9 LST retrieved by the SW5 (RTE) algorithm are more pronounced than those retrieved by the SW1 and SW3 algorithms in most regions. The 25th (75th) percentiles of the LST discrepancy between the MODIS LST and the Landsat 9 LST retrieved from the RTE, SW1, SW3, and SW5 algorithms are 0.28 (1.52), −0.66 (0.80),−0.78 (0.69), and −1.13 (0.27) K, respectively, whereas the corresponding values are −0.17 (0.94), −0.72 (0.65),−0.84 (0.54), and −0.57 (0.80) K for Figure 5e–h. The statistical results show that the SW algorithms can achieve more accurate Landsat 9 LSTs than the RTE algorithm compared with the MODIS LST, which is also consistent with previous research [11,15,19].

5. Conclusions

The LST is one of the key parameters in the Earth’s surface energy on regional and global scales. In this paper, we have adapted eleven candidate SW algorithms for LST retrieval from Landsat 9 TIRS-2 data. First, the coefficients of the SW algorithms were derived from the simulation dataset constructed using MODTRAN, emissivity spectra, and different global atmospheric profiles. Simulated data sets derived from CLAR, ERA5, GAPRI, SeeBor, and TIGR were used to calculate the coefficients, while the radiosonde profiles collected from global radiosonde stations were used independently for validation purposes. The NDVI threshold method was then used to estimate the Landsat 9 LSE to facilitate LST retrieval. Finally, the effectiveness of the SW algorithms was assessed using ground measurements collected from SURFRAS sites. The main findings are summarized below:
(1)
Any of the five atmospheric profile databases can be used for simulation when the TWV is less than 4.5 g/cm2. However, the GAPRI profile is recommended for simulation when the TWV is greater than 4.5 g/cm2. The SW5 and SW8 algorithms that do not depend on the emissivity difference exhibit the least sensitivity to emissivity uncertainty but with larger uncertainty in LST retrieval than other SW algorithms.
(2)
At the SURFRAD sites, the candidate SW algorithms exhibit comparable performance, with overall bias (RMSEs) ranging from 1.44 (3.50) to 1.61 (3.73) K compared to 2.82 K and 4.69 K for the RTE algorithm. As most of the TWVs at the SURFRAD sites were below 3.0 g/cm2, more ground measurements with higher TWVs are needed to characterize the performance of the SW algorithms under different atmospheric conditions.
(3)
The consistency of the differences between the MODIS and Landsat LSTs indicates the stability of the SW algorithms over the seasons. The cross-validation results also indicated that the SW algorithms can achieve more accurate Landsat 9 LSTs than the RTE algorithm. This study provides an alternative method for the estimation of Landsat 9 LSTs.

Author Contributions

Conceptualization, X.M.; methodology, data curation, and validation, Q.S., X.M. and L.S.; formal analysis and investigation, Q.S., X.M. and L.S.; writing—original draft preparation, Q.S., X.M. and L.S.; writing—review and editing, Q.S., X.M. and L.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Shandong Provincial Natural Science Foundation (grant numbers ZR2020MD051 and ZR2021QD055).

Data Availability Statement

The Landsat product was downloaded from https://earthexplorer.usgs.gov/ (accessed on 1 June 2024). The ground measurements were downloaded from https://www.esrl.noaa.gov/gmd/grad/surfrad/sitepage.html (accessed on 1 June 2024).

Acknowledgments

The authors wish to thank the Landsat 9 Project Teams at NASA and the USGS. We also thank the USGS EROS Landsat ground system team for producing Landsat Collection 2 data products.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV ≤ 1.5 g/cm2.
Table A1. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV ≤ 1.5 g/cm2.
TWVAlgorithmC0C1C2C3C4C5C6C7
TWV ≤ 1.5SW1−1.1491.0050.171−0.3213.2429.7883.352
SW2−1.2061.0050.171−0.3183.1689.9731.6560.017
SW3−1.1711.2091.24−0.205−0.017−0.225
SW453.5161.0153.4−57.882−2.328−90.52
SW5−1.4851.237−0.23−216.696
SW6−4.3311.0151.13657.644−87.958
SW7−4.1981.0161.12848.251−80.916
SW863.8661.040.18−74.749
SW952.0351.0151.137−56.323−83.669
SW10−4.3311.0151.13657.644−59.136
SW11−4.2631.0151.183−0.02758.247−60.984
Table A2. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV is between 1.5 and 3.0 g/cm2.
Table A2. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV is between 1.5 and 3.0 g/cm2.
TWVAlgorithmC0C1C2C3C4C5C6C7
1.5 < TWV ≤ 3.0SW12.0270.9910.162−0.2894.5024.982−0.142
SW21.5590.9930.159−0.2774.0816.371−4.2870.045
SW32.0791.191.821−0.1990.309−0.209
SW456.51713.842−57.249−2.089−91.909
SW5−3.7182.272−1.259−219.879
SW6−0.73911.81558.767−90.927
SW7−0.62511.81149.136−83.873
SW877.2911.0421.317−89.949
SW956.71511.815−57.416−86.403
SW10−0.73911.81558.767−61.544
SW11−0.71911.823−0.00258.821−61.623
Table A3. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV is between 3.0 and 4.5 g/cm2.
Table A3. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV is between 3.0 and 4.5 g/cm2.
TWVAlgorithmC0C1C2C3C4C5C6C7
3.0 < TWV ≤ 4.5SW17.0060.970.125−0.1795.8255.607−6.667
SW27.0330.9710.121−0.175.4276.546−8.6470.029
SW36.9481.1172.434−0.1473.202−0.135
SW445.4680.9767.065−40.462−4.666−65.858
SW54.4673.253−2.274−227.672
SW64.6450.9772.53150.085−65.413
SW74.7180.9772.52942.2−60.659
SW873.4071.0082.349−77.837
SW953.7890.9772.531−49.123−62.135
SW104.6450.9772.53150.085−40.371
SW114.6460.9762.585−0.00950.221−40.409
Table A4. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV < 4.5 g/cm2.
Table A4. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV < 4.5 g/cm2.
TWVAlgorithmC0C1C2C3C4C5C6C7
4.5 < TWVSW116.3030.9310.066−0.057.5497.287−12.614
SW216.6730.930.064−0.0477.2847.655−13.1980.015
SW316.2420.9883.279−0.0575.908−0.071
SW427.6550.93410.517−12.201−7.256−40.622
SW521.0314.247−3.333−235.971
SW614.6830.9343.46837.12−40.894
SW714.7630.9343.46731.48−38.14
SW867.9980.9423.431−55.77
SW951.220.9343.468−36.523−38.835
SW1014.6830.9343.46837.12−22.334
SW1114.3040.9343.596−0.01637.204−22.263
Table A5. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV is between 0.0 and 10.0 g/cm2.
Table A5. The candidate SW algorithm coefficients obtained from the SeeBor databases when TWV is between 0.0 and 10.0 g/cm2.
TWVAlgorithmC0C1C2C3C4C5C6C7
0.0 < TWV ≤ 10.0SW15.3290.980.161−0.3345.254−8.19912.475
SW2−2.0561.0090.158−0.1962.47−2.851−14.0010.243
SW35.4291.1832.229−0.204−8.078−0.251
SW462.6130.988−5.971−60.0138.151−99.067
SW57.0882.573−1.599−196.261
SW62.4190.991.91954.979−103.642
SW72.5960.991.91845.482−95.275
SW895.8571.0041.671−97.594
SW955.8940.991.919−53.433−98.498
SW102.4190.991.91954.979−76.153
SW11−3.0381.0110.9320.20850.854−48.481

Appendix B

Figure A1. Scatterplots between the air temperature at the first layer and the total water vapor of the six atmospheric profile databases. (a) CLAR, (b) ERA5, (c) GAPRI (d) SeeBor, (e) TIGR and (f) WYO.
Figure A1. Scatterplots between the air temperature at the first layer and the total water vapor of the six atmospheric profile databases. (a) CLAR, (b) ERA5, (c) GAPRI (d) SeeBor, (e) TIGR and (f) WYO.
Remotesensing 16 03633 g0a1aRemotesensing 16 03633 g0a1b

References

  1. Hulley, G.C.; Hughes, C.G.; Hook, S.J. Quantifying uncertainties in land surface temperature and emissivity retrievals from aster and modis thermal infrared data. J. Geophys. Res. Atmos. 2012, 117, D23. [Google Scholar] [CrossRef]
  2. Göttsche, F.-M.; Olesen, F.-S.; Trigo, I.; Bork-Unkelbach, A.; Martin, M. Long term validation of land surface temperature retrieved from msg/seviri with continuous in-situ measurements in africa. Remote Sens. 2016, 8, 410. [Google Scholar] [CrossRef]
  3. Li, H.; Wang, H.; Yang, Y.; Du, Y.; Cao, B.; Bian, Z.; Liu, Q. Evaluation of atmospheric correction methods for the aster temperature and emissivity separation algorithm using ground observation networks in the hiwater experiment. IEEE Trans. Geosci. Remote Sens. 2019, 57, 3001–3014. [Google Scholar] [CrossRef]
  4. Sobrino, J.A.; Julien, Y.; Jiménez-Muñoz, J.-C.; Skokovic, D.; Sòria, G. Near real-time estimation of sea and land surface temperature for msg seviri sensors. Int. J. Appl. Earth Obs. Geoinf. 2020, 89, 102096. [Google Scholar] [CrossRef]
  5. Trigo, I.F.; Ermida, S.L.; Martins, J.P.A.; Gouveia, C.M.; Göttsche, F.-M.; Freitas, S.C. Validation and consistency assessment of land surface temperature from geostationary and polar orbit platforms: Seviri/msg and avhrr/metop. ISPRS J. Photogramm. Remote Sens. 2021, 175, 282–297. [Google Scholar] [CrossRef]
  6. Valor, E.; Caselles, V. Mapping land surface emissivity from ndvi: Application to european, african, and south american areas. Remote Sens. Environ. 1996, 57, 167–184. [Google Scholar] [CrossRef]
  7. Wan, Z.; Dozier, J. A generalized split-window algorithm for retrieving land-surface temperature from space. IEEE Trans. Geosci. Remote Sens. 1996, 34, 892–905. [Google Scholar]
  8. Zeng, Q.; Cheng, J.; Yue, W. Estimating hourly all-sky surface longwave upward radiation using the new generation of chinese geostationary weather satellites fengyun-4a/agri. IEEE Trans. Geosci. Remote Sens. 2024, 62, 1–12. [Google Scholar] [CrossRef]
  9. Eon, R.; Gerace, A.; Falcon, L.; Poole, E.; Kleynhans, T.; Raqueño, N.; Bauch, T. Validation of landsat-9 and landsat-8 surface temperature and reflectance during the underfly event. Remote Sens. 2023, 15, 3370. [Google Scholar] [CrossRef]
  10. Meng, X.; Cheng, J.; Guo, H.; Guo, Y.; Yao, B. Accuracy evaluation of the landsat 9 land surface temperature product. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2022, 15, 8694–8703. [Google Scholar] [CrossRef]
  11. Wang, M.; Li, M.; Zhang, Z.; Hu, T.; He, G.; Zhang, Z.; Wang, G.; Li, H.; Tan, J.; Liu, X. Land surface temperature retrieval from landsat 9 tirs-2 data using radiance-based split-window algorithm. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2023, 16, 1100–1112. [Google Scholar] [CrossRef]
  12. Masek, J.G.; Wulder, M.A.; Markham, B.; McCorkel, J.; Crawford, C.J.; Storey, J.; Jenstrom, D.T. Landsat 9: Empowering open science and applications through continuity. Remote Sens. Environ. 2020, 248, 111968. [Google Scholar] [CrossRef]
  13. Montanaro, M.; McCorkel, J.; Tveekrem, J.; Stauder, J.; Mentzell, E.; Lunsford, A.; Hair, J.; Reuter, D. Landsat 9 thermal infrared sensor 2 (tirs-2) stray light mitigation and assessment. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–8. [Google Scholar] [CrossRef]
  14. Xu, H.; Ren, M.; Lin, M. Cross-comparison of landsat-8 and landsat-9 data: A three-level approach based on underfly images. GIScience Remote Sens. 2024, 61, 2318071. [Google Scholar] [CrossRef]
  15. Niclòs, R.; Perelló, M.; Puchades, J.; Coll, C.; Valor, E. Evaluating landsat-9 tirs-2 calibrations and land surface temperature retrievals against ground measurements using multi-instrument spatial and temporal sampling along transects. Int. J. Appl. Earth Obs. Geoinf. 2023, 125, 103576. [Google Scholar] [CrossRef]
  16. Jiménez-Muñoz, J.C.; Sobrino, J.A.; Skokovic, D.; Mattar, C.; Cristobal, J. Land surface temperature retrieval methods from landsat-8 thermal infrared sensor data. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1840–1843. [Google Scholar] [CrossRef]
  17. Rozenstein, O.; Qin, Z.; Derimian, Y.; Karnieli, A. Derivation of land surface temperature for landsat-8 tirs using a split window algorithm. Sensors 2014, 14, 5768–5780. [Google Scholar] [CrossRef]
  18. Caselles, V.; Rubio, E.; Coll, C.; Valor, E. Thermal band selection for the prism instrument 3. Optimal band configurations. J. Geophys. Res. 1998, 103, 17057–17067. [Google Scholar] [CrossRef]
  19. Ye, X.; Ren, H.; Zhu, J.; Fan, W.; Qin, Q. Split-window algorithm for land surface temperature retrieval from landsat-9 remote sensing images. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
  20. Ye, X.; Liu, R.; Hui, J.; Zhu, J. Land surface temperature estimation from landsat-9 thermal infrared data using ensemble learning method considering the physical radiance transfer process. Land 2023, 12, 1287. [Google Scholar] [CrossRef]
  21. Zhou, J.; Liang, S.; Cheng, J.; Wang, Y.; Ma, J. The glass land surface temperature product. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2019, 12, 493–507. [Google Scholar] [CrossRef]
  22. Yu, Y.; Privette, J.L.; Pinheiro, A.C. Evaluation of split-window land surface temperature algorithms for generating climate data records. IEEE Trans. Geosci. Remote Sens. 2008, 46, 179–192. [Google Scholar] [CrossRef]
  23. Dong, L.; Tang, S.; Wang, F.; Cosh, M.; Li, X.; Min, M. Inversion and validation of fy-4a official land surface temperature product. Remote Sens. 2023, 15, 2437. [Google Scholar] [CrossRef]
  24. Meng, X.; Cheng, J.; Zhao, S.; Liu, S.; Yao, Y. Estimating land surface temperature from landsat-8 data using the noaa jpss enterprise algorithm. Remote Sens. 2019, 11, 155. [Google Scholar] [CrossRef]
  25. Meng, X.; Liu, W.; Cheng, J.; Guo, H. An operational split-window algorithm for retrieving land surface temperature from fengyun-4a agri data. Remote Sens. Lett. 2023, 14, 1206–1217. [Google Scholar] [CrossRef]
  26. Becker, F.; Li, Z.-L. Towards a local split window method over land surfaces. Int. J. Remote Sens. 1990, 11, 369–393. [Google Scholar] [CrossRef]
  27. Wan, Z. New refinements and validation of the collection-6 modis land-surface temperature/emissivity product. Remote Sens. Environ. 2014, 140, 36–45. [Google Scholar] [CrossRef]
  28. Price, J.C. Land surface temperature measurements from the split window channels of the noaa 7 advanced very high resolution radiometer. J. Geophys. Res. Atmos. 1984, 89, 7231–7237. [Google Scholar] [CrossRef]
  29. Yu, Y.; Liu, Y.; Yu, P.; Wang, H. Enterprise Algorithm Theoretical Basis Document for Viirs Land Surface Temperature Production, version 1.0; NOAA: Silver Spring, MD, USA, 2017. [Google Scholar]
  30. Caselles, V.; Coll, C.; Valor, E. Land surface emissivity and temperature determination in the whole hapex-sahel area from avhrr data. Int. J. Remote Sens. 1997, 18, 1009–1027. [Google Scholar] [CrossRef]
  31. Ulivieri, C.; Castronuovo, M.; Francioni, R.; Cardillo, A. A split window algorithm for estimating land surface temperature from satellites. Adv. Space Res. 1994, 14, 59–65. [Google Scholar] [CrossRef]
  32. Vidal, A. Atmospheric and emissivity correction of land surface temperature measured from satellite using ground measurements or satellite data. Int. J. Remote Sens. 1991, 12, 2449–2460. [Google Scholar] [CrossRef]
  33. Ulivieri, C.; Cannizzaro, G. Land surface temperature retrievals from satellite measurements. Acta Astronaut. 1985, 12, 977–985. [Google Scholar] [CrossRef]
  34. Sobrino, J.A.; Li, Z.-L.; Stoll, M.P.; Becker, F. Improvements in the split-window technique for land surface temperature determination. IEEE Trans. Geosci. Remote Sens. 1994, 32, 243–253. [Google Scholar] [CrossRef]
  35. Coll, C.; Valor, E.; Schmugge, T.; Caselles, V. A procedure for estimating the land surface emissivity difference in the avhrr channels 4 and 5. Remote Sens. Appl. Valencia. Area 1997. [Google Scholar]
  36. Sobrino, J.; Li, Z.; Stoll, M.P.; Becker, F. Determination of the surface temperature from atsr data. In Proceedings of the 25th International Symposium on Remote Sensing of Environment, Graz, Austria, 4–8 April 1993; ERIM: Rotterdam, Netherlands, 1993; pp. 4–8. [Google Scholar]
  37. Emami, H.; Mojaradi, B.; Safari, A. A new approach for land surface emissivity estimation using ldcm data in semi-arid areas: Exploitation of the aster spectral library data set. Int. J. Remote Sens. 2016, 37, 5060–5085. [Google Scholar] [CrossRef]
  38. Cheng, J.; Liang, S.; Verhoef, W.; Shi, L.; Liu, Q. Estimating the hemispherical broadband longwave emissivity of global vegetated surfaces using a radiative transfer model. IEEE Trans. Geosci. Remote Sens. 2016, 54, 905–917. [Google Scholar] [CrossRef]
  39. Tang, B.H.; Shao, K.; Li, Z.L.; Wu, H.; Tang, R. An improved ndvi-based threshold method for estimating land surface emissivity using modis satellite data. Int. J. Remote Sens. 2015, 36, 4864–4878. [Google Scholar] [CrossRef]
  40. Meng, X.; Guo, H.; Cheng, J.; Yao, B. Can the era5 reanalysis product improve the atmospheric correction accuracy of landsat series thermal infrared data? IEEE Geosci. Remote Sens. Lett. 2022, 19, 7506805. [Google Scholar] [CrossRef]
  41. Meng, X.; Liu, W.; Cheng, J.; Guo, H.; Yao, B. Estimating hourly land surface temperature from fy-4a agri using an explicitly emissivity-dependent split-window algorithm. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2023, 16, 5474–5487. [Google Scholar] [CrossRef]
  42. Aires, F.; Chédin, A.; Scott, N.A.; Rossow, W.B. A regularized neural net approach for retrieval of atmospheric and surface temperatures with the iasi instrument. J. Appl. Meteorol. Climatol. 2001, 41, 144–159. [Google Scholar] [CrossRef]
  43. Mattar, C.; Durán-Alarcón, C.; Jiménez-Muñoz, J.C.; Santamaría-Artigas, A.; Olivera-Guerra, L.; Sobrino, J.A. Global atmospheric profiles from reanalysis information (gapri): A new database for earth surface temperature retrieval. Int. J. Remote Sens. 2015, 36, 5045–5060. [Google Scholar] [CrossRef]
  44. Borbas, E.E.; Seemann, S.W.; Huang, H.L.; Li, J.; Menzel; Paul, W. Global Profile Training Database for Satellite Regression Retrievals with Estimates of Skin Temperature and Emissivity. In Proceedings of the International ATOVS Study Conference-XIV, Beijing, China, 25–31 May 2005; University of Wisconsin—Madison: Madison, WI, USA, 2005; pp. 763–770. [Google Scholar]
  45. Galve, J.M.; Coll, C.; Caselles, V.; Valor, E. An atmospheric radiosounding database for generating land surface temperature algorithms. IEEE Trans. Geosci. Remote Sens. 2008, 46, 1547–1557. [Google Scholar] [CrossRef]
  46. Jiang, G.-M.; Zhou, W.; Liu, R. Development of split-window algorithm for land surface temperature estimation from the virr/fy-3a measurements. IEEE Geosci. Remote Sens. Lett. 2013, 10, 952–956. [Google Scholar] [CrossRef]
  47. Meng, X.; Cheng, J.; Liang, S. Estimating land surface temperature from feng yun-3c/mersi data using a new land surface emissivity scheme. Remote Sens. 2017, 9, 1247. [Google Scholar] [CrossRef]
  48. Tang, B.-H. Nonlinear split-window algorithms for estimating land and sea surface temperatures from simulated chinese gaofen-5 satellite data. IEEE Trans. Geosci. Remote Sens. 2018, 56, 6280–6289. [Google Scholar] [CrossRef]
  49. Ma, J.; Zhou, J.; Göttsche, F.-M.; Liang, S.; Wang, S.; Li, M. A global long-term (1981–2000) land surface temperature product for noaa avhrr. Earth Syst. Sci. Data 2020, 12, 3247–3268. [Google Scholar] [CrossRef]
  50. Li, H.; Yang, Y.; Li, R.; Wang, H.; Cao, B.; Bian, Z.; Hu, T.; Du, Y.; Sun, L.; Liu, Q. Comparison of the musyq and modis collection 6 land surface temperature products over barren surfaces in the heihe river basin, china. IEEE Trans. Geosci. Remote Sens. 2019, 57, 8081–8094. [Google Scholar] [CrossRef]
  51. Niclòs, R.; Galve, J.M.; Valiente, J.A.; Estrela, M.J.; Coll, C. Accuracy assessment of land surface temperature retrievals from msg2-seviri data. Remote Sens. Environ. 2011, 115, 2126–2140. [Google Scholar] [CrossRef]
  52. Tang, K.; Zhu, H.; Ni, P.; Li, R.; Fan, C. Retrieving land surface temperature from chinese fy-3d mersi-2 data using an operational split window algorithm. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 6639–6651. [Google Scholar] [CrossRef]
  53. Ermida, S.L.; Trigo, I.F. A comprehensive clear-sky database for the development of land surface temperature algorithms. Remote Sens. 2022, 14, 2329. [Google Scholar] [CrossRef]
  54. Zhang, G.; Li, D.; Li, H.; Xu, Z.; Hu, Z.; Zeng, J.; Yang, Y.; Jia, H. Improving hj-1b/irs lst retrieval of the generalized single-channel algorithm with refined era5 atmospheric profile database. Remote Sens. 2023, 15, 5092. [Google Scholar] [CrossRef]
  55. Meng, X.; Cheng, J. Evaluating eight global reanalysis products for atmospheric correction of thermal infrared sensor—Application to landsat 8 tirs10 data. Remote Sens. 2018, 10, 474. [Google Scholar] [CrossRef]
  56. Meng, X.; Cheng, J. Estimating land and sea surface temperature from cross-calibrated chinese gaofen-5 thermal infrared data using split-window algorithm. IEEE Geosci. Remote Sens. Lett. 2020, 17, 509–513. [Google Scholar] [CrossRef]
  57. Meng, Y.; Zhou, J.; Göttsche, F.-M.; Tang, W.; Martins, J.; Perez-Planells, L.; Ma, J.; Wang, Z. Investigation and validation of two all-weather land surface temperature products with in-situ measurements. Geo-Spat. Inf. Sci. 2023, 27, 670–682. [Google Scholar] [CrossRef]
  58. Sayler, K.; Glynn, T. Landsat 9 Data Users Handbook; Earth Resources Observation and Science (EROS) Center: Sioux Falls, SD, USA, 2022. Available online: https://www.usgs.gov/media/files/landsat-9-data-users-handbook (accessed on 1 June 2024).
  59. Zhang, Z.; He, G.; Wang, M.; Long, T.; Wang, G.; Zhang, X. Validation of the generalized single-channel algorithm using landsat 8 imagery and surfrad ground measurements. Remote Sens. Lett. 2016, 7, 810–816. [Google Scholar] [CrossRef]
  60. Duan, S.-B.; Li, Z.-L.; Li, H.; Göttsche, F.-M.; Wu, H.; Zhao, W.; Leng, P.; Zhang, X.; Coll, C. Validation of collection 6 modis land surface temperature product using in situ measurements. Remote Sens. Environ. 2019, 225, 16–29. [Google Scholar] [CrossRef]
  61. Martin, M.; Ghent, D.; Pires, A.; Göttsche, F.-M.; Cermak, J.; Remedios, J. Comprehensive in situ validation of five satellite land surface temperature data sets over multiple stations and years. Remote Sens. 2019, 11, 479. [Google Scholar] [CrossRef]
  62. Zeng, Q.; Cheng, J.; Sun, H.; Dong, S. An integrated framework for estimating the hourly all-time cloudy-sky surface long-wave downward radiation for fengyun-4a/agri. Remote Sens. Environ. 2024, 312, 114319. [Google Scholar] [CrossRef]
  63. Guillevic, P.C.; Biard, J.C.; Hulley, G.C.; Privette, J.L.; Hook, S.J.; Olioso, A.; Göttsche, F.M.; Radocinski, R.; Román, M.O.; Yu, Y.; et al. Validation of land surface temperature products derived from the visible infrared imaging radiometer suite (viirs) using ground-based and heritage satellite measurements. Remote Sens. Environ. 2014, 154, 19–37. [Google Scholar] [CrossRef]
  64. Malakar, N.K.; Hulley, G.C.; Hook, S.J.; Laraby, K.; Cook, M.; Schott, J.R. An operational land surface temperature product for landsat thermal data: Methodology and validation. IEEE Trans. Geosci. Remote Sens. 2018, 56, 5717–5735. [Google Scholar] [CrossRef]
  65. Ma, J.; Zhou, J.; Göttsche, F.-M.; Wang, Z.; Wu, H.; Tang, W.; Li, M.; Liu, S. An atmospheric influence correction method for longwave radiation-based in-situ land surface temperature. Remote Sens. Environ. 2023, 293, 113611. [Google Scholar] [CrossRef]
  66. Cheng, J.; Liang, S.; Yao, Y.; Zhang, X. Estimating the optimal broadband emissivity spectral range for calculating surface longwave net radiation. IEEE Geosci. Remote Sens. Lett. 2013, 10, 401–405. [Google Scholar] [CrossRef]
  67. Barsi, J.A.; Montanaro, M.; Thome, K.; Raqueno, N.G.; Hook, S.; Anderson, C.H.; Micijevic, E. Early Radiometric Performance of Landsat-9 Thermal Infrared Sensor. In Proceedings of the Earth Observing Systems XXVII, San Diego, CA, USA, 21–26 August 2022; SPIE: Washington, DC, USA, 2022; pp. 223–233. [Google Scholar]
  68. Du, C.; Ren, H.; Qin, Q.; Meng, J.; Zhao, S. A practical split-window algorithm for estimating land surface temperature from landsat 8 data. Remote Sens. 2015, 7, 647–665. [Google Scholar] [CrossRef]
  69. Duan, S.-B.; Li, Z.-L.; Zhao, W.; Wu, P.; Huang, C.; Han, X.-J.; Gao, M.; Leng, P.; Shang, G. Validation of landsat land surface temperature product in the conterminous united states using in situ measurements from surfrad, arm, and ndbc sites. Int. J. Digit. Earth 2021, 14, 640–660. [Google Scholar] [CrossRef]
  70. Cheng, J.; Meng, X.; Dong, S.; Liang, S. Generating the 30-m land surface temperature product over continental china and USA from landsat 5/7/8 data. Sci. Remote Sens. 2021, 4, 100032. [Google Scholar] [CrossRef]
Figure 1. Performance of the eleven candidate SW algorithms under different total water vapor (TWV) subranges: (a) 0 < TWV < 1.5 g/cm2, (b) 1.5 < TWV < 3.0 g/cm2, (c) 3.0 < TWV < 4.5 g/cm2, (d) 4.5 < TWV < 10 g/cm2, and (e) 0 < TWV < 10 g/cm2.
Figure 1. Performance of the eleven candidate SW algorithms under different total water vapor (TWV) subranges: (a) 0 < TWV < 1.5 g/cm2, (b) 1.5 < TWV < 3.0 g/cm2, (c) 3.0 < TWV < 4.5 g/cm2, (d) 4.5 < TWV < 10 g/cm2, and (e) 0 < TWV < 10 g/cm2.
Remotesensing 16 03633 g001aRemotesensing 16 03633 g001b
Figure 2. Performance of the eleven candidate SW algorithms before and after adding random noise under different total water vapor (TWV) subranges: (ac) 0 < TWV < 1.5 g/cm2, (df) 4.5 < TWV < 10 g/cm2. The random noise is added in the emissivity (a,d), brightness temperature (b,e), and TWV (c,f).
Figure 2. Performance of the eleven candidate SW algorithms before and after adding random noise under different total water vapor (TWV) subranges: (ac) 0 < TWV < 1.5 g/cm2, (df) 4.5 < TWV < 10 g/cm2. The random noise is added in the emissivity (a,d), brightness temperature (b,e), and TWV (c,f).
Remotesensing 16 03633 g002
Figure 3. The bias (a) and RMSE (b) between the retrieved Landsat 9 LST and the in situ LST for all matchups at the SURFRAD sites.
Figure 3. The bias (a) and RMSE (b) between the retrieved Landsat 9 LST and the in situ LST for all matchups at the SURFRAD sites.
Remotesensing 16 03633 g003
Figure 4. The bias (a) and RMSE (b) between the retrieved Landsat 9 LST and the MODIS LST.
Figure 4. The bias (a) and RMSE (b) between the retrieved Landsat 9 LST and the MODIS LST.
Remotesensing 16 03633 g004
Figure 5. Spatial distribution of the difference between the retrieved Landsat 9 LST and the MODIS LST. (a,e) LST retrieved from RTE algorithm, (b,f) LST retrieved from SW1 algorithm, (c,g) LST retrieved from SW3 algorithm, and (d,h) LST retrieved from SW5 algorithm.
Figure 5. Spatial distribution of the difference between the retrieved Landsat 9 LST and the MODIS LST. (a,e) LST retrieved from RTE algorithm, (b,f) LST retrieved from SW1 algorithm, (c,g) LST retrieved from SW3 algorithm, and (d,h) LST retrieved from SW5 algorithm.
Remotesensing 16 03633 g005aRemotesensing 16 03633 g005b
Table 1. Summary of the eleven candidate LST algorithms used, including names, formulas, and references.
Table 1. Summary of the eleven candidate LST algorithms used, including names, formulas, and references.
NameFormulaReferences
SW 1 T s = C 0 + ( C 1 + C 2 1 ε ε + C 3 ε ε 2 ) T 10 + T 11 2 + ( C 4 + C 5 1 ε ε + C 6 ε ε 2 ) T 10 T 11 2 [7,26]
SW 2 T s = C 0 + C 1 + C 2 1 ε ε + C 3 ε ε 2 T 10 + T 11 2 + C 4 + C 5 1 ε ε + C 6 ε ε 2 T 10 T 11 2 + C 7 ( T 10 T 11 ) 2 [27]
SW 3 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ε 10 T 10 + C 4 ( 1 ε 10 ) T 10 T 11 + C 5 T 11 ε [28]
SW 4 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ε + C 4 ε T 10 T 11 + C 5 ε [29]
SW 5 T s = C 0 + C 1 T 10 ε + C 2 T 11 ε + C 3 1 ε ε [30]
SW 6 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ( 1 ε ) + C 4 ε [31]
SW 7 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 1 ε ε + C 4 ε ε 2 [32]
SW 8 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ε [33]
SW 9 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ε + C 4 ε ε [34]
SW 10 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ( 1 ε 10 ) + C 4 ε [35]
SW 11 T s = C 0 + C 1 T 10 + C 2 T 10 T 11 + C 3 ( T 10 T 11 ) 2 + C 4 ( 1 ε 10 ) + C 5 ε [36]
Note: C 0 C 7 are the split-window algorithm coefficients; T 10 and T 11 are the brightness temperatures of channels 10 and 11, respectively; ε = 0.5 * ( ε 10 + ε 11 ) , ε = ε 10 ε 11 , ε 10 , and ε 11 are the emissivity values of channels 10 and 11, respectively.
Table 2. The regression coefficients for the emissivity calculation.
Table 2. The regression coefficients for the emissivity calculation.
Channela1a2a3a4a5a6a7
100.9766−0.10680.1524−0.0398−0.05680.0791−0.0712
110.98200.0265−0.05650.0574−0.06630.0761−0.0603
Table 3. Detailed information on seven SURFRAD sites.
Table 3. Detailed information on seven SURFRAD sites.
SiteAbbreviationLatitudeLongitudeLand CoverPath/RowLandsat PeriodPath/Row
BondvilleBND40.052−88.373cropland023/0322021/10~2024/05023/032
GoodwinCreekGWN34.255−89.873grassland023/036023/036
Penn. State Univ.PSU40.720−77.931cropland016/032016/032
Sioux FallsSXF43.734−96.623grassland029/030029/030
Fort PeckFPK48.308−105.102grassland035/026035/026
TableMountainTBL40.126−105.238grassland033/032033/032
Desert RockDRA36.623−116.020shrubland040/035040/035
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Su, Q.; Meng, X.; Sun, L. Investigation and Validation of Split-Window Algorithms for Estimating Land Surface Temperature from Landsat 9 TIRS-2 Data. Remote Sens. 2024, 16, 3633. https://doi.org/10.3390/rs16193633

AMA Style

Su Q, Meng X, Sun L. Investigation and Validation of Split-Window Algorithms for Estimating Land Surface Temperature from Landsat 9 TIRS-2 Data. Remote Sensing. 2024; 16(19):3633. https://doi.org/10.3390/rs16193633

Chicago/Turabian Style

Su, Qinghua, Xiangchen Meng, and Lin Sun. 2024. "Investigation and Validation of Split-Window Algorithms for Estimating Land Surface Temperature from Landsat 9 TIRS-2 Data" Remote Sensing 16, no. 19: 3633. https://doi.org/10.3390/rs16193633

APA Style

Su, Q., Meng, X., & Sun, L. (2024). Investigation and Validation of Split-Window Algorithms for Estimating Land Surface Temperature from Landsat 9 TIRS-2 Data. Remote Sensing, 16(19), 3633. https://doi.org/10.3390/rs16193633

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop