Warming E ﬀ ort and Energy Budget Di ﬀ erence of Various Human Land Use Intensity: Case Study of Beijing, China

: Human land use intensity a ﬀ ects the surface energy balance by changing the biogeophysical parameters. This study used Moderate Resolution Imaging Spectroradiometer remote sensing data and surface energy balance algorithms to quantify changes in surface energy budgets corresponding to changes in land use in Beijing from 2000 to 2015. Land use was reclassiﬁed by considering land use intensity. The di ﬀ erence in the latent heat ﬂux (LE) and net radiation (R n ) (LE − R n ) expressed the warming or cooling e ﬀ ect. The results showed that: (i) The increasing trend of net longwave radiation in Beijing o ﬀ set the decreasing trend of net shortwave radiation. The R n changed slightly, while the LE and LE − R n showed a signiﬁcant increase of 0.55 and 0.56 W / (m 2 · year), respectively. The ﬁndings indicated that considering only radiative forcing, or even R n , was not enough to measure the impacts of land use change on the energy budget. (ii) The order of R n , LE, and LE − R n values from high to low were natural and seminatural areas, cropland, mixed pixel areas, urban expansion areas, and old urban areas. Compared with natural and seminatural areas, the changing LE − R n trend in the other four land use types decreased with the increase in human impact intensity, indicating that human activities weakened the positive change trend of LE − R n and increased the warming e ﬀ ect. (iii) Although the temporal trend of LE increased in Beijing from 2000 to 2015, the e ﬀ ect of R n on LE − R n was greater than that of LE, especially in the four land use types a ﬀ ected by human activities. The results for surface temperature in various land use types conﬁrmed this point. This study highlights the energy budget di ﬀ erences of various land use types a ﬀ ected by human activities. It makes an important contribution to understanding the urban heat island e ﬀ ect from a biogeophysical perspective.


Introduction
Land use change (LUC), resulting from human activities, is one of the main driving factors of global climate change [1,2]. LUC changes the surface biogeochemical cycle and biogeophysical mechanisms that affect the regional and global climate systems [3][4][5][6]. At present, most studies emphasize the impact of the carbon source and sink aspect of the LUC system [7][8][9][10], but they ignore the effect of biogeophysical mechanisms on surface energy. Therefore, studying the impact of LUC on the surface energy budget is necessary to correctly understand the climate regulation effect of LUC and urban heat island effect of urban expansion [11].
The exchange process of surface radiation reflects the energy transmission in the ground-air coupling mechanism [12][13][14]. As a link between land-air energy conversion, the surface energy balance is an important intermediate mechanism for temperature changes caused by LUC [15]. LUC affects radiation and the energy budget by changing biogeophysical parameters, such as surface albedo and latent heat flux (LE), resulting in temperature changes [16][17][18][19][20][21]. As one of the important forcing factors, the albedo measures the reflection of solar radiation from different land use types, and thus can change the surface energy balance [22]. It has been found that the surface albedo can act on the planetary albedo at the top of the troposphere; if the global planetary albedo changes by 5%, the global surface temperature will change by~1 • C [23]. Additionally, soil moisture and vegetation can affect the surface LE. More vegetation may increase the absorption of solar radiation on the ground, but if the extra energy is used in evapotranspiration it may cause a cooling effect [24]. Therefore, when studying the feedback effect of LUC on the surface energy budget, it is necessary to start from albedo-led radiative forcing (or energy intake or net radiation (R n )) and LE-led energy expenditure to identify their interactions.
Land use and land management changes affect the energy budget and result in obvious surface temperature changes corresponding to various land use types [6,25]. In recent years, many researchers have studied the energy balance and constructed a series of methods for estimating surface energy budget processes under various land use types or conditions, such as the Surface Energy Algorithms for Land model and Surface Energy Balance System model [16,17,26,27]. In China, it has been confirmed that these methods can describe the surface energy balance process more accurately and quantitatively, and can highlight the interaction and role of various energy elements at a regional scale [28,29]. Some studies compared the radiative forcing corresponding to the albedo changes caused by typical land use types [30,31]; other studies analyzed the differences in energy budgets and evapotranspiration effects of different types of underlying surfaces based on local-scale urban meteorological parameterization schemes, they emphasized the differences in the energy budget under the complex urban surfaces in cities [18,32].
Using energy balance parameterization schemes combined with remote sensing to explore the energy and temperature feedbacks of LUC has also attracted research attention. Lee andGe et al. (2011, 2019) introduced remote sensing data and methods to compare the energy and temperature feedbacks of typical LUC to radiant and non-radiative energy contributions in China [33,34]. Duveiller et al. (2018aDuveiller et al. ( , 2018b, 2020) also published a series of reports analyzing the role of biogeophysical factors, combined with remote sensing, to pioneer a framework for evaluating the energy and temperature effect of LUC [35][36][37]. These studies provide us with some information, but because of the diversity of land use and the complexity of the changes of surface parameters related to land use, understanding how the changes of surface energy intake and expenditure correspond to LUC still faces some challenges [38]. There are no practical guidelines for specific regions, and in and around urban areas, in particular, the impacts of human activities on the underlying surface are variable and complex [13,[39][40][41]. The final impacts of LUC on the interaction between energy factors, under different intensities of human activities, remain inconclusive.
On the basis of previous research, we combined MODIS (Moderate Resolution Imaging Spectroradiometer) remote sensing products and land energy balance algorithms to calculate the surface energy intake and expenditure corresponding to the underlying surface changes in Beijing The results of mutual forcing between R n and LE were also analyzed. Additionally, on the basis of the study results, we explored the relationship between temperature (LST), LE, and R n .

Study Data
The land use data for 2000 and 2015 for Beijing, with a spatial resolution of 1 km used in this study, came from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resdc.cn/), as did the annual interpolated temperature data for Beijing. The water vapor pressure data were obtained from the Scientific Data Sharing Center of the China Meteorological Administration (http://data.cma.cn/). This study mainly used remote sensing products with MODIS as the source data, which have mature algorithms and theoretical supports, and have been widely used in the world. The solar radiation data were downloaded from the Environmental Ecology Laboratory of the National University of Seoul, South Korea (http://environment.snu.ac.kr/). As shown in Table 1, the remote sensing inversion data for land surface albedo (Albedo), LST, emissivity (Emissivity), and LE all came from MODIS data products provided by NASA (https://modis.gsfc.nasa.gov/), over the time span 2000-2015. We unified these remote sensing data into annual average data with a spatial resolution of 1 km to analyze the surface energy budget of Beijing from 2000 to 2015.

Land Use Data Reclassification
The land use datasets for 2000 and 2015 in Beijing were reclassified based on the effects of human activities and the natural resilience of the land use types [39]. First, three land use types with obvious differences in underlying surface cover were classed as: urban areas, cropland, and natural and seminatural areas for the two datasets. Second, areas that were always urban from 2000 to 2015 were defined as old urban areas (OU); areas that were always cultivated land were defined as cropland (CP); and unchanged natural or seminatural areas, the urban expansion areas, and areas with changes in land use type that were not urban expansion were defined as pure pixel areas (PP), urban expansion areas (UE), and mixed pixel areas (MP), respectively. According to their natural resilience, the land use types after reclassification were ordered from high to low as: PP > CP > MP > UE > OU. Conversely, based on the intensity of human activities, the order of land use types from high to low was: OU > UE > MP > CP > PP, so that the land use types reflected the strength of human activities to some extent ( Figure 1a). Overall, the area of pixels with changed land use was 789 km 2 , accounting for 4.8% of the total area. As shown in Figure 1b, the total area of the OU was 1259 km 2 ; the UE around the OU covered an area of 517 km 2 ; and the MP covered a total area of 272 km 2 . The CP and PP occupied most of the land area in Beijing, at 4619 km 2 and 9743 km 2 , respectively.

Surface Energy Balance and Warming Effect
Based on the surface energy balance algorithms of land uses, this study calculated the changes in the energy budget in Beijing and analyzed the effect of different land uses on the energy budget. On the basis of the five land use types (OU, UE, MP, CP, and PP), we explored the differences and relationships in energy intake and expenditure parameters under various land use types ( Figure 2). According to the studies of Bastiaanssen et al. (1998) and Huang et al. (2018) [16,28], the Rn was calculated as:

Surface Energy Balance and Warming Effect
Based on the surface energy balance algorithms of land uses, this study calculated the changes in the energy budget in Beijing and analyzed the effect of different land uses on the energy budget. On the basis of the five land use types (OU, UE, MP, CP, and PP), we explored the differences and relationships in energy intake and expenditure parameters under various land use types ( Figure 2).

Surface Energy Balance and Warming Effect
Based on the surface energy balance algorithms of land uses, this study calculated the changes in the energy budget in Beijing and analyzed the effect of different land uses on the energy budget. On the basis of the five land use types (OU, UE, MP, CP, and PP), we explored the differences and relationships in energy intake and expenditure parameters under various land use types ( Figure 2). According to the studies of Bastiaanssen et al. (1998) and Huang et al. (2018) [16,28], the Rn was calculated as:  [16,28], the R n was calculated as: Land 2020, 9, 280 where R ns is net shortwave radiation; R nl is net longwave radiation; R ↓ s is the downward shortwave radiation; R ↑ s is the upward shortwave radiation, and its unit of radiation is W/m 2 ; and α is the surface albedo. The net longwave radiation was calculated by the following formula: where R ↓ l is downward longwave radiation; R ↑ l is upward longwave radiation, and its unit of radiation is W/m 2 ; ε s is the surface emissivity; σ is the Boltzmann constant, which is 5.67 × 10 −8 W/m 2 /K 4 ; T a is the air temperature; and T s is the land surface temperature and its unit is K. ε a is the emissivity of air, and its calculation formula is: where e a is the saturated water vapor pressure and its unit is hpa.
In keeping with the study of Anderson-Teixeira et al. (2012), this study used −R n to denote the energy intake forcing caused by biogeophysical factor factors such as albedo [42]. According to the calculation of radiative forcing of albedo (Tang et al., 2020), the negative sign indicates that the R n decreases with the increase of albedo. Meanwhile, LE represents the forcing of energy expenditure by consuming R n . Therefore, LE−R n in this study represents the final result of the energy intake and expenditure forcing balance caused by LUC. In the forcing balance, larger values of LE−R n reduce the energy from the sensible heat flux and the soil heat flux. Therefore, less energy is available to heat the atmospheric and surface temperature, corresponding to a cooling feedback or a weak warming feedback. Conversely, lower values of LE−R n equate to a stronger warming feedback.

Spatial Analysis and Statistical Analysis
During the process of specific analysis, we also used the zonal statistical method of spatial statistics to summarize the mean values of all pixels within a land use type. We also used the linear regression, trend fitting, and t test methods (a significance test with small sample) of classical statistics to analyze the temporal changes of different energy factors for various land use types from 2000 to 2015.

Net Longwave and Shortwave Radiation
The difference between the first and the last periods of data was used to calculate the spatial variation in the net shortwave and net longwave radiation from 2000 to 2015. Overall, the change in net longwave radiation in most areas of Beijing was positive ( Figure 3a). The multi-year trend values of R nl and R ns were 0.33 W/(m 2 ·year) and −0.34 W/(m 2 ·year), and the multi-year average values of R nl and R ns were 142.94 W/m 2 and −83.91 W/m 2 .
In terms of spatial distribution, the size of the net longwave radiation increment was generally large in the northeast and smaller in the southeast (Figure 3a). The regional distribution of negative differences was concentrated in the OU and CP regions in the southeast, while the high positive differences were mainly distributed in the PP region of the northeast, where R nl rose to 50 W/m 2 .
In contrast, Figure 3b shows that the net shortwave radiation R ns in many areas of Beijing decreased, and its high value areas were mainly distributed in the OU area, with a maximum value of 61 W/m 2 .
The variation of R nl showed that PP had the largest multi-year average of 146.79 W/m 2 , followed by OU, CP, UE, and MP. Although the surface parameters had a great influence on R nl , this result showed that R nl did not correspond to obvious rules for the land use types under different intensities of human activities. For R ns , when human activities increased, the corresponding R ns decreased. As shown in Figure 3b, the multi-year mean of R ns was negative, and the order from large to small was CP > PP > MP > UE > OU. The variation of Rnl showed that PP had the largest multi-year average of 146.79 W/m 2 , followed by OU, CP, UE, and MP. Although the surface parameters had a great influence on Rnl, this result showed that Rnl did not correspond to obvious rules for the land use types under different intensities of human activities. For Rns, when human activities increased, the corresponding Rns decreased. As shown in Figure 3b, the multi-year mean of Rns was negative, and the order from large to small was CP> PP> MP> UE> OU.

Changes in Net Radiation
The net surface radiation is the net energy obtained from the net longwave and shortwave radiation processes, and therefore it is an important indicator to measure the regional energy intake. The average annual surface radiation in Beijing from 2000 to 2015 was 59.03 W/m². Although the overall Rn in 2015 was higher than that in 2000, the changes in Rn over the years were very small. Thus, the variation in Rn was not statistically significant (P = 0.751, Figure 4). Moreover, after superimposing the net longwave and shortwave radiation, the Rn changed significantly in some areas, showing increases in the northeast but decreases in the southeast (Figure 4a).

Changes in Net Radiation
The net surface radiation is the net energy obtained from the net longwave and shortwave radiation processes, and therefore it is an important indicator to measure the regional energy intake. The average annual surface radiation in Beijing from 2000 to 2015 was 59.03 W/m 2 . By superimposing the net longwave and shortwave radiation, the R n changed significantly in some areas, showing increases in the northeast but decreases in the southeast (Figure 4a). Although the overall R n in 2015 was higher than that in 2000, the changes in R n over the years were very small. Thus, the variation in Rn was not statistically significant (p = 0.751, Figure 5). Although the difference between net shortwave and net longwave radiation corresponding to different land use types was not dissimilar, a clear difference in Rn was observed between land use types (Figure 4b). Considering that the value of Rn decreases as the impact of human activities on the surface increases, the Rn of PP was much larger than that of other land use types ( Figure 5). Its multiyear average was 63.54 W/m 2 , followed by CP and MP. Moreover, the Rn of UE and OU, which were more affected by human intervention, were always smaller than other types of land use. The value of OU was the smallest, and the multi-year average was 44.19 W/m 2 . The trends of the land use types (OU, UE, and MP) that were greatly affected by human activities during the study period were all negative. The values of land use types more intensively affected by human factors were also lower than PP and CP. Compared with PP, the change trends of Rn in the four land use types with various human influences were all negative, with a decreasing trend. However, the changes were not statistically significant (P > 0.1). Although the difference between net shortwave and net longwave radiation corresponding to different land use types was not dissimilar, a clear difference in R n was observed between land use types ( Figure 4b). Considering that the value of R n decreases as the impact of human activities on the surface increases, the R n of PP was much larger than that of other land use types ( Figure 5). Its multi-year average was 63.54 W/m 2 , followed by CP and MP. Moreover, the R n of UE and OU, which were more affected by human intervention, were always smaller than other types of land use. The value of OU was the smallest, and the multi-year average was 44.19 W/m 2 . The trends of the land use types (OU, UE, and MP) that were greatly affected by human activities during the study period were all negative. The values of land use types more intensively affected by human factors were also lower than PP and CP. Compared with PP, the change trends of R n in the four land use types with various human influences were all negative, with a decreasing trend. However, the changes were not statistically significant (p > 0.1).  Although the difference between net shortwave and net longwave radiation corresponding to different land use types was not dissimilar, a clear difference in Rn was observed between land use types (Figure 4b). Considering that the value of Rn decreases as the impact of human activities on the surface increases, the Rn of PP was much larger than that of other land use types ( Figure 5). Its multiyear average was 63.54 W/m 2 , followed by CP and MP. Moreover, the Rn of UE and OU, which were more affected by human intervention, were always smaller than other types of land use. The value of OU was the smallest, and the multi-year average was 44.19 W/m 2 . The trends of the land use types (OU, UE, and MP) that were greatly affected by human activities during the study period were all negative. The values of land use types more intensively affected by human factors were also lower than PP and CP. Compared with PP, the change trends of Rn in the four land use types with various human influences were all negative, with a decreasing trend. However, the changes were not statistically significant (P > 0.1).

Analysis on the Change of the Surface Energy Expenditure
As shown in Figure 6a, the LE changes in most parts of Beijing were positive from 2000 to 2015, except in central Beijing. Expanding outward from the central city, the LE difference increased from low to high. The value range of LE corresponding to the land use type under the influence of human activity showed marked differences (Figure 6b). Most of the northeast part of region, far from Beijing city, showed high values of LE. In Beijing, the multi-year average LE was 31.8 W/m 2 . PP had the largest multi-year mean value, at 37.61 W/m 2 , followed by CP, MP, UE, and OU.

Analysis on the Change of the Surface Energy Expenditure
As shown in Figure 6a, the LE changes in most parts of Beijing were positive from 2000 to 2015, except in central Beijing. Expanding outward from the central city, the LE difference increased from low to high. The value range of LE corresponding to the land use type under the influence of human activity showed marked differences (Figure 6b). Most of the northeast part of region, far from Beijing city, showed high values of LE. In Beijing, the multi-year average LE was 31.8 W/m². PP had the largest multi-year mean value, at 37.61 W/m 2 , followed by CP, MP, UE, and OU.

Comparison of Net Radiation and Latent Heat Flux
The spatial distribution of LE−R n was almost the opposite to the LE distribution (Figure 8a). The values of LE−R n in downtown Beijing and its surroundings were negative, while the distribution of positive values was dominant in the remaining areas. However, there were scattered grid cells with negative values. The multi-year mean value of LE−R n was negative, −25.92 W/m 2 , and the order of LE−R n from large to small was PP > CP > MP > UE > OU (Figure 8b). This order was consistent with the strength of natural resilience, and contrasted with the intensity of human activities. This indicates that LE−R n will gradually decrease, and the warming feedback will become stronger, as human activities increase.

Comparison of Net Radiation and Latent Heat Flux
The spatial distribution of LE−Rn was almost the opposite to the LE distribution (Figure 8a). The values of LE−Rn in downtown Beijing and its surroundings were negative, while the distribution of positive values was dominant in the remaining areas. However, there were scattered grid cells with negative values. The multi-year mean value of LE−Rn was negative, −25.92 W/m 2 , and the order of LE−Rn from large to small was PP > CP > MP > UE > OU (Figure 8b). This order was consistent with the strength of natural resilience, and contrasted with the intensity of human activities. This indicates that LE−Rn will gradually decrease, and the warming feedback will become stronger, as human activities increase.

Comparison of Net Radiation and Latent Heat Flux
The spatial distribution of LE−Rn was almost the opposite to the LE distribution (Figure 8a). The values of LE−Rn in downtown Beijing and its surroundings were negative, while the distribution of positive values was dominant in the remaining areas. However, there were scattered grid cells with negative values. The multi-year mean value of LE−Rn was negative, −25.92 W/m 2 , and the order of LE−Rn from large to small was PP > CP > MP > UE > OU (Figure 8b). This order was consistent with the strength of natural resilience, and contrasted with the intensity of human activities. This indicates that LE−Rn will gradually decrease, and the warming feedback will become stronger, as human activities increase.  From 2000 to 2015, the LE−R n in various land types showed an upward trend (0.56 W/(m 2 ·year), indicating that the feedback effect of LUC on regional warming was weakening. However, compared with the trend of PP, the LE−R n trends in other four types were relatively small. The order from small to large was OU < MP~UE < CP (Figure 9). This also shows that the more intense the human activity, the greater the corresponding warming feedback effect.
Land 2020, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/land From 2000 to 2015, the LE−Rn in various land types showed an upward trend (0.56 W/(m 2 •year), indicating that the feedback effect of LUC on regional warming was weakening. However, compared with the trend of PP, the LE−Rn trends in other four types were relatively small. The order from small to large was OU < MP ~ UE < CP (Figure 9). This also shows that the more intense the human activity, the greater the corresponding warming feedback effect. To find out whether LE or Rn had a greater impact on net radiation (LE−Rn), we compared the three relationships among the three energy factors. Compared with Rn, the differences of LE in the five land cover types were more obvious, and the correlation between Rn and LE−Rn was greater than that between LE and LE−Rn during the study period ( Figure 10). This indicates that, in terms of the final forcing balance, the LE−Rn in land use type was more sensitive to the change of Rn indicating the importance of energy intake on warming or cooling effects. To find out whether LE or R n had a greater impact on net radiation (LE−Rn), we compared the three relationships among the three energy factors. Compared with R n , the differences of LE in the five land cover types were more obvious, and the correlation between R n and LE−R n was greater than that between LE and LE−R n during the study period ( Figure 10). This indicates that, in terms of the final forcing balance, the LE−R n in land use type was more sensitive to the change of R n indicating the importance of energy intake on warming or cooling effects.
Land 2020, 9, x; doi: FOR PEER REVIEW www.mdpi.com/journal/land From 2000 to 2015, the LE−Rn in various land types showed an upward trend (0.56 W/(m 2 •year), indicating that the feedback effect of LUC on regional warming was weakening. However, compared with the trend of PP, the LE−Rn trends in other four types were relatively small. The order from small to large was OU < MP ~ UE < CP (Figure 9). This also shows that the more intense the human activity, the greater the corresponding warming feedback effect. To find out whether LE or Rn had a greater impact on net radiation (LE−Rn), we compared the three relationships among the three energy factors. Compared with Rn, the differences of LE in the five land cover types were more obvious, and the correlation between Rn and LE−Rn was greater than that between LE and LE−Rn during the study period ( Figure 10). This indicates that, in terms of the final forcing balance, the LE−Rn in land use type was more sensitive to the change of Rn indicating the importance of energy intake on warming or cooling effects.

Discussion
Since the 20th century, Beijing has been in the stage of rapid urbanization. Correspondingly, the land use pattern of Beijing has undergone obvious changes, especially the most significant changes in cultivated land and urban land, and urban construction has taken up a lot of cultivated land [19]. Vegetation growth conditions of different land use types are different, and the differences in each land use type are reflected in biogeophysical parameters such as albedo, resulting in different energy budget effects [18,32]. Many researchers have explored the climatic effects of LUC from the perspective of radiative forcing [43,44]. Previous studies have shown that if only the radiative forcing corresponding to albedo is considered, the radiative forcing at the top of the atmosphere would increase with the intensity of human activity [39]. Simply considering albedo's radiative forcing only includes shortwave radiation, which has a limited impact on the overall energy budget and on changes in various parameters within the energy balance mechanism that regulate local energy [32].
The results of our study showed that the regularity of R ns in various land use types was not immediately obvious, but after superimposing R nl , the R n did show an obvious regularity corresponding to various land use types as the intensity of human activities increased. Other studies have also confirmed that there is a large non-radiation impact on regional energy and temperature [33,34]. Our results further found that LE−R n was more sensitive to R n than to LE, namely, the energy regulation effect of R n is greater than LE in Beijing.
Additionally, the energy budget itself will be directly affected by the LST during the calculation process (by upward longwave radiation), and the energy expenditure process will also react to the LST, including LE, which can cause changes to LST. Here, we explored the relationship between energy and temperature and found that the relationships between LST, R n , and LE were negative while the relationships between LST and LE−R n were positive in all the five land use types ( Figure 11). Overall, the relationship between R n and LST was the closest, followed by LE and LST, which passed the significance test of P = 0.1, while the relationship between LST and LE−R n did not pass the significance test. However, with the increase of natural resilience, the relationships between LE−R n and LST became weaker, and LE showed a more significant relationship with LST (p < 0.1), indicating that LE had a strong regulatory effect on LST in the areas with natural attributes. This is consistent with our previous model simulation [45].
When analyzing the urban heat island effect, the role of R n on LST or various land use types cannot be ignored. In terms of specific land use types, LE−R n has a significant impact on LST in OU, UE, and MP, which have been strongly impacted by human activities, while LE plays an important role in LST in CP and PP, which have natural or seminatural resilience. That is to say, the energy budget item that dominates the LST change is not immutable. This implies that when studying urban heat island, problems may arise with the spatial scale and extent of the study area. Additionally, the impact area of the urban heat island is far greater than the urban physical boundary [46,47]. Therefore, when we analyze the spatial differences in the urban thermal environment, the setting of its boundary extent and land use type will affect not only the heat island intensity but also the intrinsic energy budget corresponding to its LST change. When analyzing the urban heat island effect, the role of Rn on LST or various land use types cannot be ignored. In terms of specific land use types, LE−Rn has a significant impact on LST in OU, UE, and MP, which have been strongly impacted by human activities, while LE plays an important role in LST in CP and PP, which have natural or seminatural resilience. That is to say, the energy budget item that dominates the LST change is not immutable. This implies that when studying urban heat island, problems may arise with the spatial scale and extent of the study area. Additionally, the impact area of the urban heat island is far greater than the urban physical boundary [46,47]. Therefore, when we analyze the spatial differences in the urban thermal environment, the setting of its boundary extent and land use type will affect not only the heat island intensity but also the intrinsic energy budget corresponding to its LST change.

Conclusions
On the basis of the different human impacts, this study divided the land use types of Beijing into five types-OU, UE, MP, CP, and PP-and analyzed the differences in Rn and LE corresponding to various land use types and changes from 2000 to 2015. The difference in the relationship between LST and energy intake and expenditure on different underlying surfaces was also discussed. This study reached the following conclusions: From 2000 to 2015, under the impacts of different human activity intensities and natural resilience, the energy budgets of various land use types showed very obvious differences. The order of Rn, LE, and LE−Rn values were all PP > CP > MP > UE > OU. The trends of LE−Rn on the four land use types affected by human activities were all less than that of PP, showing that human activities weaken the trend of LE−Rn and increased the warming effect.
From the perspective of energy balance, the role of Rn in energy regulation for a warming or cooling effect was greater than LE. This study also found that the main energy regulation factors on LST varied in different land use types. With the increasing intensity of human activities in relation to land use, the warming/cooling effect of LE−Rn gradually increased.

Conclusions
On the basis of the different human impacts, this study divided the land use types of Beijing into five types-OU, UE, MP, CP, and PP-and analyzed the differences in R n and LE corresponding to various land use types and changes from 2000 to 2015. The difference in the relationship between LST and energy intake and expenditure on different underlying surfaces was also discussed. This study reached the following conclusions: From 2000 to 2015, under the impacts of different human activity intensities and natural resilience, the energy budgets of various land use types showed very obvious differences. The order of R n , LE, and LE−R n values were all PP > CP > MP > UE > OU. The trends of LE−R n on the four land use types affected by human activities were all less than that of PP, showing that human activities weaken the trend of LE−R n and increased the warming effect.
From the perspective of energy balance, the role of R n in energy regulation for a warming or cooling effect was greater than LE. This study also found that the main energy regulation factors on LST varied in different land use types. With the increasing intensity of human activities in relation to land use, the warming/cooling effect of LE−R n gradually increased.
This study reveals the clear differences in the energy budget correspondingly to various land use types under the impact of human activities. It makes an important contribution to understanding the urban heat island effect from the perspective of biogeophysical mechanisms. Furthermore, this study suggests that for urban planners and decision makers, urban vegetation or water bodies can be used to regulate evapotranspiration during the urban design process to mitigate the urban heat island effect. At the same time, adjusting net radiation factors may have more obvious effects on reducing the heat island effect.