1. Introduction
Land surface temperature (LST) is an important parameter for formulating surface–atmosphere interactions (Dickinson 1996; Sobrino et al. 2003a,b, 2004). It controls the upward terrestrial radiation and surface–atmosphere sensible and latent heat fluxes, and it is an indicator of the energy balance at the earth’s surface (Sellers et al. 1988). LST retrievals are more challenging than sea surface temperature (SST) retrievals because of surface heterogeneity and the spatial variability of surface emissivity.
Attempts to derive LST from satellite observations have been ongoing for several decades. The focus was on polar-orbiting systems, such as the Advanced Very High Resolution Radiometer (AVHRR) and the Moderate Resolution Imaging Spectroradiometer (MODIS), because of global coverage and high spatial resolution. Much less attention has been given to geostationary satellites. Price (1984) applied the AVHRR SST split-window algorithm to agricultural land and found the retrieval accuracy to be about 3 K. Because the land surface emissivity is not equal to 1 and has spectral dependence (Becker 1987), Becker and Li (1990) developed a local split-window method (referred to here as BL90) with coefficients dependent on surface emissivity. Wan and Dozier (1996) optimized the split window algorithm developed by Becker and Li (1990) by computing coefficients that vary with satellite zenith angle (SZA), water vapor, and lower boundary air temperature [known as the “generalized split window (GSW) algorithm”, also referred to here as WD96]. Retrievals of LST from satellites are made in atmospheric windows that have relatively small water vapor absorption. To fully avoid such effects, several studies have added atmospheric correction to the GSW algorithms (Prata 1993; Sobrino et al. 1994; Becker and Li 1995; Francois and Ottle 1996; Coll and Caselles 1997). Such modifications improved the standard error in the temperature estimates to approximately ±1 and ±1.5 K, with no significant bias (Sun et al. 2002). Qin et al. (2002a,b) used the National Oceanic and Atmospheric Administration (NOAA) AVHRR observations to study the sand-dune region across the Israel–Egypt border. They found that the 3-K across-border difference during summer is due to differences in emissivities of biogenic crust and bare sand (being 0.972 and 0.954, respectively). Bhattacharya and Dadhwal (2003) also retrieved and evaluated land (soil–vegetation complex) surface temperature from NOAA/AVHRR thermal images over India. They found that in May the retrieved LST values are closer to the soil temperature. This was attributed to the lower fractional vegetation cover as characterized by the normalized difference vegetation index. To study the effects of different surface emissivities and the three-dimensional structure of topography on LST, Chen et al. (2004) indicate that using different emissivity values for soil and leaves will produce a 1-K difference in LST estimates. Ouaidrari et al. (2002) found that quadratic form is more accurate than the linear form, especially for high water vapor contents.
Oku and Ishikawa (2004) used images from the Geostationary Meteorological Satellite Visible/Infrared Spin-Scan Radiometer to estimate diurnal variations of LST distributions over the Tibetan Plateau. LST retrieval from a Geostationary Operational Environmental Satellite (GOES-8) was addressed by Faysash and Smith (1999), Sun and Pinker (2003), and Sun et al. (2004, 2006a,b).
Evaluation of LST algorithms is complicated because LST is highly variable in space and time and ground measurements are not compatible with the spatial resolution of satellites. Observations from high-resolution satellite platforms (e.g., Landsat) are desirable; however, they are of low temporal resolution (about once every 16 days) and, as such, are of interest as a “proof of concept.” There is a strong demand for information on LST for evaluation of numerical weather prediction models at resolutions of about 40 km. As such, there is a need to evaluate the quality of the LST information derivable from satellites at scales of interest to the modeling community. The disparity between the ground and satellite observations is unavoidable, and the best possible matchup between ground and satellite observations is unable to resolve this intrinsic problem. The objective of this study is to address this complex issue, namely, How well do satellite grid averages of LST compare to available observations of LST?
We present a comprehensive evaluation of a split-window LST algorithm developed by Sun and Pinker (2003) (SP-SW) against ground observations as available from several sources and from independent algorithms. Results are compared with the NOAA/National Environmental Satellite, Data, and Information Service (NESDIS) operational daytime algorithm over the United States (Wu et al. 1999), the same algorithm implemented with an independent cloud screening approach developed at the University of Maryland for daytime and nighttime (Li et al. 2007), the previously published split-window algorithms with emissivity correction (Price 1984; Prata and Platt 1991; Ulivieri and Camizzaro 1985; Ulivieri et al. 1992; Vidal 1991; Sobrino et al. 1994, 1993; Coll et al. 1994), the widely used generalized split-window algorithms (Becker and Li 1990; Wan and Dozier 1996), and the Becker and Li (1995) algorithm with water vapor correction. The evaluation method is presented in section 2. In section 3 the data used for evaluation are described. The algorithms used are introduced in section 4. Section 5 presents the evaluation results for the different algorithms. Errors of LST estimation are addressed in section 6. Section 7 presents conclusions.
2. Method
The split-window algorithm of Sun and Pinker (2003) as well as the respective triple-widow algorithm will be evaluated against ground truth and against several existing algorithms. To make the algorithm intercomparison consistent and relevant for the satellite observations used, it is necessary to rederive relevant coefficients of the algorithms in a systematic manner using the same inputs as used for the SP algorithms. The coefficients from the previously published algorithms are rederived using the same forward simulations as performed in this study with the Moderate Resolution Atmospheric Transmission (MODTRAN) code for GOES-8.
In what follows, the simulations that were performed and the algorithms that were used are described. A table is presented with all of the algorithms and the newly derived coefficients.
a. Simulations performed
To derive regression coefficients that relate satellite-observed brightness temperature to LST, we perform forward simulations for GOES-8 using the MODTRAN4.0 radiative transfer model. To make the LST simulation results applicable on a global scale, the National Centers for Environmental Prediction (NCEP) reanalyzed global surface temperature dataset with matching atmospheric profiles with 17 vertical layers at a resolution of 2.5° is used for forward simulations. The solar zenith angle was calculated according to date, time of day, latitude, and longitude for each point in this dataset. The column water vapor (CWV) was integrated from the profiles; the air temperature was also obtained from the NCEP reanalysis. The range for CWV (precipitable water) is from 0.0 to 7.0 cm. Surface–air temperature differences range from −36.1 to 30.0 K. The Department of Geography, University of Maryland (UMD), generated a 1-km-resolution global land cover product (Hansen et al. 2000; http://glcf.umiacs.umd.edu/data/landcover/index.shtml). These products include the following 14 International Geosphere–Biosphere Programme (IGBP) classes: 1) water, 2) evergreen needle-leaf forest, 3) evergreen broadleaf forest, 4) deciduous needle-leaf forest, 5) deciduous broadleaf forest, 6) mixed forest, 7) woodland, 8) wooded grassland, 9) closed shrubland, 10) open shrubland, 11) grassland, 12) cropland, 13) bare ground, and 14) urban and built-up. In addition to these 14 IGBP classes, we add snow/ice to replace the bare ground at high latitudes, thus totaling 15 classes. The land classification was aggregated first to the resolution of the global simulations at 2.5° and subsequently to the analyzed data resolutions of 0.5° and 0.125°. In the aggregation process, the land cover type in each grid box was assigned on the basis of the dominant surface type and the surface emissivity was calculated according to a linear mixing with weighted sum of the land cover percentage times the emissivity of this surface type; it varied from 0.58 to 0.98 for the 3.9-μm channel (Fig. 1 of Sun and Pinker 2003).
From the simulations, the coefficients for the GSW algorithm (Wan and Dozier 1996) are derived as follows: 1) the atmospheric temperature profiles are aggregated according to the 11-μm brightness temperature, T11 ≤ 280 K or > 280 K; 2) the atmospheric CWV is aggregated into intervals of 0.5 cm; 3) day/night is separated using solar zenith angle θs (θs ≤ 87.5° is daytime and > 87.5° is nighttime; and 4) viewing angles are separated into intervals of 2°. The emissivity needed for the formulation of the split-window algorithm is assumed according to surface types following the approach proposed by Snyder et al. (1998).
b. Algorithms used in evaluation
1) Sun and Pinker (2003) split-window algorithm
2) Sun and Pinker (2003) triple-window algorithm (nighttime only)
3) The NOAA/NESDIS operational algorithm
Original version using channel-1 cloud detection
Original version using multichannel cloud detection
To allow the implementation of the NOAA/NESDIS operational algorithm during both daytime and nighttime, we use a four-channel coupled cloud and snow detection algorithm (CCSDA) for both daytime and nighttime (Li et al. 2007) and couple it with the original NOAA/NESDIS algorithm to rederive LST (NOAA-R algorithm). The CCSDA cloud detection method also has the capability to detect clouds over snow (Pinker et al. 2007).
4) Split-window algorithms with emissivity correction
The formulas for all of the algorithms of this type that were used in the intercomparison are presented in Table 2. Some details about the GSW algorithms are provided here.
Wan and Dozier (1996) optimized the local split-window algorithm developed by Becker and Li (1990) by computing coefficients that vary with SZA, water vapor, and lower boundary air temperature (known as the generalized split window), presented in Table 2 as No. 9. In this formulation, atmospheric attenuation with optical path is accounted for. In accord with this, the coefficients vary with cold/warm conditions, day/night, viewing angles, and column water vapor.
3. Data used
a. Satellite observations
Satellite information needed for implementing the various algorithms includes cloud cover fraction and brightness temperatures for clear sky radiances in the 11.0- and 12.0-μm channels from GOES-8. These are produced as a by-product from the operational NOAA/NESDIS GOES surface shortwave radiative flux product (Pinker et al. 2003) developed in support of the Global Energy and Water Cycle Experiment (GEWEX) Continental Scale International Project (GCIP) and the GEWEX Americas Prediction Project as archived at the University of Maryland (http://www.atmos.umd.edu/~srb/gcip/). This information is at hourly time scales over the domain area of 125°–70°W and 25°–50.5°N, at 0.5° and at 0.125° spatial resolutions. To avoid cloud effects in LST retrievals, 10% cloud thresholds are used (if cloud amount is less than 10%, LST will be retrieved). Because cloud detection criteria are very strict, the compromise of a 10% threshold was necessary to ensure that a sufficiently sizeable sample of observations is available for the retrievals (experiments were performed with different percentages before the 10% value was selected).
b. Ground observations
Ground observations of surface skin temperature needed for the evaluation of satellite-based estimates are not readily available. The observations differ from each other, and a certain adjustment is needed to make them compatible with each other. Some of the observations are based on infrared radiometers while others measure upwelling longwave radiative fluxes. Some of the radiometric observations used have been converted to skin temperatures by data providers while others are converted in this study using information on surface emissivity.
1) ARM/GCIP NESOB-97 30-min skin temperature composite
The Atmospheric Radiation Measurement Program (ARM)/GCIP Near-Surface Observation Dataset-1997 (NESOB-97) came from the ARM–Cloud and Radiation Test Bed site (34°–39°N and 94.5°–100.5°W) during the period 1 April 1997–31 March 1998. The surface skin temperature used in this paper is observed at 30-min intervals at the Central Facility (36.6°N, 97.48°W) of the Southern Great Plains (SGP) site (Table 3). The instrument used to observe the skin temperature is the multifilter radiometer (MFR). It is a sensing head of the multifilter shadow band radiometer that is mounted to look downward to measure upwelling spectral irradiances. The filters of the MFRs can drift significantly and induce errors of about ±5%. The upwelling irradiance is converted to skin temperature using a NOAA/Atmospheric Turbulence and Diffusion Division (ATDD) algorithm (http://www.arm.gov/instruments/instrument.php?id=mfrsr).
2) The Surface Radiation Budget Network
The evaluation of the various algorithms against ground observation is done at monthly time scale over the diurnal cycle.
3) The Oklahoma Mesonet
The Oklahoma Mesonet is an automated network of over 110 remote, meteorological stations across Oklahoma (Brock et al. 1995; McPherson et al. 2007) (http://www.mesonet.org). During the period of this study, 10 sites measured surface energy fluxes using the eddy correlation method. These sites are referred to as the Oklahoma Atmospheric Surface-Layer Instrumentation System (OASIS) Super Sites. The sensors at OASIS Super Sites allow estimating net radiation (shortwave and longwave components) and sensible, latent, and ground heat fluxes (Basara and Crawford 2002). Surface skin temperature is measured using an infrared temperature sensor (IRT) manufactured by Apogee Instruments, Inc., mounted at 2 m (Fiebrich et al. 2003) and LST is converted from the upwelling longwave radiation with unit emissivity while surface reflection of downwelling longwave radiation is ignored (Fiebrich et al. 2003). Sun et al. (2004) discussed the effects of these two factors and found that the total effect may be a slight underestimation of the skin temperature. Surface skin temperature is measured at an additional 79 sites designated as OASIS Standard Sites.
Observations such as skin temperature and longwave radiation are referred to as research measurements. The research measurements do not have the routine quality assurance applied to other parameters (McPherson et al. 2007). Even so, for this study, a combination of automated and manual tests was applied to the data provided before use to ensure the skin temperature observations are of research quality.
4. Results
a. Model evaluation resulting from MODTRAN simulations
To select algorithms for comparison with ground truth, simulation results for several independent algorithms were first analyzed (Figs. 2 –5). From Figs. 3 and 4 it is evident that the Ulivieri2 algorithm with emissivity difference term shows improvement when compared with the Ulivieri1 algorithm that uses only an average emissivity correction term, especially for small viewing angles (<5°). The Vidal, Coll, and Ulivieri2 algorithms show similar results; the Sobrino2 algorithm with a nonlinear term of brightness temperature difference between the two split-window channels shows improved bias error for small satellite viewing angles (<6°) while the Price algorithm performance is below the others. As evident from Fig. 4, the Becker and Li (1995) (referred to as BL95) algorithm with water vapor correction shows further improvement as compared with the Sobrino2 algorithm for satellite viewing angles of less than 5°, and the GSW algorithm that considers viewing angle effects shows improvement at large viewing angles (>5°) as compared with the Ulivieri1, Ulivieri2, Vidal, Coll, Price, and Sobrino2 algorithms. The SP-SW algorithm with viewing angle correction and nonlinear terms shows smaller errors at both small and large viewing angles as compared with the other algorithms. Because the Ulivieri2, Vidal, and Coll algorithms are similar in structure and give similar results, the later Coll algorithm was selected as a representative of this group. Included in the evaluation against ground truth are the Price, Ulivieri1, Coll, Sobrino2, GSW, BL95, and SP-SW algorithms.
b. Evaluation of algorithms
The retrieved and ground-measured LST data are stratified as daytime and nighttime according to the solar zenith angle (≤87.5° is daytime, and >87.5° is nighttime). The retrieved LST in a grid closest to the observing site is used to match with ground observations. Then the data are aggregated for each season and the root-mean-square (RMS) error, bias, and standard deviation (STD) are calculated.
1) Evaluations at the ARM SGP site
Split-window-type algorithms
Detailed error statistics for different algorithms at the ARM SGP site are presented in Table 4. When compared with other algorithms, the LST retrieval errors from the SP-SW algorithm are the smallest, followed by GSW, Coll, Ulivieri1, Price, NOAA-R, BL95, and Sobrino2 algorithms. The algorithms selected for further evaluation are GSW, NOAA-R, BL95, and SP-SW.
The cloud-screening algorithm applied to the satellite data allows the detection of snow. Under snow conditions, the SP-SW algorithm as appropriate for snow conditions used in the simulations is utilized; for the other algorithms, snow emissivity (ε11 = 0.996; ε12 = 0.962) is assigned. Figure 6 presents evaluation results for the daytime LST retrievals from the SP-SW algorithm against ARM SGP observations for different seasons of 1997. This algorithm has low bias and RMS errors in winter, summer, and autumn, with the largest errors in the spring. Figure 7 compares the RMS error differences between retrievals aggregated at 1/2° and 1/8° resolution during nighttime as evaluated against ARM SGP observations in 1997. For all four seasons, the results from the SP-SW algorithm are in best agreement with observations. Figure 8 compares the bias and residual standard errors evaluated at Bondville during daytime in winter for 1996–2000. The NOAA, GSW, and SP-SW algorithms do not show large differences in LST errors at different resolutions, and the BL95 results show larger errors at higher resolution.
Triple-window algorithm versus split-window algorithm
The Sun and Pinker triple-window (SP-TW) algorithm is applicable only during nighttime. It is of interest to compare it with the SP-SW algorithm for the same time period. As evident from Table 5, during nighttime, LST retrieval errors from the SP-TW are lower than those from the SP-SW algorithm.
2) Evaluations at the SURFRAD sites
Figure 9 shows the time series of bias and STD errors for daytime LST retrievals resulting from the different algorithms for the period from 1996 to 2000. The bias error shows overestimation of LST from the BL95 algorithm. The NOAA algorithm shows larger underestimation of LST when compared with the SP-SW algorithm (Fig. 9a). Among the tested algorithms, the SP-SW algorithm is in best agreement with ground observations, giving the lowest STD errors of about 1.0 K from 1996 to 2000 (Fig. 9b). During nighttime, the NOAA algorithm underestimates LST most (Fig. 9a). The STD errors are smallest from the SP-SW algorithm and largest from the BL95 algorithm (Fig. 9b). All the algorithms show larger underestimates of bias in spring and autumn. This may be because, as shown in Fig. 6, LST is in the range of 280–300 K during spring and autumn, when the water vapor is the highest (Sun and Pinker 2003).
The seasonal daytime and nighttime bias of LST retrievals evaluated at Fort Peck from 1996 to 2000 is shown in Fig. 10. The LST retrieval errors from the SP-SW algorithm are the lowest, and the average bias during the 5-yr period is about 1 K; the underestimate of LST from the NOAA algorithm can be as large as −8 K.
3) Evaluations at the Oklahoma Mesonet sites
Results from the evaluation of the various LST retrievals against the Oklahoma Mesonet measurements are given in Table 6. The NOAA/NESDIS algorithm using the multichannel cloud screening gives the lowest errors when compared with other algorithms, followed by the SP-SW algorithm, generalized split-window algorithm, and BL95 algorithm. The bias error (accuracy) of the SP-SW algorithm is positive, and the bias of the NOAA algorithm implemented with the CCSDA cloud screening is negative. The apparent “better” performance of the NOAA algorithm at the mesonet sites seems to be related to the fact that the measured skin temperature is converted from the energy detected by IRT using the Stefan–Boltzmann law and assumed surface emissivity of 1.0 (Fiebrich et al. 2003). The actual emissivity at the Mesonet sites should be lower than 1.0 and therefore, the IRT-measured skin temperature is lower than the actual temperature of the ground.
4) Evaluations at MODIS test sites
Limited observations are available from homogeneous sites where data were collected in support of the validation of MODIS LST products (Wan et al. 2002, 2004). To be specific, data from the following sites are available (Wan et al. 2002): two locations at Mono Lake, California; three locations at Bridgeport, California; and two rice fields in California.
The data are matched to the GOES overpass for seven cases, and errors are presented as a difference between the ground observation and the GOES estimated value. A similar pattern was followed for the SP-SW algorithm as shown in Table 7. The LST from the SP-SW algorithm show very close agreement with the ground observations.
5. Discussion of error sources
a. Observational error
The filters of the MFRs at the ARM SGP sites drift during windy weather conditions, resulting in errors in the observed upwelling irradiance that can reach ±5%. The instruments are installed at 10 and 25 m, and the skin temperature is obtained from the upwelling irradiance by applying the NOAA/ATDD algorithm. The skin temperature is converted from measured upwelling irradiance with an assumed emissivity value, which might result in some errors. Errors in the LST retrievals include scattering patterns induced by heterogeneous surfaces (Weng et al. 2004), uniform surfaces having smaller scattering effects than mountain areas. The accuracy of the IRT used by the Oklahoma Mesonet for skin temperature measurement is approximately 0.2 K from 288 to 305 K and 0.3 K from 278 to 318 K (Bugbee et al. 1998).
b. Land cover classification
DeFries et al. (1995) indicate that land cover classification accuracy of 70%–90% can be achieved, and the correct typing probability for 1 km can be expected to be 80% (Loveland and Belward 1997). Errors from land cover classification will produce errors in the coefficients relating land cover types to emissivity variability (Sun and Pinker 2003).
c. Surface emissivity
Prata et al. (1995) showed that 1% emissivity uncertainty can affect LST values for each pixel by about 0.7°C. Yu et al. (2008) tested the sensitivity of 10 different split-window algorithms to uncertainties in surface emissivity and found that the split-window algorithms are very sensitive to emissivity values. Sun and Pinker (2003) performed model (MODTRAN) simulations with input errors introduced directly to emissivity or to land cover. The results indicated that the “error” in the simulated brightness temperature increases with the increase in emissivity error, whereas the simulated brightness temperature error in the split-window channels due to land cover error was very small (less than 0.2 K), even with 50% of pixels being misclassified. The SP-SW algorithm proposes to use land cover information with land cover–dependent regression coefficients and shows improvements to the LST algorithm if emissivities are used. The surface emissivity effect is considered implicitly through land cover types.
d. Atmospheric correction error
In the atmospheric window channels, water vapor is the major absorber of upwelling radiation, especially near the surface. Most of the water vapor in the atmosphere is for conditions within the temperature range of 280–305 K; it can vary from 0.25 to 7 cm (Sun and Pinker 2003). This may be the reason why negative bias errors (underestimates due to water vapor absorption) are usually higher at the LST range of 280–305 K as compared with those below 280 K.
e. Reliability of the evaluation process
Ground observations of LST are point measurements; satellites measure at pixel level. Land surface temperature is not homogeneous within a pixel area. When we match the ground observations to the satellite retrievals, we use satellite LST at the closest grid point at the same time as the ground measurements. The spatial difference in the satellite–ground match process may introduce errors.
f. Cloud contamination
To ensure that a sufficiently sizeable sample of observations is available for the retrievals, a 10% cloud threshold is used. This may include some cloud contamination and cause underestimation of LST or a negative bias. On the other hand, ground measurements are taken under all cloud conditions and may include cloud effects. All of these error sources are listed in sections 5a to 5f and may cause LST retrieval errors.
6. Summary and discussion
Because of the difficulty in obtaining ground observations of LST, the available observations that are directly and indirectly linked to LST are used in the evaluation process of satellite retrievals. Evaluation against the ARM and SURFRAD observations shows that the LST retrievals from the SP-SW algorithms are in better agreement with the ground observations than are the other split-window-type algorithms tested here. Evaluation against the mesonet observations shows that the NOAA algorithm is in closest agreement with the observations, followed by the SP-SW algorithm, GSW algorithm, and the BL95 algorithm. The apparent “better” agreement of the NOAA algorithms with the Mesonet site observations is due to the fact that the “observed” LST used was converted from the measured energy units using a surface emissivity of 1.0 (Fiebrich et al. 2003) as assumed in the NOAA/NESDIS algorithm. Overall, the SP-SW algorithm is in better agreement with ground observations than are the other algorithms that were tested, and the triple-window (SP-TW) algorithm shows improvement over the split-window (SP-SW) algorithm for nighttime LST retrieval.
Support from Grants NA86GP0202 and NAG59916 from the NOAA Office of Global Programs is appreciated. Data were obtained from the Atmospheric Radiation Measurement Program sponsored by the U.S. Department of Energy, Office of Science, Office of Biological and Environmental Research, Environmental Sciences Division. The installation of the Oklahoma Mesonet’s skin temperature sensors was made possible, in part, by an NSF MRI Grant (ATM-9724594). Continued funding for maintenance of the network is provided by the taxpayers of Oklahoma through state taxes. In addition, support from the NOAA Office of Global Programs (NOAA Grant NA17RJ1227) was instrumental in the development of research-quality datasets. We appreciate the reviewers for their helpful comments to improve this manuscript.
REFERENCES
Basara, J. B., , and K. C. Crawford, 2002: Linear relationships between root-zone soil moisture and atmospheric processes in the planetary boundary layer. J. Geophys. Res., 107 , 4274. doi:10.1029/2001JD000633.
Becker, F., 1987: The impact of spectral emissivity on the measurement of land surface temperature from a satellite. Int. J. Remote Sens., 8 , 1508–1522.
Becker, F., , and Z. L. Li, 1990: Toward a local split window method over land surface. Int. J. Remote Sens., 11 , 369–393.
Becker, F., , and Z. L. Li, 1995: Surface temperature and emissivity at various scales: Definition, measurement and related problems. Remote Sens. Rev., 12 , 225–253.
Bhattacharya, B. K., , and V. K. Dadhwal, 2003: Retrieval and validation of land surface temperature (LST) from NOAA AVHRR thermal images of Gujarat, India. Int. J. Remote Sens., 24 , 1197–1206.
Brock, F. V., , K. C. Crawford, , R. L. Elliott, , G. W. Cuperus, , S. J. Stadler, , H. L. Johnson, , and M. D. Eilts, 1995: The Oklahoma Mesonet: A technical overview. J. Atmos. Oceanic Technol., 12 , 5–19.
Bugbee, B., , M. Droter, , O. Monje, , and B. Tanner, 1998: Evaluation and modification of commercial infrared transducers for leaf temperature measurement. Adv. Space Res., 22 , 1425–1434.
Chen, L. F., , Z-L. Li, , Q. H. Liu, , S. Chen, , Y. Tang, , and B. Zhong, 2004: Definition of component effective emissivity for heterogeneous and non-isothermal surfaces and its approximate calculation. Int. J. Remote Sens., 25 , 231–244.
Coll, C., , and V. Caselles, 1997: A split-window algorithm for land surfaces temperature from advanced very high-resolution radiometer data: Validation and algorithm comparison. J. Geophys. Res., 102 , 16697–16713.
Coll, C., , V. Caselles, , J. A. Sobrino, , and E. Valor, 1994: On the atmospheric dependence of the split-window equation for land surface temperature. Int. J. Remote Sens., 15 , 105–122.
DeFries, R. S., , M. Hansen, , and J. Townshend, 1995: Global discrimination of land cover types from metrics derived from AVHRR pathfinder data. Remote Sens. Environ., 54 , 209–222.
Dickinson, R. E., 1996: Land surface processes and climate modeling. Bull. Amer. Meteor. Soc., 76 , 1445–1448.
Faysash, D. A., , and E. A. Smith, 1999: Simultaneous land surface temperature–emissivity retrieved in the infrared split window. J. Atmos. Oceanic Technol., 16 , 1673–1689.
Fiebrich, C. A., , J. E. Martinez, , J. A. Brotzge, , and J. B. Basara, 2003: The Oklahoma Mesonet’s skin temperature network. J. Atmos. Oceanic Technol., 20 , 1496–1504.
Francois, C., , and C. Ottle, 1996: Atmospheric corrections in the thermal infrared: Global and water vapor dependent split-window algorithms—Applications to ATSR and AVHRR data. IEEE Trans. Geosci. Remote Sens., 34 , 457–470.
Hansen, M. C., , R. S. Defries, , J. R. G. Townshend, , and R. Sohlberg, 2000: Global land cover classification at 1 km spatial resolution using a classification tree approach. Int. J. Remote Sens., 21 , 1331–1364.
Hicks, B. B., , J. J. DeLuisi, , and D. Matt, 1996: The NOAA Integrated Surface Irradiance Study (ISIS)—New surface radiation monitoring network. Bull. Amer. Meteor. Soc., 77 , 2857–2864.
Li, X., , R. T. Pinker, , M. M. Wonsick, , and Y. Ma, 2007: Toward improved satellite estimates of short-wave radiative fluxes—Focus on cloud detection over snow: 1. Methodology. J. Geophys. Res., 112 , D07208. doi:10.1029/2005JD006698.
Loveland, T. R. & , and A. S. Belward The IGBP-DIS global 1 km land cover data set, DISCover: First results. Int. J. Remote Sens., 18 , 289–3295.
McClain, E. P., , W. G. Pichel, , and C. C. Walton, 1985: Comparative performance of AVHRR-based multichannel sea surface temperatures. J. Geophys. Res., 90 , 11587–11601.
McPherson, R. A., and Coauthors, 2007: Statewide monitoring of the mesoscale environment: A technical update on the Oklahoma Mesonet. J. Atmos. Oceanic Technol., 24 , 301–321.
Oku, Y., , and H. Ishikawa, 2004: Estimation of land surface temperature over the Tibetan Plateau using GMS data. J. Appl. Meteor., 43 , 548–561.
Ouaidrari, H., , S. N. Goward, , K. P. Czajkowski, , J. A. Sobrino, , and E. Vermote, 2002: Land surface temperature estimation from AVHRR thermal infrared measurements: An assessment for the AVHRR Land Pathfinder II data set. Remote Sens. Environ., 81 , 114–128.
Pinker, R. T., and Coauthors, 2003: Surface radiation budgets in support of the GEWEX Continental-Scale International Project (GCIP) and the GEWEX Americas Prediction Project (GAPP), including the North American Land Data Assimilation System (NLDAS) project. J. Geophys. Res., 108 , 8844. doi:10.1029/2002JD003301.
Pinker, R. T., , X. Li, , E. Yegorova, , and W. Meng, 2007: Toward improved satellite estimates of short-wave radiative fluxes—Focus on cloud detection over snow: 2. Results. J. Geophys. Res., 112 , D09204. doi:10.1029/2005JD006699.
Prata, A. J., 1993: Land surface temperature derived from the advanced very high resolution radiometer and the along-track scanning radiometer: 1. Theory. J. Geophys. Res., 98 , 16689–16702.
Prata, A. J., , and C. M. R. Platt, 1991: Land surface temperature measurements from the AVHRR. Proc. Fifth AVHRR Data Users Conf., Tromso, Norway, EUM P09, 443–438.
Prata, A. J., , V. Caselles, , C. Coll, , J. A. Sobrino, , and C. Ottle, 1995: Thermal remote sensing of land surface temperature from satellite: Current status and future prospects. Remote Sens. Rev., 12 , 175–224.
Price, J. C., 1984: Land surface temperature measurements from the split window channels of the NOAA 7 Advanced Very High Resolution Radiometer. J. Geophys. Res., 89 , 7231–7237.
Qin, Z. H., , P. Berliner, , and A. Karnieli, 2002a: Micrometeorological modeling to understand the thermal anomaly in the sand dunes across the Israel–Egypt border. J. Arid Environ., 51 , 281–318.
Qin, Z. H., , A. Karnieli, , and P. Berliner, 2002b: Remote sensing analysis of the land surface temperature anomaly in the sand-dune region across the Israel–Egypt border. Int. J. Remote Sens., 23 , 3991–4018.
Sellers, P. J., , F. G. Hall, , G. Asrar, , D. E. Strebel, , and R. E. Murphy, 1988: The First ISLSCP Field Experiment (FIFE). Bull. Amer. Meteor. Soc., 69 , 22–27.
Snyder, W., , Z. Wan, , Y. Zhang, , and Y. Feng, 1998: Classification-based emissivity for land surface temperature measurement from space. Int. J. Remote Sens., 19 , 2753–2774.
Sobrino, J. A., , Z. L. Li, , M. Ph Stoll, , and F. Becker, 1993: Determination of the surface temperature from ATSR data. Proc. 25th Int. Symp. on Remote Sensing of Environment, Graz, Austria, ERIM, II-19–II-109.
Sobrino, J. A., , Z. L. Li, , M. Ph Stoll, , and F. Becker, 1994: Improvements in the split-window technique for land surface temperature determination. IEEE Trans. Geosci. Remote Sens., 32 , 243–253.
Sobrino, J. A., , J. C. Jimenez-Munoz, , and L. Paolini, 2003a: Land surface temperature retrieval from LANDSAT TM 5. Remote Sens. Environ., 90 , 434–440.
Sobrino, J. A., , J. E. Kharraz, , and Z. L. Li, 2003b: Surface temperature and water vapour retrieval from MODIS data. Int. J. Remote Sens., 24 , 5161–5182.
Sobrino, J. A., , J. C. Jimenez-Munoz, , J. El-Kharraz, , M. Gomez, , M. Romaguera, , and G. Soria, 2004: Single-channel and two-channel methods for land surface temperature retrieval from DAIS data and its application to the Barrax site. Int. J. Remote Sens., 25 , 215–230.
Sun, D., , and R. T. Pinker, 2003: Estimation of land surface temperature from Geostationary Operational Environmental Satellite (GOES-8). J. Geophys. Res., 108 , 4326. doi:10.1029/2002JD002422.
Sun, D., , R. T. Pinker, , and J. B. Basara, 2004: Land surface temperature estimation from the next generation of geostationary operational environmental satellites: GOES M-Q. J. Appl. Meteor., 43 , 363–372.
Sun, D., , M. Kafatos, , R. Pinker, , and D. Easterling, 2006a: Seasonal variations in diurnal temperature range from satellites and surface observations. IEEE Trans. Geosci. Remote Sens., 44 , 2779–2785.
Sun, D., , R. Pinker, , and M. Kafatos, 2006b: Diurnal temperature range over the United States: A satellite view. Geophys. Res. Lett., 33 , L05705. doi:10.1029/2005GL024780.
Sun, Y-Y., , F-M. Gottsche, , F-S. Olesen, , and H. Fisher, 2002: Retrieval of land surface temperature from combined AVHRR data. Ann. Geophys., 20 , 1257–1259.
Ulivieri, C., , and G. Cannizzaro, 1985: Land surface temperature retrievals from satellite measurements. Acta Astronaut., 12 , 977–985.
Ulivieri, C., , M. M. Castronouvo, , R. Francioni, , and A. Cardillo, 1992: A SW algorithm for estimating land surface temperature from satellites. Adv. Space Res, 14 , 59–65.
Vidal, A., 1991: Atmospheric and emissivity correction of land surface temperature measured from satellite using ground measurements or satellite data. Int. J. Remote Sens, 12 , 2449–2460.
Wan, Z., , and J. Dozier, 1996: A generalized split-window algorithm for retrieving land-surface temperature from space. IEEE Trans. Geosci. Remote Sens., 34 , 892–905.
Wan, Z., , Y. Zhang, , Q. Zhang, , and Z. Li, 2002: Validation of the land-surface temperature products retrieved from Terra Moderate Resolution Imaging Spectroradiometer data. Remote Sens. Environ., 83 , 163–180.
Wan, Z., , Y. Zhang, , Q. Zhang, , and Z. Li, 2004: Quality assessment and validation of the MODIS global land surface temperature. Int. J. Remote Sens., 25 , 261–274.
Weng, Q., , D. Lu, , and J. Schubring, 2004: Estimation of land surface temperature–vegetation abundance relationship for urban heat island studies. Remote Sens. Environ., 89 , 467–483.
Wu, X., , W. P. Menzel, , and G. S. Wade, 1999: Estimation of sea surface temperatures using GOES-8/9 radiance measurements. Bull. Amer. Meteor. Soc., 80 , 1127–1138.
Yu, Y., , J. L. Privette, , and A. Pinheiro, 2008: Evaluation of split-window land surface temperature algorithms for generating climate data records. IEEE Trans. Geosci. Remote Sens., 46 , 179–192.
Locations of the SURFRAD and ARM observation sites.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
LST retrieval bias error for different algorithms: (a) Price, (b) Prata and Platt, and (c) BL95.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
As in Fig. 2, but for (a) Ulivieri1, (b) Ulivieri2, and (c) Vidal.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
As in Fig. 2, but for (a) Sobrino1, (b) Sobrino2, and (c) Coll.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
As in Fig. 2, but for (a) BL90 local split-window, (b) WD96 generalized split-window, and (c) Sun and Pinker (2003) split-window algorithms.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
Evaluation (bias and RMS) of daytime LST retrievals (1/8°) from the GOES-8 with the SP-SW algorithm against the ARM observations for four different seasons.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
The RMS of nighttime LST retrievals for the NOAA-R, GSW, SP-SW, and BL95 algorithms at (a) 1/2° and (b) 1/8° grid evaluated against ARM SGP observations in 1997.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
The (a) bias and (b) residual standard errors for the NOAA-R, GSW, SP-SW, and BL95 algorithms evaluated at Bondville during daytime in winter (December–February) from 1996 to 2000.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
The seasonal (spring, summer, autumn, and winter) (a) bias and (b) residual standard errors for nighttime LST retrievals as evaluated at Bondville.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
The seasonal (spring, summer, autumn, and winter) LST bias during (a) daytime and (b) nighttime at Fort Peck from 1996 to 2000.
Citation: Journal of Applied Meteorology and Climatology 48, 1; 10.1175/2008JAMC1781.1
List of selected algorithms and coefficients.
ARM and SURFRAD ground station information.
The LST retrieval errors (1/8°) evaluated against ARM SGP data.
Nighttime LST retrieval errors (1/8°) evaluated against ARM SGP data.
The LST retrieval errors (1/8°) evaluated against the Mesonet observations.
Comparison between the retrieved LST (1/8°) and in situ measured LSTs in MODIS validation sites conducted in 2000; the in situ LSTs are taken from Wan et al. (2002). The time/date information is in Pacific standard time.
