Comparison of Three Methods for Estimating Land Surface Temperature from Landsat 8-TIRS Sensor Data

Abstract

After Landsat 8 was launched in 2013, it was observed that for Thermal Infrared sensor (TIRS) bands, radiance from outside of an instrument’s field-of-view produced a non-uniform ghost signal across the focal plane that varied depending on the out-of-scene content (i.e., the stray light effect). A new stray light correction algorithm (SLCA) is currently operational and has been implemented into the United States Geological Survey (USGS) ground system since February 2017. The SLCA has also been applied to reprocess historical Landsat 8 scenes. After approximately two years of SLCA implementation, more land surface temperature (LST) validation studies are required to check the effect of correction in the estimation of LST from different retrieval algorithms. For this purpose, three different LST estimation method algorithms (i.e., the radiative transfer equation (RTE), single-channel algorithm (SCA), and split-window algorithm (SWA)) have been assessed. The study site is located on the campus of the University of Balearic Islands on the island of Mallorca (Spain) in the western Mediterranean Sea. The site is considered a heterogeneous area that is composed of different types of surfaces, such as buildings, asphalt roads, farming areas, sloped terrains, orange fields, almond trees, lawns, and some natural vegetation regions. Data from 21 scenes, which were acquired by the Landsat 8-TIRS sensor and extracted from a 100 × 100 m2 pixel, were used to retrieve the LST with different algorithms; then, they were compared with in situ LST measurements from a broadband thermal infrared radiometer located on the same Landsat 8 pixel. The results show good performances of the three methods, with the SWA showing the lowest observed RMSE (within 1.6–2 K), whereas the SCA applied to the TIRS band 10 (10 µm) was also appropriate, with a RMSE ranging within 2.0–2.3 K. The LST estimates using the RTE algorithm display the highest observed RMSE values (within 2.0–3.6 K) of all of the compared methods, but with an almost unbiased value of −0.1 K for the case of techniques applied to band 10 data. The SWAs are the preferred method to estimate the LST in our study area. However, further validation studies around the world are required.