1. Introduction
Cities are already the main human settlements since global urban population exceeded rural population in 2007, and about 68% of the world population is projected to live in cities by 2050 [
1]. The increasing urbanization trend has brought a variety of challenges to cities, such as urban climate change, limited access to open and public spaces, low-quality housing conditions and constrained access to public transport [
2]. One of the most acknowledged urban climate change phenomena is the urban heat island (UHI) effect, referring to the fact that cities are warmer than their surrounding suburban or rural areas [
3,
4]. Addressing the UHI effect is critical considering its significant impacts on urban systems, citizens’ living and ecosystems [
5]. Many studies have revealed that urban warming leads to an increase in energy and water consumption for cooling [
6,
7], deteriorates outdoor thermal comfort and air quality [
8], and thereby results in an increase in mortality and morbidity [
9]. What is worse, the heat-induced impacts can be more severe because of the interaction of UHI and heat waves that are getting more frequent, longer and more severe along with global warming [
10,
11].
An accurate measurement, assessment and identification of UHI effects is essential for the better communication of urban climate knowledge to urban planners and decision makers [
12]. Understanding surface UHI (SUHI) effects based on land surface temperature (LST) and its associated drivers has been an important theme to achieve this, apart from studies on the canopy UHI effects [
13]. Many studies have adopted SUHI intensity (SUHII), an important indicator that has significant implications for land use and land cover change changes [
14,
15], energy demand [
16] and urban living suitability [
17], to measure SUHI effects, and examined the drivers of SUHII variations. For instance, Peng, Piao [
18] examined the UHII of 419 cities around the world, pointing out that average annual daytime SUHII was higher than the nighttime one and the driving factors for the daytime and nighttime ones were different. Li et al. [
19] analyzed the SUHI effects of Berlin based on several hypothesized scenarios and pointed out that SUHII was affected by city size, urban density and compactness. Liu et al. [
15] further analyzed the SUHI of 1288 urban clusters in China, concluding that the daytime SUHII was also more prominent than nighttime one and the SUHII was a function of urban size, shape, centrality and background conditions.
UHII calculation is a relative value of urban–rural/suburban temperature differences according to its original definition, implying that the UHII can be sensitive to the selection of urban and rural/suburban sites [
20]. On the one hand, the LST in cities is strongly associated with land use/land cover (LULC), where LST of buildings and roads could be up to 10 °C higher than that of water bodies and grass land and the locations of different LULCs in both radial and circumferential directions influence the LST [
21,
22,
23]. Moreover, urban artificial landscapes exhibit a high degree of heterogeneity, making it difficult to accurately choose urban sites for SUHI assessment. On the other hand, the determination of representative rural sites is also important considering the different ecosystems of the rural surface properties [
24]. For instance, Peng et al. [
18] assessed the sensitivity of UHII to the suburban areas in Beijing, suggesting that SUHII could have a similar magnitude when the suburban areas were 50%, 100% and 150% of the urban areas. These results may indicate that the suburban/rural sites should be at least at the suburban ring–buffer zone border of the 50% of the urban areas. Liu et al. [
15] found that the sensitivity of SUHII to the rural sites decreased with the increase of ex-urban distance or area, implying that an arbitrary selection of rural sites could under-estimate SUHII magnitude [
15,
18].
To overcome such challenges in urban climate studies, the identification of the homogeneity and heterogeneity of local morphology has been prioritized as an important approach [
25,
26,
27], among which, the local climate zones (LCZ) scheme was proposed by Stewart and Oke [
20] to standardize surface structure and cover description and thereby standardize urban and suburban/rural sites for temperature comparisons. The LCZ scheme has been widely adopted as an objective framework to analyze intra-urban temperature differences in numerous cities such as Berlin [
28], Dublin [
29], Hong Kong [
30], Nagpur [
31], Nairobi [
32], Olomouc [
33], Shanghai [
34], Vancouver [
35] across various continents [
36,
37,
38,
39]. Overall, the LCZ scheme that is comprised of 10 built types and seven land cover types indicates that both air and surface temperatures reduce from compact to sparse built form and at the same time, from high-rise to low-rise built form [
20]. Nevertheless, the responses of LST to the LCZ types can vary significantly because of the combination of different built-land cover types, geographical conditions, background climates and landscape effect (e.g., topography, distance from the sea) [
31,
37], making it significant to examine the suitability of the application of the LCZ scheme for SUHI studies. The work conducted by Eldesoky, Gil [
40] confirmed that the LCZ scheme is applicable for tropical, temperate and cold climates, but not arid climates, with different levels of uncertainty.
Pending queries relevant to the LST responses to LCZ types include the seasonality or thermal anisotropy [
41]. For instance, Du, Chen [
42] analyzed LST variability of LCZ types in different seasons in Nanjing, China, concluding that the warmest or coolest zones varied with seasons, while the LST of built-up types increased with the reduction of building height. Geletič et al. [
41] analyzed the seasonal LST variations of different LCZ types in Prague, Brno and Novi Sad, reporting that SUHII of a specific LCZ type exhibited the largest difference in summer and spring and the lowest in winter. Meanwhile, the dense built-up type and industrial type had the largest SUHII, and the sparse built-up types had the smallest. Some other studies have also examined the seasonal responses of LST or SUHII to LCZ types [
43,
44]. Nevertheless, the seasonal variation of SUHI to LCZ types has not been well understood, especially in different geographic and macroclimatic contexts. Given the seasonal variability, the applicability of the LCZ scheme to assess SUHI should be examined. Moreover, the impact of landscape effect on the SUHI responses to LCZ types should be further analyzed.
To address the above-mentioned challenges and expand the application of the LCZ scheme in urban planning and design, this study aims to detect monthly variations of LST responses to different LCZ types in the city of Shenyang, China. This study will analyze the impact of landscape effects on such variations through analyzing the variations of SUHII and the LST responses to LCZ types in urban (not fully urbanized) and rural landscapes. Built upon this, the applicability of the LCZ scheme to assess SUHI will be examined. This paper is structured as follows.
Section 2 introduces the case study area and its urban heat challenges.
Section 3 presents data sources and research methods.
Section 4 analyzes the results in aspects of overall LST and SUHII variations, seasonal variations of SUHI responses to LCZ types and the applicability of LCZ schemes for SUHI studies. Following this,
Section 5 discusses the results and
Section 6 concludes this paper.
2. Study Area
This study is conducted in Shenyang (41°11′51″ N–43°02′13″ N, 122°25′09″ E–123°48′24″ E), the capital city of Liaoning Province in the southern part of the Northeast China (
Figure 1). The terrain of Shenyang is flat, gradually extending eastward into hilly areas, with an altitude of 5–441 m. Shenyang has a monsoon-influenced humid continental climate (Dwa), with four distinctive seasons. Summers are hot and humid, where July is the hottest month, with an average temperature of 24.6 °C. Winters are dry and cold, where January is the coldest month with the average temperature of −11.2 °C. The annual average temperature is 8.5 °C. Recently, subject to climate change, Shenyang is also undergoing extreme climatic conditions like many other cities, with the extreme temperature of 39.3 °C on 2 August 2018. South-dominant wind prevails in Spring, Summer and Autumn, while in winter the prevailing wind is north-dominant. The annual average wind speed is about 2.6 m/s, where the wind is the strongest in April and the weakest in August. The rainfall of Shenyang ranges between 600 and 800 mm, while the average annual rainfall of the central city (
Figure 1) is about 716.2 mm. In addition, the annual average relative humidity of Shenyang is about 55.7% [
45].
Shenyang is the geographic center of Northeast Asia. Economically, it is in the center of the Northeast Asia Economic Circle and the Bohai Rim Economic Circle. Being a major city of the Greater Shenyang Metropolitan Area, Shenyang covers an area of 12,948 km
2, where the central city is about 3495 km
2 [
45]. The city has seen rapid population increase in the past 20 years from 7.20 million in 2000 to 9.07 million in 2020 [
46]. Accordingly, Shenyang is now a megacity, the only one in Northeast China. Meanwhile, Shenyang is urbanizing rapidly with the urbanization ratio increasing from 70.33% in 2000 to 84.52% in 2020 [
47]. More than 50% of the population and infrastructure are within the three-ring area (the core city,
Figure 1), consisting of Heping, Shenhe, Dadong, Huanggu, Tiexi and some parts of other districts such as Dongling, Yuhong and Hunnan New District, with the area of about 455 km
2. Within the three-ring area, the SUHII could reach 4–5 °C [
23].
Along with the upward trend of urbanization, Shenyang is undergoing three trends including urban densification (inner city), urban sprawl (outer ring) and industrial structure change (China’s land use status remote sensing monitoring database,
http://www.resdc.cn, accessed on 19 October 2021). First, the population within the three-ring area saw an upward trend from 2.998 million in 1985 to 3.772 million in 2015 [
23]. Second, upon the three-ring area, the urbanized area of central city has been expanding towards different directions, forming a new urban pattern (built-up area A, 635.36 km
2 in area in
Figure 1). Third, Shenyang has traditionally been an old industrial base for heavy industry and manufacture, while such factories are relocating from the middle of Shenyang to the surrounding cities, counties and suburban areas to ensure environmental and living quality during urbanization. Fourth, to create the Greater Shenyang Metropolitan Area, some factories have been relocated to the middle of Shenyang and surrounding cities (e.g., Fushun), along which a new built-up area B (29.03 km
2 in area) forms (
Figure 1).
4. Results and Analysis
4.1. Monthly LST Variation
Figure 3 presents the monthly LST variation in the study area from 2018 to 2020, in which the three-ring area (inner circle) and 100% suburban area (outer circle) were marked to compare the urban–suburban–rural temperature for heat/cool island phenomenon. The results indicate an upward trend in the monthly maximum temperature from 4.35 °C in January to 45.56 °C in August and then a downward trend to 5.53 °C in December. The monthly minimum temperature increased from −32.33 °C in December to 22.52 °C in August and then decreased to −15.91 °C in November. December and January were the two coldest months of the year, with the lowest minimum and maximum temperatures, consistent with the lowest air and surface temperatures in
Table 1.
August was the hottest month with the highest minimum and maximum temperatures, consistent with the data given in
Table 1 as well. However, it should be noted that July has been statistically the hottest month of the year (
Section 2). The disparity might be caused by air temperature differences in this study for 2019 compared to the climatological values, that August was hotter than July (
Table 1). The similar scenario was applicable for the minimum temperatures in April and May.
Moreover, Shenyang is enduring a longer period of cool island phenomenon than heat island. The heat island phenomenon was obvious in July, August and September where urban–suburban temperature was generally higher than rural one. Along the urban–suburban–rural profile, there was an obvious gradient of LST reduction. In comparison, the cool island phenomenon was observed in November, December, January, February, March, April and May when rural temperature was generally higher. Along the urban–suburban–rural profile, there was a LST increase gradient, particularly in January, February, March, April and May. In both October and June, the LST distributed evenly within the study area and a clear pattern of cool/heat island phenomenon was not recognizable. Therefore, October was a month corresponding to the transition of heat island to cold island, while June was the transition month towards heat island from cold island.
4.2. Monthly Variation of Urban Heat/Cool Island Intensity
This section quantifies SUHII and analyses the monthly variation of SUHII of the study area.
Figure 4 presents the monthly LST and SUHII variations of the three-ring, the Central city, in buffer 50%, buffer 100% and buffer 0% scenarios. Overall, the maximum, mean and minimum LSTs in the central city, suburban-50%, suburban-100% and rural scenarios exhibited similar patterns. Nevertheless, some differences were observed, and the extent of such differences was dependent on the month. For instance, mean LSTs of the Central city were lower than those of rural/suburban areas, indicating cool island phenomenon, from November to May. In comparison, mean LSTs of the Central city were higher, indicating heat island phenomenon, from July to September. There were small differences among rural/suburban LSTs and central city LST in June and October. Such results were consistent with the recognition of cool/heat island phenomenon and the transition months in
Section 4.1.
Both maximum and minimum LSTs exhibited different patterns, in aspect of the temperature order of Central city, buffer 50%, buffer 100% and buffer 0% scenarios. Regarding the maximum LST, the Central city had the lowest value from February to October and the rural (buffer 0%) area had the highest value from January to October, compared with those of buffer 50% and buffer 100%. Nevertheless, the central city had the highest minimum LST throughout the year and the buffer 50%, buffer 100% and buffer 0% areas had the same LST. The minimum LST of Central city was about 2.86–13.72 °C higher than those of buffer 50%, buffer 100% and buffer 0% scenarios.
The SUHII varied seasonally and the SUHII was different depending on the definition of rural/suburban areas. The months of June, July, August and September underwent heat island phenomenon, in which the heat island phenomenon in July and August was the most intense with the highest SUHII of 2.86 and 3.06 °C according to the Buffer 0% scenario. The study area underwent cool island phenomenon from October to December, in which the cool island phenomenon was the most intense in February, March and April with an intensity of 2.78, 2.36 and 2.46 °C (Buffer 0% scenario), respectively. It should be noted that whilst June and October experienced heat and cool island phenomena, respectively, their intensities were very weak, further indicating that these two months were transition months. Furthermore, in the buffer 50% scenario, both SUHII and surface urban cool island intensities (SUCII) were the weakest, compared with those in both buffer 100% and buffer 0% scenarios. Such results indicate that a larger rural area corresponded to stronger SUHII or SUCII on the one hand, and it is important to define appropriate rural/suburban areas.
Figure 4 presents the monthly variation of temperature with respect to the three-ring area, while
Figure 5 presents monthly variation of LST of the built-up areas and those of SUHII, in buffer 50%, buffer 100% and buffer 0% scenarios. Overall, the monthly variation of maximum, mean and minimum LSTs in the case of built-up area (
Figure 5) were similar to those in the case of three-ring area (
Figure 5). For instance, the urban area (built-up area) had the highest minimum LST and buffer 50%, buffer 100% and buffer 0% areas had the same LST. From July to September, the mean LSTs of the built-up area was much higher, while the ones of the built-up area were lower from October to June. SUHII based on the built-up area also exhibited a similar pattern to that based on a three-ring area. July and August underwent the highest SUHI effects, and February, March and April had the most intense cool island impacts. However, the cool island phenomenon in May was also obvious, with a SUHII of 2.91 °C (Buffer 0% scenario). June underwent cool island effects according to the built-up area (
Figure 5), different from heat island effects according to the three-ring area (
Figure 4). Under such conditions, June, September and October were the transition months with very weak heat island or cool island phenomenon. In addition, the SUHII in the buffer 0% scenario was the strongest, followed by the buffer 100% scenario and then buffer 50% scenario.
Moreover, the heat island phenomenon according to built-up area was less intense than that according to three-ring area, and the cool island phenomenon according to the built-up area was more intense. Overall, according to the results in
Figure 4 and
Figure 5, the definition of urban area also influences the SUHII, apart from the definition of rural/suburban area. A different definition of urban area could lead to diverse results in SUHII magnitude, cool/heat island phenomenon and transition months.
4.3. Temporal Variation of Land Surface Temperature with Local Climate Zones
Table 3 presents the LST of different types of LCZs in 12 months. The results indicate that in July and August, when the SUHII was the strongest, Heavy industry areas (LCZ-10) had the highest mean LST of 34.58 and 36.87 °C, followed by Compact midrise (LCZ-2, 34.34 and 36.10 °C) and then Lightweight low-rise (LCZ-7). In comparison, the land cover of water was the strongest heat sink with the mean LST of 24.91 and 28.91 °C, respectively, in July and August. Both dense trees and low plants were also heat sinks with the second and third lowest LST.
In February and March, when the cool island phenomenon was the most obvious, the strongest heat sinks and sources were found in the land cover types of the LCZ scheme. The land cover of Low plants (LCZ-D) had the highest average LST of 11.59 and 15.82 °C, followed by Bush, scrub (LCZ-C) with the average LST of 10.26 and 15.50 °C, implying the strongest heat sources. Water (LCZ-G) exhibited the lowest average temperature of −5.04 and 6.90 °C, indicating the strongest heat sink. In addition, that built from Compact high-rise (LCZ-1) exhibited the lowest LST of 6.21 and 10.73 °C, followed by the Open high-rise (LCZ-4) with the average LST of 6.71 and 11.42 °C.
In June, a transition month from cool to heat island, the land cover of Water (LCZ-G) was the strong heat sink (24.83 °C), followed by Dense trees (LCZ-A, 28.28 °C), while Bare rock or paved was a strong heat source (LCZ-E, 32.61 °C), followed by Bush, scrub (LCZ-C, 32.24 °C) and Bare soil or sand (LCZ-F, 32.48 °C). For the built LCZs, the Sparsely built (LCZ-9) exhibited the lowest LST (29.85 °C), while the Heavy industry (LCZ-10) indicated the highest LST (33.59 °C). In October, a transition month from heat to cool island, the land cover of Bush, scrub (LCZ-C) had the highest average LST of 19.05 °C, followed by Bare soil or sand (LCZ-F, 18.86 °C), while dense trees had the lowest average LST of 14.91 °C and then the second lowest LST of Water (LCZ-G, 15.40 °C). The built LCZ of Heavy industry (LCZ-10) had the highest average LST (19.34 °C) while the Compact high-rise (LCZ-1) had the lowest LST of 17.08 °C.
Overall, LST of different types of LCZs varied temporally. In hot seasons (e.g., July, August), Heavy industry (LCZ-10) had the highest LST, forming the strongest heat sources among the built LCZs. Such results were relevant to the waste heat emissions during factory operation in hot seasons [
50] and the strong solar radiation incidence due to large sky view factor [
20]. Following this, the Compact midrise (LCZ-2), Lightweight low-rise (LCZ-7) and Large low-rise (LCZ-8) were also strong heat sources (
Figure 6), which may be because of strong solar radiation incidence. In comparison, in cold seasons (from December to March), Sparsely built (LCZ-9) and Compact low-rise (LCZ-3) had the highest LST among all built LCZs, whilst the Compact high-rise (LCZ-1) had the lowest LST (
Figure 6).
Temporal variation of LST to LCZs was also observed in land-cover LCZs. Whilst the Water had the lowest LST in almost all months, the Scattered trees (LCZ-B), Dense trees (LCZ-A) and Low plants (LCZ-D) exhibited the lowest average LST in September, October and November, respectively (
Figure 6). Both Bush, scrub (LCZ-C) and Low plants (LCZ-D) had the highest average LST from December to May, which may be because of the dual impacts of the acceptance of solar radiation and vegetation. In comparison, the Bare-rock or paved (LCZ-E) and Bare soil or sand (LCZ-F) had the highest average LST from June to September, and the Bush, scrub (LCZ-C) had the highest LST in October and November (
Figure 6).
4.4. Spatial Variation of the Responses of Land Surface Temperature to Local Climate Zones
Responses of LST to different types of LCZs within urban and rural areas is examined to analyze the impact of spatial context. The LSTs of different LCZs in urban and rural contexts are given in
Appendix A and
Appendix B.
Table 4 presents the two maximum LSTs and two minimum LSTs within urban and rural areas in 12 months. Overall, there was an obvious difference between the LST patterns of urban and rural built LCZs. In urban area, LCZ-7 (Lightweight low-rise), LCZ-8 (Large low-rise) and LCZ-10 (Heavy industry) had the highest LSTs among ten built LCZs depending on months, while the LCZ-2 (Compact midrise), LCZ-3 (Compact low-rise), LCZ-8 (Large low-rise), LCZ-9 (Sparsely built) and LCZ-10 (Heavy industry) underwent the highest LSTs within rural context. The LCZ-10 (Heavy industry) within an urban context had the highest temperatures from May to October with the average temperature ranging between 19.27 and 36.88 °C, while the LCZ-10 (Heavy industry) within a rural context showed the highest temperature throughout the year, excluding two months of January and April, with average temperature ranging from −3.43 to 36.83 °C.
Apart from the LCZ-10 (Heavy industry), within the urban area, the LCZ-7 (Lightweight low-rise) and LCZ-8 experienced the top two highest temperature throughout the year excluding July (11 times), and LCZ-8 (Large low-rise) was prominent for five times. In comparison, within a rural context, apart from LCZ-10 (Heavy industry), the top two highest temperatures were observed in different built LCZs, including LCZ-2 (Compact midrise, 4 times), LCZ-3 (Compact low-rise, 2 times), LCZ-8 (Large low-rise, 6 times) and LCZ-9 (Sparsely built, 2 times). For the minimum temperature, LCZ-1 (Compact high-rise), LCZ-4 (Open high-rise) and LCZ-9 (Sparsely built) were the built LCZs undergoing the lowest temperature within an urban context. Likewise, it was such three built LCZs that experienced the lowest LST within a rural context. It should be noted that the built LCZs of Compact high-rise, Compact midrise and Compact low-rise were not the case with the highest temperatures, different from the situation of the whole study area (
Figure 6), while the Compact high-rise underwent the lowest temperature.
Table 5 compares the LSTs of land-cover LCZs within urban and rural contexts. In both urban and rural contexts, LCZ-C (Bush, scrub), LCZ-D (Low plants) and LCZ-E (Bare rock or paved) had the highest temperature among seven land-cover LCZs. Moreover, the LCZs with peaked temperatures with urban and rural contexts were generally the same throughout a year excluding December, May and September. From January to April, LCZ-D (Low plants) had the highest temperature within both urban and rural contexts, LCZ-E (Bare rock or paved) had the highest temperature from July to September, and LCZ-C (Bush, scrub) had the highest temperature from October to November. Within urban area, the LCZ-G (Water) was a strong heat sink, where it had the lowest temperature from December to August, and it had the second lowest temperature from September to November. A similar case was observed in the rural context; LCZ-G (Water) had the lowest temperature from November to August, and it had the second lowest temperature in September and October. Likewise, the LCZ-A (Dense trees) had the minimum temperature within both urban and rural contexts.
Overall, the results indicate the urban context affected the responses of LST to LCZs. The highest temperature within built LCZs within an urban context exhibited a higher divergence compared with that within a rural context. Compared with built LCZs, the temperature of land-cover LCZs was more convergent. Moreover, there were limited differences between the land-cover LCZs with maximum (LCZ-C, -D, and -E) and minimum (LCZ-A, -G) temperatures within urban and rural contexts. Furthermore, the overall LST patterns in urban and rural contexts were different, as presented in
Appendix C and
Appendix D.
4.5. Land Surface Temperature Difference among Different Local Climate Zones
Given the spatiotemporal variations of the responses of LST to LCZs, the suitability of the LCZ scheme for LST differentiation was further examined, as shown in
Figure 7. The results indicate that LSTs of different LCZ types were generally different throughout the year, while some types of LCZs failed to differentiate LST. For instance, there was no significant difference between the LST of LCZ 1 (Compact high-rise) and LCZ-2 (Compact midrise) in December. The same results were found among other pairs including LCZ-2 (Compact midrise) and LCZ-5 (Open midrise), LCZ-4 (Open high-rise) and LCZ-5 (Open midrise), LCZ 3 (Compact low-rise) and LCZ-6 (Open low-rise), LCZ-3 (Compact low-rise) and LCZ-7 (Lightweight low-rise), LCZ-6 (Open low-rise) and LCZ-7 (Lightweight low-rise), and LCZ-8 (Large low-rise) and LCZ-9 (Sparsely built) in December. Such results indicate the suitability of the LCZ scheme was compromised in differentiating urban temperatures. The insignificant difference occurred among the pairs with the same characteristics in terms of compactness (e.g., compact, open) and building height (e.g., low-rise).
The suitability of the LCZ scheme for surface temperature differentiation varied temporally. For instance, in July, there was no significant difference between the LST of LCZ-1 (Compact high-rise) and that of LCZ-5 (Open midrise). The same case was also found in the pairs of LCZ-1 (Compact high-rise) and LCZ-8 (Large low-rise), LCZ-5 (Open midrise) and LCZ-8 (Large low-rise), and LCZ-7 (Lightweight low-rise) and LCZ-8 (Large low-rise). In August, there was no significant difference in the three pairs including LCZ-6 (Open low-rise) and LCZ-4 (Open high-rise), LCZ-5 (Open midrise) and LCZ-8 (Large low-rise), and LCZ-7 (Lightweight low-rise) and LCZ-8 (Large low-rise). In July and August, the insignificant pairs including LCZ-5 and LCZ-8, and LCZ-7 and LCZ-8 were observed. However, two pairs of LCZ-1 and LCZ-5, and LCZ-1 and LCZ-8 were found in July, but not in August.
Suitability of LCZ scheme for LST differentiation was also compromised among land-cover LCZs and it exhibited temporal variations. In December, there was no significant difference between five pairs, including LCZ-A (Dense trees) and LCZ-C (Bush, scrub), LCZ-B (Scattered trees) and LCZ E (Bare rock or paved), LCZ-D (Low plants) and LCZ-F (Bare soil or sand), LCZ-D (Low plants) and LCZ G (Water), and LCZ-F (Bare soil or sand) and LCZ-G (Water). Only temperatures of LCZ-A (Dense trees) and LCZ-G (Water) did not have significant differences in July. The pair of LCZ-C (Bush, scrub) and LCZ-F (Bare soil or sand) did not have significant differences in their surface temperatures.
Table 6 assesses the capability of the LCZ scheme in differentiating LST in different months. It shows that for the built LCZs within the whole study area, the LCZ scheme could differentiate 93.9% of the LST in August, indicating the strongest capability. This was followed by that in January, March and October, with a proportion of 91.1%. In comparison, the LCZ scheme could only differentiate 80.0% of the LST in May, indicating the weakest capability. For land-cover LCZs, the LCZ scheme could differentiate 95.2% of the temperature in July, August and October. However, the capability was the lowest in December with a proportion of 76.2%.
Figure 8 and
Figure 9 reveal the suitability of LCZ scheme in differentiating LST in urban and rural area, respectively. Compared with the whole study area, the suitability weakens within the purely urban area. For instance, in December, there were 13 pairs of LCZ types, indicating insignificant LST difference within the urban area, while there were only six pairs when analyzing the whole study area. Compared with the urban area, there were 23 pairs of LCZs indicating insignificant LST difference in December. Such results indicate the spatial variation of the suitability of LCZ scheme in dividing LST. In hot seasons when urban thermal environments were a critical problem, there were six pairs of LCZs, exhibiting insignificant LST difference in both July and August within the urban area.
In comparison, there were five and six pairs in such two months within the rural area. Nevertheless, the pairs were different. For instance, in July, the six pairs in the urban area were LCZ-1 (Compact high-rise) and LCZ-3 (Compact low-rise), LCZ-1 (Compact high-rise) and LCZ-5 (Open midrise), LCZ-1 (Compact high-rise) and LCZ-7 (Lightweight low-rise), LCZ-3 (Compact low-rise) and LCZ-5 (Open midrise), LCZ-4 (Open high-rise) and LCZ-6 (Open low-rise), LCZ-7 (Lightweight low-rise) and LCZ-8 (Large low-rise). In August, such six pairs in urban area were LCZ-1 (Compact high-rise) and LCZ-4 (Open high-rise), LCZ-1 (Compact high-rise) and LCZ-6 (Open low-rise), LCZ-2 (Compact midrise) and LCZ-7 (Lightweight low-rise), LCZ-2 (Compact midrise) and LCZ-8 (Large low-rise), LCZ-4 (Open high-rise) and LCZ-6 (Open low-rise), and LCZ-7 (Lightweight low-rise) and LCZ-8 (Large low-rise). In addition, the results indicate that many pairs of land-cover LCZs did not have significant differences in their average LSTs.
Table 6 also presents the capability of the LCZ scheme in differentiating urban temperatures of urban and rural areas. It is observed that the built LCZs had the strongest capability in March by differentiating 95.6% of the urban temperatures. Following this, in both February and November, the built LCZs could differentiate 91.1% of the urban temperatures. However, the weakest capability was observed in December where only 71.1% of the urban temperatures was distinct significantly. The capability of built LCZs was far weaker when differentiating rural temperatures, compared with urban temperature. The strongest capability of built LCZs in differentiating rural LSTs was found in July, with 89.9% of the temperature distinguishable. However, only 48.9% of the temperature was differentiated by built LCZs in December, indicating the weakest capability. Moreover, the land-cover LCZs, compared with built LCZs, showed a much weaker capability in differentiating LSTs. The strongest capability of land-cover LSTs was recognized in May, with 85.7% of (urban) LSTs differentiable, while the weakest, with only 47.6% of the (urban) LSTs distinguishable, was observed in December. In July and August, the rural land-cover LCZs exhibited the strongest capability (95.2%), while the lowest (57.1%) was found in September.
5. Discussion and Implications
5.1. Spatiotemporal Variations of Urban Thermal Environments and the Proper Month and Method Selection for Accurate Assessment
This study revealed that in July and August, Shenyang underwent the worst thermal environments with the highest LST, and both December and January were the coldest months throughout the year. Shenyang underwent an obvious heat island phenomenon in summertime, while it had a longer period of cool island phenomenon from November to May. In comparison, both June and October were transition months, implying cool–heat and heat–cool island alterations, respectively. Such results, on the one hand, indicate the dependency of LST and SUHII on month variations because of solar incidence, cooling performance and thermal inertia of vegetation and water bodies, aerodynamic properties and their interactions with urban form (e.g., materials, typology). On the other hand, LST and SUHII have been important indicators for the quantification and assessment of urban thermal environment and heat stress so that an appropriate selection of the month for study is essential, in order to avoid an underestimation of heat-related impacts and hazards. For instance, whilst June and September were also typical months in summer and autumn, the LST and SUHII were much lower than those in July and August.
Apart from temporal variations, the SUHII experienced spatial variations which is relevant to the definition of urban and suburban/rural boundaries. The adoption of economic (three-ring,
Figure 4) and administrative (built-up area,
Figure 5) border to screen urban and suburban/rural areas resulted in the upper limits of SUHII or SUCII. The definition of 50% and 100% buffering zones, adjacent to urban areas, as the suburban/rural areas, resulted in a lower SUHII or SUCII, and a smaller buffering zone corresponded to a weaker SUHI or SUCI phenomenon. This result was different from the findings, reported in a study by Peng et al. [
18] which originally introduced the buffering method for quantifying the SUHI phenomenon, that the definition of 50%, 100% and 150% buffering zones as suburban/rural areas did not make differences to SUHII in Beijing. Our results indicate that such a hypothesis was not applicable to the case study city of Shenyang, where the SUHII based on a 100% buffering zone was 0.13–0.16 °C higher than that based on a 50% buffering zone in July and August. With a larger buffering zone, the SUHII went higher until the upper SUHII limits, about two times that based on the 50% buffering zone.
The definition of the urban area also made a difference to the SUHII or SUCII because of the variation of land use/land cover included in the study area. Our results indicate that the SUHI for the built-up (urban) area was weaker than that calculated for the three-ring (urban) area, while the SUCI for the built-up area was stronger. With the change of urban area, the months that underwent cool island or heat island phenomenon varied as well so that the transition months migrated. Such results further exhibited the significance of definition of urban/suburban/rural area and the importance of determining proper months for investigating urban thermal environment and assessing heat-induced impacts.
Overall, the same definition of urban/suburban/rural area resulting in distinct SUHII within different cities may be relevant to the fact that the SUHII formation is associated with urban form (e.g., urban size, shape, density, centrality) and urban macro climatic background [
15,
19,
51,
52]. Therefore, there is a need to seek for a flexible method for suburban/rural definition that could result in stable SUHII [
15,
24,
53].
5.2. Spatiotemporal Variations of the Responses of Land Surface Temperature to Local Climate Zone
The LCZ scheme was developed to differentiate temperature across different urban zones and existing studies have indicated that Compact high-rise (LCZ-1), Compact midrise (LCZ-2), Large low-rise (LCZ-8) and Heavy industry (LCZ-10) could generally have the highest LST [
20,
54]. Nevertheless, our results revealed that the heat sinks/sources varied significantly with the month. In hot months (e.g., June, July, August), Compact midrise (LCZ-2), Heavy industry (LCZ-10), Open low-rise (LCZ-6) and Lightweight low-rise (LCZ-7) were the built LCZs with the highest LST. Such results were partially similar to the existing results [
20,
54]. However, this was not a fixed pattern. Lightweight low-rise (LCZ-7), Large low-rise (LCZ-8) and Sparsely built (LCZ-9) had the highest LST, while Compact high-rise (LCZ-1) and Open high-rise (LCZ-4) had the lowest LST among ten built LCZs from October to May. Urban greenery and water bodies and their relationship with built LCZs could show different cooling/heating influences due to thermal inertia, resulting in the change of temperature pattern of LCZs [
45,
55]. The inconsistent pattern of LST temperature of different built LCZs also indicate the variations of different combinations of heat sources and sinks, such as deciduous trees for shading and evapotranspiration (a lower urban greenery ratio in winter compared that in summer), the intermittent operation of factories, the use of heating/cooling system in cold/hot seasons and heat linkages, and the alternation of macro climate in hot–cold seasons.
For land-cover LCZs, the Water (LCZ-G) and Dense trees (LCZ-A) were the strong heat sinks throughout the year, while in other studies, the Water (LCZ-G) could be heat sources in winter because of its high inertia [
42]. Such results were consistent with previous studies that water bodies and forestry land were conducive to lowering urban temperature [
23]. Moreover, icing and snow, in combination with Water and Dense trees, fostered its lowest LST in winter. Nevertheless, the contribution of other land-cover LCZs varied temporally. The Bare rock and paved (LCZ-E) was a heat source in warm months (e.g., from June to October), while it was heat sinks from December to May (
Figure 6). Such results may be relevant to the capability of heat absorption and storage of rock and paved materials (e.g., stone, concrete, asphalt) in summer, while in winter the higher radiative capacity, compared with other land-cover LCZS, made it have a lower temperature. Likewise, the Low plants (LCZ-D) made different contributions depending on cold or warm seasons, where it was a heat source from December to May and it was a heat sink from June to November. This result might be relevant to the growth cycle of low vegetation (e.g., rice with irrigation, grass), where in warm months vegetation grows and generates cooling performance, while it changes to the land with dry grasses (e.g., straw, detritus) contributing to heat sources in cold months. Such a result may imply a change of LCZ types in different months.
The response of LST to LCZs diversified with urban and rural contexts. For instance, Lightweight low-rise (LCZ-7), Large low-rise (LCZ-8) and Heavy industry (LCZ -10) were the three built LCZs that had the highest LST throughout the year within an urban context. Different from this, the Compact midrise (LCZ-2), Compact low-rise (LCZ-3), Large low-rise (LCZ-8), Sparsely built (LCZ-9) and Heavy industry (LCZ-10) were the five built LCZs that had the highest LST within rural context. Both the hottest built LCZs within urban and urban contexts were different from those within the whole study area. Nevertheless, the Compact high-rise (LCZ-1), Open high-rise (LCZ-4) and Sparsely built (LCZ-9) were the three built LCZs that had the lowest LST within both urban and rural contexts. Moreover, Bush, scrub (LCZ-C), Low plants (LCZ-D), Bare rock or paved (LCZ-E) and Bare soil or sand (LCZ-F) had the highest LST within urban contexts, which was consistent with that within rural context. Both Dense trees (LCZ-A) and Water (LCZ-G) had the lowest LST within both urban and rural areas. Overall, such results indicate that the responses of LST to LCZs was a function of the scope of area of interest, particularly the highest urban temperature. LCZ scheme have been thought of as an effective tool [
12], to support climate-sensitive urban planning and design (e.g., outdoor thermal comfort, heat exposure, heat stress) [
56,
57,
58]. The spatiotemporal change of the hottest LCZs in our study implies the consideration of only a month or an improper selection of study area for identifying heat stress may lead to inaccurate results and mislead actions of urban heat mitigation and adaptation.
5.3. Suitability of Local Climate Zone Scheme for Urban Temperature Differentiation
Existing studies have found that the suitability of the LCZ scheme for LST differentiation can be affected by macroclimate (e.g., tropical, arid, temperate and cold) [
40] and seasons (e.g., spring, summer, autumn and winter) [
42]. It has been indicated that the LCZ scheme had the weakest capability of differentiating LST within arid climates, and it had a moderate level of capability within cold climates [
40]. The study area of Shenyang is in cold regions and the LCZ scheme could differentiate LST, with 80.0–93.3% and 76.2–95.2% of the urban and rural temperature being differentiated, respectively. Nevertheless, the capability was dependent on different months and urban context (e.g., urban and rural areas). The case study in Nanjing indicated the capability was about 86.1%, 93.1%, 77.88% and 82.0% in spring, summer, autumn and winter, respectively [
42]. In comparison, our study indicated in hot months (from July to October), the LCZ scheme had the strongest capability (89.9-93.3%) and in May it had the weakest capability (80.0%). For land-cover LCZs, our study indicated that the LCZ scheme had the strongest capability from July to October (90.5–95.2%) and the lowest from November to February (76.2–85.7%).
Our study further added a new finding that the capability of the LCZ scheme in differentiating temperature is dependent on the urban context (e.g., urban and rural), apart from macroclimate, seasons and months. Compared with the whole study area (a combination of urban and rural areas), the capability of built LCZs in either urban or rural contexts weakened. From July to October, the built LCZs could differentiate less than 90% of the urban temperature and only 71.1% and 77.8% of the temperature were differentiated in December and January, respectively. Furthermore, in a rural context, apart from 86.7% and 88.9% in August and July, less than 80% of the urban temperature was differentiated in all other months and even peaked at 48.9% in December. Such results can indicate that built LCZs had a higher level of applicability in urban context (highly urbanized area) compared with rural context (barely urbanized area). However, land-cover LCZs indicated a different result from built LCZs that within rural context the land-cover LCZs had a higher level of capability compared with that within urban context. Overall, spatiotemporal suitability of LCZ scheme for differentiating LST implies that the adoption of LCZ in urban planning and design should be pre-examined to avoid misleading results. One-month (e.g., in summer) LCZ data within a specific urban/rural context cannot fully represent and identify heat-induced impacts such as outdoor thermal comfort, heat stress and heat exposure of a city. It is essential to document several-month data based on specific contexts to overall support urban planning and design.
6. Conclusions
An accurate quantification of urban thermal environments is the premise of mitigating and avoiding the several consequences of urban heat challenges which are experienced by many cities. LST and SUHII are two important indicators, and the LCZ scheme is an important tool to differentiate surface temperatures and intra-urban temperatures. However, the landscape effects on monthly variation of SUHII and the applicability of LCZ scheme are not well understood. This study investigated the variation of SUHI effects and the suitability of the LCZ scheme for LST differentiation in Shenyang, China. The findings indicated that both the SUHII and the suitability of the LCZ scheme exhibited spatiotemporal variations. An accurate analysis of SUHII should, therefore, properly define urban and rural contexts and specify the month, and one-month analysis cannot fully represent urban thermal environments of a season. The LST of both built and land-cover LCZs could not follow a fixed order, particularly for the highest temperature that varied significantly with both month and spatial boundary. Moreover, the suitability of LCZ scheme for LST differentiation depended on both urban context and the month. In hot months such as July and August, the LCZ scheme in aspects of both built types and land-cover types had a high level of suitability while in cold months such as November, December and January the suitability weakened. The built LCZs within urban areas could have a higher level of suitability than that within rural areas, while the land-cover LCZs exhibited a reverse pattern. Overall, this paper added new findings on spatial variability and temporal heterogeneity of urban temperature and the applicability of LCZ scheme for LST differentiation. It can also provide important implications for the assessment of heat-induced impacts and supports climate-sensitive planning and design.
This paper has some possible limitations and future studies are needed to reveal the spatiotemporal variations of LST distribution, responses of LST to LCZs and the suitability of the LCZ scheme. First, built upon the landscape in August, this study only applied one suite of the LCZ scheme, which as a result cannot respect the possible LCZ variation (e.g., changes of both built and land-cover LCZs with deciduous trees and vegetation death in cold seasons) throughout the year. Therefore, it is essential to reproduce seasonal and monthly LCZ images to further verify the results and conclusions. Second, the urban and rural contexts considered in this study were differentiated within a same city, which might be not capable of representing highly urbanized and barely urbanized cities. Therefore, future studies are needed to verify the results and conclusions relevant to the urban context with the selection of metropolitans and their surrounding satellite cities. Third, this study analyzed the monthly variation of LST distribution, responses of LST to LCZs and the suitability of LCZ scheme within only one city with Dwa in cold regions and it is wise to conduct further investigations in other cities with diverse climates to gain a better understanding of the dynamics of the LCZ–LST relationships. Fourth, this study analyzed the single-year LST responses to the LCZ scheme, while the multi-year analysis could be of interest to indicate the LCZ–LST relationships in a dynamic context along with urbanization and climate change.