Retrieval of land surface albedo and temperature using data from the Indian geostationary satellite: a case study for the winter months

The shortwave and longwave radiation budget at land surfaces is largely dependent on two fundamental quantities, the albedo and the land surface temperature (LST). A time series (November 2005 to March 2006) of daily data from the Indian geostationary satellite Kalpana‐1 Very High Resolution Radiometer (K1VHRR) sensor in the visible (VIS), water vapour (WV) and thermal infrared (TIR) bands from noontime (0900 GMT) observations were processed to retrieve these quantities in clear skies for five winter months. Cloud detection was carried out using bispectral threshold tests (in both VIS and TIR bands) in a dekadal time series. Surface albedo was retrieved using a simple atmospheric transmission model. K1VHRR albedo was compared with Moderate Resolution Imaging Spectroradiometer (MODIS) AQUA noontime albedo over different land targets (agriculture, forest, desert, scrub and snow) that showed minimum differences over agriculture and forest. The comparison of spatial albedo over different landscapes yielded a root mean square deviation (RMSD) of 0.021 in VHRR albedo (9% of MODIS albedo). A mono‐window algorithm was implemented with a single TIR band to retrieve the LST. Its accuracy was also verified over different land targets by comparison with aggregated MODIS AQUA LST. The maximum RMSD was obtained over agriculture. Spatial comparison of VHRR and AQUA LSTs over homogeneous and heterogeneous landscape cutouts revealed an overall RMSD of 2.3 K. An improvement in the retrieval accuracy is expected to be achieved with atmospheric products from the sounder and split thermal bands in the imager of future INSAT 3D missions.     

Preliminary applications of a land surface temperature retrieval method to IASI and AIRS data

Land surface temperature (LST) is one of the key state variables for many applications. This article aims to apply our previously developed LST retrieval method to infrared atmospheric sounding interferometer (IASI) and atmospheric infrared sounder (AIRS) data. On the basis of the opposite characteristics of the atmospheric spectral absorption and surface spectral emissivity, a ‘downwelling radiance residual index’ (DRRI) has been recalled and improved to obtain LST and emissivity. To construct an efficient DRRI, an automatic channel selection procedure has been proposed, and 11 groups of channels have been selected within the range 800–1000 cm^?1. The DRRI has been tested with IASI and AIRS data. For the IASI data, the radiosonde data have been used to correct for atmospheric effects and to retrieve LST, while the atmospheric profiles retrieved from AIRS data have been used to perform the atmospheric corrections and subsequently to estimate LST from AIRS data. The differences between IASI- and Moderate Resolution Imaging Spectroradiometer (MODIS)-derived LSTs are no more than 2 K, while the differences between AIRS- and MODIS-derived LSTs are less than 5 K. Even though an exceptionally problematic value occurred (–12.89 K), the overall differences between AIRS-estimated LST and the AIRS L2 LST product are no more than 5 K. Although the IASI-derived LST is more accurate than the AIRS-derived one, the convenient retrieval of AIRS atmospheric profile made this method more applicable. Limitations and uncertainties in retrieving LST using the DRRI method are also discussed.     

Derivation of land surface temperatures from MODIS data using the general split‐window technique

Fast Atmospheric Signature Code (FASCODE), a line‐by‐line radiative transfer programme, was used to simulate Moderate Resolution Imaging Spectroradiometer (MODIS) data at wavelengths 11.03 and 12.02 µm to ascertain how accurately the land surface temperature (LST) can be inferred, by the split‐window technique (SWT), for a wide range of atmospheric and terrestrial conditions. The approach starts from the Ulivieri algorithm, originally applied to Advanced Very High Resolution Radiometer (AVHRR) channels 4 and 5. This algorithm proved to be very accurate compared to several others and takes into account the atmospheric effects, in particular the water vapour column (WVC) amount and a non‐unitary surface emissivity. Extended simulations allowed the determination of new coefficients of this algorithm appropriate to MODIS bands 31 and 32, using different atmospheric conditions. The algorithm was also improved by removing some of the hypothesis on which its original expression was based. This led to the addition of a new corrective term that took into account the interdependence between water vapour and non‐unitary emissivity values and their effects on the retrieved surface temperature. The LST products were validated within 1 K with in situ LSTs in 11 cases.     

A C++ program for retrieving land surface temperature from the data of Landsat TM/ETM+ band6

A C++ language-based software tool for retrieving land surface temperature (LST) from the data of Landsat TM/ETM+ band6 is developed. It has two main functional modules: (1) Three methods to compute the ground emissivity based on land use/cover classification image, NDVI image and the ratio values of vegetation and bare ground and (2) Converting digital numbers (DNs) from TM/ETM+ band6 to LST. In the software tool, Qin et al.'s mono-window algorithm and Jiménez-Muňoz and Sobrino's single channel algorithm are programmed to retrieve LST. It will be a useful software tool to study the thermal environment of ground surface or the energy balance between the ground and the bottom atmosphere by using the thermal band of Landsat TM/ETM+.     

Validation of land surface temperature derived from MSG/SEVIRI with in situ measurements at Gobabeb,Namibia

Land surface temperature (LST) derived from Meteosat Second Generation/?Spinning-Enhanced Visible and Infrared Imager MSG/SEVIRI data is an operational product of the Land Surface Analysis Satellite Applications Facility (LSA SAF). The LST has a temporal resolution of 15 minutes, a sampling distance of 3 km at nadir, and a targeted accuracy of better than 2 K. Gobabeb (Namibia) is one of Karlsruhe Institute of Technology's (KIT's) four dedicated stations for LST validation. In March 2010, a field survey was performed to characterize the Gobabeb site more closely. SAF LST and in situ LST obtained over a period of 3 days from additional measurements with a telescopic mast on the Namib gravel plains were in good agreement with each other (bias 1.0 K). For the same period, the bias between SAF LST and Gobabeb main station LST was even smaller (0.4 K). A mobile measurement system was set up by fixing the telescopic mast to a four-wheel drive. Around solar noon, LST from in situ measurements along a 40 km track and LST from Gobabeb main station had a bias of 0.4 K and a standard deviation of 1.2 K, which means that in situ LSTs at Gobabeb main station are representative for large parts of the gravel plains. Exploiting this relationship, 2 years of LST from MSG/SEVIRI were compared with in situ LST from Gobabeb main station. The magnitude of the monthly biases between the two data sets was generally less than 1.0 K and root mean square errors were below 1.5 K. Furthermore, the bias appears to exhibit a seasonality, which could be accounted for in future validation work.     

Surface temperature and water vapour retrieval from MODIS data

This paper gives operational algorithms for retrieving sea (SST), land surface temperature (LST) and total atmospheric water vapour content (W) using Moderate Resolution Imaging Spectroradiometer (MODIS) data. To this end, the MODTRAN 3.5 radiative transfer program was used to predict radiances for MODIS channels 31, 32, 2, 17, 18 and 19. To analyse atmospheric effects, a simulation with a set of radiosonde observations was used to cover the variability of surface temperature and water vapour concentration on a worldwide scale. These simulated data were split into two sets (DB1 and DB2), the first one (DB1) was used to fit the coefficients of the algorithms, while the second one (DB2) was used to test the fitted coefficients. The results show that the algorithms are capable of producing SST and LST with a standard deviation of 0.3 K and 0.7 K if the satellite data are error free. The LST product has been validated with in situ data from a field campaign carried out in the Mississippi (USA), the results show for the LST algorithm proposed a root mean square error lower than 0.5K. Regarding water vapour content, a ratio technique is proposed, which is capable of estimating W from the absorbing channels at 0.905, 0.936, and 0.94,µm, and the atmospheric window channel at 0.865,µm, with a standard deviation (in the comparison with radiosonde observations) of 0.4 g cm^?2.     

Land surface temperature retrieval from LANDSAT TM 5

In this paper, three methods to retrieve the land surface temperature (LST) from thermal infrared data supplied by band 6 of the Thematic Mapper (TM) sensor onboard the Landsat 5 satellite are compared. The first of them lies on the estimation of the land surface temperature from the radiative transfer equation using in situ radiosounding data. The others two are the mono-window algorithm developed by Qin et al. [International Journal of Remote Sensing 22 (2001) 3719] and the single-channel algorithm developed by Jiménez-Muñoz and Sobrino [Journal of Geophysical Research 108 (2003)]. The land surface emissivity (LSE) values needed in order to apply these methods have been estimated from a methodology that uses the visible and near infrared bands. Finally, we present a comparison between the LST measured in situ and the retrieved by the algorithms over an agricultural region of Spain (La Plana de Requena-Utiel). The results show a root mean square deviation (rmsd) of 0.009 for emissivity and lower than 1 K for land surface temperature when the Jiménez-Muñoz algorithm is used.     

Accuracy assessment of land surface temperature retrievals from MSG2-SEVIRI data

The accuracy of the Land Surface Temperature (LST) product generated operationally by the EUMETSAT Land Surface Analysis Satellite Applications Facility (LSA SAF) from the data registered by the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on board the geostationary METEOSAT Second Generation 2 (MSG2, Meteosat 9) satellite was assessed on two test sites in Eastern Spain: a homogeneous, fully vegetated rice field and a high-plain, homogeneous area of shrubland. The LSA SAF LSTs were compared with ground LST measurements in the conventional temperature-based (T-based) method. We also validated the LSA SAF LST product by using an alternative radiance-based (R-based) method, with ground LSTs calculated from MSG-SEVIRI channel 9 brightness temperatures (at 10.8 μm) through radiative transfer simulations using atmospheric temperature and water vapor profiles together with surface emissivity data. Two lakes were also used for validation with the R-based method. Although the LSA SAF LST algorithm works mostly within the uncertainty expectation of ± 2 K, both validation methods showed significant biases for the LSA SAF LST product, up to 1.5 K in some cases. These biases, with the LSA SAF LST product overestimating reference values, were also observed in previous studies. Nevertheless, the present work points out that the biases are related to the land surface emissivities used in the operational generation of the product. The use of more appropriate emissivity values for the test sites in the LSA SAF LST algorithm led to better results by decreasing the biases by 0.7 K for the shrubland validation site. Furthermore, we proposed and checked an alternative algorithm: a quadratic split-window equation, based on a physical split-window model that has been widely proved for other sensors, with angular-dependent coefficients suitable for the MSG coverage area. The T-based validation results for this algorithm showed LST uncertainties (robust root-mean-squared-errors) from 0.2 K to 0.5 K lower than for the LSA SAF LST algorithm after the emissivity replacement. Nevertheless, the proposed algorithm accuracies were significantly better than those obtained for the current LSA SAF LST product, with an average accuracy difference of 0.6 K.     

FY3A/MERSI地表温度反演

MERSI是我国第二代极轨气象卫星上的重要传感器,可获取高空间分辨率和高时间分辨率的对地观测影像。为使Jimènez-Mu珘nozSobrino算法更适用于FY3A/MERSI传感器通道特性,更新了大气函数的估算系数,并引入观测角度因子,以获取更为精确像元间更为平滑的地表温度。用MODTRAN4模拟验证该算法精度,得引入角度因子后反演精度显著提升,所有角度下平均误差为-0.6±2.2K。用实测的敦煌戈壁地表温度和MODIS地表温度产品评价MERSI反演结果,显示MERSI地表温度的空间分布准确,结果精度也较高。与实测温度对比,平均误差为1.74K,均方根误差小于1.9K。研究区域与MODIS地表温度间差异平均为2.6307K。虽然会受云检测精度和观测亮温偏高的影响,由MERSI反演的高精度地表温度在相关科研和业务方面仍然具有极好的应用前景。     

Evaluation of 10 year AQUA/MODIS land surface temperature with SURFRAD observations

As the 10 year Moderate Resolution Imaging Spectroradiometer Land Surface Temperature MODIS LST becomes available, it is significant to perform a comprehensive evaluation on the long-term product before downstream users use it for climate studies and atmospheric models. In this study, a validation is carried out using observations from the US Surface Radiation budget (SURFRAD) network. Strict quality control removes cloud-contaminated samples from MODIS LST collection and decreases noise information from SURFRAD measurements, thereby making the validation more persuasive. With analysis on 19,735 valid samples, Aqua/MODIS LST from a split-window algorithm shows retrieval errors from –14 K to 17 K with a bias of –0.93 K, an RMSE of 2.65 K, and a standard deviation of 2.48 K. The errors also show strong seasonal signals. With correlation tests between LST errors and several other factors, it is disclosed that LST retrieval errors mainly come from atmospheric effects and surface emissivity uncertainties, which are closely related to relative air humidity, absolute air humidity, sensor zenith angle, wind speed, normalized difference vegetation index (NDVI), and soil moisture. In addition, the impacts from these factors may not be independent. These impact factors suggest a deficiency of the split-window algorithm in dealing with atmospheric and surface complexity and variety.     

A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region

Remote sensing of land surface temperature (LST) from the thermal band data of Landsat Thematic Mapper (TM) still remains unused in comparison with the extensive studies of its visible and near-infrared (NIR) bands for various applications. The brightness temperature can be computed from the digital number (DN) of TM6 data using the equation provided by the National Aeronautics and Space Administration (NASA). However, a proper algorithm for retrieving LST from the only one thermal band of the sensor still remains unavailable due to many difficulties in the atmospheric correction. Based on thermal radiance transfer equation, an attempt has been made in the paper to develop a mono-window algorithm for retrieving LST from Landsat TM6 data. Three parameters are required for the algorithm: emissivity, transmittance and effective mean atmospheric temperature. Method about determination of atmospheric transmittance is given in the paper through the simulation of atmospheric conditions with LOWTRAN 7 program. A practicable approach of estimating effective mean atmospheric temperature from local meteorological observation is also proposed in the paper when the in situ atmospheric profile data is unavailable at the satellite pass, which is generally the case in the real world especially for the images in the past. Sensitivity analysis of the algorithm indicates that the possible error of ground emissivity, which is difficult to estimate, has relatively insignificant impact on the probable LST estimation error i T, which is sensible to the possible error of transmittance i 6 and mean atmospheric temperature i T a . Validation of the simulated data for various situations of seven typical atmospheres indicates that the algorithm is able to provide an accurate LST retrieval from TM6 data. The LST difference between the retrieved and the simulated ones is less than 0.4°C for most situations. Application of the algorithm to the sand dunes across the Israel-Egypt border results in a reasonable LST estimation of the region. Based on this LST estimation, spatial variation of the interesting thermal phenomenon has been analysed for comparison of LST difference across the border. The result shows that the Israeli side does have significantly higher surface temperature in spite of its denser vegetation cover than the Egyptian side where bare sand is prevalent.     

Intercomparison of methods for estimating land surface temperature from a Landsat-5 TM image in an arid region with low water vapour in the atmosphere

This study compares the methods for retrieving the land surface temperature (LST) (T _s) from Landsat-5 TM (Thematic Mapper) data, including the radiative transfer equation (RTE) method, the mono-window algorithm (MWA) and the generalized single-channel (GSC) method in an arid region with low atmospheric water vapour content. In addition, T _s calculated without atmospheric correction of TM band 6 is also assessed. The intercomparison is divided into two parts. The first part is applying the methods at the Biandukou site (100° 58′ E, 38° 16′ N, elevation?=?2690 m) and the second part is applying them at Binggou (100° 13′ E, 38° 42′ N, elevation?=?3400 m) and Arou (100° 27′ E, 38° 36′ N, elevation?=?2960 m) sites. Results demonstrate that these methods provide acceptable accuracies at the Biandukou site. At this site, GSC generates nearly the same accuracy as RTE; MWA estimations are slightly less accurate than RTE and GSC; estimations without atmospheric correction of TM band 6 exhibit the largest errors. On the other hand, MWA is a good choice for retrieving the LST at Binggou and Arou sites. In cases where the meteorological parameters are unavailable, it is an alternative option to calculate T _s directly from TM band 6 image without atmospheric correction at these two sites.     

Developed split window method to retrieve land surface temperature from HYTES thermal hyperspectral images using genetic algorithm

This article aims to establish a new method to retrieve land surface temperature (LST) from hyperspectral thermal emission spectrometer (HYTES) data with split window (SW) algorithm. First, the optimal bands of HYTES sensor were selected with the genetic algorithm and then were used in the SW algorithm. In the SW algorithm, its coefficients were obtained based on several subranges of atmospheric column water vapours (CWVs) and view zenith angle (VZA) under various land surface conditions, in order to remove the atmospheric effect and improve the retrieval accuracy. Results showed that the root-mean-square error (RMSE) varies for different CWV and VZA, and with the increasing CWV and VZA, the RMSE value also increases. The emissivity, CWV, and VZA were also obtained for pixels. The sensitive analysis of LST retrieval to instrument noise and uncertainty of pixel emissivity and water vapour demonstrated the good performance of the proposed algorithm. Finally, the new algorithm was applied to HYTES sensor data, and the LST was validated using LST product of HYTES sensor obtained by NASA. The results showed that the RMSE of the LST retrieval with the proposed algorithm and the LST product of sensor for data 1 and data 2 is 1.3 K and 1.6 K, respectively.     

基于Landsat-5 TM的洪河湿地地表温度估算方法对比研究

利用洪河湿地2008年5月15日过境的Landsat/TM图像和实测地面数据以及MODIS 地表发射率数据,分别运用大气辐射传输模型、覃志豪的单窗算法和Jimenez\|Munoz & Sobrino 的单波段算法估算洪河湿地的地表温度,并且对比了大气校正前后的NDVI、LSE以及各种算法估算地表温度的差异。分析估算结果表明,覃志豪的单窗算法与实测地面数据估算结果非常一致。指出在没有实时探空数据的情况下,应用只有一个热红外通道的Landsat/TM数据源,采用覃志豪的单窗算法估算的精度是可以接受的。     

地表温度遥感中大气平均作用温度估算模型研究

大气平均作用温度Ta是地表温度遥感单窗算法中一个关键的参数,利用2008~2011年全国123个探空站点资料,针对大气水汽量的垂直分布特征,分析了利用近地层气温T0估算大气有效平均温度的可行性;进一步分析了T0和Ta之间的相关性,建立了适合我国地区大气平均温度估算的最佳模型Ta=44.97098+0.80512 T0,模型的决定系数R²为0.859,均方根误差为4.198 K。通过对44幅HJ\|1B/IRS热红外图像地温反演的敏感性分析,结果表明:模型估算的Ta用于地表温度反演时的误差为1.734 K;当大气透射率τ很小时,模型估算的Ta误差对地温反演很敏感,较小的估算误差会给地温反演带来很大的误差;随着大气透射率τ的增加,Ta的估算误差对地温反演的敏感性逐渐降低。
     

A three-channel algorithm for retrieving night-time land surface temperature from MODIS data under thin cirrus clouds

Thin cirrus clouds are dominated by non-spherical ice crystals with an effective emissivity of less than 0.5. Until now, the influences of clouds were not commonly considered in the development of algorithms for retrieving land-surface temperature (LST). However, numerical simulations showed that the influence of thin cirrus clouds could lead to a maximum LST retrieval error of more than 14 K at night if the cirrus optical depth (COD) at 12 μm was equal to 0.7 (cirrus emissivity equivalent to 0.5). To obtain an accurate estimate of the LST under thin cirrus using satellite infrared data, a nonlinear three-channel LST retrieval algorithm was proposed based on a widely used two-channel algorithm for clear-sky conditions. The variations in the cloud top height, COD, and effective radius of cirrus clouds were considered in this three-channel LST retrieval algorithm. Using Moderate Resolution Imaging Spectroradiometer (MODIS) channels 20, 31, and 32 (centred at 3.8, 11.0, and 12.0 μm, respectively) and the corresponding land surface emissivities (LSEs), the simulated data showed that this algorithm could obtain LSTs with root mean square errors (RMSEs) of less than 2.8 K when the COD at 12 μm is less than 0.7 and the viewing zenith angle (VZA) is less than 60°. In addition, a sensitivity analysis of the proposed algorithm showed that the total LST errors, including errors from the uncertainties in input parameters and algorithm error, were nearly the same as the algorithm error itself. Some lake surface water temperatures measured in Lake Superior and Lake Erie were used to test the performance of the proposed LST retrieval algorithm. The results showed that the proposed nonlinear three-channel algorithm could be used for estimating LST under thin cirrus with an RMSE of less than 2.8 K.     

环境一号B星热红外波段单通道算法温度反演

文中在考虑环境一号B星(HJ-1B)热红外波段(infrared scanner,IRS4)光谱响应函数和有效波长的基础上,通过MODTRAN4模型模拟,对Jimenez-Munoz和Sobrino(JM&S)单通道算法中的大气函数进行改进,重新计算得到了适合HJ-1B星IRS4地表温度(land surface temperature,LST)反演的3个大气函数公式,并反演了福州地区的地表温度.采用基于星上辐亮度法对反演的地表温度进行精度评价,并将反演的地表温度与JM&S算法、段四波等修正的JM&S算法反演的地表温度进行对比分析.结果表明:使用文中改进后的大气参数对HJ-1B星IRS4进行地表温度反演,可取得较好结果.     

利用TM6数据反演陆地表面温度新算法研究

陆地表面温度(LST)反演一直是热红外遥感研究中的一大难题。虽然TM 6数据具有较高的空间分辨率(120 m),但由于只有一个热通道,要得到地表真实温度,原来需要利用辐射传输方程的方法,实时资料的缺乏限制了该方法的应用。因而由TM 6数据得到的通常都是星上亮度温度,而星上亮度温度与实际地表温度差距较大,因此,其反演的温度精度不高。而单窗算法和普适性单通道算法的提出为从TM 6数据较高精度地反演陆地表面温度提供了可能。分析和研究了这两个新的单通道温度反演算法,并针对北京市的实际情况,利用2005年5月6日的TM数据对北京市的陆地表面温度进行了反演,并用实地测量数据进行了比较验证。结果表明这两种温度反演算法都取得了较高的精度,它们的rm sd值分别为1.38°和2.18°。     

基于IDL的遥感影像大气与地形校正方法实现

介绍了利用交互式数据语言（Interactive Data Language,IDL）开发TM/ETM遥感影像大气与地形校正模型的详细过程,以2000年4月30日密云ETM影像为例,对大气与地形校正方法的有效性和实用性进行了验证。结果表明,该方法有效地消除了大气与地形影响,提高了地表反射率等地表参数的反演精度和数据质量,为进一步开展定量遥感研究提供了数据质量保障。     

Comparison of satellite-derived and in-situ observations of ice and snow surface temperatures over Greenland

The most practical way to get spatially broad and continuous measurements of the surface temperature in the data-sparse cryosphere is by satellite remote sensing. The uncertainties in satellite-derived LSTs must be understood to develop internally-consistent decade-scale land surface temperature (LST) records needed for climate studies. In this work we assess satellite-derived “clear-sky” LST products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), and LSTs derived from the Enhanced Thematic Mapper Plus (ETM+) over snow and ice on Greenland. When possible, we compare satellite-derived LSTs with in-situ air temperature observations from Greenland Climate Network (GC-Net) automatic weather stations (AWS). We find that MODIS, ASTER and ETM+ provide reliable and consistent LSTs under clear-sky conditions and relatively-flat terrain over snow and ice targets over a range of temperatures from ? 40 to 0 °C. The satellite-derived LSTs agree within a relative RMS uncertainty of ~ 0.5 °C. The good agreement among the LSTs derived from the various satellite instruments is especially notable since different spectral channels and different retrieval algorithms are used to calculate LST from the raw satellite data. The AWS record in-situ data at a “point” while the satellite instruments record data over an area varying in size from: 57 × 57 m (ETM+), 90 × 90 m (ASTER), or to 1 × 1 km (MODIS). Surface topography and other factors contribute to variability of LST within a pixel, thus the AWS measurements may not be representative of the LST of the pixel. Without more information on the local spatial patterns of LST, the AWS LST cannot be considered valid ground truth for the satellite measurements, with RMS uncertainty ~ 2 °C. Despite the relatively large AWS-derived uncertainty, we find LST data are characterized by high accuracy but have uncertain absolute precision.