Spatiotemporal Differentiation of Land Surface Thermal Landscape in Yangtze River Delta Region, China

: Advancements in the integrated development of the Yangtze River Delta are changing the structure and function of the surface thermal landscape and triggering a series of ecological and environmental problems. Therefore, examining the spatiotemporal differentiation characteristics and evolution laws of this land surface thermal landscape has great theoretical and practical signiﬁcance in the context of optimizing functional zoning and realizing the harmonious development of the economy, society and nature. The paper uses the LST (land surface temperature) data retrieved by MODIS (MOD11A2) remote sensing satellites in 2007, 2010, 2013, 2016 and 2019 to extract a land surface thermal rating map of the Yangtze River Delta region, and to analyze the spatiotemporal differentiation in the land surface thermal landscape, combining of the land surface thermal landscape strip proﬁle and thermal landscape pattern indices. The results show that the LST in the Yangtze River Delta region has increased in the past 12 years, the proportion of middle-, sub-high- and high-temperature zones increased by 33.42%, and the high-temperature zone has gradually extended into inland areas. The high-temperature zones in the areas surrounding core cities such as Shanghai, Nanjing, and Hangzhou have expanded. The corridor effect of thermal changes on the surface is obvious. The degree of aggregation in the lower-temperature areas has gradually decreased. The degree of aggregation in the higher-temperature regions has increased. The patch types of thermal landscape pattern increase, and the distribution of landscape area among various types tends to be even. this trend is most signiﬁcant in optimized development region.


Introduction
The deteriorating thermal environment in urban spaces has become one of the most significant features of global urban climate change [1,2]. Urban expansion will inevitably lead to significant changes in the spatial patterns of land use/cover change (LUCC) in urban areas [3][4][5], causing the land surface temperature (LST) to rise [6][7][8], which in turn affects the material structure and energy flow of urban ecosystems, and thus the quality of urban human settlements [9].
The Yangtze River Delta displays the highest level and fastest speed of urbanization of any region in China, and it is also core to international competition. Its development and evolution are the epitome of the integration and urbanization of China's developed regions [10,11]. The evolution of this area reflects the complex processes by which human activities transform surface structure and function, which in turn affect the spatiotemporal evolution of the surface thermal landscape [12]. Therefore, reviewing the mechanisms of the spatiotemporal changes in this region's surface thermal landscape is very important to formulating measures to alleviate the urban heat island effect [13][14][15], and to thus improve the modern urban ecological environment [16,17].
Urban space thermal environment research concerns the two key concepts of urban atmospheric heat island, and urban surface heat island [18]. With the continuous improvement of remote sensing and geographic information system (GIS) technology [19], many scholars have used remote sensing images to assess the LST and to characterize the urban thermal environment and urban heat island effect. Domestic scholars such as Chen et al. have used the research methods of landscape ecology for the first time, and in doing so proposed the concept of the thermal landscape, established a thermal landscape evaluation system, obtained urban surface temperature based on TM (land thematic mapper) images, and studied the spatial pattern of the thermal environment in Shanghai [20]. Qiao et al. used MODIS (moderate resolution imaging spectroradiometer) surface temperature data to characterize the evolution of Beijing's surface thermal landscape from the perspective of quantity, shape, and structure, using the landscape index, a surface thermal grade change map, thermal landscape centroid migration, and other methods [21]. Yue used SPOT images to interpret the types of urban land use, while also using ETM+ (Enhanced Thematic Mapper Plus) images from a Landsat7 satellite to retrieve LST. By analyzing the relationship between land use spatial patterns and LST, they revealed the internal mechanism of the urban heat island effect [22]. Gebrekidan studied the relationship between LUCC, normalized difference vegetation index (NDVI), surface emissivity and surface temperature for a city in Ethiopia. The results showed that LUCC is the main factor that triggers increases in surface temperature [23]. There have also been many studies on the thermal environment in the Yangtze River Delta. For example, Xia analyzed the urban expansion of Wuxi and its heat island response [24]; Zhang et al. systematically described the spatiotemporal change trajectory and the characteristics of the thermal environment pattern in Nanjing, and revealed the driving factors of these changes [25]; Han et al. studied the changing characteristics of the thermal environment security patterns of 16 cities in the Yangtze River Delta [26]. Overall, most of the existing studies have had a small scale, or have focused on a single city and thus have shown less consideration of differences within a region. Therefore, from a meso-scale perspective, this paper uses the thermal landscape index to analyze the spatiotemporal differentiation characteristics of the surface thermal landscape in the Yangtze River Delta, and further identifies the relationship between the spatial structures of urban agglomerations (core area-peripheral area; core city-node city-general city; corridor-hinterland) in the region and the thermal landscape of its surface.

Regional Overview of the Study
The Yangtze River Delta is located in the lower reaches of the Yangtze River, east of the Yellow Sea and East China Sea. It is an alluvial plain formed by the Yangtze River entering the sea. It has a subtropical monsoon climate, with high temperatures and rain in the summer, mild rain in the winter, and small annual temperature differences. Its geographical position is highly favorable-the terrain is flat, and transportation is convenient, which is conducive to the development of external relations and imparts a broad geographical advantage on the hinterland. The Yangtze River Delta region studied in this paper incorporates Shanghai, Jiangsu Province, Zhejiang Province and Anhui Province (three provinces and one city), with 41 prefecture-level cities and a total area of 358,000 km 2 ( Figure 1).
The Yangtze River Delta's superior geographical location means its total economic output ranks among the highest in the country, with a solid economic foundation, developed manufacturing and high-tech industries, and a rapidly developing service industry. It only accounts for 2.5% of the country's total area, but it contains 16% of the country's population and contributes 24% of the total economic output (2019). The Yangtze River Delta has been integrated since 1982, and its scope of integration has increased since then (Table 1). In 2018, the integrated development of the Yangtze River Delta became a national strategy. From 2007 to 2019, the population density, power consumption, and proportion of built-up area in the Yangtze River Delta region showed upward trends ( Figure 2). Changes The Yangtze River Delta's superior geographical location means its total economic output ranks among the highest in the country, with a solid economic foundation, developed manufacturing and high-tech industries, and a rapidly developing service industry. It only accounts for 2.5% of the country's total area, but it contains 16% of the country's population and contributes 24% of the total economic output (2019). The Yangtze River Delta has been integrated since 1982, and its scope of integration has increased since then (Table 1). In 2018, the integrated development of the Yangtze River Delta became a national strategy. From 2007 to 2019, the population density, power consumption, and proportion of built-up area in the Yangtze River Delta region showed upward trends ( Figure  2). Changes in the urban environment have caused significant changes in the surface temperature of the region.

2010
The Regional Plan for the Yangtze River Delta Region defines the scope of the Yangtze River Delta as 25 prefecture-level cities in Shanghai, Jiangsu Province and Zhejiang Province. The plan still lists 16 cities as the "core areas" of the Yangtze River Delta regional development plan.

2010
The Regional Plan for the Yangtze River Delta Region defines the scope of the Yangtze River Delta as 25 prefecture-level cities in Shanghai, Jiangsu Province and Zhejiang Province. The plan still lists 16 cities as the "core areas" of the Yangtze River Delta regional development plan.

2014
The Guiding Opinions of the State Council on Relying on the Golden Waterway to Promote the Development of the Yangtze River Economic Belt stated that the Yangtze River Delta City Group should build an urban agglomeration centered on Shanghai and Nanjing, and with Hangzhou and Hefei as the sub-centers. Anhui was officially included in the Yangtze River Delta region.

2016
The National Development and Reform Commission issued the Development Plan for the Yangtze River Delta Urban Agglomeration, which included 26 cities.

2019
The Outline of the Yangtze River Delta Regional Integration Development Plan was given, and its planning scope covered the whole area of Shanghai, Jiangsu Province, Zhejiang Province, and Anhui Province, with 27 cities constituting the central area.

2019
The Outline of the Yangtze River Delta Regional Integration Development Plan was given, and its planning scope covered the whole area of Shanghai, Jiangsu Province, Zhejiang Province, and Anhui Province, with 27 cities constituting the central area.

Remote Sensing MODIS Date
MODIS is the abbreviation of the moderate resolution imaging spectroradiometer. It is a detection instrument mounted on Terra and Aqua satellites. The ground resolution of MODIS data is 250, 500 and 1000 m [27], the scan width is 2330 km, and the observation data contain 36 discrete spectral bands [28]. This multi-wavelength observational equipment can closely scrutinize the dynamic changes of various features on the earth, such as changes in surface temperature [29], vegetation coverage [30], aerosol and atmospheric properties [31][32][33][34], etc. This paper employs the LST product retrieved by the land standard product MOD11A2 on the Aerra satellite (data source: https://ladsweb.modaps.eosdis.nasa.gov/search/ (accessed on 15 October 2020)). The spatial resolution is 1000 m, and the product contains 8d synthetic data. On average, four images are taken per month. As the weather in August is clear and cloudless, the ground features are clear and the quality of remote sensing images is good. Studies have shown that summer temperatures in urban areas are high, the heat island is strong, and its area is large; the calculation results are easily analyzable [34,35]. Therefore, this article uses the same four days (217th, 225th, 233rd and 241st day) for each year (2007,2012,2015 and 2019 year). We referred these to a MODIS China regional rank number, and selected the three areas of H27V05, H28V05, and H28V06 to cover the Yangtze River Delta.

MODIS Data Preprocessing
First, we used the MRT (MODIS Reprojection Tool) to splice the data of the three scenes, selected the band of LST_Day_1 km, and added the geographic coordinate system

Remote Sensing MODIS Date
MODIS is the abbreviation of the moderate resolution imaging spectroradiometer. It is a detection instrument mounted on Terra and Aqua satellites. The ground resolution of MODIS data is 250, 500 and 1000 m [27], the scan width is 2330 km, and the observation data contain 36 discrete spectral bands [28]. This multi-wavelength observational equipment can closely scrutinize the dynamic changes of various features on the earth, such as changes in surface temperature [29], vegetation coverage [30], aerosol and atmospheric properties [31][32][33][34], etc. This paper employs the LST product retrieved by the land standard product MOD11A2 on the Aerra satellite (data source: https://ladsweb.modaps.eosdis. nasa.gov/search/ (accessed on 15 October 2020)). The spatial resolution is 1000 m, and the product contains 8d synthetic data. On average, four images are taken per month. As the weather in August is clear and cloudless, the ground features are clear and the quality of remote sensing images is good. Studies have shown that summer temperatures in urban areas are high, the heat island is strong, and its area is large; the calculation results are easily analyzable [34,35]. Therefore, this article uses the same four days (217th, 225th, 233rd and 241st day) for each year (2007,2012,2015 and 2019 year). We referred these to a MODIS China regional rank number, and selected the three areas of H27V05, H28V05, and H28V06 to cover the Yangtze River Delta.

MODIS Data Preprocessing
First, we used the MRT (MODIS Reprojection Tool) to splice the data of the three scenes, selected the band of LST_Day_1 km, and added the geographic coordinate system (WGS84); secondly, the band calculator (Band Math) was used to synthesize the data for the 217th, 225th, 233rd, and 241st days to derive the August surface temperature using the maximum value synthesis method, and we cropped the synthesized remote sensing data according to the vector boundary of the Yangtze River Delta In ENVI5.3. Thirdly, the image was reprojected (Reprojection coordinates: WGS1984 UTM Zone 50N); then, we converted the LST data from Kelvin to Celsius using the description of MODIS data [36]. The specific formula is shown in (1). After the temperature conversion, we derived the DN (Digital where MPVK is MODIS pixel value in Kelvin, and MPVC is MODIS pixel value in Celsius. The function of data normalization is to summarize the statistical distribution of a unified sample. As such, the range of land temperature values is distributed between 0 and 1 [37] for further research and statistical analyses. We used the band calculator [38] to enter the Formula (2) into ENVI5.3: where T * i is the value after the normalization of the i pixel, T i is the actual land temperature of the i pixel, T max is the maximum LST in all the pixels, and T min is the minimum LST in all the pixels.
After normalization, using the reclassification tool [39] in ArcGIS 10.2, the standardized LST was divided into five levels ( Table 2), based on which the classification map of the surface thermal landscape of the Yangtze River Delta region from 2007 to 2019 was produced ( Figure 3). Table 2. Classification of surface thermal rating.

Thermal Landscape Rating
Temperature Zone Range High-temperature zone [0.8, 1.0]

Strip Profile of the Surface Thermal Landscape
The strip profile extraction method uses the principle of mathematical statistics to obtain elevation information for a buffer area in a direction perpendicular to the profile line [40]. This paper employs the strip profile extraction method to further analyze the changes in the surface thermal landscape of different regions in the Yangtze River Delta. The Development Plan for the Yangtze River Delta Urban Clusters proposes a networked spatial pattern of "one core, five circles and four belts", wherein the four belts are the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt, the Yangtze River Development Belt, the Coastal Development Belt, and the Shanghai-Hangzhou-Jinhua Development Belt. These belts link 27 cities (Table 3). We used ArcGIS10.2 to extract surface thermal landscape strip profiles of the four development zones in 2007 and 2019, based on the classification map of the surface thermal landscape.

Strip Profile of the Surface Thermal Landscape
The strip profile extraction method uses the principle of mathematical statistics to obtain elevation information for a buffer area in a direction perpendicular to the profile line [40]. This paper employs the strip profile extraction method to further analyze the changes in the surface thermal landscape of different regions in the Yangtze River Delta. The Development Plan for the Yangtze River Delta Urban Clusters proposes a networked spatial pattern of "one core, five circles and four belts", wherein the four belts are the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt, the Yangtze River Development Belt, the Coastal Development Belt, and the Shanghai-Hangzhou-Jinhua Development Belt. These belts link 27 cities (Table 3). We used ArcGIS10.2 to extract surface thermal landscape strip profiles of the four development zones in 2007 and 2019, based on the classification map of the surface thermal landscape.

Thermal Landscape Pattern Indices
The thermal landscape pattern indices draw lessons from the research methods of landscape ecology and analyzes the evolution characteristics of thermal landscape through the calculation of landscape indices. The landscape pattern indices contain a wealth of landscape pattern information, which reflects the number, shape, and spatial features of each landscape area [21,41,42]. We used Fragstats 4.2 to calculate each landscape pattern indices based on the classification map of the surface thermal landscape (Table 4).

1.
The AI concerns the aggregation of patches of a certain surface heat level [41]. It calculates the length of the common boundary between pixels in a certain level of the surface thermal landscape. When the common boundary value for all patches in the same level is largest, the aggregation index is largest, which means that the level of surface thermal integration is good. In some studies, the landscape distribution is relatively concentrated, while in others it is more scattered [43,44].

2.
The CONTAG represents the trend in the extension of different surface heat levels throughout the entire landscape, and the discrete relationship between different patches [45]. When the CONTAG index increases, the connections between the dominant patches in the landscape are better. On the contrary, there are many types of patches in the landscape, and they are densely distributed, and the degree of fragmentation in the landscape is high. 3.
The SHDI reflects the richness and complexity of the landscape pattern. When SHDI = 0, the landscape pattern is composed of patches of the same type. When the SHDI index increases, this means that the number of patch types increases, or the patch types are more evenly distributed in the landscape [46,47]. 4.
The SHEI reflects the uniformity of the distribution of patches in the landscape. In the study of surface heat, when the SHEI index increases, this indicates that heat levels of different types are homogeneously distributed throughout the entire landscape, and when it decreases, they are concentrated in some areas [47]. Table 4. Landscape indices calculation and feature description.

Landscape Indices Calculation Formula Variable Interpretation
Aggregation Index (AI) g ij is the connection between all surface thermal regions in the i level thermal landscape, ax→g ii is the maximum number of connections between all pixels in the i level thermal landscape, and P i is the ratio of the i level thermal landscape pixels to the entire landscape area.
m is the number of thermal levels in the entire landscape, P ij represents the number of nodes between the i level surface heat and the j level surface heat. Shannon's Diversity Index (SHDI) [43,44] n is the number of all pixels in the landscape, and P i is the ratio of i level surface heat to the total landscape area.
Simpson's Diversity n is the number of all pixels in the landscape, and P i is the ratio of i level surface heat to the total landscape area.

The Universality of LST Rise
In 2007, about 85.54% of the Yangtze River Delta region was covered by sub-low-and middle-temperature zones; the sub-low-temperature zone accounted for about 42.60%, and the middle-temperature zones accounted for about 42.94%. The sub-high-and hightemperature zones accounted for only 14.46% (Figures 3 and 4), mainly distributed in Shanghai, Nanjing, Hangzhou, Hefei, Suzhou (Jiangsu), Wuxi, Ningbo and other urban areas, forming scattered and independent patches; the sub-low-and low-temperature zones were mainly distributed in Jiangsu and northern Anhui, the hilly and mountainous areas in western and southern Anhui, the hills of western and eastern Zhejiang, the mountains of southern Zhejiang, and the regions with rivers and lakes, which are distributed continuously over a large area.
In 2010, the temperatures of the low-temperature zones decreased slightly, mainly located in northern Jiangsu, the thermal rating of which changed from low-temperature to sub-low-temperature, and the area of sub-low-temperature zones decreased significantly (to 28.09%). The sizes of middle-, sub-high-and high-temperature zone have increased compared with 2007, and their proportions have increased to 53.83%, 12.02%, and 2.05%, respectively. Most of the sub-low-temperature zones surrounded by middletemperature zones has been upgraded to middle-temperature zones. The sub-high-and high-temperature zones are continuing to expand towards the peripheries of the cities, and a trend of continuous development has appeared.
In 2013, the areas of the sub-low-temperature zones continued to decrease, while the areas of the middle-, sub-high-and high-temperature zones continued to increase. The middle-temperature zones were more widely distributed, covering 58.19% of the total area of the Yangtze River Delta, and the mid-temperature area in southern Zhejiang has increased the most. The sub-high-and high-temperature zones in the western Hangzhou and Taihu Lake areas increased significantly. temperature zones accounted for only 14.46% (Figures 3 and 4), mainly distributed Shanghai, Nanjing, Hangzhou, Hefei, Suzhou (Jiangsu), Wuxi, Ningbo and other urb areas, forming scattered and independent patches; the sub-low-and low-temperatu zones were mainly distributed in Jiangsu and northern Anhui, the hilly and mountaino areas in western and southern Anhui, the hills of western and eastern Zhejiang, the mou tains of southern Zhejiang, and the regions with rivers and lakes, which are distribut continuously over a large area. In 2010, the temperatures of the low-temperature zones decreased slightly, main located in northern Jiangsu, the thermal rating of which changed from low-temperatu to sub-low-temperature, and the area of sub-low-temperature zones decreased sign cantly (to 28.09%). The sizes of middle-, sub-high-and high-temperature zone have creased compared with 2007, and their proportions have increased to 53.83%, 12.02%, a 2.05%, respectively. Most of the sub-low-temperature zones surrounded by middle-te perature zones has been upgraded to middle-temperature zones. The sub-high-and hig temperature zones are continuing to expand towards the peripheries of the cities, and trend of continuous development has appeared.
In 2013, the areas of the sub-low-temperature zones continued to decrease, while t areas of the middle-, sub-high-and high-temperature zones continued to increase. T middle-temperature zones were more widely distributed, covering 58.19% of the to area of the Yangtze River Delta, and the mid-temperature area in southern Zhejiang h In 2016, the proportion of low-temperature zones dropped slightly to 4.11%; the sublow-temperature zones continued to decrease in size; the proportion of middle-temperature zones remained largely unchanged, while the sub-high-and high-temperature zones increased significantly, and the proportion of sub-low-temperature zones increased from 13.77% to 20.52%. Most areas in northern Jiangsu have changed from low-and sub-lowtemperature zones to middle-temperature zones, but some areas in the mountainous region of southern Zhejiang came to show a lower thermal rating. The thermal rating of most areas around Taihu Lake increased, and the sub-high-and high-temperature zones were further strengthened. Around cities such as Nantong, Jinhua, Quzhou, Wenzhou and Taizhou (Zhejiang), high-temperature patches have gradually appeared. The warming effect in cities along the Yangtze River in Anhui is obvious, such as in Ma'anshan, Wuhu, Xuancheng, Tongling, Chizhou, Anqing, etc.
In 2019, the areas of the low-temperature zones remained unchanged, the areas of the sub-low-and middle-temperature zones decreased slightly, and the areas of the subhigh-and high-temperature zones continued to increase, to 22.0 and 6.43%, respectively. Compared with 2016, the area occupied by sub-high-temperature zones decreased to the west of Taihu Lake and along the Yangtze River in Anhui, but there was a significant increase in northern Anhui.
Overall, from 2007 to 2019, the proportions of low-and sub-low-temperature zones gradually decreased, and the proportions of middle-, sub-high-and high-temperature zones gradually increased. The low-and sub-low temperature zones were mainly distributed in the hills and mountains of northern Jiangsu, western and southern Anhui, western and eastern Zhejiang, southern Zhejiang, and around rivers and lakes. Northern Jiangsu is an important part of China's coastal economic belt, an important node area of the "the Belt and Road", an important region for economic growth in the Yangtze River Delta, and one of the fastest growing regions in China. Therefore, the effect of rising surface heat levels is relatively significant. The spatial pattern of urban development in the middle-, subhigh-and high-temperature zones exhibits a "Z"-shaped arrangement (Nanjing-Shanghai-Hangzhou-Ningbo) in the Yangtze River Delta; this has gradually expanded to peripheral areas and displays a continuous distribution trend.

The Urban Hierarchical Order of LST Changes
There is an obvious relationship between LST changes and urban hierarchy: node cities around core cities display the most obvious temperature rise. The high-calorific value areas of high-level cities are broad and contiguous, while the high-calorific value areas of low-level cities are small and isolated.
From 2007 to 2019, among the 27 cities in the Yangtze River Delta region, the proportions of low-and sub-low-temperature zones showed a decreasing trend, except for in Jinhua ( Figure 5). In total, 11 cities showed a decrease in their proportion of middletemperature zones, including one super city (Shanghai), one mega city (Nanjing), one type I big city (Suzhou), four type II big cities (Wuxi, Changzhou, Shaoxing, Taizhou (Zhejiang)), and four medium-sized cities (Jiaxing, Zhenjiang, Huzhou, Jinhua). On the other hand, 16 cities showed increasing proportions, including two type I big cities (Hangzhou, Hefei), seven type II big cities (Wenzhou, Ningbo, Nantong, Wuhu, Yancheng, Yangzhou, Taizhou (Jiangsu)), three medium-sized cities (Ma'anshan, Anqing, Zhoushan), and four type I small cities (Tongling, Chuzhou, Xuancheng, Chizhou). The proportions of lowand sub-low-temperature zones showed a generally increasing trend, except for in Ningbo and Zhoushan.
Sustainability 2021, 13, x FOR PEER REVIEW 10 of 21 (33.01%); except for Nanjing, which is a mega city, all of these are large and medium-sized cities. They are all distributed in peripheral areas around super cities and mega cities (Shanghai, Nanjing, Hangzhou), around Taihu Lake, or along the Yangtze River.

Corridor Effect of LST Changes
As shown in Figure 6a, in 2007, the high-calorific value areas of Shanghai and Hangzhou occupied a large stretch of the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt. As the main urban areas of Nanjing expanded along the Yangtze River, the high-calorific value areas of Nanjing along the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt slightly smaller, while the high-calorific value areas of Hefei and Changzhou were smaller than those of Shanghai and Hangzhou. The high-calorific value areas of cities such as Suzhou (Jiangsu), Wuxi, Ningbo, Shaoxing, Zhenjiang, and Jiaxing showed a narrower distribution, and each contain two or more high-calorific value centers. From 2007 to 2019, the surface thermal rating and the spatial distribution of highcalorific value areas along the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt showed an increasing trend, especially in the Nanjing-Shanghai-Hangzhou section, and the proportion of high-temperature areas increased significantly, which is con- We defined areas with a surface thermal rating >3 as high-calorific value areas. From 2007 to 2019, the nine cities with highest high-calorific value area ratio were Zhenjiang (71.14%), Jiaxing (53.70%), Changzhou (51.70%), Wuxi (48.44%), Suzhou (Jiangsu) (44.57%), Nantong (41.78%), Huzhou (37.90%), Nanjing (36.36%) and Taizhou (Jiangsu) (33.01%); except for Nanjing, which is a mega city, all of these are large and medium-sized cities. They are all distributed in peripheral areas around super cities and mega cities (Shanghai, Nanjing, Hangzhou), around Taihu Lake, or along the Yangtze River.

Corridor Effect of LST Changes
As shown in Figure 6a, in 2007, the high-calorific value areas of Shanghai and Hangzhou occupied a large stretch of the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt. As the main urban areas of Nanjing expanded along the Yangtze River, the high-calorific value areas of Nanjing along the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt slightly smaller, while the high-calorific value areas of Hefei and Changzhou were smaller than those of Shanghai and Hangzhou. The high-calorific value areas of cities such as Suzhou (Jiangsu), Wuxi, Ningbo, Shaoxing, Zhenjiang, and Jiaxing showed a narrower distribution, and each contain two or more high-calorific value centers. From 2007 to 2019, the surface thermal rating and the spatial distribution of highcalorific value areas along the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt showed an increasing trend, especially in the Nanjing-Shanghai-Hangzhou section, and the proportion of high-temperature areas increased significantly, which is consistent with the data for other cities where the high-calorific value has increased.
As shown in Figure 6b, relying on the Yangtze River Golden Waterway, the Yangtze River Development Belt is the main transportation corridor by which the Yangtze River Delta can expand inland, the main area of river-sea combined transportation below Nanjing, and the demonstration area of Wanjiang City Belt industrial transfer. In 2007, the cities with the largest areas with high calorific value were Shanghai and Nanjing, while around Yangzhou, Ma'anshan, Wuhu, Tongling, Chizhou, Anqing, etc., small areas of highcalorific value were formed. By 2019, Nanjing and the cities along the Yangtze River, such as Taizhou (Jiangsu), Yangzhou, and Ma'anshan, were warming significantly, forming a continuous region of high-temperature areas. Only inland areas such as Tongling, Chizhou, and Anqing were low-temperature areas.
As shown in Figure 6c, the Coastal Development Belt is relatively weak in the four regions, and especially northern Jiangsu and southern Fujian, that lack both deep-water ports and core cities; therefore, the north and south ends of the Belt are relatively deprived areas, except for the Hangzhou Bay coast. As reflected in the LST, in 2007, the large-scale high-calorific value areas were mostly found in Shanghai and Ningbo, followed by the type II big cities Yancheng, Nantong, Taizhou (Zhejiang), and Wenzhou. By 2019, the highcalorific value areas of Yancheng and Tongzhou significantly expanded, and the surface temperatures in coastal areas of northern Jiangsu and southern Fujian increased to varying degrees, although these areas still dominated by middle-temperature zones.
As shown in Figure 6d, the Shanghai-Hangzhou-Jinhua Development Belt relies on the Shanghai-Kunming Channel, connecting Shanghai, Jiaxing, Hangzhou, Jinhua, Quzhou, and other cities. The Hangzhou-Jinhua-Quzhou section is mainly located in the Jinqu Basin and proceeds northeast-southwest, and the high-calorific value areas in this region are large. In 2007, the large-scale high-calorific value areas were concentrated mostly in Jinhua, followed by Zhuji and Quzhou. By 2019, the sizes of the Zhuji, Jinhua, and Quzhou high-calorific value areas had increased, and high-calorific patches appeared between Zhuji and Yiwu and between Jinhua and Quzhou.
On the whole, the high-calorific value areas of cities in the Yangtze River Delta have continuously expanded over time, and expanded along four corridors, forming a trend of the continuous spread of high-temperature zones with a significant corridor effect. Out of the four corridors, along the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt and the Yangtze River Development Belt, which both have many cities connected in series, the surface heat level is generally higher, and the surface warming effect is obvious. In contrast, the Coastal Development Belt and the Shanghai-Hangzhou-Jinhua Development Belt contain fewer cities, with low levels and small scales of urban expansion, so the general surface heat here is low, and the surface warming effect is not significant. As shown in Figure 6b, relying on the Yangtze River Golden Waterway, the Yangtze River Development Belt is the main transportation corridor by which the Yangtze River Delta can expand inland, the main area of river-sea combined transportation below Nanjing, and the demonstration area of Wanjiang City Belt industrial transfer. In 2007, the cities with the largest areas with high calorific value were Shanghai and Nanjing, while around Yangzhou, Ma'anshan, Wuhu, Tongling, Chizhou, Anqing, etc., small areas of high-calorific value were formed. By 2019, Nanjing and the cities along the Yangtze River, such as Taizhou (Jiangsu), Yangzhou, and Ma'anshan, were warming significantly, forming a continuous region of high-temperature areas. Only inland areas such as Tongling, Chizhou, and Anqing were low-temperature areas.

Regional Differences of LST Changes
In the Development Plan for the Yangtze River Delta Urban Agglomeration, based on the main functional plan for the zone, the land development intensity, the development direction, and the suitability of the population agglomeration and urban-rural construction, the land is divided into three types: optimized, key, and restricted development regions (Figure 7). Different development regions manifest different trends in the Aggregation Index (AI), Contagion Index (CONTAG), Shannon's Diversity Index (SHDI), and Simpson's Diversity Index (SHEI).
upward trend year by year. The reason is that this type is mainly distributed in the areas around Taihu Lake, such as Suzhou, Wuxi, Jiaxing, and Huzhou. With the rapid development of the Taihu Lake area, especially the cities on the east side, the sub-low temperature patches in the marginal areas of the original cities have gradually been invaded by the high temperature patches. the distribution area of sub-low-temperature is becoming less and less, and the middle-or even the sub-high-temperature zone invades the east side of Taihu Lake, and the AI rises. As shown in Figure 9, the CONTAG is lowest in the optimized development regions and highest in the restricted development regions; the SHDI and SHEI are highest in the optimized development regions and lowest in the restricted development regions. This shows that the thermal landscape of the optimized development region is highly fragmented, while the overall thermal landscape of the restricted development zone is better, and the complexity and richness of the landscape type are not as good as the optimized development region. The key development zone is between the former two. Over the past 12 years, the CONTAG of the optimized development region first increased and then decreased, while the SHDI and SHEI first decreased and then increased, indicating that the degree of thermal landscape pattern fragmentation in the optimized development regions first decreased and then increased. The fragmentation of thermal landscape in key development regions fluctuated greatly. In 2013, the middle-temperature zone accounted for As shown in Figure 8, the AI of the high-and sub-high temperature zone showed restricted, key and optimized development regions, from high to low; the overall AI index in the middle-and low-temperature zones does not change significantly; the AI of the sub-low-temperature zones showed primarily optimized, restricted and key development regions, from high to low. Judging from the differences in changes in different temperature zone, the AI in the high-and sub-low-temperature zones has changed significantly in the three development regions in the past 12 years, while the AI in the other three temperature types has not changed significantly. Among them, the AI of the three development regions in the high-temperature zone has shown an upward trend year by year. The area where the high-temperature zone is located is generally the city center. With the rapid progress of urbanization in the Yangtze River Delta and the continuous expansion of urban dense areas, discrete high temperature patches are gradually connected through corridors, making the AI increase year by year. The performance of this trend optimized development zone is the most significant. The AI of sub-low-temperature zones in key and restricted development region shows a downward trend. This is mainly due to the fact that the key and restricted development regions have undertaken many enterprises that have been transferred from the original optimized development regions. As industries, population, and transportation continue to gather in cities and towns, the urban space continues to expand. Urban hightemperature grade patches are broken, leading to a downward trend in the AI, while the optimized development region as a whole shows an upward trend year by year. The reason is that this type is mainly distributed in the areas around Taihu Lake, such as Suzhou, Wuxi, Jiaxing, and Huzhou. With the rapid development of the Taihu Lake area, especially the cities on the east side, the sub-low temperature patches in the marginal areas of the original cities have gradually been invaded by the high temperature patches. the distribution area of sub-low-temperature is becoming less and less, and the middle-or even the sub-high-temperature zone invades the east side of Taihu Lake, and the AI rises.  71.31% of the area of the key development region, and the degree of fragmentation of the thermal landscape was minimized. After 2013, as the area of the middle-temperature zone decreased, its leading role in the thermal landscape decreased year by year, and the thermal landscape is gradually developing in the direction of diversification and fragmentation. The CONTAG of the restricted development region generally showed a downward trend, while the SHDI and SHEI show an upward trend, and the landscape pattern tended towards fragmentation. It shows that the types of urban thermal landscape patches are increasing, and the distribution of landscape area among various types tends to be even. This trend is most significant in optimized development zones.

Discussion
Global climate change and the intensification of urbanization have had a strong impact on the urban environment [9]. As an important aspect of these impacts, urban land surface thermal landscape has changed urban climate and air quality, and seriously threatened the health of urban residents. Therefore, the research on land surface thermal As shown in Figure 9, the CONTAG is lowest in the optimized development regions and highest in the restricted development regions; the SHDI and SHEI are highest in the optimized development regions and lowest in the restricted development regions. This shows that the thermal landscape of the optimized development region is highly fragmented, while the overall thermal landscape of the restricted development zone is better, and the complexity and richness of the landscape type are not as good as the optimized development region. The key development zone is between the former two. Over the past 12 years, the CONTAG of the optimized development region first increased and then decreased, while the SHDI and SHEI first decreased and then increased, indicating that the degree of thermal landscape pattern fragmentation in the optimized development regions first decreased and then increased. The fragmentation of thermal landscape in key development regions fluctuated greatly. In 2013, the middle-temperature zone accounted for 71.31% of the area of the key development region, and the degree of fragmentation of the thermal landscape was minimized. After 2013, as the area of the middle-temperature zone decreased, its leading role in the thermal landscape decreased year by year, and the thermal landscape is gradually developing in the direction of diversification and fragmentation. The CONTAG of the restricted development region generally showed a downward trend, while the SHDI and SHEI show an upward trend, and the landscape pattern tended towards fragmentation. It shows that the types of urban thermal landscape patches are increasing, and the distribution of landscape area among various types tends to be even. This trend is most significant in optimized development zones.

Discussion
Global climate change and the intensification of urbanization have had a strong impact on the urban environment [9]. As an important aspect of these impacts, urban land surface thermal landscape has changed urban climate and air quality, and seriously threatened the health of urban residents. Therefore, the research on land surface thermal

Discussion
Global climate change and the intensification of urbanization have had a strong impact on the urban environment [9]. As an important aspect of these impacts, urban land surface thermal landscape has changed urban climate and air quality, and seriously threatened the health of urban residents. Therefore, the research on land surface thermal landscape based on the process of urbanization has attracted extensive attention of scholars [48].
Urban agglomeration has become a distinctive feature worldwide over recent decades [49][50][51], it are an inevitable outcome of industrialization and urbanization. Compared with a single city, it has more complex functions and structures. An urban agglomeration can cause a significant change in the thermal environment over a wide region compared to that generated by individual cities [52][53][54]. However, most studies are limited to the perspective of a single city, ignoring the impact of urban agglomeration [55]. In previous research studies, it has been determined that urbanization is one of the significant drivers of LST alteration across the world [56][57][58].
However, so far, an investigation of the spatiotemporal variation of regional thermal environment in urban agglomerations has not been carried out in detail, little is known about the spatiotemporal patterns and evolution of the LST. Therefore, a detailed research of the spatiotemporal evolution of the LST in urban agglomerations are crucial for the government to take more environment-friendly action in order to address adapted plans for mitigation and adaptation to these issues [53]. It has been widely acknowledged that urban agglomeration significantly enhances the LST [59][60][61]. However, the effect of LST involves all stages of development: single city, mega city and city agglomeration. Studies concerning this phenomena have so far largely focused on large cities, or on a single area within a city [62].

Spatial Differentiation of Thermal Environment in Urban Agglomerations
Compared with a single city, urban agglomeration involves more cities and larger regions. Therefore, LST effect is not only reflected in cities, but also in regions. The spatiotemporal changes of LST effect are more complex. In urban agglomeration areas, not only the high-temperature zones around core cities and important node cities are expanding, but also the high-temperature zones of adjacent cities tend to integrate into one; Along the complex traffic corridor, the high-temperature zone is also gradually integrated, and the corridor effect is very obvious. In the Pearl River Delta Metropolitan Region (PRDR), the isolated urban heat islands were gradually connected and interacted with each other, forming the regional heat island [63]. This is consistent with the research results of this paper ( Figure 6). In Chengdu-Chongqing urban agglomeration, the heat island intensity in megacities was greater than other cities, and a significant correlation was observed between the built-up area and the daytime heat island intensity [64]. However, the high temperature center is not necessarily located in the city center, may be distributed in industrial areas and commercial intensive areas [65,66]. In Shanghai, except the city center, the intensive industrial parks at the urban fringe, docks and harbors, and major satellite towns outside the outer rings exhibited higher LST and remarkable variation in LST.
Many results show that LST generally decreases from the urban center to the suburb (or surrounding rural areas), and the low temperature is generally distributed in the area where a large area of green space and water body are located [67], this study also draws the same conclusion (Figures 3 and 4). Initial observations of LST in recent years have indicated that LST has tended to increase in both urban and suburban areas [68].
There are also spatial differences in thermal landscape pattern changes. The degree of aggregation of high-level thermal landscape types is on the rise, the distribution of thermal landscape area among various types gradually tends to be even [69]. Similar conclusions were reached in this study (Figures 8 and 9).

The Causes of Changes in Urban Thermal Environment
LST is a very important parameter in surface energy balance which associated with some dominant LUCC types. Urban agglomeration development has significantly transformed natural and semi-natural surfaces into impervious structures, which has disturbed the balance of the earth's surface radiation and energy [63,70,71], leads to large-scale LUCC that result in changes to physical properties [72]. Many results confirmed the most obvious dependence of the LST on different land cover types [73]. The surface thermal distribution is seriously affected by the impervious surface coverage and vegetation coverage [74].
Urbanization, which is a form of land cover, occupies less than 2% of the global land surface, evidence shows that human disturbance due to urbanization has significantly altered the natural landscape all over the world. Due to urbanization and industrialization, the number of built-up areas are increasing and replacing the vegetation and land cover area. Human activity has led to an increase in impervious surfaces, a decrease in vegetated areas, and a rise in greenhouse gas emissions [75][76][77], which has resulted in an increase in the occurrence of urban heat island globally. The urban heat island is a condition wherein the temperature within an urban agglomeration is higher than that of the surrounding rural localities. LUCC are significant factor for the spatiotemporal variation of LST [78][79][80]. Natural land surfaces are gradually replaced by impervious surfaces [81], this change leads to different exchange processes of water, substances and energy between the land surface and the atmosphere. A larger urban area generally has more impervious surfaces, which can absorb more heat and raise the urban temperature [82]. Larger urban areas always mean a developed regional economy and higher energy consumption, which can expand the temperature gap between urban and suburban areas and increase the LST intensity.
LST is closely related to land cover types, meteorological conditions, anthropogenic heat sources, the size of the urban area, urban morphology [83]. Urban centers and industrial areas are generally high-temperature centers; on the other hand, the river, surface water reservoirs, forests, bigger city parks and wooded cemeteries are generally characterized by much lower temperatures than the average. Land cover type pattern are the consequence of different thermal properties of surfaces (conductivity and capacity); The vegetation cycle causes seasonal changes of thermal properties of particular land cover classes; Land relief and building structure insolated slopes and parts of buildings are warmer whereas shadowed surfaces are colder, creating a positive daytime energy balance of street canyon structure [73].
With the rapid urban expansion, many neighboring cities are socially and economically connected, and the infrastructure networks are interconnected, thereby forming a huge urban agglomeration [49]. The special micro-climate effect forms and introduces higher temperatures in the urban area than its surroundings. Higher levels of energy consumption indicate higher LST intensity [84]. The LST are not only influenced by land cover patterns and process but also significantly dominated by the specific land conversion processes [63].
In the Pearl River Delta Metropolitan Region, the cooling effects of ecological land loss and gain are significantly different, which provides evidence that the existing natural ecosystems (especially forests) are valuable for climatic adaptation because newly constructed ecological land does not provide the same cooling effect, the relationship between urbanization rate and temperature may conform to the Environmental Kuznets Curve [63]. In the Yangtze River Delta Urban Agglomeration, urban areas were significantly positively correlated with LST intensity, while the average wind speed and precipitation was significantly negatively correlated with LST intensity [52]. In the Central Guizhou urban agglomeration, about 96.8% of the regional area NDVI is negatively correlated with LST, and about 3.2% of the regional area is positively correlated [85]. In the Chang-Zhu-Tan urban agglomeration, 80.825% of all the pixels has changed to impervious areas [86].
In addition, the active factory may influence the changed pattern of LST. It was well documented that with the rapid urbanization and population growth, increase in built-up land and replace of natural lands with artificial buildings altered surface energy budgets, and may affect local, regional, and global climate by changing the surface physical properties [87]. LST has been usually regarded as the local or meso-scale climatic phenomenon associated with human activities, there were significant relationships between LST and fractional vegetation cover, population density, and road density [66]. In Shenyang, the built-up area has doubled in 1986-2007, most of which is converted from croplands around the urban fringe area. The built-up area has close association with the population increase, the gross domestic production, and fixed asset investments [74]. During this period, large amounts of impervious surfaces have replaced the semi-natural and natural landscape, especially in the high-rise city proper and newly emerging industrial parks as well as intensive residences along the outer traffic ring [66]. Over the past one decade many residents in CBD removed to the urban fringe and the surrounding satellite towns, owing to the recent urban regeneration. Thus, population density in the CBD significantly decreased and the former residences gave the way to newly emerging green spaces. Old residences in the city center were demolished, this resulted in significant population agglomeration and enhancement of infrastructures in the urban fringe and suburban towns [88]. Therefore, the former boundary between suburban satellite towns and the city center became unclear, and the urban form has shifted from the single-core to the hybrid of multi-nuclei and multiple axis-based sprawl.

Conclusions
This paper has analyzed the spatiotemporal differentiation characteristics of the surface thermal landscape in the Yangtze River Delta for the past 12 years, using MODIS surface temperature data. The research conclusions are as follows: From 2007 to 2019, the LST in the Yangtze River Delta region generally increased (the proportion of middle-, sub-high-and high-temperature zones increased by 33.42%), and the higher temperature areas gradually expanded into their peripheral areas in accordance with the "Z"-shaped spatial pattern (Nanjing-Shanghai-Hangzhou-Ningbo) of urban development in the Yangtze River Delta. The relationship between LST changes and city hierarchy was apparent, in that the node cities around the core cities underwent a large temperature increase (Zhenjiang (71.14%), Jiaxing (53.70%), Changzhou (51.70%), Wuxi (48.44%), Suzhou (Jiangsu) (44.57%), Nantong (41.78%), Huzhou (37.90%), Nanjing (36.36%) and Taizhou (Jiangsu) (33.01%)), the high-calorific value areas of high-level cities were wide and contiguous, and the high-calorific value areas of low-level cities were small and isolated. The corridor effect of thermal changes on the surface was very significant, forming a trend of reciprocal integration and development. Among said corridors, the Shanghai-Nanjing-Hefei-Hangzhou-Ningbo Development Belt and the Yangtze River Development Belt showed obvious warming effects, followed by the Coastal Development Belt and the Shanghai-Hangzhou-Jinhua Development Belt. The changes in the surface thermal landscape index were subject to regional differences, The more significant change is the decrease in the degree of agglomeration in the sub-low-temperature zone (The AI of the sub-low-temperature zone in key development region dropped from 80.35 to 69.63); the degree of agglomeration in the high-temperature zone gradually increases (The AI of the high-temperature zone in the optimized development region increased from 75.94 to 83.68). The patch type of thermal landscape pattern increases, and the distribution of landscape area among various types tends to be even. This phenomenon is most significant in optimized development region (In 2019, the SHDI and SHEI were as high as 1.3015 and 0.8087 respectively).