Abstract

The aim of this study is to evaluate the aerosol influence on LST retrieval with two algorithms (split-window (SW) method and a four-channel based method) using simulated data under typical conditions. The results show that the root mean square error (RMSE) decreases to approximately 2.3 K for SW method and 1.5 K for four channel based method when VZA = 60° and visibility = 3 km; an RMSE would be increased by approximately 1.0 K when visibility varies from 3 km to 23 km. Moreover, a detailed sensitivity analysis under a visibility of 3 km and 23 km is performed in terms of uncertainties of land surface emissivity (LSE), water vapor content (WVC), and instrument noise, respectively. It is noted that the four-channel based method is more sensitive to LSE than SW method, especially for dry atmosphere; LST error caused by a WVC uncertainty of 20% is within 1.5 K for SW method and within 0.8 K for four-channel based method; the instrument noise would introduce LST error with a maximum standard deviation of 0.5 K and 0.04 K for the four-channel based method and SW method, respectively.

1. Introduction

As the direct driving force in the exchange of long-wave radiation and turbulent heat fluxes at the surface-atmosphere interface, land surface temperature (LST) is a key parameter in the research of land surface processes at global or regional scale, energy, and water cycle [1, 2]. LST has high spatial and temporal variation [3] and the traditional in situ measurement methods can only provide LST at point or local scales [4], so that it is unreliable to acquire LST at global scale by interpolation method. Remote sensing is a unique way of measuring LST at regional and global scales [5], and in the last twenty years, various infrared sensors onboard satellite have been launched successfully, such as Moderate Resolution Imaging Spectroradiometer (MODIS), Advanced Very High Resolution Radiometer (AVHRR), Advanced Along-Track Scanning Radiometer (AASTR), and Spinning Enhanced Visible and Infrared Imager (SEVIRI), among others. Meanwhile, the LST retrieval technology from remote sensing data made a great progress, and various methods have been proposed, such as single channel method [6, 7], split-window (SW) method [812], temperature emissivity separation (TES) method [13], and multichannel method [14, 15], among others. However, these algorithms just provide a multichannel correction for attenuation arising from molecular water vapor absorption under clear skies, and aerosol effect on LST retrieval is overlooked [16].

Atmospheric aerosols comprise a dispersed system of small solid and liquid particles suspended in air for varying periods of time and transported by vertical and horizontal wind currents, frequently to great distances. Even for clear skies, the particulates affect the heating or cooling of the atmosphere through scatter and absorption of energy, depending on the chemical and physical properties of the particles [17]. Therefore, there are some previous investigations into the problem on sea surface temperature [16, 1823]. However, because the uncertainty of accurately determining the land surface emissivity (LSE) is still crucial for improving LST retrieval algorithms, few studies have considered reducing the influences of different aerosol types and loadings on the retrieval of LST from remote sensing data.

The transmittance for the middle infrared (MIR) channel in atmosphere window is more stable with less change to the surface temperature, and most values are above 0.8 [15]. In consideration of this, a four-channel based method for LST algorithm using 2 thermal infrared (TIR) channels and 2 MIR channels has been developed [15]. Since the aerosol size range covers more than five orders of magnitude, from 10 nm to several hundred micrometers [17] and the aerosol effect decreases with the increase of wavelength, MIR channels are more sensitive to aerosol particles than TIR channels. Therefore, in this study, a detailed analysis on the performance of this method under typical conditions is investigated in terms of aerosol, LSE, water vapor content (WVC), instrument noise, using simulated airborne hyperspectral scanner (AHS) data, which has multiple TIR and MIR channels providing an excellent opportunity for evaluating the algorithms’ performances; meanwhile, an intercomparison with a two-TIR channel (11 μm and 12 μm) based SW method is performed. In addition, in this study, to avoid the interference of solar direct radiation, only the algorithms for nighttime LST retrieval are analyzed.

In this paper, Section 2 shows the basic theory in infrared spectrum. Section 3 describes the data characteristics, atmospheric and aerosol characteristics in AHS channels for LST inversion. Section 4 presents the simulated data. Section 5 gives the results and discussions. Finally, the conclusions and summary are drawn in Section 6.

2. Methodology

2.1. Radiative Transfer Model

Based on the radiative transfer theory and assuming a cloud-free atmosphere under local thermodynamic equilibrium, the infrared spectral radiance at the top of atmosphere (TOA) for a viewing zenith angle (VZA), , can be expressed as [24]withwhere is the TOA brightness temperature (BT), is the BT for surface-leaving radiance, and and are the directional spectral emissivity and the total atmospheric spectral transmittance at a VZA of , respectively; and represent the upwelling atmospheric spectral radiance and the upwelling diffusion radiance that results from the scattering of solar radiance at a VZA of , respectively; is the LST; and represent the downwelling hemispheric atmospheric spectral irradiance and the downwelling diffusion irradiance that results from the scattering of solar radiation divided by , respectively; is the bidirectional spectral reflectivity; is the solar zenith angle and is the relative azimuth angle between the view azimuth angle and the solar azimuth angle; and is the direct solar spectral irradiance at ground level. For the spectral radiance in the TIR channels and the nighttime measurements in the MIR channels, , , and can be neglected without loss of accuracy.

2.2. LST Retrieval Method

The TOA radiance measured by infrared sensors is influenced by surface properties (LST and LSE) and atmosphere. To retrieve LST physically, the LSE for each channel must be known and is acquired by classification. Under a given VZA, the SW method and the four-channel based method for each surface type can be expressed as follows [15]:(i)two-TIR channel based SW method (denoted as Algorithm 1):where is the LST for the th IGBP surface types and and are the TOA BT of the adjacent thermal channels. These parameters , , , and are the regression coefficients for the th surface type for nighttime LST retrievals.(ii)The four-channel based method (denoted as Algorithm 2):where and are the TOA BT of the two MIR channels in atmospheric window, respectively. These parameters are the regression coefficients for the th surface type for nighttime LST retrievals.

In this study, we mainly analyze the performances of the two LST retrieval methods affected by aerosol for two typical surface types (soil and grass surfaces) under dry (WVC: 0–1.5 g/cm2) and wet atmosphere (WVC: 4.0–5.5 g/cm2), respectively.

3. Optimal Channels Selection

The main problem in surface temperature retrieval from space is the modification of the surface radiation on its way through the atmosphere and the contribution of the atmosphere itself to the final signal measured by the sensor [25]. The content of atmospheric constituents as water vapor, ozone, and uniformly mixed gases (CO2, CO, CH4, N2O, and O2) as well as the temperature profile strongly affecting the accuracy of LST retrieval. In this section, the atmospheric characteristics in AHS channels are analyzed for retrieving LST with optimal channels.

3.1. AHS Description

The AHS is an imaging line-scanner radiometer, acquiring images in 80 spectral channels covering the visible and near infrared (VNIR, 20 channels), short wave infrared (SWIR, 43 channels), MIR (7 channels), and TIR (10 channels) spectral ranges. AHS has an instantaneous field of view (IFOV) of 2.5 mrad and a FOV of ±45° [26]. The spectral characteristics for the MIR and TIR channels are shown in Figure 1.

3.2. Atmosphere Characteristics in AHS Channels

Absorption by atmospheric molecules is a key parameter impacting the radiance measured by sensors; therefore, it is essential to analyze the absorption characteristics of atmospheric molecular in MIR and TIR spectrum for selecting the optimal channels for LST retrieval efficiently. As an example, the absorption properties of several atmospheric gases (H2O, CO2, O3, CO, N2O, and CH4) for an US standard atmosphere in MIR and TIR spectrum including urban aerosol are displayed in Figure 2. It is noted that the MIR range is divided by the several absorption subranges mainly induced by N2O and CO2 [25]. The ozone absorption splits the thermal infrared region into two parts such that one is more transparent than the other [25] (see Figure 2). In both infrared ranges, water vapor is the primary absorber and therefore substantially responsible for the shape of the total transmission. For MIR range, it is obvious that there exist two atmospheric windows of approximately 3.3~4.2 μm and 4.8~5.1 μm, and a strong absorption of CO2 is presented approximately from 4.2 μm to 4.4 μm, which is almost blind to land surface; CH4 molecular appears strong absorption approximately from 3.2 μm to 3.5 μm, and less absorption can be found in other spectrum; there exists strong absorption of N2O approximately from 4.4 μm to 4.6 μm; CO molecular has a relatively weak absorption approximately from 4.5 μm to 4.8 μm; other atmospheric molecules have relatively high transmittance in MTR spectrum. For TIR range, besides O3 absorption band ranging from 9 to 10 μm, H2O mainly contributes to the total transmittance.

In consideration of the locations of AHS spectral channels (see Figure 1) and the absorption characteristics of atmospheric molecular (see Figure 2), it is noted that AHS channel 64 (CH64) and channel 65 (CH65) are affected seriously by the strong absorption of CH4; AHS channel 67 (CH67) is influenced by a strong absorption by CO2 and N2O; AHS channel 69 (CH69) and channel 70 (CH70) are strongly affected by H2O. Furthermore, since AHS CH64 and AHS CH70 are too close to the edges of 3~5 μm atmosphere window and are usually rather noisy [26], the two channels AHS CH66 and CH68 weakly by atmospheric absorption are selected.

In addition, for TIR spectrum, AHS channels 71 (CH71) and 80 (CH80) show the highest atmospheric absorption; channels 75 (CH75) to 79 (CH79) are located in the atmospheric window 10~12.5 μm; channel 74 (CH74) is located in the region of the ozone absorption. In consideration of the combined effects of atmospheric molecular absorption and previous studies on LST retrieval for AHS data using SW method [27], AHS CH75 and CH79 are excellent channels for LST retrieval.

4. Data Simulation

In order to determine the coefficients and , a serial of AHS data under different atmospheric and surface conditions are simulated with the atmospheric radiative transfer code MODTRAN 4. In this study, the atmospheric profiles are extracted from the Thermodynamic Initial Guess Retrieval (TIGR) dataset at Laboratoire de Météorologie Dynamique [28, 29], and 705 clear-sky profiles with relative humidity at one of levels lower than 90% are extracted from TIGR dataset, covering a wide range of bottom atmospheric temperature (: 250~310 K) and atmospheric WVC (0.06~5.39 g/cm2).

To make the simulation more representatives, for analyzing the aerosol effect on LST retrieval over urban area, three types of aerosol, namely, rural, urban, and tropospheric aerosol [30], with visibility varying from 3 to 23 km with a step of 5 km are included, and the Mie theory is employed for the aerosol scattering. LST varies from to in step of 5 K for from 250 K to 310 K ( is the atmospheric temperature in the first boundary layer of the selected atmosphere profiles). In addition, the emissivity from ASTER Spectral library is used in this study, which includes spectra of rocks, minerals, lunar soils, terrestrial soils, manmade materials, meteorites, vegetation, snow, and ice covering the visible to thermal infrared wavelength region (0.4~15.4 μm) (http://speclib.jpl.nasa.gov/). As an example, the spectrum for soil and grass, used in this simulation, is shown in Figure 3. Furthermore, to account for angular effects, a total of 6 VZAs at a height of 100 km (0°, 33.56°, 44.42°, 51.32°, 56.25°, and 60°) are used in this simulation. Therefore, for a certain type of land surface, with these qualified atmospheric profiles, various aerosols, and VZAs as MODTRAN input parameters, the spectral atmospheric parameters, such as atmospheric transmittance, downwelling radiance, and upwelling radiance, can be acquired, and then the TOA radiance or BT is subsequently generated with the aid of land surface parameters (LST and LSE) and the channel spectral response function (see Figure 4).

Similar to LST retrieval schemes proposed by Wan and Dozier [31], for each VZA, the WVC and LST are divided into several tractable subranges for improving LST retrieval accuracy. In this paper, the LST is divided into five subranges with an overlap of 5 K, that is, <280 K, 275~295 K, 290~310 K, 305~325 K, and >320 K, and WVC is divided into five subranges with an overlap of 0.5 g/cm2, that is, 0~1.5 g/cm2, 1.0~2.5 g/cm2, 2.0~3.5 g/cm2, 3.0~4.5 g/cm2, and 4.0~5.5 g/cm2. For each land surface, the corresponding coefficients and for each subrange can therefore be determined by the statistical regression (Levenberg-Marquardt method) with these simulated data.

5. Results

5.1. LST Retrieval under Aerosol Contaminated Condition

For LST retrieval, it is customary to express channel-integrated TOA radiance as BT, which is a coupling of atmosphere and land surface information and could be attenuated by the aerosol existing in the atmosphere, producing an effect on LST retrieval. In this paper, we begin exploring the effects of aerosol on TOA BT using radiative transfer simulations in terms of aerosol optical depth (AOD) (it is also expressed as visibility in this paper) and aerosol types; then the LST retrieval accuracy and uncertainty analysis under different AOD conditions for the two methods are analyzed systematically, with an aim of investigating the performances of the methods, especially under the aerosol contaminated skies.

5.1.1. Aerosol Effect on TOA BT

The magnitude of aerosol-induced attenuation in the IR depends on a number of factors, the more important ones being AOD (or visibility), the effective temperatures (vertical profile) of the gaseous and particulate absorbing layers, the particle size distribution, and particle chemical composition. The latter two provide information about the aerosol attenuation mechanism which cannot be inferred objectively from AOD value and is depicted as aerosol types.

To analyze the influence of AOD on TOA BT, the simulated data with various visibilities under the six VZAs is used. For a certain aerosol type and surface type (grass or soil), the differences between the aerosol-contaminated TOA BT and aerosol-free TOA BT are calculated under dry (WVC: 0–1.5 g/cm2) and wet (WVC: 4–5.5 g/cm2) atmospheres, respectively. Generally, we can expect such factors to affect TOA BT in different channels. As an example, for urban aerosol type, the ΔBT caused by different visibilities under wet and dry atmosphere for soil surface are given in Figures 5 and 6, respectively. It is noted that the ΔBT is highly variable, and the AOD effect is much more serious in MIR channels (AHS CH66 and CH68) than that in TIR channels (AHS CH75 and CH79) due to the wavelength dependence of the aerosol scattering. Moreover, AOD under dry atmosphere causes much more TOA BT variation than that under wet atmosphere, since in case of a wet atmosphere the spectrum is dominated by water vapour and hence the aerosol does not play such a dominant role. It is noted that the ΔBT for AHS CH66 and CH68 is proportional to the increase of VZA, while inversely for AHS CH75 and CH79; for urban aerosol, AOD mostly produce positive ΔBT under wet or dry atmosphere.

In addition, the effect of aerosol type on TOA BT is investigated in terms of three typical aerosol types, namely, rural, urban, and tropospheric aerosols, under six VZAs with a visibility of 23 km, and the TOA BT differences under two different aerosol types are calculated. It is found that the TOA BT variation caused by aerosol type is much lower than that caused by AOD, especially for dry atmosphere. As an example, BT differences for soil surface under wet atmosphere is given in Figure 7, and an attenuation within 1.0 K is presented in channel 66, while the aerosol type almost produce slight attenuation on TOA BT in channel 79 among these three aerosol type. Therefore, as following, its effect on LST retrieval is not discussed, and only AOD is focused.

5.1.2. The Effect of Aerosol on LST Retrieval

As described above, aerosol could produce a great influence on TOA BT, especially for MIR channels, which will directly affect LST retrieval accuracy. However, few studies have addressed the problem of aerosol effect in the existing LST retrieval algorithms, even if the aerosol loading retrieval algorithms from the remote sensing data have been developing for many years. In this section, the LST retrieval accuracies of Algorithm 1 and Algorithm 2 under various atmosphere and visibility conditions are investigated with simulated data.

For urban aerosol type, the root mean square errors (RMSEs) of LST retrieval with Algorithm 1 and Algorithm 2 are calculated under various VZAs, five visibilities (3, 8, 13, 18, and 23 km), two types of atmosphere (wet and dry), and two land surfaces (grass and soil) (see Figures 8 and 9). It is worth noting that both Algorithm 1 and Algorithm 2 are affected by aerosol when visibility is lower than 8 km; Algorithm 2 still presents lower RMSE than Algorithm 1, especially for wet atmosphere; when the VZA is equal to 60° and visibility is equal to 3 km, Algorithm 1 presents a maximal RMSE of approximately 2.3 K, while Algorithm 2 presents 1.5 K; RMSE of approximately 1.0 K is increased when visibility varies from 3 km to 23 km; the RMSE for soil is 0.2 K higher than that for grass.

5.2. Uncertainty Analysis under Aerosol Contaminated Condition

LSE, atmospheric correction error, the instrument noise, and the accuracy of the algorithm itself, and so forth affect the accuracy of LST retrieval. To analyze the uncertainty of LST retrieval under aerosol contaminated condition, a sensitivity analysis in terms of the uncertainties of LSE, WVC, and instrument noise under 3 km and 23 km visibility is performed for Algorithm 1 and Algorithm 2 using the simulated data with a VZA of 0° to eliminate VZA effect.

5.2.1. The Sensitivity to Uncertainty of LSE

The LST retrieval algorithms depend on an accurate knowledge of land surface types since the LSE is acquired by classification in this paper. The largest LST uncertainty would be induced by misclassifying soil into grass, and the LST with error-added LSE for soil (denoted as ) is retrieved with the coefficients for grass. Under dry and wet atmospheres, the LST error caused by LSE uncertainty is expressed by the differences between and the error-free LST (denoted as ) retrieved with the right coefficients (see Figure 10) when the visibility is equal to 3 km and 23 km, respectively. It is demonstrated that the LST error caused by the uncertainty of LSE under wet atmosphere is larger than that under dry atmosphere, and Algorithm 2 is more sensitive to the misclassification than Algorithm 1 due to large variation of LSE in MIR channels. Furthermore, the standard deviation of LST error under 3 km visibility is larger than that under 23 km visibility; under 3 km and 23 km visibility, for Algorithm 1, the LST error caused by the misclassification is approximately within 1.0 K under dry atmosphere, while more than 2.0 K under wet atmosphere; for Algorithm 2, this error is more than 2.0 K under various atmosphere conditions.

5.2.2. The Sensitivity to Uncertainty of WVC

For some channels in the infrared atmospheric window, the main atmospheric effect comes from water vapor absorption and has important effect to LST retrieval error. In this study, the atmospheric WVC is used to select the coefficients in (3) and (4). The WVC error may result in wrong subrange selection and subsequently lead to a large error in the retrieved LST. This is also one of the reasons for two consecutive WVC subranges overlapping by 0.5 g/cm2 in the algorithm development. For example, the falling into two subranges and ; thus, the retrieved LSTs using the coefficients of the subrange with and are different from the retrieved LSTs using the coefficients of the subrange with , . In order to characterize the WVC error on the LST retrieval, an error of 20% is added to original WVC to generate the error-added WVC, which is then used to select the coefficients in (3) and (4), and the simulated data with the overlap WVC are used to estimate the LST with different coefficients of the two adjacent subranges. The histograms of the LST difference between them for dry and wet atmospheres with a visibility of 3 km and 23 km, respectively, are shown in Figure 11. It is noted that Algorithm 1 is more sensitive to the uncertainty of WVC than Algorithm 2; under 3 km or 23 km visibility, LST error caused by a WVC uncertainty of 20% is within 1.5 K for Algorithm 1, and within 0.8 K for Algorithm 2; in most cases, the LST errors under 3 km visibility is a little larger than those under 23 km visibility.

5.2.3. The Sensitivity to Instrument Noise

The expected radiometric noise of AHS CH66, 68, 75, and 79 is 0.33 K, 0.33, 0.2 K, and 0.2 K, respectively. In order to analyze the instrument noise effect on LST retrieval, the corresponding instrument radiometric noises for these channels are added to their TOA BT; then, the noised LST is estimated again using Algorithm 1 and Algorithm 2, respectively. The differences between noise-free and are calculated for dry and wet atmosphere with a visibility of 3 km and 23 km, respectively, and the corresponding histograms are shown in Figure 12. It can be found that the for Algorithm 2 is a little higher than that for Algorithm 1 since the higher instrument noise would produce when more channels are used to retrieved LST; the maximum standard deviation of for Algorithm 2 is approximately 0.5 K, while it is 0.04 K for Algorithm 1 under various atmosphere and visibility conditions.

6. Conclusion and Summary

This paper has addressed an analysis of aerosol effect on LST retrieval using two methods with the data simulated with MODTRAN model derived by cloud-free TIGR atmosphere profiles, various land surface conditions (LSE and LST), VZAs, and spectral response functions. The analysis results show that, under various AOD conditions, Algorithm 2 presents higher LST retrieval accuracy than Algorithm 1, with a maximal RMSE of approximately 2.3 K for Algorithm 1 and 1.5 K for Algorithm 2 when the VZA is equal to 60° and visibility is equal to 3 km. Besides, LST retrieval accuracy would be reduced by approximately 1.0 K when visibility varying from 3 km to 23 km.

Furthermore, an analysis under two typical AOD (visibility = 3 km and 23 km) is performed in terms of uncertainties of LSE, WVC and instrument noise with simulated data. The analysis results show that the Algorithm 2 is more sensitive to LSE than that for Algorithm 1 due to large variation of LSE in MIR channel. When soil is misclassified into grass, the LST error for Algorithm 1 is within 1.0 K, and more than 2.0 K under dry and wet atmosphere, respectively; this error for Algorithm 2 is more than 2.0 K under various atmosphere conditions. In addition, a WVC uncertainty of 20% would introduce a LST error within 1.5 K for Algorithm 1 and within 0.8 K for Algorithm 2. The maximum standard deviation of the LST error caused by the instrument noise is approximately 0.5 K for Algorithm 2, while it is 0.04 K for Algorithm 1.

In summary, aerosol indeed produces influence on LST retrieval, and Algorithm 2 integrating MIR and TIR channels presents lower RMSE than Algorithm 1, especially under low visibility condition, which would be benefit to developing an aerosol correction algorithm for LST retrieval. However, due to the larger variation in LSE for MIR channels, Algorithm 2 is more sensitive to LSE than Algorithm 1, so Algorithm 2 could be used effectively with the accurate knowledge of LSE.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work has been supported by the National Natural Science Foundation of China under Grants 41301387 and 41371353. The authors would like to thank anonymous reviewers for the valuable comments and suggestions.