4.1. Characteristics of sUHI Spatiotemporal Variation
The LST
R and other similar LST normalization methods are widely used to analyze the sUHI effect [
42,
43,
44,
45]. Based on the pixel number trend of sUHI intensity grades, the characteristics of sUHI areas and intensity development in the area within the 6th Ring Road area of Beijing in 2000–2017 were investigated. The results showed that there were different stages in the sUHI development. In order to get the actual physical value of sUHI intensity and verify the results based on LST
R, this study also estimated the sUHI intensity based on the LCZ scheme, which was calculated as 2–6 K.
In terms of the evolution of sUHI intensity grades based on LSTR, for the sUHI area, 2006 was the turning point. From 2000 to 2006, the sUHI area increased, but it decreased after 2006. However, the turning point of sUHI intensity was in 2009. The sUHI intensity increased from 2000 to 2009, after which it remained stable. Nevertheless, the turning point of 2006, which was different from the LSTR scheme, was observed in the time series of sUHI intensity based on the LCZ scheme. The decreasing trend of sUHI intensity from 2006 to 2009 was consistent with that calculated based on the LSTR method. However, there was still a declining trend of sUHI intensity based on the LCZ scheme from 2009 to 2015. The different results indicate that the change in sUHI intensity might get lost in the division of sUHI intensity grade and therefore cannot accurately reflect the change in sUHI intensity. The rising trend of sUHI intensity in 2000–2006 was contrary to that of the estimation based on the LSTR method. To understand this, the pixel numbers of the medium and strong sUHI in 2000–2006 were further counted. The overall trend of the sum of the two sUHI grades increased, which is consistent with the trend of the LCZ-based sUHI intensity variation for the same period. Therefore, it can be concluded that when the LSTR method was used to obtain the medium and strong sUHI sequences, part of the increasing signal was masked by the overall downward trend in the average process.
According to the turning points of sUHI area and intensity, the development of sUHI in Beijing can be divided into three stages. (1) 2000–2006: both area and intensity of sUHI increased; (2) 2006–2015: both area and intensity of sUHI decreased; (3) 2015–2017: sUHI intensity increased (which was uncertain because of the short period), but the area decreased. The overall trend of an initial increase in sUHI intensity followed by a decrease presented in our study is consistent with other research on sUHI intensity variation in Beijing. Liu et al., (2020) investigated the magnitude and spatial patterns of sUHI from 2003 to 2018 in a similar area of Beijing [
46]. The turning year of 2009, before which the sUHI intensity increased and after which it decreased, was also found in their study. Meng et al., (2018) also found an upward trend of sUHI intensity from 2003 to 2008 [
47]. Combining all the studies above with ours, the sUHI intensity in Beijing increased from 2000 and then presents a downward trend. The turning year may be around 2006–2009. Furthermore, the sUHI intensity increased from 2016 to 2018 in the study by Liu et al., (2020) [
46], which is similar to our finding. Therefore, the sUHI intensity may present upward trend recently, after 2015 or 2016, which needs further research.
For the full research period (2000–2017), the spatiotemporal pattern of sUHI change shows that areas with increased sUHI intensity were mainly concentrated in a parabola belt area of the westward opening outside the 4th Ring Road, where significant urbanization occurred [
48], especially in the southeast area. In addition, the area in the southwest corner and the Capital International Airport Area in the northeast corner were also areas with significantly increased sUHI, while the sUHI intensity in the urban area within the 4th Ring Road showed a downward trend, with a slight increase in some areas. Since 2000, the sUHI intensity in the central area of Beijing decreased, or increased insignificantly, but it increased in peripheral areas. The spatial distribution of the Theil–Sen regression slope of LST
R in different periods reflected the spatiotemporal heterogeneity of sUHI intensity change. There was a general trend for sUHI enhancement in space, which was presented as the sUHI intensity increase from north to south, from east to west, and from periphery to center. The sUHI pattern changed before and after 2009. Before 2009, the sUHI mainly increased in the suburbs, developing from the inner to outer suburbs, and it decreased or slightly increased in the central urban area. After 2009, the sUHI intensity enhancement area moved to the central urban area, and sUHI decreased in a large area outside the 4th Ring Road. Nevertheless, the spatial patterns of sUHI change identified by Liu et al., (2020) [
46] are opposite from ours. The main reason may be the difference in the definition of sUHI intensity used in their study and ours, which is based on the LST difference between urban and rural area identified by the proportion of impervious surface, and LST
R grade, respectively. Moreover, a similar increasing trend is identified in the rural area, similar to the area detected in our study in the research of Li et al., (2020) [
49]. The research investigates the spatial pattern variation of Beijing from 2000 to 2013 using the difference between the LST of each pixel and the averaged of background areas as the sUHI intensity, which is similar to our LST
R based sUHI intensity definition. Therefore, the different method of calculating sUHI may lead to different or even opposite results for sUHI spatiotemporal patterns.
4.2. Relationship between LCZ, Population Density, and sUHI
The thermal differential of LCZs has been widely recognized [
50,
51,
52,
53,
54]. According to the analysis of the relationship between LCZ type and LST
R, there are significant differences in the sUHI intensity of different urban forms. The order of sUHI intensity of built types showed that the sUHI intensity was closely related to the height of buildings. In general, the lower the building height, the stronger the sUHI intensity, which is consistent with previous studies [
54,
55,
56]. For the same height level, the sUHI intensity in open zones was significantly lower than those in compact zones. For different heights, the difference of sUHI intensities in open zones was larger, which demonstrates that the change in sUHI intensity according to height is more obvious in an open environment.
The LCZ classification maps for different years allowed for analysis of the relationship between the change in urban spatial form and sUHI change. Wang et al., (2019) found the warming trend in the Pearl River Delta from 2000 to 2015 is related to conversion between and within LCZ types [
57]. Based on the sUHI intensity change caused by the LCZ conversion in 2000–2009, 2009–2017, and the overall performance in 2000–2017, specific LCZ changes would lead to a stable sUHI intensity change trend. If other built types were transformed into LCZ 4, or LCZ 3 was transformed into other built types (except LCZ 8), the sUHI intensity would be effectively reduced, and the sUHI intensity would significantly increase if other LCZ types were converted to LCZ 8. Research by Wang et al., (2019) also shows that the expansion of LCZ 8 lead to LST increase [
57]. An instability of sUHI intensity change caused by the LCZ transformation was also observed, especially in the transformation between built types. For example, the transformation of LCZ 2 into other built types (except LCZ 4) led to the decrease in sUHI intensity in 2000–2009 and the increase in 2009–2017.
Population has close relationship with UHI [
58,
59]. Zhang et al., (2013) used multiple linear regression to investigate the relationship between population density, which is extracted from statistical yearbooks, and other factors, and UHI in Shanghai [
60]. A positive effect of population density on UHI was found in this study. However, Du et al., (2016) found that there was no significant correlation between population density and UHI intensity in the Yangtze River Delta Urban Agglomeration based on population data from statistical yearbooks and the linear regression method [
61]. The statistical population data and the simple regression method may have led to the inconsistent results. In our study, based on the relationship between LST
R and the population density distribution map, increasing population aggravated the sUHI effect, and the impact varied with different population density ranges and regions. The influence was most significant when the population density was at a low level and began to increase. In that stage, the sUHI intensity increased sharply, but the sUHI was still weak (LST
R < 0.1). With the increase in population density, the direct impact of human activities on sUHI tended to be small, and the sUHI intensity approached or reached medium sUHI. After reaching a certain threshold (approximately 300 people/hm
2), the increase in population density no longer increased the sUHI intensity effectively. In terms of space, the relationship between sUHI intensity and human activities was the most significant in the suburbs outside the 4th Ring Road, with a high positive correlation, but it was not significant in the central urban area.
The change in population density affected the sUHI differently according to different stages and regions. The reason for this may be that during the increase in population density from extremely low to the first turning point, not only the artificial heat increased, but the underlying surface also changed dramatically as the vegetation was converted to construction land [
35]. These factors induced the LST to increase sharply, thus leading to the sUHI effect. From the first to the second turning point, urbanization occurred continuously. With the increase in population density, anthropogenic heat emissions, such as heat emissions from air conditioning in summer and transportation, gradually increase, which aggravates the sUHI effect [
62]. However, when the population density reached the second turning point, the load of energy consumption reached its peak, and the anthropogenic heat emissions do not increase as they are limited by the urban carrying capacity. In addition, the process of increasing population density is often accompanied by the process of land urbanization, with a consequent increase in the area of building surfaces and impervious surfaces, especially in the early urbanization process, i.e., the early stage of rapid population density increase. During this stage, a large amount of natural surface, which has a large specific heat capacity and can remove heat from the surface through evaporation or transpiration, is replaced by building surfaces and impervious surfaces. In contrast, the lower specific heat capacity and lack of evaporation and transpiration on building surfaces and impervious surfaces can exacerbate the sUHI. However, as land urbanization reaches a mature stage, the increase in population will no longer significantly increase building surfaces and impervious surfaces due to the limitations of plot ratio and land development, which means that the area and ratio of building surfaces and impervious surfaces will reach a relatively stable condition. Consequently, its enhancement of sUHI will also reach a relatively stable condition. That is, when the population density reaches a certain value, its effect of enhancing the sUHI effect by influencing the increase in building surfaces and impervious surfaces will not change significantly after the population density reaches a certain value. Therefore, the increase in population density no longer affects the sUHI intensity significantly.
Recently, GeoDetector has increasingly been used in the analysis of driving factors of UHI [
63,
64,
65], indicating that GeoDetector is appropriate for driving factor analysis for UHI. The quantitative method of GeoDetector helps in revealing the interpretability of LCZ, population, and their interaction to sUHI. Our driving factors analysis based on GeoDetector showed that both LCZ type and population density can effectively explain the spatial differentiation of sUHI, but the explanatory ability of LCZ type was significantly higher than that of population density. The interaction effect of the two factors can significantly improve the explanatory ability. Compared to population density, LCZ, which characterizes the urban space form, was the main factor influencing sUHI. Population density is an important factor in human activities and has a significant impact on sUHI, but its ability to explain the spatial differentiation of sUHI independently was poor. The interaction effect of the two factors was presented as enhancement and bivariate, thus indicating that the combination of urban space form and human activities can enhance the ability to explain the spatial differentiation of sUHI.
The dynamic analysis results show that LCZ conversion can contribute 30–40% of the spatial structure of sUHI intensity change, and the change in population density had a weak ability to explain the spatial structure of sUHI intensity change. This result shows that, compared to population density, LCZ conversion was the key factor leading to the sUHI intensity change. The phenomenon of urban spatial form being the main influencing factor, rather than the change in population density, has also been found in the center area of Kuala Lumpur [
66]. Nevertheless, the results of interaction detection showed that when considering both LCZ conversion and change in population density, the q values of the three time periods were significantly improved. The explanatory power was about 10% higher than that of LCZ variation only, and significantly higher than that of population density variation only. The results indicate that despite the weak ability of population density change to interpret sUHI change independently, population density was indispensable to explain the changes in the sUHI effect.