Validation and Analysis of Long-Term AATSR Land Surface Temperature Product in the Heihe River Basin, China

The Advanced Along-Track Scanning Radiometer (AATSR) land surface temperature (LST) product has a long-term time series of data from 20 May 2002 to 8 April 2012 and is a crucial dataset for global change studies. Accuracy and uncertainty assessment of satellite derived LST is important for its use in studying land–surface–atmosphere interactions. However, the validation of AATSR-derived LST products is scarce in China, especially in arid and semi-arid areas. In this study, we evaluated the accuracy of the AATSR LST product using ground-based measurements from 2007 to 2011 in the Heihe River Basin (HRB), China. The AATSR-derived LST results over Yingke site are closer to ground measurements than those over A’rou site for both daytime and nighttime temperatures. For nighttime, the averaged bias, STD, RMSE and R2 over both sites are 0.67 K, 3.03 K, 3.13 K and 0.93 K, respectively. Based on the accuracy assessment, we analyzed the AATSR-derived annual LST variations both in the HRB region and the two validation sites for the period of 2003 to 2011. The results at the A’rou site show an obvious increasing trend for daytime from 2003 to 2011. For the whole HRB region, the warming trend is clearly shown in the downstream of HRB.


Study Area
The study area is located in the semiarid HRB of northwest China between the coordinates 97 • 24 ~102 • 10 E, 37 • 41 ~42 • 42 N. The simultaneous satellite and ground-based remote sensing WATER experiments were conducted in summer 2008, and the ground based measurements were acquired from 2007 to 2011. The semiarid basin is mainly controlled by the distribution of water resources from south to north, distributed with diverse land cover types from the upstream, midstream to downstream including glacier (the southern Qilian Mountains), frozen soil, alpine meadow, forest, irrigated crops, riparian ecosystem, desert and Gobi. The basin is influenced by Asian monsoon, and most rainfall occurs from May to September in the rainy season, and the rainfall decreases from south to north. Figure 1 shows the geography and land cover encompassing the HRB with the two validation sites as well [27].
Remote Sens. 2017, 9,152 3 of 18 the methods. Section 3 describes the corresponding results. A discussion and a conclusion are provided in Sections 4 and 5, respectively.

Study Area
The study area is located in the semiarid HRB of northwest China between the coordinates 97°24′~102°10′E, 37°41′~42°42′N. The simultaneous satellite and ground-based remote sensing WATER experiments were conducted in summer 2008, and the ground based measurements were acquired from 2007 to 2011. The semiarid basin is mainly controlled by the distribution of water resources from south to north, distributed with diverse land cover types from the upstream, midstream to downstream including glacier (the southern Qilian Mountains), frozen soil, alpine meadow, forest, irrigated crops, riparian ecosystem, desert and Gobi. The basin is influenced by Asian monsoon, and most rainfall occurs from May to September in the rainy season, and the rainfall decreases from south to north. Figure 1 shows the geography and land cover encompassing the HRB with the two validation sites as well [27].  [27]). The right image is the AATSR LST image on 23 December 2008.
As indicated by Wan et al. [25], the major requirements of an ideal LST validation site should be large enough to cover at least dozens of satellite pixels and homogeneous. There are ten sites in WATER campaign in HRB [27]. In order to guarantee the consistency between in situ observations with remote sensing data, homogeneous in situ sites should be chosen. In this study, relatively homogeneous ground measurements from two of the ten sites were chosen: A'rou and Yingke.
A'rou, is located in the source of Heihe River, a relatively homogeneous area with a height of 3032.8 m. The land cover type is Alpine meadow. A'rou is an important station on the northeast edge of the Tibetan Plateau with a lot of hydrological and ecological in situ observations, which could both benefit the future studies of HRB and Tibetan Plateau as well. It should be noted that the seasonally frozen soil is widely distributed in A'rou, which may hinder the wide application of satellite derived Figure 1. Geographical/Land cover map encompassing the HRB and the two in situ sites position (for details see [27]). The right image is the AATSR LST image on 23 December 2008.
As indicated by Wan et al. [25], the major requirements of an ideal LST validation site should be large enough to cover at least dozens of satellite pixels and homogeneous. There are ten sites in WATER campaign in HRB [27]. In order to guarantee the consistency between in situ observations with remote sensing data, homogeneous in situ sites should be chosen. In this study, relatively homogeneous ground measurements from two of the ten sites were chosen: A'rou and Yingke.
A'rou, is located in the source of Heihe River, a relatively homogeneous area with a height of 3032.8 m. The land cover type is Alpine meadow. A'rou is an important station on the northeast edge of the Tibetan Plateau with a lot of hydrological and ecological in situ observations, which could both benefit the future studies of HRB and Tibetan Plateau as well. It should be noted that the seasonally frozen soil is widely distributed in A'rou, which may hinder the wide application of satellite derived LST product due to the uncertainty of the algorithm in the freeze/thaw condition. The aims of the validation have also focused on the LST uncertainty effect under freeze/thaw status [27]. Yingke site is located at 38 • 51 N, 100 • 25 E with an elevation of 1519 m. The dominant land cover at Yingke site is homogeneous cropland which turns to bare soil in winter months [27]. The geographical encompass and site characteristics are shown in Figure 1 and Table 1. The land cover map in Figure 1 is acquired from the yearly MODIS Land Cover Type product (MCD12Q1) in 2009 with a spatial resolution of 500 m.

AATSR LST Product
AATSR onboard the Envisat satellite is the third sensor following ATSR-1 and ATSR-2 designed to obtain the high levels of Sea Surface Temperature (SST) and LST for global climate change studies. AATSR measures seven channels of reflected and emitted radiation at 0.55 µm, 0.66 µm, 0.87 µm, 1.6 µm, 3.7 µm, 11 µm and 12 µm with both a nadir view resolution of 1 × 1 km 2 and a forward view resolution of 1.5 × 2 km 2 . The AATSR swath width is 512 km and able to provide a global LST coverage about every three days [22,26]. AATSR-derived LST products have been widely used in climate change, land-atmosphere feedbacks, modeling studies, land cover changes and many other applications [17].
Only the cloud-free two channel brightness temperatures of 11 µm (T 11 ) and 12 µm (T 12 ) with nadir view dataset were used for the AATSR LST product retrievals to avoid the angular effects on temperature and emissivity [18,25]. The basic split-window algorithm was used for the LST retrieval with the following mathematics: where a 0 , b 0 and c 0 are regression coefficients determined by simulated data, which are dependent on the land surface conditions, including the land cover types (27 biomes), vegetation fraction, season, time of day (daytime or nighttime), as well as atmospheric conditions and satellite viewing zenith angle [18]. The RMSE of LST is less than 1 K in ideal conditions [22,23]. The AATSR LST data acquisition period is from 20 May 2002 to 8 April 2012, which is available from the AATSR GlobTemperature Level-2 LST product [26] provided by ESA DUE GlobTemperature Project [28]. The sensor is in a sun-synchronous orbit and the overpass time is 10 a.m. local solar time. The spatial resolution of AATSR LST product is 0.01 • in NetCDF (Network Common Data Form) format in this study. The AATSR LST product is estimated using the nadir view observations, and the viewing angle at nadir is less than 23.5 • [19,29]. The GlobTemperature Level-2 AATSR LST products include a primary LST dataset and an additional auxiliary dataset. The primary LST dataset is comprised by Julian Date, pixel center latitude and longitude, time difference from a reference time, LST, LST uncertainty, quality control flags, satellite zenith viewing angle and satellite azimuth angle. The auxiliary dataset includes channel description, emissivity, brightness temperature, land cover classification, fractional vegetation cover, water vapor content, normalized difference vegetation index, solar zenith angle, solar azimuth angle, land-water mask, etc. To obtain accurate at-sensor observation time, the dataset also includes the "dtime", i.e., the difference between acquisition time and product time stamp.

In Situ Measurements
To validate AATSR-derived LST product, the 10-min interval in situ quality controlled four-component radiation measurements at A'rou and Yingke were acquired from the automatic meteorological observations on the specific HRB experiments WATER [27,30,31] in the watershed from 2007 to 2011. The information on two validation sites is listed in Table 1. A'rou is chosen as a fully equipped homogeneous superstation among the HRB comprehensive hydrometeorological observation network. To be noted that there are more rainfall/cloudy days at A'rou than that at Yingke. The rainy weather in A'rou may affect the final validation results. Yingke site is located in large areas around the oasis which is much more homogeneous and suitable for LST validation.
The in-situ LST was derived from the in-situ upwelling (Rlu, W·m −2 ) and downwelling (Rld, w W·m −2 ) longwave radiation fluxes from the traditional equation: where ε is the broadband surface emissivity; and σ is the Stefan-Boltzman constant (5.67 × 10 −8 W·m −2 ·K −4 ). In this study, the broadband surface emissivity at each site was obtained using a broadband emissivity method from Cheng et al. [31] as the following equation: where ε 29 and ε 31 are the MODIS emissivities for channels 29 and 31, for the spectral domain of 8.400-8.700 µm and 10.780-11.280 µm, respectively. The broadband surface emissivity accuracy and stability will not heavily affect the ground-based LST calculation since the right side of the equation has a fourth root which makes the emissivity (0.8-0.99) magnitude affect less.
The in situ measurements comparison was conducted with cloud free AATSR-derived LST data with the quality control flags. Due to the different acquisition time of in situ data, the period of the validation process is from August 2007 to November 2011 at A'rou and November 2007 to November 2011 at Yingke, respectively.

Methods
This research includes two parts: First, we compared the AATSR-derived LST with in situ measurements over two sites in HRB for 2007-2011. This evaluation is intended to show the accuracy of the AATSR-derived LST and determine the feasibility of using it in studies involving temperature data. Second, based on the evaluation, we analyzed the time series LST of the validation sites and the whole HRB for the period of 2003-2011 and gave a preliminary conclusion for the LST dynamics.

Assessment of AATSR-Derived LST with in Situ Measurements over A'rou and Yingke for 2007-2011
The reliability of the long-term accuracy of AATSR-derived LST data was first evaluated. The corresponding in situ measurements both collocation in space and concurrent in time with the AATSR-derived LST data were chosen. The in situ data were collected every 10 min. The in-situ measurements which have nearest collection time matching with AATSR data were selected and used in the comparison. Nearly all the cloud-free satellite data from 2007 to 2011 available for the in-situ sites were used. It should be noted that the AATSR-derived LST has a spatial resolution of 0.01 • at nadir, while the in situ LST measurements were calculated from the four-component radiations at the local scale. The validation results, to some extent, might be affected by the different observation scales.
The assessment processes are shown as follows: First, the simple linear regression model was composed between AATSR-derived LST and in situ measurements. The linear relationship between them is shown in the scatterplot. In this regression model, the coefficient of determination, i.e., R-squared, which shows how well the data are explained by the best-fit line, is also presented as the standard regression output. This model composition was done at each site for both daytime and nighttime observations, respectively. In addition, the histograms of the differences between the satellite-derived LST and in situ measurements were also shown.
Next, three more error measures were calculated considering all the validation points for the two validation sites of A'rou and Yingke: bias (AATSR-derived LST minus in situ LST), STD (standard deviation) and RMSE (root-mean-square error). The annual and seasonal statistics between AATSR-derived LST and in situ LST were both computed.

Long-Term Variations of AATSR-Derived LST for 2003-2011
For the two validation sites, the annual cloud-free AATSR-derived LST variations at A'rou and Yingke site were first analyzed for 2003-2011. Then, the daytime, nighttime and overall time averaged cloud-free AATSR-derived LSTs were computed. In addition, the linear regression was also conducted to show the LST variation trend for nine years. The linear trend is often considered by climatologists rather than a polynomial trend, especially for long-term trend analysis.
For the whole HRB region, the cloud-free statistical analysis includes three parts: First, the nine-year averaged daytime, nighttime and the difference between daytime and nighttime LSTs were shown. Second, the AATSR-derived annually maximum, minimum and averaged LSTs for each pixel of the HRB were calculated for both daytime and nighttime from 2003 to 2011. Third, the corresponding variation trends of the maximum, minimum and averaged LSTs were presented to express the LST changes over this nine years period. In this paper, the trend is presented by the slope calculated from the simple linear regression model between the maximum (minimum and averaged) LST with the period.

Annual Analysis
The daytime and nighttime cloud-free validation results between AATSR-derived LST and in situ LST over A'rou and Yingke site for 2007-2011 were analyzed and summarized in Table 2. The scatterplots of the comparisons between the averaged AATSR-derived and ground based LSTs are given in Figure 2, as well as the difference histograms. The annual phase comparisons are shown in Figure 3.
The overall number of cloud-free satellite and in situ concurrent observations for 2007-2011 are 461 at A'rou and 477 at Yingke, respectively. The numbers of cloud-free nighttime observations (244, 250) are greater than that of daytime observations (217, 227) due to the stable weather conditions at night. GlobTemperature Project quality control information was used for the cloud mask in this research. The detail results at each site are discussed as follows.
During the daytime, the bias, STD, RMSE and R 2 over A'rou site are 1.57 K, 5.60 K, 5.80 K and 0.82 K for 2007-2011, respectively. At Yingke site, the corresponding statistics are 1.77 K, 3.12 K, 3.58 K and 0.96 K, respectively. During the nighttime, the AATSR-derived LST performs better: at A'rou site, 1.19 K, 3.32 K, 3.52 K and 0.90 K; and, at Yingke site, 0.15 K, 2.74 K, 2.74 K and 0.96 K. The daytime results at A'rou site are not good for at least two potential reasons: first, the A'rou site is located in the Qilian Mountains, and the mountain topography and windy weather condition will affect the satellite LST retrieval; second, the A'rou site has a freeze/thaw condition, especially in the daytime winter/spring when the freeze/thaw transition status happens frequently. When this transition happens, the LSE of mixed ice, water and soil are more difficult to be determined.
For daytime, the averaged bias, STD, RMSE and R 2 over both sites is 1.67 K, 4.36 K, 4.69 K and 0.89 K for 2007-2011, respectively. For nighttime, the averaged bias, STD, RMSE and R 2 over both sites is 0.67 K, 3.03 K, 3.13 K and 0.93 K for 2007-2011, respectively. The nighttime atmospheric water vapor is less than that in daytime. In addition, when at night the surface will not be heated by direct solar radiation resulted in the minimized thermal contrast between sunlit spots and shadows [32].
Accordingly, AATSR-derived nighttime LST performs better than the daytime LST for both validation sites in Table 2 and Figure 2. For the nighttime comparisons, the AATSR-derived LST at Yingke site with relatively more homogeneous cropland surface performs better than that at A'rou site.    For most sensors, clearly degradation has been shown during the satellite lifetime. In this study, the annual comparisons between AATSR-derived LST and in situ LST for 2008-2011 were also conducted, and the corresponding statistics are shown in Table 2 (Figure 3d), the AATSR-derived LST meets perfectly with ground based measurements within the whole study period.

Seasonal Analysis
The seasonal validations are shown in Table 3. The results suggest that in summer the AATSR-derived LST leads to the largest bias over both validation sites. There are two possible reasons for this. First, the number of co-located data pairs in summer is fewer than those in other seasons, except for the daytime comparison in Yingke site. The small sample size is more vulnerable for abnormal values. Second, in summer, the cloudy sky is much more common than in other seasons. The cloud contamination would also cause more bias in the retrieval results even though researchers have tried the best to mask clouds.
For the daytime comparisons at A'rou site, the error statistics in winter and autumn are smaller than those in spring and summer. In addition to the mentioned cloudiness problem, the LST has shown more dynamic changes in spring and summer than in winter and autumn due to the not equally distributed solar radiation. When the winter and autumn comparisons are analyzed, the autumn LST data give better matching with in situ measurements. As mentioned before, in winter some of the soils over the two sites are in the freeze/thaw transition conditions, especially over A'rou site.
For the nighttime bias comparisons at Yingke site, which are the best results in the seasonal comparison, the AATSR-derived LSTs were overestimated in summer (1.60 K) and autumn (1.17 K) and underestimated in winter (−1.06 K) and spring (−0.65 K). These results could strengthen the seasonal variation cycle conclusion made by Noyes et al. [21].
In conclusion, the AATSR-derived LST results over homogeneous Yingke site are closer to ground measurements than those over A'rou site both for daytime and nighttime. The discrepancies over the two sites are higher for daytime than nighttime. The cloud mask uncertainties, land surface emissivity instabilities associated with land cover types, freeze/thaw transition conditions and acquisition times are the major reasons for the differences in validation results.

Two Sites Analysis
The long-term time series of AATSR-derived LST at each site providing annual and seasonal variations are shown in Figure 4. The low values in each figure might be cloud contaminated pixels not screened by the quality flags. As shown in the figures, due to the different elevations, AATSR-derived LST at Yingke site is higher than that at A'rou site, especially at nighttime. For the daytime results at Yingke site, the highest temperature is around 1 June 2008, both in and Figures 3c and 4c during the nine years AATSR-derived LST data.    The annual averaged cloud-free AATSR-derived LSTs were then computed for A'rou and Yingke over 2003-2011 for daytime, nighttime and overall time, respectively (shown in Figure 5). Table 4 shows the number of available cloud-free AATSR data for each site used in Figure 5. For daytime, the linear regression model at A'rou shows a warming trend with a slope of 0.41 K/year for the AATSR-derived LST over the nine years, while at Yingke the model shows a cooling trend with a slope of −0.11 K/year. For nighttime, at A'rou and Yingke the linear regression slopes are −0.05 K/year (cooling) and 0.17 K/year (warming), respectively. Although the temporal variations of AATSR-derived nighttime averaged LST in A'rou show a little decrease, the temporal variations of AATSR-derived daytime (overall) averaged LST at A'rou site show an increase of nearly 4 K (2 K) in cloud-free conditions. Meanwhile, the warming trend at Yingke is not obvious due to the low value for daytime of 2011.
Remote Sens. 2017, 9, x FOR PEER REVIEW 12 of 18 The annual averaged cloud-free AATSR-derived LSTs were then computed for A'rou and Yingke over 2003-2011 for daytime, nighttime and overall time, respectively (shown in Figure 5). Table 4 shows the number of available cloud-free AATSR data for each site used in Figure 5. For daytime, the linear regression model at A'rou shows a warming trend with a slope of 0.41 K/year for the AATSR-derived LST over the nine years, while at Yingke the model shows a cooling trend with a slope of −0.11 K/year. For nighttime, at A'rou and Yingke the linear regression slopes are − 0.05 K/year (cooling) and 0.17 K/year (warming), respectively. Although the temporal variations of AATSR-derived nighttime averaged LST in A'rou show a little decrease, the temporal variations of AATSR-derived daytime (overall) averaged LST at A'rou site show an increase of nearly 4 K (2 K) in cloud-free conditions. Meanwhile, the warming trend at Yingke is not obvious due to the low value for daytime of 2011. Figure 6 shows the HRB nine-year (2003-2011) AATSR-derived LST average. The nine-year averaged LST shows clearly the mountain (grassland, blue color on Figure 6a, upstream HRB), the lake (blue color on Figure 6c) and the flat area (barren, red color, downstream HRB). In addition, the Heihe River contour and the around cropland could be easily obtained in the middle reaches. Figure  6c shows clearly the averaged LST difference of water with the barren area. Figure 7 shows the trend of the maximum, minimum and averaged AATSR-derived LST dynamics for 2003-2011. The averaged daytime (Figure 7c) and nighttime (Figure 7f) LST trends are generally mild except for a warmer trend around the upstream (yellow color on the southeast of the image) of HRB for daytime averaged LST. It should be pointed out that, in the downstream (north) of the HRB region, there are some lakes which were dried up before 2004. After 2004, the government makes the HRB water diversion to these areas. The efforts have significant effects on the lake environment which are shown clearly in Figure 7c,f: the daytime averaged LST trend is cooling (the blue color on the north of Figure 7c), while the nighttime averaged LST trend is remarkable warming  Table 4. The number of available cloud-free AATSR data for each site used in Figure 5.  2003  48  56  104  60  64  124  2004  54  61  115  61  64  125  2005  62  56  118  66  56  122  2006  50  51  101  58  63  121  2007  53  64  117  62  67  129  2008  47  67  114  63  69  132  2009  52  47  99  61  58  119  2010  53  65  118  58  62  120  2011  59  53  112  56  70  126 3.2.2. HRB Region Analysis Figure 6 shows the HRB nine-year (2003-2011) AATSR-derived LST average. The nine-year averaged LST shows clearly the mountain (grassland, blue color on Figure 6a, upstream HRB), the lake (blue color on Figure 6c) and the flat area (barren, red color, downstream HRB). In addition, the Heihe River contour and the around cropland could be easily obtained in the middle reaches. Figure 6c shows clearly the averaged LST difference of water with the barren area.     (Figure 7f) LST trends are generally mild except for a warmer trend around the upstream (yellow color on the southeast of the image) of HRB for daytime averaged LST. It should be pointed out that, in the downstream (north) of the HRB region, there are some lakes which were dried up before 2004. After 2004, the government makes the HRB water diversion to these areas. The efforts have significant effects on the lake environment which are shown clearly in Figure 7c,f: the daytime averaged LST trend is cooling (the blue color on the north of Figure 7c), while the nighttime averaged LST trend is remarkable warming (the red color on the north of Figure 7f).

HRB Region Analysis
The most standout result is implied in Figure 7b from the trend of the daytime minimum LST which shows a large number of LST decrease in the midstream of HRB region, while the downstream of HRB region shows an intensively increasing trend during the nine years. The increasing trend is clearly shown in the downstream of HRB by the nine-year statistics. Table 5 shows the region-averaged values of annual AATSR-derived LST departure from the nine-year averaged LST in the HRB. The anomaly of averaged values is tagged red in 2005 and 2010 when the region-averaged anomaly is larger than that of other years. These results should also be obtained in Figure 7. However, for the whole HRB region, the LST trend could hardly be concluded for 2003-2011 because of some unexpected anomalies which might be due to the moderate La Niña (cool) events during 2010-2011 [33]. In addition, the results are limited by satellite observation uncertainties, if the climate anomalies occur near the beginning or the end of the observations, the trend estimates may be biased [34,35]. Table 5. Averaged values of annual AATSR-derived LST departure from the nine-year (2003-2011) average shown as Figure 6 in the HRB.

Satellite Derived LST Validation
The AATSR-derived LST product has a wide variety of future applications, especially with its ancestor ATSR-2 and its successor SLSTR (Sea and Land Surface Temperature Radiometer) onboard Sentinel-3 satellite. Combing these sensors derived LST products together, twenty years globally data could be obtained and could give a substantial LST data for studies of climate change, land-atmosphere feedbacks, land cover change, crop management, and other applications [17].
A worse performance of A'rou was found in this paper. The potential possibility might be the high elevation (over 3000 m), freeze/thaw problems and the cloud and rainy weather condition in the mountain area.
Another issue affecting the validation results is the scaling mismatch problems caused by the different representativeness between the ground-based and satellite-based LSTs. Many efforts have been made to address this issue in the last decades [19,25,36]. Most of the previous studies of LST validation have focused on the homogeneous and flat areas in order to avoid the thermal heterogeneity. In fact, the ideal validation area is not easy to find. Thus, the satellite-based LST product validation is often restricted by the limited sites, especially in China. However, in order to better use satellite-based LST products in different applications, the assessments of the satellite-based LSTs for all kinds of land surface types are necessary.

Temperature Trend in HRB
Some previous efforts have made on the analysis of HRB temperature warming trend from satellite based LST and ground based measurements to simulated models. Jin et al. [35] used 11 years (2000-2010) of MODIS monthly observation data to show a slight warming trend over the Tibetan Plateau, in which HRB is located. Luo et al. [37] used a Mann-Kendall trend test to analyze the dynamics of atmospheric temperature in HRB and concluded that annual mean temperature increased significantly in the entire HRB during 1980-2009. They also concluded that the downstream HRB contribute largely to the whole HRB warming trend. Li and Wang [38] found that on the annual scale, time series of averaged atmospheric temperature at all the selected stations from HRB presented remarkable increasing trends.
In this work, the nine years AATSR-derived LST at the A'rou site showed an increase from 2003 to 2011. The increasing trend at Yingke site is also apparent, although a little decrease occurs for 2011. For the whole HRB region, the increasing trend was not obvious for 2003-2011. The warming trend is clearly shown in the downstream of HRB. This confirmed the conclusion made by Luo et al. [37] as the downstream make the largest contribution to the warming trend in HRB. The warming trend results in this work provide another independent evidence for the temperature change studies over HRB region.
The long-term variation analysis at A'rou and Yingke sites also demonstrated the need for the further cloud pixel detection processes in the LST product. As mentioned by Kogler et al. [14], the cloud contamination might affect the long-term LST warming/cooling analysis. In addition, the trend analysis might also be impacted by multi-sensors calibration biases, as well as the degradation of the sensors themselves. Therefore, improved and strict cloud-detection algorithms and multi-sensor calibration from long-term LST data are quite necessary for the LST dataset utilizations.

Conclusions
The updated long-term AATSR-derived LST products were compared with ground based LST measurements from 2007 to 2011 at the two validation sites in HRB, China. The comparison was conducted with simultaneous in situ temperature measurements at areas with the homogeneous land use but with different elevations. Besides the cloud effect and other data accuracy problems in the LST retrievals, i.e., the time and viewing angle differences and observation footprint differences, the results from the two validation sites indicate that the AATSR-derived LST products yield a reasonable accuracy, especially at nighttime. For nighttime, the averaged bias, STD, RMSE and R 2 over both sites are 0.67 K, 3.03 K, 3.13 K and 0.93 K, respectively. In addition, the Yingke site with an elevation of 1519 m performs better than the A'rou site at 3032.8 m. The corresponding nighttime results from Yingke site are 0.15 K, 2.74 K, 2.74 K and 0.96 K.
Next, the time series of AATSR-derived LST over HRB and two individual sites were analyzed for 2003-2011. The long term series of the sites showed some cloudy possibility in A'rou and Yingke sites which is very useful for the further cloud pixel detection processes. Moreover, the nine years AATSR-derived LST from the A'rou site for daytime showed an obvious increase from 2003 to 2011. Due to the low value for daytime of 2011, the LST trend is not obvious. For the whole HRB region, the warming trend is clearly shown in the downstream of HRB.
The future inter comparison between different satellites derived LST should be done to validate the AATSR-derived LST and make an inter-calibration results avoiding the different scale representatives between satellite derived LST and in situ measurements. Further researches extending or updating the satellite observed LST to the past (e.g., the 1980s) will make a better understanding of the local temperature change trend, as well as the global/regional warming.