Abstract
The urban heat island (UHI) effect exerts a great influence on the Earth’s environment and human health and has been the subject of considerable attention. Landscape patterns are among the most important factors relevant to surface UHIs (SUHIs); however, the relationship between SUHIs and landscape patterns is poorly understood over large areas. In this study, the surface UHI intensity (SUHII) is defined as the temperature difference between urban and suburban areas, and the landscape patterns are quantified by the urban-suburban differences in several typical landscape metrics (ΔLMs). Temperature and land-cover classification datasets based on satellite observations were applied to analyze the relationship between SUHII and ΔLMs in 332 cities/city agglomerations distributed in different climatic zones of China. The results indicate that SUHII and its correlations with ΔLMs are profoundly influenced by seasonal, diurnal, and climatic factors. The impacts of different land-cover types on SUHIs are different, and the landscape patterns of the built-up and vegetation (including forest, grassland, and cultivated land) classes have the most significant effects on SUHIs. The results of this study will help us to gain a deeper understanding of the relationship between the SUHI effect and landscape patterns.
Similar content being viewed by others
Introduction
Increasingly strong impacts resulting from accelerating urbanization pose a threat to the Earth’s environment1. Among these effects, the urban heat island (UHI) effect, i.e., the phenomenon of higher temperatures in urban areas relative to the surrounding areas2, has been the subject of considerable attention in recent years3,4,5,6,7,8. The UHI effect has potential influences on energy consumption9, vegetation growth10, 11, air and water quality1, 9, and can even cause harm to human health9, 12, 13. Therefore, a better understanding of the UHI effect and its related factors is obviously of critical importance.
UHIs can be broadly divided into atmospheric UHIs calculated from weather station networks, and surface UHIs (SUHIs) evaluated from thermal infrared remote sensing techniques. Compared with the atmospheric UHIs, the SUHIs have the advantages of easier access, wall-to-wall continuous coverage, and direct linkage with surface conditions5, 14. Hence, the SUHI effect has been extensively applied for heat island study3,4,5,6, 15,16,17.
A number of studies have signified that urban landscape patterns, quantified by landscape metrics18 and comprised of two fundamental aspects19, 20 (i.e. landscape composition and configuration), are among the most important factors relevant to SUHIs4, 5, 9, 21, 22. However, to date, almost all the studies that have investigated the relationship between SUHIs and landscape patterns have been conducted only in a single city23,24,25,26,27,28,29,30 or a few cities21, 31, 32. The conclusions drawn from these studies may not be comprehensive or may even be mutually contradictory, because of the limitations of the regional climate conditions, geographical locations, development levels, or some other factors of the target city/cities. For instance, some studies have indicated that a higher density of urban development results in a more obvious SUHI effect21, 26. However, a quite different result was reported in another investigation, where it was concluded that more sprawling city patterns contribute to the formation of surface heat islands29. The conclusion made in some studies that water can significantly alleviate the SUHI phenomenon28, 32 has not been supported in other research26. Moreover, numerous studies have revealed that SUHIs and the landscape patterns of cities vary obviously with seasonal and diurnal changes3, 4, 21, 33. However, most studies that have investigated the relationship between SUHIs and landscape patterns have ignored the seasonal and diurnal factors24, 26,27,28. There is therefore a strong need to understand the correlation between SUHIs and landscape patterns, considering seasonal and diurnal factors, for different cities located in diverse climatic zones.
China has witnessed a rapid urbanization process over the past few decades34, and the trend is expected to continue35. Sustained and rapid urbanization has not only led to significant UHI effects in many cities of China3, 4, 22, but has also profoundly changed China’s urban landscape and land-use patterns36. Moreover, China covers diverse climatic conditions, varying from tropical to subarctic/alpine, and changing from arid to humid37. Meanwhile, located in the East Asian monsoon region, China experiences obviously seasonal dynamics and environmental changes38. These factors make China an ideal study area to investigate the relationship between SUHIs and landscape patterns. Although studies of UHIs and landscape patterns have been carried out in a number of Chinese cities21, 23, 26, 28, few studies have systematically quantified the relationship between SUHIs and landscape patterns at the national scale.
The purpose of this study was thus to investigate the SUHI effect and its relationship to the patterns of the urban landscape in 332 cities/city agglomerations distributed in five different climatic zones of China (Fig. 1). Seasonal (i.e. summer and winter) and diurnal (i.e. day and night) factors were considered in the research. The SUHI intensity (SUHII) was defined as the average land surface temperature (LST) difference between an urban area and its surrounding suburban area (Fig. 2), and the landscape patterns were quantified by the urban-suburban differences in several typical landscape metrics (ΔLMs) (Table 1). The Moderate Resolution Imaging Spectroradiometer (MODIS) LST product (version 5) and China’s Land-Use/Cover Datasets (CLUDs) were used to evaluate the SUHII and ΔLMs, respectively (see Methods). The results of this study will help to deepen our understanding of how landscape patterns affect the UHI effect. In addition, new and important insights could be provided to urban planners and managers on how to mitigate the UHI effect from the perspective of urban landscape design, taking into account the climatic conditions of different cities in China.
Results
SUHIs and landscape patterns across climatic zones
The SUHI intensities (i.e. SUHII) of 332 Chinese cities/city agglomerations were calculated in this study, and the average SUHII values across five climatic zones and China are shown in Table 2. In the winter daytime, the average SUHII of the EW climatic zone (equatorial climate and warm and fully humid temperate climate), situated in the southeast of China, is clearly higher than that of any other climatic zone. Meanwhile, in the winter nighttime, the average SUHII in the EW climatic zone is the lowest, and the most intense average SUHII occurs in climatic zone A (arid steppe and desert climate), which is mainly located in the northwest of China. Interestingly, in the summer daytime, climatic zone S (snow climate with dry winter) and climatic zone TS (tundra climate and snow climate with cool summer and cold winter) experience the highest and lowest SUHII, respectively. Seasonally, the SUHII in summer is mostly higher than that in winter. Diurnally, the daytime SUHII is higher than the nighttime SUHII in summer, but the opposite phenomenon occurs in winter (except for climatic zone EW).
The landscape metrics applied in this study were: (1) percentage of landscape (PLAND); (2) Shannon’s diversity index (SHDI); (3) patch density (PD); (4) mean shape index (MSI); (5) clumpiness index (CI); and (6) contagion index (CONTAG). These landscape metrics were calculated at both the class and landscape levels. Class-level metrics describe the characteristics of each single type of land cover, while landscape-level metrics examine the spatial structure in multi-class patch mosaics (see Methods and Table 1 for more details). At the landscape level, ΔPD (“Δ” means the difference between urban and suburban) and ΔCONTAG are both negative, while ΔSHDI and ΔMSI are both positive in nearly all the climatic zones (Supplementary Table S1). At the class level, the signs of the ΔLMs depend on both the climatic zones and the land-cover types. For instance, the ΔPD values are all no more than zero, while ΔPLAND is negative in almost all the climatic zones for every land-cover type, except for the built-up class (Supplementary Table S1). On average, the percentages of built-up and cultivated land are both above 30% for the urban areas in nearly all the climatic zones (except for the climatic zone situated at the Qinghai–Tibet Plateau (TS)), generally followed by forest, grassland, bare land, wetland, and waterbody (Supplementary Table S2).
The relationship between SUHIs and landscape patterns
Spearman’s rank correlation coefficients of the SUHII and ΔLMs across climatic zones and China were applied to quantify the relationship between SUHIs and landscape patterns. At the class level, we investigated the correlations between SUHIs and landscape patterns for seven different types of land cover (see Methods), respectively.
For the built-up class, SUHII is significantly (p < 0.05) and positively correlated with ΔPLAND for most cases (Table 3). In contrast, ΔCI is mostly negatively correlated with SUHII. It should be noted that, in the winter daytime, the correlation between SUHII and ΔPD is significantly positive in the EW (r = 0.40, p < 0.001) and W (r = 0.58, p < 0.001) climatic zones, but is not significant in the A (r = 0.18, p > 0.05), S (r = 0.08, p > 0.05), and TS (r = −0.02, p > 0.05) climatic zones. Similarly, the winter daytime SUHII has a significant and positive correlation with ΔMSI in the EW (r = 0.24, p < 0.01) and W (r = 0.39, p < 0.001) climatic zones, but it is not significantly correlated with ΔMSI in the other climatic zones.
The correlations between SUHII and the ΔLMs, in most circumstances, are not significant for bare land, waterbody, and wetland (Supplementary Tables S3–S5). At the national scale, for bare land, ΔPLAND is generally positively correlated with the daytime SUHII (r = 0.14, p < 0.05 for annual; r = 0.19, p < 0.05 for summer; r = 0.06, p > 0.05 for winter), and ΔPLAND shows a weak-and-negative correlation with the nighttime SUHII (r = −0.18, p < 0.05 for annual; r = −0.18, p < 0.05 for summer; r = −0.20, p < 0.05 for winter). In contrast, for waterbody and wetland, ΔPLAND is generally negatively correlated with daytime SUHII and positively correlated with nighttime SUHII across China.
Although the forest, grassland, and cultivated land classes all belong to vegetation, the relationships between their landscape patterns and SUHIs appear quite different when seasonal, diurnal, and climatic factors are taken into consideration (see Table 4 and Supplementary Tables S6–S8). For the forest class, the daytime SUHII usually shows negative correlations with ΔPLAND, and most of these negative correlations are significant in summer (Table 4). In comparison, the nighttime SUHII usually shows positive correlations with the ΔPLAND of forest (except for the EW climatic zone), but the majority of these positive correlations are not significant (Table 4). From summer to winter, the daytime alleviation effect of ΔPLAND becomes weaker in the north of China (e.g. for climatic zone A, r = −0.46, p < 0.05 in summer daytime, and r = −0.05, p > 0.05 in winter daytime), and the nighttime enhancement effect of ΔPLAND becomes stronger in many climatic zones (e.g. for climatic zone W, r = 0.01, p > 0.05 in summer nighttime, and r = 0.48, p < 0.001 in winter nighttime) (Table 4). Specifically, the ΔPD of forest is mostly negatively correlated with daytime SUHII and positively correlated with nighttime SUHII (Supplementary Table S6). In addition, there is no significant correlation between SUHII and ΔMSI or ΔCI, in most cases (Supplementary Table S6).
For the grassland class, the increase of ΔPLAND generally has no significant mitigating effect on SUHII in the daytime, and ΔPLAND is even significantly and positively correlated with summer daytime SUHII in the A (r = 0.42, p < 0.05) and S (r = 0.44, p < 0.05) climatic zones (Table 4). Similar to the forest class, SUHII is not significantly correlated with the ΔMSI or ΔCI of grassland (Supplementary Table S7). However, it should be noted that, in the TS climatic zone, SUHII is significantly negatively correlated with ΔPLAND in the nighttime and ΔPD in the daytime, respectively (Table 4 and Supplementary Table S7).
For the cultivated land class, ΔPLAND is significantly correlated with SUHII in nearly all the climatic zones (except for climatic zone TS) (Table 4). In summer, the correlations between SUHII and ΔPLAND are significantly negative in the A and S climatic zones, but are rarely significant in the EW or W climatic zones. In winter daytime, SUHII is significantly and positively correlated with ΔPLAND in the EW and W climatic zones. In addition, the relationship between the ΔCI of cultivated land and SUHII is negative in winter daytime, but positive in winter nighttime (except for climatic zone TS), and most of the correlations between ΔMSI and SUHII are not significant (Supplementary Table S8).
At the landscape level, SUHII is generally positively correlated with ΔCI and negatively correlated with ΔSHDI in all the climatic zones, except for TS (where the opposite correlation occurs). Furthermore, there is no significant relationship between SUHII and ΔPD or ΔMSI, in most cases (Table 5).
Discussion
The necessity of taking seasonal, diurnal, and climatic factors into account
The spatio-temporal heterogeneity of the UHI effect in China is explicitly delineated in this study (Table 2). Firstly, Seasonal differences of SUHII are observed in this study, which might be relevant to the change of the physical and biochemical properties of the land cover3, 4. For instance, through calculating the urban-suburban differences in the Enhanced Vegetation Index (ΔEVI) of the different vegetation types, we can clearly see the seasonal variation of vegetation (Fig. 3). Secondly, we found obvious differences between daytime and nighttime SUHII (Table 2), and these diurnal differences are probably attributed to the fact that the mechanism of daytime heat island formation is different from that at night3, 4, 17, 22. The daytime SUHII was widely considered to be the result of an increase in sensible heat flux and a reduction in latent heat flux due to large areas of vegetated and evaporating soil surfaces are encroached by impervious surface2,3,4, 39, while the release of the more stored energy in the urban zone compared to surrounding area contributes to nighttime heat island4, 17. Thirdly, the intensities of the SUHIs in different climatic zones usually show differences (Table 2). Similar to the findings in a previous investigation4, the southeast of China (EW) tends to experience stronger SUHII in the daytime and weaker SUHII in the nighttime. Many cities located in the northwest of China (A) and the Qinghai–Tibet Plateau (TS) were included in this study, and the average SUHII in the two climatic zones (i.e. A and TS) also shows unique characteristics. Finally, from the perspective of landscape patterns, the values of the landscape metrics in the urban areas and the urban-suburban difference of the landscape metrics (i.e. ΔLMs) both vary with the change of land-cover types and climatic zones (Supplementary Tables S1 and S2). Therefore, it is important to take seasonal, diurnal, and climatic factors into account when investigating the quantitative correlations between SUHIs and landscape patterns.
Theoretical implications of the relationship between SUHIs and landscape patterns
The results of Spearman’s rank correlation analysis show that the relationships between ΔLMs of different land-cover types and SUHII are different, and these relationships are closely related with seasonal, diurnal and climatic factors (Tables 3–5, and Supplementary Tables S3–S8). Generally speaking, three reasons account for this. Firstly, the SUHII itself is related with seasonal, diurnal and climatic factors, and ΔLMs also vary with climatic regions and land-cover types. Accordingly, from a statistical point of view, the relationships between SUHIs and landscape patterns shall be likely connected with seasonal, diurnal and climatic factors. Secondly, different land-cover types cause different feedbacks to LST and further SUHII, due to the variation of their physical and biochemical properties21, 31, 40. For instance, the results of Wang et al. showed the variation and differences of albedo, surface roughness and aerodynamic resistance across the cultivated land, built-up and grassland31. Thirdly, local climatic conditions, including precipitation, temperature, soil moisture and air humidity, have notable influences on the correlation between urban landscapes and heat island effects41, 42. The detailed discussions of the relationships between SUHII and ΔLMs of each type of land cover are given separately below.
The built-up class, as one of the most important land-cover types in urban areas, can be expected to have a great effect on SUHII. Significantly positive correlations between ΔPLAND of the built-up class and SUHII were observed in this study (Table 3), which indicates that the larger percentage difference of the built-up class between urban and suburban areas has the potential to result in more intense SUHII in the cities. This can be directly explained by the physical characteristics of the built-up class. Built-up areas tend to have lower albedo, higher thermal conductivity, and larger heat capacity, which will result in higher surface heat storage3, 28, 43. Besides, with the increase of built-up areas, the heat loss ability in urban areas will be weaker due to denser and taller buildings21. Furthermore, some other SUHI-induced factors, including population, energy consuming and air pollution, are closely relevant to the percent of built-up class43,44,45. These factors will indirectly enhance SUHI effects as the ΔPLAND of the built-up class increases. Several studies, focusing on big cities in China such as Beijing, Shanghai, and Wuhan, have reported that the increase of the patch density, aggregation, and irregularity of the built-up class results in higher LST25, 26, 28. However, the results in this study indicate that the relationships between SUHII and the characteristics of the built-up class vary with the climatic zones. For instance, the correlations between the winter daytime SUHII and the urban-suburban difference in patch density (i.e. ΔPD) are significantly positive in the EW and W climatic zones, but are not significant and even weakly negative across climatic zones A and S (Table 3). The following two reasons may contribute to this. Firstly, the average patch density of the built-up class in the urban areas of cities located in the EW and W climatic zones is higher than that of the cities located in the A and S climatic zones (Supplementary Table S2). Thus, the built-up patches in the urban area of a city located in climatic zones EW or W connect more easily with each other than those in the urban area of a city located in climatic zones A or S when the patch density of the built-up class in the urban area keeps growing. A previous study concluded that the contiguity of urban buildings, regardless of their density, was the critical factor influencing the magnitude of the UHI effect46. Therefore, a positive relationship between SUHII and ΔPD for the built-up class occurs in the EW and W climatic zones, rather than the A and S climatic zones. Secondly, the difference of geographic locations may be another possible reason for the different characteristics of the correlations between SUHII and ΔPD for the built-up class in these two groups of climatic zones. Climatic zones A and S are located in the north of China and have a higher latitude than climatic zones EW and W. As a result, less solar radiation shall be received because of the shorter hours of sunlight and lower solar altitude in winter daytime. Furthermore, in winter daytime, shadows are more easily formed by buildings in climatic zones A and S due to their lower solar altitude21, which reduces the amount of shortwave radiation, and thus decreases the LST in urban areas.
Because of the high thermal capacity and inertia of water21, the wetland and waterbody classes are able to maintain a low temperature in the day, but a relatively high temperature in the night. Therefore, the increase of the urban-suburban difference in the percentage (ΔPLAND) of wetland and waterbody results in a weaker SUHII in the daytime (13:30 pm), but a more intense SUHII in the nighttime (01:30 am) (Supplementary Tables S4 and S5). In contrast, the ΔPLAND of the bare land class is positively correlated with daytime SUHII, but is negatively correlated nighttime SUHII (Supplementary Table S3). The lower water retention rate of bare land is a possible reason for this opposite correlation. However, most of the correlations between SUHII and the ΔLMs of wetland, waterbody, and bare land are not significant (Supplementary Tables S3–S5), which might be due to their very low percentages (1.88 ± 2.66% for wetland, 0.04 ± 0.05% for waterbody, and 3.76 ± 10.9% for bare land in urban areas at the national level, Supplementary Table S2).
Vegetation has a particular functionality for the increase of the latent heat flux and the decrease of the sensible heat flux via transpiration. It can therefore be expected to have a cooling effect on LST, and thus relieve the SUHI effect3,4,5. However, our results indicate that the relationships between landscape patterns of vegetation and SUHIs are very complex, which need in-depth and detailed analyses.
In summer daytime, the increase of the percentage difference of forest between urban and suburban areas (i.e. the ΔPLAND of forest) significantly mitigates the SUHII in many climatic zones except for EW and TS (Table 4). In climatic zone EW, the insignificant mitigation effect of the ΔPLAND of forest is probably due to the aerodynamically smoother in urban areas than surrounding suburban areas31, 41. The cities in climatic zone EW are suited at the southeast of China and have a humid and warm climate in summer. Accordingly, the increase of the ΔPLAND of forest should seem to be able to significantly mitigate the SUHII in EW climatic zone, because the denser vegetation in humid climatic regions has a higher evaporation rate than that in dry climatic regions41, 47. However, at the same time, urban areas will have higher aerodynamically resistance to sensible heat diffusion due to the denser vegetation (especially forest) in surrounding suburban areas41. As a result, the convection of dissipating heat from urban areas is less efficient than from the surrounding suburban areas in cities located in the EW climatic zone, which will offset evaporation cooling effect to a great degree. This may be the reason of the insignificant effect of the ΔPLAND of forest in the EW climatic zone. In comparison, the insignificant correlation between SUHII and the ΔPLAND of forest in the TS climatic zone (Table 4) is possibly due to the particular climatic characteristics of the Qinghai–Tibet Plateau. In winter daytime, SUHII is significantly and negatively correlated with the ΔPLAND of forest in the EW and W climatic zones, but is insignificantly correlated in the other climatic zones (Table 4). This is primarily due to the lower ΔEVI of forest in climatic zones A, S, and TS than in climatic zones EW and W (Fig. 3). However, the ΔPLAND of forest is insignificantly correlated with summer nighttime SUHII in nearly all the climatic zones (except EW), and is positively correlated with winter nighttime SUHII in the W and TS climatic zones (Table 4), which can likely be attributed to the absence of transpiration in the nighttime3, 4.
Grassland, another important type of vegetation, cannot significantly relieve the SUHI effect in most cases, and even enhances the summer daytime SUHII effect in the A and S climatic zones (Table 4). There are two possible reasons for this. Firstly, compared with forest, grassland generally maintains lower rates of evapotranspiration, and thus has a weaker mitigating effect on SUHII47. Secondly, the percentage of grassland is usually less than other types of vegetation (13.57 ± 15.27% for forest, 38.62 ± 14.27% for cultivated land, and only 8.58 ± 13.38% for grassland in urban areas at the national scale, Supplementary Table S2). However, the percentage of grassland in the TS climatic zone is relatively high (36.41 ± 25.12%), and several ΔLMs of Grassland are significantly and negatively correlated with SUHII in the TS climatic zone (Table 4 and Supplementary Table S7).
The effect of cultivated land on SUHII cannot be ignored, because of its relatively high proportion in urban areas across China (except for TS) (Supplementary Table S2). In summer, the increase of the ΔPLAND of cultivated land helps to mitigate SUHII in the A and S climatic zones, but it enhances nighttime SUHII in the EW and W climatic zones (Table 4). The uneven distribution of the two main types of cultivated land, i.e., most of the paddy fields are found in south and central China and the dry farming land is mainly located in the north of China38, might be able to explain this difference. From summer daytime to winter daytime, the mitigating effect of ΔPLAND on SUHII weakens and even turns into an enhancing effect (Table 4), which is probably due to harvesting resulting in a land-cover change from plants to bare land31.
At the landscape level, both the proportion and distribution of the different land-cover types have a pronounced effect on SUHII (Table 5). We found that the increase of the urban-suburban difference in SHDI and CONTAG decrease and increase SUHII, respectively, which is consistent with previous findings26, 28. This suggests that more diverse land-use types and more homogeneous mixing of them in the urban area could help to relieve the SUHI effect, which is important information for urban planning. For instance, in order to reduce the urban surface temperature, we could increase the degree of mutual mixing between buildings and vegetation25.
Significance and uncertainties
In contrast to previous studies on the relationship between the UHI effect and landscape characteristics, which have generally been limited to a single city or a small number of cities, we extended the study area to a national scale, with 332 cities/city agglomerations distributed in various climatic zones. In addition, seasonal and diurnal factors were also taken into consideration, and both class and landscape levels were considered in the study. Our research and findings can therefore be considered to be more comprehensive. For instance, the conclusion that the patch density of the built-up class can significantly enhance UHI intensity has been supported in many studies confined to certain big cities25, 26, 28. However, our study found that it is only in regions with a relatively high building density that the increase of the patch density of the built-up class in the urban area significantly increases SUHII, but not in regions with a low building density. Furthermore, we found that although forest, grassland, and cultivated land classes belong to vegetation, the impacts of their landscape patterns on SUHII show different characteristics, which has been not revealed in other studies3, 4. Overall, the results of this study will help us to gain deeper insights into the relationship between the SUHI effect and landscape patterns.
Some possible uncertainties remain in this research. Firstly, because of the enormous temperature variation of cities in the different regions of China, the direct assessment of the relationship between LST and landscape patterns was not possible in this research. Therefore, in this study, we adopted another way of analyzing the correlation between the urban-suburban difference in surface temperature (i.e. SUHII) and the corresponding difference in landscape patterns (i.e. ΔLMs). Therefore, the difference in the analysis methods needs to be considered when conducting a comparative discussion. Secondly, the numbers of cities in the different climatic zones are different, which is due to the difference in the city densities resulting from the regionally unbalanced development in China. Thirdly, we found that the relationship between SUHII and the ΔLMs in the TS climatic zone, located in the Qinghai–Tibet Plateau, is dissimilar to those in the other climatic zones. The particular geographic location, the special surface conditions, and the lower level of urban development of cities in the Qinghai–Tibet Plateau may be possible reasons for this phenomenon, but the actual causes of this result need to be further investigated.
Methods
Study areas
We divided China into five climatic zones (Fig. 1) characterized by different climatic conditions (i.e. temperature and precipitation) according to the Köppen-Geiger climate classification system48. In total, 336 cities are distributed in the different climatic zones, including 305 prefecture-level cities, 23 autonomous prefectures, and four municipalities (i.e. Beijing, Shanghai, Tianjin, and Chongqing). Several cities (e.g. Jiangmen and Zhongshan in Guangdong province) were aggregated as one city agglomeration due to the urban areas of these cities being spatially contiguous as a city-cluster. Finally, the analysis addressed 332 cities/city agglomerations (Fig. 1).
Datasets
The 2015 China Land-Use/Cover Datasets (CLUDs) with a resolution of 30 m were applied to extract the cities’ urban areas and calculate the landscape metrics. The CLUDs were provided by the Chinese Academy of Sciences, and the overall accuracy of the 25 categories of these datasets has been reported as being as high as 90%49, 50. We reclassified the categories into eight classes: bare land, built-up, cultivated land, grassland, forest, waterbody, wetland, and permanent snow and ice. Due to the inexistence of permanent snow and ice in most of the cities, only the first seven classes were taken into consideration when computing the landscape metrics. LST was obtained from the EOS-Aqua-MODIS 8-day composite product (version 5) with a spatial resolution of 1 km (MYD11A2) during the period from 2014 to 2016. The Aqua MODIS LST data were acquired in both the daytime (13:30 pm) and nighttime (01:30 am), using a split-window algorithm51. The retrieval of LST was further improved by correcting the noise from cloud contamination, topographic differences, and zenith angle changes, obtaining an accuracy of better than 1 K51, 52. Considering the seasonal and diurnal factors, we calculated the daytime and nighttime average LST both annually and for summer (from June to August) and winter (from December to February). The Shuttle Radar Topography Mission (SRTM) 3 arc-second (approximately 90 m) digital elevation model (downloaded from http://earthexplorer.usgs.gov/) was also utilized in order to exclude the altitude effect. The MODIS EVI products are able to reflect the spatial distribution and seasonal variation of vegetation. Seasonally average EVI data (MYD13A3) from 2014 to 2016 were employed to estimate the difference of vegetation between urban and suburban areas.
Landscape metrics
Numerous landscape metrics have been developed to characterize landscape patterns, for both composition and configuration18, 19. For this study, we selected six landscape metrics (Table 1) according to the following principles28, 36, 53: (1) commonly used; (2) minimal redundancy; and (3) interpretable. The landscape composition was characterized by PLAND and SHDI. PLAND is the most frequently applied composition metric25, and SHDI, a popular measurement of a landscape’s diversity, increases as the number of different patch types increases. PD, MSI, CI, and CONTAG were used to describe the landscape configurational features54. PD is equal to the number of patches of the corresponding patch type divided by the entire landscape area, and a high PD usually signifies a more fragmented landscape. MSI is an effective indicator to characterize the complexity of landscapes, and a high MSI is generally due to the more irregular and complex shape of patches. CI and CONTAG are both aggregation metrics, and more aggregated landscapes usually correspond to higher values of these metrics. These landscape metrics were computed with the FRAGSTATS public domain software54. SHDI and CONTAG were calculated only at the landscape level, CI was computed only at the class level, and the other metrics were calculated at both levels. Considering the definition of SUHII, we calculated the average value differences of all the landscape metrics between urban and suburban areas (i.e. ΔLMs), and our main objective was to comprehensively and systematically analyze the relationship between SUHII and the ΔLMs.
Analysis
In this study, we defined SUHII as the LST difference between an urban area and its surrounding suburban area3, 4, 21. The delineation of urban and suburban areas was based on the method proposed by Zhou et al.4, according to the land-cover map (i.e. CLUDs). The suburban area (excluding water pixels) within a ring zone around the urban area covered the same area as the urban area (excluding water pixels) (Fig. 2). In order to reduce bias, suburban pixels that satisfied one of the following two conditions were excluded when calculating SUHII: (1) suburban pixels with elevations more than 50 m higher than the highest point or more than 50 m lower than the lowest point in the urban area5, 6; and (2) suburban pixels in a city falling into the urban area of a neighboring city5. The Spearman’s rank correlation coefficient, a non-parametric method, was applied to quantity the correlation between SUHII and the ΔLMs. The Spearman’s rank correlation coefficient ranges from −1 to 1, and a higher absolute value indicates a more powerful correlation. A positive correlation coefficient means a positive correlation between SUHII and the ΔLMs, and vice versa. The significance test was performed by a two-tailed t-test, and the standard 0.05 significance level was adopted in the analysis.
Data availability
All the land surface temperature data can be downloaded on the MODIS product website (https://lpdaac.usgs.gov/data_access/data_pool). The land-cover classification datasets are provided by the Chinese Academy of Sciences. Other relevant data in this study are available from the authors.
References
Grimm, N. B. et al. Global Change and the Ecology of Cities. Science. 319, 756–760 (2008).
OKE, T. R. The energetic basis of the urban heat island. Quarterly Journal of the Royal Meteorological Society. 108, 1–24 (1982).
Peng, S. et al. Surface urban heat island across 419 global big cities. Environ Sci Technol. 46, 696–703 (2012).
Zhou, D., Zhao, S., Liu, S., Zhang, L. & Zhu, C. Surface urban heat island in China’s 32 major cities: Spatial patterns and drivers. Remote Sensing of Environment. 152, 51–61 (2014).
Imhoff, M. L., Zhang, P., Wolfe, R. E. & Bounoua, L. Remote sensing of the urban heat island effect across biomes in the continental USA. Remote Sensing of Environment. 114, 504–513 (2010).
Zhou, D., Zhao, S., Zhang, L., Sun, G. & Liu, Y. The footprint of urban heat island effect in China. Sci Rep. 5, 11160 (2015).
Feng, H., Zhao, X., Chen, F. & Wu, L. Using land use change trajectories to quantify the effects of urbanization on urban heat island. Advances in Space Research. 53, 463–473 (2014).
Zhu, R., Guilbert, E. & Wong, M. S. Object-oriented tracking of the dynamic behavior of urban heat islands. International Journal of Geographical Information Science. 31, 405–424 (2016).
Rizwan, A. M., Dennis, L. Y. C. & Liu, C. A review on the generation, determination and mitigation of Urban Heat Island. Journal of Environmental Sciences. 20, 120–128 (2008).
Zhao, S., Liu, S. & Zhou, D. Prevalent vegetation growth enhancement in urban environment. Proc Natl Acad Sci USA 113, 6313–6318 (2016).
Zhou, D., Zhao, S., Zhang, L. & Liu, S. Remotely sensed assessment of urbanization effects on vegetation phenology in China’s 32 major cities. Remote Sensing of Environment. 176, 272–281 (2016).
Chen, K. et al. Spatial analysis of the effect of the 2010 heat wave on stroke mortality in Nanjing, China. Sci Rep. 5, 10816 (2015).
Xu, Z. et al. Assessment of the temperature effect on childhood diarrhea using satellite imagery. Sci Rep. 4, 5389 (2014).
Voogt, J. A. & Oke, T. R. Thermal remote sensing of urban climates. Remote Sensing of Environment. 86, 370–384 (2003).
Clinton, N. & Gong, P. MODIS detected surface urban heat islands and sinks: Global locations and controls. Remote Sensing of Environment. 134, 294–304 (2013).
Shen, H., Huang, L., Zhang, L., Wu, P. & Zeng, C. Long-term and fine-scale satellite monitoring of the urban heat island effect by the fusion of multi-temporal and multi-sensor remote sensed data: A 26-year case study of the city of Wuhan in China. Remote Sensing of Environment. 172, 109–125 (2016).
Cao, C. et al. Urban heat islands in China enhanced by haze pollution. Nat Commun. 7, 12509 (2016).
O’Neill, R. V. et al. Indices of landscape pattern. Landscape Ecology. 1, 153–162 (1988).
J., G. E. Quantifying Landscape Spatial Pattern: What Is the State of the Art? Ecosystems. 1, 143–156 (1998).
Turner, M. G. Landscape Ecology: What Is the State of the Science? Annual Review of Ecology, Evolution, and Systematics. 36, 319–344 (2005).
Du, H. et al. Influences of land cover types, meteorological conditions, anthropogenic heat and urban area on surface urban heat island in the Yangtze River Delta Urban Agglomeration. Sci Total Environ. 571, 461–470 (2016).
Tran, H., Uchihama, D., Ochi, S. & Yasuoka, Y. Assessment with satellite data of the urban heat island effects in Asian mega cities. International Journal of Applied Earth Observation and Geoinformation. 8, 34–48 (2006).
Peng, J., Xie, P., Liu, Y. & Ma, J. Urban thermal environment dynamics and associated landscape pattern factors: A case study in the Beijing metropolitan region. Remote Sensing of Environment. 173, 145–155 (2016).
Coseo, P. & Larsen, L. How factors of land use/land cover, building configuration, and adjacent heat sources and sinks explain Urban Heat Islands in Chicago. Landscape and Urban Planning. 125, 117–129 (2014).
Zhou, W., Huang, G. & Cadenasso, M. L. Does spatial configuration matter? Understanding the effects of land cover pattern on land surface temperature in urban landscapes. Landscape and Urban Planning. 102, 54–63 (2011).
Li, J. et al. Impacts of landscape structure on surface urban heat islands: A case study of Shanghai, China. Remote Sensing of Environment. 115, 3249–3263 (2011).
Connors, J. P., Galletti, C. S. & Chow, W. T. L. Landscape configuration and urban heat island effects: assessing the relationship between landscape characteristics and land surface temperature in Phoenix, Arizona. Landscape Ecology. 28, 271–283 (2012).
Wu, H., Ye, L.-P., Shi, W.-Z. & Clarke, K. C. Assessing the effects of land use spatial structure on urban heat islands using HJ-1B remote sensing imagery in Wuhan, China. International Journal of Applied Earth Observation and Geoinformation. 32, 67–78 (2014).
Jr., B. S. & Rodgers, M. O. Urban Form and Thermal Efficiency: How the Design of Cities Influences the Urban Heat Island Effect. Journal of the American Planning Association. 67, 186–198 (2001).
Feng, H., Liu, H. & Wu, L. Monitoring the relationship between the land surface temperature change and urban growth in Beijing, China. IEEE Journal of Selected Topics in Applied Earth Observations & Remote Sensing. 7, 4010–4019 (2014).
Wang, X. et al. Quantifying the contribution of land use change to surface temperature in the lower reaches of Yangtze River. Atmospheric Chemistry and Physics Discussions 1–13 (2016).
Chen, X.-L., Zhao, H.-M., Li, P.-X. & Yin, Z.-Y. Remote sensing image-based analysis of the relationship between urban heat island and land use/cover changes. Remote Sensing of Environment. 104, 133–146 (2006).
Zhou, W., Qian, Y., Li, X., Li, W. & Han, L. Relationships between land cover and the surface urban heat island: seasonal variability and effects of spatial and thematic resolution of land cover data on predicting land surface temperatures. Landscape Ecology. 29, 153–167 (2013).
Gong, P. et al. Urbanisation and health in China. The Lancet. 379, 843–852 (2012).
Seto, K. C., Guneralp, B. & Hutyra, L. R. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proc Natl Acad Sci USA. 109, 16083–16088 (2012).
Fang, C., Li, G. & Wang, S. Changing and Differentiated Urban Landscape in China: Spatiotemporal Patterns and Driving Forces. Environ Sci Technol. 50, 2217–2227 (2016).
Wu, S. Aridity/humidity status of land surface in China during the last three decades. Science in China Series D. 48, 1510 (2005).
Chen, X., Hu, B. & Yu, R. Spatial and temporal variation of phenological growing season and climate change impacts in temperate eastern China. Global Change Biology. 11, 1118–1130 (2005).
Arnfield, A. J. Two decades of urban climate research: a review of turbulence, exchanges of energy and water, and the urban heat island. International Journal of Climatology. 23, 1–26 (2003).
Jin, M. S., Kessomkiat, W. & Pereira, G. Satellite-Observed Urbanization Characters in Shanghai, China: Aerosols, Urban Heat Island Effect, and Land–Atmosphere Interactions. Remote Sensing. 3, 83–99 (2011).
Zhao, L., Lee, X., Smith, R. B. & Oleson, K. Strong contributions of local background climate to urban heat islands. Nature. 511, 216–219 (2014).
Peng, S. S. et al. Afforestation in China cools local land surface temperature. Proc Natl Acad Sci USA 111, 2915–2919 (2014).
Huang, Q. & Lu, Y. Urban heat island research from 1991 to 2015: a bibliometric analysis. Theoretical and Applied Climatology (2017).
Larkin, A., van Donkelaar, A., Geddes, J. A., Martin, R. V. & Hystad, P. Relationships between Changes in Urban Characteristics and Air Quality in East Asia from 2000 to 2010. Environ Sci Technol. 50, 9142–9149 (2016).
Liao, W., Liu, X., Wang, D. & Sheng, Y. The Impact of Energy Consumption on the Surface Urban Heat Island in China’s 32 Major Cities. Remote Sensing. 9, 250 (2017).
Debbage, N. & Shepherd, J. M. The urban heat island effect and city contiguity. Computers, Environment and Urban Systems. 54, 181–194 (2015).
Bonan, G. B. Forests and climate change forcings, feedbacks, and the climate benefits of forests. Science. 320, 1444–1449 (2008).
Rubel, F. & Kottek, M. Observed and projected climate shifts 1901–2100 depicted by world maps of the Köppen-Geiger climate classification. Meteorologische Zeitschrift. 19, 135–141 (2010).
Liu, W., Liu, J., Kuang, W. & Ning, J. Examining the influence of the implementation of Major Function-oriented Zones on built-up area expansion in China. Journal of Geographical Sciences. 27, 643–660 (2017).
Yao, R. et al. Investigation of Urbanization Effects on Land Surface Phenology in Northeast China during 2001–2015. Remote Sensing. 9, 66 (2017).
Wan, Z. New refinements and validation of the MODIS Land-Surface Temperature/Emissivity products. Remote Sensing of Environment. 112, 59–74 (2008).
Rigo, G., Parlow, E. & Oesch, D. Validation of satellite observed thermal emission with in-situ measurements over an urban surface. Remote Sensing of Environment. 104, 201–210 (2006).
Li, H. & Wu, J. Use and misuse of landscape indices. Landscape Ecology. 19, 389–399 (2004).
McGarigal, K., Cushman, S. A. & Ene, E. Spatial pattern analysis program for categorical and continuous maps. Computer software program produced by the authors at the University of Massachusetts, Amherst. FRAGSTATS v4. See h ttp://wwwumassedu/landeco/research/fragstats/fragstatshtml (2012).
Acknowledgements
This work was supported in part by the China National Science Fund for Excellent Young Scholars under Grant 41522110, in part by the National Key Research and Development Program of China under Grant 2016YFB0501403, and in part by the Foundation for the Author of National Excellent Doctoral Dissertation, China, under Grant 201348.
Author information
Authors and Affiliations
Contributions
Q.Y. and X.H. designed the study; Q.Y. and J.L. performed the data analysis; and all the authors contributed to the interpretation of the results and approved the final version of the manuscript.
Corresponding authors
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yang, Q., Huang, X. & Li, J. Assessing the relationship between surface urban heat islands and landscape patterns across climatic zones in China. Sci Rep 7, 9337 (2017). https://doi.org/10.1038/s41598-017-09628-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-09628-w
This article is cited by
-
Impact of COVID-19 Restrictions on the Urban Thermal Environment of Edmonton, Canada
Environmental Management (2023)
-
Variations in Surface Urban Heat Island and Urban Cool Island Intensity: A Review Across Major Climate Zones
Chinese Geographical Science (2023)
-
Spatial-temporal dynamics of alpine grassland coverage and its response to climate warming in Mt. Qomolangma Nature Preserve during 2000–2019
Journal of Mountain Science (2022)
-
The effects of human movements on urban climate over Eastern China
npj Urban Sustainability (2021)
-
Temporal-Spatial Analysis of the Warming Effect of Different Cultivated Land Urbanization of Metropolitan Area in China
Scientific Reports (2020)