Retrieval of land surface temperature (LST) from landsat TM6 and TIRS data by single channel radiative transfer algorithm using satellite and ground-based inputs

The present study proposes land surface temperature (LST) retrieval from satellite-based thermal IR data by single channel radiative transfer algorithm using atmospheric correction parameters derived from satellite-based and in-situ data and land surface emissivity (LSE) derived by a hybrid LSE model. For example, atmospheric transmittance (τ) was derived from Terra MODIS spectral radiance in atmospheric window and absorption bands, whereas the atmospheric path radiance and sky radiance were estimated using satellite- and ground-based in-situ solar radiation, geographic location and observation conditions. The hybrid LSE model which is coupled with ground-based emissivity measurements is more versatile than the previous LSE models and yields improved emissivity values by knowledge-based approach. It uses NDVI-based and NDVI Threshold method (NDVI^THM) based algorithms and field-measured emissivity values. The model is applicable for dense vegetation cover, mixed vegetation cover, bare earth including coal mining related land surface classes. The study was conducted in a coalfield of India badly affected by coal fire for decades. In a coal fire affected coalfield, LST would provide precise temperature difference between thermally anomalous coal fire pixels and background pixels to facilitate coal fire detection and monitoring. The derived LST products of the present study were compared with radiant temperature images across some of the prominent coal fire locations in the study area by graphical means and by some standard mathematical dispersion coefficients such as coefficient of variation, coefficient of quartile deviation, coefficient of quartile deviation for 3rd quartile vs. maximum temperature, coefficient of mean deviation (about median) indicating significant increase in the temperature difference among the pixels. The average temperature slope between adjacent pixels, which increases the potential of coal fire pixel detection from background pixels, is significantly larger in the derived LST products than the corresponding radiant temperature images.     

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.     

陆地卫星TM6波段范围内地表比辐射率的估计

地表比辐射率是用热红外波段遥感数据反演地表温度的关键参数。目前,应用陆地卫星TM6波段数据反演地表温度共有3种算法,即大气校正法、单窗算法和单通道算法。这3种算法都需要TM6波段范围内的地表比辐射率作为地表温度反演参数。本文首先简介这3种反演算法;然后着重探讨TM6波段地表比辐射率估计方法;最后,利用这一方法对山东省陵县附近农田地区进行地表比辐射率估计和地表温度反演。结果表明,该方法能获得较合理的地表温度反演结果。     

基于ASTER GED产品的地表发射率估算

地表发射率是地表温度反演的重要输入参数,为了解决现有地表发射率估算方法在裸露地表精度较差的问题,本文基于最新的ASTER全球地表发射率产品(ASTER GED)和基于植被覆盖度的方法(VCM),提出了一个改进的地表发射率估算方法。首先,利用ASTER GED产品求解裸土发射率,然后,利用ASTER波谱库中的植被发射率和植被覆盖度结合VCM方法计算地表发射率。利用张掖地区2012年11景ASTER TES算法反演的地表发射率产品和实测地表发射率数据进行了验证,同时利用一景Landsat 8 TIRS数据分析了对地表温度反演精度的影响。结果表明该方法估算的地表发射率整体精度较高,可以有效改进裸露地表的发射率估算精度,用于支持利用多种热红外传感器数据生产高精度的地表温度产品。     

Characterizing spring phenology of temperate broadleaf forests using Landsat and Sentinel-2 time series

Vegetation phenology has a great impact on land-atmosphere interactions like carbon cycling, albedo, and water and energy exchanges. To understand and predict these critical land-atmosphere feedbacks, it is crucial to measure and quantify phenological responses to climate variability, and ultimately climate change. Coarse-resolution sensors such as MODIS and AVHRR have been useful to study vegetation phenology from regional to global scales. These sensors are, however, not capable of discerning phenological variation at moderate spatial scales. By offering increased observation density and higher spatial resolution, the combination of Landsat and Sentinel-2 time series might provide the opportunity to overcome this limitation.In this study, we analyzed the potential of combined Sentinel-2 and Landsat time series for estimating start of season (SOS) of broadleaf forests across Germany for the year 2018. We tested two common statistical modeling approaches (logistic and generalized additive models using thin plate splines) and the two most commonly used vegetation indices, the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI).We found strong agreement between SOS estimates from logistic and spline models (r_EVI = 0.86; r_NDVI = 0.65), whereas agreement was higher for EVI than for NDVI (RMSD_EVI = 3.07, RMSD_NDVI = 5.26 days). The choice of vegetation index thus had a higher impact on the results than the fitting method. The EVI-based SOS also showed higher correlation with ground observations compared to NDVI (r_EVI = 0.51, r_NDVI = 0.42). Data density played an important role in estimating land surface phenology. Models combining Sentinel-2A/B, with an average cloud-free observation frequency of 12 days, were largely consistent with the combined Landsat and Sentinel-2 models, suggesting that Sentinel-2A/B may be sufficient to capture SOS for most areas in Germany in 2018. However, in non-overlapping swath areas and mountain areas, observation frequency was significantly lower, underlining the need to combine Landsat and Sentinel-2 for consistent SOS estimates over large areas. Our study demonstrates that estimating SOS of temperate broadleaf forests at medium spatial resolution has become feasible with combined Landsat and Sentinel-2 time series.     

热红外遥感地表温度与发射率地面验证进展

地表温度与发射率是地表—大气系统长波辐射和潜热通量交换的直接驱动力,是描述区域和全球尺度上地表能量平衡与水平衡的重要参数,其时空变化信息在气象预测、气候变化、水循环、地质勘探、农林监测和城市热环境等诸多领域具有广泛的应用。热红外遥感作为当前获取区域或全球尺度上地表温度和发射率的最有效手段之一,相较于传统的地面点位测量方法,具有空间覆盖范围大和重复观测等优势。对热红外遥感定量反演的地表温度与发射率产品进行地表真实性验证,有利于发现遥感数据自身或其反演算法的缺陷,确定产品的精度与不确定度,便于遥感产品的应用与推广。本文首先回顾了地表温度和发射率的定义,阐述了热红外遥感可反演、地面可测量的地表温度和发射率的科学内涵,并对利用热红外遥感数据反演地表温度和发射率的理论和方法作了概述;对地表温度和发射率地面验证的框架体系、验证指标进行总结,建立了基于精度、精确度、不确定度、完整性和稳定性的验证评价指标体系;总结了地表温度和地表发射率的地面验证方法、地面测量方法、辅助数据的获取方法、地表温度地面测量的采样方法,以及在验证异质非同温地表时从点到像元尺度的地表温度尺度转换方法等,分析了地面验证过程的主要误差来源;归纳了目前地表温度和地表发射率主要验证站点、观测网络及其空间分布特征;最后,本文讨论了地表温度与发射率地面验证存在的若干问题,并对地表温度与发射率验证工作的发展前景和趋势进行了相关展望。     

A sub-pixel analysis of urbanization effect on land surface temperature and its interplay with impervious surface and vegetation coverage in Indianapolis,United States

This study developed an analytical procedure based upon a spectral unmixing model for characterizing and quantifying urban landscape changes in Indianapolis, Indiana, the United States, and for examining the environmental impact of such changes on land surface temperatures (LST). Three dates of Landsat TM/ETM+ images, acquired in 1991, 1995, and 2000, respectively, were utilized to document the historical morphological changes in impervious surface and vegetation coverage and to analyze the relationship between these changes and those occurred in LST. Three fraction endmembers, i.e., impervious surface, green vegetation, and shade, were derived with an unconstrained least-squares solution. A hybrid classification procedure, which combined maximum-likelihood and decision-tree algorithms, was developed to classify the fraction images into land use and land cover classes. Correlation analyses were conducted to investigate the changing relationships of LST with impervious surface and vegetation coverage. Results indicate that multi-temporal fraction images were effective for quantifying the dynamics of urban morphology and for deriving a reliable measurement of environmental variables such as vegetation abundance and impervious surface coverage. Urbanization created an evolved inverse relationship between impervious and vegetation coverage, and brought about new LST patterns because of LST's correlations with both impervious and vegetation coverage. Further researches should be directed to refine spectral mixture modeling by stratification, and by the use of multiple endmembers and hyperspectral imagery.     

Retrieval of leaf water content spanning the visible to thermal infrared spectra

The objective of this study was to investigate the entire spectra (from visible to the thermal infrared; 0.390–14.0 μm) to retrieve leaf water content in a consistent manner. Narrow-band spectral indices (calculated from all possible two band combinations) and a partial least square regression (PLSR) were used to assess the strength of each spectral region. The coefficient of determination (R²) and root mean square error (RMSE) were used to report the prediction accuracy of spectral indices and PLSR models. In the visible-near infrared and shortwave infrared (VNIR–SWIR), the most accurate spectral index yielded R² of 0.89 and RMSE of 7.60%, whereas in the mid infrared (MIR) the highest R² was 0.93 and RMSE of 5.97%. Leaf water content was poorly predicted using two-band indices developed from the thermal infrared (R² = 0.33). The most accurate PLSR model resulted from MIR reflectance spectra (R² = 0.96, RMSE = 4.74% and RMSE cross validation RMSE_CV = 6.17%) followed by VNIR–SWIR reflectance spectra (R² = 0.91, RMSE = 6.90% and RMSE_CV = 7.32%). Using thermal infrared (TIR) spectra, the PLSR model yielded a moderate retrieval accuracy (R² = 0.67, RMSE = 13.27% and RMSE_CV = 16.39%). This study demonstrated that the mid infrared (MIR) and shortwave infrared (SWIR) domains were the most sensitive spectral region for the retrieval of leaf water content.     

遥测地表温度与比辐射率的迭代反演方法──理论推导与数值模拟

首先在地表比辐射率为已知的条件下,提出一个非线性迭代温度反演模型,我们对不同的地表和大气条件进行了模拟计算,结果表明当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±２０％时的温度均方根误差为０．４８Ｋ。当大气模式误差一个模式时反演的温度均方根误差为０．８５Ｋ。在此基础上,引人相邻像元的概念,相邻像元的大气状况可以认为是相同的,应用两个时相的遥感影像数据,假定在两个相近时相之间地表比辐射率值不变,建立地表比辐射率与温度的反演模型。我们对不同的地表和大气条件进行了模拟计算,结果表明当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±２０％时地表温度均方根误差小于１．５Ｋ,地表比辐射率均方根误差小于０．０２,地表辐射均方根误差为１％;当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±１０％时,地表温度均方根误差小于１．０Ｋ,地表辐射均方根误差小于０．６％。     

Parameterization of air temperature in high temporal and spatial resolution from a combination of the SEVIRI and MODIS instruments

Some applications, e.g. from traffic or energy management, require air temperature data in high spatial and temporal resolution at two metres height above the ground (T_2m), sometimes in near-real-time. Thus, a parameterization based on boundary layer physical principles was developed that determines the air temperature from remote sensing data (SEVIRI data aboard the MSG and MODIS data aboard Terra and Aqua satellites). The method consists of two parts. First, a downscaling procedure from the SEVIRI pixel resolution of several kilometres to a one kilometre spatial resolution is performed using a regression analysis between the land surface temperature (LST) and the normalized differential vegetation index (NDVI) acquired by the MODIS instrument. Second, the lapse rate between the LST and T_2m is removed using an empirical parameterization that requires albedo, down-welling surface short-wave flux, relief characteristics and NDVI data. The method was successfully tested for Slovenia, the French region Franche-Comté and southern Germany for the period from May to December 2005, indicating that the parameterization is valid for Central Europe. This parameterization results in a root mean square deviation RMSD of 2.0 K during the daytime with a bias of −0.01 K and a correlation coefficient of 0.95. This is promising, especially considering the high temporal (30 min) and spatial resolution (1000 m) of the results.     

光谱指数趋势面的城市地表温度降尺度转换

针对城市及周边区域建造区和自然地表交织分布的特点,探讨了利用归一化植被指数(NDVI)和归一化建造指数(NDBI)构造趋势面的地表温度(LST)降尺度方法,以北京市市区及周边较平坦区域为例实现了LST自960 m向120 m的降尺度转换。分析了LST空间分布特征及NDVI、NDBI对地物的指示性特征;以北京市四至六环为界分析NDVI、NDBI趋势面对地表温度的拟合程度及各自的适用区域;在120 m、240 m、480 m和960 m 4个尺度上评价了NDVI、NDBI和NDVI+NDBI趋势面对LST的拟合程度和趋势面转换函数的尺度效应;对NDVI、NDBI和NDVI NDBI等3种方法的降尺度结果分覆盖类型、分区域对比评价。实验结果表明结合两种光谱指数的NDVI NDBI方法降尺度转换精度有所改善,改善程度取决于地表覆盖类型组合。     

Quantification of LST and CO2 levels using Landsat-8 thermal bands on urban environment

Land surface temperature (LST) is an important aspect in global to regional change studies, for control of climate change and balancing of high temperature. Urbanization is one of the influencing factors increasing land surface and atmospheric temperature, by the emission of greenhouse gases (e.g. CO₂, NO and methane). In the present study, LST was derived from Landsat-8 of multitemporal data sets to analyse the spatial structure of the urban thermal environment in relation to the urban surface characteristics and land use–land cover (LULC). LST is influenced by the greenhouse gases i.e. CO₂ plays an important role in increasing the earth’s surface temperature. In order to provide the evidence of influence of CO₂ on LST, the relationship between LST, air temperature and CO₂ was analysed. Landsat-8 satellite has two thermal bands, 10 and 11. These bands were used to accurately to calculate the temperature over the study area. Results showed that the strength of correlation between ground monitoring data and satellite data was high. Based on correlation values of each month April (R² = 0.994), May (R² = 0.297) and June (R² = 0.934), observed results show that band 10 was significantly correlating with air temperature. Relationship between LST and CO₂ levels were obtained from linear regression analysis. band 11 was correlating significantly with CO₂ values in each of the months April (R² = 0.217), May (R² = 0.914) and June, (R² = 0.934), because band 11 is closer to the 15-micron band of CO₂. From the results, it was observed that band 10 can be used for calculating air temperature and band 11 can be used for estimation of greenhouse gases.     

Retrieving spatial variations of land surface temperatures from satellite data–Cairo region,Egypt

Cairo region is characterized by a range of physiographic features, including: flat agricultural lands, bare sandy deserts, highlands, calcareous terrains and urban land use. A time series data-set (300 images) acquired from the Moderate Resolution Imaging Spectroradiometer for the period July 2002–June 2015 were utilized to retrieve the spatial variations in the mean land surface temperature (LST) for the above-mentioned surface features. Results showed that vegetation, topography and surface albedo have negative correlations with LST. Vegetation/LST correlation has the maximum regression coefficient (R² = 0.68) and albedo/LST has the minimum (R² = 0.03). Cultivated lands reveal the lowest mean LST (<32 °C), whereas industrial lands exhibit the highest LST (>40 °C) of Cairo region. There is a considerable urban heat island formed at Helwan south of Cairo, where heavy industries are settled. Industrial activities raised the mean LST of the region by at least 4 °C than the surrounding urban lands.     

Soil respiration mapped by exclusively use of MODIS data for forest landscapes of Saskatchewan,Canada

Soil respiration (R_s) is of great importance to the global carbon balance. Remote sensing of R_s is challenging because of (1) the lack of long-term R_s data for model development and (2) limited knowledge of using satellite-based products to estimate R_s. Using 8-years (2002–2009) of continuous R_s measurements with nonsteady-state automated chamber systems at a Canadian boreal black spruce stand (SK-OBS), we found that R_s was strongly correlated with the product of the normalized difference vegetation index (NDVI) and the nighttime land surface temperature (LSTn) derived from Moderate Resolution Imaging Spectroradiometer (MODIS) imagery. The coefficients of the linear regression equation of this correlation between R_s and NDVI × LSTn could be further calibrated using the MODIS leaf area index (LAI) product, resulting in an algorithm that is driven solely by remote sensing observations. Modeled R_s closely tracked the seasonal patterns of measured R_s and explained 74–92% of the variance in R_s with a root mean square error (RMSE) less than 1.0 g C/m²/d. Further validation of the model from SK-OBS site at another two independent sites (SK-OA and SK-OJP, old aspen and old jack pine, respectively) showed that the algorithm can produce good estimates of R_s with an overall R² of 0.78 (p < 0.001) for data of these two sites. Consequently, we mapped R_s of forest landscapes of Saskatchewan using entirely MODIS observations for 2003 and spatial and temporal patterns of R_s were well modeled. These results point to a strong relationship between the soil respiratory process and canopy photosynthesis as indicated from the greenness index (i.e., NDVI), thereby implying the potential of remote sensing data for detecting variations in R_s. A combination of both biological and environmental variables estimated from remote sensing in this analysis may be valuable in future investigations of spatial and temporal characteristics of R_s.     

Land surface temperature as potential indicator of burn severity in forest Mediterranean ecosystems

Forest fires are one of the most important causes of environmental alteration in Mediterranean countries. Discrimination of different degrees of burn severity is critical for improving management of fire-affected areas. This paper aims to evaluate the usefulness of land surface temperature (LST) as potential indicator of burn severity. We used a large convention-dominated wildfire, which occurred on 19–21 September, 2012 in Northwestern Spain. From this area, a 1-year series of six LST images were generated from Landsat 7 Enhanced Thematic Mapper (ETM+) data using a single channel algorithm. Further, the Composite Burn Index (CBI) was measured in 111 field plots to identify the burn severity level (low, moderate, and high). Evaluation of the potential relationship between post-fire LST and ground measured CBI was performed by both correlation analysis and regression models. Correlation coefficients were higher in the immediate post-fire LST images, but decreased during the fall of 2012 and increased again with a second maximum value in summer, 2013. A linear regression model between post-fire LST and CBI allowed us to represent spatially predicted CBI (R-squared_adj > 85%). After performing an analysis of variance (ANOVA) between post-fire LST and CBI, a Fisher's least significant difference test determined that two burn severity levels (low-moderate and high) could be statistically distinguished. The identification of such burn severity levels is sufficient and useful to forest managers. We conclude that summer post-fire LST from moderate resolution satellite data may be considered as a valuable indicator of burn severity for large fires in Mediterranean forest ecosytems.     

Diurnal and seasonal impacts of urbanization on the urban thermal environment: A case study of Beijing using MODIS data

Beijing has experienced rapid urbanization and associated urban heat island effects and air pollution. In this study, a contribution index was proposed to explore the effect of urbanization on land surface temperature (LST) using Moderate-Resolution Imaging Spectroradiometer (MODIS)-derived data with high temporal resolution. The analysis indicated that different zones and landscapes make diurnally and seasonally different contributions to the regional thermal environment. The differences in contributions by the three main functional zones resulted from differences in their landscape compositions. The roles of landscapes in this process varied diurnally and seasonally. Urban land was the most important contributor to increases in regional LSTs. The contributions of cropland and forest varied distinctly between daytime and nighttime owing to differences in their thermal inertias. Vegetation had a notable cooling effect as the normalized vegetation difference index (NDVI) increased during summer. However, when the NDVI reached a certain value, the nighttime LST shifted markedly in other seasons. The results suggest that urban design based on vegetation partitions would be effective for regulating the thermal environment.     

Land surface temperature and its impact factors in Western Sichuan Plateau,China

The understanding influence of multiple factors variations on land surface temperature (LST) remains elusive. LST was retrieved by the atmospheric correction algorithms. Based on the correlation coefficients, stepwise regression analysis was developed to examine how multiple factors variability led to LST variations. The differences in LST between impact factors vary depending on time in a day. The elevation and land use types significantly affect the LST in sunny slope or shadow areas has a significantly quadratic curve correlation or a negative linear correlation with it, the influence of slope and aspect is not very significant. LST for forestland, grassland and bare land in the sunny slope and shadow area was the cubic polynomial related to its elevation. Normalized difference vegetation index (NDVI) and normalized difference moisture index (NDMI) effectively express LST in mountainous. LST and NDMI or NDVI have a significantly negative correlation, NDMI is more effective and more applicable for the expression of LST.     

Developing a temporally land cover-based look-up table (TL-LUT) method for estimating land surface temperature based on AMSR-E data over the Chinese landmass

The land surface temperature (LST) is an important parameter when studying the interface between the atmosphere and the Earth's surface. Compared to satellite thermal infrared (TIR) remote sensing, passive microwave (PMW) remote sensing is better able to overcome atmospheric influences and to estimate the LST, especially in cloudy regions. However, methods for estimating PMW LSTs at the country and continental scales are still rare. The necessity of training such methods from a temporally dynamic perspective also needs further investigations. Here, a temporally land cover based look-up table (TL-LUT) method is proposed to estimate the LSTs from AMSR-E data over the Chinese landmass. In this method, the synergies between observations from MODIS (Moderate Resolution Imaging Spectroradiometer) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS), which are onboard the same Aqua satellite, are explored. Validation with the synchronous MODIS LSTs demonstrates that the TL-LUT method has better performances in retrieving LSTs with AMSR-E data than the method that uses a single brightness temperature in 36.5 GHz vertical polarization channel. The accuracy of the TL-LUT method is better than 2.7 K for forest and 3.2 K for cropland. Its accuracy varies according to land cover type, time of day, and season. When compared with the in-situ measured LSTs at four sites without urban warming in the Tibet Plateau, the standard errors of estimation between the estimated AMSR-E LST and in-situ measured LST are from 5.1 K to 6.0 K in the daytime and 3.1 K to 4.5 K in the nighttime. Further comparison with the in-situ measured air temperatures at 24 meteorological stations confirms the good performance of the TL-LUT method. The feasibility of PMW remote sensing in estimating the LST for China can complement the TIR data and can, therefore, aid in the generation of daily LST maps for the entire country. Further study of the penetration of PMW radiation would benefit the LST estimations in barren and other sparsely vegetated environments.     

Retrieving aboveground biomass of wetland Phragmites australis (common reed) using a combination of airborne discrete-return LiDAR and hyperspectral data

Wetland biomass is essential for monitoring the stability and productivity of wetland ecosystems. Conventional field methods to measure or estimate wetland biomass are accurate and reliable, but expensive, time consuming and labor intensive. This research explored the potential for estimating wetland reed biomass using a combination of airborne discrete-return Light Detection and Ranging (LiDAR) and hyperspectral data. To derive the optimal predictor variables of reed biomass, a range of LiDAR and hyperspectral metrics at different spatial scales were regressed against the field-observed biomasses. The results showed that the LiDAR-derived H_p99 (99th percentile of the LiDAR height) and hyperspectral-calculated modified soil-adjusted vegetation index (MSAVI) were the best metrics for estimating reed biomass using the single regression model. Although the LiDAR data yielded a higher estimation accuracy compared to the hyperspectral data, the combination of LiDAR and hyperspectral data produced a more accurate prediction model for reed biomass (R² = 0.648, RMSE = 167.546 g/m², RMSE_r = 20.71%) than LiDAR data alone. Thus, combining LiDAR data with hyperspectral data has a great potential for improving the accuracy of aboveground biomass estimation.