Estimation of Land Surface Temperature in an Agricultural Region of Bangladesh from Landsat 8: Intercomparison of Four Algorithms.

Abstract

The presence of two thermal bands in Landsat 8 brings the opportunity to use either one or both of these bands to retrieve Land Surface Temperature (LST). In order to compare the performances of existing algorithms, we used four methods to retrieve LST from Landsat 8 and made an intercomparison among them. Apart from the direct use of the Radiative Transfer Equation (RTE), Single-Channel Algorithm and two Split-Window Algorithms were used taking an agricultural region in Bangladesh as the study area. The LSTs retrieved in the four methods were validated in two ways: first, an indirect validation against reference LST, which was obtained in the Atmospheric and Topographic CORection (ATCOR) software module; second, cross-validation with Terra MODerate Resolution Imaging Spectroradiometer (MODIS) daily LSTs that were obtained from the Application for Extracting and Exploring Analysis Ready Samples (AEEARS) online tool. Due to the absence of LST-monitoring radiosounding instruments surrounding the study area, in situ LSTs were not available; hence, validation of satellite retrieved LSTs against in situ LSTs was not performed. The atmospheric parameters necessary for the RTE-based method, as well as for other methods, were calculated from the National Centers for Environmental Prediction (NCEP) database using an online atmospheric correction calculator with MODerate resolution atmospheric TRANsmission (MODTRAN) codes. Root-mean-squared-error (RMSE) against reference LST, as well as mean bias error against both reference and MODIS daily LSTs, was used to interpret the relative accuracy of LST results. All four methods were found to result in acceptable LST products, leaving atmospheric water vapor content (w) as the important determinant for the precision result. Considering a set of several Landsat 8 images of different dates, Jiménez-Muñoz et al.’s (2014) Split-Window algorithm was found to result in the lowest mean RMSE of . Du et al.’s (2015) Split-Window algorithm resulted in mean RMSE of . The RTE-based direct method and the Single-Channel algorithm provided the mean RMSE of and , respectively. For Du et al.’s algorithm, the w range of to g cm−2 was considered, whereas for the other three methods, w values as retrieved from the NCEP database were considered for corresponding images. Land surface emissivity was retrieved through the Normalized Difference Vegetation Index (NDVI)-threshold method. This intercomparison study provides an LST retrieval methodology for Landsat 8 that involves four algorithms. It proves that (i) better LST results can be obtained using both thermal bands of Landsat 8; (ii) the NCEP database can be used to determine atmospheric parameters using the online calculator; (iii) MODIS daily LSTs from AEEARS can be used efficiently in cross-validation and intercomparison of Landsat 8 LST algorithms; and (iv) when in situ LST data are not available, the ATCOR-derived LSTs can be used for indirect verification and intercomparison of Landsat 8 LST algorithms.