Exploring Adaptive UHI Mitigation Solutions by Spatial Heterogeneity of Land Surface Temperature and Its Relationship to Urban Morphology in Historical Downtown Blocks, Beijing

: Heat stress brought on by the intensiﬁcation of urban heat island (UHI) has caused many negative effects on human beings, which were found to be more severe in highly urbanized old towns. With the inconsistent ﬁndings on how urban spatial morphological characteristics inﬂuence land surface temperature (LST) and gaps between design practices being found, we chose Beijing Old Town (BOT) as the study area and took the basic planning implementation module “block” as a study to reveal the spatial heterogeneity of LST and its relationship to multiple urban morphological characteristics with higher spatial resolution calculated via WorldView3. Our results have shown that (1) UHI effect was signiﬁcant and spatially heterogeneous in BOT, and signiﬁcant hot areas with high LST value and small LST differences were found, as cold areas were the exact opposite. (2) The proportion of vegetated area, water, impervious surface, and urban spatial structure indicators i.e., building coverage ratio, mean height, highest building index, height ﬂuctuation degree, space crowd degree and sky view factor were identiﬁed as signiﬁcantly affecting the LST of blocks in BOT. (3) The effects of GBI components and conﬁguration on LST varied within different block types; generally, blocks with GBI with larger patches that were more complex in shape, more aggregated, and less fragmented were associated with lower LST. Finally, in the context of integrating our study results with relevant planning and design guidelines, a strategy sample of adaptive GBI planning and vegetation design for blocks with different morphological features was provided for urban planners and managers to make a decision on UHI mitigation in the renewal process of BOT. coefﬁcients less than Thus, all above morphological indicators identiﬁed to be signiﬁcant correlated with LST mean were selected to classify block types in BOT. The sum of the squared errors method (SEE) was applied in Python v3.8.5 to determine the optimal number of clustering “K” [79]. All indicators related were standardized in SPSS 25.0, and K-means clustering was performed afterwards. Partial correlation analysis was performed between the GBI pattern indicators and LST mean & LST std within different block types. Below, Figure shows the ﬂow chart of the methods and data pre-processing for the present study.


Introduction
It has long been recognized that urbanization changes spatial morphology and brings about population aggregation of cities, thus resulting in various urban problems [1]. One of the most reported examples is that urban areas are warmer than the surrounding countryside, a phenomenon known as the urban heat island (UHI) [2]. The UHI effect has been observed worldwide, contributing to a suite of negative impacts such as increasing energy consumption [3,4] and urban smog formation [5], reducing thermal comfort of urban dwellers, and increasing health problems [6][7][8] and mortality rates [9,10]. Consequently, the UHI effect of urban spatial morphology and how to mitigate it via Nature-based Solution (NbS) has attracted increasing attention from scientists and urban managers [11][12][13][14]. Lately, the UHI effect has been generally studied by remote-sensed land surface temperature (LST) or station-based air temperature [12,13,15]. Due to convenience and effectiveness in acquiring spatial variations of temperature, LST is widely applied to analyze the spatial pattern of UHI and its relationship to urban morphology [16][17][18].

Planning Orientation of Urban Thermal Environment Improvement Studies
As the basic unit of urban construction in China, block is the main unit used to manage land use, control development intensity, allocate public facilities, etc. [53]. It is also the basic scale for shaping urban morphology [54]. Therefore, it is very meaningful to study the relationship between LST and urban morphology at this scale, and to explore the applicable thermal mitigation strategies for urban planning (UP). Some researchers have conducted studies on block level in Chinese cities such as Wuhan and Xi'an [22,24]. They revealed the significant impact of urban characteristics on LST and proposed some directions for optimizing urban 2D or 3D characteristics; however, these findings can be further refined to form a research strategy and integrate the results with local planning documents to enhance the value of practical application.
In order to offer practical, operational guidelines for urban managers, the selection of morphological indicators and the research process should be more comprehensively considered. Some researchers mentioned that the indicators selected for the study should be those that can be directly used or provide specific directions for urban planning; e.g., building height, floor area ratio, building coverage ratio, and so on [28,29]. As for indicators that are not used in UP practices but characterize the urban form in a comprehensive way, their significance for guiding UP should be explored [24]. In addition, urban planning and design has a comprehensive process and a clear sequence, in which after the spatial layout of the city is clearly defined and a regional infrastructure planning framework is formed, various important urban construction projects can be drawn up; i.e., the construction at the three-dimensional level is restricted by the planning at the two-dimensional level [54]. Most of the current studies have analyzed different levels of indicators simultaneously during the analysis [30,55], which may lead to interference in the analysis results.
Therefore, this study aims to explore the adaptive UHI mitigation solutions based on the correlations between fine-texture morphological indicators and LST in old downtown areas by using imagery with high spatial resolution. We conducted this research in the historical downtown blocks in Beijing City, China. Specifically, it addresses three questions as follows: (1) what are the LST characteristics of blocks and how do these LST characteristics differ among different blocks in Beijing Old Town? (2) Which urban morphological indicators influence the LST significantly in these blocks? (3) What are the differences in relationships between the LST and GBI pattern among blocks with different morphological features?
Land 2022, 11, 544 4 of 24 Important insight of adaptive GBI planning and vegetation design can be provided for urban planners and managers to make decisions about mitigating the UHI effect in the old city renewal process.

Study Area and Block Units
Beijing Old Town (BOT), located in the center of Beijing City (39 • 28 ~41 • 50 N, 115 • 25 ~117 • 30 E), with a total area of about 62.5 km 2 , was defined as the area within the Second-Ring Road by the Urban Master Plan of Beijing (2016-2035) [56] (Figure 1). BOT has been developed for 3000 years, and has been the capital city for 800 years, where not only cultural relics and historic sites are richer, but where a greater variety of buildings exist and a larger population aggregates. BOT acts as the most key area for the preservation of urban historical and cultural heritage in Beijing. Nevertheless, there are both significant hot and cool spots of UHI observed in this area [57]. Consequently, efforts to mitigate the UHI effect should strictly implement the requirements of "the old town can no longer be demolished". old city renewal process.

Study Area and Block Units
Beijing Old Town (BOT), located in the center of Beijing City (39°28′~41°50′ N, 115°25′~117°30′ E), with a total area of about 62.5 km 2 , was defined as the area within the Second-Ring Road by the Urban Master Plan of Beijing (2016Beijing ( -2035 [56] (Figure 1). BOT has been developed for 3000 years, and has been the capital city for 800 years, where not only cultural relics and historic sites are richer, but where a greater variety of buildings exist and a larger population aggregates. BOT acts as the most key area for the preservation of urban historical and cultural heritage in Beijing. Nevertheless, there are both significant hot and cool spots of UHI observed in this area [57]. Consequently, efforts to mitigate the UHI effect should strictly implement the requirements of "the old town can no longer be demolished".
According to the Detail Regulatory Plan in the Functional Core Area of Beijing (Block Level) [58], BOT involves 25 districts that are composed of 123 blocks (including 47 historical and cultural blocks, HCB) ( Figure 1a). HCB refers to blocks with historical value that need to be protected, in which they are dotted with many traditional housing quadrangles (siheyuan) and alleyways (called as hutongs). Some of the most famous hutongs such as Dashilar Commercial Street (Figure 1b

Data Source and Processing
We obtained the LST and urban morphological indicators by remote sensing image data and building vector data. The former, including Landsat-8 OLI Imagery (September, 2017) from the Geographical Spatial Data Cloud (http://www.gscloud.cn (accessed on 9 According to the Detail Regulatory Plan in the Functional Core Area of Beijing (Block Level) [58], BOT involves 25 districts that are composed of 123 blocks (including 47 historical and cultural blocks, HCB) ( Figure 1a). HCB refers to blocks with historical value that need to be protected, in which they are dotted with many traditional housing quadrangles (siheyuan) and alleyways (called as hutongs). Some of the most famous hutongs such as Dashilar Commercial Street (Figure 1b

Data Source and Processing
We obtained the LST and urban morphological indicators by remote sensing image data and building vector data. The former, including Landsat-8 OLI Imagery (September, 2017) from the Geographical Spatial Data Cloud (http://www.gscloud.cn (accessed on 9 September 2020)) and WorldView-3 Imagery (October 2018) with a high resolution of 0.30 m, are used to extract the LST value and landscape characteristic indicators [46]. They are high-quality images in which cloud cover in the study area is less than 3%, and the weather on these days was mainly sunny with light wind; that is, the thermal environment was not affected by rainfall or wind [59]. Based on the ENVI5.3, the remote sensing data was pre-processed through radiometric calibration, atmospheric correction, and orthorectification, as well as image fusion, cropping, and enhancement [60]; the latter are the Amap building vector data (2016). After spatial correction, we applied it to calculate the spatial structure indicators of built-up.

Characteristic Indicators of UHI Based on LST
The radiative transfer equation (RTE) method was applied to retrieve the LST in the thermal infrared (TIR) band (10.60~11.90 µm) of pre-processed Landsat-8 OLI using ENVI 5.3, the specific calculation steps of which are described below [60]: where B(T S ) is the blackbody radiance, L λ is the thermal infrared radiance brightness value; ε is the surface emissivity. The atmospheric profile parameters include τ, L ↑ , and L ↓ , where τ is the atmospheric transmittance in the thermal infrared band, L ↑ is the atmospheric upward radiance, and L ↓ is the atmospheric downward radiance; parameters were obtained from media files of USGS ((https://www.usgs.gov (accessed on 3 March 2021)). ε was calculated differently for different regions with following formulas [61]: where ε wt , ε ua , ε ns are surface emissivity of water body, urban area and natural surface, respectively. The equation for P v (vegetation coverage) is [62]: where NDV I is the normalized vegetation index of the whole BOT; N Soil is the NDVI of the completely bare, soil-covered area; N Vegetaion is the NDVI value of areas with complete vegetation covered. Based on empirical values, N Soil equals 0.05 and N Vegetaion equals 0.70 [59]. Planck's formula was further applied to calculate the land surface temperature T s [63]: where T s is the surface brightness temperature; constants K 1 = 774.89 W/(m 2 µm sr), K 2 = 1321.08 K [22]). The mean value (LST mean ), standard deviation (LST std ), and value range (LST range ) of LST were used to quantify the variation of thermal environment in BOT. Among them, LST mean characterized the thermal intensity of urban blocks, the larger its value, the higher the LST, the stronger the UHI intensity. LST std and LST range characterized the degree of thermal variation of urban blocks, the smaller its value, the smaller the LST difference, the lower the thermal dispersion intensity [64]. Meanwhile, the Getis-Ord Gi* local statistics of LST mean was also conducted to demonstrate the spatial heterogeneity of UHI (hot, warm, Land 2022, 11, 544 6 of 24 and cold area) based on ArcGIS 10.6, as it was considered the most appropriate method to identify hot spots [65]. The equations for calculating the Getis-Ord Gi* are [66]: where i, j are designations of different raster cells/elements; w i,j is the spatial weight between elements i and j; n is the number of raster cells; x is the raster cell attribute; X is the average of attribute values [66]. The cohesion index and the aggregation index of hot, warm, and cold areas were calculated with the software of Fragstats 4.2 [67].

Spatial Morphological Indicators of Urban Blocks
Here, characteristic indicators of land cover and spatial structure indicators of builtups (Table 1) were chosen to quantitatively describe urban morphology in BOT. The land cover characteristic indicators were calculated by the land cover map based on WorldView-3 Imagery. Based on the eCognition Developer 9.0, remote sensing interpretation was carried out using the pre-processed WorldView-3 images [68]. Basic units were divided using multi-scale segmentation, then the object-oriented classification method was applied to firstly distinguish the water body and non-water body areas based on values of NIR band, after which NDVI was defined in the software to distinguish the vegetated and non-vegetated areas [69]. Additionally, combining with the human-computer interaction method, grassland areas were distinguished from trees and shrubs. Final results were further corrected by visual interpretation. The accuracy was verified by the field accuracy check with sample points. Results showed that the overall interpretation accuracy was 93.7% with the kappa coefficient up to 0.92, which met the requirements of the experimental data analysis [70]. The Normalized Difference Vegetation Index (NDVI) and the Impervious Surface Area (ISA) were selected to represent the vegetation quality and the impervious status of the land surface. The urban spatial structure indicators were calculated based on modified Amap building vector data [71]. Formulas of these indicators are detailed in Table A1. Block-based zonal statistics were conducted by using images after interpretation in ArcGIS10.6. Table 1. Description of land cover and urban spatial structure indicators.

Indicators Definition Unit
Land cover characteristics Impervious land proportion (IP) The ratio of impervious area to block area. % Vegetated land proportion (VP) The ratio of vegetated area to block area. % Water proportion (WP) The ratio of water area to block area. % Bare soil proportion (SP) The To find the mitigation measures of UHI in urban blocks, spatial pattern indicators of GBI were further selected in Table 2. The composition indicators of TP, GP and WP were calculated from interpretation results of WorldView3. By merging trees, shrubs, grass and water in ArcGIS10.6, we obtained the GBI raster map and calculated its spatial configuration within blocks by the software of Fragstats4.2. Formulas of the configuration indicators were detailed in Table A2.

Relationship Analysis Methods
The Pearson correlation analysis was conducted between urban morphology and the LST with the help of SPSS25.0. Pearson correlation test was also conducted in order to avoid the effect of the multicollinearity of indicators. Variables with coefficients larger than 0.6 have a strong correlation to each other and one of them should be excluded from the following analysis [77]. By using the absolute values of correlation coefficients, we explained the influence and dominance of the morphological indicators on the thermal environment. Indicators with significance coefficients less than 0.05 were considered being significantly correlated with the thermal environment [78].
Urban green space planning is dependent on block morphological characteristics. Thus, all above morphological indicators identified to be significant correlated with LST mean were selected to classify block types in BOT. The sum of the squared errors method (SEE) was applied in Python v3.8.5 to determine the optimal number of clustering "K" [79]. All indicators related were standardized in SPSS 25.0, and K-means clustering was performed afterwards. Partial correlation analysis was performed between the GBI pattern indicators and LST mean & LST std within different block types. Below, Figure 2 shows the flow chart of the methods and data pre-processing for the present study.

EER REVIEW 8 of 24
explained the influence and dominance of the morphological indicators on the thermal environment. Indicators with significance coefficients less than 0.05 were considered being significantly correlated with the thermal environment [78]. Urban green space planning is dependent on block morphological characteristics. Thus, all above morphological indicators identified to be significant correlated with LSTmean were selected to classify block types in BOT. The sum of the squared errors method (SEE) was applied in Python v3.8.5 to determine the optimal number of clustering "K" [79]. All indicators related were standardized in SPSS 25.0, and K-means clustering was performed afterwards. Partial correlation analysis was performed between the GBI pattern indicators and LSTmean & LSTstd within different block types. Below, Figure 2 shows the flow chart of the methods and data pre-processing for the present study. Flowchart of this study. Note: Getis-Ord Gi*, a hotspot analysis tool, was applied to identify those spatially aggregated high and low LST points with statistically significance.

Spatial Heterogeneity of UHI in Historical Downtown Blocks
Beijing Old Town (BOT) was located in central Beijing, where high LST aggregated ( Figure 3a). The LSTmean of BOT was 33.06 C higher than 27.09 °C of the whole city and 31.15 °C of the central urban area, while LSTstd was 2.09 °C and LSTrange was 25.25 °C, respectively, lower than those of the whole city and central urban area of Beijing (Table  3). The thermal characteristics of higher mean temperature value and lower temperature difference indicated that this historical downtown area was under more severe heat stress.
The thermal environment characteristics within BOT were spatially differentiated among blocks, while the hot and cold areas had conjugated each other (Figure 3a,b). Integrated with the LST value and spatial characteristics of hot, warm and cold areas, we found that the LSTmean of hot areas reached as high as 35.61 °C, while the LSTstd and LSTrange was 0.86 °C and 4.95 °C, respectively (Table 3). Those typical heat islands with higher cohesion and lower connectivity value (Table 3) were mainly aggregated in the Dashilar Figure 2. Flowchart of this study. Note: Getis-Ord Gi*, a hotspot analysis tool, was applied to identify those spatially aggregated high and low LST points with statistically significance.

Spatial Heterogeneity of UHI in Historical Downtown Blocks
Beijing Old Town (BOT) was located in central Beijing, where high LST aggregated ( Figure 3a). The LST mean of BOT was 33.06 C higher than 27.09 • C of the whole city and 31.15 • C of the central urban area, while LST std was 2.09 • C and LST range was 25.25 • C, respectively, lower than those of the whole city and central urban area of Beijing (Table 3). The thermal characteristics of higher mean temperature value and lower temperature difference indicated that this historical downtown area was under more severe heat stress. and connectivity value (Table 3) consisted of the matrix of LST distribution pattern in BOT ( Figure 3c). The LSTmean of warm areas was 33.19 °C that was lower than that of hot areas, while the LSTstd (0.86 °C) was equal to that of the hot area and even the LSTrange (3.88 °C) was lower than that of hot areas (Table 3). This indicated that the warm area also had a feature with high average temperature and smaller temperature differences but was not as significant as hot areas.

Relationship between UHI and Block Morphology
To avoid the multicollinearity of indicators in Table 1, the Pearson's correlation coefficient among initial morphological features is calculated and shown in Figure 4. Here, the NDVI was highly correlated with the vegetation proportion (VP), ISA was highly correlated with the impervious surface proportion (ISP), the building surface area (BSA) was highly correlated with the mean height (MH), the average volume (AV) and the building structural index (BSI), BSI was highly correlated with MH and AV and AV was highly correlated with FAR, BSI and BSA, all of which had a correlation coefficient over 0.6. Because VP, ISP and MH are commonly used in urban planning and design, we excluded  The thermal environment characteristics within BOT were spatially differentiated among blocks, while the hot and cold areas had conjugated each other (Figure 3a,b). Integrated with the LST value and spatial characteristics of hot, warm and cold areas, we found that the LST mean of hot areas reached as high as 35.61 • C, while the LST std and LST range was 0.86 • C and 4.95 • C, respectively (Table 3). Those typical heat islands with higher cohesion and lower connectivity value (Table 3) were mainly aggregated in the Dashilar Block, Xinjiekou Block, Tiyuguanlu Block and Jianguomen Block (Figure 3c(i-iv)). Conversely, the LST mean of cold areas was lower than 30 • C, with the LST std reaching 1.37 • C and LST range up to 17.47 • C. Those areas also had higher cohesion and lower connectivity value (Table 3) and mainly clumped in blocks of Shichahai, Tiantan and Longtan (Figure 3c(v-vii)). Different from the above two types of areas, the warm area with highest cohesion and connectivity value (Table 3) consisted of the matrix of LST distribution pattern in BOT (Figure 3c). The LST mean of warm areas was 33.19 • C that was lower than that of hot areas, while the LST std (0.86 • C) was equal to that of the hot area and even the LST range (3.88 • C) was lower than that of hot areas (Table 3). This indicated that the warm area also had a feature with high average temperature and smaller temperature differences but was not as significant as hot areas.

Relationship between UHI and Block Morphology
To avoid the multicollinearity of indicators in Table 1, the Pearson's correlation coefficient among initial morphological features is calculated and shown in Figure 4. Here, the NDVI was highly correlated with the vegetation proportion (VP), ISA was highly correlated with the impervious surface proportion (ISP), the building surface area (BSA) was highly correlated with the mean height (MH), the average volume (AV) and the building structural index (BSI), BSI was highly correlated with MH and AV and AV was highly correlated with FAR, BSI and BSA, all of which had a correlation coefficient over 0.6. Because VP, ISP and MH are commonly used in urban planning and design, we excluded NDVI and ISA from land cover characteristic indicators, as well as BSI, BSA and AV from spatial structure indicators.

ER REVIEW 10 of 24
NDVI and ISA from land cover characteristic indicators, as well as BSI, BSA and AV from spatial structure indicators. Then, we further conducted the correlation analysis between screened block morphological indicators and the LST. As shown in Table 4, ISP reaching 67.86% of the total area of BOT is significantly positively correlated with LSTmean and negatively correlated with both LSTstd and LSTrange, while it is the opposite for GBI-related indicators such as VP and WP. Nevertheless, the presence of bare soil did not have a significant effect on either LSTmean or LSTstd, nor on LSTrange. Meanwhile, the absolute value of correlation coefficient between ISP and LSTmean was significantly higher than that of VP and WP, indicating that the positive influence of gray components (characterized by ISP) on LSTmean was higher than the negative influence of the blue-green components characterized by VP and WP. The absolute value of correlation coefficient between ISP and LSTstd (as well as LSTrange) was a bit lower than that of WP, but significantly higher than VP, showing that the influence of increasing the temperature difference from WP was greater than that from VP, as well as the influence of reducing the temperature difference from ISP.  Then, we further conducted the correlation analysis between screened block morphological indicators and the LST. As shown in Table 4, ISP reaching 67.86% of the total area of BOT is significantly positively correlated with LST mean and negatively correlated with both LST std and LST range , while it is the opposite for GBI-related indicators such as VP and WP. Nevertheless, the presence of bare soil did not have a significant effect on either LST mean or LST std , nor on LST range . Meanwhile, the absolute value of correlation coefficient between ISP and LST mean was significantly higher than that of VP and WP, indicating that the positive influence of gray components (characterized by ISP) on LST mean was higher than the negative influence of the blue-green components characterized by VP and WP. The absolute value of correlation coefficient between ISP and LST std (as well as LST range ) was a bit lower than that of WP, but significantly higher than VP, showing that the influence of increasing the temperature difference from WP was greater than that from VP, as well as the influence of reducing the temperature difference from ISP.
Considering blocks with higher ISP tended to have higher LST, smaller temperature differences and more gathered hot areas, the relationship between spatial structure indicators of built-up and the LST was further analyzed by partial correlation analysis ( Table 5). The building coverage ratio (BCR), the space crowd degree (SCD) and the sky view factor (SVF) were all positively significant with LST mean , with the absolute value of correlation coefficient compared as BCR > SCD > SVF, while the building height related indicators including the mean height (MH), the highest building index (HBI) and the height fluctuation degree (HFD) were all negatively significant with LST mean , with the absolute value of correlation coefficient compared as MH > HFD > HBI. Nevertheless, BCR was the only indicator with significantly negative correlation with LST std . The above results indicate that the increase of buildings proportion in blocks will aggravate the increasing of LST value and the reducing of LST variation, while the increase of average building height, maximum building index and fluctuation can decrease the LST value but not affect the LST difference. Overall, the BCR and MH were the primary two indicators of built-up's spatial structure affecting the LST in BOT blocks according to the absolute value comparison. Note: C is the coefficient and S is the significance, and ** means at the 0.01 level (two-tailed), the correlation is significant. Note: C is the coefficient and S is the significance, and ** means at the 0.01 level (two-tailed), the correlation is significant. * means at the 0.05 level (two-tailed), the correlation is significant.

Relationships between UHI and GBI in Different Spatial Morphological Blocks
Morphological indicators with significant correlation with LST mean obtained from the above were used to classify blocks with different urban form in BOT. After applying the sum of the squared errors method (SSE) to identify K-value (K = 3 in this study), we classified those blocks into three categories by K-means clustering, which were high GBI proportion and mid-rise blocks (HMB), mid GBI proportion and high-rise blocks (MHB), and low GBI proportion and low-rise blocks (LLB), with a quantity of 30, 45 and 48, respectively (Figure 5a). There existed significant differences of thermal environment and urban morphology among three categories of block types (Figure 5b,c). The LST mean of HMB was 31.57 • C, lower than 32.82 • C for MHB and 34.42 • C for LLB. HMB also had the highest LSTstd and LST range (1.61 • C, 7.02 • C), both higher than that for MHB (1.07 • C, 0.90 • C) and LLB (5.26 • C, 4.26 • C). As for block morphology, HMB had the highest GBI proportion (that is the sum of VP and WP) and the lowest ISP, while LLB was exactly opposite and MHB was in-between ( Figure 5c). LLB has the highest proportion of buildings with the highest spatial congestion and lowest building height, which was significantly different from the other two block types (Figure 5c). Then, we further conducted Pearson correlation statistics to analyze the relationship between the LST and GBI composition and configuration characteristic in three categories of blocks. Here, considering that LSTstd and LSTrange showed similar feature differences among three block types (Figure 5b), only LSTstd and LSTmean was selected for the following correlation analysis. As shown in Table 6, the tree cover proportion (TP) in HMB and LLB while the grass cover proportion (GP) in MHB were presented significantly correlated with LSTmean, but no correlation with LSTstd, and the water cover proportion (WP) in HMB had significant correlation both with LSTmean and LSTstd, and in MHB only had significant correlation with LSTmean, but no correlation with LSTmean and LSTstd in LLB. The absolute values of correlation coefficients in HMB were obviously larger than those in MHB and LLB. The above showed that the relationship between GBI composition and LST varied depending on block types, and the mean LST value might be reduced by increasing the proportion of trees, shrubs, water bodies or grassland according to local conditions.
As for the GBI configuration indicators, their correlations with the LSTmean and LSTstd varied among different blocks, too. Overall, only the AREA_MN was negatively corre- Then, we further conducted Pearson correlation statistics to analyze the relationship between the LST and GBI composition and configuration characteristic in three categories of blocks. Here, considering that LST std and LST range showed similar feature differences among three block types (Figure 5b), only LST std and LST mean was selected for the following correlation analysis. As shown in Table 6, the tree cover proportion (TP) in HMB and LLB while the grass cover proportion (GP) in MHB were presented significantly correlated with LST mean , but no correlation with LST std , and the water cover proportion (WP) in HMB had significant correlation both with LST mean and LST std , and in MHB only had significant correlation with LST mean , but no correlation with LST mean and LST std in LLB. The absolute values of correlation coefficients in HMB were obviously larger than those in MHB and LLB. The above showed that the relationship between GBI composition and LST varied depending on block types, and the mean LST value might be reduced by increasing the proportion of trees, shrubs, water bodies or grassland according to local conditions. Note: TP-trees and shrubs proportion; GP-grass proportion; WP-water proportion; LPI-largest patch index; AREA_MN-mean patch size; LSI-landscape shape index; FRAC_AM-area-weighted fractal dimension index; NP-number of patches; ENN_AM-Mean Euclidean nearest neighbor distance; AI-aggregation index. ** means at the 0.01 level (two-tailed), the correlation is significant. * means at the 0.05 level (two-tailed), the correlation is significant.
As for the GBI configuration indicators, their correlations with the LST mean and LST std varied among different blocks, too. Overall, only the AREA_MN was negatively correlated with LST mean in all block types; the LPI, FRAC_AM, and AI were negatively correlated with LST mean , while the ENN_AM was positively correlated with LST mean significantly in HMB and LLB, and no significant correlation was found in MHB; the LSI and NP were positively correlated with LST mean in MHB. Fewer indicators were identified to be significantly correlated with LST std , which were the LPI and AI in HMB, and AREA_MN and AI in MHB.
By comparing the absolute values of correlation coefficients between LST mean and GBI configuration indicators, we further found that the absolute values of coefficients for indicators that characterized GBI patches' area and their aggregation, including the LPI, AREA_MN and AI, were generally higher than those for other indicators, and the absolute values of coefficients in HMB were also larger than those in LLB ( Table 5). As for the FRAC_AM, the indicator characterizing the shape complexity of GBI patches, the absolute value of coefficient in HMB was found to be higher than that in LLB. Nevertheless, for the indicator of ENN_AM representing the connectivity of GBI, the absolute value of coefficient in HMB were smaller than that in LLB.
As a whole, both the correlation significance between LST and GBI's pattern indicators and the absolute value of their correlation coefficients were significantly different; thus, it is important to emphasize the spatial morphological dependence of the GBI pattern influence on the UHI effect in BOT blocks.

Potential Adaptive UHI Mitigation Solutions for Historical Downtown Blocks with Different Spatial Morphology
As a typical historical downtown area, Beijing Old Town (BOT) has been strictly protected to preserve the original build style and cultural features; any planning and design acts must be initiated on the premise of protecting the historical and cultural heritage [56,58,80]. Thence, a comprehensive and feasible GBI planning design for thermal mitigation in BOT might consider at least the following three aspects: urban morphology characteristics, the composition and configuration of GBI, as well as the policy terms and control intensity. Integrated with a broad review of requirements and guidelines of urban planning and management in BOT, we extracted the GBI measures permitted in BOT, and provided a set of thermal mitigation measures by GBI, as Table 7 showed. Based on this measure set, we provided an adaptive thermal mitigation solution using GBI for each block (Figures 6-8). Table 7. Summary of GBI measures to mitigate the UHI in Beijing Old Town.

Key Points Requirements Sources
Promoting street shading Improve the construction of boulevard system, promote the transformation of street shading facilities. Regulatory plan, 2020.
Optimize plant communities and improve the quality of street shade. Regulatory plan, 2020.
Increasing the area and accessibility of green space Promote greening coverage, green space per capita, etc.
Regulatory plan, 2020. Constructing both community parks and small pocket green spaces on marginal land, unused land and etc.

Key Points Requirements Sources
Promoting the tree protection Integrate trees preservation with street greening or public space design programs, strictly implement the protection of famous trees, large trees, old and historic trees.
Carry out the "one tree for one courtyard" replanting program. Regulatory plan, 2020. Promoting shading transformation of checkerboard-shaped street network Enhance the proportion of boulevards, improve greening rate of street network by increasing the green space on both sides of streets.

Strengthening preservation and maintenance of traditional buildings and alleys
Strictly implement the requirements of historical buildings or buildings with traditional features can no longer be demolished. • High GBI proportion and mid-rise blocks (HMB) These blocks are mainly located in the central and southern regions of Beijing Old Town; e.g., West Chang'an Block, Donghuamen Block, Tiantan Block, Taoranting Block, etc. There are larger bodies of water, parks and historical gardens concentrated in some blocks, and the land surface temperature is relatively low. Spaces available for greening in such blocks are relatively limited and scattered due to historical and cultural heritage protection regulations. Hence, the priority in the GBI planning and design is the optimization of GBI composition such as increasing trees, shrubs and water bodies in HMB blocks.
As for the optimization of GBI configuration, priority is firstly given to greening the paved, vacant or unused lands in the residential clusters where heat islands exist (Figure 5a). New green spaces also might be integrated with the renovation of shanty towns, dangerous houses, economical buildings and public spaces, which is helpful to increase the average area of GBI patches in HMB. Most of the existing parks, gardens and green patches in HMB should be further enlarged or connected by building greenbelts and shaded walkways if possible. Higher green view rate (GVR) often indicating abundant vegetation and better plant communities is better for cooling. The boulevard with comfortable street shading might be considered by urban planners, which is a GBI type with complex shape and high connectivity. However, it is notable for areas around royal palaces and gardens, ancient government offices, altars and temple architecture clusters that the greening form, including plant selection and layout, should be consistent with their own spatial style. average area of GBI patches in HMB. Most of the existing parks, gardens and green patches in HMB should be further enlarged or connected by building greenbelts and shaded walkways if possible. Higher green view rate (GVR) often indicating abundant vegetation and better plant communities is better for cooling. The boulevard with comfortable street shading might be considered by urban planners, which is a GBI type with complex shape and high connectivity. However, it is notable for areas around royal palaces and gardens, ancient government offices, altars and temple architecture clusters that the greening form, including plant selection and layout, should be consistent with their own spatial style. These blocks are more concentrated in the east and west sides of Beijing Old Town, mainly dominated by residential and commercial buildings (e.g., Wangjing SOHO and Chaonai Block and Fenghuiyuan Block), with almost no historical and cultural blocks inside. There is a high proportion of tall buildings that provide both ground spaces for greening and shade, which contributes to the fact that the area and intensity of heat islands in MHB is not as high as in HMB. However, a certain proportion of bare lands exists in MHB. Thus, the planning and design of GBI to mitigate thermal stress gives priority to retrofitting these bare lands into community parks and small pocket parks, or part playgrounds into lawns. These renewal measures could be integrated with the Urban Master Plan of Beijing (2016-2035) that requires such neighborhoods (non-historical and cultural blocks, the main components of MHB) to achieve a 500 m service radius coverage of green spaces, and to achieve the goal of park area per capita >6.8 m 2 by 2035.
Considering that the shape complex and fragmentation of GBI has a significant impact on thermal environment in MHB, it is also encouraged to enlarge isolated green patches and to connect them with boulevards into a network. There are many feasible ways to enhance the GBI network such as building linear greenbelts along streets, enlarging small gardens around buildings and improving vertical greening of buildings. Here, green stormwater infrastructure and multilayer plantation could be integrated with greenbelt planning, as well as green roof or roof gardens with energy-saving building transformation. All of the above can make for great contributions to the formation of a green, cool, shared nature space system for all.  • Low GBI proportion and low-rise blocks (LLB) These blocks with a high percentage of historic and cultural blocks, such as Dashilar Block and Andingmen Block, are mainly located in the core and north of Beijing Old Town. Dense bungalows and two-storied houses, narrow public spaces and stricter protection requirements are among the top three challenges posed to thermal mitigation by way of GBI in LLB. It is most urgent to increase the GBI proportion and improve the GBI connectivity as much as possible; it is more feasible to greening streets and courtyards in LLB.

Discussion
It is increasingly acknowledged that urban morphology has an important influence on the formation and distribution of UHI in cities [16,32,33,81]. Previous studies have paid more attention to the spatial heterogeneity of UHI and its relationship with spatial morphology in a city, although blocks are taken as the analysis units. Our study extends these studies by focusing on the finer spatial morphology and its influence on UHI in the highly urbanized downtown blocks, using Worldview3 remote sensing data with spatial resolution of 0.3 m.
We found there was a significant UHI effect in Beijing Old Town, a famous and typical historical downtown region of Beijing in China, as reflected by the higher average surface temperature (LSTmean of 33.06 °C) and smaller temperature difference (LSTstd and LSTrange were 2.09 °C and 25.25 °C, respectively). It was demonstrated as the contiguous warm matrix and concentrated hot patches in Figure 3c. The combined application of LSTmean, LSTstd, LSTrange and Hot Spot Analysis (Getis-Ord Gi *) can better reflect the spatial heterogenous characteristics of UHI in highly urbanized blocks.
The correlation analysis between the LST indicators and morphological characteristics revealed that the impervious surface, vegetation and water cover proportion were the main land cover indicators influencing the average surface temperature and temperature differences in the study area, which was consistent with the results of most related studies [7,47,82]. Unlike previous studies which rarely associated the correlation analysis with the process of urban planning, this study conducted partial correlation analysis between the LST indicators and spatial structure feature of built-up by controlling the impervious surface, vegetation and water cover proportion.
Results showed that the building coverage ratio (BCR) and sky view factor (SVF) were positively correlated within the LSTmean, which were consistent with findings of most studies [24,34,35,83]. The space crowd degree (SCD)-a less used indicator-was also positively correlated with LSTmean, which reflected the difference of outdoor space encroachment by buildings in blocks of Beijing Old Town; thus, this indicator was very applicable in this study area. Meanwhile, all height-related indicators, including MH, HBI and HFD, were negatively correlated with LSTmean, which was consistent with the findings of most studies [30,[84][85][86][87]. This probably resulted from increasing and differentiating building height in one area that can lead to shadows provision, wind corridors enhancement, building land saving and blue or green spaces reservation, thus contributing to UHI mitigation. A comprehensive comparison of the absolute values of correlation coefficients showed These blocks are more concentrated in the east and west sides of Beijing Old Town, mainly dominated by residential and commercial buildings (e.g., Wangjing SOHO and Chaonai Block and Fenghuiyuan Block), with almost no historical and cultural blocks inside. There is a high proportion of tall buildings that provide both ground spaces for greening and shade, which contributes to the fact that the area and intensity of heat islands in MHB is not as high as in HMB. However, a certain proportion of bare lands exists in MHB. Thus, the planning and design of GBI to mitigate thermal stress gives priority to retrofitting these bare lands into community parks and small pocket parks, or part playgrounds into lawns. These renewal measures could be integrated with the Urban Master Plan of Beijing (2016-2035) that requires such neighborhoods (non-historical and cultural blocks, the main components of MHB) to achieve a 500 m service radius coverage of green spaces, and to achieve the goal of park area per capita >6.8 m 2 by 2035.
Considering that the shape complex and fragmentation of GBI has a significant impact on thermal environment in MHB, it is also encouraged to enlarge isolated green patches and to connect them with boulevards into a network. There are many feasible ways to enhance the GBI network such as building linear greenbelts along streets, enlarging small gardens around buildings and improving vertical greening of buildings. Here, green stormwater infrastructure and multilayer plantation could be integrated with greenbelt planning, as well as green roof or roof gardens with energy-saving building transformation. All of the above can make for great contributions to the formation of a green, cool, shared nature space system for all.

•
Low GBI proportion and low-rise blocks (LLB) These blocks with a high percentage of historic and cultural blocks, such as Dashilar Block and Andingmen Block, are mainly located in the core and north of Beijing Old Town. Dense bungalows and two-storied houses, narrow public spaces and stricter protection requirements are among the top three challenges posed to thermal mitigation by way of GBI in LLB. It is most urgent to increase the GBI proportion and improve the GBI connectivity as much as possible; it is more feasible to greening streets and courtyards in LLB.
Greening priority is given to hutongs with more serious heat islands. Greening corridors by planting street trees, bioswales or flower stands along wider hutongs is more advocated because of the better connectivity and shape complexity. Street trees in LLB are mainly tall species with larger canopies such as acacia and elm, mostly old or famous. So, the protection of large trees should be firmly implemented to improve their maintenance, the plant community structure of trees and shrubs and grasses that can be enriched around large, old trees.
For the narrow streets and alleys, the potential spaces for greening should be various and the greening form should be flexible due to the constraint of spaces. The building setback spaces, unused paved frontiers, reserving micro-courtyards, paved spaces around public transport stations in LLB could be greened into corner gardens, "one-meter gardens" or hedgerows, etc. Integrated with pavilions and planters, the planting ponds, flower beds or flower stands could also be used to green streets, walls and courtyards. All plants should give priority to persimmon, acacia, magnolia, begonia and other native plants with traditional cultural characteristics and ease to grow natively. For example, traditional plants such as Parthenocissus tricuspidata can be used on the greening of walls or houses, which can promote cooling and humidification benefits at the pedestrian height. Here, it is worth emphasizing that public participation is very important to courtyard and hutong greening in LLB.
Moreover, asphalt pavements are most used in hutongs; these pavements with dark color and low albedo will lead to deterioration of the thermal environment. Those materials might be replaced by high albedo materials in an urban renewal solution, giving priority to the traditional pavement materials with light color and high reflectivity such as strips of stone, stone slabs and gray processed bricks.

Discussion
It is increasingly acknowledged that urban morphology has an important influence on the formation and distribution of UHI in cities [16,32,33,81]. Previous studies have paid more attention to the spatial heterogeneity of UHI and its relationship with spatial morphology in a city, although blocks are taken as the analysis units. Our study extends these studies by focusing on the finer spatial morphology and its influence on UHI in the highly urbanized downtown blocks, using Worldview3 remote sensing data with spatial resolution of 0.3 m.
We found there was a significant UHI effect in Beijing Old Town, a famous and typical historical downtown region of Beijing in China, as reflected by the higher average surface temperature (LST mean of 33.06 • C) and smaller temperature difference (LST std and LST range were 2.09 • C and 25.25 • C, respectively). It was demonstrated as the contiguous warm matrix and concentrated hot patches in Figure 3c. The combined application of LST mean , LSTstd, LST range and Hot Spot Analysis (Getis-Ord Gi*) can better reflect the spatial heterogenous characteristics of UHI in highly urbanized blocks.
The correlation analysis between the LST indicators and morphological characteristics revealed that the impervious surface, vegetation and water cover proportion were the main land cover indicators influencing the average surface temperature and temperature differences in the study area, which was consistent with the results of most related studies [7,47,82]. Unlike previous studies which rarely associated the correlation analysis with the process of urban planning, this study conducted partial correlation analysis between the LST indicators and spatial structure feature of built-up by controlling the impervious surface, vegetation and water cover proportion.
Results showed that the building coverage ratio (BCR) and sky view factor (SVF) were positively correlated within the LST mean , which were consistent with findings of most studies [24,34,35,83]. The space crowd degree (SCD)-a less used indicator-was also positively correlated with LST mean , which reflected the difference of outdoor space encroachment by buildings in blocks of Beijing Old Town; thus, this indicator was very applicable in this study area. Meanwhile, all height-related indicators, including MH, HBI and HFD, were negatively correlated with LST mean , which was consistent with the findings of most studies [30,[84][85][86][87]. This probably resulted from increasing and differentiating building height in one area that can lead to shadows provision, wind corridors enhancement, building land saving and blue or green spaces reservation, thus contributing to UHI mitigation. A comprehensive comparison of the absolute values of correlation coefficients showed that BCR and MH were the most important built-up factors influencing the LST mean . Since these two indicators are often used in urban planning, future consideration for buildings management should be mainly from these two aspects to make adaptive UHI mitigation measures.
This study tried to meet the needs of urban planning in terms of analysis process setting and indicators and research unit selection. It will be meaningful and suitable for urban planning, as it can offer practical, operational guidelines for urban managers. We therefore focused on the correlation between LST and the composition and configuration of GBI in blocks with different morphological characteristics, for GBI had been recognized as the nature-based solution (NbS) to mitigate the UHI effect in cities. It was reflected that both the GBI composition (trees/shrubs, water bodies, or grasslands) and the GBI patch size and aggregation (AREA_MN and AI) all contributed to reducing the LST mean , but partly contributed to increasing the LST std , and the significance of the coefficients differed in three types of blocks. Greater cooling potential can be provided by changing the patch area and aggregation degree among three types of blocks in Beijing Old Town, which is similar to the findings with Masoudi [13]. Besides, in the configuration section, the relationships were more consistent in HMB and LLB, and significantly different in MHB. The relatively consistent in terms of architectural structure between HMB and LLB might be one possible reason, as they are both dominated by hutongs and bungalows. It is also found that the absolute value of the correlation coefficient between GBI configuration and LST mean was higher in HMB and LLB, where low or medium-rise buildings dominate, than in MHB, where there are more high-rise buildings, suggesting that the contribution of building shadows on heat reduction may make the benefits of GBI on cooling less remarkable. Consequently, cooling strategies should be differentiated by block types.
There are also some aspects that can be further improved in our study. We applied the Pearson correlation analysis method, which is a simpler and commonly used method (which was applied in 68% of the relevant studies [78]). However, some studies pointed to the fact that the UHI phenomenon should be modelled locally instead of having an aggregated model for an entire area, as thermal environment characteristics are context sensitive, i.e., it varies significantly over space [63,88]. Although the block types were divided and conducted the correlation analysis separately, which avoided the impact of spatial context differences to some extent, some methods can provide more precise results of morphological indicators and LST; e.g., curve fit linear regression. We will take these methods into account in future studies [89]. Additionally, improving the ventilation environment has an important role in enhancing urban thermal environment, and the element of wind is also closely related to some important indicators in UP. Therefore, conducting the air path studies at block level by assessing current ventilation environment and proposing potential and existing ventilation corridors should be taken into consideration in further studies. Finally, in the description of urban 3D morphology, we used spatial structure indicators of buildings. However, as vegetated areas are an important composition of urban landscape, different plant types also have differences in 3D characteristics, which may have an impact on thermal environment. In the future, we can combine high resolution remote sensing data with LiDAR data to describe urban 3D morphology characteristics more comprehensively.

Conclusions
As the UHI effect and its associated consequences are expected to be more severe, mitigation of UHI has attracted increasing attention in urban studies. We followed a common but practice-oriented framework of UHI research to guide thermal mitigation solutions in this study. The spatial heterogeneity of UHI in blocks of Beijing Old Town (BOT) are analyzed, with an emphasis on the relationship between thermal environment and urban morphological characteristics, including basic land cover proportions and spatial structure indicators of built-up, as well as landscape composition and configuration of green and blue infrastructure within different spatial formed blocks.
Results showed that an obvious UHI effect was found in BOT, thermal environment within was block heterogenous and significant hot and cold areas with significant high or low land surface temperature, respectively, were all found, both of which were highly aggregated spatially. Proportions of vegetated area and water body were found to have a significant negative impact on block LST, while the proportion of impervious surface (ISP) was the opposite, with the highest absolute value of correlation coefficient among land cover characteristic indicators, reflecting its dominance in LST influencing. Bare soil area, on the other side, the only land cover type of which proportion was not significant with LST. By excluding the effect from land cover indicators, we revealed the significant correlation between LST and several spatial structure indicators, i.e., building coverage ratio, mean height, highest building index, height fluctuation degree, space crowd degree and sky view factor. It also showed that coverage ratio and the mean height of buildings were the two primary indicators having the greatest impact on LST in BOT through the comparation of absolute values; thus, more attention can be paid to those two aspects in further urban renewal as they are commonly used in urban planning.
The above laid a certain foundation for classifying block types in BOT by their morphological characteristics, after which we analyzed the impact of a comprehensive selection of indicators characterizing the composition and configuration of GBI in different block types. It showed that the effectiveness of different indicators for guiding the construction of GBI in each block varies, but generally by optimizing blue-green space composition, increasing patch size, aggregation, dominance and connectivity, all therefore have potential for promoting a block thermal environment. It is hoped that adaptive GBI planning and design for block thermal mitigation could be conducted based on this correlation study in the future renewal of old downtown.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Category Indicators Formula Note
Impervious surface area (ISA) In this paper, ISA of urban area was extracted by linear spectral hybrid image element decomposition model. It mainly includes the minimum noise separation, pure image element processing, end element collection, linear spectral separation, result checking and correction of the pre-processed images [73,74] Urban spatial structure characteristics Building coverage ratio (BCR) Hi , Fi refer to the height, volume, footprint and perimeter of the building No.i respectively, n is the number of buildings, A is the area of the block, C = 3.0 m is a constant, H max is the maximum height of the buildings in the block, H min is the minimum height of the buildings in the block [71].  Aggregation index (AI) AI = g ii max−g ii g ii is the number of similar neighboring patches of GBI in a block.
Note: Methematical formulas of above indicators are refering to McGarigal and Marks [67].