1. Introduction
With increasing socio-economic development globally, the process of urbanization has accelerated, land types such as forests and farmlands have been gradually replaced, and the proportion of impervious surfaces in cities has increased. Coupled with the large amount of pollutants and greenhouse gases emitted from industrial production and inhabitants’ lives, this has resulted in a general increase in urban temperatures, which has led to the even more severe phenomenon of the surface urban heat island (SUHI) [
1,
2]. According to the United Nations Intergovernmental Panel on Climate Change (IPCC), the land surface temperature (LST) continues to rise in almost all regions of the globe and it will continue to increase over the next 70 years [
3]. It has been shown that the SUHI can be observed in 98.9% of Chinese cities at night during the summer months [
4], seriously threatening the ecological environment and human health. With the progress of science and technology, high-resolution remote sensing images provide great help for SUHI research [
5,
6], which helps reveal the characteristics of spatial and temporal changes in the urban thermal environment as well as the mechanism of evolution, promotes the development of measures to mitigate the urban heat island effect, and has a far-reaching impact on the improvement of the ecological environment and enhancement of the well-being of human beings.
There are a variety of methods for quantifying SUHIs [
7,
8,
9,
10,
11], the most commonly used being the urban–rural dichotomy [
12,
13,
14], in which the difference between the average urban and rural LST is defined as the SUHI. According to published research, the urban–rural dichotomy method has the highest frequency of application in studies related to surface urban heat islands, which is approximately 1.5 times higher than other methods, and it is currently widely used worldwide [
15]. The critical issue in the urban–rural dichotomy is the selection of appropriate urban and rural areas. Many scholars have proposed ideas and methods for defining urban and rural areas. Peng et al. [
16] used a city-clustering algorithm (CCA) to define urban areas and calculated the heat island intensity for equal suburbs, smaller suburbs, and larger suburbs (100%, 50%, and 150% of the urban areas). The difference between the determined LST of the urban areas and the average LST within the surrounding n-km buffer zone was taken as the SUHI intensity [
12,
17]. The advantage of this type of method is that it can make the heat island intensity of different periods horizontally comparable; therefore, it is widely used in SUHI research, but its disadvantage is that it has not yet explored a universal buffer width for the time being. The question that needs to be considered is the buffer width that should be chosen for different study areas to better analyze SUHIs. Zhou et al. [
17] formed an urban boundary by aggregating 2 km areas with >50% building intensity. Subsequently, a buffer zone of the same size as the urban area was constructed as a rural area in order to quantify the SUHI. Many researchers have classified ranges corresponding to urban and rural areas based on administrative units [
18,
19]. Fields, forests, and water bodies can be used as representative pixels of rural areas [
13,
20,
21], but they can be subjective. ISA, MCD12Q1, and DMSP/OLS data were used, similar to those in the SUHI study [
22,
23,
24,
25,
26].
Various SUHI characterization methods provide a wide range of options for study; however, there are several challenges [
27,
28]. The rapid development of cities, rapid expansion of built-up land, and conversion of different land use types have made the boundaries between urban and rural areas unclear, and the urban–rural differences between cities are also quite different [
29]. Therefore, for a large-scale study involving multiple cities, each of which varies in size, population density, building intensity, and land cover type, the selection of a single research method is not necessarily applicable to all cities. There have been a few studies conducted to confirm this. Yao et al. [
30] analyzed how different methods and data affected the assessment of SUHI intensity in 31 Chinese cities. This study suggested that ignoring the impact of water bodies and elevation will overestimate the daytime SUHI intensity. In contrast, using nearby suburban areas largely underestimates the SUHI intensity. Liu et al. [
31] assessed SUHIs using seven methods for delineating non-urban areas and found that using surrounding buffers as suburbs resulted in a significant increase in the downward trend of SUHIs, which seriously underestimates the need to mitigate SUHIs in the future. In contrast, non-urban areas, defined according to administrative units, had an overestimated SUHI trend in cities with large altitudinal differences and better green areas. Schwarz et al. [
32] compared the 11 most commonly used SUHI quantitative indicators and found weak correlations. In this regard, he indicated that the variability and instability of quantitative indicators should be considered when conducting SUHI studies and that multiple indicators should be used simultaneously to describe SUHIs.
In summary, the choice of urban and rural areas and SUHI quantification methods significantly affected the SUHI results. However, no uniform standards or systems are available for SUHI studies. Moreover, few studies have compared different heat island quantification methods, and only a few researchers have discussed them [
33,
34,
35]. However, by scientifically and accurately defining the most representative urban and rural areas, analyzing the results of the study using different research methods according to the degree of development and population size of different cities, and screening out the SUHI research methods that were most suitable for the study area, the SUHI results were comparable. Taking the Yangtze River Delta (YRD) urban agglomeration as an example, we have applied GIS spatial statistics, the Mann–Kendall method, and Sen’s slope estimation method based on LST data to explore the influence of different urban and rural delineation methods on the SUHI results, providing references for choosing appropriate SUHI methods to accurately and scientifically quantify SUHIs.
2. Study Area and Data Sources
2.1. Study Area
The YRD urban agglomeration, which was defined in the Yangtze River Delta Urban Agglomeration Development Plan released in 2016, was selected as the study area (
Figure 1). It is comprised of 26 cities: Shanghai, Nanjing, Wuxi, Changzhou, Suzhou, Nantong, Yancheng, Yangzhou, Zhenjiang, and Taizhou in Jiangsu Province; Hangzhou, Ningbo, Jiaxing, Huzhou, Shaoxing, Jinhua, Zhoushan, and Taizhou in Zhejiang Province; and Hefei, Wuhu, Ma’anshan, Tongling, Anqing, Chuzhou, Chizhou, and Xuancheng in Anhui Province. The total build-up area of the YRD increased from 5419.24 km
2 in 2012 to 6112.84 km
2 in 2016. Among them, there was an increase of 445.14 km
2 from 2012 to 2014, and an increase of 248.46 km
2 from 2014 to 2016.
Bordering the Yellow Sea and East China Sea, the YRD is located at the confluence of the river and the sea, with many coastal ports along the river, and an alluvial plain formed before the Yangtze River enters the sea, covering an area of 211,700 km2. In terms of climate, it is mainly a subtropical monsoon climate with sufficient light, four distinct seasons, rainy summers, and annual precipitation of 1000–1400 mm. In terms of topography, the northern part of the study area is dominated by plains, whereas the southern part has a higher terrain with a high majority of mountains and hills. In terms of water resources, the study area is rich in water resources, with numerous rivers and lakes, making it the most densely networked area in China, but also leading to frequent flooding in the region. In terms of economic development, the YRD is one of the most economically active urban agglomerations in China, with a well-developed economy and a high degree of town clustering. In 2016, the GDP amounted to 14.7 trillion RMB and the GDP growth averaged over 8.4%, higher than the national average by 1.7 percentage points. The YRD region has a resident population of over 150 million as of 2016, accounting for 11% of the country’s population, with a difference between the resident population and the household population of approximately 20 million, making it one of the regions with the largest population inflow in the country and strong demographic support.
2.2. Data Sources
The data used in this study included LST, urban boundary, land use, and administrative district data.
The LST product used in this study is the Daily 1 km all-weather land surface temperature dataset for the China’s landmass and its surrounding areas (TRIMS LST) [
36,
37], which was obtained from the National Tibetan Plateau Data Center (
https://data.tpdc.ac.cn, accessed on 29 May 2024), with a temporal resolution of 4 times a day, a spatial resolution of 1 km, and a time span of 2000–2022. The time span of 2003–2022 was selected for this study. Because the Aqua satellite overpasses at 13:30 and 1:30, and the Terra satellite overpasses at 10:30 and 22:30, the Aqua satellites are more representative of day and night temperatures at the surface. Therefore, a dataset prepared using the Aqua satellite was selected for this study.
Urban boundary data were derived from the Global Urban Boundary (GUB) dataset developed by the Peng Cheng Laboratory (
https://data-starcloud.pcl.ac.cn, accessed on 29 May 2024) [
38]. This dataset is based on the developed Global Artificial Impervious Area (GAIA) mapping product with a high resolution (30 m). The time range included 1990, 1995, 2000, 2005, 2010, 2015, and 2018, and in this study, the years 2000–2018 were selected.
Administrative district and land-use data (CNLUCC) were obtained from the Resource and Environmental Science Data Platform (
https://www.resdc.cn/, accessed on 29 May 2024). The years selected for the CNLUCC dataset were 2000, 2005, 2010, 2015, and 2018; the spatial resolution of the dataset was 30 m. In order to satisfactorily calculate the SUHI intensity, this study extracted land types such as farmland, forest, and water body for multiple periods of time.
5. Discussion
The analyses performed in this study revealed that different research methods yielded different results in terms of heat island intensity, heat island areas, and trends over time. In particular, there were significant differences between the urban–buffer and urban–water methods, whereas there were similarities in the SUHI results obtained from the municipal–nonmunicipal, urban–field, and urban–forest methods. In addition, this study quantified and analyzed the SUHI intensities and SUHI areas at three different scales: annual, summer, and winter. This is an advantage of this paper over other previous studies, which focused more on “summer” and “a certain year,” while the content analysis in this paper is more detailed in terms of time and provides a comprehensive description of the YRD SUHI in terms of two indicators, which reveals the similarities between different research methods. This study is more detailed in terms of time and provides a comprehensive description of the YRD SUHI in terms of two indicators, finding that there should be different scopes of application for different research methods, which will be discussed in detail in the following section with the hope of providing scientific references for the study of SUHI.
5.1. Applicability of the Urban–Buffer Method
In the urban–buffer method, the extent of the buffer zone was set based on the periphery of the city, which was equal to the city’s area. In addition to removing the influence of water image elements on the SUHI results, buffer zones beyond the city’s administrative boundaries were excluded to ensure that the selected area was within the study area (
Figure 20). However, for cities with urban cores close to administrative boundaries, such as Shanghai, Hangzhou, Ningbo, Taizhou, and Nantong, the extent of the buffer zone is likely largely excluded. This will result in rural reference areas not being well covered, thus affecting SUHI research results. For cities where urban construction and development are concentrated in the geographic center, the range within the buffer zone is preserved with a higher degree of integrity, such as Hefei, Taizhou, and Xuancheng. At this time, it is more objective, reasonable, and accurate to use the buffer zone as a rural reference area for the study of urban–rural thermal environment differences. In addition, buffer zones are close to core urban areas, may be more susceptible to surrounding urban development, and may not be representative of rural areas undergoing rapid urban development, which can lead to an overestimation of rural background LST. For example, the “Shanghai Metropolitan Area” in the eastern part of the YRD, centered on Shanghai, is the strongest core area in the YRD region, and the buffer zone of the Suzhou–Wuxi-Changsha area, which is adjacent to Shanghai, may underestimate the SUHII of the study area when used as a rural reference. Therefore, this method is not applicable to cities with urban cluster development, and the LST of the buffer zone in this type of urban center is influenced by the development intensity of the surrounding cities, resulting in an increase in the LST of the buffer zone driven by them. In summary, the choice of buffer distance is crucial in the current era of rapid social and economic development, especially for first-tier and second-tier cities, to reasonably grasp the boundaries of “urban,” “suburban,” and “rural,” and to consider the impact of the development of neighboring cities on the study area, which is of great significance to the accuracy of SUHI research.
5.2. Applicability of the Municipal–Nonmunicipal Method
With continuous changes in administrative divisions, there are now cities in the YRD, such as Shanghai and Nanjing, which are fully categorized as municipal districts without nonmunicipal districts. Therefore, the municipal–nonmunicipal method cannot be applied to these cities. In some cities, the municipal district includes a large part of the study area, not just its core area. In reality, the municipal district covers rural areas adjacent to the city center, such as urban villages and non-building land types, such as fields and forests in the subsurface, which have lower surface temperatures, resulting in an underestimation of the average LST value of the municipal district in this study, and thus an underestimation of the SUHII. The results of the division of administrative areas change over time, not only between individual cities but they may also be transformed into each other in different cities. For example, Chaohu City in Anhui Province was abolished in 2011, and the areas belonging to it were assigned to the municipal jurisdictions of Hefei City, Wuhu City, and Ma’anshan City, which led to a certain degree of bias in the subsequent division of the urban and rural areas based on administrative units. Therefore, in the study of SUHIs for a long time series, it is recommended to objectively consider the division of administrative districts of the study area during different periods so that the scope of urban and rural areas can be more in line with the law of urban development and change. Although municipal districts are established in each city based on differences in economic and social development levels, population size structures, etc., there are still differences in the divisions of administrative units in different cities. Therefore, this method is preferred for the study of SUHI in a single city; however, it may have limitations when comparing SUHI differences between different cities.
5.3. Applicability of the Urban–Field Method
The urban–field method considers all of the field image pixels in the study area. In the YRD, the field area is smaller in the south and larger in the north, and the field LST is affected by the planting system, structure, and other factors. This method is applicable to areas where the field area is far from the urban build-up area and where the planting system and area are more stable. This method is suitable for areas where the field is far from urban build-up areas, the image pixels are more dispersed, and the cultivation system and area are more stable. For the YRD, the method is more applicable to the northern region because the northern field covers a wide area, whereas in the southern part, owing to the topography, the flat fields are mostly concentrated around the cities, which may lead to an overestimation of the rural LST.
5.4. Applicability of the Urban–Forest Method
The urban–forest method considers all of the forested land image pixels in the study area. According to the terrain and geomorphology of the YRD, most of the forest is distributed in the south, with less distribution in the north. Taizhou has a very small area of forested land, and as of 2018, there is no forested land. For this type of city, the use of forests as a rural reference range is not representative. It can be seen that the urban–forest and urban–field methods have similar conclusions in most studies. Therefore, if the selected study area contains both large areas of field and large areas of forest land, then either of them can be selected as the rural reference area, or fields can be selected as the rural reference context in plains, and forests can be selected as the rural reference context in mountainous and hilly areas.
5.5. Applicability of the Urban–Water Method
The urban–water method includes all of the water body image pixels in the study area. However, for cities with large lakes, this method may seriously overestimate daytime SUHI results and underestimate nighttime SUHI results because of the specific heat capacity properties of water bodies, which have lower daytime and higher nighttime temperatures. Therefore, when used as a rural reference background, the SUHI has a large diurnal variation and does not accurately reflect actual urban–rural thermal environment differences.
5.6. Shortcomings and Prospects
In addition, this study proposed to strengthen the importance of selecting methods for urban and rural reference areas and to consider various aspects of SUHI research methods applicable to the selected study area in urban heat island related research. However, this study has several shortcomings that should be addressed in subsequent studies. For example, this study compared only five SUHI research methods, which fall into the category of the urban–rural dichotomy and help explore urban–rural thermal environment differences. However, the research methods for SUHIs are not limited to these and other methods are also widely used in China and abroad, such as the local climate zone (LCZ), Gaussian surface model (GSM), urban heat island ratio index (URI), urban thermal field variance index (UTFVI), and some non-parametric models that do not consider urban and rural areas. In the future, we can choose various methods, not only focusing on the urban–rural dichotomy, to study urban heat islands by providing more comparability and selectivity for research methods related to SUHIs and can carry out a comparative study of SUHI among multiple cities in a more reasonable way. Moreover, in this study, an area of equal size on the periphery of the urban core area was designated as the urban buffer zone; however, many current studies consider that the area around the city belongs to the suburbs, whereas the areas farther away from the city are the rural areas, and the buffer zone adjacent to the urban area is still within the scope of the SUHI, resulting in the underestimation of the SUHI [
30]. In contrast, in this study, only equal-sized buffer zones neighboring urban areas were considered, and no comparison was made of the impact on the SUHI results when ranges farther or closer to urban areas were used as rural reference areas. Future studies should consider the extent of the study area more comprehensively and examine how the distance of the buffer zone from the urban center affects SUHIs. Third, in this study, all field, forest, and water bodies in each city were selected, and the cropping system, cropping structure, terrain, slope, distance to the urban area, and stability of land types were not considered; therefore, there is a certain degree of subjectivity in this selection, which can affect the results of the SUHI study. In an actual study, it is recommended that the terrain conditions, geographic location, and cropping structure of the studied areas be included in order to define rural areas more effectively. This is because the type of subsurface planting and crop canopy can affect the monitored LST. Fourth, the urban, buffer zone, field, forest, municipal districts, and nonmunicipal districts selected in this study excluded water body pixels to avoid the influence of water bodies on the surrounding SUHI study but did not take into account the influence of elevation on the SUHI of the study area. The southern part of the YRD has high terrain and complex topographic conditions, and differences in elevation may have affected the study results. In a subsequent study, the effects of the water bodies and elevation on SUHIs should be considered comprehensively in order to determine whether the effects of these factors need to be excluded to minimize the error in the study results.
6. Conclusions
In this study, we have examined the SUHI of the Yangtze River Delta urban agglomeration by selecting different urban and rural reference ranges, exploring the differences in SUHI intensity and SUHI area produced using different urban and rural delineation methods, and examining their impacts on SUHI trend changes. The main conclusions are as follows:
(1) There were significant differences between the results of the different methods in terms of the SUHI intensity. Specifically, the daytime SUHI intensity can be ranked from highest to lowest as follows: urban–water, urban–forest, urban–field, urban–buffer, and municipal–nonmunicipal methods. The nighttime SUHI intensity can be sorted from highest to lowest as follows: urban–forest, urban–field, urban–buffer, municipal–nonmunicipal, and urban–water methods. There were also significant differences in the seasonal and diurnal variations. The SUHI intensity was high in summer and low in winter. The seasonal and diurnal variations in SUHI intensity calculated by the urban–buffer and municipal–nonmunicipal methods were small, whereas the seasonal and diurnal variations calculated by the other methods were large. The SUHI intensity results obtained by the urban–buffer and municipal–nonmunicipal methods were related to the degree of urban development, whereas the SUHI intensity results obtained by the urban–field, urban–forest, and urban–water methods were related to the area of the corresponding type of subsurface and also the degree of agglomeration. The selection of different urban and rural delineation ranges has a significant impact on the SUHI intensity results, and the choice of research method is crucial when conducting heat island studies in different cities.
(2) We found significant differences between the SUHI area values generated using the different methods. In terms of spatial distribution, the results obtained from the urban–buffer method showed that the heat island area was smallest during the daytime, whereas the results obtained from the urban–water method showed that the heat island area was largest during the daytime and smallest during the nighttime. The heat island area was mainly distributed in the southeastern part of the YRD and built-up areas of cities, and the heat island effect was stronger in the southern and eastern parts than in the northern part. With respect to spatial changes, the heat island area was mainly concentrated in the metropolitan area centered on Shanghai and the southern cities of the Zhejiang Province in the early stage. With the passage of time, as well as the policy support of the YRD, the heat island areas gradually expanded outside the cities. The west is represented by the city of Hefei, which had the most significant degree of expansion of the heat island effect. The expansion of the heat island area during the day was mainly in the west and north, whereas the expansion of the heat island area at night was mainly concentrated in the center and south.
(3) By analyzing the annual daytime and nighttime trends as well as the seasonal trends, it was found that the results of the trend changes obtained by the different heat island calculation methods were significantly different. When applying the urban–buffer and municipal–nonmunicipal methods, most cities showed an upward trend, whereas most cities showed a downward trend when applying the other methods. In summer, the selection of different research methods resulted in greater differences in the results, whereas in winter, the results were less different. Changes in the rural reference areas of Nanjing and Tongling produced smaller changes in the study results. Compared to the western cities of the YRD, the eastern cities showed greater differences in the results when using different research methods.
(4) Different methods have different advantages and disadvantages as well as different scopes of application. The urban–buffer method is more suitable for urban core areas that are geographically located in the center of the city, so that the buffer zone can be preserved to a greater extent and is more representative. The municipal–nonmunicipal method is more applicable to studies of a single city, and municipal districts can be selected to consider years of administrative change. The urban–field method is more applicable to cities with a flat topography. The urban–forest method is more applicable to cities with higher elevations. The urban–water method is not applicable to areas with larger water bodies in the study area, which may result in an overestimation of SUHI intensity during the day and underestimation at night.