Ground surface temperature (GST) is an integrated result of heat and moisture exchanges through the ground surface.1 In permafrost regions, GST can serve as upper boundary conditions in permafrost models for modeling the freeze-thaw depths of the active layer and the thermal properties of permafrost, and is, therefore, a crucial parameter for the thermal state of the active layer and permafrost.1 Thus, accurate estimation of GST in permafrost regions spatially and temporally could provide valuable dataset in the research on not only permafrost, but also hydrology, ecology, and climate.23.–4 Unfortunately, GST is not a standard meteorological measurement and is rarely available for a specific case.
Accompanied with the development of remote sensing technology, satellite data make it possible to measure land surface temperature (LST) over the entire globe with sufficiently high temporal resolution and with complete spatial coverage rather than only point values.5 Currently, LST can be derived from sensors, including the Landsat Thematic Mapper/Enhanced Thematic Mapper Plus, the advanced spaceborne thermal emission and reflection radiometer, the advanced very high-resolution radiometer, the advanced along-track scanning radiometer, the spinning enhanced visible and infrared imager, and the moderate-resolution imaging spectroradiometer (MODIS),220.127.116.11.11.–12 and especially with improved moderate-resolution satellite sensors, such as the MODIS instruments onboard the Terra and Aqua platforms.
The MODIS Terra and Aqua satellites are in Sun-synchronous orbit near polar orbits, with Aqua in an ascending orbit and Terra in a descending orbit and equatorial crossings at 10:30 a.m. for Terra and at 1:30 p.m. for Aqua (local solar time). MODIS Terra data became available in February 2000, and Aqua data became available in July 2002. The generalized split-window algorithm is used to derive LST from bands 31, 32 and known emissivities.13 Comparisons between the MODIS LST products and in situ measurements indicate that the accuracy of the MODIS LST products is better than 1 K for a given observation time and angle. Therefore, the MODIS LST products have been widely used in various studies,14,15 including research on the cryosphere,16 which is sensitive to climate changes.
Integrated mean daily values from Aqua and Terra daytime and nighttime LST acquisitions in Quebec and the northern slope of Alaska were performed by Hachem et al.16,17 The results showed good agreement with air temperature and were used to study the spatial distribution of permafrost. However, some research found that the performance of MODIS LST over the regions of wet polygonal tundra in Siberia1 and high Arctic tundra in Svalbard18 was affected by cloud cover, water bodies, snow cover, and soil properties.
Despite the fact that MODIS LST has been widely used in the circum-Arctic regions, there are few reports to evaluate the LST in permafrost regions of high-altitude areas. Permafrost on the Tibetan Plateau (TP) occupies in area1920.–21 and exerts great influences on regional and global climates through thermal forcing mechanisms.22 Many studies have analyzed climate change in the TP using observational data from meteorological stations,23 and the permafrost is considered as one of the most sensitive indicators. For low-statured vegetation, GST is likely a better indicator than air temperature when used to evaluate the thermal condition of permafrost regions.24 Wu et al.25 analyzed the changes of the mean annual GST and the surface freezing and thawing indices over the period of 1980 to 2007 based on the GST records from 16 meteorological stations, and the results indicated that the dramatic ground surface warming could have significant influence on the change of the permafrost thermal regime in the study region. However, meteorological stations are sparse in the west of TP, and a spatial expansion of temperature combined with remote sensing data is essential to assess the regional ground thermal regime in cold regions.
The MODIS sensors observed the instantaneous LST during the satellite overpass times, while the data that we need to analyze or drive the permafrost model mostly are the mean daily GST. Previous studies focused on the method to derive spatial distribution of mean daily air temperature in the TP, and most of these studies aimed at using MODIS LST combined with the statistical approach or the temperature-vegetation index.26,27 There are only few studies on deriving the mean daily GST from MODIS LST, which is mainly due to the scarcity of meteorological stations and lack of consecutive and high temporal resolution instrumental records.
It is well known that single daytime MODIS observations, which are significantly influenced by solar radiation and cloud cover, are used for the model. The higher observation frequency (for example, combined daytime and nighttime LST) used in models would generate more accurate spatial GST products.28 Meanwhile, the accuracy is influenced by the number and period of the available pixels of MODIS LST.
This study estimates mean daily GST using empirical models combining single daily tiled daytime and nighttime MODIS LST observations (considering the acquisition time during the day). A key difference from the previous approaches is that the model training datasets consisted of the automatic weather station (AWS) observations at the satellite overpass times. Three models were employed in this study according to different compositions of daytime/nighttime MODIS LST as predictors. Data from 2012 at three AWS were used to establish the empirical models, and the models were assessed with cross-validation.
In order to establish the suitable methods to estimate the GST in permafrost regions of TP, the study area was arbitrarily selected within the permafrost region of TP. It is located in central TP with an area of . Three AWS were set up in the study area along the Qinghai-Tibet Highway (Fig. 1 and Table 1). These stations provide diurnal temperature curves, which are used in conjunction with MODIS LST to predict mean daily GST.
Description and location of automatic weather stations in the central Tibetan Plateau.
|Station||Elevation (m)||Lat./Long.||Vegetation type||Measurement frequency|
|Xidatan (XDT)||4538||35.72°N; 94.13°E||Alpine meadow||Half hour average of 5 min readings|
|Wudaoliang (WDL)||4783||35.22°N; 93.08°E||Alpine desert||Half hour average of 5 min readings|
|Tanggula (TGL)||5100||32.58°N; 91.86°E||Alpine meadow||Half hour average of 5 min readings|
Data and Methods
Automatic Weather Stations
Three AWS provide GST observation data as half-hour averages of measurements recorded at 5-min intervals. The Xidatan station (XDT) and Tanggula station (TGL) are located in alpine meadow, and the Wudaoliang station (WDL) is in alpine desert. The distances from XDT to WDL, WDL to TGL, and XDT to TGL are 110, 313, and 406 km, respectively. GST was measured with an SI-111 precision infrared radiometer (Campbell Scientific Inc., Logan, Utah, USA). The infrared radiometer was typically installed 2 m above ground, with an effective observation area of and an accuracy of . The dataset used in this study covered the year 2012. Table 1 provides information regarding the location, elevation, and other characteristics of the three AWS.
MODIS LST Data
The satellite data used in this study are the 1 km gridded clear-sky reprocessing version 5 MOD11A1 and MYD11A1 h26 v05 tile, which have been downloaded from the Level 1 and Atmosphere Archive and Distribution System.29 For each of the 366 days of the year 2012, daytime and nighttime LST, the corresponding quality control (QC), and the coincident local solar view time (VT) were retrieved from the downloaded LST products at pixels identified by the geographical locations of the three AWS. LST values with low-quality QC values due to cloud cover or other processing failures, were omitted. After quality filtering, to 56.1% of LST values remained for the three AWS. Table 2 shows the clear-sky days at four satellite overpass times and complete observation days (both the daytime and nighttime are clear-sky) of Terra MODIS and Aqua MODIS. The complete observation days are to 34.7% for only Terra MODIS or Aqua MODIS and to 20.5% for both of them in 2012. The filtered LST values were then converted to Celsius from the recorded scaled temperature values in Kelvin with a scale factor of 0.02. The ground observations of GST and the MODIS LST ranged from to 13.19°C and to 18.17°C, respectively. The minimum and maximum LST values are biased since only clear-sky pixels were considered.
Statistics of clear-sky days at three automatic weather stations.
|Clear-sky days||Complete observation daysa|
|T_day||T_night||A_day||A_night||Terra||Aqua||Terra and Aqua|
|XDT||177 (48.4%)||187 (51.1%)||153 (41.8%)||220 (60.1%)||93 (25.4%)||109 (29.8%)||52 (14.2%)|
|WDL||190 (51.9%)||217 (59.3%)||142 (38.8%)||244 (66.7%)||122 (33.3%)||102 (27.9%)||54 (14.8%)|
|TGL||212 (57.9%)||222 (60.7%)||151 (41.3%)||237 (64.8%)||127 (34.7%)||113 (30.9%)||75 (20.5%)|
Note: T_day, Terra_daytime; T_night, Terra_nighttime; A_day, Aqua daytime; A_night, Aqua nighttime.
Both the daytime and nighttime are clear skies.
Data Processing and Model Establishment
Due to the intrinsic scanning characteristics of the MODIS instrument onboard the polar-orbiting satellites, the differences in local solar VT for the same pixel on different days in one revisit period (16 days for reaching the same ground track) may reach up to 2 h. Figure 2 shows MODIS Terra satellite daytime overpass times at the TGL station within one revisit period (from January 1 to January 16, 2012), and the other three overpass times (Terra satellite nighttime, Aqua satellite daytime, and Aqua satellite nighttime overpass) have corresponding and similar fluctuations. To match the AWS observation time, the local solar VT was converted to local standard time (Beijing time).
The models establishing procedure was composed of three steps, outlined in detail below.
Step 1: Select the AWS half-hour GST observations closest to satellite overpass times (hereafter MODIS-GST) as predictor(s) for each of the three AWS. Then, we explore the linear correlations between the mean daily GST (hereafter meteo-GST) and different predictor(s) combinations of daytime/nighttime MODIS-GST with data for all 366 days of the year 2012. The meteo-GST were averages of the 48 half-hour observation values in one day. The candidate predictor combinations include , , , and ; and ; and ; and , , , and , as shown in Table 3.
Step 2: By comparing determination coefficients () and root mean square errors (RMSE) of the linear regressions, proper predictor(s) and statistic correlations were selected as empirical models. Then, the models established by observations from one of the three AWS were applied to MODIS LST data to estimate meteo-GST, and observations from the remaining two were used to validate the models.
Step 3: The models’ validation was performed by using statistical parameters, such as the mean error (ME), mean absolute error (MAE), and RMSE. The ME, MAE, and RMSE were calculated as follows:
Description of the terminology.
|GSTmoddt (°C)||GST observed by AWS at Terra daytime overpass time|
|GSTmodnt (°C)||GST observed by AWS at Terra nighttime overpass time|
|GSTmyddt (°C)||GST observed by AWS at Aqua daytime overpass time|
|GSTmydnt (°C)||GST observed by AWS at Aqua nighttime overpass time|
|GSTmean (°C)||Mean daily GST observed by AWS|
|LSTmoddt (°C)||Daytime LST obtained by Terra|
|LSTmodnt (°C)||Nighttime LST obtained by Terra|
|LSTmyddt (°C)||Daytime LST obtained by Aqua|
|LSTmydnt (°C)||Nighttime LST obtained by Aqua|
|LSTmean (°C)||Mean daily LST estimated by the models|
Note: GST, ground surface temperature; AWS, automatic weather station; LST, land surface temperature; moddt, Terra daytime; modnt, Terra nighttime; myddt, Aqua daytime; mydnt, Aqua nighttime.
Relationships Between MODIS-GST and Meteo-GST
The statistics of linear regressions between the MODIS-GST and the meteo-GST at three AWS are shown in Table 4. The determination coefficients () were for the relationship between only daytime MODIS-GST and meteo-GST. When only nighttime MODIS-GST were used as predictors of the regression equations, the were greater than 0.88 and more strongly correlated than daytime MODIS-GST. The RMSE were 4.38 to 5.43°C and 2.77 to 3.52°C for regressions with daytime and nighttime MODIS-GST, respectively. Both indicators in Table 4 show that the nighttime MODIS-GST was better than daytime MODIS-GST for deriving meteo-GST.
Statistics of linear regressions between meteo-GST and MODIS-GST with data from 366 days of the year 2012 at three AWS.
|GSTmoddt (°C), GSTmodnt (°C)||732||0.95||2.13||0.96||2.03||0.95||2.29|
|GSTmyddt (°C), GSTmydnt (°C)||732||0.97||1.6||0.98||1.46||0.97||1.9|
|GSTmoddt (°C), GSTmodnt (°C), GSTmyddt (°C), GSTmydnt (°C)||1464||0.99||1.06||0.99||0.95||0.99||1.17|
Note: XDT, statistic with Xidatan AWS observation; WDL, statistic with Wudaoliang AWS observation; TGL, statistic with Tanggula AWS observation.
When employing both daytime and nighttime MODIS-GST to estimate meteo-GST using regression equations, the statistical reached 0.95, and the RMSE decreased to (Table 4). Compared to estimates derived using daytime or nighttime MODIS-GST alone, the combined daytime and nighttime MODIS-GST data significantly improved the performance of the regression equations for estimating . In Table 4, the MODIS-GST of Aqua satellite, and combination, yielded of 0.97 to 0.98, and the RMSE was 1.46 to 1.9°C; corresponding indicators in Table 4 for and combination of Terra satellite were 0.95 to 0.96 and 2.03 to 2.29°C, respectively. The Aqua satellite provides a better estimation of meteo-GST compared to the Terra satellite in regression equations.
Further improvement was observed when employing all four MODIS-GST in the regression equation (Table 4), yielding of 0.99 and RMSE of 0.95 to 1.17°C. Compared to the performance of the combined daytime and nighttime MODIS-GST data, three models were created with and , and , and , , , and as independent predictor combinations for meteo-GST estimations (Table 5). The regression coefficients of the equations are listed in Table 5.
Empirical models and their coefficients.
|III||LSTmoddt, LSTmoddt, LSTmyddt, LSTmyddt||LSTmean=a1×LSTmoddt+a2×LSTmodnt+a3×LSTmyddt+a4×LSTmydnt+b||a1=0.144||a1=0.156||a1=0.18|
Note: moddt, Terra daytime; modnt, Terra nighttime; myddt, Aqua daytime; mydnt, Aqua nighttime.
Results of Model Validation
Estimations of mean daily LST () were made by the three models with the MODIS LST data (Table 5). Figure 3 demonstrates the models’ estimated against the ground observations of ; the slopes of the regression lines of the scatter plots were 0.97 to 1.12 at three AWS. Compared to models I and II, model III yielded a higher (0.93, 0.92, and 0.91 for XDT, WDL, and TGL, respectively) (Fig. 3).
Table 6 lists the statistics of the estimation residuals. The XDT observations in Table 6 indicate the empirical models established by the observations at XDT station and validated with observations at WDL and TGL stations (so to with the WDL and TGL observations). The annual averages of the ME of the three models were near zero and ranged from to 1, with slightly negative values with WDL observations and positive values with XDT and TGL observations. Table 6 showed that model III yielded the lowest MAE and RMSE, 2.28 to 2.42°C and 2.96 to 3.05°C, of all the three AWS observations, respectively. The WDL observations performed best in model III, with ME, MAE, and RMSE of , 2.28, and 2.97°C, respectively. For model III, there was no obvious difference in the MAE and RMSE of the three AWS observations, which were and 3°C, respectively. Although the GST estimated from model III at each of the three AWS provided the best agreement with meteo-GST, the number of MODIS LST available was reduced to to 20.5%, which indicates that model III was not always feasible for the study area. This is a reminder that model I and model II should be supplemental methods. Compared to model II, model I had a lower MAE (2.48 to 2.71°C) and RMSE (3.24 to 3.53°C), which indicates that the Terra MODIS LST performed better than the Aqua MODIS LST when estimating meteo-GST in the study area. Generally, those models performed well in spite of some possible bias caused by the clouds.
Statistics of the estimation residuals of the three models.
|Model no.||XDT observations||WDL observations||TGL observations|
Note: XDT observations, the empirical model established by the observation data at XDT station and validated with WDL and TGL stations; WDL observations, the empirical model established by the observation data at WDL station and validated with XDT and TGL stations; TGL observations, the empirical model established by the observation data at TGL station and validated with XDT and WDL stations; ME, mean error; MAE, mean absolute error; RMSE, root mean squared error.
The mean daily GST (meteo-GST) estimated by the empirical models combining MODIS LST showed a good agreement with AWS observations. There are two major reasons leading to this. First is the similar observation principle between infrared radiometer and sensors on board satellites though MODIS sensors considered the atmospheric transmission. Second, the training datasets we used in following statistical analysis and models establishment with mean daily GST are the AWS observations at the satellite overpass times rather than the MODIS LST.
To compare the accuracy with routine modeling methods that have been used in previous study,17 MODIS LST data have also been applied as training datasets to establish the empirical models. Table 7 shows the resulting estimation residuals. Almost all the ME, MAE, and RMSE of the three models in Table 7 were larger than those in Table 6, which indicates that the models established by the AWS observations at the satellite overpass times were better than the models established by MODIS LST. The MAE and RMSE for model III improved to 0.85°C and 0.36 to 0.83°C, respectively. Though the models established by the AWS observations at the satellite overpass times cannot solve the mismatch of point and 1 km per pixel scales, the increased number of training data theoretically still improves the accuracy of estimated meteo-GST.
Statistics of the estimation residuals of the three models by training data of MODIS LST.
|Model no.||XDT observations||WDL observations||TGL observations|
Note: XDT observations, the empirical model established by the observation data at XDT station and validated with WDL and TGL stations; WDL observations, the empirical model established by the observation data at WDL station and validated with XDT and TGL stations; TGL observations, the empirical model established by the observation data at TGL station and validated with XDT and WDL stations.
During the daytime, the difference between MODIS LST and AWS observations is mainly controlled by the surface energy balance, which is a complex system dependent on information not easily available (such as solar radiation, cloud cover, wind speed, soil moisture, and surface roughness).30 During the night, solar radiation does not affect the thermal infrared signal. Therefore, nighttime LST is a better predictor than daytime LST for estimating meteo-GST. This also can be proved from the regression coefficient in Table 5. Models coupling nighttime and daytime LST were more accurate than those using nighttime or daytime LST alone.28 Correlation analysis between the mean daily GST and the MODIS-GST showed that the Terra MODIS performed better than the Aqua MODIS, while Aqua MODIS outperformed Terra MODIS when applied the models to estimated meteo-GST. The model results using nighttime and daytime LST reflect real LST because they reduce the possible bias caused by the clouds. Further analysis demonstrates that the coupled nighttime and daytime LST of Terra and Aqua will perform better than using Terra or Aqua alone. For regions with extensive cloud cover that cannot obtain the daytime and nighttime LST at the same day, Williamson et al.24 proposed the interpolated curve mean daily surface temperature method, which requires only a single daytime observation.
There is no significant difference in the MAE and RMSE of the three AWS observations, and 3°C, respectively, which indicates that the deviation error is no different when a specific the AWS observation is selected to establish the models. However, whether the ME would be negative or positive depends on the difference of meteo-GST between the modeling station and validation stations being lower or higher. The morphology of areas where the three AWS located are flat and the underlying surfaces are alpine meadow and alpine desert, which can represent the general characteristic of the central TP. Using an averaged result or dividing the area into several subregions (such as utilizing the Thiessen polygon)31 when used for regional GST estimates will be a better method.
The difference between estimated and observed mean daily GST was likely due to the combination of several factors, including temperature differences due to time offset between AWS half-hour observation and satellite overpass time and cloud contamination and cloud shadow causing differences in LST across spatial and temporal scales. Although the MODIS cloud-detecting algorithm takes as many as 22 of the 36 MODIS bands to maximize the reliability of cloud detection, uncertainties still exist in the MODIS cloud mask product applications.7 For some instances, the residual clouds contaminate the data, and in others, the overdetection of clouds diminishes the availability of clear observations, which leads to many fill-values in the downstream products (such as LST); these effects could affect data availability. For example, the false-positive cloud flags in hot and bright pixels likely lead to a summertime cooling bias in LST over urban areas and the presence of residual clouds would cause misidentification of vegetation changes.3233.–34 Future work should be devoted to the identification of cloud-contaminated LST data and the degree to which contaminated LST values are affected. By quantifying the degree of contamination, data volume will be preserved and quality improved. Furthermore, mismatches in spatial and temporal scale should also be investigated. Statistical techniques should be improved to make the most use of limited ground station data. The implementation of empirical models in this study could provide a powerful technique to obtain the spatial distribution of GST using the MODIS LST products.
To assess the potential of remotely sensed data to estimate GSTs, an empirical model was established based on three AWS in the central TP. The empirical models were established by the correlation between the mean daily ground surface temperature (meteo-GST) of the three AWS observations and the AWS observations at the satellite overpass times, with three different combinations of daytime/nighttime MODIS LST (Terra daytime and nighttime, Aqua daytime and nighttime, and Terra and Aqua daytime and nighttime). Employing both daytime and nighttime from Terra or Aqua MODIS as predictors performed better than using daytime or nighttime alone. Among the three models, the model established by observations of Terra and Aqua daytime and nighttime, MODIS LST, yielded the highest and the lowest MAE and RMSE, but the number of available pixels was substantially reduced. There was no obvious difference in , MAE, and RMSE of three AWS observations, , 2.4°C, and 3°C, respectively. Terra MODIS LST performed better than Aqua MODIS LST to estimate meteo-GST. Models established by the AWS observation values at the MODIS overpass times were better than the models established by the MODIS LST observations. These three AWS are located in the central TP and, therefore, indicate the thermal conditions of permafrost, so the estimation of mean daily GST could be a powerful dataset for the monitoring and modeling of permafrost.
The authors thank the Level 1 and Atmosphere Archive and Distribution System for providing MODIS land products. This research received financial support from the National Major Scientific Project of China “Cryospheric Change and Impacts Research” (2013CBA01803), the National Natural Science Foundation of China (Grant Numbers: 41271086, 41271081, and 41101069), and the Hundred Talents Program of the Chinese Academy of Sciences granted to Tonghua Wu (51Y251571).
Defu Zou is pursuing his PhD in the Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou, China. He received his bachelor’s degree in 2009 and his master’s degree in 2012 from Lanzhou University, Lanzhou, China. His main research interests focus on estimation of land surface variables from satellite observations, studies on Tibetan Plateau permafrost distribution mapping, and the active layer thickness modeling.
Lin Zhao is the group leader of the Cryosphere Research Station on Qinghai-Tibetan Plateau, State Key Laboratory of Cryospheric Science, Cold and Arid Regions Environmental and Engineer Research Institute, Chinese Academy of Sciences, where he got his PhD in 2004. He also is the national correspondent of the international permafrost monitoring network for China. His research interests include general geocryology, heat and moisture exchange processes within active layers, boundary layer processes, and ecology in permafrost regions.
Tonghua Wu studied permafrost and climate change at the Cryosphere Research Station on Qinghai-Tibetan Plateau, State Key Laboratory of Cryospheric Science, Cold and Arid Regions Environmental and Engineer Research Institute, Chinese Academy of Sciences, where he received his PhD degree in 2006. His interests as a research scientist included general geocryology, application of geographic information system and remote sensing techniques in permafrost-related studies, and climate changes in cold regions. At present, he performs his research on the Qinghai-Tibetan and Mongolian Plateau.
Xiaodong Wu is a research scientist in the Cold and Arid Regions Environmental and Engineering Research Institute. He received his PhD degree from the Nanjing Institute of Limnology and Geography, Chinese Academy of Sciences, in 2008. He is currently an associate professor in the State Key Laboratory of Cryospheric Sciences, Chinese Academy of Sciences. His research interests cover global change, ecology of permafrost regions, and soil carbon.
Qiangqiang Pang is a research assistant in Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou, China. He received his PhD degree from Cold and Arid Regions Environmental and Engineering Research Institute, Chinese Academy of Sciences, Lanzhou, China, in 2009. His current research interests include permafrost and environment in cold regions, and active layer hydrothermal processes.
Zhiwei Wang is a PhD student at the Cold and Arid Regions Environmental and Engineer Research Institute, Chinese Academy of Sciences. He received his BS degree from Inner Mongolia University for the Nationalities in 2007 and his MS degree from Lanzhou University in 2010. For more than four years, he has been studying the vegetation dynamics by NDVI and decision tree in permafrost zone of Qinghai-Tibet Plateau using MODIS, GIMMS, and ASTER GDEM v002.