Land surface temperature and emissivity retrieval by integrating MODIS data onboard Terra and Aqua satellites

This paper presents an algorithm to retrieve land surface temperature (LST) and emissivity by integrating MODIS (Moderate Resolution Imaging Spectroradiometer) data onboard Terra and Aqua satellites. For a study area, there will be four pairs of day and night observations by MODIS onboard two satellites every day. Solar zenith angle, view zenith angle, and atmospheric water vapour have first been taken as independent variables to analyse their sensitivities to the same infrared channel measurements of MODIS on both Terra and Aqua satellites. Owing to their similar influences on the same MODIS band from Terra and Aqua satellites, four pairs of MODIS data from Terra and Aqua satellites can be thought of as MODIS measurement on a satellite at different viewing angles and viewing time. Comparisons between the retrieved results and in-situ measurements at three test sites (Qinghai Lake, Poyang Lake and Luancheng in China) indicate that the root mean square (rms) error is 0.66 K, except for the sand in Poyang Lake area. The rms error is less than 0.7 K when the retrieved results are compared with Earth Observing System (EOS) MODIS LST data products using the physics-based day/night algorithm. Emissivities retrieved by this algorithm are well compared to EOS MODIS emissivity data products (V5). The proposed algorithm can therefore be regarded as complementary and an extension to the EOS physics-based day/night algorithm.     

Evaluating NDVI-based emissivities of MODIS bands 31 and 32 using emissivities derived by Day/Night LST algorithm

Surface emissivity estimation is a significant factor for the land surface temperature estimation from remotely sensed data. For fully vegetated surfaces, the emissivity estimation is performed in a simple manner since the emissivity is relatively uniform. However, for arid land with sparse vegetation, the estimation is more complicated since the emissivity of the exposed soil and rock is highly variable. In this study, mean and difference emissivity for bands 31 and 32 of MODIS sensor have been derived based on NDVI values. First, the NDVI thresholds have been determined to separate bare soil, partially vegetated soil and fully vegetated land. Then regression relations have been derived to estimate mean and difference emissivity of the bare soil samples and partially vegetated surfaces. A constant emissivity is also used for fully vegetated area. Along with the correlations, standard deviations of the regression relations have been examined for a set of representative soil types. Standard deviations smaller than 0.003 in mean emissivity and smaller than 0.004 in difference emissivity are resulted in regression linear relations. Evaluation of the NDVI derived regression relations has been performed using the results of MODIS Day/Night Land Surface Temperature (LST) algorithm on a pair of MODIS images. Using around 45,500 pixels with different soil and land cover types, emissivity of each pixel in bands 31 and 32 have been estimated. The calculated emissivities have been compared with emissivities calculated by MODIS Day/Night LST algorithm. Biases and standard deviations of NDVI-based relations show relatively high agreement for mean and difference emissivity relations with Day/Night method results. It may be concluded that the proposed algorithm can be used as a rather simple alternative to complex emissivity estimation algorithms.     

中国江淮、黄淮地区陆面微波比辐射率的变化特征

陆面微波比辐射率较高且易变,造成陆面上反演降水以及其它大气参数较为困难。对于地表特征复杂的中国,陆面微波比辐射率的研究还很有限。通过利用Tropical Rainfall Measuring Mission (TRMM)卫星上同步扫描的VIRS（红外和可见光）与TMI（微波）资料以及微波辐射传输模式反演了中国江淮、黄淮地区陆面微波比辐射率。然后,结合MODIS提供的地表类型数据,分析了江淮、黄淮地区不同地表微波比辐射率的时空变化特征。 结果表明该地区的农作物地表比辐射率最小,垂直与水平比辐射率极化差最大;而森林地表比辐射率最大,极化差最小。此外,不同地表的微波比辐射率昼夜变化明显,季节变化不明显。比辐射率估算误差中,地表温度、微波亮温和大气相对湿度3因子的准确计算对22 GHz和85 GHz的影响较为明显,对其它通道影响较小。对于小于85 GHz的通道,比辐射率估算精度受微波亮温的影响最为明显,地表温度其次,相对湿度最小;对于高频85 GHz,相对湿度的影响最明显,其次是微波亮温,最后是地表温度。     

针对MODIS 数据的地表温度非线性迭代反演方法

地表温度是气象、水文、生态等研究领域中的一个重要参数。构建了MODIS31/ 32 波段的热辐射传输方程, 讨论了方程的数值迭代解法, 提出了针对MODIS 数据地表温度的非线性迭代反演方法, 并介绍了大气透过率和地表比辐射率这两个中间参数的估计方法。误差及敏感性分析表明,提出的方法对大气透过率和地表比辐射率都不敏感, 反演精度优于传统的线性分裂窗算法。     

Ground measurements for the validation of land surface temperatures derived from AATSR and MODIS data

An experimental site was set up in a large, flat and homogeneous area of rice crops for the validation of satellite derived land surface temperature (LST). Experimental campaigns were held in the summers of 2002-2004, when rice crops show full vegetation cover. LSTs were measured radiometrically along transects covering an area of 1 km². A total number of four thermal radiometers were used, which were calibrated and inter-compared through the campaigns. Radiometric temperatures were corrected for emissivity effects using field emissivity and downwelling sky radiance measurements. A database of ground-based LSTs corresponding to morning, cloud-free overpasses of Envisat/Advanced Along-Track Scanning Radiometer (AATSR) and Terra/Moderate Resolution Imaging Spectroradiometer (MODIS) is presented. Ground LSTs ranged from 25 to 32 °C, with uncertainties between ± 0.5 and ± 0.9 °C. The largest part of these uncertainties was due to the spatial variability of surface temperature. The database was used for the validation of LSTs derived from the operational AATSR and MODIS split-window algorithms, which are currently used to generate the LST product in the L2 level data. A quadratic, emissivity dependent split-window equation applicable to both AATSR and MODIS data was checked as well. Although the number of cases analyzed is limited (five concurrences for AATSR and eleven for MODIS), it can be concluded that the split-window algorithms work well, provided that the characteristics of the area are adequately prescribed, either through the classification of the land cover type and the vegetation cover, or with the surface emissivity. In this case, the AATSR LSTs yielded an average error or bias of − 0.9 °C (ground minus algorithm), with a standard deviation of 0.9 °C. The MODIS LST product agreed well with the ground LSTs, with differences comparable or smaller than the uncertainties of the ground measurements for most of the days (bias of + 0.1 °C and standard deviation of 0.6 °C, for cloud-free cases and viewing angles smaller than 60°). The quadratic split-window algorithm resulted in small average errors (+ 0.3 °C for AATSR and 0.0 °C for MODIS), with differences not exceeding ± 1.0 °C for most of the days (standard deviation of 0.9 °C for AATSR and 0.5 °C for MODIS).     

Evaluation and improvement of the MODIS land surface temperature/emissivity products using ground‐based measurements at a semi‐desert site on the western Tibetan Plateau

Current MODerate‐resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST, surface skin temperature)/emissivity products are evaluated and improvements are investigated. The ground‐based measurements of LST at Gaize (32.30° N, 84.06° E, 4420 m) on the western Tibetan Plateau from January 2001 to December 2002 agree well (mean and standard deviation of differences of 0.27 K and 0.84 K) with the 1‐km Version 004 (V4) Terra MODIS LST product (MOD11A1) generated by the split‐window algorithm. Spectral emissivities measured from surface soil samples collected at and around the Gaize site are in close agreement with the landcover‐based emissivities in bands 31 and 32 used by the split‐window algorithm. The LSTs in the V4 MODIS LST/emissivity products (MYD11B1 for Aqua and MOD11B1 for Terra) from the day/night LST algorithm are higher by 1–1.7 K (standard deviation around 0.6 K) in comparisons to the 5‐km grid aggregated values of the LSTs in the 1‐km products, which is consistent with the results of a comparison of emissivities. On average, the emissivity in MYD11B1 (MOD11B1) is 0.0107 (0.0167) less than the ground‐based measurements, which is equivalent to a 0.64 K (1.25 K) overestimation of LST around the average value of 285 K. Knowledge obtained from the evaluation of MODIS LST/emissivity retrievals provides useful information for the improvement of the MODIS LST day/night algorithm. Improved performance of the refined (V5) day/night algorithm was demonstrated with the Terra MODIS data in May–June 2004.     

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 generalized split-window algorithm for land surface temperature estimation from MSG-2/SEVIRI data

This paper aims to determine land surface temperature (LST) using data from a spinning enhanced visible and infrared imager (SEVIRI) on board Meteosat Second Generation 2 (MSG-2) by using the generalized split-window (GSW) algorithm. Coefficients in the GSW algorithm are pre-determined for several overlapping sub-ranges of the LST, land surface emissivity (LSE), and atmospheric water vapour content (WVC) using the data simulated with the atmospheric radiative transfer model MODTRAN 4.0 under various surface and atmospheric conditions for 11 view zenith angles (VZAs) ranging from 0° to 67°. The results show that the root mean square error (RMSE) varies with VZA and atmospheric WVC and that the RMSEs are within 1.0 K for the sub-ranges in which the VZA is less than 30° and the atmospheric WVC is less than 4.25 g cm^?2. A sensitivity analysis of LSE uncertainty, atmospheric WVC uncertainty, and instrumental noise (NEΔT) is also performed, and the results demonstrate that LSE uncertainty can result in a larger LST error than other uncertainties and that the total error for the LST is approximately 1.21 and 1.45 K for dry atmosphere and 0.86 and 2.91 K for wet atmosphere at VZA = 0° and at VZA = 67°, respectively, if the uncertainty in the LSE is 1% and that in the WVC is 20%. The GSW algorithm is then applied to the MSG-2 – SEVIRI data with the LSE determined using the temperature-independent spectral indices method and the WVC either determined using the measurements in two split-window channels or interpolated temporally and spatially using European Centre for Medium Range Weather Forecasting (ECMWF) data. Finally, the SEVIRI LST derived in this paper (SEVIRI LST1) is evaluated through comparisons with the SEVIRI LST provided by the land surface analysis satellite applications facility (LSA SAF) (SEVIRI LST2) and the Moderate Resolution Imaging Spectroradiometer (MODIS) LST product (MOD11B1 LST product). The results show that more than 80% of the differences between SEVIRI LST1 and SEVIRI LST2 are within 2 K, and approximately 70% of the differences between SEVIRI LST1 and MODIS LST are within 4 K. Furthermore, compared to MODIS LST, for four specific areas with different land surfaces, our GSW algorithm overestimates the LST by up to 1.0 K for vegetated surfaces and by 1.3 K for bare soil.     

Quality assessment and validation of the MODIS global land surface temperature

This paper presents an evaluation of the Earth Observing System (EOS) Moderate Resolution Imaging Spectroradiometer (MODIS) thermal infrared bands and the status of land surface temperature (LST) version-3 standard products retrieved from Terra MODIS data. The accuracy of daily MODIS LST products has been validated in more than 20 clear-sky cases with in situ measurement data collected in field campaigns in 2000–2002. The MODIS LST accuracy is better than 1°C in the range from ?10 to 50°C. Refinements and improvements were made to the new version of MODIS LST product generation executive code. Using both Terra and Aqua MODIS data for LST retrieval improves the quality of the LST product and the diurnal feature in the product due to better temporal, spatial and angular coverage of clear-sky observations.     

Evaluation of land surface temperature and emissivities retrieved from MSG/SEVIRI data with MODIS land surface temperature and emissivity products

Land surface temperature (LST) and land surface emissivity (LSE) are two key parameters in global climate study. This article aims to cross-validate LST/LSE products retrieved from data of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on board the first geostationary satellite, Meteosat Second Generation (MSG), with Moderate Resolution Imaging Spectroradiometer (MODIS) LST/LSE version 5 products over the Iberian Peninsula and over Egypt and the Middle East. Besides time matching, coordinate matching is another requirement of the cross-validation. An area-weighted aggregation algorithm was used to aggregate SEVIRI and MODIS LST/LSE products into the same spatial resolution. According to the quality control (QC) criterion and the view angle, the cross-validation was completed under clear-sky conditions and within a view angle difference of less than 5° for the two instruments to prevent land surface anisotropic effects. The results showed that the SEVIRI LST/LSE products are consistent with MODIS LST/LSE products and have the same trend over the two study areas during both the daytime and the night-time. The SEVIRI LST overestimates the temperature by approximately 1.0 K during the night-time and by approximately 2.0 K during the daytime compared to MODIS products over these two study areas. The SEVIRI LSE underestimates by about 0.015 in 11 μm and by about 0.025 in 12 μm over the Iberian Peninsula. However, both LSEs agree and show a difference of less than 0.01 over Egypt and the Middle East.     

MODIS数据反演地表温度的参数敏感性分析

在利用MODIS卫星遥感数据进行地表温度反演过程中，有两个基本参数需要确定，即地表比辐射率和大气透过率，尽管采用了比较合理的参数估计方法，但仍会有一些不可避免的因素导致误差的产生。为了进一步研究可能的参数误差对地表温度反演精度的影响，我们对该算法的两个参数进行敏感性分析。结果表明，当31、32两个波段的参数估计都有中等误差时，可能的地表温度误差对大气透过率和地表比辐射率都不敏感，所引起的地表温度误差大约为0．6～0．8℃，算法能够得到较高精度的地表温度反演结果。     

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.     

A practical split‐window algorithm for retrieving land‐surface temperature from MODIS data

This paper presents a practical split‐window algorithm utilized to retrieve land‐surface temperature (LST) from Moderate‐resolution Imaging Spectroradiometer (MODIS) data, which involves two essential parameters (transmittance and emissivity), and a new method to simplify Planck function has been proposed. The method for linearization of Planck function, how to obtain atmosphere transmittance from MODIS near‐infrared (NIR) bands and the method for estimating of emissivity of ground are discussed with details. Sensitivity analysis of the algorithm has been performed for the evaluation of probable LST estimation error due to the possible errors in water content and emissivity. Analysis indicates that the algorithm is not sensitive to these two parameters. Especially, the average LST error is changed between 0.19–1.1°C when the water content error in the simulation standard atmosphere changes between ?80 and 130%. We confirm the conclusion by retrieving LST from MODIS image data through changing retrieval water content error. Two methods have been used to validate the proposed algorithm. Results from validation and comparison using the standard atmospheric simulation and the comparison with the MODIS LST product demonstrate the applicability of the algorithm. Validation with standard atmospheric simulation indicates that this algorithm can achieve the average accuracy of this algorithm is about 0.32°C in LST retrieval for the case without error in both transmittance and emissivity estimations. The accuracy of this algorithm is about 0.37°C and 0.49°C respectively when the transmittance is computed from the simulation water content by exponent fit and linear fit respectively.     

New refinements and validation of the MODIS Land-Surface Temperature/Emissivity products

This paper discusses the lessons learned from analysis of the Moderate Resolution Imaging Spectroradiometer (MODIS) Land-Surface Temperature/Emissivity (LST) products in the current (V4) and previous versions, and presents eight new refinements for V5 product generation executive code (PGE16) and the test results with real Terra and Aqua MODIS data. The major refinements include considering surface elevation when using the MODIS cloudmask product, removal of temporal averaging in the 1 km daily level-3 LST product, removal of cloud-contaminated LSTs in level-3 LST products, and the refinements for the day/night LST algorithm. These refinements significantly improved the spatial coverage of LSTs, especially in highland regions, and the accuracy and stability of the MODIS LST products. Comparisons between V5 LSTs and in-situ values in 47 clear-sky cases (in the LST range from − 10 °C to 58 °C and atmospheric column water vapor range from 0.4 to 3.5 cm) indicate that the accuracy of the MODIS LST product is better than 1 K in most cases (39 out of 47) and the root of mean squares of differences is less than 0.7 K for all 47 cases or 0.5 K for all but the 8 cases apparently with heavy aerosol loadings. Emissivities retrieved by the day/night algorithm are well compared to the surface emissivity spectra measured by a sun-shadow method in two field campaigns. The time series of V5 MODIS LST product over two sites (Lake Tahoe in California and Namco lake in Tibet) in 2003 are evaluated, showing that the quantity and quality of MODIS LST products depend on clear-sky conditions because of the inherent limitation of the thermal infrared remote sensing.     

Evaluation of ASTER and MODIS land surface temperature and emissivity products using long-term surface longwave radiation observations at SURFRAD sites

Land surface temperature (LST) and emissivity are key parameters in estimating the land surface radiation budget, a major controlling factor of global climate and environmental change. In this study, Terra Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and Aqua MODerate resolution Imaging Spectroradiometer (MODIS) Collection 5 LST and emissivity products are evaluated using long-term ground-based longwave radiation observations collected at six Surface Radiation Budget Network (SURFRAD) sites from 2000 to 2007. LSTs at a spatial resolution of 90 m from 197 ASTER images during 2000-2007 are directly compared to ground observations at the six SURFRAD sites. For nighttime data, ASTER LST has an average bias of 0.1 °C and the average bias is 0.3 °C during daytime. Aqua MODIS LST at 1 km resolution during nighttime retrieved from a split-window algorithm is evaluated from 2002 to 2007. MODIS LST has an average bias of − 0.2 °C. LST heterogeneity (defined as the Standard Deviation, STD, of ASTER LSTs in 1 × 1 km² region, 11 × 11 pixel in total) and instrument calibration error of pyrgeometer are key factors impacting the ASTER and MODIS LST evaluation using ground-based radiation measurements. The heterogeneity of nighttime ASTER LST is 1.2 °C, which accounts for 71% of the STD of the comparison, while the heterogeneity of the daytime LST is 2.4 °C, which accounts for 60% of the STD. Collection 5 broadband emissivity is 0.01 larger than that of MODIS Collection 4 products and ASTER emissivity. It is essential to filter out the abnormal low values of ASTER daily emissivity data in summer time before its application.     

Development of a daily long term record of NOAA-14 AVHRR land surface temperature over Africa

We developed a new 6-year daily, daytime and nighttime, NOAA-14 AVHRR based land surface temperature (LST) dataset over continental Africa for the period 1995 through 2000. The processing chain was developed within the Global Inventory Modeling and Mapping System (GIMMS) at NASA's Goddard Space Flight Center. This paper describes the processing methodology used to convert the Global Area Coverage Level-1b data into LST and collateral data layers, such as sun and view geometries, cloud mask, local time of observation, and latitude and longitude. We used the Ulivieri et al. [Ulivieri, C., M.M. Castronuovo, R. Francioni, and A. Cardillo (1994), A split window algorithm for estimating land surface temperature from satellites, Adv. Space Research, 14(3):59-65.] split window algorithm to determine LST values. This algorithm requires as input values of surface emissivity in AVHRR channels 4 and 5. Thus, we developed continental maps of emissivity using an ensemble approach that combines laboratory emissivity spectra, MODIS-derived maps of herbaceous and woody fractional cover, and the UNESCO FAO soil map. A preliminary evaluation of the resulting LST product over a savanna woodland in South Africa showed a bias of < 0.3 K and an uncertainty of < 1.3 K for daytime retrievals (< 2.5 K for night). More extensive validation is required before statistically significant uncertainties can be determined. The LST production chain described here could be adapted for any wide field of view sensor (e.g., MODIS, VIIRS), and the LST product may be suitable for monitoring spatial and temporal temperature trends, or as input to many process models (e.g., hydrological, ecosystem).     

Estimating afternoon MODIS land surface temperatures (LST) based on morning MODIS overpass,location and elevation information

The MODerate Resolution Imaging Spectroradiometer (MODIS) instrument on‐board the Terra and Aqua satellites is a critical tool for providing daily estimates of land surface temperature (LST). Terra launched in late 1999 has a morning (AM) overpass, whereas Aqua launched in early 2002 has an afternoon (PM) overpass. Generally, LST is expected, under cloudless conditions, to be warmer in the early afternoon than the morning due to the link between maximum skin temperature and solar insolation peak time, therefore the Aqua PM LST is likely to be closer to the maximum daily LST than that acquired from Terra. This letter investigated differences between the Aqua MODIS PM and Terra MODIS AM LST estimates over a range of land cover classes, locations, and dates, across Canada. The aim was to develop a simple adjustment which can be applied to Terra AM LST estimates to approximate a “synthetic” Aqua PM LST product from 2000 to mid‐2002, thereby providing a seamless afternoon MODIS LST product from 2000 to 2006. Results indicate that there are statistically significant differences between the AM and PM LST ranging from 0.3°C to 3.2°C depending on cover type, and between 1.2° and 5.0° depending on time of year. On average, over 90% of the variation observed in the PM record can be explained by the AM LST, land cover types and location.     

AMSR-E亮温数据与MODIS陆表分类产品结合反演全球陆表温度

AMSR-E被动微波传感器获取的亮温数据与MODIS陆表分类产品(MOD12)相结合,将全球陆表分为16类,并假设每种类型的地表在各个被动微波通道具有较一致的发射率,在此基础上针对每种陆表类型分别建立了陆表温度反演算法。在算法的建立过程中,为了避免混合像元以及冻土、积雪发射率不确定性带来的影响,仅对单一地表类型占90%以上以及MODIS陆表温度产品高于273K的被动微波像元进行回归。同时,考虑到降雨对回归结果的影响,在数据选择中加入了降雨判识,在被动微波亮温数据中除去了降雨像元。利用上述算法,用2004年1～10月的全球部分地区AMSR-E数据在MODIS陆表分类产品的基础上对每种地表类型分别进行了陆表温度反演,并与MODIS陆表温度产品进行对比,结果显示相关性较好,均方根误差为2～4 K。     

基于MODIS 影像数据的劈窗算法研究及其参数确定

劈窗算法是目前由热红外遥感数据获取陆面温度的主要方法。在介绍劈窗算法的一般表现形式的基础上, 我们推导出适合于MOD IS 影像数据的劈窗算法。大气透过率和地表比辐射率是求解地表温度的两个关键参数。由于MOD IS 图像分辨率较低,MOD IS 像元主要由水面、植被和裸土3种地物类型构成, 故可依据这3 种地物的构成比例确定地表比辐射率。从遥感影像上反演大气的水汽含量, 再根据大气水汽含量与大气透过率的关系计算出大气透过率。最后将文中推导的劈窗算法用于江苏省地表温度的反演。反演出来的地表温度图显示出明显的地表温度空间差异、城市热岛效应和不同的地物类型。     

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.