A Calculation Model for Ground Surface Temperature in High-Altitude Regions of the Qinghai-Tibet Plateau, China

: As a major parameter in the energy balance of the ground surface, temperature represents the level of exchange of energy and moisture between the ground and air. The Qinghai-Tibet Plateau (QTP) has the permafrost region with the highest altitude and the largest area in low–middle latitude of the world. The variation in ground surface temperature has an impact on the existence and development of the permafrost. Therefore, the analysis of the ground surface temperature in the QTP is signiﬁcant to reﬂect the energy exchange in permafrost regions. This paper collected solar radiation data and calculated the conversion coefﬁcient from total solar radiation to long-wave radiation of the ground surface on different underlying surfaces. The ground surface temperature was inversely calculated and modiﬁed based on the reception of solar radiation on different underlying surfaces. A simpliﬁed calculation model of ground surface temperature was built to reﬂect the ground surface temperature on different underlying surfaces of the QTP. The calculation results were compared with MODIS and showed good ﬁtness, providing a systematic and reliable method for calculating the ground surface temperature on the QTP. The above model plays a signiﬁcant role in the estimation of soil moisture, ground surface energy and water balance.


Introduction
As a typical interior plateau in Asia, the Qinghai-Tibet Plateau (QTP) covers 2.5 million km 2 and the average altitude is higher than 4000 m; it is the largest plateau in China and the highest plateau in the world. The solar radiation controls the balance of ground surface radiation, energy exchange and the climate system. The influence factors of total solar radiation include water vapor content, visibility, aerosol, cloud, precipitation, latitude and altitude. The QTP has the largest area of permafrost in low-middle latitude of the world, in which the variation in the permafrost will affect the formation, change and development of the climate in Eastern China and East Asia. As the key parameter of the energy balance in the ground surface, the ground surface temperature reflects the exchange of energy and moisture between the ground surface and air [1]. Under the scenario of climate warming, the mean ground surface temperature has increased by 0.85 • C, which affects the existence and development of permafrost [2]. The thawing of the permafrost will release the methane and CO 2 stored in the permafrost, conversely facilitating the increase in climate warming [3][4][5]. Otherwise, the ground surface temperature is a key parameter of the ground heating field. The diurnal variation of the heating field in the QTP will cause the diurnal variation of the air flow field and affects the weather process in adjacent regions via the diurnal variation of local circulation in the QTP and surrounding regions [6].
The ground surface temperature was explored by several researchers, with promising achievements. The annual variation in ground surface temperature can be calculated by a theoretical model including solar radiation and the sun angle. Dumais et al. proposed a simplified model for calculating the road surface temperature. The model was inversely related to the surface temperature by inputting the data of road albedo, air temperature, wind speed and solar radiation [7]. Wang et al. analyzed the influence of solar radiation on the pavement type and underlying permafrost. The pavement surface temperature was inversely calculated by total radiation and atmospheric temperature, which provided a parameterized method for the calculation of the pavement surface temperature [8]. Hu et al. provided another parameterized method for the calculation of the surface temperature of an embankment and its slope in the Qinghai-Tibet Railway based on the absorbed solar radiation by the ground surface [9].
Moreover, the ground surface temperature can be inversely calculated by remote sensing technology. Weng et al. provided an integrated calculation method of ground surface temperature in permafrost regions of the QTP based on multi-source remote sensing [10]. Hu et al. built a down-scaling method to calculate the ground surface temperature based on the diurnal temperature cycle of the ground surface, which can obtain the ground surface temperature under the conditions of a sunny and cloudy sky [11]. Wu et al. proposed a model to calculate the glacier surface temperature combined with the Jiménez-Muoz (JM) model and Landsat 8 data [12]. Shen et al. explored the influence of temporal and spatial scale on the ground surface temperature (MODISLST) and air temperature, to seek an optimal temporal-spatial scale to improve the accuracy of calculation of the ground surface temperature using a MODIS dataset [13]. Zhou et al. proposed an improved temperature and emissivity separation (TES) algorithm, and the improved TES algorithm can be used to accurately derive LST&E from AHI data [14]. The existing temperature calculation model mainly focuses on the Qinghai-Tibet Engineering Corridor, which is based on monitoring data. The calculation results represent the given condition and not all the underlying surfaces.
As a key boundary condition for numerical simulation, the ground surface temperature can be fitted by field-measured data. The temperature variation can be precisely obtained by the above fitting method; however, the temperature distribution in the study region cannot be reflected. Hu et al. have fitted the regression equation of the temperature boundary based on measured temperature data along the Qinghai-Tibet Highway at the mileage of K2968 + 200 in the Qinghai-Tibet Plateau [15]. The temperature variation of the ground surface on both sides of an embankment was analyzed along the QTH and QTR. The correlation between the surface temperature of the embankment slope and the theoretical radiation value was revealed based on the field measurement results. Chou et al. built a regression equation to connect the ground temperature at 0.5 m depth and the radiation difference [16]. Huang et al. built a model to reflect the relationship between air temperature, wind speed, soil moisture and ground surface temperature. It is known from the above model that the ground surface temperature has an obvious response to air temperature throughout the whole year, has an obvious response to wind speed in winter and summer and has an obvious response to soil moisture in summer [17].
The incidence distribution of solar radiation is significant to the hydrology, ecology and carbon balance of a region. The atmospheric transmissivity has an influence on the energy budget on the ground surface by solar radiation incidence. The strong solar radiation on the Qinghai-Tibet Plateau, with high altitude, has a more obvious influence on the energy budget than in other regions. The present research on the inversion of the ground surface temperature mainly focuses on the Qinghai-Tibet Engineering Corridor, which can only reflect the specific engineering condition due to the lack of field-monitored data within a large range [18]. The boundary layer theory can calculate the underlying surface temperature based on air temperature. However, the temperature increase value cannot systematically represent the condition of each underlying surface due to the complexity of the field and radiation differences on the embankment slope surface [19]. This paper extracted the ground surface temperature from MODIS data on the different underlying surfaces of the Qinghai-Tibet Plateau (QTP) and analyzed the temperature characteristics on six types of underlying surface. The coefficient of conversion rate and correction was obtained in the above underlying surface, which provides a reliable calculation method for the temperature inversion in different underlying surfaces on the QTP.

Study Region
The study region in this paper is the QTP (26 • 00 12 N~39 • 46 50 N, 73 • 18 52 E~104 • 46 59 E), which is located in the south-west of China. The altitude of the QTP is above 4000 m, and it is referred to as "the roof of the world". The unique local environment has created different underlying surfaces, including bare land, wetland, grassland, forests, shrub land, water bodies, building land and glaciers ( Figure 1). cannot systematically represent the condition of each underlying surface due to the plexity of the field and radiation differences on the embankment slope surface [19 paper extracted the ground surface temperature from MODIS data on the differe derlying surfaces of the Qinghai-Tibet Plateau (QTP) and analyzed the temperature acteristics on six types of underlying surface. The coefficient of conversion rate an rection was obtained in the above underlying surface, which provides a reliable ca tion method for the temperature inversion in different underlying surfaces on the Q

Study Region
The study region in this paper is the QTP (26°00′12″N~39°46 73°18′52″E~104°46′59″E), which is located in the south-west of China. The altitude QTP is above 4000 m, and it is referred to as "the roof of the world". The uniqu environment has created different underlying surfaces, including bare land, we grassland, forests, shrub land, water bodies, building land and glaciers ( Figure 1).

Data
The digital elevation model (DEM) was provided by the Geospatial Data Clou Computer Network Information Center, Chinese Academy of Sc (http://www.gscloud.cn (accessed on 1 December 2021.)). The land cover dataset wa the GlobeLand30 (http://www.globallandcover.com (accessed on 1 December which comprises datasets collected at 30-m resolution, more than ten times that of ous datasets. The above datasets are valuable for monitoring environmental chang for resource management at global, regional and local scales. The overall accuracy GlobeLand30 V2020 and the Kappa coefficient are 85.72% and 0.82, respectively. Th cover dataset is shown in Table 1.

Data
The digital elevation model (DEM) was provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http://www. gscloud.cn (accessed on 1 December 2021)). The land cover dataset was from the Glo-beLand30 (http://www.globallandcover.com (accessed on 1 December 2021)), which comprises datasets collected at 30-m resolution, more than ten times that of previous datasets. The above datasets are valuable for monitoring environmental changes and for resource management at global, regional and local scales. The overall accuracy of the GlobeLand30 V2020 and the Kappa coefficient are 85.72% and 0.82, respectively. The land cover dataset is shown in Table 1. The MODIS MOD11A2 dataset of ground surface temperature was provided by the Geospatial Data Cloud site, Computer Network Information Center, Chinese Academy of Sciences (http://www.gscloud.cn (accessed on 1 December 2021)). The dataset was merged, projected, clipped and converted to obtain the ground surface temperature in the QTP with a 1 km resolution.

Ground Surface Temperature Calculation Model
Based on the Stefan-Boltzmann Law, the radiant temperature on different underlying surfaces was calculated by the long-wave radiation of the ground surface. The conversion scale from the total solar radiation to long-wave radiation of the ground surface is different on different underlying surfaces, among which the heat exchange is different and is represented as the value of ground surface temperature.

Direct Solar Radiation on Different Underlying Surfaces
The following equations governing the net energy budget of the Earth system were used in this study: where G is the ground heat flux, H is the sensible heat flux, E is the evapotranspiration rate and λ is the latent heat of vaporization.
where Tr is the air temperature within the lowest layer of the atmospheric model, γa is the aerodynamic resistance between the surface and the lowest layer of the atmosphere, ρ and cp are the density and specific heat of air, β is the moisture availability function (0 ≤ β ≤ 1), e* (Ts) is the saturated vapor pressure at temperature T s , e r is the vapor pressure within the lowest layer of the atmospheric model and g is the psychrometric constant. Solar radiation is the most direct factor that results in the variation in ground temperature and air temperature. The absorbent solar radiation on a level ground surface can be expressed as follows: in which R s is the absorbent solar radiation on a level ground surface; S 0 is the solar constant, S 0 = 1353.3 W·m −2 , without consideration of atmospheric influences; Z is the solar zenith angle ( Figure 2). Solar radiation on the Qinghai-Tibet Plateau is negatively correlated with latitude and positively correlated with altitude and relative sunshine duration. The sunshine duration, latitude and altitude have a dominant influence on the solar radiation of the ground surface [20]. Therefore, the influence of atmosphere was ignored in this study.
zenith angle ( Figure 2). Solar radiation on the Qinghai-Tibet Plateau is neg lated with latitude and positively correlated with altitude and relative sunsh The sunshine duration, latitude and altitude have a dominant influence on t ation of the ground surface [20]. Therefore, the influence of atmosphere w this study. Z is calculated as follows: cos sin sin cos cos cos Z ϕ δ φ δ ψ = + in which φ is the latitude, is the solar declination and ψ is the time angle. The solar zenith angle is closely related to the regional longitude, latitu The variation in the zenith angle over time in 15°N, 30°N, 45°N and 60°N Figure 3.  Z is calculated as follows: cos Z = sin ϕ sin δ + cos φ cos δ cos ψ (5) in which ϕ is the latitude, δ is the solar declination and ψ is the time angle. The solar zenith angle is closely related to the regional longitude, latitude and time. The variation in the zenith angle over time in 15 • N, 30 • N, 45 • N and 60 • N is shown in Figure 3.
in which Rs is the absorbent solar radiation on a level ground surface; S0 is the solar stant, S0 = 1353.3 W·m −2 , without consideration of atmospheric influences; Z is the s zenith angle (Figure 2). Solar radiation on the Qinghai-Tibet Plateau is negatively co lated with latitude and positively correlated with altitude and relative sunshine dura The sunshine duration, latitude and altitude have a dominant influence on the solar r ation of the ground surface [20]. Therefore, the influence of atmosphere was ignore this study. Z is calculated as follows: in which φ is the latitude, is the solar declination and ψ is the time angle. The solar zenith angle is closely related to the regional longitude, latitude and t The variation in the zenith angle over time in 15°N, 30°N, 45°N and 60°N is show Figure 3.  The zenith in winter is larger than that in summer, and the solar altitude in summer is larger than that in winter. The solar zenith changes with time in a one-day cycle, which has the minimum value, and the solar altitude has the maximum value at noon. The heating area on the ground surface is small and the heat is concentrated with a large solar altitude, which has more solar radiation. With the increase in latitude, the difference in the solar zenith angle in spring and winter is increasingly obvious. The zenith angle at the summer solstice is smaller than that at the winter solstice. In 24 h during the winter solstice, the order of the solar zenith in different latitudes from high to low is 30 The declination δ (Figure 4) was calculated as follows: The time angle ψ (Figure 4) was calculated as follows: in which N is the data order: common year, 1 January is and 31 December is 365; leap year, 1 January is and 31 December is 366; ψ is 0 based on the local true solar time at noon.
The zenith in winter is larger than that in summer, and the solar altitude in summer is larger than that in winter. The solar zenith changes with time in a one-day cycle, which has the minimum value, and the solar altitude has the maximum value at noon. The heating area on the ground surface is small and the heat is concentrated with a large solar altitude, which has more solar radiation. With the increase in latitude, the difference in the solar zenith angle in spring and winter is increasingly obvious. The zenith angle at the summer solstice is smaller than that at the winter solstice. In 24 h during the winter solstice, the order of the solar zenith in different latitudes from high to low is 30°N > 15°N > 45°N > 60°N.
The declination (Figure4) was calculated as follows: The time angle ψ (Figure 4) was calculated as follows: ( ) 15 12 t ψ = × − (7) in which N is the data order: common year, 1 January is and 31 December is 365; leap year 1 January is and 31 December is 366; ψ is 0 based on the local true solar time at noon.

Long-Wave Radiation on Different Underlying Surfaces
The long-wave radiation on ground surface Rl can be expressed as follows: l s l n R vR R = + (8) in which v is the conversion ratio from total solar radiation to long-wave radiation on the ground surface. The conversion ratio on different underlying surfaces in the QTP is shown in Table 2. Rl|n is the equilibrium value of long-wave radiation at night; Rs is the direct solar radiation on the ground surface.

Long-Wave Radiation on Different Underlying Surfaces
The long-wave radiation on ground surface R l can be expressed as follows: in which v is the conversion ratio from total solar radiation to long-wave radiation on the ground surface. The conversion ratio on different underlying surfaces in the QTP is shown in Table 2. R l | n is the equilibrium value of long-wave radiation at night; Rs is the direct solar radiation on the ground surface. The equilibrium value of long-wave radiation on the ground surface at night is the radiation value that reaches an equilibrium status whereby the ground surface receives solar radiation during the day and releases effective radiation at night with the thermal insulation of the atmosphere [15]. The sine function can precisely calculate the equilibrium value of long-wave radiation on the ground surface at night every day, and the fitted equation of measured data is shown as follows [26]: in which, A, B and C are coefficients fitted by the measured results, A is relative with altitude and latitude; C is π/3. A can be expressed as follows: in which H is the altitude, ϕ is the latitude and H a0 is the average thickness of the atmosphere, which is 8072.5 m in the QTP [9]. B can be calculated via A as follows [26]:

Thermal Boundary Model
According to the Stefan-Boltzmann Law, the radiation temperature on different underlying surfaces can be calculated based on the long-wave radiation of the ground surface as follows: in which T g is the radiation temperature, ε is the emissivity of a gray body, which is 1; σ is the Stefan-Boltzmann constant, which is 5.67 × 10 −8 W·m −2 ·K −4 . However, the above equation did not consider the influence of atmosphere and precipitation under an ideal condition. The calculation result has some errors because the radiation temperature is different from the real temperature on the ground surface. The calculation result has a linear relationship with the real ground surface temperature. Therefore, a linear equation was adopted to modify the ground surface radiation temperature as follows: in which k and b are the fitted coefficient; T g is the ground surface temperature. The k and b were fitted by the MODIS data of the ground surface temperature. The values of k and b of different underlying surfaces on the QTP are shown in Table 3.

Process of Surface Temperature Calculation
Based on the determination of latitude, altitude and solar radiation, the hour angle, declination angle and zenith angle were calculated in the QTP. Then, the ground surface temperature could be approximately solved by the modified and simplified model based on the conversion coefficient and correction coefficient of different underlying surfaces. The surface temperature calculation process utilized in this study is shown in Figure 5.
Therefore, a linear equation was adopted to modify the ground surface radiation temperature as follows: ' g g T kT b = + (13) in which k and b are the fitted coefficient; Tg is the ground surface temperature. The k and b were fitted by the MODIS data of the ground surface temperature. The values of k and b of different underlying surfaces on the QTP are shown in Table 3.

Process of Surface Temperature Calculation
Based on the determination of latitude, altitude and solar radiation, the hour angle, declination angle and zenith angle were calculated in the QTP. Then, the ground surface temperature could be approximately solved by the modified and simplified model based on the conversion coefficient and correction coefficient of different underlying surfaces. The surface temperature calculation process utilized in this study is shown in Figure 5.

The Mean Monthly Ground Surface Temperature (MMGST)
The variation in mean monthly ground surface temperature (MMGST) calculated from the radiant ground surface temperature is shown in Figure 6. The variation in MMGST has a similar trend on different underlying surfaces; it reaches the maximum during May to August and reaches the minimum during December to January. The

The Mean Monthly Ground Surface Temperature (MMGST)
The variation in mean monthly ground surface temperature (MMGST) calculated from the radiant ground surface temperature is shown in Figure 6. The variation in MMGST has a similar trend on different underlying surfaces; it reaches the maximum during May to August and reaches the minimum during December to January. The MMGST from MODIS can precisely validate the results of the modified calculation model on the ground surface temperature of different underlying surfaces.  The variation in MMGST for bare land is similar to that of grassland. The calculated MMGST is smaller than the MODIS data value. The MMGST of bare land in January, November and December is lower than 0 °C, and it is higher than 0 °C from February to October. The maximum value of MMGST reaches 20.74 °C in July, and the minimum reaches −4.29 °C in December. MMGST before the modification is higher than the MODIS result; the maximum difference is 7.79 °C in July, of which the minimum difference value is 2.49 °C after the modification.
The calculation results from the simplified model are larger than those of the MODIS data before correction. The minimum difference between the MODIS data and model results is 0.87 °C in July before correction, and it is 0.50 °C in November. The maximum temperature difference for water bodies before correction is 5.68 °C, and it is 1.0 °C after correction, with an increase of 82.39% in precision. The precision of the corrected temperature of wetland and forest increased by 95.35% and 87.07%, respectively. The study region has a high altitude; therefore, the influence of the latitude and underlying surface on MMGST is more significant than that of altitude. The variation in MMGST for forests, water bodies and glaciers throughout the whole year is smaller than that of wetland, grassland and bare land. The order of MMGST is forest > wetland > water body > glacier.

The Mean Annual Radiant Ground Surface Temperature (MARGST)
The comparison of the calculation results of the mean annual radiant ground surface temperature (MARGST) and modified MARGST on different underlying surfaces is shown in Figures 7-9. The correlation coefficient R 2 , root-mean-square error (RMSE) and deviation bias statistics were obtained during the fitting. The calculation results of MARGST before modification were under the line y = x, which means that the calculation results are small, with the reason being that the bare land is relatively dry (Figure 6). The temperature after modification scatters uniformly on the two sides of the line y = x. The bias before and after the modification is 6.36 and 6.25, respectively, increasing by 1.73%. The variation in MMGST for bare land is similar to that of grassland. The calculated MMGST is smaller than the MODIS data value. The MMGST of bare land in January, November and December is lower than 0 • C, and it is higher than 0 • C from February to October. The maximum value of MMGST reaches 20.74 • C in July, and the minimum reaches −4.29 • C in December. MMGST before the modification is higher than the MODIS result; the maximum difference is 7.79 • C in July, of which the minimum difference value is 2.49 • C after the modification.
The calculation results from the simplified model are larger than those of the MODIS data before correction. The minimum difference between the MODIS data and model results is 0.87 • C in July before correction, and it is 0.50 • C in November. The maximum temperature difference for water bodies before correction is 5.68 • C, and it is 1.0 • C after correction, with an increase of 82.39% in precision. The precision of the corrected temperature of wetland and forest increased by 95.35% and 87.07%, respectively. The study region has a high altitude; therefore, the influence of the latitude and underlying surface on MMGST is more significant than that of altitude. The variation in MMGST for forests, water bodies and glaciers throughout the whole year is smaller than that of wetland, grassland and bare land. The order of MMGST is forest > wetland > water body > glacier.

The Mean Annual Radiant Ground Surface Temperature (MARGST)
The comparison of the calculation results of the mean annual radiant ground surface temperature (MARGST) and modified MARGST on different underlying surfaces is shown in Figures 7-9. The correlation coefficient R 2 , root-mean-square error (RMSE) and deviation bias statistics were obtained during the fitting. The calculation results of MARGST before modification were under the line y = x, which means that the calculation results are small, with the reason being that the bare land is relatively dry (Figure 6). The temperature after modification scatters uniformly on the two sides of the line y = x. The bias before and after the modification is 6.36 and 6.25, respectively, increasing by 1.73%. The RMSE before and after the modification is 4.66 and 1.27, respectively, decreasing by 72.75%. It is shown that the calculation results of MARGST after modification can precisely reflect the real temperature condition. The RMSE before and after the modification is 4.66 and 1.27, respectively, decreasing by 72.75%. It is shown that the calculation results of MARGST after modification can precisely reflect the real temperature condition. The calculation results of MARGST on wetland are large before modification. The RMSE before and after modification is 26.99 and 0.91, respectively, decreasing by 96.63% (Figure 7). The bias before and after modification is 14.65 and 3.82, decreasing by 73.92%. The radiant ground surface temperature for forests is higher than that on wetland before correction. The RMSE and bias for forests decreased by 97.52% and 85.54% after correction, respectively. The RMSE and bias of water bodies before correction were 5.82, which was induced by high humidity on the water surface. After correction, the RMSE and bias decreased by 17.01% and 77.32%, respectively. The simplified calculation model increased the accuracy of the temperature calculation results. The calculation results of MARGST on wetland are large before modification. The RMSE before and after modification is 26.99 and 0.91, respectively, decreasing by 96.63% (Figure 7). The bias before and after modification is 14.65 and 3.82, decreasing by 73.92%. The radiant ground surface temperature for forests is higher than that on wetland before correction. The RMSE and bias for forests decreased by 97.52% and 85.54% after correction, respectively. The RMSE and bias of water bodies before correction were 5.82, which was induced by high humidity on the water surface. After correction, the RMSE and bias decreased by 17.01% and 77.32%, respectively. The simplified calculation model increased the accuracy of the temperature calculation results.

Results Validation
The altitude, longitude and latitude of different underlying surfaces on the Qinghai-Tibet plateau in 2015-2019 were obtained and inputted into the modified, simplified calculation model of MARGST. The calculation results were compared with the MODIS data on different underlying surfaces, which showed high goodness of fit ( Figure 10). The R 2 on different underlying surfaces is larger than 0.86, in the order of bare land > grassland >

Results Validation
The altitude, longitude and latitude of different underlying surfaces on the Qinghai-Tibet plateau in 2015-2019 were obtained and inputted into the modified, simplified calculation model of MARGST. The calculation results were compared with the MODIS data on different underlying surfaces, which showed high goodness of fit ( Figure 10). The R 2 on different underlying surfaces is larger than 0.86, in the order of bare land > grassland > forest > wetland > glacier > water body. The RMSE on different underlying surfaces ranges from 0.97 to 1.42, and bias is lower than 5.89, which represents the good fitting results and high precision of model calculation. forest > wetland > glacier > water body. The RMSE on different underlying surfaces ranges from 0.97 to 1.42, and bias is lower than 5.89, which represents the good fitting results and high precision of model calculation.

Discussion
In the Qinghai-Tibet Plateau, the absorbing capacity of solar radiation is different on different underlying surfaces, and the conversion rate from total solar radiation to longwave radiation on different underlying surfaces is different as well. The different underlying surfaces on the Qinghai-Tibet Plateau differ in optical, hydro-thermal, chemical and

Discussion
In the Qinghai-Tibet Plateau, the absorbing capacity of solar radiation is different on different underlying surfaces, and the conversion rate from total solar radiation to long-wave radiation on different underlying surfaces is different as well. The different underlying surfaces on the Qinghai-Tibet Plateau differ in optical, hydro-thermal, chemical and biochemical properties, as well in the process of ground surface radiation and energy exchange, which results in heterogeneity of land surface processes [27]. On different underlying surfaces, the total solar radiation is almost the same; however, the long-wave radiation is different, which is similar to the temperature distribution. The long-wave radiation was dominated by the canopy temperature, soil temperature and ground surface color [28].
The equilibrium value of long-wave radiation at night indicates a mechanism that absorbs solar radiation during the daytime and releases the effective radiant energy at night, which reaches an equilibrium state under the thermal insulation of the atmosphere. The equilibrium value of long-wave radiation at night varies with seasonal change and is related to the regional condition and altitude. In this paper, the order of the conversion rate coefficient on different underlying surfaces is grassland > wetland > forest > bare land > water body > glacier. Moreover, the ground surface temperature in high-altitude regions is influenced by ecological and environmental factors. The degradation of permafrost, decrease in seasonal snow and melting of glaciers will absorb thermal energy and lower the ground surface temperature. However, albedo decreased and the absorption capacity of solar radiation increased after the melting of snow and glaciers, which increased the ground surface temperature (Figure 11). The variation in MMGST on different underlying surfaces showed a normal distribution in the Qinghai-Tibet Plateau. The radiant energy had the maximum in the southwest of China throughout the year. The solar radiation reaches the maximum value in summer, and the ground surface temperature has the maximum from June to August [29]. The conversion ratio coefficient is different on different underlying surfaces; therefore, the conversion ratio from total solar radiation to long-wave radiation on the ground surface and the order of ground surface temperature is grassland > wetland > forest > bare land > water body > glacier. The MMGST obtained by the simplified calculation model has a good fit with that of MODIS data for the entire year.
The calculation results of MARGST from the simplified model have a correlation coefficient higher than 0.9 with remote sensing data. After the modification, the calculation results reduce the deviation and mean square root. The ground surface temperature on The variation in MMGST on different underlying surfaces showed a normal distribution in the Qinghai-Tibet Plateau. The radiant energy had the maximum in the south-west of China throughout the year. The solar radiation reaches the maximum value in summer, and the ground surface temperature has the maximum from June to August [29]. The conversion ratio coefficient is different on different underlying surfaces; therefore, the conversion ratio from total solar radiation to long-wave radiation on the ground surface and the order of ground surface temperature is grassland > wetland > forest > bare land > water body > glacier. The MMGST obtained by the simplified calculation model has a good fit with that of MODIS data for the entire year.
The calculation results of MARGST from the simplified model have a correlation coefficient higher than 0.9 with remote sensing data. After the modification, the calculation results reduce the deviation and mean square root. The ground surface temperature on different underlying surfaces can be inverted by the conversion ratio coefficient and modification coefficient, which are determined by the regional latitude and altitude. The horizontal ground surface temperature was calculated by solar radiation; however, the absorbance of solar radiation differs with the slope and aspect of different underlying surfaces. In this paper, the calculation results of the ground surface temperature have a limitation regarding the slope and aspect. For further study, other environmental factors, such as air temperature, wind speed, soil moisture, precipitation, land cover and transpiration, should be considered in the calculation model to obtain better results.

Conclusions
In this paper, the radiant ground surface temperature on different underlying surfaces was calculated based on the energy balance of ground surface radiation. A modified calculation model was built to precisely reflect the ground surface temperature on different underlying surfaces. Some conclusions are summarized as follows.
(1) In the Qinghai-Tibet Plateau, the absorbance capacity of solar radiation differs on different underlying surfaces, which represents the ground surface temperature. The reasons are mainly associated with the difference in the conversion ratio from total solar radiation to ground surface long-wave radiation on different underlying surfaces.
Based on the calculation model, the order of conversion ratio coefficient is distributed as follows: grassland > wetland > forest > bare land > water body > glacier. (2) A simplified calculation model of ground surface temperature was built on different underlying surfaces, which considered latitude, altitude and solar radiation. The ground surface temperature can be approximately solved by the conversion coefficient and correction coefficient on different underlying surfaces. From the calculated results of MMGST and MARGST, the modified model has a high correlation with the MODIS data, which is higher than 0.85, and the RMSE and deviation were improved by the modified model. The inter-annual variation in MMGST has a unimodal curve trend with rising and falling. (3) The simplified and modified model presented in this paper can calculate the ground surface temperature at every hour of every day. By comparison with the MODIS data, the calculated results are found to have high precision. The application of the model will be expanded to other regions with different altitudes. The existing temperature calculation model mainly focuses on the Qinghai-Tibet Engineering Corridor, which is based on monitoring data. The calculation results represent the given condition and not all the underlying surfaces. The model presented in this manuscript has a higher resolution and wider application potential in ground surface analysis.