Landsat-Based Monitoring of the Heat Effects of Urbanization Directions and Types in Hangzhou City from 2000 to 2020

Rapid urbanization has produced serious heat effects worldwide. However, the literature lacks a detailed study on heat effects based on the directions and types of urban expansion. In this work, a typical city with an extremely hot summer climate, Hangzhou, was selected as a case study to determine the relationships between the urban heat-effect dynamics and spatiotemporal patterns of impervious surface expansion. Based on long-term Landsat imagery, this study characterized the spatiotemporal patterns of urban expansion and normalized surface temperatures in Hangzhou City from 2000 to 2020 using object-based backdating classification and a generalized single-channel algorithm with the help of a land-use transfer matrix, expansion index, and spatial centroids. Relevant policies, industries, and traffic networks were discussed to help explain urban expansion and thermal environment changes. The results demonstrated that in 2020, the area of impervious surfaces covered 1139.29 km2. The majority of the gains were in farmland, water, and forests, and the annual growth rate was 32.12 km2/year beginning in 2000. During the expansion of impervious surfaces, the city warmed at a slower rate, and more thermal contributions came from sub-urban areas. The southeast-oriented expansion of impervious surfaces was the key reason for the spatiotemporal dynamics of the urban heat effects. The dominant urban edge expansion intensified the local heat effects. This research provides a Landsat-based methodology for better understanding the heat effects of urban expansion.


Introduction
Urban areas comprise approximately 5% of the total surface of the Earth [1]. However, more than half of the human population lives in cities [2]. This arrangement has produced a series of global environmental problems due to the replacement of natural land cover with impervious surfaces [3]. Meanwhile, urbanization introduces a number of local modifications, such as anthropogenic heat, emissions, roughness and canyons, multiple solar inter-reflections in packed areas, and diversified shading, most of which alter the radiative forcing and air circulation [4,5]. Heat stress due to these local modifications on land-use types, radiative forcing, and air circulation is one of the most severe challenges facing urban residents [6,7]. This phenomenon is traditionally called the urban heat island effect, and induces the growth of energy consumption, especially in summer. This phenomenon has a strong impact on the local air quality and greenhouse gas emissions [8,9]. Therefore, subjects related to the spatiotemporal patterns of urbanization and its heat effects have attracted wide attention and extensive investigations [10][11][12].
Satellite remote sensing data, such as imagery from Landsat, the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), the Advanced Very High Resolution Radiometer (AVHRR), and the Moderate Resolution Imaging Spectroradiometer (MODIS), provide thermal infrared bands and spatially explicit earth observations over long-term periods in a comparable and repeatable manner [13][14][15]. Due to their relatively high spatial resolution and long time-sequenced data, Landsat series have become a common measure for characterizing urbanization and its heat effects [11,16,17]. Lu et al. [18] summarized methods to extract impervious surface areas from satellite images, concluding that object-based approaches are preferable for reducing the heterogeneity of the urban landscape, and that Landsat images are valuable for long-term comparisons. Aina et al. [19] used time-series Landsat images to map impervious surfaces and analyzed the urban landscape through change analysis and determining expansion types. As a direct surface phenomenon, the land-use change from natural land cover to impervious surfaces was a common perspective of remote sensing-based monitoring of the urbanization. Namely, the classification of impervious surfaces from Landsat imagery commonly represents the first step. Moreover, an object-based backdating approach based on existing mapping efforts, with a consistent classification scheme, was able to improve the accuracy by decreasing pseudo changes and heterogeneity [18,20]. Then, urbanization patterns were conventionally characterized via the spatiotemporal analysis of impervious surfaces. This type of analysis involves studying areal changes, spatial distribution, and conversion from the land-use transfer matrix; expansion types from expansion index calculations; and directions based on spatial centroid movements [14,[21][22][23]. After that, the linkage between land-use changes and local thermal patterns was historically determined in one of two ways: (1) a direct temperature comparison between impervious surfaces and other land-use types or a comparison between temperature at different distances from the city; and (2) indirect analysis based on the temporal phenological changes of urban plants, with the former method being more popular [11,24,25]. Du et al. [26] adopted Landsat-based impervious surface mapping and land surface temperature (LST) retrieval and analyzed urban expansion and its heat-island effects using a transfer matrix and buffer analysis, along with a discussion of driving forces. Estoque et al. [27] conducted a correlation analysis between the density of impervious surfaces or green space and LST and compared the mean LSTs of different sizes, shape complexities, and aggregations of patches of impervious surfaces and green space to quantify the urban heat effect based on Landsat images. Roy et al. [17] assessed land-use dynamics based on Landsat images using areal changes and type conversions and linked the spectral indices of urban expansion to LST through linear regressions. The generalized single-channel algorithm is a mature method with wide applicability that is used to accurately map fine-resolution LSTs [26,27]. Overall, previous studies have closely analyzed the relevant classification methods, the spatiotemporal patterns of urban expansion, the surface urban heat island effects, and the driving forces [24,28,29]. However, there remains a lack of studies relating urban expansion types to local heat effects.
Hangzhou, an emerging "stove city" (which implies an extremely hot summer climate) experienced rapid urban development and an increasing urbanization rate, from 22.96% in 1978 to 61.81% in 2015, thereby becoming the second-hottest metropolitan area in China [14,30]. Urban expansion and its heat effects in Hangzhou City have drawn additional attention due to government planning efforts related to polycentric urban structures and market forces. Nevertheless, prior research was limited to separately studying urban expansion patterns and the land/air temperature changes before 2015. However, the precise up-to-date influences of urban expansion patterns on the thermal environment of Hangzhou City require further exploration. Most notably, the detailed heat effects of different urbanization types and directions in Hangzhou have been insufficiently researched but remain crucial for supporting government planning to achieve urban heat mitigation. Hence, in this study, the spatiotemporal patterns of urban expansion and their heat effects were quantified using object-based backdating classification and the generalized singlechannel algorithm based on Landsat imagery between 2000 and 2020. Then, the explicit spatiotemporal effects of urban expansion compared to LST patterns were analyzed with the help of a land-use transfer matrix, an expansion index, and spatial centroids. The specific objectives were to: (1) characterize the spatiotemporal patterns of urban expansion; (2) quantify the spatiotemporal variations in the LST; and (3) analyze the heat effects of urbanization directions and types in Hangzhou City from 2000 to 2020.

The Study Area
Hangzhou is the capital of Zhejiang Province in China and the second largest metropolis on the Yangtze River delta. It is located in the lower reaches of the Qiantang River ( Figure 1). Based on previous studies, the study area (119 • 40 38 -120 • 42 50 E, 30 • 33 54 -29 • 50 48 N), Hangzhou City, was selected from the main metropolitan area, which includes the Shangcheng, Xiacheng, Jianggan, Gongshu, Xihu, Binjiang, Xiaoshan, and Yuhang districts, covering 3345.91 km 2 [14]. This region has a humid subtropical monsoon climate with distinct urban heat island effects in summer [31]. The study area has an annual average temperature of 17.83 • C and annual precipitation of 1941.23 mm [32,33]. The summer daytime and nighttime temperatures of the area can reach nearly 40 and 20 • C, respectively. In the past two decades, characterized by its famous landscape and rapid economic growth, this area has experienced rapid urbanization, with the urban population rising from 3.73 million in 2000 to 8.13 million in 2019 [32,33]. Moreover, this region has experienced intense land-use changes during urban expansion, and consequently suffered from environmental effects [34].
heat effects were quantified using object-based backdating classification and the generalized single-channel algorithm based on Landsat imagery between 2000 and 2020. Then, the explicit spatiotemporal effects of urban expansion compared to LST patterns were analyzed with the help of a land-use transfer matrix, an expansion index, and spatial centroids. The specific objectives were to: (1) characterize the spatiotemporal patterns of urban expansion; (2) quantify the spatiotemporal variations in the LST; and (3) analyze the heat effects of urbanization directions and types in Hangzhou City from 2000 to 2020.

The Study Area
Hangzhou is the capital of Zhejiang Province in China and the second largest metropolis on the Yangtze River delta. It is located in the lower reaches of the Qiantang River ( Figure 1). Based on previous studies, the study area (119°40′38″-120°42′50″E, 30°33′54″-29°50′48″N), Hangzhou City, was selected from the main metropolitan area, which includes the Shangcheng, Xiacheng, Jianggan, Gongshu, Xihu, Binjiang, Xiaoshan, and Yuhang districts, covering 3345.91 km 2 [14]. This region has a humid subtropical monsoon climate with distinct urban heat island effects in summer [31]. The study area has an annual average temperature of 17.83 °C and annual precipitation of 1941.23 mm [32,33]. The summer daytime and nighttime temperatures of the area can reach nearly 40 and 20 °C, respectively. In the past two decades, characterized by its famous landscape and rapid economic growth, this area has experienced rapid urbanization, with the urban population rising from 3.73 million in 2000 to 8.13 million in 2019 [32,33]. Moreover, this region has experienced intense land-use changes during urban expansion, and consequently suffered from environmental effects [34].

Data and Preprocessing
The 16 remote sensing images adopted in this study included Landsat Thematic Mapper (TM) imagery from 2000 and 2010 and Landsat Operational Land Imager (OLI) imagery from 2020 ( Table 1). The level 2 surface reflectance images (atmospherically corrected from level 1 images) with the smallest cloud cover used for LST retrieval in summer (June to August) were downloaded from the United States Geological Survey (USGS) Center (http://glovis.usgs.gov/, accessed on 23 January 2021). These three images were preprocessed and converted from the digital number (DN) values of the thermal bands (band 6 for Landsat-5 and band 10 for Landsat-8) to the at-sensor brightness temperatures [25,35]. For land-use classification, mosaic images from 1 May to 31 October were downloaded from the Google Earth Engine (GEE) cloud-computing platform after pre-processing ( Figure 2). Based on GEE, a series of 16 low-cloud-cover (less than 20%) Landsat surface reflectance tier 1 images were selected and processed to composite median values of the multi-spectral bands at 30 m spatial resolution by remove various types of clouds and noise, with a compositing function [36][37][38]. The acquired mosaic image of the Landsat-8 OLI is shown in Figure 1. GlobeLand30 2020, a 30-m resolution global land-cover data product, was developed by China, with a total accuracy of 85.72% and a Kappa coefficient of 0.82 [39]. The names of the mapsheets of N50_30, N51_25, and N51_30 from GlobeLand30 2020 were applied and downloaded from the National Geomatics Center of China (http://www.globallandcover.com/, accessed on 4 February 2021) as the baselines for the classification.

Methods
As shown in Figure 2, the workflow of this study contained three major parts: (1) extracting spatiotemporal characteristics of impervious surfaces using the object-based backdating classification and the land-use transfer matrix (Sections 2.3.1 and 2.3.4); (2) quantifying the spatiotemporal patterns of the thermal environment through normalized land surface temperature (NLST) retrieval, and calculation of the surface urban heat island intensity and spatial centroids (Sections 2.3.2, 2.3.3 and 2.3.6); and (3) analyzing the heat effects of the urban expansion types and directions by calculating the expansion index, spatial centroids, and NLSTs of different expansion types (Sections 2.3.5 and 2.3.6).

Object-Based Backdating Classification
Based on the classification system of GlobeLand30 products, land-use types in the study area contained farmland, grassland, forests (including forest and shrub-wood), wetlands, water, and impervious surfaces [40]. To achieve change analysis effectively with a consistent classification scheme, while recognizing existing mapping efforts, the object-  Figure 2, the workflow of this study contained three major parts: (1) extracting spatiotemporal characteristics of impervious surfaces using the object-based backdating classification and the land-use transfer matrix (Sections 2.3.1 and 2.3.4); (2) quantifying the spatiotemporal patterns of the thermal environment through normalized land surface temperature (NLST) retrieval, and calculation of the surface urban heat island intensity and spatial centroids (Sections 2.3.2, 2.3.3 and 2.3.6); and (3) analyzing the heat effects of the urban expansion types and directions by calculating the expansion index, spatial centroids, and NLSTs of different expansion types (Sections 2.3.5 and 2.3.6).

Object-Based Backdating Classification
Based on the classification system of GlobeLand30 products, land-use types in the study area contained farmland, grassland, forests (including forest and shrub-wood), wetlands, water, and impervious surfaces [40]. To achieve change analysis effectively Remote Sens. 2021, 13, 4268 6 of 18 with a consistent classification scheme, while recognizing existing mapping efforts, the object-based backdating method was adopted to acquire land-use maps of 2000, 2010, and 2020. As a common time interval for temporal comparisons, a decade was used in this study. The accuracy assessment of land-use classification was based on randomly selected samples from Google Earth images, presenting overall accuracy of over 85% (Table 2). Classification involved the following five steps [20,41]: Object-based image analysis: The GlobeLand30 map of 2020, as Phase I data, and Landsat-8 OLI mosaic image of 2020, as Phase II data, were segmented using the eCognition Developer 9.0 software with 50, 0.1, and 0.5 as values of the scale, shape, and compactness, respectively. 2.
Stratified image classification: The classes of objects from Phase II data were assigned by that from Phase I data.

3.
Visual interpretation and manual modification: The changed objects from the Phase II data were modified manually with the help of finer resolution images from Google Earth.

4.
Dissolve was applied to obtain the final land-use map of 2020 at a 30 m spatial resolution.

5.
Backdating: The final land-use map of 2020 was used as the Phase I data and the Landsat-5 TM mosaic image of 2010 was used as the Phase II data. The steps from 1 to 4 were conducted to acquire a land-use map for 2010. In the same way, based on the land-use map of 2010 (as Phase I data) and the Landsat-5 TM mosaic image of 2000 (as Phase II data), the steps from 1 to 4 were conducted to acquire a land-use map for 2000.

Normalized Land Surface Temperature Retrieval
The LST was retrieved based on the Landsat imagery and the generalized singlechannel algorithm. In detail, the radiative transfer equation was utilized to estimate the surface temperature from the thermal infrared bands of Landsat-5 TM (band 6) and Landsat-8 (band 10). The retrieved LSTs of 2000 and 2010 were consistent with the temperature products from the DATABANK Remote Sensing Data Engine (http://databank.casearth. cn/, accessed on 13 January 2021), the errors of which were below 2 K [42]. Due to the lack of LST products in the summer of 2020 and insufficient weather station data in the study area, the LST in 2020 was not validated in this study.
The LST was retrieved using Equation (1) [42]: where T s is the LST; L i is the radiance of band i; ε i is the land-surface emissivity calculated by the normalized difference vegetation index (NDVI) threshold method [14,25]; Ψ 1 , Ψ 2 , and Ψ 3 are the atmospheric functional parameters calculated by the atmosphere water vapor content (w); w is calculated from the MODIS water-vapor product (MOD05, with a similar time of passing territory with used Landsat images) based on the ratio technique (http://modis-atmos.gsfc.nasa.gov/, accessed on 20 January 2021) [42]; γ and δ are the Planck' parameters, which were calculated by Equations (2) and (3), respectively: where T i is the at-sensor brightness temperature; b γ was 1256 K for Landsat-5 band 6 and 1324 K for Landsat-8 band 10.
To compare spatiotemporal variations, the NLST was adopted using the following equation: where T sj , T smin , and T smax are the LST values of pixel j and the minimum and maximum of all pixels, respectively. Here, larger values of NLST denote relatively high LSTs at spatiotemporal scales. The smaller values of NLST variations between impervious surfaces and other land-use types represent weaker heat effects.

Analysis of Surface Urban Heat Island Intensity
To quantify the spatiotemporal patterns of heat effects from urbanization, the surface urban heat island intensity (SUHII) was calculated. In this study, the SUHII was defined as the NLST difference between areas of high and sub-high intensity impervious surfaces, which is different from other studies that used LST as the base [8,43]. Our SUHII definition can reduce the uncertainties associated with the land surface temperature retrieval and site-specific rural conditions. The SUHII was acquired through following steps:

1.
Impervious surface intensity maps were generated from land use data using moving windows of 120 × 120 m for 2000 and 2010 and 100 × 100 m for 2020 to match the pixel sizes of Landsat-5 band 6 and Landsat-8 band 10, respectively.

2.
A 50% threshold of the impervious surface intensity was used as a criterion to separate the maps into high-and low-intensity impervious surfaces [44].

3.
The buffer areas outside the high intensity impervious surfaces with distances of 240 m for 2000 and 2010 and 200 m for 2020 were used as the sub-high areas. Finally, the NLST difference between the high-intensity area and its surrounding sub-high areas was taken as the SUHII.

Land-Use Transfer Matrix
The land-use transfer matrix is widely adopted to characterize the structure of landuse dynamics. It describes areal conversions across different land-use types in a certain region during a particular period [45]. To quantify the contributions of other land-use types to the growth of impervious surfaces in Hangzhou City, the land-use transfer matrix was calculated using the tabulate area tool in the ArcGIS software: where S ij is the area of land-use type i transferred to land-use type j, and n is the number of land-use types. In this study, n was 6. i and j (1, 2, . . . , 6) are land-use types before and after a certain transfer process, respectively. Each element on the main diagonal of the matrix denotes the area of each land-use type remaining unchanged. The contributions of other land-use types to the growth of impervious surfaces represent the values from one of the columns of this matrix.

The Expansion Index
The expansion index (E) was chosen to determine three main urban expansion types: infilling, edge-expansion, and outlying [46,47]. A newly-grown patch expanding as a gap or hole among or within old patches is defined here as an infilling expansion, whereas uniaxial spreading in parallel strips from one edge is an edge-expansion type. The expansion is Remote Sens. 2021, 13, 4268 8 of 18 referred to as outlying expansion when a newly-grown patch is found isolated from an old patch [48]. E was acquired by using Equation (6) via the analysis tool in the ArcGIS software: where P new is the perimeter of a newly-developed patch of impervious surfaces, L com is the length of the common edge of the expansion patch of impervious surfaces, and L com and P new have the same units (such as kilometers). In this method, the patch is identified as infilling when E > 0.5, edge-expansion when 0 < E ≤ 0.5, and outlying when E = 0.

Spatial Centroids
The spatial centroid is defined as the location determined by the coordinates of the weighted geometric center of a polygon or multiple polygons, as applied in the analysis of land-use type evolution and thermal environmental patterns [49,50]. The spatial centroids of impervious surfaces and the surface temperatures in 2000, 2010, and 2020 were mapped to delineate the direction of changes as follows: where X t and Y t are the longitude and latitude of a centroid of impervious surfaces and the surface temperature in year t, respectively; C ti is the area of patch i in year t when calculating the centroid of impervious surfaces; C ti is is the NLST values of the patch i in year t for the calculation of centroids of the temperature; X i and Y i are the longitude and latitude of patch i used for impervious surfaces or NLST, respectively; and N is the total patch number of impervious surfaces or NLST. When impervious surfaces or the surface temperature have the same growth or decline rates in every direction, the spatial centroid remains invariant; otherwise, its centroid moves toward the direction in which impervious surfaces (or the surface temperature) gain or reduce more [51].

Spatiotemporal Distribution of Impervious Surfaces
Impervious surfaces were spread widely across Hangzhou City from 2000 to 2020 and were mainly located to the north of the Qiantang River ( Figure 3). In 2020, the area of impervious surfaces totaled 1139.29 km 2 , which represented 34.05% of the total study area ( Figure 4). Xiacheng, Gongshu, Jianggan, Shangcheng, and Binjiang were the top five districts with the greatest percentages of impervious surfaces, with values of 96.32%, 77.64%, 69.65%, 69.09%, and 67.77%, respectively (Figure 4). In the last 20 years, although expansion slowed in Shangcheng, Xiacheng, Jianggan, Gongshu and Binjiang districts, impervious surfaces expanded with an increasing growth rate of 6.11% in Hangzhou City (Figure 4). From 2000 to 2010, the expansion rate of impervious surfaces in Hangzhou City was 31.16 km 2 /year, with Xiaoshan, Yuhang and Jianggan being the top three districts. In the following 10 years, this rate increased to 33.07 km 2 /year. Xiaoshan, Yuhang, and Xihu districts had the highest growth speeds of impervious surfaces (Figure 4).
Over the 20 years under analysis, impervious surfaces mostly occupied farmland and water, which reached above 90% of the contribution proportion for the growth of impervious surfaces ( Figure 5). During the two phases, the contributions of farmland and water notably decreased, while that of grassland reduced slightly and that of forests increased by 3.15%. Most reductions in farmland and water due to urban expansion occurred in Jianggan, Binjiang, Xiaoshan, and Yuhang districts (Figure 3). Impervious surfaces occupied forests and grassland principally in the Xiaoshan and Yuhang districts.  Over the 20 years under analysis, impervious surfaces mostly occupied farmland and water, which reached above 90% of the contribution proportion for the growth of impervious surfaces ( Figure 5). During the two phases, the contributions of farmland and water notably decreased, while that of grassland reduced slightly and that of forests increased by 3.15%. Most reductions in farmland and water due to urban expansion occurred in Jianggan, Binjiang, Xiaoshan, and Yuhang districts (Figure 3). Impervious surfaces occupied forests and grassland principally in the Xiaoshan and Yuhang districts.   Over the 20 years under analysis, impervious surfaces mostly occupied farmland and water, which reached above 90% of the contribution proportion for the growth of impervious surfaces ( Figure 5). During the two phases, the contributions of farmland and water notably decreased, while that of grassland reduced slightly and that of forests increased by 3.15%. Most reductions in farmland and water due to urban expansion occurred in Jianggan, Binjiang, Xiaoshan, and Yuhang districts (Figure 3). Impervious surfaces occupied forests and grassland principally in the Xiaoshan and Yuhang districts.

Spatiotemporal Dynamics of LST
The NLST showed dramatic spatiotemporal variations. The increasing distribution of high values indicated the warming environment of Hangzhou City ( Figure 6). Overall,

Spatiotemporal Dynamics of LST
The NLST showed dramatic spatiotemporal variations. The increasing distribution of high values indicated the warming environment of Hangzhou City ( Figure 6). Overall, irrespective of the outlying smallest values of NLST in 2020 due to the cloud cover, especially in the northwestern and southern forests, spatial variations in surface temperature decreased during the two phases of 2000-2010 and 2010-2020 (Figures 6 and 7). Impervious surfaces had the largest mean values of NLST, followed by farmland (Figure 7). However, the smallest mean values of NLST were observed for the wetlands in 2000 and for forests in 2010 and 2020.

Spatiotemporal Dynamics of LST
The NLST showed dramatic spatiotemporal variations. The increasing distribution of high values indicated the warming environment of Hangzhou City ( Figure 6). Overall, irrespective of the outlying smallest values of NLST in 2020 due to the cloud cover, especially in the northwestern and southern forests, spatial variations in surface temperature decreased during the two phases of 2000-2010 and 2010-2020 (Figures 6 and 7). Impervious surfaces had the largest mean values of NLST, followed by farmland (Figure 7). However, the smallest mean values of NLST were observed for the wetlands in 2000 and for forests in 2010 and 2020.  The SUHII decreased in the last 20 years (Table 3). This result was consistent with the afore-mentioned spatiotemporal distribution of NLST, through which spatial variations in the heat effects of urban expansion reduced from 2000 to 2020. However, the NLST

Spatiotemporal Dynamics of LST
The NLST showed dramatic spatiotemporal variations. The increasing distribution of high values indicated the warming environment of Hangzhou City ( Figure 6). Overall, irrespective of the outlying smallest values of NLST in 2020 due to the cloud cover, especially in the northwestern and southern forests, spatial variations in surface temperature decreased during the two phases of 2000-2010 and 2010-2020 (Figures 6 and 7). Impervious surfaces had the largest mean values of NLST, followed by farmland (Figure 7). However, the smallest mean values of NLST were observed for the wetlands in 2000 and for forests in 2010 and 2020.  The SUHII decreased in the last 20 years (Table 3). This result was consistent with the afore-mentioned spatiotemporal distribution of NLST, through which spatial variations in the heat effects of urban expansion reduced from 2000 to 2020. However, the NLST The SUHII decreased in the last 20 years (Table 3). This result was consistent with the afore-mentioned spatiotemporal distribution of NLST, through which spatial variations in the heat effects of urban expansion reduced from 2000 to 2020. However, the NLST changes yielding impervious surfaces with different intensities varied over time. From 2000 to 2010, the high and sub-high intensity impervious surfaces featured a decreasing NLST. Nonetheless, in the next 10 years, the NLST of high and sub-high intensity areas grew. This result indicated that the spatial variations in heat effects in high and sub-high intensity urban areas weakened and then increased, but were still weaker than the initial values. The centroid of NLST moved to the southeast from 2000 to 2020, and was mainly located in the north of the Qiantang River, which was in line with the major distribution of impervious surfaces (Figures 3 and 8). Movements of the NLST centroid revealed that the heat effects of urban expansion in the southeastern part of Hangzhou city experienced a larger change rate than the effects in other directions. The centroid movement was more strongly distinguished during the first 10 years, similarly to the changes in SUHII. In other words, the change in the heat effects of urban expansion was faster from 2000 to 2010 than that during the last 10 years. intensity urban areas weakened and then increased, but were still weaker than the initial values. The centroid of NLST moved to the southeast from 2000 to 2020, and was mainly located in the north of the Qiantang River, which was in line with the major distribution of impervious surfaces (Figures 3 and 8). Movements of the NLST centroid revealed that the heat effects of urban expansion in the southeastern part of Hangzhou city experienced a larger change rate than the effects in other directions. The centroid movement was more strongly distinguished during the first 10 years, similarly to the changes in SUHII. In other words, the change in the heat effects of urban expansion was faster from 2000 to 2010 than that during the last 10 years.

Directions, Types, and Temperature Variations of Urban Expansion
The centroid of impervious surfaces moved to the east from 2000 to 2010, and then turned to the south in the next decade ( Figure 8). This result reveals that in the first decade, urban expansion in the eastern part of Hangzhou City was more notable than that in other directions, and the growth pole transformed to the south. This result agrees with NLST temporal dynamics showing that the spatial variations of urban expansion in Hangzhou

Directions, Types, and Temperature Variations of Urban Expansion
The centroid of impervious surfaces moved to the east from 2000 to 2010, and then turned to the south in the next decade ( Figure 8). This result reveals that in the first decade, urban expansion in the eastern part of Hangzhou City was more notable than that in other directions, and the growth pole transformed to the south. This result agrees with NLST temporal dynamics showing that the spatial variations of urban expansion in Hangzhou City were more remarkable in the first 10 years according to the smaller movement of centroids from 2010 to 2020.
The distribution of different expansion types illustrates that the edge-expansion was the uppermost type and widely distributed across Hangzhou City in the last two decades (Figure 9). The outlying expansion was mainly clustered in Jianggan, Xihu, Xiaoshan and Yuhang districts. Infilling expansion was the least common and scattered in Hangzhou City, with the majority located in the east. The edge-expansion increased from 2000 to 2020, while outlying and infilling decreased. Firstly, edge-expansion and the outlying urban patches featured relatively close proportions to those of the former aggregations in Jianggan and Binjiang districts, with the latter typically located in Xiaoshan District. Nonetheless, the edge-expansion became dominant from 2010 to 2020, especially in Xiaoshan and Yuhang districts. At the same time, outlying urban patches were primarily located in Xihu, Xiaoshan, and Yuhang districts. The heat effects of the three types of urban expansion displayed different variations ( Table 4). The edge-expansion had the largest values of maximum, average, and variations of NLST during the last 20 years. Conversely, infilling urban expansion possessed the smallest values of the maximum, average, and variation of NLST. However, the minimum value of NLST belonged to the outlying expansion from 2000 to 2010, which appeared in the edge-expanded urban patch during the next 10 years. Generally, the heat effects of infilling growth were slightly weaker than those of outlying growth, and those of the edge-expansion were strongest. City, with the majority located in the east. The edge-expansion increased from 2000 2020, while outlying and infilling decreased. Firstly, edge-expansion and the outlying u ban patches featured relatively close proportions to those of the former aggregations Jianggan and Binjiang districts, with the latter typically located in Xiaoshan Distri Nonetheless, the edge-expansion became dominant from 2010 to 2020, especially Xiaoshan and Yuhang districts. At the same time, outlying urban patches were primari located in Xihu, Xiaoshan, and Yuhang districts. The heat effects of the three types of u ban expansion displayed different variations ( Table 4). The edge-expansion had the lar est values of maximum, average, and variations of NLST during the last 20 years. Co versely, infilling urban expansion possessed the smallest values of the maximum, averag and variation of NLST. However, the minimum value of NLST belonged to the outlyin expansion from 2000 to 2010, which appeared in the edge-expanded urban patch durin the next 10 years. Generally, the heat effects of infilling growth were slightly weaker tha those of outlying growth, and those of the edge-expansion were strongest.

Spatiotemporal Patterns of Urban Expansion
We observed that the expansion of Hangzhou City was mainly concentrated in the northeast of the city with flat terrain (Figures 3 and 8). This concentration was due to the limitations of the mountainous regions in the south and west, such as Tianmu and Chuanwu Mountains [52]. Given the constraints of land availability, Xiaoshan and Yuhang districts were merged into Hangzhou City in 2001, which provided these two administrative units the highest urban expansion rates (Figure 4). Further, Xiaoshan and Yuhang districts featured expanses of available farmland, which made farmland the greatest contributor to urban sprawl (Figures 3 and 5). The centroid movement revealed that urban expansion was more rapid in the southeast (Figure 8), which resulted from the extensive urban growth in the flat terrain of southeastern Hangzhou City, such as in Xiaoshan District and the east of Yuhang District. This phenomenon also agreed with the urban planning emphasized in the Hangzhou Ninth Party Congress involving "active east, quiet west, new south, beautiful north, and prosperous center" and "eastward urban expansion, westward tourism and industry, and development along the two sides of the Qiantang River" [53].
The results illustrated that edge-expansion was the dominant type of urban growth from 2000 to 2020, while infilling expansion was sparsely distributed (Figure 9). This result agrees with prior studies on urban expansion from 2000 to 2015 and can be explained by strip development with subjection to the southwestern mountains and the Qiantang River [54]. A previous study concluded that outlying growth, followed by edge-expansion, was the primary type of urban sprawl from 1995 to 2005 in districts that were built earlier [55]. It can be inferred that the urban expansion of Hangzhou City transformed from outlying-oriented to leading edge-expansion (Figure 9a). Infilling expansion was lowest because of the high land use efficiency in Hangzhou City, which featured few gaps or holes across urban areas.

Spatial Characteristics of Heat Effects from the Urban Expansion
Unlike previous studies, the present study focused on the relative values of LST across Hangzhou City, and their spatiotemporal dynamics based on NLST. It was reported that the absolute values of the air and land temperature in urban areas increased during the urbanization process in the study area, which was shown in this study as high increases in NLST ( Figure 6) [14,30,31,56]. It was found that the spatial variations in the heat effects of urban expansion decreased during the last two decades, which was principally due to the growing distribution of high values of NLST ( Figure 6 and Table 3). In other words, the whole city has heated during the process of urban expansion. Although the whole city became warmer, the cooling effects of forests and wetlands were more distinguished from 2000 to 2020, as shown in Figure 7, which was explained in past studies as being caused by carbon sequestration and the climate-regulation functions of urban forests and wetlands [25,57,58].
The SUHII values and NLST centroid movements revealed that the spatial variations in the heat effects of urban expansion were significant in the eastern part of Hangzhou City from 2000 to 2010 (Figure 8 and Table 3). Specifically, the sprawl of impervious surfaces in the east was a key factor for the decrease in the spatial variation of urban heat effects. However, in the last 10 years, despite the NLST growth in high and sub-high intensity impervious surfaces, spatial variations in urban heat effects have decreased southwards at a slower rate. This observation indicates that changes in heat effects were larger in sub-urban than core-urban areas and that the south-oriented expansion of impervious surfaces largely induced warming of the city. Meanwhile, the warming rate slowed, and sub-urban areas contributed more from 2010 to 2020.
The increase in the edge-expansion of impervious surfaces also intensified urban heat effects (Figure 9). The new edge-built urban areas had the warmest relative temperatures (Table 4) because edge-expansion was extensive across Hangzhou City. Meanwhile, edge-expansion conventionally occurred in areas where human activities were intense accompanied by economies of scale and the movement of centers, which hampered the air and radiative circulation [55,59]. These factors simultaneously transformed the urban heat effect from a monocentric to polycentric pattern [14]. These results demonstrated that edge-built impervious surfaces were the hot spots of Hangzhou City with abundant conversions from cooler areas, such as farmland, water, and forests. The declining spatial variations in urban heat effects resulted from decreasing NLST differences among impervious surfaces and these cooler land-use types. The increasing of edge-expansion and decreasing of infilling also verified the growth of contributions from sub-urban areas to urban heat effects in the last decade.

Spatial Characteristics of Heat Effects from Urban Expansion
The driving forces of urban expansion patterns in Hangzhou City included topographic factors, urban planning policies and other socioeconomic determinants, while the influence of natural drivers weakened gradually [52,60,61]. Topographic factors such as the elevation, slope, and natural barriers of mountains and rivers directly resulted in urban expansion patterns as reported in 4.1. Based on the urban master plan of Hangzhou (2001-2020), the urban expansion pattern of this city was designated as a polycentric structure with a main center, three sub-centers, and six clusters [62]. Moreover, this plan emphasized urban expansion along the Qiangtang River and the six ecological belts connecting suburban and exurban ecological, green, and agricultural areas [63]. The theme of urban expansion thus changed from the "Xihu (West Lake) Era" to the "Qiantang River Era" [64]. This urban planning promoted urban expansion across the city and mitigated urban heat effects, with the NLST values of impervious surfaces reduced from 2000 to 2020 (Figures 3 and 7). Although the NLST values of impervious surfaces slightly increased from 2000 to 2010, new polices were put forward with consideration of the landscape composition and configuration and took effects from 2010 to 2020. In this way, starting in 2011, the municipal government took measures to propel industrial transformation and shut down outdated production facilities [65]. Moreover, the metropolis developed a strategy to enhance the optimization of inner-city functions and structures to improve land-use efficiency and restore lost green spaces, resulting in a relatively lower warming rate and greater contributions from sub-urban areas in recent years [66].
Industrial suburbanization was articulated in previous studies as a vital socioeconomic factor for underlying the spatiotemporal dynamics of urban expansion from 2000 to 2015 [62,64,67]. The development core of industries moved from Jianggan District to Xiaoshan District, which was consistent with the movements of the NLST centroid ( Figure 8). This agreement was because industries increased the anthropogenic heat and emissions in southeastern areas of Hangzhou City. In recent years, the development of technologyand knowledge-intensive industries, and science and technology based industrial zones located in Binjiang, Xiaoshan, and Yuhang districts, has entailed increasing labor, land, and production value [32,68,69], thus supporting polycentric urban development and thermal environment and decreasing the spatial variations of urban heat effects (Figures 3 and 6, and Table 3). This development of technology-and knowledge-intensive industries also accelerated roughness changes by the construction of tall buildings, which modified radiative forcing and air circulation by leading to multiple solar inter-reflections in packed areas and diversified shading. Traffic network construction was reported to be another socioeconomic determinant [52,61]. Two subway lines completed in 2012 facilitated the urban sprawl into urban-rural frontiers with a massive physical edge-expansion along major road corridors toward the east and south, which also promoted a boom in urban heat effects on the south shore of the Qiantang River (Figures 3 and 6). Population and economic conditions, as additional common socioeconomic factors that are easily quantified, were considerably influenced by policies in the short term and indirectly influenced urban expansion [55,70].
Thus, the important roles of policies, industries, and traffic networks in the urban expansion of Hangzhou City and its thermal environment were widely recognized in previous studies and further highlighted in this study.

Conclusions
Urban land-use changes hinder ecosystem sustainability through the loss of habitats and carbon storage and an increase in warming the environment. It is urgent further explore the heat effects of different urbanization types and directions in Hangzhou to support urban heat mitigation strategies.
The present research selected Hangzhou City as a case study to determine the spatiotemporal patterns of urban expansion and their heat effects using Landsat imagery between 2000 and 2020. Policies, industries, and traffic networks were discussed to help explain the urban expansion of Hangzhou City and the changes in its thermal environment. The following conclusions were drawn from this study: (1) The impervious surface expansion of Hangzhou City was mainly concentrated in the northeast of the city with flat terrain. In 2020, the area of impervious surfaces totaled 1139.29 km 2 , with the majority being in farmland, water, and forests. The annual growth rate was 31.16 km 2 /year from 2000 to 2010 and 33.07 km 2 /year in the last 10 years. (2) The city has heated over time with decreases in spatial variation during the urbanization process. Impervious surfaces were cores of heat, and forests and wetlands had the lowest surface temperatures. Urban areas expanded with decreases in the warming rate and featured greater contributions from sub-urban areas.