Abstract

The thermal environment is closely related to human well-being. Determinants of surface urban heat islands (SUHIs) have been extensively studied. Nevertheless, some research fields remain blank or have conflicting findings, which need to be further addressed. Particularly, few studies focus on drivers of SUHIs in massive cities with different sizes under various contexts at large scales. Using multisource data, we explored 11 determinants of surface urban heat island intensity (SUHII) for 1449 cities in different ecological contexts throughout China in 2010, adopting the Spearman and partial correlation analysis and machine learning method. The main results were as follows: (1) Significant positive partial correlations existed between daytime SUHII and the differences in nighttime light intensity and built-up intensity between cities and their corresponding villages except in arid or semiarid western China. The differences in the enhanced vegetation index were generally partially negatively correlated with daytime and nighttime SUHII. The differences in white sky albedo were usually partially negatively correlated with nighttime SUHII. The mean air temperature was partially positively correlated with nighttime SUHII in 40% of cases. Only a few significant partial relationships existed between SUHII and urban area, total population, and differences in aerosol optical depth. The explanation rates during daytime were larger than during nighttime in 72% of cases. The largest and smallest rates occurred during summer days in humid cold northeastern China (63.84%) and in southern China (10.44%), respectively. (2) Both the daytime and nighttime SUHII could be well determined by drivers using the machine learning method. The RMSE ranged from 0.49°C to 1.54°C at a national scale. The simulation SUHII values were always significantly correlated with the actual SUHII values. The simulation accuracies were always higher during nighttime than daytime. The highest accuracies occurred in central-northern China and were lowest in western China during both daytime and nighttime.

1. Introduction

In 2018, 55% of the world’s population lived in urban areas and this proportion is expected to reach 68% in 2050 [1]. Rapid, high-intensity urbanization has obviously changed the ecosystems and environment at urban, regional, and even global scales [2]. One of the most obvious urban climate features is the urban heat island effect [3]. This effect can directly and indirectly affect regional climate [4], energy use [5], air quality [6], urban hydrology [7], soil physicochemical properties [8], creature distribution and activities [9], and human health, comfort, and quality of life [10].

The approaches used for studying urban heat island (UHI) include weather station observations [11], fixed-point field measurements [12], mobile belt transect surveys [13], numerical modelling [14], and remote sensing monitoring, which have been widely accepted and adopted [15, 16]. Although land surface temperatures (LSTs) derived by remote sensing are not identical to above-ground air temperatures, they are closely related [17, 18]. Nonetheless, UHIs are called “surface urban heat islands (SUHIs)” when derived from remote sensing data in order to distinguish them from traditional UHIs analyzed using air temperatures [15, 1921]. Many determinants of SUHIs have been extensively studied and can be divided into seven categories [22]: land use and land cover types, surface biophysical conditions, landscape components and configurations, manners and intensities of human activity, meteorological conditions and geographical location, policy factors, and synthetic analyses of the abovementioned factors. The analytical methods applied to assess the relationship between drivers and SUHI intensity (SUHII) have mainly included Pearson’s correlation analysis, Spearman correlation analysis, comparative analysis, ordinary least squares regression, geographically weighted regression, regression tree model [16], and machine learning models [23, 24].

Nevertheless, certain issues remain unresolved. First, most previous studies only used simple correlation analysis [21, 2528] or linear [2730] or nonlinear regression [30] methods to analyze the relationships between the drivers and SUHIs. More effective and innovative methods should be considered, particularly when synthetically analyzing multiple factors. Additionally, several influencing factors should be considered together rather than simply analyzing the relationships between SUHI and a single influencing factor, such as vegetation index [20, 25, 30, 31], albedo [20, 25], built-up intensity [30, 32], total population [31, 32], and urban area [33, 34]. Second, the research has concentrated on large and mega cities [2528, 34, 35], while few studies focus on medium-sized and small cities. Third, certain disadvantages are associated with the use of previously proposed monitoring indices of SUHII [36]. For instance, certain studies defined SUHIIs as the differences in LSTs between a city or central city and its adjacent regions within a certain radius [27, 28, 3739]. Nonetheless, rapid urbanization usually occurs in areas surrounding cities, particularly large or mega cities in certain developing countries [40, 41]. Thus, biases can be introduced when monitoring SUHIIs. In addition, certain previous studies utilized differences in LSTs between cities or central cities and rural regions. In addition, rural regions were defined as zones within a certain distance from a city. However, the distance values in previous studies were too large for monitoring cities of different sizes [29, 30, 42], specifically because only regions that were sufficiently far away could be defined as reference villages. Therefore, few regions remained due to the dense distribution of cities in certain well-developed regions with dense populations. Moreover, certain studies have used the Gaussian fitting method [32, 43, 44], but this approach is only effective when the left field of the LST fits the Gaussian distribution after extraction of the background LST. However, this condition does not always exist [41]. For example, the method will fail when cold heat islands occur [28, 37] or when the outskirts or village regions have not been limited by the terrain [2729, 37], water bodies [29], or satellite towns. In addition, the impervious surface percentage should be considered when defining rural regions [29]. Fourth, although the ecological context has an obvious impact on SUHII, including water and heat combinations, vegetation characteristics, landforms, and human activities, most previous studies at large scale have not considered this factor [25, 28, 34, 35, 37, 45] but simply use administrative boundary data to determine regions of different contexts [27, 42] or identify these regions too coarsely [20, 26], with the exception of a few studies [29, 30, 46, 47]. Finally, certain research areas remain unstudied or have conflicting findings [22]. For instance, to our knowledge, unlike during the daytime, the impacts of vegetation on SUHIs during the nighttime are still in debate [21, 28, 42, 4851]. The relationships between albedo and daytime SUHIs are complex [27, 28, 52, 53]. Very few previous studies have analyzed impacts of aerosol optical depth (AOD) on SUHIs combined with other factors [38]. The effects of city size and total population on SUHIs have always been controversial [31, 33, 38, 54]. Different findings have been reported on the correlation between air temperature and SUHIs [2628, 34, 38, 42, 55]. Therefore, the objective of this study is to accurately study 11 associated determinants of SUHII for 1449 cities of different sizes in different ecological contexts throughout China in 2010.

2. Materials and Methods

2.1. Data

Monthly LST data at a 1 km resolution and enhanced vegetation index (EVI) data at 250 m resolution from Dec. 2009 to Nov. 2010 from MODIS/Aqua were downloaded from the Geospatial Cloud of Computer Network Information Center, Chinese Academy of Sciences website. Each pixel value of LST or EVI is a simple average of all the corresponding MYD11A1 LST pixels or maximum of all MYD09GQ EVI pixels collected within that month. The MYD11A1 product uses the split-window technique to derive LSTs. From the official MODIS website, we downloaded 8-day albedo data (MCD43B3) and quality data (MCD43B2) at 1 km resolution and daily aerosol data and quality data (MCD04_3K) at 3 km resolution from Dec. 2009 to Nov. 2010.

Digital elevation data with a 1 km resolution were taken from the Cold and Arid Regions Sciences Data Center at Lanzhou. Land use (Supplementary material Text S1) and population data from 2010, which had 1 km resolutions, were provided by the Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences. Population data are produced using the spatial interpolation method on statistical population data from each county. The 500 m resolution impervious percentage data for 2010 were provided by the Beijing City Lab (Supplementary material Text S2). The calibrated 1 km nighttime light data for 2010 were downloaded from the National Centers for Environmental Information website. Monthly 0.5° × 0.5° resolution temperature and precipitation data for Dec. 2009 to Nov. 2010 were downloaded from the Meteorological Data Website of China (Supplementary material Text S3).

Ecological function recognition data were provided by the Data Sharing Infrastructure of Urban and Regional Ecological Science. China is divided into three first-level, 50 second-level, 206 third-level, and 1434 fourth-level ecological function regions based on their landforms, water and heat combinations, and vegetation characteristics.

2.2. Methods
2.2.1. Regionalization of the Background Environment in China for SUHII Analysis

China was divided into four regions based on the ecological function recognition results for the first, second, and third levels (Figure 1) to analyze the potential different influencing factors of SUHII in each region. Region I mainly corresponded to the vast majority of northeastern China, which has climatic zones that range from cold temperate to midtemperate. The region is affected by monsoons and has four distinct seasons—warm and rainy summers and cold and dry winters—with long periods of melting ice and snow. The climate of region II is similar to that of region I. However, winters there are warmer and shorter than those in region I. Region III is the only zone that is not affected by monsoons in China. It is a typical arid temperate climate, with a mean annual precipitation (MAP) below 400 mm. Region III is subjected to cold winters and hot or warm summers. The vast majority of region IV belongs to the subtropical zone, and the rest is in the tropical zone. This region has a typical rainy and hot temperate climate, with MAP exceeding 800 mm and evergreen vegetation.

2.2.2. Definitions of Urban and Rural Regions

The urban land polygons were first aggregated at a distance of 1.1 km, which is sufficient to include most adjacent and scattered urban land polygons into the urban class and is able to separate the main city zones and satellite cities or two closely adjacent cities that should be considered separately according to normal human perceptions. Subsequently, only the urban land regions in these separate aggregated polygons remain. These polygons with remaining areas of larger than 6 km2 were considered, which included the vast majority of cities in the eastern regions in China, those with the densest populations and most well-developed economies. The corresponding rural zones were defined as the buffer zones of cities, with buffer distances between 5 and 10 km [37] and impervious percentages below 5% [29, 30]; however, the rural zones did not include waterbodies, regions with slopes exceeding 7.5°, or elevations that are 50 m greater than or less than the maximum or minimum elevations of the urban zones, respectively [30, 33, 42]. We found that 99.70% of the cities in China were located in plains regions with slopes less than 7.5° in 2010.

2.2.3. Calculation of SUHII and Its Drivers

The study timeframes were defined based on the definitions of meteorological seasons [56]. Winter ranges from December of the previous year to February of the current year, spring from March to May, summer from June to August, and autumn from September to November. The annual and seasonal mean values of daytime and nighttime LST, EVI, white sky albedo (WSA), and aerosol optical depth (AOD) in China were first calculated, and the differences in the abovementioned indices, population density (PD), nighttime light intensity (NL), and built-up intensity (BI)—defined as the percentage of the landscape taken by built-up land—in between the urban and rural regions were further calculated. Moreover, the urban area size (UAS) and total population (TP) values for these cities as well as the mean air temperature (MAT), total precipitation (P), and humidity index (HI) values at the centroids of the cities were calculated.

2.2.4. Exploration of the Relationships between SUHII and Its Drivers

The Spearman and partial correlation coefficients between SUHII and its drivers in China, as well as within its four regions, were calculated during daytime and nighttime for all four seasons and the entire year. Moreover, the explanation rates for the daytime and nighttime SUHII in each region and timeframe were derived using the stepwise regression method.

To simulate the SUHII based on drivers using the machine learning method, the cities were first randomly divided into training and testing sample sets at ratios of 4 : 1 for China and its four regions. Then, several simulation models were built, adopting four intelligence algorithms (Supplementary material Text S4), including the general regression neural network (GRNN), support vector machine based on genetic algorithm (SVM_GA), particle swarm optimization (SVM_PSO), and grid search method (SVM_GS), using three types of indices (Supplementary material Text S5), including all drivers, the drivers that were found to be significantly partially correlated, and those that were simply correlated with SUHII. Finally, the model was used to estimate the SUHII of the test samples, and the root-mean-square errors (RMSEs), correlation coefficients, and significance levels were calculated to indicate the simulation accuracy. The parameter optimization and simulation processes were repeated 49 times in order to avoid uncertainties during the process and to obtain stable results.

3. Results

3.1. Effects of Single Associated Determinants of SUHII

The differences in NL and BI between cities and villages (ΔNL and ΔBI) were both partially and positively correlated with SUHII in many cases () (Table 1). Nonetheless, they were almost insignificantly partially correlated with daytime SUHII in region III. Moreover, they were almost insignificantly partially correlated with daytime and nighttime SUHII in region I and daytime SUHII in region IV. UAS was only significantly positively partially correlated with daytime SUHII in region II in summer () and autumn () and with nighttime SUHII in region I in winter and spring. However, significant positive correlation coefficients generally existed, except in winter in northern China. Similarly, significant partial correlation did not exist or rarely existed between TP and SUHII during daytime and nighttime, respectively. In addition, positive correlations generally existed during daytime, except in winter in each region and in spring in region III. During nighttime, the differences in PD between cities and villages (ΔPD) exhibited significant positive partial correlations with nighttime SUHII in most cases in region IV and at the national scale although a significant positive Spearman correlation always existed, except in region III during winter, spring, and autumn and in region IV in summer. During daytime, significant partial correlations were minimal.

The differences in the EVIs between cities and villages (ΔEVI) were generally partially or negatively correlated with daytime and nighttime SUHII in many cases (). ΔEVI had a predominant cooling effect on daytime SUHII in region III. This was the only determinant that was significantly partially correlated with SUHII in summer (r = −0.79, ), autumn (r = −0.69, ), and the entire year (r = −0.66, ). However, during nighttime, ΔEVI and SUHII were not significantly partially correlated in region III, with the exception of summer (r = −0.50, ), in region I, with the exception of winter (r = −0.44, ), and the entire year (r = −0.48, ). ΔWSA was generally partially and negatively correlated with nighttime SUHII (). Both the partial and Spearman correlation coefficients were larger in northern China than in southern China. Nevertheless, the relationships were complex during daytime. For the most part, significant partial or Spearman correlations did not exist in southern China. ΔWSA was significantly partially negatively correlated with SUHII in winter in northern China (). ΔWSA was the only driver that was significantly partially correlated with daytime SUHII in winter in region I (r = −0.39, ). During nighttime, the differences in AOD between cities and villages (ΔAOD) were only negatively significant relative to SUHII at the national scale in spring () and autumn () and in region II in spring (). During daytime, significant positive partial relationships existed in spring in region I () and in autumn in region IV (). Nevertheless, significant negative partial relationships existed in summer at the national scale () and in autumn in region II ().

During nighttime, HI was significantly partially negatively correlated with SUHII in all four seasons at a national scale and in region II (), but this correlation was not observed or was seldom observed in other regions. In addition, significant negative correlations existed in all four seasons at a national scale and in region I and during certain seasons in other regions (). During daytime, significant positive partial correlations existed in region IV, with the exception of spring, which seldom existed in other regions. In addition, significant negative correlations did not exist in region III but did occur in all four seasons at a national scale () and during certain timeframes in other regions, especially II and IV. The patterns for P were similar to those for HI. However, the timeframes in which significant relationships existed in those regions may have differed. MAT was insignificantly partially correlated with nighttime SUHII in region I but was positively correlated at certain times at a national scale and in other regions, especially in region II. Only one partial negative correlation existed during summer nights in region IV (r = −0.16, ). In addition, significant negative and positive Spearman correlations existed in many cases at the national scale and in region IV, respectively. Few or complex significant Spearman correlations were identified for the other regions. During daytime, no significant partial correlation coefficients existed at the national scale, but positive Spearman correlation coefficients existed, except during summer. At the regional scale, significant partial or Spearman correlations seldom existed or were complex.

3.2. Combined Effects of Associated Determinants of SUHII

The explanation rates during daytime were larger than those during nighttime in 72% of the cases (Table 2). During daytime, they were larger in summer than in winter in northern China and smaller in southern China. The largest and smallest rates were found for summer days in regions I (63.84%) and IV (10.44%), respectively.

The SUHII values was considered to be best simulated by adopting the SVM-GA algorithm using the drivers that were found to be significantly partially correlated with SUHII values, after comparing four intelligence algorithms (Supplementary material Text S4 and Table S1) and three types of indices (Supplementary material Text S5 and Table S1). Both daytime and nighttime SUHII values were estimated at the national scale (Table 3). The simulation accuracy was always higher during nighttime than during daytime. The most accurate results were obtained for annual nighttime measurements (RMSE = 0.49°C, r = 0.71), while the least accurate results were obtained during daytime in summer (RMSE = 1.54°C, r = 0.55). Significant correlations always existed between the estimated and actual SUHII ().

At the regional scale, SUHII could also be simulated in each region for each timeframe. The simulation accuracy was always higher during nighttime than during daytime. The most accurate results were obtained for annual nighttime measurements in region I (RMSE = 0.38°C, r = 0.80, and ), while the least accurate measurements were recorded in region IV during summer days (RMSE = 2.18°C, r = 0.53, and ). Significant correlations existed between the estimated and actual SUHII () in 87.50% of cases. In northern China, the daytime SUHII results for autumn and summer were most accurate, while the results for spring and winter were less accurate. In southern China, the most accurate results also occurred in autumn, but the least accurate results were measured in summer. Seasonal variability was complex for nighttime SUHII in different regions. Spatially, the most accurate results for daytime SUHII usually occurred in region II although the results for winter in region IV also exhibited high accuracy. However, the least accurate results usually occurred in region III, with poor results also evident during summer in region IV and autumn in region I. Similarly, the most accurate nighttime SUHII values were associated with region II, except in winter and autumn, during which the most accurate results were associated with regions IV and I, respectively. The least accurate results also usually occurred in region III, except in summer, when they occurred in region IV.

4. Discussion

4.1. Associated Determinants for SUHII
4.1.1. Anthropogenic Heat Emissions

ΔNL is considered to be an indicator of anthropogenic heat flux. This factor can be converted into sensible heat flux or other energy components [57], thereby indirectly enhancing SUHII during daytime [28, 58] or nighttime [27, 28]. These findings are supported by the results of the current study. ΔNL exhibited a positive partial correlation with SUHII in most cases, except during daytime in regions I and III and nighttime in region I. These exceptions were due to the influence of ΔNL, which was much smaller than the influences of the other factors. For example, the key determinants were the cooling effects caused by vegetation in arid and semiarid western China and the low ΔALB in cold, humid northeastern China due to the presence of ice and snow over long periods in the villages there (Table 1). In addition, ΔNL was collinearly correlated with other factors in the regression, such as ΔBI [5961]. Previous studies also found that ΔNL was not significantly partially correlated with daytime SUHII in 419 large cities globally [28] and 32 large cities in China [27].

4.1.2. Built-Up Intensity

Increases in BI can cause reductions in latent heat flux and increases in sensible heat [30, 62]; decreases in sky views that reduce outgoing long-wave radiation [63]; the emergence of urban topographies, including the “canyon effect,” which increases the absorption of short-wave radiation [63]; and increases in surface roughness, which reduces boundary layer winds and hinders sensible loss [63], among other impacts. In addition, BI can cause the shade effect [64, 65]. Together, these factors generally enhance SUHII, except during daytime in cities in desert or xeric shrubland biomes [27, 30, 32, 65, 66], especially after stored heat is released at night [27]. The rules are complex for desert cities during daytime. Different conclusions have been drawn. In five cities in desert environments across the continental USA, uncharacteristic “U-shaped” horizontal gradients were found for the response of LST to BI [30]. In other words, LST decreased from the urban core to the outskirts and then increased again in the suburban to rural zones. A quadratic polynomial fit was used in that study because LST was not significantly linearly correlated with BI. A nonlinear but positive relationship was found in the central region of the city of Phoenix in the USA [66]. LST increased more rapidly with BI when its values were less than 0.4 and increased more slowly thereafter. Our findings support those results. ΔBI generally exhibited a positive partial correlation with nighttime SUHII, except in southern China in winter and autumn. During daytime, positive partial correlations occurred much less frequently and the coefficients were smaller than during nighttime in most cases. No significant partial correlations were found for arid and semiarid western China.

4.1.3. City Size and Total Population

UAS and TP were generally significantly positively correlated with daytime SUHII and foot-print area [29, 30, 32, 42], with the strongest correlations found in summer [29, 30, 42]. Nevertheless, these partial correlations were quite weak [28]. The increases in amplitude slowed as UAS increased, especially when the area was larger than 100 km2. Our work supports these findings. Significant positive correlations were generally found, with the exception of northern China, and the strongest relationships in each region occurred in summer. Nonetheless, in region I, significant partial correlations only existed in summer (r = 0.25, ) and autumn (r = 0.17, ). Nevertheless, different conclusions were also found in previous studies. Insignificant relationships have been found for 32 major cities in China [27]. The positive correlation coefficients associated with winter daytime SUHII in 419 large global cities were found to be quite similar to summer values [28]. These differences may have occurred because the correlation method only considers two variables, while their relationship may be influenced or masked by other factors that play more important roles [28, 54].

UAS and TP were generally significantly positively correlated with summer and winter nighttime SUHII and foot-print area for 32 major cities in China [27, 33], premonsoon summer and winter nighttime SUHII for 84 major cities in India [31], summer nighttime SUHII for 38 major cities in the continental US [30], mean nighttime SUHII during cooling degree days for many global cities [37], and nighttime SUHII and foot-print area for 8 Asian mega cities [32]. Nevertheless, insignificant correlations were found relative to annual SUHII for 39 major cities in mainland China [38] and annual, winter, and summer SUHII for 419 global cities [28]. Moreover, UAS exhibited insignificant, 2% significant positive, and 5% significantly negative correlations with annual, winter, and summer SUHII for 419 global cities, respectively [28]. In our study, significant positive correlation coefficients generally existed, with weak positive partial correlations found in 16% of cases ().

4.1.4. Population Density

The effects of population density were found to be quite small relative to daytime and nighttime SUHII in 419 world cities [28]. In our study, ΔPD and nighttime SUHII were commonly significantly positively Spearman correlated (). In addition, they exhibited extremely significant positive partial correlations in 32% of cases, especially in southern China (). During daytime, ΔPD and SUHII commonly exhibited extremely significant positive Spearman correlations in humid cold region I (). Significant partial correlations were rarely found. The different findings may be due to different study objects or biases in the spatial distribution of the population. Furthermore, the gridded population data used in Peng et al. (2011) had a spatial resolution of 2.5′, whereas the resolution for this study was 1 km.

4.1.5. Vegetation Activity

During daytime, the transpiration of vegetation can increase latent heat fluxes. Therefore, this factor can produce cooling effects on LST and mitigate SUHII [21, 2628, 31, 32, 65], as verified here. ΔEVI exhibited a negative partial correlation with daytime SUHII in all cases except winter and spring in region I, winter in region II, and summer in southern China. This was mainly due to the minimal influence of low-activity vegetation in these regions at these times [18, 27, 28], such as in winter in regions I and II; it may also have been due to the effects of ΔEVI, which was collinearly correlated with other factors in the regression [27]. For example, ΔEVI was correlated with ΔNT in southern China in summer and in region I for the entire year as well as with TP in spring in region I. Note that vegetation has a predominant cooling effect on daytime SUHII in arid and semiarid western China. ΔEVI was the only driver that was partially correlated with SUHII in summer (r = −0.79, ), autumn (r = −0.69, ), and annually (r = −0.66, ).

The rules were more complex during nighttime than daytime. The LST of grass was generally found to be less than that for other vegetation, and short grasses are among the coolest urban elements [48, 49, 67]. Long-wave radiation from grass surfaces is released very quickly after sunset, and those surfaces continue to cool throughout the night due to their low heat capacity, soil canopy, high SVF values, and mixture of bare soil, gravel, and other inert materials [48]. Regarding forests, all-day changes in LST have been observed for various surface features using thermal video radiometers in the summer in Beijing and in the winter in Tel Aviv, Israel. Researchers have found that differences in LST between forests and other surface features are small at night and are less than those in daytime [48, 49]. Nevertheless, the mean forest LST was found to be lower than that of concrete pavement in the summer in Beijing [49] but higher than that of most other urban elements in the winter in Tel Aviv [48]. These differences may have been caused by vegetation types, density, coverage, and surface materials of man-made structures. Moreover, forest LST has been found to be higher than that for grasses or crops throughout the vast majority of China during nighttime [50]. Forests generally experience higher LSTs than most surface features for a number of reasons. First, penetrated solar energy can be stored within and beneath canopies because the water content and high foliage density of vegetation can decrease radiative heat loss from the ground [48]. Therefore, the observed surface temperatures of vegetation reflect a combination of re-emitted long-wave radiation from the vegetation canopy and the underlying ground [48]. Second, vegetation can decrease the sky view factor (SVF), reducing the amount of radiation emitted into the open sky [48]. However, the shading effects of vegetation, especially trees can reduce the amount of heat stored during the daytime and thereby reduce the nighttime SUHII [45, 53, 68]. In our work, ΔEVI generally exhibited a negative partial correlation with SUHII in China, except in arid and semiarid western China. These negative relationships likely exist because the corresponding rural regions of the vast majority of cities were mainly covered by croplands, after limiting the elevation. This conclusion needs to be confirmed using a thermal video radiometer to monitor the LSTs of urban man-made structures and croplands at the same time. The poorest vegetation conditions existed in western China and could not have efficiently produced an effect. Therefore, negative relationships disappeared there. Weak negative contributions were also found for nighttime SUHII in 419 world cities [28]. In addition, negative correlations existed between SUHII and EVI in certain global cities during the growing season [21], and in certain cities in China during both summer and winter [26, 45], between LST and vegetation cover in the entire Phoenix metropolitan area in the USA during summer [69], and in the Yangon region of Myanmar regardless of seasonality [51]. However, insignificant correlations were found between EVI and nighttime SUHII of 32 major cities in China during summer, winter, and the whole year [27, 42].

EVI is commonly used to indicate vegetation activity. Nevertheless, it cannot represent differences in vegetation types. However, its values can correctly reflect the degree of vegetation density when only similar vegetation types are considered. If this is not accounted for, incorrect interpretations and decisions may result. For example, humans should plant trees more densely to mitigate nighttime SUHII, given the results of this study.

4.1.6. Albedo

Materials with lower albedos can absorb more energy. Stored heat can be emitted upward into the atmosphere during the night, contributing to higher surface and air temperatures in the vast majority of studies [27, 28, 38, 62, 70], which is supported by our results. ΔWSA exhibited a negative partial correlation with nighttime SUHII except in southern China in winter and summer. These exceptional cases occurred when the effects of ΔWSA were collinearly correlated with other factors in the regression. For example, ΔWSA was correlated with ΔEVI, P, and T in southern China during the winter (). Both the partial and Spearman correlation coefficients were larger in northern China than in the south. This weaker relationship for southern cities was primarily related to the presence of evergreen forests and continuous croplands surrounding them, which resulted in small positive ΔWSA values [27]. This was also consistent with previous findings that cities surrounded by croplands and deciduous trees exhibited larger contrasts in albedo relative to urban areas adjacent to suburban areas [70].

During daytime, the relationship was complex. Both significant positive and negative relationships have been identified in the literature [28, 52]; however, insignificant cases were also found [27]. WSA was found to weakly explain daytime SUHII for 419 world cities [28], while it was found to be the primary factor affecting seasonal amplitude in Beijing [53]. In our study, ΔWSA was found to exhibit both positive and negative partial correlations with daytime SUHII at the national scale and in the northern regions. The strongest correlation was negative and always occurred in winter. In particular, ΔWSA was the only factor that was partially correlated with daytime SUHII in winter in the humid, cold region I (r = −0.39, ).

4.1.7. Aerosol Optical Depth

Aerosols can decrease LST by reducing the amount of short-wave radiation that reaches the ground surface [38]. They can also increase LST because they are more effective at absorbing and emitting radiation than are water vapour and greenhouse gases in the long-wave atmospheric spectrum under specific conditions [38, 71]. The combined effect depends on initial particle size and growth due to ageing and absorption of water vapour [38, 72]. ΔAOD has been found to positively impact annual mean nighttime SUHII in 39 major cities in mainland China () [38]. The impacts were stronger in semiarid cities than in cities in humid climates because coarser aerosols result in stronger long-wave radiative forcing [38]. Our work is consistent with these findings to some degree. ΔAOD was significantly positively correlated with SUHII in 56.00% of cases. In addition, the correlation coefficients were obviously greater in arid and semiarid western China than they were in other regions of the country. An exception was the negative relationship in winter at the national scale. Moreover, ΔAOD exhibited a positive partial correlation with SUHII in spring (r = 0.11, ) and autumn (r = 0.07, ) in China and in spring in region II (r = 0.15, ). No significant contributions were found for the 39 major cities in mainland China for annual daytime SUHII. ΔAOD was poorly correlated with annual and summer SUHII. In this study, ΔAOD was found to be significantly positively correlated with SUHII in 48.00% of cases. In addition, the largest correlation coefficients (r = 0.60, ) were also found for annual mean SUHII in arid and semiarid western China. Weak and negative partial correlations were only found in summer in China (r = −0.07, ) and autumn in region II (r = −0.14, ), while positive correlations were found in spring in region I (r = 0.34, ) and autumn in southern China (r = 0.12, ).

4.1.8. Climate

Large soil moisture differences between urban and rural regions can occur in areas with high HI or P. These drivers can mitigate nighttime SUHII and enhance daytime SUHII through increases in thermal inertia and soil moisture heat capacity, respectively [26, 27, 38, 42, 54, 73]. This conclusion was supported by our findings. In some cases, HI and P exhibited significant negative and positive partial correlations with nighttime and daytime SUHII, respectively. It is worth noting that farmlands are a highly prominent land cover type in rural regions. In China, 48% of farmland receives water from irrigation rather than from rain [38]. Irrigation causes domesticated plants to behave differently than vegetation in natural ecosystems [38]. This difference can impact the effects of HI and P on SUHII to a large extent [38].

Higher air temperatures can affect LST in both urban and rural regions [55]. This factor can also affect factors such as vegetation activity, evapotranspiration from soils and plants, and anthropogenic heat emissions [27, 42, 55], all of which can directly or indirectly affect SUHII. During nighttime, significant negative correlations have been found between MAT and SUHII in winter, summer, and annually in most studies [27, 28, 42]. Our study verified these findings to some degree. Negative relationships existed at the national scale in all seasons except in autumn. However, significant positive relationships existed in southern China during all seasons except summer, and both significant positive and negative relationships occurred in some regions of northern China during certain timeslots. Moreover, insignificant MAT contributions have been found relative to annual mean SUHII for 39 cities in mainland China [38]. MAT exhibited a negative partial correlation with SUHII during winter and summer during the night in 419 world cities. Nonetheless, significant positive partial correlations existed in some cases, with one exception during summer in southern China.

MAT has been found to be significantly positively correlated with daytime SUHII in winter and annually, but this correlation was insignificant in summer for 32 major cities in China [27, 42]. The results of our study were consistent with these findings. Nevertheless, no significant MAT contributions were found to impact annual mean SUHII for 39 major cities in mainland China [38], and MAT was negatively correlated with SUHII during not only winter but also summer [26]. Moreover, we found significant positive correlations in spring and autumn at the national scale. Both significant positive and negative correlations existed in certain regions in China. In previous studies, MAT has also been found to be significantly negatively correlated with SUHII in summer and annually, while it was found to be quite weakly positively correlated in winter in 419 world cities [28]. SUHII was also found to be approximately 5°C higher during heat waves than on non-heat-wave days in Taipei [55]. The differences in SUHII for both large and medium cities can be as small as 1.0°C, suggesting that SUHII is enhanced in both large and medium-sized cities during heat waves [55]. During an extensive heat-wave event in Cyprus, SUHII was enhanced in two administrative districts but was less affected in two other districts [74]. Moreover, MAT was found to be insignificantly partially correlated with daytime SUHII for 419 world cities [28]. Similarly, no significant partial relationships were found at the national scale in our study. Nevertheless, both positive and negative significant partial relationships were found in certain regions during certain timeslots.

In addition, Yao and Wang [20] found that the LSTs in urban surface were less sensitive to climate variability than in rural areas in northern China. However, this phenomenon was less evident in hot and humid southern China.

4.1.9. Combined Effects

The factors studied here better explained daytime values than nighttime values in winter, summer, and the whole year for 419 global cities in a previous study [28] and 72% of the cases in our work. Nevertheless, opposite finding existed in winter and summer for 32 large cities in China [27]. The summer explanation rates were 4% and 18% smaller than the winter rates for daytime and nighttime SUHII in 32 cities in China, respectively [27]; however, they were 8% and 13% larger for 419 world cities, respectively [28]. Our study found that SUHII was found to be 2.65% greater and 1.68% smaller during summer than winter at the national scale, but this finding showed obvious spatial heterogeneity. During daytime, the explanation rates for summer were higher than for winter in northern China, especially in northeastern China. However, the patterns were reversed in northern China. During nighttime, the explanation rates for summer were slightly smaller than in winter. Nevertheless, different rules existed in the other two regions, especially in western China, where the explanation rates in summer were 30.50% higher than in winter. The largest and smallest explanation rates were found for summer days in humid and cold northeastern China (63.84%) and in southern China (10.44%), respectively. The above results show that different background environmental conditions can have obvious and differing impacts on SUHII. The explanation rates in this study are generally lower than those reported in previous studies, mainly because our study objects are many cities of various sizes, whereas other studies only focused on large or mega cities [27, 28]. In addition, more years of data should be analyzed in order to obtain more reliable conclusions in future [27, 28].

To our knowledge, only one study has developed estimates of SUHII in several cities at large scale simply by using drivers and not adopting a climatologic model [39]. Daytime and nighttime summer and winter SUHII values were estimated for 67 major cities in northern and southern China [39]. In our work, we simulated the SUHII for the vast majority of cities in five ecological function zones and for the entirety of China for each season and annually (Table 3). The simulation accuracy was always higher during the night, which is consistent with previous findings [39]. Nevertheless, the variables analyzed better explained daytime values than nighttime values in 72% of cases. This conflict needs to be explored and discussed further in future. Multisubject knowledge and methods should be combined. To our knowledge, similar studies have not been reported. Nowadays, it is complex to clearly understand the mechanism of SUHIs and accurately estimate the distribution and intensity of SUHIs. The highest estimation accuracies were determined for central-northern China (region II) and were lowest in arid and semiarid western China for both nighttime and daytime SUHII. The RMSE for southern China in this work was 0.05°C and 0.76°C lower than in a previous study of nighttime and daytime summer SUHII in southern China, respectively. However, the present study simulated the nighttime SUHII for 517 cities and the daytime SUHII for 514 cities, rather than both for 27 major cities. Moreover, data on eleven variables related to SUHII were collected in a statistical yearbook in a previous study [39], whereas all indices employed in the present study used common remotely sensed or meteorological data, with the exception of TP and ΔPD, neither of which contributed significantly to daytime SUHII. Weak TP contributions were only found relative to 16% of nighttime SUHII (), whereas ΔPD impacted 32% of cases. In fact, population data were easier to obtain than most other social or economic indices. Therefore, in order to improve the generality of our method, neither can be excluded in practice.

4.2. Uncertainties and Limitations

Although MODIS LST data have been validated and have been found to be highly accurate in many cases and are widely accepted and used [25, 26, 37, 42, 43], uncertainties and issues remain, especially relative to urban environments and rural areas with complex terrain or vegetation caused by the highly complicated spatial heterogeneity of landscape components, greater air pollution in cities, and the known anisotropy issue [37]. The ready-made composite monthly mean LSTs data of China were used rather than high-volume raw daily LST data stored in blocks. The processing of monthly data mainly included mosaicking, slicing, projection transformation, and mean value calculations. Whether quality control information was considered was not determined. Moreover, quality control was not performed on the albedo and aerosol data in order to retain a sufficient number of pixels for this variable. These issues for quality control of LST, EVI, and ALB may have introduced bias or uncertainty to our results [75], especially when cloudy conditions prevail. However, the impact of the MODIS quality on long-term SUHI intensity is weaker than that on the spatial pattern of LST [75]. Moreover, population data are generally collected based on administration cells at the county level. In fact, it is not difficult to accurately identify discrete spatial patterns in the data, especially for internal administrative units within which natural conditions change sharply.

In this study, the spatial extents of cities were derived from the aggregation results of urban land polygons rather than the raw data. Therefore, the raster cells within a certain distance can be treated together as a whole city. The extracted urban land patches are treated as independent cities, if they are far from other urban land patches. After the validation check, the distributions can be well extracted for overwhelming majority cities as the public understanding of humans. Nevertheless, very few exceptional cases existed, such as Panzhihua, whose western and eastern urban areas are more than 4 km apart.

Certain factors were not considered or needed to be modified for SUHII or LST. For example, BI cannot be used to derive detailed information on the two- and three-dimensional characteristics of man-made structures that have strong correlations with LST, including building materials, fraction of buildings or flat impervious surfaces, building volume, density, fraction index, building shadows, mean and maximum height, type, color, floor area ratio, sky view factor, road length, area, and number of nodes [22, 48, 65, 76, 77]. With the exception of EVI, LST is also correlated with the fraction of trees, grass, or vegetation combinations [22, 78], fractional coverage of tree canopy [23], vegetation height [79], canopy configuration (e.g., canopy area and edge) [23], etc. It has also been found to correlate with surface wetness or bareness [22]. Moreover, LST is commonly influenced not only by landscape components but also by landscape configuration and their interactions [22]. The NL data failed to support more detailed social or economic information, such as the proportions of different industries, races or consumption of various resources, and unemployment rates. Moreover, we did not consider certain cropland management activities that may affect SUHII, including irrigation [38] and crop harvest. Finally, more ecological process indices should be introduced using the combination of remotely sensed, observed, and numerical model methods. For example, researchers could attempt to calculate absorbed energy differences by considering AOD, albedo, and incoming solar energy outside the atmosphere in a closed experiment. Obviously, it is unrealistic to consider so many factors at a large scale. A feasible method may be to choose typical cities in various climate zones based on observed spatiotemporal patterns and then study them thoroughly.

5. Conclusions

We studied the associated determinants of SUHII and estimated SUHII values across 1449 cities in China. Some important conclusions are summarized as follows:(1)In arid and semiarid western China, significant positive partial correlations were not found between daytime SUHII and the differences in NL and BI between cities and their corresponding villages. Although UAS and TP were generally positively correlated with daytime and nighttime SUHII, significant partial correlations only existed in 8% and 16% of cases. ΔPD was positively partially correlated with nighttime SUHII in 32% of cases, particularly in southern China. However, only two weak negative partial correlations were found.ΔEVI was generally negatively partially correlated with daytime and nighttime SUHII, except under certain poor vegetation conditions or when disturbed by other collinear factors included in the regression. ΔWSA was usually negatively partially correlated with nighttime SUHII. During daytime, this relationship was complex. ΔAOD was significantly positively correlated with nighttime and daytime SUHII in 56% and 48% of cases, respectively. The correlation coefficients were obviously larger in arid and semiarid western China. Moreover, ΔAOD was weakly positively partially correlated with nighttime and daytime SUHII in three and two cases, respectively. Moreover, weak and negative partial correlations only existed in two cases.In certain cases, HI and P were partially negatively correlated with nighttime SUHII and partially positively correlated with daytime SUHII, respectively. These relationships could be disturbed by common irrigation activities on farmlands. MAT was partially positively correlated with nighttime SUHII in 40% of cases. During daytime, no significant partial correlations existed at the national scale. Nevertheless, both positive and negative relationships were found in certain regions and during certain timeslots.In total, the variables analyzed better explained daytime values than nighttime values in 72% of cases. The best and worst explanation rates occurred during summer days in humid, cold northeastern China (63.84%) and in southern China (10.44%), respectively.(2)Both daytime and nighttime SUHII can be accurately estimated based on drivers using machine learning. The simulation accuracy was always higher during nighttime than during daytime. The highest estimation accuracies for both daytime and nighttime SUHII occurred in central-northern China, while the lowest occurred in arid and semiarid western China. At the national scale, RMSE ranged from 0.49°C to 1.54°C. Significant correlations always existed between the estimated and actual SUHII (). At the regional scale, RMSE ranged from 0.38°C to 2.18°C. Significant correlations existed in 87.50% of cases ().

Data Availability

Both LST and EVI data were downloaded from http://www.gscloud.cn/. Both albedo and aerosol data were provided by https://ladsweb.modaps.eosdis.nasa.gov/search/. The digital elevation data were taken from http://westdc.westgis.ac.cn/data/92bb3089-cc0c-46d2-908d-aae810ef064e. Both land use and population data were provided by http://www.resdc.cn/. The impervious percentage data were downloaded from https://www.beijingcitylab.com/. The nighttime light data were provided by the National Centers for Environmental Information website https://ngdc.noaa.gov/eog/dmsp/download_radcal.html. The temperature and precipitation data were taken from the Meteorological Data Website of China http://data.cma.cn/. The other data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (41701501 and 31600375), Social Science Fund of China (General Projects) (grant no. 17BJL065), Key Scientific and Technological Project of Henan Province (192102310003), Key Scientific Research Projects of Henan Colleges and Universities (18A170014), and Philosophy and Social Science Planning Project of Henan Province (2017BJJ005).

Supplementary Materials

The following supporting information is available as part of the online article: Text S1: introduction of the land use data. Text S2: introduction of the impervious percentage data. Text S3: introduction of the temperature and precipitation data. Text S4: comparison of four intelligence algorithms for prediction of SUHII. Text S5: comparison of three types of indices for prediction of SUHII. Table S1: accuracy statistics derived from the adoption of different indices and intelligence algorithms for daytime and nighttime SUHII simulation in each region of China during the whole year and each season. (Supplementary Materials)