Retrieval of Land Surface Temperature over the Heihe River Basin Using HJ-1 B Thermal Infrared Data

The reliable estimation of spatially distributed Land Surface Temperature (LST) is useful for monitoring regional land surface heat fluxes. A single-channel method is developed to derive the LST over the Heihe River Basin in China using data from the infrared sensor (IRS) onboard the Chinese “Environmental and Disaster Monitoring and Forecasting with a Small Satellite Constellation” (HJ-1B for short for one of the satellites), with ancillary water vapor information from Moderate Resolution Imaging Spectroradiometer (MODIS) products (MOD05) and in situ automatic sun tracking photometer CE318 data for the first time. In situ LST data for the period from mid-June to mid-September 2012 were acquired from automatic meteorological stations (AMS) that are part of Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project. MOD05-based LST and CE318-based LST are compared with in situ measurements at 16 AMS sites with land cover types of vegetable, maize and orchards. The results show that the use of the MOD05 product could achieve a comparable accuracy in LST retrieval with that achieved using the CE318 data. The largest difference between the MOD05-based LST and CE318-based LST is 0.84 K throughout the study period over the Heihe River Basin. The standard deviation (STD), root mean square error (RMSE), and correlation coefficient (R) of HJ-1B/IRS vs. the in situ measurements are 2.45 K, 2.78 K, and 0.67, respectively, whereas those for the MODIS 1 km LST product vs. the in situ measurements are 4.07 K, 2.98 K, and 0.79, respectively. The spatial pattern of the HJ-1B/LST over the study area OPEN ACCESS Remote Sens. 2015, 7 301 in the Heihe River Basin generally agreed well with the MODIS 1 km LST product and contained more detailed spatial textures.


Introduction
Land Surface Temperature (LST), one of the most important key parameters in the water balance and energy exchange between the surface and atmosphere, is widely used in many applications [1][2][3][4].To obtain spatial and temporal records of LST, retrievals from satellite observations, such as geostationary and polar-orbit satellites, are of fundamental importance in global and regional studies [5,6].Due to technical constraints, it is difficult to acquire remotely sensed LST with high spatial, temporal and spectral resolutions simultaneously from a single sensor.For example, although the new Landsat 8 Thermal Infrared Sensor (TIRS) provides two thermal infrared bands, the revisit frequency of 16 days makes it difficult to apply in environmental studies that require more frequent information [6].Low spectral information, which usually has only one thermal band for sensors with high spatial resolution, impedes the accurate atmospheric correction of satellite data because single-channel sensors observe much less information on atmospheric conditions than multi-channel sensors [7][8][9][10].
The Chinese "Environmental and Disaster Monitoring and Forecasting with a Small Satellite Constellation" (HJ-1B for short for one of the satellites) have been widely used in various studies [11,12], including those that use agricultural production analyses, vegetation indexes and evapotranspiration estimates.In addition, HJ-1B satellite could also be used in Urban Heat Island (UHI) studies and fragile ecosystem studies due to its relatively high temporal frequency and high spatial resolution, which can provide global coverage every two days with a spatial resolution of 300 m.However, the HJ-1B satellite has only one infrared sensor (IRS), and only single-channel methods can be used in LST estimation.
Many studies of single-channel methods have been published.To solve the atmospheric correction problems associated with single-channel sensors, Qin and Karnieli [7] developed a mono-window algorithm for LST retrieval from Landsat TM6 data.Jiménez-Muñoz and Sobrino [13] developed a generalized single-channel method and applied to the Landsat Thematic Mapper (TM) for LST retrieval via Water Vapor Content (WVC) rather than atmospheric vertical profiles of air temperature and relative humidity obtained by radiosonde.Sobrino and Jiménez-Muñoz [14] also noted that the single-channel method provides similar or better results than the split-window for low atmospheric WVC conditions.The accuracy of atmospheric WVC is crucial because of its large impacts on atmospheric downwelling and upwelling radiances.Li [15] developed a single-channel algorithm based on a parametric model, which is more suitable for HJ-1B/IRS due to its better accuracy.Zhao et al. [16] performed an error analysis of LST retrieval from HJ-1B/IRS data and confirmed that the LST retrieved from HJ-1B/IRS could not be higher than 1 K unless accurate atmospheric profile and land surface emissivity (LSE) data were applied.Coll et al. [17] used different sources of atmospheric profiles to retrieve LST with the single-channel method, and the result indicated that the retrieved LSTs based on radiosonde atmospheric profiles were better than that based on re-analysis profiles and satellite sounder products.Yet it is still difficult to obtain radiosonde data routinely.Zhou et al. [18] improved a mono-window method for HJ-1B/IRS retrievals by considering the viewing zenith angle, yielding an accuracy better than 1 K.
Although studies have been performed on HJ-1B LST retrieval, further studies of HJ-1B LST retrieval are still needed, especially in semi-arid regions where accurate land surface emissivity estimates are usually difficult to obtain.The objective of this paper is to develop and validate a single-channel algorithm by using different ancillary WVC data and single-channel thermal infrared data from the HJ-1B/IRS with a spatial resolution of 300 m in semi-arid regions in Heihe River Basin.
Section 2 describes the principle of the single-channel method developed for HJ-1B/IRS observations.Section 3 introduces the data used in the present study.Section 4 presents the comparison of the results with in situ measurements.Finally, the conclusions are presented in Section 5.

Methodology
When the sky is clear and a local thermodynamic equilibrium exists in the troposphere, the LST can be calculated from the Top Of Atmosphere (TOA) Brightness Temperature (BT) observations on the basis of the radiative transfer theory as follows Equation ( 1) [1]: where , εi is the land surface emissivity, and τi is the total atmospheric transmittance along the path between the land surface and the satellite sensor.
The LSE can be estimated either based on the vegetation cover method [19], or on the classification-based method [15,20,21].The effective emissivity of the land surface is calculated as follows Equation (2) [4]: (1 ) where εv is the vegetation emissivity, εg is the bare soil emissivity, and FVC is the fractional vegetation cover of a pixel and is calculated from the MODIS Normalized Difference Vegetation Index (NDVI) product (MOD13Q1) as follows Equation (3) [4]: where NDVIv and NDVIs are the NDVI values of fully covered vegetation and bare soil, respectively.In this paper, NDVIv and NDVIs are taken as the maximum and minimum values of the image, respectively.Directional hemispherical reflectance values for different land cover types are provided by the Johns Hopkins University (JHU) spectral library, Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) spectral library and MODIS spectral library (University of California, Santa Barbara (UCSB)) [22].The channel land surface emissivities of HJ-1B/IRS are calculated by convolving the Spectral Response Function (SRF) with the spectral library emissivities as follows Equation (4) [15]: where f(λ) is the HJ-1B/IRS SRF (Figure 1), λ1 to λ2 is the wavenumber range in the SRF, and λ is the wavenumber.Another essential procedure in the LST retrieval is to obtain the atmospheric parameters (including the atmospheric transmittance, atmospheric upward radiation and atmospheric downward radiation) for atmospheric correction.In the simulation using atmospheric transfer model MODTRAN 4.0 to acquire these parameters, Cloudless Land Atmosphere Radiosounding (CLAR), including 382 equally distributed global atmospheric profiles, is used in the method proposed by Li [15].The relationship between the simulated atmospheric parameters (transmittance τ, upwelling atmospheric radiance L↑and downwelling atmospheric radiance L↓) and water vapor content (w) are established by statistical regression from the simulation above.
However, Li [15] notes that the quadratic relationship between atmospheric parameters and WVC proposed by Jiménez-Muñoz and Sobrino [13] could lead major errors when the WVC is lower than 0.5 cm.Therefore, these parameters are derived as the cubic function of the WVC in our study as shown in Equations ( 5) to ( 7) [15].The regression correlation coefficient (R) for each equation is 0.995, 0.990 and 0.992, respectively.Consequently, the atmospheric parameters could also be obtained with WVC when detailed relatively high accuracy atmospheric profiles are not available.

HJ-1B/IRS Data
The "Environmental and disaster monitoring and forecasting with a small satellite constellation" [23] comprises two small optical satellites (HJ-1A and HJ-1B) that were launched in 2008.The sensor of HJ-1B/IRS consists of four bands, i.e., Near Infrared band (NIR, 0.75-1.10microns), Shortwave Infrared band (SWIR, 1.55-1.75microns), Middlewave Infrared band (MWIR, 3.50-3.90microns) and Longwave Infrared band (LWIR, 10.5-12.5 microns).The spatial resolution of the NIR, SWIR and MWIR bands is 150 m and that of the LWIR band is 300 m.HJ-1B/IRS can provide global coverage every two days, and the revisit frequency is between the values of the MODIS and Landsat TM frequencies.The SRF of LWIR of HJ-1B/IRS is shown in Figure 1.Although HJ-1B/CCD suffers from a lens anomaly due to sensor design and atmospheric conditions [24], no lens anomalies have been found in HJ-1B/IRS data yet.
The HJ-1B/IRS data used in this paper cover 13 cloud-free days in 2012: 19 and 30 June; 8, 23 and 27 July; 14, 15, 22, and 29 August; and 2, 6, 14 and 17 September.The results are then compared with the ground data.

Ground Data
The ground data are from the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) experiment in 2012.The experiment area is located in the middle reaches of the Heihe River Basin, a typical inland river basin with a semi-arid climate in Northwest China (Figure 2).The land surface is heterogeneous and is dominated by maize, vegetables, orchards, and residential areas.The spatial distribution of in situ Water Vapor Content measurements sites is also shown in Figure 2, including Wuxing Village (including Building 7 and Gouqu) and the Gaoya hydrological station.The details are shown on Section 3.3.2.Additional details can be found in Li et al. [25] and Liu et al. [26].
The in situ LST is derived from the land surface (Upwelling) Longwave Radiation (Rlu) and atmospheric (Downwelling) Longwave Radiation (Rld), which were observed using four-component radiation measurements.For the radiometer intercomparison in HiWATER, the four-component radiation measurements exhibit good agreement.The average RMSE and Mean Relative Error (MRE) are 6.43 W•m −2 and 0.62%, respectively, for Rld and 3.16 and −0.34%, respectively, for Rlu [25,27].The in situ measurements are stored in 10-min intervals, and the cloud-free measurements concurrent with the satellite overpasses (approximately 3:40 UTC, 11:40 Beijing Time both for HJ-1B and MODIS) are used for comparison with the retrieved LST.
Because of the large spatial variations, daytime in situ measurement data used in previous validation studies are scarce [2,28,29].Although the spatial and temporal variability of the LST is smaller at nighttime, satellite-derived LST products should also be validated during the daytime because the LST data are mostly required during the day, e.g., for monitoring surface heat fluxes [30].For these reasons, only the daytime values at the homogeneous sites (14 maize sites, 1 vegetable site, and 1 orchard site) were used in this study.

MODIS Data
MODIS, which is on board the Terra and Aqua satellites, has 36 spectral channels ranging from 0.405 to 14.385 μm.The MODIS data products used in this study are given in Table 1 and, are downloaded from "Reverb", a next generation metadata and service discovery tool [31].It is worth noting that the MOD05 product is affected by the MODIS band 5 stripe noise.In this study, the noise is removed using median filtering with a 9 × 9 pixel filter kernel.Moreover, the International Geosphere-Biosphere Program (IGBP) global vegetation classification from the land cover product MCD12Q1 is used.The in situ automatic sun tracking photometer CE318, which is provided by HiWATER from 19 June to 15 September 2012, includes two detectors that cover the entire spectrum from 340 to 1640 nm with eight bands.The detectors determine atmospheric parameters, such as aerosols and WVC (retrieved from the measurements at 936 nm).Other than the ground data, only 26 cloud-free samples of the CE318 data were acquired in this study.These data are obtained at the three sites in Zhangye and one site in Arou; these areas are shown in Figure 2. In Zhangye, the sites include Wuxing Village (Building 7), Wuxing Village (Gouqu) and the Gaoya hydrological station.It is worth noting that the CE318 data acquired at Building 7 and at Gouqu are from different days.As shown in Table 2, the data on June 19 are obtained at Building 7 and the other data are obtained at Gouqu.

Comparisons between MOD05 and In Situ Measurements of Water Vapor Content
The determination of atmospheric WVC is of vital importance in the retrieval of LST from satellite data.Although the usage of in situ WVC measurements can provide more accurate results than other satellite-based WVC data, the in situ WVC measurements are acquired for many reasons, such as human resources, time and other irresistible things.Therefore, the satellite-based WVC data are more economic and flexible.In this paper, the performance of the WVC from MODIS (MOD05), which are used for the HJ-1B/IRS atmospheric correction, is first assessed with the WVC from the in situ automatic sun-tracking photometer CE318 (referred to CE318-WVC).
The results from the cloud-free observations (Figure 3a) demonstrate that the values of the RMSE, MRE, and the determinant coefficient R 2 between the MOD05-WVC and CE318-WVC are 0.34 cm, 18%, and 0.93, respectively.Figure 3b clearly shows that the MOD05-WVC is slightly overestimated compared with the in situ CE318-WVC.When the WVC is approximately less than 1.5 cm (which usually occurs in semi-arid area), the WVC differences between the MOD05 product and the CE318 data are lower than 0.4 cm, which is acceptable in the LST retrievals.This result makes MOD05-WVC data reasonable for LST estimation in semi-arid area of the Heihe River Basin.
Among all the CE318 sites, Wuxing Village (Figure 2) provides the relatively complete observations: both CE318 data and in situ LST data are obtained here, and more cloud-free days occur.In this study, Table 2 presents the retrieved LST using the two WVC datasets for a maize site (AMS 5) in Wuxing Village.The largest LST difference between the MOD05-based LST and CE318-based LST appeared on June 19 as 0.84 K (with a WVC difference of 0.38 cm).The large error may be caused by the different atmospheric conditions at Building 7. The RMSE and R 2 between the MOD05-WVC and CE318-WVC in Wuxing Village are 0.12 cm and 0.98, respectively.Meanwhile, the RMSE and R 2 between the MOD05-based LST and in situ LST are 1.27 K and 0.53, respectively, and those between the CE318-based LST and in situ LST are 1.35 K and 0.71, respectively.The results demonstrate that use of the MOD05 product ensures a LST retrieval accuracy that is comparable to that produced when using the CE318 data.

HiWATER-MUSOEXE Results
Statistics of the comparisons between the HJ-1B/IRS LST, MODIS 1 km LST, and in situ AMS measurements for each of the HiWATER sites shown in Figure 2 are summarized in Table 3.Although there are fewer days with HJ-1B/IRS data than with MODIS data, both comparisons with in situ data yield similar results: the HJ-1B/IRS standard deviation (STD), RMSE, and linear correlation coefficient (R) values of HJ-1B/IRS vs. AMS are 2.45 K, 2.78 K, and 0.67 (with 204 samples), respectively, whereas those for the MODIS 1 km LST values are 4.07 K, 2.98 K, and 0.79 (with 886 samples), respectively.
A single Taylor [32] diagram is shown in Figure 4 to summarize multiple aspects of the method performance.The labels are assigned to six cases of the available data listed in Table 3.The position of each label on the plot quantifies how closely the estimation matches the in situ measurements.The RMSE (blue dotted semi-circle) between the estimated and in situ measurements is the difference from the "observed" values, which are defined as the in situ measurement STDs.The labels with a small distance the "observed" points (i.e., smaller values of "observed" contours) correspond to the estimates that agree well with the in situ measurements.Most of these estimates have relatively high correlations and low RMSEs.The estimates that lie on or near the dashed arc have the correct STD.In Figure 4, it can be observed that label "a" (vegetables, HJ-1B/IRS) generally agrees better with the in situ measurements and exhibits lower RMSE and higher correlation.For the other labels (maize and orchards), the results are not as good, as the RMSEs are larger than 1.5 K (outside the dotted semi-circle of 1.50) and the correlations are lower than 0.9 (outside the green line of 0.9).The vegetable site (AMS 1, labeled "a" in Figure 4) is a relatively flat area with an elevation of 1552.75 m. Figure 5a indicates good agreement between the site observations and the satellite estimates at the vegetable site, which are also shown in Figure 4 and Table 3.
The 14 maize sites (AMS 2-AMS 3, AMS 5-AMS 16) are typical maize croplands in the HiWATER project.These sites are analyzed as a group in this study.The RMSE of 2.83 K in the maize area is acceptable due to the intercropping of certain croplands.The MODIS 1 km LST produces significant underestimations (the red ellipse in Figure 5b).These underestimations may partly be caused by the effect of clouds.It is difficult to efficiently exclude the clouds, especially in the partly cloudy pixels.
The orchard site (AMS 17) is planted with apple trees and intercropping maize.Unlike the vegetable and maize sites, the landscape at this site is heterogeneous.However, the results at the orchard site are similar to those at the maize sites (Table 4 and Figure 4).This result is confirmed in Figure 5c. Figure 5c also shows that for the pixels of low temperatures, the estimation by HJ-1B/LST is much more consistent with the in situ values.

Spatial Distribution
Because the purpose of this study is to retrieve the spatial distribution of LST that can be used for monitoring regional surface heat fluxes, it is important to assess the spatial distribution of the retrieved LST from the HJ-1B data.Figure 6 shows examples of the spatial distribution of the HJ-1B/IRS and MODIS 1 km LST.Both satellite images were acquired on 22 August 2012, and the MODIS 1 km LST data were resampled to 300 m.The HJ-1B/IRS data were acquired at approximately 3:50 UTC, Figure 7 shows the input parameters for the HJ-1B/IRS LST retrieval.The TOA radiance is between 4 and 11 W•m −2 •sr −1 •μm −1 .In the solid square (the irrigated oasis), the TOA radiance is approximately 9 W•m −2 •sr −1 •μm −1 when including the 16 AMS sites.The LSE is between 0.96 and 0.99, and the mean LSE is approximately 0.99.The WVC is approximately 1.5 cm.The WVC of the irrigated oasis is higher than that in the mountainous areas (the DEM shown in Figure 7d is higher than 4200 m) and desert (or sparsely vegetated in Figure 7e) because in the irrigated oasis, the irrigation provided more moisture.As mentioned (Figure 3b), when the WVC is approximately 1.5 cm, the WVC differences between the MOD05 product and the CE318 data are lower than 0.4 cm; thus, the LST differences between the MOD05-based LST and CE318-based LST are less than 1 K. Figure 7d shows the Digital Elevation Data (DEM) obtained from the Shuttle Radar Topography Mission (SRTM) of the study area.The irrigated oasis is a nearly flat plain.Figure 7e shows the IGBP land cover type from MCD12Q1 product.These data are also the source data for LSE production.It should be mentioned that the spatial resolution of land cover data from MCD12Q1 is 500 m, which may have certain scaling effects on the LSE production that is used as a spatial resolution of 300 m in LST retrieval.
The spatial patterns of the retrieved HJ-1B/IRS LST over part of the middle stream and upstream regions of the Heihe River Basin generally agree well with the MODIS 1 km LST product (Figure 6).Generally, the HJ-1B/IRS LST product (Figure 6a) provides a more detailed texture for the study area (and other complex landscapes) than the MODIS product (Figure 6b) because of the higher spatial resolution of the HJ-1B data.Although the cloud mask (pixels in white color) in the MODIS 1 km product is not used in the HJ-1B/IRS product, the actual low temperatures could be discriminated from the image, such as the dark blue color in the image of Figure 6a.In the solid square, the mean LST is approximately 300 K for both images.Figure 6c shows the density comparison between HJ-1B/IRS and MODIS derived LST for the 10,000 pixels in the solid square.In Figure 6c, the cloud mask data (80 pixels in the solid square) in MODIS are excluded.To analyze the difference between HJ-1B/IRS and MODIS derived LST, the difference image between them is also obtained in Figure 6d.In the solid square and other flat and homogeneous places (cropland and grassland in Figure 7e), the differences between the two sensors derived LST are clearly shown to be within 1 K.The large differences occur both in mountainous areas and desert in the southwest and north of the image, respectively.On 16 May 2012, the plants were seeded in the Zhangye oasis (in the solid square) when the landscape was covered by sparse vegetation.Therefore, the LST image on this day is not as smooth as that of the fully covered vegetation days, such as 19 June 2012 and 2 September 2012.However, the difference between the results from two sensors is not as obvious.Furthermore, the large difference shown in Figure 8b may be heavily affected by the cloud on the retrieval images from two sensors.
In the solid square, the view zenith angles (VZAs) of HJ-1B/IRS are all approximately 30 degrees.Meanwhile, the VZAs of MODIS vary from 10 degrees to 50 degrees in the same area.The VZA comparisons are shown in Table 4. Apparently, the minimum VZA differences between MODIS and HJ-1B/IRS occurred on 22 August 2012, in which the LST differences (Figure 7c) of MODIS and HJ-1B/IRS are also smaller than those on other days (Figure 8).Wan et al. [2] demonstrated that different VZAs between different sources of measurements may result in different land surface emissivity values, as well as different LST retrieval results.In this study, although the VZA differences between HJ-1B and MODIS are expected to be small for fully covered croplands [17,33], the differences are still an important source for bias in the LST retrievals and impede the validation between different sensors.Additionally, low spectral information impedes accurate atmospheric correction of satellite data because single-channel sensors observe much less information on atmospheric conditions than multi-channel sensors.The proposed single-channel method in this paper is mainly suitable in semi-arid regions with low atmospheric WVC.

Conclusions
The reliable estimation of spatially distributed Land Surface Temperature (LST) is particularly crucial for monitoring regional surface heat fluxes.For the first time, a feasible single-channel algorithm has been developed to derive spatially distributed LST from the infrared sensor (IRS) onboard the Chinese "Environmental (HuanJing in Chinese) and disaster monitoring and forecasting with a small satellite constellation" (HJ-1B) with different ancillary data, i.e., the MODIS water vapor content (WVC) product and in situ automatic sun tracking photometer CE318 data.The results show that the retrieved HJ-1B LST data, through use of the MOD05 WVC product, can achieve an accuracy comparable to that of the CE318 data for WVC.The largest LST difference between the MOD05-based LST and CE318-based LST was 0.84 K (associated with a WVC difference of 0.38 cm).
Additionally, the Standard Deviation (STD), Root Mean Square Error (RMSE), and correlation coefficient (R) of HJ-1B/IRS vs. the in situ measurements were 2.45 K, 2.78 K, and 0.67, whereas those for the MODIS 1 km LST vs. the in situ measurements were 4.07 K, 2.98 K, and 0.79, respectively.
The spatial pattern of the HJ-1B/LST over the study area in the Heihe River Basin generally agreed well with the MODIS 1 km LST product and contained more detailed spatial textures.In this study, although the view zenith angle (VZA) difference between HJ-1B and MODIS are expected to be small for fully covered croplands, these differences are still important source of bias in LST retrievals and the validation between different sensors.Low spectral information impedes the accurate atmospheric correction of satellite data because single-channel sensors observes much less information on atmospheric conditions than multi-channel sensors.The proposed single-channel method in this paper is mainly suitable for semi-arid regions with low atmospheric WVC.
To improve the validation process, future work should be dedicated to conducting careful and quantitative validations between different remotely sensed estimations and should eliminate the sensor attribute differences.Additionally, more surface sensors should be deployed to obtain in situ measurements during a field campaign.

Figure 2 .
Figure 2. Location of the study area and distribution of the ground sites (including in situ LST measurements sites and WVC sites) used in this study.The location of the HiWATER-MUSOEXE area is indicated with the small red quadrangle (in the left image) and a solid black quadrangle (in the right image).The black square in the left image is the Heihe Matrix study area, which will be fully analyzed in the LST spatial distribution analyses.

Figure 3 .
Figure 3. Scatter plots of (a) the MOD05-WVC and in situ CE318-WVC for cloud-free observations; (b) the difference between the MOD05-WVC and in situ CE318-WVC compared with the CE318-WVC at the HiWATER sites from 19 June to 15 September 2012.The number of valid samples is 26.

Figure 4 .
Figure 4. Taylor diagram that describes the performance of the LST retrieval from HJ-1B/IRS and MODIS over the HiWATER experiment sites in 2012.

( 11 :Figure 6 .
Figure 6.(a) The retrieved LST images with a 300 m resolution on 22 August 2012 (3:50 UTC, 11:50 Beijing Time); (b) the MODIS LST (3:30 UTC, 11:30 Beijing Time), where the white color is the MODIS cloud mask; (c) the density of the HJ-1B/IRS and MODIS derived LST for the 10000 pixels in the solid square; and (d) the difference between HJ-1B/IRS and MODIS derived LST.

Figure 8
shows other comparison examples between the HJ-1B/IRS and MODIS derived LST in the growing season: (a) 16 May 2012 (in situ data were not available on this day in the former discussion); (b) 19 June 2012; and (c) 2 September 2012.The result with the obvious color changes indicates that the LST greatly changes during the growing season.

Table 1 .
Specifications of the MODIS products used in this study.

Table 2 .
Comparison between MOD05-based LST, CE318-based LST and in situ LST at HiWATER AMS 5 in Wuxing Village.

Table 3 .
Comparison of the HJ-1B/IRS and MODIS 1 km LST algorithm performance at the HiWATER sites (from June to September 2012).

Table 4 .
Comparison of the HJ-1B/IRS and MODIS view zenith angles (VZAs) in the solid square from Figure 6 and Figure 8 in 2012.