Multi-Time Scale Analysis of Urbanization in Urban Thermal Environment in Major Function-Oriented Zones at Landsat-Scale: A Case Study of Hefei City, China

: Urbanization and increasing demand for natural resources and land have affected the urban thermal environment. This is an important hot topic in urban climate research. In this study, we obtained multi-time scale land surface temperatures (LST) at the Landsat scale in Hefei, China, from 2011 to 2020. The evolution of the surface urban heat island (SUHI) was analyzed, and the contribution index (CI), urban thermal ﬁeld variation index (UTFVI), and landscape pattern were evaluated to analyze the thermal environment mechanism of a major function-oriented zone (MFOZ). In addition, we explored the role and mechanism of different MFOZs in a thermal environment. Our results show that the multi-time scale differences in the SUHI were obvious, with the phenomenon of heat islands being concentrated in the main city zone. There are signiﬁcant multi-time scale differences in the CI of different landscapes under the MFOZ. The UTFVI analysis of the MFOZ shows that the livability of the cities in the core optimization zone (COZ) and modern urbanization and industrialization cluster development zone (IDZ) is poor. MFOZ planning moderately alleviated the urban thermal environment of the entire study area, especially in the agricultural development zone (ADZ) and ecological conservation zone (ECZ). This study can guide the planning of the MFOZ and guide decision-makers in selecting governance zones when planning policies or dividing the key restoration areas of the thermal environment.


Introduction
Urbanization in China has accelerated in recent decades. By the end of 2020, the urbanization rate of China's resident population exceeded 60% [1]. Urbanization has reshaped the local geometry and ecological characteristics of urban environments [2]. Among these, land surface change is one of the most direct results of urbanization. One of the direct changes in land conditions due to urbanization is radiative forcing with the alteration of land cover and surface albedo [3]. Changes in land use and cover formation, along with urban canopy structure, usually lead to an increase in urban heat storage and heat within the atmosphere [4], which shapes the local climate and results in the surface urban heat island effect (SUHI).
The SUHI effect was discovered in the early nineteenth century [5]. Subsequently, many researchers have studied this by observing the air temperatures of urban and suburban areas in cities of different latitudes and types [6]. With the development of remote sensing technology, SUHI effects are typically estimated from thermal infrared remote sensing techniques [7,8]. In general, current research mainly studies the impact of urbanization on SUHI from the aspects of surface biophysical parameters, socio-economic factors and landscape components, diversity, and configuration factors [9]. SUHI analysis helps rethink urban planning and execution based on observable environmental indices [10]. This could enlighten us to leverage these legacies and contribute to sustainable and harmonious future planning [11]. Owing to the high heterogeneity of urban areas, it is difficult to capture the representative values of urban temperature using standard weather stations [12]. In addition, dense meteorological networks, which are rarely available in cities, are necessary to capture intra-urban temperature variability, which affects heat exposure [13,14]. Land surface temperatures derived from satellite thermal infrared sensors have the advantages of high efficiency, low cost, and broad spatial coverage that traditional field measurement methods do not provide [15,16]. Therefore, this method has been used to explore urban thermal environment variations on regional and global scales. Previous studies have clarified the typical research methods of urban thermal environment research as follows [17]: (1) the direct observation of the temporal and spatial evolution of land surface temperature (LST) to study the problems of surface urban heat islands (SUHI) [18,19]; (2) the analysis of the direct spatial relationship between various factors from the perspective of vertical dimensions and models [20][21][22], mainly including LST, the relationship between specific land use types [23], normalized difference vegetation index (NDVI) [24], normalized difference built-up index (NDBI) [25], socio-economic factors [26], and meteorological factors [27]; and (3) research on the thermal environment and the relationship between urban land system architecture, such as the shape, size, and pattern of local climate zones (LCZs), and connectivity [28][29][30][31]. Recently, attention has been paid to urban forms, such as the 2D/3D urban morphology [32] and spatial patterns [33] of SUHIs, concerning possible mitigation strategies.
In summary, the existing research has fully explored the urban thermal environment, and research methods are increasingly varied, providing many useful ideas for an in-depth understanding of the evolution of the thermal environment and an exploration of the response factors [34,35]. Many methods for the quantitative assessment of urban thermal environments across varying spatial or temporal scales have been rapidly developed with the availability of thermal infrared remote sensing data from satellites. However, the trade-off between the spatial and temporal resolutions of TIR sensors limits the study of urban thermal environments in both spatial and temporal domains [36,37]. Moreover, owing to the failure to fully consider the spatial functional differences in the study areas within the scope of implementation, the prevention and control effects of policies are often limited [38,39]. Major function-oriented zone (MFOZ) planning with spatial governance as the basic concept can fully consider the spatial differences in the thermal environment, emphasize the spatial and overall nature of prevention and control, and provide a way to accurately control the thermal environment effect through differentiated policies. A better understanding and monitoring of the thermal environment effect of urban functional zones is critically important in improving the quality of life and the environment of urban residential areas and could also help achieve sustainable development goals, especially "sustainable cities and communities" [39,40]. To fulfill these research goals, we used the spatiotemporal reconstruction and fusion method used in our previous study [41] to obtain a multi-time scale (interannual, seasonal, and diurnal) LST at the Landsat scale. Using the time series, Landsat-scale LST derived from MODIS and Landsat to simulate the urban thermal environment of the MFOZ, diurnal and seasonal changes and spatial details can be captured by partially compensating for different data availability. Then, the contribution index (CI), the urban thermal field variation index (UTFVI), and landscape patterns were used to analyze the multi-time scale evaluation of the long-term thermal environmental effects on the MFOZ scale during urbanization.
The city of Hefei, China, was chosen as a representative case study because of its thermal environmental changes due to urbanization. We believe Hefei could become a representative case, considering that it is a comprehensive national science center in conjunction with Shanghai and Beijing, a strategic dual-node city of "One Belt One Road" [42] and in the "Yangtze River Economic Belt" [43]. By 2020, the resident population of Hefei exceeded 9 million, the city's GDP was 1004.572 billion CNY, and the urbanization rate of the resident population had reached 76%, making it a regional megacity with a strong influence on the Yangtze River Economic Belt [43]. In this assessment, we aimed to (1) reveal the diurnal and seasonal spatiotemporal evolution characteristics of the SUHI effect in Hefei, (2) analyze the evolution characteristics of spatial details of CI and estimation of UTFVI under different MFOZs to understand the regulatory role of the urban thermal environment in the process of urbanization, (3) analyze the influence of landscape patterns under different MFOZs on the urban thermal environment and explore effective ways to alleviate urban heat islands in the process of rapid urbanization, and (4) explore the impact of MFOZ planning on the urban thermal environment and provide theoretical guidance for future MFOZ planning schemes. The results from this study provide crucial feedback to urban thermal environment investigations, especially long-term trends, and help develop urban thermal environment strategies under different MFOZs.

Study Area
The study was conducted in Hefei, the capital city of Anhui Province, with a geographic location of 116 • 41 -117 • 58 E and 30 • 57 -32 • 37 N ( Figure 1). It has a subtropical humid monsoon climate, with distinct monsoons and four distinct seasons. The average annual temperature is 15.7 • C. The temperature is higher in June, July, and August and lower in December, January, and February. The average temperature in July is the highest, about 28.1 • C, and the temperature in January is the lowest, at approximately 2.6 • C. The average annual precipitation is approximately 1000 mm. The month with the highest precipitation is July, with an average precipitation of approximately 174 mm. The month with the least precipitation is January, with an average precipitation of approximately 25 mm. Precipitation during spring, summer, autumn, and winter accounts for 25%, 44%, 20%, and 11% of the total annual precipitation, respectively. The major types of vegetation are evergreen and deciduous forests, shrubs, and short grass. The main cultivated land types are rice, wheat fields, rape, and cotton. Most of the built areas do not use space heating for winter but cool it for summer. Hefei comprises four counties (Feidong, Feixi, Changfeng, and Lujiang), one county-level city (Chaohu), and four districts (Yaohai, Luyang, Shushan, and Baohe). According to the statistical data of climate resources released by the Hefei Meteorological Bureau, the seasons are defined as spring (March, April, and May), summer (June, July, and August), autumn (September, October, and November), and winter (December, January, and February).
According to the differences between resource and environment carrying capacity, resource endowment, development foundation and development potential, and considering the strategic pattern of urban space, agricultural space, and ecological space into an overall consideration, the MFOZ planning of Hefei divides the land space into five MFOZs: the core optimization zone (COZ), new-type urbanization and industrialization cluster development zone (IDZ), agricultural development zone (ADZ), ecological conservation zone (ECZ), and prohibited development zone. The prohibited development zone consists of 123 broken patches, which are difficult to count. Therefore, the map of Hefei's MFOZ planning in the official documents includes only four MFOZs, as shown in Figure 1. A map of the prohibited development zones is provided separately. This study did not explore prohibited development zones.

Materials and Methods
The methods used in this study comprised several elements, including Landsat and MODIS LST processing, UHI intensity, the land-use/land-cover (LULC) interpretation, contribution index, and the urban thermal field variation index. The general approach is summarized in the flowchart shown in Figure 2.

Materials and Methods
The methods used in this study comprised several elements, including Landsat and MODIS LST processing, UHI intensity, the land-use/land-cover (LULC) interpretation, contribution index, and the urban thermal field variation index. The general approach is summarized in the flowchart shown in Figure 2.

Satellite Data Acquisition and Pre-Processing
Cloud-free Landsat images (Landsat 5/7/8) covering Hefei were obtained from 2011 to 2020 (the acquisition time of Landsat was approximately 10:30 a.m. local time) from the

Satellite Data Acquisition and Pre-Processing
Cloud-free Landsat images (Landsat 5/7/8) covering Hefei were obtained from 2011 to 2020 (the acquisition time of Landsat was approximately 10:30 a.m. local time) from the official website, USGS GLOVIS (http://glovis.usgs.gov/, accessed in 12 May 2021). The spatial resolutions for the different Landsat thermal bands were 120 m (Landsat-5 thematic mapper TM), 60 m (Landsat-7 enhanced thematic mapper plus ETM+), and 100 m (Landsat-8 thermal infrared sensor TIRS). The radiative transfer equation (RTE) method [44,45] was applied for LST retrieval using these images. A series of diurnal MODIS LST products (MODIS/Terra LST/Emissivity Daily L3 Global 1 km SIN Grid V006; MOD11A1) from 2011 to 2020 were obtained from the USGS (acquisition times of MOD11A1 were approximately 10:30 am and 10:15 pm local time). The Landsat and MODIS data were projected onto the Universal Transverse Mercator (zone 50N) projection system and resampled to a spatial resolution of 30 m. Diurnal and seasonal Landsat-scale LSTs were obtained using the spatiotemporal reconstruction and fusion method used in our previous study [41]. The data are presented in Table 1. Table 1. Data used in the study (The dates of the fusion data are marked in red).

Data Dates
Landsat

Intensity of Surface Urban Heat Island
The average standard deviation method uses a combination of the average surface temperature and different standard deviation multiples to segment the surface temperature field to effectively describe the urban heat island [46]. In the evolution of SUHI in different phases, it can avoid the time phase difference to a certain extent [47]. This method indicates that areas with above-average temperatures are more likely to be classified as urban heat islands [47]. We divided the LSTs into five categories using the mean standard deviation method: high-temperature, sub-high-temperature, medium-temperature, sub-medium- temperature, and low-temperature zones ( Table 2). The high-temperature and sub-hightemperature zones are urban heat island zones. Table 2. The classification of the surface temperature and the used formula (T s is the zone of different temperature levels, µ is the average temperature of the study area, and std is the standard deviation of the LST in the study area).

Levels Formula
High The intensity of surface urban heat islands (SUHII) [48,49] is calculated as follows: where T H is the average temperature of the heat island zone, and T O is the average temperature corresponding to the other temperature levels, except for the heat island zone.

The Contribution Index (CI) of MFOZ
The COZ is the most densely populated area in Hefei. The IDZ has obvious location and transportation advantages and a good industrial foundation, as it is the main manufacturing agglomeration zone in the city. The ADZ has a good foundation for agricultural development and is an important grain, cotton, and oil-producing region. The ECZ is an important ecological barrier and ecological product supply area in the city. The ecological environment of the ECZ is fragile and sensitive. In this study, a contribution index [38,39,50] was used to quantify the contribution of each MFOZ or landscape to the regional thermal environment. The contribution index is calculated as follows: where CI is the contribution of the MFOZ or landscape to the regional LST, T a is the average surface temperature of zone a or landscape a, T Avg is the average surface temperature of the region, S a is the area of zone a, and S is the area of the region.

The Urban Thermal Field Variation Index
The urban thermal field variation index (UTFVI) was used to evaluate the effect of the SUHI considering the ecological aspects of the municipality, according to the methodology proposed by Zhang [51]. The retrieved surface temperature was used to calculate the UTFVI [52]. The UTFVI can be calculated using Equation (3) as follows: where UTFVI is the urban thermal field variation index, T s is the temperature of a pixel in the study area, and T m is the average temperature over the entire study area. The UTFVI was used to quantitatively describe the urban heat island effect [51]. In this study, we used the UTFVI as an ecological indicator to evaluate the urban environmental quality [53]. To intuitively describe the change in the urban thermal field, the threshold method was used to divide the urban thermal field into six levels [53]. The threshold divisions of the ecological indicators and their corresponding ecological significance levels are presented in Table 3.   Figure 3 shows the day and night SUHI images of Hefei during the different seasons from 2011 to 2020. The fused Landsat images presented more detailed information (e.g., image texture structure) and clearer boundaries than the MODIS images. The urban heat island phenomenon occurs in spring and summer during the day and becomes more obvious in summer, whereas the SUHI in other seasons is not obvious. Apparently, SUHI has been concentrated in the main city zone of Hefei over the past 10 years. The built-up areas of Changfeng, Chaohu, and Lujiang were small, and the effect of SUHI was not prominent. Heat islands appear on cultivated land during the spring, autumn, and winter. Because part of the cultivated land is not cultivated during these seasons, the surface of the land is exposed and rapidly heated by solar radiation. Thus, the surface temperature was higher. An urban heat island exists in all seasons at night, which is relatively evident in spring and summer. In addition, a heat island also appears in water bodies at night, mainly because the larger specific heat capacity of the water area slows down the cooling rate at night and, hence, increases the temperature. The temperature at night in areas covered by vegetation during winter was relatively high. Owing to the high thermal inertia of vegetation, dense vegetation has a heat preservation effect at night, and the temperature is higher than in other areas [54].
Overall, the heat island area in the main city zone of Hefei has increased with the expansion of the main city zone from 2011 to 2020 and is highly consistent with the spatial distribution and direction of urban expansion. It has mainly expanded to the southwest and northeast of the main city zone and is comprised of high and sub-high zones. The urban heat island phenomena in Changfeng, Chaohu, and Lujiang did not change significantly. In the last 10 years, the main city zone of Hefei has expanded significantly, the underlying surface has changed significantly, and the impervious surface area has increased, prominently changing the thermal environment of the study area.
Overall, during the day, the area ratio of the medium zone was higher than that of the sub-high zone, whereas the area ratios of the high, sub-low, and low zones were similar from 2011 to 2020 (Figure 4). During the day, the area ratio of the medium zone was higher than that of the sub-low zone, whereas the area ratios of the high, sub-high, and low zones were relatively similar. Seasonal differences in the area ratio of the temperature level zones were evident, and the area ratio of the medium zone was the highest in all seasons. In general, the SUHII showed significant diurnal and seasonal differences ( Figure 5). The SUHII in summer was higher during the day than at night. In spring, SUHII during the day was similar to that at night. Overall, the intensity and fluctuation range of SUHI were higher in spring and summer than in autumn and winter. This is because the vegetation is lush in spring and summer, which are growing seasons. Compared with impervious surfaces, such as concrete buildings, areas with high vegetation coverage usually have lower sensible heat flux and higher latent heat flux owing to the strong effect of evaporation, which causes the temperature to be lower than that of impervious surfaces, such as buildings [55,56]. The seasonal difference in SUHII was higher during the day than at night, and the fluctuation range in different seasons at night was more concentrated than that during the day (especially when excluding water bodies). In addition, the urban heat island at night in autumn and winter was larger than that in the daytime, but this phenomenon weakened after excluding water bodies. This is mainly due to the more effective radiative cooling in vegetated land compared to urban canopy structures with complex geometry and large heat capacity. The larger specific heat capacity of water slows the cooling rate at night. Overall, during the day, the area ratio of the medium zone was higher than that of the sub-high zone, whereas the area ratios of the high, sub-low, and low zones were similar from 2011 to 2020 (Figure 4). During the day, the area ratio of the medium zone was higher than that of the sub-low zone, whereas the area ratios of the high, sub-high, and low zones were relatively similar. Seasonal differences in the area ratio of the temperature level zones were evident, and the area ratio of the medium zone was the highest in all seasons. In general, the SUHII showed significant diurnal and seasonal differences ( Figure  5). The SUHII in summer was higher during the day than at night. In spring, SUHII during the day was similar to that at night. Overall, the intensity and fluctuation range of SUHI were higher in spring and summer than in autumn and winter. This is because the vegetation is lush in spring and summer, which are growing seasons. Compared with imper-  than at night, and the fluctuation range in different seasons at night was more concentrated than that during the day (especially when excluding water bodies). In addition, the urban heat island at night in autumn and winter was larger than that in the daytime, but this phenomenon weakened after excluding water bodies. This is mainly due to the more effective radiative cooling in vegetated land compared to urban canopy structures with complex geometry and large heat capacity. The larger specific heat capacity of water slows the cooling rate at night.

Diurnal and Seasonal Contribution of MFOZ
In this study, LULC data from 2011, 2014, 2017, and 2020 were selected as examples to further calculate the percentage of each landscape area under the MFOZ. The interpretation of remote sensing images refers to land use data of 30 m released by the Data Center of the Institute of Geographic Sciences and Natural Resources Research of the Chinese Academy of Sciences (http://www.resdc.cn/, accessed in 28 May 2021). Referring to the National Land Use Cover Classification System of remote sensing monitoring, the research results of land use in Hefei, and the current status of land use in the study area, the four types of LULC extracted in the study were defined as cultivated land, construction land, water, and vegetation [57,58]. According to the accuracy evaluation, the kappa coef- Referring to the National Land Use Cover Classification System of remote sensing monitoring, the research results of land use in Hefei, and the current status of land use in the study area, the four types of LULC extracted in the study were defined as cultivated land, construction land, water, and vegetation [57,58]. According to the accuracy evaluation, the kappa coefficients of the remote sensing images for 2011, 2014, 2017, and 2020 were 0.89, 0.88, 0.87, and 0.89, respectively.
For the entire study area, with the increasing urbanization process, the cultivated land occupied by construction land decreased from 2011 to 2020. This may have led to the occupation of agricultural water. Water and vegetation were relatively stable, with little change. The four MFOZs showed significant differences in terms of area and landscape composition. The ratios of the four MFOZs (COZ, IDZ, ADZ, and ECZ) in Hefei were 1%, 46%, 28%, and 25%, respectively ( Figure 6). The level of urbanization rapidly decreased from the COZ to the ECZ. Although only 0.01% of the study area is covered, according to Hefei's MFOZ planning (2016-2020), the development intensity of the COZ is more than 90%, with approximately 6700 people per square kilometer, which is between five and ten times the average level of the whole city. From 2011 to 2020, construction land was dominant in the COZ, accounting for more than 94% of the area. In the IDZ, the main landscape is cultivated land, which accounts for more than 61% of the study area. Construction land was another important component of the IDZ, accounting for more than 23% of the total. Cultivated land accounted for approximately 80% of the ADZ, whereas construction land accounted for approximately 14%. The ECZ had the lowest level of urbanization among the four MFOZs. The coverage rates of cultivated land, construction land, water, and vegetation were approximately 53%, 9%, 28%, and 8%, respectively. The average LST differences between each MFOZ were calculated (Figure 7), and the day and night data for the four seasons from 2011 to 2020 were selected to calculate the contributions of the four MFOZs ( Figure 8). Overall, during the day, the average temperatures in the COZ and IDZ were higher than those in the ADZ and ECZ from 2011 to 2020 ( Figure 6). The average temperature in the study area was very close to that of the ADZ. Compared with spring and summer, the average temperature regularity of the MFOZ in autumn and winter was not obvious. This is because, in autumn and winter, part of the cultivated land is in the non-cultivation period, the landscape type is greatly affected by solar radiation, and the temperature is unstable. At night, the average temperatures of the COZ and ECZ were higher than those of the IDZ and ADZ. The average temperature in the study area was between that of the ECZ and IDZ. This is because the ECZ contains a large amount of water, and the large specific heat capacity of the water slows down the cooling rate at night, so the temperature in the ECZ is higher at night. During the day or night, the average temperature of the COZ was much higher than those of the other MFOZs. When comparing the difference between day and night, LST was relatively stable and more regular during the night. Because of the different thermal inertia of the landscape during the day, the influence of solar radiation is also different, whereas the artificial The average LST differences between each MFOZ were calculated (Figure 7), and the day and night data for the four seasons from 2011 to 2020 were selected to calculate the contributions of the four MFOZs ( Figure 8). Overall, during the day, the average temperatures in the COZ and IDZ were higher than those in the ADZ and ECZ from 2011 to 2020 ( Figure 6). The average temperature in the study area was very close to that of the ADZ. Compared with spring and summer, the average temperature regularity of the MFOZ in autumn and winter was not obvious. This is because, in autumn and winter, part of the cultivated land is in the non-cultivation period, the landscape type is greatly affected by solar radiation, and the temperature is unstable. At night, the average temperatures of the COZ and ECZ were higher than those of the IDZ and ADZ. The average temperature in the study area was between that of the ECZ and IDZ. This is because the ECZ contains a large amount of water, and the large specific heat capacity of the water slows down the cooling rate at night, so the temperature in the ECZ is higher at night. During the day or night, the average temperature of the COZ was much higher than those of the other MFOZs. When comparing the difference between day and night, LST was relatively stable and more regular during the night. Because of the different thermal inertia of the landscape during the day, the influence of solar radiation is also different, whereas the artificial heat resources during the night are relatively concentrated and stable. Overall, the average temperature of the four MFOZs after 2017 was lower than that before 2017, indicating that the implementation of MFOZ planning alleviated the urban heat island effect to a certain extent.
Land 2022, 11, x FOR PEER REVIEW 12 of 20 ranges from 0.01 to 0.05. The CI of the IDZ mostly ranges from −0.01 to −0.25. The CI of the ADZ mostly ranges from −0.03 to −0.51. This is because the larger specific heat capacity of water slows the cooling rate at night.

Contribution of Landscape Types under MFOZ
Overall, the absolute CI values of the landscape types in the IDZ and ECZ were much higher than those in the other functional zones (Figure 9). In the COZ, construction land

Contribution of Landscape Types under MFOZ
Overall, the absolute CI values of the landscape types in the IDZ and ECZ were much higher than those in the other functional zones (Figure 9). In the COZ, construction land had a significant positive impact on the thermal environment during the daytime in  Figure 8 shows the thermal environment contribution index of the MFOZ to the study area during the day and night in the different seasons. The thermal environment contribution index is defined as the product of the average temperature difference between the MFOZ and the study area and the area percentage; thus, the value of the area and the average temperature have a greater influence. During the day, IDZ contributed the most to the urban thermal environment. The CI of IDZ ranges from 0.03-0.62. The CI of the ECZ ranges from −0.69 to −0.03. The CI values of the ADZ and COZ fluctuate around zero. The ADZ had a constant negative contribution in the summer because of the cooling effect of the lush crops. Although the COZ had the highest LST, its contribution to the urban thermal environment was not the greatest because it had the smallest area, accounting for only 0.01% of the study area. The ECZ had a constant negative contribution, mainly because it had a large amount of vegetation and water, alleviating the effect of the urban heat island. However, the contribution of the ECZ changed significantly following the changes in water and vegetation during the four seasons. In spring, the cooling effect of water and vegetation was the most obvious in the ECZ. At night, ECZ contributed the most to the urban thermal environment. The CI of the ECZ mostly ranges from 0.01 to 0.60. The CI of COZ ranges from 0.01 to 0.05. The CI of the IDZ mostly ranges from −0.01 to −0.25. The CI of the ADZ mostly ranges from −0.03 to −0.51. This is because the larger specific heat capacity of water slows the cooling rate at night.

Contribution of Landscape Types under MFOZ
Overall, the absolute CI values of the landscape types in the IDZ and ECZ were much higher than those in the other functional zones (Figure 9). In the COZ, construction land had a significant positive impact on the thermal environment during the daytime in spring, summer, and autumn. Simultaneously, the cultivated land and water alleviated the regional thermal environment. At night, construction land had a greater positive impact on the four-season regional thermal environment, while cultivated land and water had a negative impact. Overall, the COZ, with construction land as the main source, contributes the most heat to the regional thermal environment during the day and night. Because the zone covers more than 94% of the impervious surface, there is less energy loss when latent heat evaporates.

Estimation of UTFVI
Thermal comfort in the environment is a way for the UTFVI to characterize urban health and ecological quality. In this study, the UTFVI evaluation was conducted in summer during the day when the most obvious heat island effect occurred. According to the analysis in Figure 10, the spatial distribution characteristics of the UTFVI in Hefei were significantly different from those of SUHI and LULC. By comparing and analyzing the In the IDZ, regardless of the time of day, construction land played a positive role in the regional thermal environment in all seasons. In particular, construction land contributed the most to the summer regional thermal environment during the daytime, at 0.61. At the same time, the average LST of construction land is more than 0.24 • C higher than the regional average LST. Cultivated land, vegetation, and water play negative roles during the day, particularly in spring and summer. At night, cultivated land played a negative role, and the CI of water and vegetation fluctuated around zero. Cultivated land, vegetation, and water can alleviate this regional thermal environment. The three landscape types have similar contributions to the urban thermal environment, but there are some differences. Because cultivated land is lush in summer and accounts for more than 61% of the area, the cooling effect in summer is the strongest during the daytime.
In the ADZ, cultivated land accounted for approximately 80% of the total area. During the daytime, cultivated and construction land had a positive impact on the regional thermal environment, and the impact of construction land was higher than that of cultivated land. Water and vegetation negatively affect regional thermal environments. At night, the impact of the landscape in different seasons on the regional thermal environment varied greatly, and the regularity was not strong. Cultivated land has a small positive impact in summer but a negative impact at night in other seasons. Vegetation had a positive impact. This is because the thermal inertia of cultivated land is low, and its ability to resist changes in the surface temperature is poor. However, vegetation has high thermal inertia. Therefore, the LST of cultivated land changes significantly during the day and night, and the vegetation maintains a higher temperature at night.
In the ECZ, cultivated land and water play important roles in the regional thermal environment, accounting for approximately 53% and 28% of the total area, respectively. During the day, cultivated land had a positive impact on the regional thermal environment, whereas water had a negative impact. At night, the impacts of cultivated land and water are opposite to those during the daytime. Construction land and vegetation had little impact on the thermal environment because of their small proportions. At night, in summer, the CI of cultivated land and water was close to zero. The seasonal difference was caused by the difference in thermal inertia. Cultivated land was the main contributor to the increase in the regional thermal environment, whereas water was the main source of cooling during the day. At night, the roles of cultivated land and water were opposite to those during the day. This is because the larger specific heat capacity of water slows down the cooling rate at night; therefore, the temperature is higher in water than in other terrains.

Estimation of UTFVI
Thermal comfort in the environment is a way for the UTFVI to characterize urban health and ecological quality. In this study, the UTFVI evaluation was conducted in summer during the day when the most obvious heat island effect occurred. According to the analysis in Figure 10, the spatial distribution characteristics of the UTFVI in Hefei were significantly different from those of SUHI and LULC. By comparing and analyzing the distribution of LULC and SUHI in the study area, we found that the spatial distribution pattern of the heat island zone in Hefei is consistent with construction land, whereas the distribution of water, cultivated land, and vegetation corresponds to the non-heat island zone. However, it can be seen from the UTFVI images that there are two extremes in the eco-environmental quality of Hefei: the zones with excellent and the worst ecological evaluation indicators are the most extreme; the zones with good, normal, bad, and worse ecological evaluation indicators are less extreme. The zones with the worst UTFVI were mainly distributed in the COZ and IDZ from 2011 to 2020. The zones with excellent UTFVI were mainly distributed in the ADZ and ECZ. The high coverage of water and vegetation in these zones indicates that they play a role in alleviating the thermal environment and improving the ecological environment. The UTFVI of the Dashu Mountain, Dafangying and Dongpu reservoirs, and Baiyan Lake in the main city zone (Figure 1) are excellent, showing that the large ecological parks in the main city zone have played a role in lowering the temperature, alleviating the urban heat island effect, and improving the ecological environment. Although we can see from the LULC and SUHI images since 2014 that the heat island effects of Chaohu, Changfeng, and Lujiang have not been significant, it can be seen from the UTFVI that these zones are under the worst UTFVI phenomenon, indicating that their ecological environment is gradually deteriorating.
Land 2022, 11, x FOR PEER REVIEW 15 of 20 that the ecological environment improved to a certain extent after the implementation of MFOZ planning, and the standard of living in cities in the ADZ and ECZ zones was better.

Discussion
Monitoring the dynamics of the urban thermal environment in the MFOZ generally requires both high spatial and temporal resolution. It is difficult to monitor the regional thermal environment because of the trade-off between spatial and temporal resolutions. In this study, we first used spatiotemporal reconstruction fusion [41] to obtain multi-time scale LSTs and then combined this with the CI, UTFVI, and landscape patterns to analyze the urban thermal environment and mechanism of the MFOZ. To understand the specific results after the implementation of MFOZ planning (Table 4), we analyzed the characteristics of the ecological indicators in each MFOZ. From 2011 to 2020, the zones with the most severe UTFVI accounted for 39.70%, 37.84%, 38.27%, and 36.72% of Hefei City, respectively, and the zones with excellent UTFVI accounted for 45.20%, 46.89%, 48.49%, and 50.96%, respectively. Overall, the thermal comfort in the environment has improved over the past decade. We further analyzed the UTFVI of each MFOZ and found that thermal comfort was improved for all MFOZs, except the COZ, which decreased thermal comfort. However, because its area percentage was only 1%, its impact on the UTFVI in the study area was not obvious. The analysis of the UTFVI showed that the ecological environment improved to a certain extent after the implementation of MFOZ planning, and the standard of living in cities in the ADZ and ECZ zones was better.

Discussion
Monitoring the dynamics of the urban thermal environment in the MFOZ generally requires both high spatial and temporal resolution. It is difficult to monitor the regional thermal environment because of the trade-off between spatial and temporal resolutions. In this study, we first used spatiotemporal reconstruction fusion [41] to obtain multi-time scale LSTs and then combined this with the CI, UTFVI, and landscape patterns to analyze the urban thermal environment and mechanism of the MFOZ.
According to the SUHI results, the fused LST image series at the Landsat scale reveals clear spatial details in the MFOZ, which is of great significance for the accurate classification and dynamic monitoring of the urban thermal environment. A remarkable SUHI core was found in the main city zone of Hefei, and the progress of urbanization can also be recognized. The SUHI have obvious differences, and the differences in the degree of urbanization determine the different characteristics of LST [59]. During the daytime, the zone with dense buildings in the main city zone belongs to the heat island zone, whereas the zones with ecological parks belong to the non-heat island zone. This is because the compact block has a larger proportion of hard ground, absorbs more heat, and has poor ventilation, whereas ecological parks have better air circulation and higher vegetation coverage [60,61]. Both transpiration and shade effects of vegetation during the day help alleviate the heat island effect [6]. Although the surrounding surface temperatures of water and vegetation near or within towns are relatively high, their thermal inertia is high; therefore, their surface temperatures are relatively low [4,54]. At night, areas where buildings are sparse and belong to non-heat island zones expand. This is because the heat island effect is primarily caused by the release of stored heat at night. Areas where buildings are sparse store less heat during the day and have better ventilation conditions; therefore, the cooling rate is faster than in dense areas at night [62]. The results showed that the natural surface area (vegetation and water) interspersing or encircling an area has an obvious cooling effect, which efficiently alleviates the nighttime heat island effect.
According to the results of the CI and UTFVI, the role of each MFOZ can be determined through the analysis of the contribution of different MFOZs to the urban thermal environment. During the daytime, the zones of the COZ, IDZ, and ADZ contributed significantly to the regional thermal environment in different seasons. At night, the COZ and IDZ zones contributed significantly. In general, the COZ and IDZ were the main sources of LST rise. In the ADZ and ECZ, the roles of cultivated land and vegetation change between day and night. During the daytime, cultivated land was the main contributor to regional thermal energy. In contrast, forests had a significant negative effect on regional LST, which was reversed at night. This is because vegetation has higher thermal inertia than cultivated land, causing it to maintain a higher temperature at night [54]. Although the water area is limited, it has high specific heat capacity and stable thermal environment properties, which can effectively isolate the heat island and regulate the temperature in local areas [22]. Large ecological parks effectively cool the surfaces of the COZ and IDZ [54,63]. To find an appropriate balance between sustainable urban development and the increase in urban green space, urban planners should optimize the allocation of green space patches in selected zones by increasing the size of existing green space patches rather than building new small patches [64,65]. In addition, the estimation of UTFVI further revealed zones with high thermal risk in the different MFOZs. Comparing the analysis results of the whole city and different main function zones, we found that research on different main function zones can reveal spatial differences in the influence of different regional factors [47]. For research areas with strong landscape heterogeneity, the study of different MFOZs better reflects the impact of various factors on the thermal environment so that decision-makers can select governance zones when planning policies or dividing the key restoration areas of the thermal environment [32,47].
However, several limitations and constraints should be addressed in future studies. First, many factors affect the LST, such as climate, precipitation, and wind speed. There were differences in LST during the same season or even in the same month. Our study selected only one LST scenario for each month, which may be one-sided. Second, there are certain errors in the LST inversion, reconstruction, fusion, and LULC classification. Third, there are only four land-use types in our study, which leads to less refinement of the role of more land types in the thermal environment. Fourth, the UTFVI is the spatial temperature difference normalized by the average temperature, which strongly depends on the background air temperature and has certain limitations. We acknowledge that these uncertainties limited the accuracy of this study. Nevertheless, our research provides a feasible method for exploring the multi-time scale evolution characteristics and influencing factors of the thermal environment in different MFOZs on a fine scale, fills the knowledge gap in urban thermal environment estimation, and provides a reference for effective planning of MFOZs in the future. In our future work, we will focus on solving the following problems. First, we will continue to increase the density of the Landsat-scale LST and analyze the changes and impacts of the urban thermal environment in different months in the same season or even at different times in the same month. Second, the role of vegetation in the estimation of UTFVI is unclear. In the next stage, we will further refine the study zone and land-use types to better explore their roles in developing a more effective urban planning scheme. Third, we will use a variety of methods to define SUHI and its intensity and explore their differences and accuracy.

Conclusions
In this study, Landsat and MODIS data were integrated using spatiotemporal reconstruction and fusion to derive high spatiotemporal resolution data, which were used to capture accurate fine-spatial resolution details. By analyzing the diurnal and seasonal CI of different landscape types under the MFOZ, we determined the role of each MFOZ in the thermal environment, and the landscape components in different MFOZs directly determined the heat absorption and heat release processes, especially in the COZ and IDZ. The estimation of UTFVI further revealed that the COZ and IDZ are zones with high thermal risk. For the entire study area, the thermal environment was alleviated after the implementation of MFOZ planning, among which the planning results of the ADZ and ECZ were relatively obvious. Therefore, to alleviate the heat island effect of different MFOZs, urban forests, ecological parks, and artificial lakes should be centrally planned at a certain scale in the COZ and IDZ to break up continuous and concentrated construction land. Protective measures should be taken for the original forestland, lakes, rivers, and other water in the ADZ and ECZ, and the expansion and infiltration of construction land to cultivated land should be restricted to maintain its role in regulating temperature. Although the study results were derived from Hefei, the methodology illustrated here can be readily applied to many other cities to explore whether the urban thermal environment is related to different geographical environments and development models.