Investigation and validation of algorithms for estimating land surface temperature from Sentinel-3 SLSTR data

Land surface temperature (LST) is an important indicator of global ecological environment and climate change. The Sea and Land Surface Temperature Radiometer (SLSTR) onboard the recently launched Sentinel-3 satellites provides high-quality observations for estimating global LST. The algorithm of the official SLSTR LST product is a split-window algorithm (SWA) that implicitly assumes and utilizes knowledge of land surface emissivity (LSE). The main objective of this study is to investigate alternative SLSTR LST retrieval algorithms with an explicit use of LSE. Seventeen widely accepted SWAs, which explicitly utilize LSE, were selected as candidate algorithms. First, the SWAs were trained using a comprehensive global simulation dataset. Then, using simulation data as well as in-situ LST, the SWAs were evaluated according to their sensitivity and accuracy: eleven algorithms showed good training accuracy and nine of them exhibited low sensitivity to uncertainties in LSE and column water vapor content. Evaluation based on two global simulation datasets and a regional simulation dataset showed that these nine SWAs had similar accuracy with negligible systematic errors and RMSEs lower than 1.0 K. Validation based on in-situ LST obtained for six sites further confirmed the similar accuracies of the SWAs, with the lowest RMSE ranges of 1.57–1.62 K and 0.49−0.61 K for Gobabeb and Lake Constance, respectively. While the best two SWAs usually yielded good accuracy, the official SLSTR LST generally had lower accuracy. The SWAs identified and described in this study may serve as alternative algorithms for retrieving LST products from SLSTR data

A global historical twice-daily (daytime and nighttime) land surface temperature dataset produced by Advanced Very High Resolution Radiometer observations from 1981 to 2021

Land surface temperature (LST) is a key variable for monitoring and evaluating global long-term climate change. However, existing satellite-based twice-daily LST products only date back to 2000, which makes it difficult to obtain robust long-term temperature variations. In this study, we developed the first global historical twice-daily LST dataset (GT-LST), with a spatial resolution of 0.05∘, using Advanced Very High Resolution Radiometer (AVHRR) Level-1b Global Area Coverage (GAC) data from 1981 to 2021. The GT-LST product was generated using four main processes: (1) GAC data reading, calibration, and preprocessing using open-source Python libraries; (2) cloud detection using the AVHRR-Phase I algorithm; (3) land surface emissivity estimation using an improved method considering annual land cover changes; (4) LST retrieval based on a nonlinear generalized split-window algorithm. Validation with in situ measurements from Surface Radiation Budget (SURFRAD) sites and Baseline Surface Radiation Network sites showed that the overall root-mean-square errors (RMSEs) of GT-LST varied from 1.6 to 4.0 K, and nighttime LSTs were typically better than daytime LSTs. Intercomparison with the Moderate Resolution Imaging Spectroradiometer LST products (MYD11A1 and MYD21A1) revealed that the overall root-mean-square difference (RMSD) was approximately 3.0 K. Compared with MYD11A1 LST, GT-LST was overestimated, and relatively large RMSDs were obtained during the daytime, spring, and summer, whereas the significantly smaller positive bias was obtained between GT-LST and MYD21A1 LST. Furthermore, we compared our newly generated dataset with a global AVHRR daytime LST product at the selected measurements of SURFRAD sites (i.e., measurements of these two satellite datasets were valid), which revealed similar accuracies for the two datasets. However, GT-LST can additionally provide nighttime LST, which can be combined with daytime observations estimating relatively accurate monthly mean LST, with an RMSE of 2.7 K. Finally, we compared GT-LST with a regional twice-daily AVHRR LST product over continental Africa in different seasons, with RMSDs ranging from 2.1 to 4.3 K. Considering these advantages, the proposed dataset provides a better data source for a range of research applications. GT-LST is freely available at https://doi.org/10.5281/zenodo.7113080 (1981–2000) (Li et al., 2022a), https://doi.org/10.5281/zenodo.7134158 (2001–2005) (Li et al., 2022b), and https://doi.org/10.5281/zenodo.7813607 (2006–2021) (J. H. Li et al., 2023).</p

Land Surface Temperature Measurements form EOS MODIS Data

We have developed a physics-based land-surface temperature (LST) algorithm for simultaneously retrieving surface band-averaged emissivities and temperatures from day/night pairs of MODIS (Moderate Resolution Imaging Spectroradiometer) data in seven thermal infrared bands. The set of 14 nonlinear equations in the algorithm is solved with the statistical regression method and the least-squares fit method. This new LST algorithm was tested with simulated MODIS data for 80 sets of band-averaged emissivities calculated from published spectral data of terrestrial materials in wide ranges of atmospheric and surface temperature conditions. Comprehensive sensitivity and error analysis has been made to evaluate the performance of the new LST algorithm and its dependence on variations in surface emissivity and temperature, upon atmospheric conditions, as well as the noise-equivalent temperature difference (NE(Delta)T) and calibration accuracy specifications of the MODIS instrument. In cases with a systematic calibration error of 0.5%, the standard deviations of errors in retrieved surface daytime and nighttime temperatures fall between 0.4-0.5 K over a wide range of surface temperatures for mid-latitude summer conditions. The standard deviations of errors in retrieved emissivities in bands 31 and 32 (in the 10-12.5 micrometer IR spectral window region) are 0.009, and the maximum error in retrieved LST values falls between 2-3 K. Several issues related to the day/night LST algorithm (uncertainties in the day/night registration and in surface emissivity changes caused by dew occurrence, and the cloud cover) have been investigated. The LST algorithms have been validated with MODIS Airborne Simulator (MAS) dada and ground-based measurement data in two field campaigns conducted in Railroad Valley playa, NV in 1995 and 1996. The MODIS LST version 1 software has been delivered

Supplement of the radiance-based method to validate satellite-derived land surface temperature products over heterogeneous land surfaces

Land surface temperature (LST) retrieved from satellite remote sensing data has become a key parameter in research on global environmental change; therefore, the acquisition of accurate satellite-derived LST information is crucial for the diagnosis and analysis of global change. However, it is relatively difficult to obtain the true value of a pixel due to the scale mismatch between in situ measurements and satellite-based observations, especially for commonly heterogeneous and nonisothermal land areas, which greatly increases the difficulty in estimating pixel-representative LST values from in situ measurements for validation of satellite-based LST products. In this study, a supplemented radiance-based (SR-based) validation method was developed to evaluate the latest moderate resolution imaging spectroradiometer (MODIS) Collection 6 Level 2 daily LST/land surface emissivity (LSE) products over a heterogeneous and nonisothermal region of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project, West China. In the SR-based framework, pixel-representative LST values are simulated by the MODTRAN model from the corresponding in situ measurements, such as LSE and atmospheric profile measurements, to evaluate the MODIS LST products. The validation results show that the MODIS daytime LST products from the Aqua satellite (MYD11_L2) have a greater accuracy than those from the Terra satellite (MOD11_L2). Analyses of the effect factors indicate a strong correlation between the errors in the MOD11_L2 LST product and the corresponding difference in the MODIS brightness temperature between bands 31 and 32. Although the requirement of synchronous or quasisynchronous in situ measurements for the validated LST products may limit the applicability of the SR-based method, it is still an effective and simple method for validating satellite-derived LST products over mixed pixels. Our method is an indispensable supplement for the validation methods of satellite-derived LST products, and it can be applied in West China and other areas with heterogeneous land surfaces

Remote Sensing Monitoring of Land Surface Temperature (LST)

This book is a collection of recent developments, methodologies, calibration and validation techniques, and applications of thermal remote sensing data and derived products from UAV-based, aerial, and satellite remote sensing. A set of 15 papers written by a total of 70 authors was selected for this book. The published papers cover a wide range of topics, which can be classified in five groups: algorithms, calibration and validation techniques, improvements in long-term consistency in satellite LST, downscaling of LST, and LST applications and land surface emissivity research

Comparison of methods to estimate lake-surface-water temperature using Landsat 7 ETM+ and MODIS imagery : case study of a large shallow subtropical lake in southern Brazil

Water temperature regulates many processes in lakes; therefore, evaluating it is essential to understand its ecological status and functioning, and to comprehend the impact of climate change. Although few studies assessed the accuracy of individual sensors in estimating lake-surface-water temperature (LSWT), comparative analysis considering different sensors is still needed. This study evaluated the performance of two thermal sensors, MODIS and Landsat 7 ETM+, and used Landsat methods to estimate the SWT of a large subtropical lake. MODIS products MOD11 LST and MOD28 SST were used for comparison. For the Landsat images, the radiative transfer equation (RTE), using NASA’s Atmospheric Correction Parameter Calculator (AtmCorr) parameters, was compared with the single-channel algorithm in different approaches. Our results showed that MOD11 obtained the highest accuracy (RMSE of 1.05 C), and is the recommended product for LSWT studies. For Landsat-derived SWT, AtmCorr obtained the highest accuracy (RMSE of 1.07 C) and is the recommended method for small lakes. Sensitivity analysis showed that Landsat-derived LSWT using the RTE is very sensitive to atmospheric parameters and emissivity. A discussion of the main error sources was conducted. We recommend that similar tests be applied for Landsat imagery on different lakes, further studies on algorithms to correct the cool-skin effect in inland waters, and tests of different emissivity values to verify if it can compensate for this effect, in an effort to improve the accuracy of these estimates

Linking thermal variability and change to urban growth in Harare Metropolitan City using remotely sensed data.

Doctor of Philosophy in Environmental Science. University of KwaZulu-Natal. Pietermaritzburg, 2017.Urban growth, which involves Land Use and Land Cover Changes (LULCC), alters land surface thermal properties. Within the framework of rapid urban growth and global warming, land surface temperature (LST) and its elevation have potential significant socio-economic and environmental implications. Hence the main objectives of this study were to (i) map urban growth, (ii) link urban growth with indoor and outdoor thermal conditions and (iii) estimate implications of thermal trends on household energy consumption as well as predict future urban growth and temperature patterns in Harare Metropolitan, Zimbabwe. To achieve these objectives, broadband multi-spectral Landsat 5, 7 and 8, in-situ LULC observations, air temperature (Ta) and humidity data were integrated. LULC maps were obtained from multi-spectral remote sensing data and derived indices using the Support Vector Machine Algorithm, while LST were derived by applying single channel and split window algorithms. To improve remote sensing based urban growth mapping, a method of combining multi-spectral reflective data with thermal data and vegetation indices was tested. Vegetation indices were also combined with socio-demographic data to map the spatial distribution of heat vulnerability in Harare. Changes in outdoor human thermal discomfort in response to seasonal LULCC were evaluated, using the Discomfort Index (DI) derived parsimoniously from LST retrieved from Landsat 8 data. Responses of LST to long term urban growth were analysed for the period from 1984 to 2015. The implications of urban growth induced temperature changes on household air-conditioning energy demand were analysed using Landsat derived land surface temperature based Degree Days. Finally, the Cellular Automata Markov Chain (CAMC) analysis was used to predict future landscape transformation at 10-year time steps from 2015 to 2045. Results showed high overall accuracy of 89.33% and kappa index above 0.86 obtained, using Landsat 8 bands and indices. Similar results were observed when indices were used as stand-alone dataset (above 80%). Landsat 8 derived bio-physical surface properties and socio-demographic factors, showed that heat vulnerability was high in over 40% in densely built-up areas with low-income when compared to “leafy” suburbs. A strong spatial correlation (α = 0.61) between heat vulnerability and surface temperatures in the hot season was obtained, implying that LST is a good indicator of heat vulnerability in the area. LST based discomfort assessment approach retrieved DI with high accuracy as indicated by mean percentage error of less than 20% for each sub-season. Outdoor thermal discomfort was high in hot dry season (mean DI of 31oC), while the post rainy season was the most comfortable (mean DI of 19.9oC). During the hot season, thermal discomfort was very low in low density residential areas, which are characterised by forests and well maintained parks (DI ≤27oC). Long term changes results showed that high density residential areas increased by 92% between 1984 and 2016 at the expense of cooler green-spaces, which decreased by 75.5%, translating to a 1.98oC mean surface temperature increase. Due to surface alterations from urban growth between 1984 and 2015, LST increased by an average of 2.26oC and 4.10oC in the cool and hot season, respectively. This decreased potential indoor heating energy needed in the cool season by 1 degree day and increased indoor cooling energy during the hot season by 3 degree days. Spatial analysis showed that during the hot season, actual energy consumption was low in high temperature zones. This coincided with areas occupied by low income strata indicating that they do not afford as much energy and air conditioning facilities as expected. Besides quantifying and strongly relating with energy requirement, degree days provided a quantitative measure of heat vulnerability in Harare. Testing vegetation indices for predictive power showed that the Urban Index (UI) was comparatively the best predictor of future urban surface temperature (r = 0.98). The mean absolute percentage error of the UI derived temperature was 5.27% when tested against temperature derived from thermal band in October 2015. Using UI as predictor variable in CAMC analysis, we predicted that the low surface temperature class (18-28oC) will decrease in coverage, while the high temperature category (36-45oC) will increase in proportion covered from 42.5 to 58% of city, indicating further warming as the city continues to grow between 2015 and 2040. Overall, the findings of this study showed that LST, human thermal comfort and air-conditioning energy demand are strongly affected by seasonal and urban growth induced land cover changes. It can be observed that urban greenery and wetlands play a significant role of reducing LST and heat transfer between the surface and lower atmosphere and LST may continue unless effective mitigation strategies, such as effective vegetation cover spatial configuration are adopted. Limitations to the study included inadequate spatial and low temporal resolution of Landsat data, few in-situ observations of temperature and LULC classification which was area specific thus difficult for global comparison. Recommendations for future studies included data merging to improve spatial and temporal representation of remote sensing data, resource mobilization to increase urban weather station density and image classification into local climate zones which are of easy global interpretation and comparison

Flood Resilience Assessment of New Orleans after Hurricane Katrina based on Thermal and Vegetation Index Image Time Series

Resilience is a concept with increasing importance in modern risk management because of its role in reducing risks of unpreventable disasters. Previous resilience assessment studies often require extensive surveys of various social, economic, and psychological data or incorporate remote sensing data as one of the complicated physical and social parameters for assessment models. Limited data accessibility to such data due to funding, time, and labor intensity is a major challenge for their wider applications. Therefore, this study proposes the hypothesis that the overall resilience of an urban area to disturbances of natural disasters can be reflected through the time series change sequences of thermal and vegetation index from satellite images. This is because the vegetation index reflects the recoverability of vegetated areas, and thermal change pattern is a reflection of land-cover and land-use changes and energy consumption, which is the end result of various impacts such as social, economic, and physical factors. Specifically, this study introduced a rapid and objective flood resilience assessment method through time series classification based on thermal feature and vegetation index. The method first used unsupervised classification methods to identify potential flood impact levels and conducted supervised classification to obtain a more accurate classification result. Finally, the derived impact levels were classified as flood resilience levels, which are beneficiary for flood preparation and resource allocation for the local and federal government