Assessing the Distribution of Urban Green Spaces and its Anisotropic Cooling Distance on Urban Heat Island Pattern in Baotou, China

An essential part of urban natural systems, urban green spaces play a crucial role in mitigating the urban heat island effect (UHI). The UHI effect refers to the phenomenon where the temperature within a city is higher than that of the surrounding rural areas. The effects of the spatial composition and configuration of urban green spaces on urban land surface temperature (LST) have recently been documented. However, few studies have examined the effects of the directionality and distribution of green spaces on LST. In this study, we used a landscape index to describe the change in pattern of heat island intensity for the city of Baotou, China. We then used a semi-variable function and nearest neighbor algorithm to analyze the cooling effects of green spaces. We found that: (1) the cooling distance of an urban green space was not only influenced by its size, vegetation cover, and shape, but also showed anisotropy. In general, the larger the area of the urban green space and the higher the value of Normalized Difference Vegetation Index (NDVI; a measure of plant photosynthetic activity), the larger the cooling distance within a certain threshold. Green spaces with more regular shapes displayed higher LST mitigation; however, the cooling distance was directional, and cooling effects depended on the semi-major axis and semi-minor axis of the green space. (2) The distribution of the urban green space within the landscape played a key role in mitigating the UHI effect. Within a certain area, the cooling effect of green spaces that are evenly distributed was greater than that which was associated with either green spaces that were large in area or where greens spaces were aggregated in the landscape. Therefore, within urban areas, where space is limited, urban planning should account for green spaces that are relatively scattered and evenly distributed to maximize cooling effects. The results of this study have key implications for sustainable urban planning and development; to mitigate urban heat island effects it is important to not only increase canopy cover or the size of urban green spaces, but also to optimize their spatial configuration.


Introduction
Urban areas are becoming increasingly important and the United Nations projects that more than 65% of the global population will be urban dwellers by the year 2050 [1]. It is therefore essential to work towards a better understanding of the complex processes that are involved with

Study Area
Baotou, the largest city of the Inner Mongolia Autonomous Region and an important industrial center of China, is located within western Inner Mongolia between 109˝25 1 E and 110˝25 1 E longitude and 40˝25 1 N and 40˝50 1 N latitude (Figure 1). The administrative area of Baotou covers 11 districts and a total area of 2.776ˆ10 10 m 2 , and has a metropolitan area covering 2.2ˆ10 8 m 2 , a population of 2,799,200 and includes four districts; Kundulun, Qingshan, Donghe, and Jiuyuan. Baotou has a temperate continental climate with a mean annual temperature of 6.4˝C, an average total annual rainfall of 310 mm, and a total of 110-143 frost-free days.

Study Area
Baotou, the largest city of the Inner Mongolia Autonomous Region and an important industrial center of China, is located within western Inner Mongolia between 109°25′E and 110°25′E longitude and 40°25′N and 40°50′N latitude (Figure 1).The administrative area of Baotou covers 11 districts and a total area of 2.776 × 10 10 m 2 , and has a metropolitan area covering 2.2 × 10 8 m 2 , a population of 2,799,200 and includes four districts; Kundulun, Qingshan, Donghe, and Jiuyuan. Baotou has a temperate continental climate with a mean annual temperature of 6.4 °C, an average total annual rainfall of 310 mm, and a total of 110-143 frost-free days.

Data and Processing
We used 30 m resolution remote sensing images from Landsat 5 TM taken in 2000, 2004, 2007, and 2011, and Landsat 8 OLI from 2014. We downloaded these images, which were captured during the growing season (July through August), from the US Geological Survey website. Images were only selected if they had cloud cover of 2% or below. We conducted radiometric calibration, atmospheric correction, geometric correction (single-point precision below 1 pixel), masking and other treatments on image data from all five time periods using ENVI5.1 and ArcGIS10.3. We then extracted urban green spaces found within Baotou by screen visual interpretation. We conducted an accuracy assessment of our urban green space classification using ground truthing from an independent field study, and our classification accuracy was 95% or higher.

Land Surface Temperature (LST)
We adopted the mono-window algorithm to retrieve the LST with a spatial resolution of 30 m from Landsat TM/OLI TIR channel [34]. The formulas are as follows:         6 6 6 6 6 6 6 6 6 6 6

Data and Processing
We used 30 m resolution remote sensing images from Landsat 5 TM taken in 2000, 2004, 2007, and 2011, and Landsat 8 OLI from 2014. We downloaded these images, which were captured during the growing season (July through August), from the US Geological Survey website. Images were only selected if they had cloud cover of 2% or below. We conducted radiometric calibration, atmospheric correction, geometric correction (single-point precision below 1 pixel), masking and other treatments on image data from all five time periods using ENVI5.1 and ArcGIS10.3. We then extracted urban green spaces found within Baotou by screen visual interpretation. We conducted an accuracy assessment of our urban green space classification using ground truthing from an independent field study, and our classification accuracy was 95% or higher.

Land Surface Temperature (LST)
We adopted the mono-window algorithm to retrieve the LST with a spatial resolution of 30 m from Landsat TM/OLI TIR channel [34]. The formulas are as follows: T s " ra 6 p1´C 6´D6 q`rb 6 p1´C 6´D6 q`C 6`D6 s T 6´D6 T a s {C 6 (1) C 6 " ε 6 τ 6 (3) where, T s is the real ground temperature (K); T 6 is brightness temperature (K); L λ is radiance (mW/(cm 2¨s r¨µm)), K 1 and K 2 are constants (mW/(cm 2¨s r¨µm) and K) obtained from the Landsat image header; ε 6 , τ 6 , and T a are surface emissivity, atmospheric transmittance and average atmospheric temperature (K), respectively. The surfaces were mainly divided into three types, i.e., water, urban, and natural surfaces. The ε 6 of water was taken as 0.995 and the urban and natural surface were estimated from vegetation coverage. The τ 6 was calculated using the LOWTRAN7 atmospheric model, and T a was calculated using the standard atmosphere section average temperature [34]: a 6 =´67.355, b 6 = 0.458 (temperature range of 0-70˝C).

Spatial Pattern of UHI Intensity
To account for variability in the temperature obtained from remote sensing images over several time periods, we used the following Equation (5) to calculate the relative temperature (˝C) and then divided the UHIs intensity into different grades using this threshold value (Table 1): where, T ui is brightness of the i th study area, and T is the average ground temperature of Baotou. Landscape metrics were also selected to describe the pattern of UHI intensity within each grade, such as Class Area (CA), Number of Patches (NP), Landscape Shape Index (LSI), Aggregation Index (AI), Shannon's Evenness Index (SHEI), Mean Patch Fractal Dimension (FRAC_MN), and Perimeter-Area Ratio (PARA). These indices were all calculated using the software Fragstats 4.2. The CA is equal to the sum of the area of all patches with the corresponding patch type (m 2 ) and the NP refers to the number of patches of the same patch type. The LSI provides a standardized measure of total edge or edge density of a patch while adjusting for the size of the overall landscape, and, when LSI > 1, the higher the value the larger the degree of landscape fragmentation. The AI is a measure of the relative degree of aggregation of patches in the landscape, ranging from 0 to 100 %, in which a greater value represents a higher degree of patches that are clumped. The SHEI considers an even distribution of area among patch types to be a landscape with maximum evenness; when the value is 0, it means that there is only one kind of patch and when the value is 1, it suggests the uniform distribution of different patch types. The FRAC_MN ranges from 1 to 2, and provides a measure of patch size and the degree of convolution of its border, in which the greater the value, the more complex the shape. The PARA is a simple measure of shape complexity, in which the greater the value, the more complex the shape. SPSS 20.0 statistical software was used to analyze the relationship between the isotropic cooling distance and urban green space area, and between the mean NDVI, with a spatial resolution of 30 m, and the shape index of green spaces (PARA).
Six typical sampling areas were also selected in urban green spaces which were defined as aggregates. We generated ten equal interval (30 m interval) buffers, from 0 m-300 m around the urban green space aggregates, and we calculated the mean surface temperature of each for the five time periods of Landsat images.

Cooling Distance and Direction of Urban Green Space
We detected urban green space cooling distance and direction using a semi-variance function [35], which is often used when there is a spatially-correlated distance or directional bias in the data. Our data are inherently spatially biased, which makes them appropriate for this function. The semi-variance function is defined as: where, γ(h) is the semi-variance, Z(x) is the regional variable, N(h) is of the sample size, and h is the separation distance. The spatial variation of the regional variable shows a stable value as the distance increases to a certain degree; that is, a partial sill. This distance is the maximum spatially-correlated interval, which, for this study, is the maximum cooling distance. The semi-variance function of the two-dimensional regional variable is also directional; when the ratio of two different directions is equal or close to 1, it shows isotropy, otherwise it will show anisotropy. The cooling distances of the different directions were derived from the semi-variance model parameters of ArcGIS10.3 geostatistical analyst toolkit.
Since the semi-variance function cannot be applied to aggregated urban green spaces (green space patches within a distance of less than 150 m from each other) with mutual cooling effects, we generated ten buffer rings of equidistant division from 0 m-300 m at 30 m intervals (corresponding to the resolution of the raster image) around the selected aggregated urban green spaces and calculated the average ground temperature of each interval. The point along the 0-300 m buffer distance at which the temperature semi-variance curve became stable, is approximately the cooling distance of that aggregate of urban green spaces.

Mitigating Effect of Green Space Spatial Distribution
To discover the relationship between the spatial distribution of urban green spaces and the heat island pattern, we used the k-nearest neighbor algorithm. This algorithm calculates a nearest neighbor index based on the mean distance from each feature to its nearest neighboring feature. Next, it identifies whether the distribution pattern of urban green spaces is aggregated or scattered. We calculated the k-nearest neighbor values using the Spatial Statistics Tools in ArcGIS10.3 and divided the nearest neighbor indices into different grades using the threshold value (Table 2). In addition, we calculated the percentage of area for each urban green space that is within each district to determine the influence of district on average temperature. We obtained urban green space aggregation degree and the area proportion of four districts in Baotou, including Kundulun, Qingshan, Donghe, and Jiuyuan, and the established the relationships between the aggregation degree, area proportion, and the mean temperature in each district. Table 2. k-nearest neighbor index classification.

Nearest Neighbor Index
Type More aggregated than random distribution 1.0 Random distribution 1.0-1. 5 More scattered than random distribution >1.5 Scattered

Spatiotemporal Variations of UHI Intensity
From the temperature retrieval and heat island intensity classification (Equation (5) and Table 1), we found that the heat island intensities of Baotou were mostly classified as moderate, strong, or very strong UHIs. These moderate to very strong UHIs were distributed in the downtown area, specifically in areas with relatively high building density or on unused land. In contrast, urban green spaces and their surrounding areas were classified as either no, or as weak heat islands. Extremely strong heat islands only accounted for a relatively small area, mainly located in the industrial area of Baotou ( Figure 2). spaces and their surrounding areas were classified as either no, or as weak heat islands. Extremely strong heat islands only accounted for a relatively small area, mainly located in the industrial area of Baotou ( Figure 2).

Figure 2.
Intensity map of heat islands in Baotou. The UHI intensity was classified into six categories, with the colors from blue to red representing the intensity of UHIs from weak to strong.
Seven landscape indices were selected to measure UHI intensity, targeting four types of heat islands; moderate, strong, very strong, and extremely strong, reflecting their dynamic change within the five time periods (Figure 3). We found that the areas of moderate, strong, and very strong heat islands together accounted for a relatively large proportion of the total area (each accounting for about 25%). Areas of moderate and strong heat islands increased (from 13% and 21.7% to 20.9% and 35.7%, respectively, Figure 3a,b) before 2011, and then began to decrease (to 19.6% and 21.8%, respectively) by 2014. There was an initial increase in the number of moderate and strong heat island patches from 2000 to 2007, followed by a decrease in 2011, with the number of strong heat island patches increasing again in 2014. The areas of very strong and extremely strong heat islands first decreased (from 37.2% and 1.8% to 23.4% and 0.8%, respectively) and then increased (to 27.7% and 2.7%, respectively). However, the number of patches within the very strong and extremely strong heat island classifications showed little change over the study period. The total variation of heat island landscape shape indices were not obvious (Figure 3c), but the indices of moderate and strong heat islands were relatively high and diffuse in distribution, which is supported by the AI values ( Figure 3e). However, overall, the pattern across all categories was considered to be aggregated in distribution, especially for the very strong heat island classification. There was a gradual decreasing trend (Figure 3d), ranging from 1.004 to 1.092, on the mean patch fractal dimensions of the four heat island types, shown by the total FRAC_MN value for Baotou (Figure 3f). This indicates that these Seven landscape indices were selected to measure UHI intensity, targeting four types of heat islands; moderate, strong, very strong, and extremely strong, reflecting their dynamic change within the five time periods (Figure 3). We found that the areas of moderate, strong, and very strong heat islands together accounted for a relatively large proportion of the total area (each accounting for about 25%). Areas of moderate and strong heat islands increased (from 13% and 21.7% to 20.9% and 35.7%, respectively, Figure 3a,b) before 2011, and then began to decrease (to 19.6% and 21.8%, respectively) by 2014. There was an initial increase in the number of moderate and strong heat island patches from 2000 to 2007, followed by a decrease in 2011, with the number of strong heat island patches increasing again in 2014. The areas of very strong and extremely strong heat islands first decreased (from 37.2% and 1.8% to 23.4% and 0.8%, respectively) and then increased (to 27.7% and 2.7%, respectively). However, the number of patches within the very strong and extremely strong heat island classifications showed little change over the study period. The total variation of heat island landscape shape indices were not obvious (Figure 3c), but the indices of moderate and strong heat islands were relatively high and diffuse in distribution, which is supported by the AI values ( Figure 3e). However, overall, the pattern across all categories was considered to be aggregated in distribution, especially for the very strong heat island classification. There was a gradual decreasing trend (Figure 3d), ranging from 1.004 to 1.092, on the mean patch fractal dimensions of the four heat island types, shown by the total FRAC_MN value for Baotou (Figure 3f). This indicates that these types of heat islands are relatively evenly-spaced in Baotou, and across all types, represented 84%-85% in the aggregation model (Figure 3g,h). types of heat islands are relatively evenly-spaced in Baotou, and across all types, represented 84%-85% in the aggregation model (Figure 3g,h).

Urban Green Spaces Cooling Distance
Nine independently-distributed urban green spaces, whose surrounding temperatures had not been affected by adjacent green spaces, were selected for isotropic and anisotropic cooling distance analysis within the five periods of our remote sensing image data. There was a logarithmic relationship between both area and shape index with cooling distance (R 2 = 0.6 for both relationships) (Figure 4a,c), and a linear relationship between mean NDVI and cooling distance (a significance of 10%; Figure 4b).
Urban green spaces with more area corresponded to a larger cooling distance; however, the influence distance tended to be stable with the increase of area. The shape of the urban green space patches also affected cooling distance, and cooling distance decreased with shape complexity. The internal cooling effect of high NDVI of urban green space was significant, and presented a positive correlation with cooling distance on the whole. However, this effect was not obvious because of the influence of other factors.

Urban Green Spaces Cooling Distance
Nine independently-distributed urban green spaces, whose surrounding temperatures had not been affected by adjacent green spaces, were selected for isotropic and anisotropic cooling distance analysis within the five periods of our remote sensing image data. There was a logarithmic relationship between both area and shape index with cooling distance (R 2 = 0.6 for both relationships) (Figure 4a,c), and a linear relationship between mean NDVI and cooling distance (a significance of 10%; Figure 4b). Urban green spaces with more area corresponded to a larger cooling distance; however, the influence distance tended to be stable with the increase of area. The shape of the urban green space patches also affected cooling distance, and cooling distance decreased with shape complexity. The internal cooling effect of high NDVI of urban green space was significant, and presented a positive correlation with cooling distance on the whole. However, this effect was not obvious because of the influence of other factors.  Analysis of the urban green space aggregates identified that the surrounding temperature of the six selected sampling areas increased with distance and displayed a law of gradient variation ( Figure  5). The variation tended to stable out around 180 m, before which the cooling magnitude was higher and after which the cooling magnitude was lower (<1 °C). Within 300 m, the cooling temperature was between 1.9 and 3.1 °C, and the maximum cooling distance was between 120 and 300 m. . Relationships between urban green cooling distance, area, NDVI mean, and PARA. (a) represents the relationship between urban green space area and cooling distance in five time periods remote sensing image data; (b) represents the relationship between the mean NDVI of each urban green space and cooling distance in the same images; and (c) represents the relationship between shape index of each urban green space and the cooling distance in the same images.
Analysis of the urban green space aggregates identified that the surrounding temperature of the six selected sampling areas increased with distance and displayed a law of gradient variation ( Figure 5). The variation tended to stable out around 180 m, before which the cooling magnitude was higher and after which the cooling magnitude was lower (<1˝C). Within 300 m, the cooling temperature was between 1.9 and 3.1˝C, and the maximum cooling distance was between 120 and 300 m.

The Directivity of Urban Green Spaces Cooling Distance
The anisotropic cooling distance of the selected urban green space (Table 3) was associated with green space shape. The maximum cooling distance was on the semi-major axis of the green space, while the minimum distance was on the semi-minor axis. The cooling distance from the centroid of green space increased with green space area. However, there was no clear trend to follow for the cooling distance away from the green space border. The maximum directional cooling distance was between 150 and 454 m, while the minimum directional cooling distance was between 106 and 333 m.

Cooling Effect of Urban Green Space Distribution
We found that urban green space distribution in Baotou was generally random or more scattered than random, and that green area accounted for 10%-40% of each district ( Figure 6). Additionally, there was no significant correlation between the degree of aggregation of urban green spaces and mean temperature; however, the mean temperature did decrease as the urban green space distribution became dispersed. The ratio of urban green space area in each district presented a negative effect on mean temperature, and mean temperature tended to be lower with an increase of green space to area ratio. If the area ratio was greater than 50%, a more significant cooling effect could be identified.

The Directivity of Urban Green Spaces Cooling Distance
The anisotropic cooling distance of the selected urban green space (Table 3) was associated with green space shape. The maximum cooling distance was on the semi-major axis of the green space, while the minimum distance was on the semi-minor axis. The cooling distance from the centroid of green space increased with green space area. However, there was no clear trend to follow for the cooling distance away from the green space border. The maximum directional cooling distance was between 150 and 454 m, while the minimum directional cooling distance was between 106 and 333 m.

Cooling Effect of Urban Green Space Distribution
We found that urban green space distribution in Baotou was generally random or more scattered than random, and that green area accounted for 10%-40% of each district ( Figure 6). Additionally, there was no significant correlation between the degree of aggregation of urban green spaces and mean temperature; however, the mean temperature did decrease as the urban green space distribution became dispersed. The ratio of urban green space area in each district presented a negative effect on mean temperature, and mean temperature tended to be lower with an increase of green space to area ratio. If the area ratio was greater than 50%, a more significant cooling effect could be identified.  Figure 6. Cooling effect diagram of urban green space distribution. (a) Represents the relationship between aggregation of urban green space of each district and mean temperature of each district in five time periods of remote sensing imagery data. When the value is 1, the urban green space distribution is random, when the value is less than 1, the distribution is aggregated, and when the value is more than 1 the distribution is scattered; and (b) represents the relationship between the proportion of area of urban green space in each district and mean temperature of each district in the same images.

Relationship between the Urban Heat Island and Green Space
We found that areas with high building density have more UHI, in contrast to urban green spaces and their surrounding areas where no, or only weak, heat islands were found. The amount of area of urban green space was a strong predictor for mitigating UHI effects, suggesting that larger green spaces have larger cooling distances [13,32,36]. However, a strong positive logarithmic relationship between area and the cooling distance (R 2 = 0.62) was also found, which suggested that with the increase of urban green space area the cooling effect plateaus. Mean NDVI had a weak linear correlation with cooling distance (R 2 = 0.07). Diverse land cover across the urban green space was likely the main cause for this, i.e., species composition [17], community structure [18], internal water bodies, and permeable surface [19,20] all likely contribute to the cooling effect. Accumulative or area weighted index, such as NDVI or percent cover of woody vegetation [27], might be a good candidate for analysis of green space mitigating effects of the urban heat island effects in further studies. We also found that the spatial configuration of green space had a significant effect on UHI. This result is consistent with findings from previous studies [36], in which it was reported that the more regularly spaced the urban green spaces are, the lower the surface temperatures. However, Zhou et al. (2011) suggested that LST decreased through the increase in the amount of green space edge, which enhances the thermal exchange between the green space and its surrounding areas [27]. We suggest the relatively small size of green spaces in this study was insufficient to identify significant cooling effects [37].
The directivity of urban green spaces cooling distance is anisotropic, and does have an effect on LST pattern. The maximum cooling distances (150 to 454 m) and minimum cooling distances (106 to 333 m), respectively, corresponded to the directions of the radius of the urban green spaces. The cooling effects of urban green spaces were related to their shape; green spaces in simple shapes showed isotropy while those with complex shapes showed the contrary. Interestingly, the different directions did not have the same cooling distances. However, the semi-variance function was only suitable for analyzing individual green spaces rather than concentrated spaces, which could influence the result. Due to the limited capacity of the function itself, and because the surface temperature of the surrounding areas was affected by the similarity of aggregated green spaces with their surrounding area, the cooling distance was enlarged. For this reason, we used the buffer zone method to calculate the cooling effect of different urban green spaces in concentrated areas. We found that Figure 6. Cooling effect diagram of urban green space distribution. (a) Represents the relationship between aggregation of urban green space of each district and mean temperature of each district in five time periods of remote sensing imagery data. When the value is 1, the urban green space distribution is random, when the value is less than 1, the distribution is aggregated, and when the value is more than 1 the distribution is scattered; and (b) represents the relationship between the proportion of area of urban green space in each district and mean temperature of each district in the same images.

Relationship between the Urban Heat Island and Green Space
We found that areas with high building density have more UHI, in contrast to urban green spaces and their surrounding areas where no, or only weak, heat islands were found. The amount of area of urban green space was a strong predictor for mitigating UHI effects, suggesting that larger green spaces have larger cooling distances [13,32,36]. However, a strong positive logarithmic relationship between area and the cooling distance (R 2 = 0.62) was also found, which suggested that with the increase of urban green space area the cooling effect plateaus. Mean NDVI had a weak linear correlation with cooling distance (R 2 = 0.07). Diverse land cover across the urban green space was likely the main cause for this, i.e., species composition [17], community structure [18], internal water bodies, and permeable surface [19,20] all likely contribute to the cooling effect. Accumulative or area weighted index, such as NDVI or percent cover of woody vegetation [27], might be a good candidate for analysis of green space mitigating effects of the urban heat island effects in further studies. We also found that the spatial configuration of green space had a significant effect on UHI. This result is consistent with findings from previous studies [36], in which it was reported that the more regularly spaced the urban green spaces are, the lower the surface temperatures. However, Zhou et al. (2011) suggested that LST decreased through the increase in the amount of green space edge, which enhances the thermal exchange between the green space and its surrounding areas [27]. We suggest the relatively small size of green spaces in this study was insufficient to identify significant cooling effects [37].
The directivity of urban green spaces cooling distance is anisotropic, and does have an effect on LST pattern. The maximum cooling distances (150 to 454 m) and minimum cooling distances (106 to 333 m), respectively, corresponded to the directions of the radius of the urban green spaces. The cooling effects of urban green spaces were related to their shape; green spaces in simple shapes showed isotropy while those with complex shapes showed the contrary. Interestingly, the different directions did not have the same cooling distances. However, the semi-variance function was only suitable for analyzing individual green spaces rather than concentrated spaces, which could influence the result. Due to the limited capacity of the function itself, and because the surface temperature of the surrounding areas was affected by the similarity of aggregated green spaces with their surrounding area, the cooling distance was enlarged. For this reason, we used the buffer zone method to calculate the cooling effect of different urban green spaces in concentrated areas. We found that the variation of temperature was higher before 180 m, after which the variation of temperature was lower (<1˝C), and the cooling temperature was between 1.9 and 3.1˝C.
The ground surface temperature, and its aggregation degree, showed a positive correlation with the ground surface temperature (Figure 6a). This meant that the temperature was lowered when green spaces were scattered. This result was contrary to the conclusion drawn by Wang et al. (2014), that there was a negative correlation between aggregation degree and urban ground surface temperature [38]. On one hand, it is related to the size of green spaces and the mean distance (or arrangement) between each of them. In Li et al. (2012), a similar result to Wang et al. (2014) was observed in Beijing, a city which is dominated by small patches [32]. Located near paved roads, small green spaces that are scattered across a city usually have a high thermal load and are unlikely to produce significant cool island effects [39], with the mean distance of those green spaces larger than its cooling effects. As a result, the cooling effect of aggregated urban green spaces is limited. However, in our study, the mean distance between each green space and average green space size is 151.2 m and 2.1ˆ10 5 m 2 respectively, providing more area and greater distances of urban green spaces to analyze. On the other hand, the relationship between mean temperature and the green spaces area ratio of each district showed negative correlation (Figure 6b). When the ratio was greater than 50%, the influence might become clearer. Thus, during the planning of urban green spaces, we should consider their spatial distribution and evenly distribute them as much as possible. We feel that a rational design of urban green space area would not only mitigate the UHI effect, but also play a role of leisure and recreation, earthquake preparedness, and disaster reduction, etc.

Implications for Urban Planning
Vegetation management, such as planting more trees, should be considered an effective approach to cooling the heated city [27]. This has already been done for Beijing, where they planned the wedge green belt to increase tree canopy cover to cool down the metropolitan area [40]. The size of an urban green space is an important factor that effect UHI dynamics that is widely accepted. However, it is infeasible to account for the building expense that green spaces come with in some urban area. We found that not only spatial configuration but also distribution of green space affected LST, which also supported by other researches [27,32,36]. When public space is limited, different spatial arrangements of land cover features can significantly decrease LST [27], e.g., encouraging private "greening" areas rather than building new expensive small-medium sized parks [39]; increasing patch density of the green space, rather than expanding the green space cover [36]; and covering roofs or facades with plants mitigating urban thermal environment [39,41], are all feasible options.

Conclusions
We analyzed five time periods of remote sensing image data to extract surface temperature and urban green spaces of built-up areas in Baotou using the mono-window algorithm. We also analyzed five time periods of UHI intensity change in Baotou with the landscape index, based on the effect of urban green spaces on UHI pattern using the methods of semi-variance function, buffer zone method, and aggregation degree. The conclusions as follow: (1) The cooling distance of an urban green space varied with its size and shape, and the cooling distance showed anisotropy. Within a certain threshold, the greater urban green space area and mean NDVI clearly corresponded to a larger cooling distance. With increasing of shape complexity of the green space, the cooling distance declined in certain green space areas. According to the semi-variance function, the cooling distance was directional, the maximum cooling distance was on the semi-major axis of the green space, and the cooling distance was between 150 and 454 m, while the minimum distance was on the semi-minor axis and the cooling distance was between 106 and 333 m.
(2) The distribution of urban green space plays a key role in mitigating the UHI effect. Within a certain space, in terms of the cooling effect, evenly-distributed green spaces have greater effects than aggregated distributed green spaces. Therefore, within the limited urban space, relatively scattered and evenly-distributed green spaces are a potential solution for urban planning to mitigate the UHI effect.