Examining Land Use/Land Cover Change and the Summertime Surface Urban Heat Island Effect in Fast-Growing Greater Hefei, China: Implications for Sustainable Land Development

In this study, a retrospective analysis of the relationship between the land use/land cover (LULC) change and associated surface urban heat island (SUHI) effect in fast-growing Greater Hefei between 1995 and 2016 was performed. Our results reveal the heterogeneous patterns of LULC change. The concentric buffer-based urban–rural gradient analysis reveals that most of the newly emerging developed land occurred within downtown Hefei. In contrast, in three suburban municipality/county jurisdictions, the overall area change in the non-developed land was much lower, but the net increase in developed land is remarkable. Meanwhile, the spatiotemporal patterns of SUHI are in good agreement with that of the developed land, as evidenced by the notable increase in SUHI intensity (SUHII) levels and SUHI spatial extent (SUHISE) in response to the rapid urban expansion, particularly along transportation corridors. In addition, partial least square regression (PLSR) models indicate that the buffer-based predictors/independent variables are significantly related to the responses (SUHII and SUHISE), explaining approximately 61.3% of the variance in the SUHII and 79.8% of the variance in the SUHISE, respectively. Furthermore, the relative strength of the independent variables in determining the relationship was quantitatively examined. The findings of this study provide clear evidence for decision making for sustainable land development and mitigation of the SUHI effect.


Introduction
During the past two centuries, urban areas have been the dominant human settlements globally, and approximately 68% of urban dwellers will live in the cities in 2050 [1]. However, among the consequences due to artificially altered hydrological circulation, climate regulation, soil conservation, and biodiversity preservation [2][3][4][5][6][7][8][9], the urban heat island (UHI) effect, which is associated with human activities and climate change [10][11][12], has attracted increasing concern. Urban areas, which are subject to a local warming effect, or combined local and global warming effects, are the fundamental units for climate change adaptation and mitigation [13]. Thus, how cities cope with the UHI effect, particularly optimizing land use patterns towards sustaining regional ecosystems, will significantly affect UHI mitigation and adaption to climate change [14][15][16][17][18][19][20].
Since the late 1970s, China has been the largest urbanizing country in the world. The UHI effect has been one of the pronounced environmental problems, fueling growing concern among the public. Although the occurrence of the UHI effect is not necessarily proportional to urban size [21], numerous studies on the UHI effect have been focused on the leading megacities along the eastern coastal zones, which accommodate millions of migrants and have experienced dramatic urban growth. By comparison, case studies on the less-developed western and central cities are relatively scarce [22]. Recently, with the implementation of China's national strategies towards reducing the unbalanced development between the developed and less-developed regions [23,24], for the fast-growing big cities in the less-developed western and central regions, the linkage between their booming economic activities and the pronounced UHI effect has attracted more attention [25][26][27].
Two types of UHI effects have been investigated in previous research: the canopy UHI (CUHI) effect, which measures the rural-urban air temperature difference from in situ weather stations; and the surface UHI (SUHI) effect, which measures the rural-urban land surface temperature (LST) difference. Monitoring of the UHI effect has been undertaken with data retrieved from satellite remotely sensed thermal infrared (TIR) imagery [28,29]. The former approach has the advantage of a relatively long time-series of recorded data and high temporal resolution. Unfortunately, given the complicated urban morphology and changeable microclimatic factors, the sparse distribution of weather stations with low spatial coverage cannot sufficiently represent the CUHI effect of the whole city. Thus, this approach can be useful for retrospective analysis and prediction of the fine-scale CUHI effect but is biased at the city level. Alternatively, the satellite TIR imagery, which can capture the TIR signals of land surface with broad spatial coverage, has been widely used to examine the impacts of land use/land cover (LULC) change on the SUHI effect at city and regional levels [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44]. The existing studies mostly emphasize the relationships between the LULC change and SUHI effect, using qualitative correlation analysis or quantitative methods (e.g., bivariate linear regression and multiple regression). Undoubtedly, the quantitative methods can provide more persuasive and convincing results. However, some questions arise from the application of the quantitative methods and the practical demands of decision making. Firstly, because most of the popular quantitative methods are based on the ordinary least-square (OLS) model, it remains a challenge to reveal the quantitative relationships between LULC change and the SUHI effect, considering these simple linear regression models may suffer from multicollinearity between the variables . Under such circumstance, stepwise regression cannot be a cure-all because removing the necessary independent variables may cause the loss of statistical power. Secondly, for a given regression model that contains multiple variables measured with different units, the issue of quantifying the relative importance of each independent variable and its influence on the SUHI effect has yet not been satisfactorily addressed. Thus, these debates must be clarified before the quantitative LULC change versus SUHI effect relationship can be adopted by decision makers engaged in sustainable urban land development and climatic adaption.
In this study, Greater Hefei, one of the star cities in central China and known as one of China's ten representative cities suffering from a pronounced summertime UHI effect [45], is suitable to be used to quantitatively examine the relationship between LULC change and summertime SUHI effect. The objective of the study is to achieve a better understanding of the relationship between LULC change and SUHI effect, and further provide valuable experience for urban planners and decision makers who are engaged in mitigating the SUHI effect and climate adaption for sustainable land development.

Study Region
Greater Hefei, the capital city of Anhui province, is situated in central China between latitudes 30°56′N-32°32′N and longitudes 116°40′E-117°58′E ( Figure 1). This region has a northern subtropical humid monsoon climate, featuring mild temperatures in spring and autumn; a hot summer; and a cold winter. Its long-term recorded extreme temperatures are a maximum of 37.8~39.1 °C during the summer (from late June to early September) and a minimum of −9.9~−13.5 °C during the winter (from December to the following February), with an average annual temperature of 15.7 °C. Its precipitation averages 1000 mm/year, most of which occurs during May and August. Topographically, the city is located on the hilly terrain between the Yangtze River and the Huaihe River. The elevation ranges between 12 and 45 m above sea level.
Hefei is an important city for its geographical advantage in central China, with a recorded history of 2200 years. It was one of the three cities that were first listed as China's National Garden Cities due to its natural landscape and splendid cultural heritage [46,47]. In 2011, the Greater Hefei region reshaped its boundary by merging with the former Chaohu city. Its administrative area almost doubled and the population increased by nearly 2 million. Its present administrative boundary consists of four districts (Luyang, Shushan, Baohe, and Yaohai), one municipality (Chaohu), and four counties (Feidong, Feixi, Changfeng, and Lujiang). It covers an area of approximately 11,430 km 2 , with a total population of 7,965,000 [48]. Over the past decade, this region experienced an unprecedented rural-urban transition. With the implementation of industrial upgradation and transprovincial relocation, intensive industrial parks and new settlements were developed at the periphery. However, the issue of coordinating land development plans, population growth, and mitigation of the UHI effect is absent from the official urban master planning documentation, given the extreme summertime heat events.

Data Sources
This study relies on the Landsat series images for retrieving LST data that indicate the summertime SUHI effect in Greater Hefei. Due to the limitation of weather conditions, the available Landsat series images with high quality and evenly annual intervals for this study are scarce. Therefore, multitemporal Landsat serial images (Level-1T, path 121/row 38) acquired during the summer period, including three Landsat Thematic Mapper (TM) images dated August 1, 1995, July 11, 2002, and August 2, 2007, and one Landsat Operational Land Imager (OLI)/Thermal Infrared Sensor (TIRS) image dated July 25, 2016, were used as the major dataset. All of the Landsat series images were free of cloud/haze contamination and released with terrain correction [49]. The auxiliary datasets include the commercial 30 m LULC maps (required on demand in 1995, 2002, 2007, and 2016) purchased from Beijing Digital Space Technology Co. Ltd, local socioeconomic data, and urban planning documents [46,47,50,51].

Methods
To clearly present the essential procedures for this study, an overall working flowchart, including data processing, LULC classification, LULC change measurement, retrieval of LST, and statistical analysis, is shown in Figure 2.

Processing of Satellite Imagery
To obtain the overall LULC information for the study area, bands 5, 4, and 3 of Landsat TM were combined to generate false-color images. For Landsat OLI/TIRS image, bands 6, 5, and 4 were combined. Based on prior knowledge of Hefei urban planning materials and our field surveys, according to the guidance for national LULC classification [52], the LULC categories of this study region, namely developed land, forest, cropland, water body, and bare land, were adopted (Table 1). By overlapping commercial 30 m LULC maps (1995-2016), for each LULC type, at least 50 featured training sites were selected according to the random stratified sampling method. Thereafter, the supervised maximum likelihood classification method was used to perform LULC classification [53]. The 30 m commercial LULC maps acquired in 1995-2016 and field survey data in 2012-2016 were used for validating the classification results. To assess the accuracy of the classified results, for each image 250 samples were randomly selected. The overall accuracy of LULC maps over the study

LULC Category Description
Developed land Urban and rural settlements, mainly including residences, commercial centers, economic development zone, industrial zones, college towns, railways, and highways Forest Mainly including natural and artificial woodlands, cropland shelterbelts, forest nurseries, and scrublands Cropland Paddy fields, fallow lands after harvest, drylands, and orchards Water body Lakes, rivers, ponds, reservoirs, fishponds, dikes, permanent and seasonal wetlands Bare land Bare rocks, quarries, mines, and vacant lands for urban development.

LULC Change Measurement
In this study, detection of LULC change was performed via overlapping time series of LULC maps and calculating the changed/unchanged pixels falling within each LULC category. To measure the LULC change under the pressure of rapid urbanization, particularly that attributable to human activities and socio-economic development along the rural-urban gradient, a series of concentric zonal buffers with different radii from the city center were drawn, based on our prior knowledge of the rural-urban transition and population growth of the study area. By analyzing the historical trends in population growth and urban encroachment in Greater Hefei, zonal buffers were drawn at 0-1.  (Table 2). Table 2. Description of concentric buffers drawn along the urban-rural gradient.

0-1.5
The city core of downtown Hefei.

1.5-3.0
The urban area between the inner and the first ring roads.
3.0-6.0 The urban area between the first and second ring roads.

Retrieval of LST
In this study, the image-based method, which assumes uniform atmospheric condition and does not require thermal correction of the TIR band [55,56], was employed for retrieval of LST, given all the satellite images were free of cloud contamination. The detailed procedures for retrieval of LST included the following steps. Firstly, for Landsat TM images, prior to the retrieval of LST, a quadratic model was used to convert the digital number (DN) of Landsat TM band 6 into brightness temperatures [57]: For Landsat OLI TIR bands, because band 11 is subject to high uncertainty due to telescope stray light disturbance [58,59], band 10 was used to generated radiant temperatures as follows: where Lλ is the radiance of thermal band pixels in W/(m 2 ster μm), gain is the slope of the radiance/DN conversion function, and offset is the slope of the radiance/DN conversion. Secondly, for Landsat TM and Landsat 8 bands, the spectral radiance was then converted to atsatellite brightness temperature under the assumption of uniform emissivity [60].
where TB is the effective at-satellite temperature in Kelvin (K), and both K1 and K2 are pre-launch calibration constants. Thirdly, since the values of TB were referenced to a black body rather than real earth objects, correction of spectral emissivity must be performed to retrieve the LST of true ground. It is currently a challenge to measure the emissivity values of the LULC types with remarkable variations. To reflect the spatially heterogeneous land surface emissivity, the emissivity of surface materials (ε) was calculated as follows [61]: where NDVI denotes the normalized difference vegetation index, and ρNIR and ρRed denote the reflectance of the near-infrared band and the red band, respectively. Pv denotes the vegetation proportion, NDVImin denotes the minimum value of NDVI, and NDVImax denotes the maximum value of NDVI [62]. v  denotes vegetation emissivity and s  denotes soil emissivity. The values of m and n were referenced to the results of Sobrino et al. [61].

Measuring Summertime UHI Effect Indicators
In rapidly urbanizing areas, measurement of the routine UHI intensity based on the air temperature difference between urban and rural areas is problematic, largely due to the dramatic rural-urban transformation that sharply altered the background air temperature [64,65]. Similarly, SUHI intensity measured by the LST difference between the urban and rural areas is also problematic. Alternatively, inspired by the concept of a local climate zone [65], measurement of the LST difference between the cooling surface (vegetated land and water bodies) and the impervious surface can better indicate the SUHI intensity and spatial extent influenced by the summertime SUHI effect. Therefore, based on the retrieved LST, two important summertime SUHI effect indicators, namely the SUHI intensity (SUHII) and spatial extent influenced by the SUHI effect (SUHISE), were measured as follows: where SUHII is indicated with ℃, LSTIS is the mean LST of the impervious surface within a given zonal buffer, and LSTCS is the mean LST of the cooling surface within a given zonal buffer.
where SUHISE is indicated with km 2 , P is the total number of pixels (30 m resolution) with SUHII ≥ 0, and A is the pixel-specific area. Subsequently, based on our overall understanding of the spatiotemporal pattern of SUHI indicators of the study region, the SUHII was divided into six levels according to the thresholds of percentiles pinpointed on the normal distribution curve, that is, level 1: extreme low (<=5%); level 2: slightly low (5~25%); level 3: low to slightly high (25~50%); level 4: mediate high (50~75%); level 5: very high (75~95%); and level 6: extreme high(>=95%). Similarly, the SUHISE was divided into six levels according to the thresholds of percentiles pinpointed on the distribution curve.

Statistical Analysis
Procedures for data analysis, including descriptive statistics, normality test, and Pearson's product-moment correlation, were employed before quantifying the assumed relationship between the SUHI effect indicators and independent variables. The multicollinearity between the variables was detected according to the result of correlation analysis (see Appendix Table 1). Therefore, to address the problem of multicollinearity that may cause misleading results, a partial least squares regression (PLSR) model, which combines features from and generalizes principal component analysis (PCA) and multiple linear regression [66], was employed to quantify the assumed relationship between the SUHI effect indicators and independent variables. Furthermore, to avoid the over-fitting problem and determine a reasonable model with the appropriate number of components that has good predictive ability, the leave-one-out (LOO) method was used for crossvalidation, according to the criteria for selecting the models with the highest average predicted R 2 and the lowest average prediction sum of squares (PRESS) [66]. The validated PLSR models are written as follows: where α1 is the intercept; β1~β6 are the partial coefficients; and X1~X6 are the area proportions of Developed land, Forest, Waterbody, Cropland, and Bare land within a given buffer zone, respectively. X6 is the mean population density (Pop_density, unit: 1,000 persons per km 2 ) within a given zonal buffer; ε1 is the error term.
where α2 is the intercept, β'1~ β'6 are the partial coefficients; and the definitions of X1~X6 are the same as for Equation (12). ε2 is the error term.
The statistical processes employed in this study were performed with R 3.6 [67] and the pls library [68]. In contrast, water bodies show a relatively slow but stable changing rate (7.82%), increasing from 1082.01 to 1166.65 km 2 during the study period. At the regional level, the dramatic increases in developed land and forest, and slight increases in bare land and water bodies, accompanied by a sharp decrease in cropland, should be of particular concern.

Synoptic Analysis of LULC Change at the Regional Level
In addition, the heterogeneity of LULC change between the administrative zones of Greater Hefei is notable. As can be seen, most of the newly emerging developed land occurred within the extent of downtown Hefei. With the exception of the nearly unchanged old downtown within the first and the second buffers, the new extent of downtown Hefei, which falls between the third and fifth buffers, experienced a dramatic LULC change, as evidenced by 31.75% of newly emerging developed land, 29.59% of forest, and 20.28% of water bodies from other LULC categories. This result can be explained by the importance of administrative zones and their distance to downtown Hefei, since the former city proper of Hefei dominated the outward expansion of developed land, and thus caused urban encroachment into the intensively developed land in neighboring Feidong and Feixi. In contrast, for the exurban counties and districts falling between the sixth and ninth buffers, the junior administrative importance, terrain constraints, and the far distance to downtown Hefei limit the intensive growth of developed land within their boundaries. Thus, the overall area changes in the non-developed land were considerably lower.  Table 2). Figures 4 and 5 show that the spatiotemporal patterns of summertime SUHI are in good agreement with that of the developed land (see Figure 2). The remarkable increase in SUHII levels and SUHISE can be attributed to the changing urban form and spatial linkage with neighboring county jurisdictions. As shown, in 1995-2000, the city proper dominates the spatial pattern of the intensively developed land at the regional level. The city proper exhibits much higher SUHII levels (4-6 levels within the 0-6 buffer zones) than the rest of the study region (1-6 levels). It is noteworthy that, since 2007, a mixture of in-filling and axial expansion occurred at the fringe of the city proper and eventually resulted in downtown Hefei, accompanied by increasing SUHII levels and SUHISE within the buffer zones (6-21 km). Furthermore, in marginal municipality/county jurisdictions (e.g., Chaohu, Lujiang, and Changfeng) and towns and villages within the transportation corridors, the significant increases in summertime SUHII levels and SUHISE can be observed. This indicates that the spatially close linkage between urban and rural developed land exerts a surface warming effect at regional and local levels. Besides, Figure 4 shows the overall trends in zonal SUHII and SUHISE along the urban-rural gradient. As shown, the urbanized and urbanizing areas within the 6-21km buffer zones exhibit relatively higher summertime SUHII and SUHISE than the more distant buffer zones (>21 km). It appears the range of SUHII along the urban-rural gradient is smaller than that of the SUHISE. The spatial pattern of developed land can explain this result because most of the developed land is located within the 6-21km buffer zones; in contrast, the area proportion of developed land in neighboring county jurisdictions is relatively small, and the dominant natural and semi-natural lands can help cool the land surface.  Table 3 shows the overall relationships between the buffer-based predictors/independent variables and the responses (SUHII and SUHISE) in 1995-2016. The results of PLSR models indicate the independent variables are significantly related to the responses, explaining approximately 61.30% of the variance in the SUHII and 79.80% of the variance in the SUHISE within the study region. However, in addition to the annual and seasonal variations of the regional climate, the different units between the buffer-based LULC area proportions and population density make it difficult to directly compare the role of predictors in determining the variances of the responses. Alternatively, the standardized regression coefficients can better indicate the relative strength of the relationships between the buffer-based predictors/independent variables and the responses. Herein, the relative importance of the predictors was interpreted from two perspectives.

Driving Factors Analysis of the Relationship between the LULC Change and SUHI Effect Indicators
When focusing on the relationships between the buffer-based predictors/independent variables and the summertime SUHII, for the predictors with positive coefficients, the sequential order of the relative importance is Developed land > Bare land> Pop_density. As indicated, when controlling for the other independent variables, with each one standard deviation increase in Developed land, Bare land, and Pop_density, the resultant increases in the level of the SUHII standard deviations were 0.627, 0.221, and 0.098, respectively. Given the SUHII was measured with the LST difference between the LULC categories, the relative importance of Pop_density cannot be directly embodied via its standardized coefficient, although there is a highly significant correlation between developed land and Pop_density. In contrast, for the predictors with negative coefficients, the sequential order of the relative importance is cropland > Water body > Forest. When controlling for the other independent variables, with each one standard deviation increase in Cropland, Water body, and Forest, the resultant decreases in the level of the SUHII standard deviations were 1.106, 0.16, and 0.104, respectively. Overall relationships between the buffer-based predictor variables and SUHISE were highly similar to the abovementioned buffer-based predictors versus SUHII relationships, except that the relative importance of Pop_density slightly surpasses that of Bare land. Such a change in the relative importance of Bare land can be explained by its thermal feature but small area proportion in the LULC categories.

On the Relationship Between LULC Change and Summertime SUHI Effect
On local and regional scales, the occurrence of the SUHI effect in response to LULC change should be seen in the broad context of human activities rather than the LULC change itself. Both the biophysical features of the land surface and socioeconomic levels of the study areas play a critical role in causing LULC change and the associated SUHI effect. Most existing studies have focused on the relationship between the biophysical features of the land surface (e.g., NDVI, MNDWI, and NDBI) and associated SUHI effect [69,70]. Usually, due to the unavailability and uncertainty of socioeconomic variables, few studies have focused on the role of socioeconomic variables in determining the SUHI effect [44,71,72]. In some existing studies, ordinary least square (OLS)-based univariate and multiple linear regression models and non-linear regression models were employed to examine the simple and sophisticated relationships between LULC change and the SUHI effect, depending on the data availability and different research purposes. As previously noted, the univariate linear/non-linear regression models sometimes performed better, but the results of the multiple linear regression models were poorly interpreted when considering all of the potential variables that were assumed to be reasonable in the building the models [15,73]. However, this does not mean the failure of comprehensive models in this domain. PLSR is an explicit approach capable of removing multicollinearity among the variables by extracting from the predictors a set of orthogonal factors called latent variables which have the best predictive power [66]. In this study, we used PLSR models to examine the statistically significant relationships between the buffer-based independent variables and the dependent variables (SUHII and SUHISE), in addition to the relative importance of the independent variables in determining the complicated relationships. Nevertheless, the relative importance of Pop_density in contributing to the summertime SUHI effect indicators cannot be directly detected and thus may be underestimated, because of the simply measured LST difference between the impervious surface and cooling surface, and the lack of consideration of the interactions of the LULCs and Pop_density. Thus, we argue that a comprehensive model considering the joint influences of the biophysical features of LULCs and socioeconomic variables can help better understand the sophisticated relationships between LULC change and the SUHI effect, given the existence of the noise of potential multicollinearity in the regression models.

Implications for Sustainable Land Development of a Better Urban Thermal Environment
Similar to the case of other big cities in central China, since the late 1970s, the socioeconomic development of Hefei has lagged behind that of the coastal megacities. This situation changed in 2002 when Greater Hefei implemented its development strategies for joining with the regional economic integration of the Yangtze river basin, through spatial development planning, urban agglomeration, industrial restructuring, and ecological conservation [74]. To date, most of these goals have been successfully achieved and laid the foundation for the expansion of Greater Hefei in the Yangtze River basin. In this study, as shown in Figure 2, the dramatic change in the expansion of downtown Hefei, which featured dense transportation corridors and an in-filling pattern of intensively developed land, largely reflected the LULC change in response to the land development pressure for the newly emerging high-tech parks, upgrading industrial clusters, supporting facilities, and residential communities.
Based on our results, it appears that the imbalance of the land development and population aggregation between the downtown Hefei and the marginal municipality/county jurisdictions should be responsible for the heterogeneity of LULC change and the summertime SUHI effect of the study region. Therefore, we argue that the relationship between the LULC change and associated SUHI effect should be stressed in the local authority's long-term plan actions and developmental strategies. For instance, the local authority can take advantage of the unprecedented opportunity provided by national strategies such as "Rise of central China", "Integration of the Yangtze river basin", and "Made in China 2025", which will favor internationally competitive industrial clusters, such as integrated circuits, new energy, intelligent connected vehicles, and quantum communication, in addition to the tertiary industry in Greater Hefei. However, it is notable that the local authority's recently issued plan (2018-2035) for these leading industrial clusters [75] prefers a site selection (491 km 2 ) within the 6-21 km buffer zones of downtown Hefei, given the convenience of livability and management. Nonetheless, no specific concerns on LULC change and the SUHI effect were embodied in the local authority's planning report, nor was there mention of a practical approach for suppressing the possible adverse thermal effect. It can be anticipated that such an intensive land developmental pattern will not only exacerbate land use competition for commercial development and ecological conservation, but will also not be helpful in shaping a balanced land development pattern between downtown Hefei and the marginal municipality/county jurisdictions. Consequently, the intensely developed land and over-dense population in downtown Hefei make it difficult to exploit sufficient space for spatially configured green spaces and water bodies for effective modification of the urban climate.
Thus, to alleviate the summertime SUHI effect in extreme heat events, future urban planning should reshape the urban form to suppress the over-intensive development in downtown Hefei and make room for cooling surfaces, by propagating the practices of ecological resilience and recovery to rebuild blue-green spaces across the urban-rural corridors. Based on our analysis of the present ecological land layout, we propose the planning of wetland parks and forest parks in the southeast and northwest parts of Greater Hefei, particularly to strengthen the protection zone of Chaohu Lake and optimize the low-impact development of the region. In addition, major ecological corridors, such as the Nanfei River ecological corridor and the Pai River ecological corridor, and the vast area of primary farmland area (for which urban development is permanently forbidden) should be incorporated into the ecological network to control the growth of construction land (see Figure 6).

Limitations of This Study and Remarks
Several minor limitations exist in this study. Theoretically, the time series of available satellite images may provide sufficient information for better understanding the process of LULC change and resultant levels of SUHII and SUHISE, which are highly important to decision making for sustainable land development, climate adaption, and ecosystem management. Unfortunately, due to cloud contamination and the 16-day revisiting interval of Landsat satellites, the availability of satellite images with satisfactory quality is limited. Thus, the results of this study can only partially reveal the relationship between the LULC change and the associated daytime SUHI effect. Furthermore, the site-specific air temperature recorded at the weather stations should also be referenced for better representation of the spatiotemporal patterns of the summertime SUHI effect associated with LULC change. Usually, urban encroachment in suburban and rural areas causes bias in air temperatures at local weather stations, which is known as the "city-entering" phenomenon of stations [64]. However, when overlapping the LULC and summertime SUHI indicator maps with the weather stations map, and then reclassifying these weather stations to their correct suburban or rural area, the quantitative relationship between LULC change and SUHI effect discussed in this study will be more convincing.
Another limitation is the relatively coarse resolution of the SUHI indicator products. The retrieved 30 m SUHI indicator products are suitable to depict the thermal effect of LULC change at medium and large scales with simple landscapes (e.g., forest, water bodies, cropland, and their mosaic), but they are insufficient to characterize the divergence in thermal effect of intensively developed land, especially in the relatively small-scale urban built environment with complicated LULC components and sophisticated urban morphology. Although no very high-resolution TIR band products are provided by the Landsat series satellite platforms, a few trial studies of thermal sharpening technology using a combination of Landsat TIR bands data and high-resolution satellite images have proven effective in monitoring the fine-scale SUHI effect in urban areas [76][77][78]. Moreover, the in situ measured air temperature and the CUHI effect cannot be directly compared with the SUHI effect. However, the combination of thermally sharpened SUHI effect indicators and in situ measured microclimatic factors can be used to generate the simulated CUHI indicators via computational fluid dynamics (CFD), which can be compared with the real CUHI indicators. Therefore, to enhance the applicability and feasibility of satellite remotely sensed TIR products in urban planning for UHI mitigation and climate adaption, more case studies using this thermal sharpening technology or similar methods should be encouraged.

Conclusions
In this study, the major findings can be summarized as follows.
(1) Overall, in Greater Hefei, the heterogeneous patterns of LULC change, particularly the noticeable differences in the growth patterns of the developed land between downtown Hefei and the marginal municipalities/counties, are observed. The centric buffer-based urban-rural gradient analysis reveals that most of the newly emerging developed land occurred within downtown Hefei, evidencing this city's transition from an ordinary provincial capital to a star megacity of central China during the past two decades. Consequently, the dramatic increase in developed land and slight increases in forest and water bodies caused a sharp decrease in cropland. In contrast, in three marginal municipality/county jurisdictions, the overall area change in non-developed land was much lower.
(2) The spatiotemporal patterns of summertime SUHI are in good agreement with that of the developed land. Rapid expansion of the developed land resulted in the remarkable increase in summertime SUHII levels and SUHISE. These were characterized by a mixture of in-filling and axial expansion along transportation corridors, or solely by axial expansion. The results of PLSR models indicate the zonal buffer-based predicators/independent variables are significantly related to the responses (SUHII and SUHISE), explaining approximately 61.3% of the variance in the SUHII and 79.8% of the variance in the SUHISE within the study region. In addition, the relative strength of the independent variables in determining the above-mentioned relationships was quantitatively examined. Therefore, optimization of the land development patterns is the key to mitigating the summertime SUHI effect and climatic adaption.
In summary, the findings of this study provide clear evidence for decision making toward sustainable land development and mitigating the pronounced summertime UHI effect.
Author Contributions: Hao Zhang and Rui Zhou conceived the central idea and designed the technical framework for this study. Ying-ying Li and Hao Zhang commanded data processing and analysis. Hao Zhang performed statistical analysis and wrote the manuscript. Manjula Ranagalage, Hao Zhang and Rui Zhou edited the draft manuscript. Yu Liu., M.R., and Rui Zhou participated in discussing the findings of this study and helped enhance the research design, analysis, and manuscript writing. Specially, the authors are indebted to Prof. Conghe Song, who reviewed this manuscript and made rigid comments on enhancing the quality of this manuscript. All authors have read and agreed to the published version of the manuscript.