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.     

Separating surface emissivity and temperature using two-channel spectral indices and emissivity composites and comparison with a vegetation fraction method

The temperature-independent thermal infrared spectral indices (TISI) method is employed for the separation of land surface temperature (LST) and emissivity from surface radiances (atmospherically corrected satellite data). The daytime reflected solar irradiance and the surface emission at ∼3.8 μm have comparable magnitudes. Using surface radiances and a combination of day-night 2-channel TISI ratios, the ∼3.8 μm reflectivity is derived. For implementing the TISI method, coefficients for NOAA 9-16 AVHRR channels are obtained. A numerical analysis with simulated surface radiances shows that for most surface types (showing nearly Lambertian behavior) the achievable accuracy is ∼0.005 for emissivity (AVHRR channel-5) and ∼1.5 K for LST. Data from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used for calculation of atmospheric attenuation. Comparisons are made over a part of central Europe on two different dates (seasons). Clouds pose a major problem to surface observations; hence, monthly emissivity composites are derived. Additionally, using TISI-based monthly composites of emissivities, a normalized difference vegetation index (NDVI)-based method is tuned to the particular study area and the results are intercompared. Once the coefficients are known, the NDVI method is easily implemented but holds well only for vegetated areas. The error of the NDVI-based emissivities (with respect to the TISI results) ranges between −0.038 and 0.032, but for vegetated areas the peak of the error-histogram is at ∼0.002. The algorithm for retrieving emissivity via TISI was validated with synthetic data. Due to the different spatial scales of satellite and surface measurements and the lack of homogeneous areas, which are representative for low-resolution pixels and ground measurements, ground-validation is a daunting task. However, for operational products ground-truth validation is necessary. Therefore, also an approach to identify suitable validation sites for meteorological satellite products in Europe is described.     

SEBASS hyperspectral thermal infrared data: surface emissivity measurement and mineral mapping

This study focuses on mapping surface minerals using a new hyperspectral thermal infrared (TIR) sensor: the spatially enhanced broadband array spectrograph system (SEBASS). SEBASS measures radiance in 128 contiguous spectral channels in the 7.5- to 13.5-μm region with a ground spatial resolution of 2 m. In September 1999, three SEBASS flight lines were acquired over Virginia City and Steamboat Springs, Nevada. At-sensor data were corrected for atmospheric effects using an empirical method that derives the atmospheric characteristics from the scene itself, rather than relying on a predicted model. The apparent surface radiance data were reduced to surface emissivity using an emissivity normalization technique to remove the effects of temperature. Mineral maps were created with a pixel classification routine based on matching instrument- and laboratory-measured emissivity spectra, similar to methods used for other hyperspectral data sets (e.g. AVIRIS). Linear mixtures of library spectra match SEBASS spectra reasonably well, and silicate and sulfate minerals mapped remotely, agree with the dominant minerals identified with laboratory X-ray powder diffraction and spectroscopic analyses of field samples. Though improvements in instrument calibration, atmospheric correction, and information extraction would improve the ability to map more pixels, these hyperspectral TIR data nevertheless show significant advancement over multispectral thermal imaging by mapping surface materials and lithologic units with subtle spectral differences in mineralogy.     

Air temperature estimation with MSG-SEVIRI data: Calibration and validation of the TVX algorithm for the Iberian Peninsula

Air temperature can be estimated from remote sensing by combining information in thermal infrared and optical wavelengths. The empirical TVX algorithm is based on an estimated linear relationship between observed Land Surface Temperature (LST) and a Spectral Vegetation Index (NDVI). Air temperature is assumed to be equal to the LST corresponding to the effective full vegetation cover, and is found by extrapolating the line to a maximum value of NDVI_max. The algorithm has been tested and reported in the literature previously. However, the effect of vegetation types and climates and the potential variation in NDVI of the effective full cover has not been subject for investigation. The present study proposes a novel methodology to estimate NDVI_max that uses observed air temperature to calibrate the NDVI_max for each vegetation type. To assess the validity of this methodology, we have compared the accuracy of estimates using the new NDVI_max and the previous NDVI_max that have been proposed in literature with MSG-SEVIRI images in Spain during the year 2005. In addition, a spatio-temporal assessment of residuals has been performed to evaluate the accuracy of retrievals in terms of daily and seasonal variation, land cover, landscape heterogeneity and topography. Results showed that the new calibrated NDVI_max perform well, with a Mean Absolute Error ranging between 2.8 °C and 4 °C. In addition, vegetation-specific NDVI_max improve the accuracy compared with a unique NDVI_max.     

Effect of surface wind speed and sensor view zenith angle dependence of emissivity on SST retrieval from thermal infrared data: ATSR

The Along Track Scanning Radiometer (ATSR) data are currently being processed at various places within the European community including the Rutherford Appleton Laboratory (RAL) in the U.K. In generating an atmospherically corrected sea-surface temperature (SST) field, the emissivity of the sea surface is assumed to be independent of the sensor view zenith angle, sea state and wavelength (Ian Barton, RAL, personal communication). The sensor view zenith angle dependence of the emissivity is generally not known because of the complications introduced by the surface wind speed. This paper attempts to evaluate the uncertainties introduced in the SST due to the variation of emissivity with the sensor view zenith angle and surface roughness generated by the wind speed.

Using the Cox and Munk formulation, Takashima and Takayama have simulated the sea water emissivities as a function of wind speed of up to 15ms^-1 and a range of the sensor view zenith angles. Their emissivity data for 11 and 12μm channels corresponding to the viewing geometry of the ATSR have been used in this work. It is shown that, depending on the value of the SST, there can be significant errors due to the sensor view zenith angle and sea surface roughness dependence of the emissivity. For example, if the SST is 10°C and the wind speed is 0ms^-1, then the errors due to the sensor view zenith angle dependence of the emissivity are shown to be 0·77°C and 0·55°C in 11 and 12μm channels, respectively, and at 15 ms^-l the respective errors reach about 1·ldeg K and 0·86 deg K. The errors due to the deviations of the emissivities from unity for nadir view in calm conditions are about 2·1°C and 3·5°C, respectively, in the 11 and 12μm channels. All of these errors are additive, indicating the importance of the calibration and validation.     

基于静止气象卫星数据的地表温度遥感估算

基于分裂窗算法和地表温度日周期变化模型,探讨了利用多时相热红外遥感数据反演地表温度的方法。首先,利用分裂窗算法及地表温度日周期变化形式,推导了多时相遥感数据反演地表温度的方法。其次,利用辐射传输模型(MODTRAN),以2006年夏季在禹城观测的3 d地表温度、气温及大气水汽数据做为输入参数、变化观测角及比辐射率,模拟了一日多个时刻与风云二号(F-2D)波谱响应函数一致的亮温数据,基于此,模拟数据库对所提算法进行了检验。最后,利用2010年9月30日FY-2D多时相热红外数据对新疆区域地表温度进行了反演,并与相应时刻的MODIS地表温度产品进行了比较。结果表明：利用模拟遥感数据反演地表温度,模拟值与估算值的相关系数达0.9,均方根误差在1.5 K以内;利用在轨FY-2D热红外数据反演得到的地表温度与MODIS温度产品趋势基本一致,两者的相关性达到了0.5,均方根误差为4.4 K。需要说明的是,此方法仅满足于晴朗无云的条件。     

High spatial resolution sea surface climatology from Landsat thermal infrared data

High spatial and temporal resolution maps of sea surface temperature (SST) have numerous applications in coastal and estuarine systems. A climatology map, tracking SST as a function of year-day, was produced at Southern New England using 53 Landsat TM and ETM+ thermal infrared data. A recursive curve-fitting algorithm was used to fit these data and eliminate cloud contamination, resulting in an average daily temperature at every 60-m pixel. The climatology was validated against long-term in situ records that were analyzed with the same techniques. The results show, as expected, that isolated and shallow water bodies undergo more extreme temperature variation (−2 to 25 °C) than deeper, well-connected embayments (1 to 21 °C) or the coastal ocean (4 to 18 °C). The coastal ocean is shown to lag insolation and shallow lakes by up to 44 days, with embayments showing a gradation between these extremes. Despite the subtle temperature range variation, there is rich detail in the spatial patterns which are relevant to the applied sciences of coastal and estuarine systems. The spatial pattern of the climatology reveals anomalous patterns, such as occur where anthropogenic forcing alters climatological patterns. The heat budget of Mount Hope Bay in northeast Narragansett Bay has anthropogenic thermal input from a large power plant, and this input is reflected in the climatology. From the results, it is seen that Narragansett Bay has, on average, a mean annual temperature of 11.86±0.41 °C, while the Mount Hope Bay system is consistently warmer at 12.30±0.21 °C and shows a delayed response to autumn cooling. The long history of Landsat data acquisition can be used to create a climatology of coastal and estuarine scale dynamics at an order of magnitude finer scale resolution than AVHRR climatologies.     

Bare surface soil moisture retrieval from the synergistic use of optical and thermal infrared data

Land surface soil moisture (SSM) is crucial to research and applications in hydrology, ecology, and meteorology. To develop a SSM retrieval model for bare soil, an elliptical relationship between diurnal cycles of land surface temperature (LST) and net surface shortwave radiation (NSSR) is described and further verified using data that were simulated with the Common Land Model (CoLM) simulation. In addition, with a stepwise linear regression, a multi-linear model is developed to retrieve daily average SSM in terms of the ellipse parameters x₀ (horizontal coordinate of the ellipse centre), y₀ (vertical coordinate of the ellipse centre), a (semi-major axis), and θ (rotation angle), which were acquired from the elliptical relationship. The retrieval model for daily average SSM proved to be independent of soil type for a given atmospheric condition. Compared with the simulated daily average SSM, the proposed model was found to be of higher accuracy. For eight cloud-free days, the root mean square error (RMSE) ranged from 0.003 to 0.031 m³ m^?3, while the coefficient of determination (R²) ranged from 0.852 to 0.999. Finally, comparison and validation were conducted using simulated and measured data, respectively. The results indicated that the proposed model showed better accuracy than a recently reported model using simulated data. A simple calibration decreased RMSE from 0.088 m³ m^?3 to 0.051 m³ m^?3 at Bondville Companion site, and from 0.126 m³ m^?3 to 0.071 m³ m^?3 at the Bondville site. Coefficients of determination R² = 0.548 and 0.445 were achieved between the estimated daily average SSM and the measured values at the two sites, respectively. This paper suggests a promising avenue for retrieving regional SSM using LST and NSSR derived from geostationary satellites in future developments.     

Land Surface Temperature Retrieval and Its Spatial Heterogeneity in Summer in Chongqing

Land Surface Temperature(LST)is an important parameter in land surface energy budget.In order to improve the accuracy of LST retrieval by remote sensing methods in summer in urban districts of Chongqing with hot and humid atmosphere condition.an improved methodology was presented with the improved atmospheric transmittance estimated on MODTRAN software using the atmospheric profile data of MERRA in urban districts of Chongqing.LST was retrieved from Landsat 8 TIRS band 10 data using single-window algorithm which apply the improved and unimproved atmospheric transmittance respectively.Then the retrieved LST was compared with the 0 cm soil temperatures observedby 4 meteorological stations.Finally,the spatial heterogeneity of LST was analyzed.The result indicated that：(1)The scheme proposed in this paper can improve LST retrieval in summer in urban districts of Chongqing.The Mean Absolute Error(MAE)decrease from 4.89 K to 1.73 K.(2)The retrieved LST has spatial heterogeneity with different terrain factors.Its lapse rate is about 1.17 K/100 m.It decreases with the increase of slope.Moreover,it has obvious differences with aspect.The flat slope>sunny slope>semi-sunny slope>semi-shady slope>shady slope.There also existed highly significant correlation between the LST and hill shade.The LST increases with the decrease of hill shade.(3)Influenced by land cover,the spatial distribution of LST showed significant differences.The average LST inthe built\|up area is highest,while the wet land is lowest.The difference of average LST in other land cover types is little.     

Comparison of land surface emissivity and radiometric temperature derived from MODIS and ASTER sensors

This study compared surface emissivity and radiometric temperature retrievals derived from data collected with the MODerate resolution Imaging Spectroradiometer (MODIS) and Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) sensors, onboard the NASA's Earth Observation System (EOS)-TERRA satellite. Two study sites were selected: a semi-arid area located in northern Chihuahuan desert, USA, and a Savannah landscape located in central Africa. Atmospheric corrections were performed using the MODTRAN 4 atmospheric radiative transfer code along with atmospheric profiles generated by the National Center for Environmental Predictions (NCEP). Atmospheric radiative properties were derived from MODTRAN 4 calculations according to the sensor swaths, which yielded different strategies from one sensor to the other. The MODIS estimates were then computed using a designed Temperature-Independent Spectral Indices of Emissivity (TISIE) method. The ASTER estimates were derived using the Temperature Emissivity Separation (TES) algorithm. The MODIS and ASTER radiometric temperature retrievals were in good agreement when the atmospheric corrections were similar, with differences lower than 0.9 K. The emissivity estimates were compared for MODIS/ASTER matching bands at 8.5 and 11 μm. It was shown that the retrievals agreed well, with RMSD ranging from 0.005 to 0.015, and biases ranging from −0.01 to 0.005. At 8.5 μm, the ranges of emissivities from both sensors were very similar. At 11 μm, however, the ranges of MODIS values were broader than those of the ASTER estimates. The larger MODIS values were ascribed to the gray body problem of the TES algorithm, whereas the lower MODIS values were not consistent with field references. Finally, we assessed the combined effects of spatial variability and sensor resolution. It was shown that for the study areas we considered, these effects were not critical.     

Near-surface thermal structure and surface diurnal warming in the Adriatic Sea using satellite and drifter data

An extensive set of in situ temperature data collected by surface drifters is combined with satellite-derived sea surface temperature images to study the difference between the pseudo-bulk and bulk temperatures (ΔT_pb-b) in the Adriatic Sea for the period 21 September 2002-31 December 2003. The variations of this temperature difference are described as a function of local wind speed and incoming solar radiation provided by a local area atmospheric model. The daily sea surface temperature variability is also assessed by computing the temperature difference between the daily maximal and minimal values (ΔT_day-night). The data show that the smaller the wind speed and the larger the solar radiation, the larger ΔT_pb-b. The temperature difference reached the highest value (∼5 °C) on a hot day (more than 600 W/m²) of May 2003 in weak wind condition (around 3 m/s). For strong winds (speed > 6 m/s) the dependence on both the wind and solar radiations vanishes as the temperature difference approaches zero because the near-surface water becomes thermally homogenous due to the wind-induced vertical mixing. Strong diurnal warming of the sea surface, as derived by the pseudo-bulk estimates, and a strong near-surface stratification were found during the spring/summer season. Monthly mean statistics show that the diurnal cycle of the pseudo-bulk and bulk temperature starts to become significant already in February and March. Subsequently (from April to August) both the diurnal warming and the stratification are maximal (monthly means of ΔT_day-night ∼1-2 °C and of ΔT_pb-b ∼0.5 °C ), while in fall and early winter the ΔT_pb-b values are quite small (monthly means near 0 °C) and the ΔT_day-night monthly means are bounded by 0.5-1.5 °C. Maximal amplitudes of the diurnal cycle can exceed 4 °C (mostly in spring-summer) for both the pseudo-bulk and bulk temperatures. However, the monthly means of ΔT_day-night is generally twice as large for the pseudo-bulk estimates (∼2 °C) with respect to the bulk layer (∼1 °C). The diurnal warming of the sea surface, as derived by the pseudo-bulk temperature, occurs at about 14:30 local time, that is more than 2 h after the maximal sun elevation and an hour earlier than the bulk temperature maximum at 20-40 cm depth.     

Coupled retrieval of aerosol optical thickness, columnar water vapor and surface reflectance maps from ENVISAT/MERIS data over land

An algorithm for the derivation of atmospheric parameters and surface reflectance data from MEdium Resolution Imaging Specrometer Instrument (MERIS) on board ENVIronmental SATellite (ENVISAT) images has been developed. Geo-rectified aerosol optical thickness (AOT), columnar water vapor (CWV) and spectral surface reflectance maps are generated from MERIS Level-1b data over land. The algorithm has been implemented so that AOT, CWV and reflectance products are provided on an operational manner, making no use of ancillary parameters apart from those attached to MERIS products. For this reason, it has been named Self-Contained Atmospheric Parameters Estimation from MERIS data (SCAPE-M). The fundamental basis of the algorithm and applicable error figures are presented in the first part of this paper. In particular, errors of ± 0.03, ± 4% and ± 8% have been estimated for AOT, CWV and surface reflectance retrievals, respectively, by means of a sensitivity analysis based on a synthetic data set simulated under a usual MERIS scene configuration over land targets. The assumption of a fixed aerosol model, the coarse spatial resolution of the AOT product and the neglection of surface reflectance directional effects were also identified as limitations of SCAPE-M. Validation results are detailed in the second part of the paper. Comparison of SCAPE-M AOT retrievals with data from AErosol RObotic NETwork (AERONET) stations showed an average Root Mean Square Error (RMSE) of 0.05, and an average correlation coefficient R² of about 0.7-0.8. R² values grew up to more than 0.9 in the case of CWV after comparison with the same stations. A good correlation is also found with the MERIS Level-2 ESA CWV product. Retrieved surface reflectance maps have been successfully compared with reflectance data derived from the Compact High Resolution Imaging Spectrometer (CHRIS) on board the PRoject for On-Board Autonomy (PROBA) in the first place. Reflectance retrievals have also been compared with reflectance data derived from MERIS images by the Bremen AErosol Retrieval (BAER) method. A good correlation in the red and near-infrared bands was found, although a considerably higher proportion of pixels was successfully processed by SCAPE-M.     

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

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

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+.     

Land surface phenology of North American mountain environments using moderate resolution imaging spectroradiometer data

Monitoring and understanding plant phenology is becoming an increasingly important way to identify and model global changes in vegetation life cycle events. High elevation biomes cover twenty percent of the Earth's land surface and provide essential natural resources. These areas experience limited resource availability for plant growth, development, and reproduction, and are one of the first ecosystems to reflect the harmful impact of climate change. Despite this, the phenology of mountain ecosystems has historically been understudied due to the rough and variable terrain and inaccessibility of the area. In addition, although numerous studies have used synoptically sensed data to study phenological patterns at the continental and global scales, relatively few have focused on characterizing the land surface phenology in mountainous areas. Here we use the MODIS/Terra + Aqua satellite 8-day 500 m Nadir BRDF Adjusted Reflectance product to quantify the land surface phenology. We relate independent data for elevation, slope, aspect, solar radiation, and temperature as well as longitude and latitude with the derived phenology estimates. We present that satellite derived SOS can be predicted based on topographic and weather variables with a significant R²adj between 0.56 and 0.62 for the entire western mountain range. Elevation and latitude exhibit the most significant influences on the timing of SOS throughout our study area. When examined at both the local and regional scales, as well as when accounting for aspect and temperature, SOS follows closely with Hopkins' Bioclimatic Law with respect to elevation and latitude.     

MODIS数据地表温度反演劈窗算法比较

针对MODIS数据,分析比较了QIN和Wan-Dozier两种劈窗算法地表温度(LST)反演精度和误差分布。首先利用辐射传输模型MODTRAN4.0,结合TIGR大气廓线数据,评价两种算法绝对精度,然后基于误差传递理论分析评价二者的总精度,最后对两种算法的LST反演结果进行比较。研究表明针对所有廓线数据,两种算法绝对精度相差不大,但Wan-Dozier算法绝对精度受地表温度和水汽含量变化的影响程度要大于QIN算法;两种算法总精度相差不大,且主要误差源均为算法绝对精度和地表比辐射率精度,QIN算法反演结果对地表比辐射率的敏感性要略高于Wan-Dozier算法;两种算法得到研究区LST分布情况基本一致,均可表现空间LST分布差异,其中水体和裸土的LST反演结果差异较大,城镇和植被平均温度差异在0.5 K以内。     

Atmospheric correction of optical imagery from MODIS and Reanalysis atmospheric products

In this paper we analyze the differences obtained in the atmospheric correction of optical imagery covering bands located in the Visible and Near Infra-Red (VNIR), Short-Wave Infra-Red (SWIR) and Themal-Infrared (TIR) spectral regions when atmospheric profiles extracted from different sources are used. In particular, three sensors were used, Compact High Resolution Imaging Spectrometer (CHRIS), Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) and Landsat5 Thematic Mapper (TM), whereas four atmospheric profiles sources were considered: i) local soundings launched near the sensor overpass time, ii) Moderate Resolution Radiometer (MODIS) atmospheric profiles product (MOD07), iii) Atmospheric Correction Parameter Calculator (ACPC) generated by the National Center for Environmental Prediction (NCEP) and iv) Modified Atmospheric Profiles from Reanalysis Information (MAPRI), which includes data from NCEP and National Center of Atmospheric Research (NCAR) Reanalysis project but interpolated to 34 atmospheric levels and resampled to 0.5° × 0.5°. MODIS aerosol product (MOD04) was also used to extract Aerosol Optical Thickness (AOT) values at 550 nm. Analysis was performed for three test dates (12th July 2003, 18th July 2004 and 13th July 2005) over an agricultural area in Spain. Results showed that air temperature vertical profiles were similar for the four sources, whereas dew point temperature profiles showed significant differences at some particular levels. Atmospheric profiles were used as input to MODTRAN4 radiative transfer code in order to compute atmospheric parameters involved in atmospheric correction, with the aim of retrieving surface reflectances in the case of VNIR and SWIR regions, and Land Surface Temperature (LST) in the case of the TIR region. For the VNIR and SWIR region, significant differences depending on the atmospheric profile used were not found, particularly in the Visible region in which the AOT content is the main parameter involved in the atmospheric correction. In the case of TIR, differences depending on the atmospheric profile used were appreciable, since in this case the main parameter involved in the atmospheric correction is the water vapor content, which depends on the vertical profile. In terms of LST retrieval from ASTER data (2004 test case), all profiles provided satisfactory results compared to the ones obtained when using a local sounding, with errors of 0.3 K for ACPC and MAPRI cases and 0.7 K for MOD07. When retrieving LST from TM data (2005 test case), errors for MOD07 and MAPRI were 0.6 and 0.9 K respectively, whereas ACPC provided an error of 2 K. The results presented in this paper show that the different atmospheric profile sources are useful for accurate atmospheric correction when local soundings are not available. In particular, MOD07 product provides atmospheric information at the highest spatial resolution, 5 km, although its use is limited from 2000 to present, whereas MAPRI provides historical information from 1970 to present, but at lower spatial resolution.     

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).     

通过AVHRR数据研究中国陆面温度分异规律

近年来国内外利用遥感方法在陆面温度精确反演中开展了大量的研究工作,采用了一个在大区域上适用的由NOAA？AVHRR数据反演陆面温度的方法,反演中国晴空条件下各月和全年平均陆面温度,分析了中国陆面温度的分异规律,并与气温的分异规律作了对比,同时对中国土地利用／土地覆盖变化研究（LUCC）样带上的陆面温度变化进行了分析,这项工作从晴空条件重新认识了地面温度场的空间分异,对于研究中国陆地土壤蒸发,植物光合作用,土地覆盖的分布具有重要的指示意义。     

Scaling of land surface temperature using satellite data: A case examination on ASTER and MODIS products over a heterogeneous terrain area

Land surface temperature (LST) is a key parameter in numerous environmental studies. Surface heterogeneity induces uncertainty in pixel-wise LST. Spatial scaling may account for the uncertainty, however, different approaches lead to differences in scaled values. Satellite-retrieved LST may be representative of the pixel-wise LST and useful for scaling analysis, but the limited accuracy of retrieved values adds uncertainty into the scaled values. Based on the Stefan-Boltzmann (S-B) law, this study proposed scaling approaches for LST over flat and relief areas to explore the combined uncertainties in scaling using satellite-retrieved data. To take advantage of simultaneous, multi-resolution observations at coincident nadirs by the Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and the MODerate-resolution Imaging Spectroradiometer (MODIS), LST products from these two sensors were examined for part of the Loess Plateau in China. 90-m ASTER LST data were scaled up to 1 km using the proposed approaches, and variation in the LST was generally reduced after scaling. Amongst the sources of uncertainties, surface heterogeneity (emissivity) and different scaling approaches resulted in very minor differences, with a maximum difference of 0.2 K for the upscaled LST. Terrain features, taken as an areal weighting factor, had negligible effects on the upscaled value. Limited accuracy of the retrieved LST was the major uncertainty. The overall LST increased 0.6 K on average with correction for terrain-induced angular effect and 0.4 K for both angular and adjacency effects over the study area. Accounting for terrain correction in scaling is necessary for rugged areas. With terrain correction, the upscaled ASTER LST achieved an agreement of − 0.1 ± 1.87 K and a root mean square error (RMSE) of 1.87 K overall with the 1-km MODIS LST rectified by Wan et al.'s approach [Wan, Z., Zhang, Y., Zhang Q., Li, Z.-L. (2002b), Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data. Remote Sensing of Environment, 83, 163-180]. Refining the rectification approach resulted in a better agreement of − 0.2 ± 1.57 K and a RMSE of 1.58 K.