Numerical Simulation of Local Climate Zone Cooling Achieved through Modiﬁcation of Trees, Albedo and Green Roofs—A Case Study of Changsha, China

: By exploring the cooling potential of tree quantity, ground albedo, green roofs and their combinations in local climate zone (LCZ)-4, LCZ-5, and LCZ-6, this study focuses on the optimum cooling level that can be achieved in open residential regions in Changsha. It designs and models 39 scenarios by integrating in situ measurement and ENVI-met numerical simulation and further compares cooling e ﬀ ects of various combinations of the cooling factors. The results show that (1) an increased number of trees and higher albedo are more e ﬀ ective compared to green roofs in reducing summer potential temperatures at street level (2 m high) in three LCZs. Negative correlations are observed in the pedestrian air temperature with trees and ground albedo; (2) the e ﬀ ects of cooling factors vary among di ﬀ erent LCZ classes, with the increased 60% more trees leading to lower outdoor temperatures for LCZ-4 (0.28 ◦ C), LCZ-5 (0.39 ◦ C), and LCZ-6 (0.54 ◦ C), while higher albedo of asphalt surface (increased by 0.4) is more e ﬀ ective in LCZ-4 (reaches to 0.68 ◦ C) 14:00, compare to LCZ-5 (0.49 ◦ C) and LCZ-6 (0.38 ◦ C); (3) applying combined cooling methods can provoke air temperature reduction (up to 0.96 ◦ C), especially when higher levels of tree quantities (increased by 60%) are coupled with cool ground materials (albedo increased by 0.4). The results can contribute useful information for improving thermal environment in existing residential regions and future residential planning.


Introduction
Rapid urbanization leads to the transformation of vegetation area, water bodies, and cultivated land into urban buildings and impervious ground. These transformations reduce the evapotranspiration of vegetation and increase the absorption of solar radiation by impervious materials, resulting in higher air and surface temperatures in urbanized city center areas than in surrounding rural regions, thus forming the urban heat island (UHI) phenomenon [1]. Extreme heat events (EHE) [2] caused by the UHI not only lead to climatological problems [3] but also have an impact on public health [4][5][6] and even contribute to mortality among the urban population [7][8][9]. Recently, urban planning strategies have gained importance when UHI mitigation and sustainable development of the cities are considered [10]. In fact, a number of studies have been conducted on how the land surface properties (vegetation, material, and buildings) contribute to the urban microclimate conditions [11][12][13][14][15][16][17] as well as how the optimized urban physical properties facilitate energy reduction [18][19][20] to maintain the comfortable urban thermal environment and make the city more livable and sustainable.
Among the many UHI mitigation strategies, greenery has been regarded as one of the valid measures to have significant cooling effects. By increasing the quantity of vegetation [21,22] or

Target Region
Changsha is the capital city of Hunan Province, which is located in the southeast of China. The regional range of Changsha is 111 • 53 -114 • 15 E and 27 • 51 -28 • 41 N. Changsha City has a subtropical monsoon humid climate, which is extremely hot in summer and cold in winter. The dominant wind is southeast in summer and northwest in winter. From late May, there are 85 days when the daily average temperature is above 30 • C in summer, and 30 days when the temperature is above 35 • C [63].
Three residential areas were selected as the research object in Changsha, which are an open arrangement of high-rise buildings (LCZ-4), mid-rise buildings (LCZ-5), and low-rise buildings (LCZ-6) ( Figure 1). In order to avoid the interference brought by geographical location, the research objects are located in the adjacent blocks of Yuelu District of Changsha. However, the area of the three residential areas occupied are quite different, which may influence the comparison results of the three LCZs. In order to cover the same area in each LCZ, an area of 120 × 120 m of each residential area is defined (Figure 2), which means the areas outside the simulation domain also feature similar surface Sustainability 2020, 12, 2752 4 of 23 properties to the simulated areas. The three classes contain the characteristic of open LCZs, which have an abundance of pervious land cover (low plants, scattered trees), as well as concrete, steel, stone, and glass construction materials [54]. Table 1 shows the comparison of geometric and surface cover properties' values between standard statistics and the study cases. Although three figures do not align with the surface property values of LCZ-4 and LCZ-6, the selected three sites are identified and described with the best-fit LCZ class, which is allowed in the LCZ classification [54]. In order to avoid the interference brought by geographical location, the research objects are located in the adjacent blocks of Yuelu District of Changsha.
Sustainability 2020, 12, x FOR PEER REVIEW 4 of 24 steel, stone, and glass construction materials [54]. Table 1 shows the comparison of geometric and surface cover properties' values between standard statistics and the study cases. Although three figures do not align with the surface property values of LCZ-4 and LCZ-6, the selected three sites are identified and described with the best-fit LCZ class, which is allowed in the LCZ classification [54]. In order to avoid the interference brought by geographical location, the research objects are located in the adjacent blocks of Yuelu District of Changsha. To further calculate the parameters, in the center of each site, the fisheye camera was used at 2 m to calculate the sky view factor, which indicates the ratio of the amount of sky hemisphere visible from ground level to that of an unobstructed hemisphere. The aspect ratio was calculated through mean height-to-width ratio of street canyons, while the building surface fraction is the ratio of building plan area to total plan area. Additionally, the impervious surface fraction is the ratio of impervious plan area (paved, rock) to total plan area, and the pervious surface fraction refers to the ratio of pervious plan area (bare soil, vegetation, water) to total plan area. steel, stone, and glass construction materials [54]. Table 1 shows the comparison of geometric and surface cover properties' values between standard statistics and the study cases. Although three figures do not align with the surface property values of LCZ-4 and LCZ-6, the selected three sites are identified and described with the best-fit LCZ class, which is allowed in the LCZ classification [54]. In order to avoid the interference brought by geographical location, the research objects are located in the adjacent blocks of Yuelu District of Changsha. To further calculate the parameters, in the center of each site, the fisheye camera was used at 2 m to calculate the sky view factor, which indicates the ratio of the amount of sky hemisphere visible from ground level to that of an unobstructed hemisphere. The aspect ratio was calculated through mean height-to-width ratio of street canyons, while the building surface fraction is the ratio of building plan area to total plan area. Additionally, the impervious surface fraction is the ratio of impervious plan area (paved, rock) to total plan area, and the pervious surface fraction refers to the ratio of pervious plan area (bare soil, vegetation, water) to total plan area. To further calculate the parameters, in the center of each site, the fisheye camera was used at 2 m to calculate the sky view factor, which indicates the ratio of the amount of sky hemisphere visible from ground level to that of an unobstructed hemisphere. The aspect ratio was calculated through mean height-to-width ratio of street canyons, while the building surface fraction is the ratio of building plan area to total plan area. Additionally, the impervious surface fraction is the ratio of impervious plan area (paved, rock) to total plan area, and the pervious surface fraction refers to the ratio of pervious plan area (bare soil, vegetation, water) to total plan area.

Model Settings
The simulation area was 120 × 120 m, with the number of horizontal grids being 49 × 49 (horizontal resolution is 2.5 m). A nested grid was set around the site, with ground type loam, which is a buffer zone that can be created outside the core area to solve problems with the model not working reliably at and near the borders [17]. On the basis of the report of Hunan Meteorological Bureau, China, in 2019 [64], there were 22 days of high temperature from 12 July to 20 August in Changsha City, and the highest average air temperature appeared on 1 August 2019. Therefore, the simulation started at 7:00 on 1 August 2019 and ended at 7:00 on 2 August 2019 for 24 hours, with the initial temperature 33.12 • C, the wind speed 2.2 m/s, the wind direction southeast 150 • , and the humidity 56.50% [63]. In the simulated models, nine receivers were arranged proximately every 30 m ( Figure 2).
The full forcing function developed by ENVI-met 4.0 was used to drive the simulation with external meteorological data as boundary conditions. The data come from the Changsha meteorological station (No. 57687, location 112 • .91 E, 28 • .21 N), which is approximately 3.38 kilometers from the block where the three sites are located (112 • .93 E, 28 • .24 N) ( Figure 1). The meteorological station of Changsha is one of the national surface stations that provides hourly data, including hourly observational values of such weather elements as temperature, pressure, relative humidity, moisture pressure, wind, and precipitation. Other than that, it also needs vegetation data, soil data, building data, and ground data. The plant height and crown width are estimated according to the actual measurement from the site. The leaf area density (LAD) is based on the measured values of Cinnamomum camphora and Photinia rubra. The LAD of grassland is set as the default value from ENVI-met. The thermal properties of ground materials are set according to "Code for thermal design of civil building (GB 50176-2016)" [65]. The initial soil temperature and humidity are set based on general values [66]. Table 2 lists the main input parameters of the ENVI-met simulation. The accuracy of ENVI-met was tested by the data measured at test point #5 in the open high-rise model, open middle-rise model, and open low-rise model in the summer. The air temperature was recorded by the HOBO Temperature/Humidity Data Logger (model H 08-032-08) at a height of 2 m above the ground at the #5 receiver location in each LCZ scenario ( Figure 3). Since humidity is not considered in this study, the data of humidity were not recorded. The HOBO data logger has a temperature operating range of −30 to 50 • C and a resolution of 0.02 • C at 21 • C. The on-site Ta measurement was applied on 22 July 2019, which was the day after Changsha Meteorological Observatory issued the first high-temperature yellow early warning signal, in order to test the accuracy of ENVI-met under the extreme hot weather condition. Therefore, the validation measurement duration was 24 hours from 10:00 on 22 July to 10:00 on 23 July 2019, and the output data interval was set at every 6 minutes. Since the first three hours of simulation in ENVI -met is unstable, it was discarded. Thus, the simulations were from 7:00 of 22 July to 10:00 of 23 July 2019. low-rise model in the summer. The air temperature was recorded by the HOBO Temperature/Humidity Data Logger (model H 08-032-08) at a height of 2 m above the ground at the #5 receiver location in each LCZ scenario ( Figure 3). Since humidity is not considered in this study, the data of humidity were not recorded. The HOBO data logger has a temperature operating range of −30 to 50 ℃ and a resolution of 0.02 ℃ at 21 ℃. The on-site Ta measurement was applied on July 22, 2019, which was the day after Changsha Meteorological Observatory issued the first high-temperature yellow early warning signal, in order to test the accuracy of ENVI-met under the extreme hot weather condition. Therefore, the validation measurement duration was 24 hours from 10:00 on July 22 to 10:00 on July 23, 2019, and the output data interval was set at every 6 minutes. Since the first three hours of simulation in ENVI -met is unstable, it was discarded. Thus, the simulations were from 7:00 of July 22 to 10:00 of July 23, 2019. The sensitivity of measured and simulated data was analyzed with correlation measures and difference measures. The correlation evaluation index (R 2 ) was used to quantify the correlation degree between the simulation value and the observation value, while the root mean square error (RMSE) evaluated the error between the two indexes. The RMSE expression is as follows,

Model validation
From Figure 4, the RMSE and R 2 of the two sets of data are 1.39 ℃ and 0.89, respectively, in LCZ-4. For LCZ-5, the RMSE is the lowest (0.67 ℃), and R 2 is 0.97, while for LCZ-6, RMSE is 1.44 ℃, and R 2 is 0.88. Similarities were found in the three scenarios, which are, from 12:00 to 18:00, the actual air temperature is significantly underestimated in the simulation results, while from 19:00 to 2:00, the simulated temperature data are higher than the actual measurement data. The sensitivity of measured and simulated data was analyzed with correlation measures and difference measures. The correlation evaluation index (R 2 ) was used to quantify the correlation degree between the simulation value and the observation value, while the root mean square error (RMSE) evaluated the error between the two indexes. The RMSE expression is as follows, (1)

Model Validation
From Figure 4, the RMSE and R 2 of the two sets of data are 1.39 • C and 0.89, respectively, in LCZ-4. For LCZ-5, the RMSE is the lowest (0.67 • C), and R 2 is 0.97, while for LCZ-6, RMSE is 1.44 • C, and R 2 is 0.88. Similarities were found in the three scenarios, which are, from 12:00 to 18:00, the actual air temperature is significantly underestimated in the simulation results, while from 19:00 to 2:00, the simulated temperature data are higher than the actual measurement data.
Previous studies have reported simulation errors, which could be caused by many reasons [60,67,68]. Firstly, the anthropogenic heat release is not taken into consideration [69,70], leading to underprediction of daytime Ta. Secondly, the ENVI-met simulation maintains a single wind direction for the entire simulation period, which introduces uncertainty into the model results. Overall, it can be concluded that ENVI-met simulations provide acceptable accuracy regarding outdoor air temperatures in this study.  Previous studies have reported simulation errors, which could be caused by many reasons [60,67,68]. Firstly, the anthropogenic heat release is not taken into consideration [69,70], leading to underprediction of daytime Ta. Secondly, the ENVI-met simulation maintains a single wind direction for the entire simulation period, which introduces uncertainty into the model results. Overall, it can be concluded that ENVI-met simulations provide acceptable accuracy regarding outdoor air temperatures in this study.

Vegetation and Material Variable Settings
The numerical simulations were conducted to examine the cooling potential of different modifications on urban microclimates by considering the original form of LCZ-4 (open high-rise buildings), LCZ-5 (open middle-rise buildings), and LCZ-6 (open low-rise buildings) as three Baseline scenarios. There are other 36 scenarios for this study, which are 12 cooling strategy scenarios applied for LCZ-4, LCZ-5, and LCZ-6 ( Table 3).
In the outdoor thermal environment, the ability of plant types to reduce the air temperature of pedestrian height is trees > lawn > shrub [29] according to the order of cooling effect from strongest to weakest, so for the tree factor, on the basis of the Baseline scenario, the study designed and simulated the situation of adding different numbers of trees (30% and 60%) ( Table 4). In Scenario A and Scenario B, trees were added along the road ( Figure 5). There are two reasons for this modification: First, in order to explore how to reduce the temperature with minimum modification to the current situation, it only increases the number of trees, and thereby the existing evenly spaced trees along the road are maintained. Second, as found in the simulation of all the Baseline scenarios at 14:00, the hard paved area presents higher temperature. Therefore, the study assumes that the arrangement of trees along the road is more effective in cooling. A green roof typically is a vegetative layer with a growing medium over a waterproofing membrane [71]. It usually has two distinctive types, intensive and extensive. Intensive green roofs come with very high maintenance costs and heavy structural support because they employ deep growth media and high growth media for trees and shrubs [72], while extensive green roofs are cost-effective and easy to maintain, with a thin layer of vegetation and shallow soil [41]. Therefore, this study only considers the extensive green roof, which is more widely used. In this case, the Sedum emarginatum Migo species was considered in the green roof modelling, due to its long life cycle and great drought tolerant ability [73]. Meanwhile, the height of the soil layer and vegetation layer were both set to 0.10 m. In Scenario C, 50% roofs were replaced with green roofs, while in Scenario D, 100% roofs were covered by green roofs. The layouts are shown in Figure 6.
Replacing urban surfaces with high-albedo materials is a common strategy to mitigate UHI [74]. In this study, asphalt is the most dominant material for the main streets in the Baseline scenario, which is also the ground material taking up large space in each LCZ. The parking lot grounds in LCZ-4 and LCZ-6, as well as the walkway in LCZ-4, are made of brick, and the basketball court in LCZ-4 is plastic, which only takes up a small percentage. Therefore, the parking lot and basketball court sections remain unchanged in the variable settings ( Figure 2). Since the current street is covered with asphalt, where solar radiation is absorbed due to the low albedo [50], the albedo of the asphalt pavement on the streets was increased by 0.2 and 0.4 respectively in Scenario E and Scenario F, that is, to change the value of Sustainability 2020, 12, 2752 9 of 23 reflection in the database manager in ENVI-met. Although the research shows that roof surfaces and building walls can also reflect radiation [14], this is not considered in this study.
In this study, asphalt is the most dominant material for the main streets in the Baseline scenario, which is also the ground material taking up large space in each LCZ. The parking lot grounds in LCZ-4 and LCZ-6, as well as the walkway in LCZ-4, are made of brick, and the basketball court in LCZ-4 is plastic, which only takes up a small percentage. Therefore, the parking lot and basketball court sections remain unchanged in the variable settings ( Figure 2). Since the current street is covered with asphalt, where solar radiation is absorbed due to the low albedo [50], the albedo of the asphalt pavement on the streets was increased by 0.2 and 0.4 respectively in Scenario E and Scenario F, that is, to change the value of reflection in the database manager in ENVI-met. Although the research shows that roof surfaces and building walls can also reflect radiation [14], this is not considered in this study.

Results
The results of the study focus on three hours, which are 8:00, 14:00, and 20:00. The main reason is that 8:00 is usually the morning rush hour for daily commuting, while the summer temperature 4 and LCZ-6, as well as the walkway in LCZ-4, are made of brick, and the basketball court in LCZ-4 is plastic, which only takes up a small percentage. Therefore, the parking lot and basketball court sections remain unchanged in the variable settings ( Figure 2). Since the current street is covered with asphalt, where solar radiation is absorbed due to the low albedo [50], the albedo of the asphalt pavement on the streets was increased by 0.2 and 0.4 respectively in Scenario E and Scenario F, that is, to change the value of reflection in the database manager in ENVI-met. Although the research shows that roof surfaces and building walls can also reflect radiation [14], this is not considered in this study.

Results
The results of the study focus on three hours, which are 8:00, 14:00, and 20:00. The main reason is that 8:00 is usually the morning rush hour for daily commuting, while the summer temperature

Results
The results of the study focus on three hours, which are 8:00, 14:00, and 20:00. The main reason is that 8:00 is usually the morning rush hour for daily commuting, while the summer temperature reaches the highest at 14:00. Meanwhile, the urban heat island effect is obvious during the first 2-3 h after the sun sets [60]; therefore, 20:00 is the chosen hour in the evening. In the following, the air temperature refers to the average value of nine receivers ( Figure 2).

Baseline Scenario
It can be seen from Figure 7 that during the morning of 8:00-10:00, the temperature of LCZ-4 appears to be the highest, while LCZ-5 and LCZ-6 are slightly lower. The three types of LCZ all reach their peaks at 14:00. During the high temperature hours (12:00-14:00), the average air temperature of LCZ-6 reaches the highest Ta, and the LCZ-5 air temperature remains the lowest. After reaching the peak temperature, there is a downward trend in the air temperature for all three LCZs until 22:00-23:00. During the nighttime (20:00 to 5:00), the air temperature of LCZ-4, LCZ-5, and LCZ-6 are substantially coincident. reaches the highest at 14:00. Meanwhile, the urban heat island effect is obvious during the first 2-3 hours after the sun sets [60]; therefore, 20:00 is the chosen hour in the evening. In the following, the air temperature refers to the average value of nine receivers ( Figure 2).  reaches the highest at 14:00. Meanwhile, the urban heat island effect is obvious during the first 2-3 hours after the sun sets [60]; therefore, 20:00 is the chosen hour in the evening. In the following, the air temperature refers to the average value of nine receivers (Figure 2). It can be seen from Figure 7 that during the morning of 8:00-10:00, the temperature of LCZ-4 appears to be the highest, while LCZ-5 and LCZ-6 are slightly lower. The three types of LCZ all reach their peaks at 14:00. During the high temperature hours (12:00-14:00), the average air temperature of LCZ-6 reaches the highest Ta, and the LCZ-5 air temperature remains the lowest. After reaching the peak temperature, there is a downward trend in the air temperature for all three LCZs until 22:00-23:00. During the nighttime (20:00 to 5:00), the air temperature of LCZ-4, LCZ-5, and LCZ-6 are substantially coincident. Figure 8 shows the air temperature distribution of the Baseline scenario at the selected three hours (8:00, 14:00, and 20:00). At 14:00, the northwest parts of LCZ-4 are cooler than the rest of the areas due to geometrical shading of streets by buildings. In LCZ-4, the street areas between buildings and the open areas in the southeast have higher temperatures, which can be explained by the low ground albedo materials and the lack of greenery within those areas. The low Ta areas of LCZ-5 (open middle-rise) are mainly concentrated in interspaced areas between buildings on the east. In contrast, the Ta layout of LCZ-6 at 14:00 indicates the asphalt has higher temperature compared with other pavement materials, mainly because the low buildings in LCZ-6 cannot create large areas of shadows, causing the air temperature increase at 14:00.

Tree quantity
All LCZs' Ta values are reduced after the addition of trees. A 30% increase of tree covers (Scenario A) is sufficient to cut the incoming solar radiation and reduce average air temperature by 0.10-0.30 ℃, while the 60% increase (Scenario B) shows an approximate drop of 0.10-0.60 ℃ during the day, indicating that more trees provide more effective cooling. However, the temperature drop is insignificant (0.08-0.18 ℃) for every 30% tree increase by 20:00. This may due to an increase in tree block long wave radiation escaping to the night sky, which prevents the function of nocturnal cooling [75].
According to Figure 9, the greatest reductions are achieved along the street in LCZ-4, LCZ-5, and LCZ-6. LCZ-4 and LCZ-5 have relatively tall buildings, which are already well-shaded most of the day, which causes the increased tree quantity to be less effective. It is noted in Figure 9 that the Ta reduction around buildings is lower than in other areas because these zones are already cooler than the street area in the Baseline scenario due to geometrical shading by buildings. In LCZ-4, even though trees that are along the streets keep an obvious cooling effect at the peak hour 14:00 ( Figure  9a), the impact is not significant enough to reduce the temperature of the entire site, especially in the areas that are not covered by tree shade and building shadows. For example, the temperature at measurement point #9 increases by 0.06 ℃ at 14:00, which is attributed to the direct exposure to

Tree Quantity
All LCZs' Ta values are reduced after the addition of trees. A 30% increase of tree covers (Scenario A) is sufficient to cut the incoming solar radiation and reduce average air temperature by 0.10-0.30 • C, while the 60% increase (Scenario B) shows an approximate drop of 0.10-0.60 • C during the day, indicating that more trees provide more effective cooling. However, the temperature drop is insignificant (0.08-0.18 • C) for every 30% tree increase by 20:00. This may due to an increase in tree block long wave radiation escaping to the night sky, which prevents the function of nocturnal cooling [75].
According to Figure 9, the greatest reductions are achieved along the street in LCZ-4, LCZ-5, and LCZ-6. LCZ-4 and LCZ-5 have relatively tall buildings, which are already well-shaded most of the day, which causes the increased tree quantity to be less effective. It is noted in Figure 9 that the Ta reduction around buildings is lower than in other areas because these zones are already cooler than the street area in the Baseline scenario due to geometrical shading by buildings. In LCZ-4, even though trees that are along the streets keep an obvious cooling effect at the peak hour 14:00 (Figure 9a), the impact is not significant enough to reduce the temperature of the entire site, especially in the areas that are not covered by tree shade and building shadows. For example, the temperature at measurement point #9 increases by 0.06 • C at 14:00, which is attributed to the direct exposure to sunlight, and thereby increases the evaporation of water in the soil. Likely, the prediction in LCZ-6 reveals a very similar Ta distribution, with larger Ta reduction along the street (marked by the red color in Figure 9c), because the trees could shade the street by intercepting the solar radiation as well as cooling the air by evapotranspiration, which results in relatively low temperatures for LCZ-6.
Sustainability 2020, 12, x FOR PEER REVIEW 12 of 24 sunlight, and thereby increases the evaporation of water in the soil. Likely, the prediction in LCZ-6 reveals a very similar Ta distribution, with larger Ta reduction along the street (marked by the red color in Figure 9c), because the trees could shade the street by intercepting the solar radiation as well as cooling the air by evapotranspiration, which results in relatively low temperatures for LCZ-6.

Green roof
Comparing and contrasting among Baseline scenarios, Scenario C (50% green roof) and Scenario D (100% green roof) adding 100% green roof coverage had some cooling effect in certain regions for LCZ-6, mainly along the southeast wind direction ( Figure 10). However, the air temperature

Green Roof
Comparing and contrasting among Baseline scenarios, Scenario C (50% green roof) and Scenario D (100% green roof) adding 100% green roof coverage had some cooling effect in certain regions for LCZ-6, mainly along the southeast wind direction ( Figure 10). However, the air temperature reduction for LCZ-4 and LCZ-5 was insignificant. Additionally, the Ta reduction caused by 50% green roof was not obvious in all three LCZ classes. This result is consistent with the conclusion that as the building height increases, the microclimate impact of the green roof at pedestrian level decreases [76]. Figure 10 indicates that the elevation air temperature difference between Scenario D and the Baseline scenario in LCZ-6. The vertical average temperature drop is 0.09 • C within four meters below the green roof. As the vertical height decreases, the air temperature reduction also drops. The air temperature reduction at the pedestrian height is only 0.04 • C (Figure 10b).

Green roof
Comparing and contrasting among Baseline scenarios, Scenario C (50% green roof) and Scenario D (100% green roof) adding 100% green roof coverage had some cooling effect in certain regions for LCZ-6, mainly along the southeast wind direction ( Figure 10). However, the air temperature reduction for LCZ-4 and LCZ-5 was insignificant. Additionally, the Ta reduction caused by 50% green roof was not obvious in all three LCZ classes. This result is consistent with the conclusion that as the building height increases, the microclimate impact of the green roof at pedestrian level decreases [76]. Figure 10 indicates that the elevation air temperature difference between Scenario D and the Baseline scenario in LCZ-6. The vertical average temperature drop is 0.09 ℃ within four meters below the green roof. As the vertical height decreases, the air temperature reduction also drops. The air temperature reduction at the pedestrian height is only 0.04 ℃ (Figure 10b).

Ground albedo
In Scenario E and Scenario F, the ground albedo was increased from 0.1 to 0.3 and 0.5, respectively. As can be observed in Table 5, the maximum Ta reduction of Scenario E is 0.12-0.49 ℃ during the day. Although the cooling effect in Scenario F is more significant than that in Scenario E (0.20-0.68 ℃), similar Ta distributions are observed in LCZ-4, LCZ-5, and LCZ-6 ( Figure 11).
Compared to the Baseline scenario at 14:00 (Figure 11), the Ta reduction of Scenario E and Scenario F in three LCZs is comparatively larger in the hard-paved area, while the air temperature drops in other areas are smaller. Specifically, it can be seen that the 0. Overall, the air temperature reduction in LCZ-4 is larger than that of LCZ-5 and LCZ-6 at 8:00 and 14:00 with higher albedo ( Table 5). The air temperature reduction reaches a peak at 14:00, while it becomes negligible at 20:00 in all LCZs. The reason may due to the ground materials with higher

Ground Albedo
In Scenario E and Scenario F, the ground albedo was increased from 0.1 to 0.3 and 0.5, respectively. As can be observed in Table 5, the maximum Ta reduction of Scenario E is 0.12-0.49 • C during the day. Although the cooling effect in Scenario F is more significant than that in Scenario E (0.20-0.68 • C), similar Ta distributions are observed in LCZ-4, LCZ-5, and LCZ-6 ( Figure 11).
Compared to the Baseline scenario at 14:00 (Figure 11), the Ta reduction of Scenario E and Scenario F in three LCZs is comparatively larger in the hard-paved area, while the air temperature drops in other areas are smaller. Specifically, it can be seen that the 0.20 albedo increase case exhibits temperatures up to 0.     From the simulation results we can see that the tree quantities and albedo have a negative correlation with air temperature, which confirms the cooling effect on pedestrian air temperature modification. However, the temperature drops of green roofs remain within 0.10 ℃, even at 14:00 when the solar radiation is strong. Additionally, the cooling effects related to green roofs are not obvious regardless of whether 50% or 100% green roofs are covered, according to the linear regression in Figure 12 (c).
As can be observed in Figure 12, there is nearly a 0.14-0.40 ℃ reduction in air temperature per 30% increase in tree quantity and an approximate 0.10-0.50 ℃ decrease for every 0.20 increase in albedo at 14:00. Specifically, the slopes of the relationships are different depending on the LCZ classes. For LCZ-4, the cooling effect of increasing albedo outperforms the tree quantities and green roofs at 8:00 and 14:00. The decline in cooling capacity provided by increasing albedo 0.20 (from 0.40 ℃ reduce to less than 0.15 ℃) is larger than that of increasing tree quantity. On the contrary, the cooling effect of albedo is better than tree quantities for LCZ-5 and LCZ-6 at 8:00 and 14:00. Overall, the air temperature reduction in LCZ-4 is larger than that of LCZ-5 and LCZ-6 at 8:00 and 14:00 with higher albedo ( Table 5). The air temperature reduction reaches a peak at 14:00, while it becomes negligible at 20:00 in all LCZs. The reason may due to the ground materials with higher albedo cutting down ground surface temperature and promoting favorable evapotranspiration processes during the day.
From the simulation results we can see that the tree quantities and albedo have a negative correlation with air temperature, which confirms the cooling effect on pedestrian air temperature modification. However, the temperature drops of green roofs remain within 0.10 • C, even at 14:00 when the solar radiation is strong. Additionally, the cooling effects related to green roofs are not obvious regardless of whether 50% or 100% green roofs are covered, according to the linear regression in Figure 12c.
As can be observed in Figure 12, there is nearly a 0.14-0.40 • C reduction in air temperature per 30% increase in tree quantity and an approximate 0.10-0.50 • C decrease for every 0.20 increase in albedo at 14:00. Specifically, the slopes of the relationships are different depending on the LCZ classes. For LCZ-4, the cooling effect of increasing albedo outperforms the tree quantities and green roofs at Sustainability 2020, 12, 2752 14 of 23 8:00 and 14:00. The decline in cooling capacity provided by increasing albedo 0.20 (from 0.40 • C reduce to less than 0.15 • C) is larger than that of increasing tree quantity. On the contrary, the cooling effect of albedo is better than tree quantities for LCZ-5 and LCZ-6 at 8:00 and 14:00.

Two variables
To further investigate the cooling method combinations, six different scenarios were simulated for LCZ-4-6. The 24-hour air temperature reduction is shown for the integrated scenarios in Figure  13. Comparing among Scenario A+C, Scenario A+E, Scenario B+D, Scenario B+F, Scenario C+E, and Scenario D+F, the air temperatures of the scenarios adopted with two variables are always lower than that in the Baseline scenario. Overall, the ΔTa-max during 24 h occurs at 14:00 in the afternoon, and the ΔTa-min during 24 h occurs during 1:00 to 4:00 in the morning for all the simulated scenarios. Moreover, two-variable scenarios have a more significant average cooling effect (reduction range from 0.02 to 0.96 ℃) in LCZ-4, LCZ-5, and LCZ-6 compared to one-variable scenarios during the day (from 7:00 to 18:00). At nighttime (from 19:00 to 6:00), the temperature reductions of all scenarios are much less.
According to the previous single-variable simulation results, within the whole 24 h periods, the highest reduction of air temperature is originally obtained from large tree quantities combined with ground albedo. Even though less of an effect is found from the changes of green roof, it is easily believed that the two-variables combination has cooling potential, especially when tree quantities and higher ground albedo are integrated. Figure 12. (a) Correlation between average air temperature and tree quantities, albedo, and green roof at 8:00; (b) Correlation between average air temperature and tree quantities, albedo, and green roof at 14:00; (c) Correlation between average air temperature and tree quantities, albedo, and green roof at 20:00.

Two Variables
To further investigate the cooling method combinations, six different scenarios were simulated for LCZ-4-6. The 24-hour air temperature reduction is shown for the integrated scenarios in Figure 13. Comparing among Scenario A+C, Scenario A+E, Scenario B+D, Scenario B+F, Scenario C+E, and Scenario D+F, the air temperatures of the scenarios adopted with two variables are always lower than that in the Baseline scenario. Overall, the ∆Ta-max during 24 h occurs at 14:00 in the afternoon, and the ∆Ta-min during 24 h occurs during 1:00 to 4:00 in the morning for all the simulated scenarios. Moreover, two-variable scenarios have a more significant average cooling effect (reduction range from 0.02 to 0.96 • C) in LCZ-4, LCZ-5, and LCZ-6 compared to one-variable scenarios during the day (from 7:00 to 18:00). At nighttime (from 19:00 to 6:00), the temperature reductions of all scenarios are much less.
According to the previous single-variable simulation results, within the whole 24 h periods, the highest reduction of air temperature is originally obtained from large tree quantities combined with ground albedo. Even though less of an effect is found from the changes of green roof, it is easily believed that the two-variables combination has cooling potential, especially when tree quantities and higher ground albedo are integrated.
In fact, the combined-variable scenarios have a more significant cooling effect. For instance, Scenario A+E (increasing trees by 30%, increasing albedo by 0.2) provides higher reduction (0.96 • C), when compared to Scenario A (0.32 • C) and Scenario E (0.41 • C) in LCZ-4 at 14:00. Additionally, for the value of two-variable doubles (compare Scenario A+E and Scenario B+F), temperature reduction is also enhanced during the day. Particularly, Scenario B+F has the most significant air temperature reduction (0.30-0.96 • C) during the daytime. At the peak hour (14:00), the air temperature differences between Baseline scenario and Scenario B+F for LCZ-4, LCZ-5, and LCZ-6 are 0.96 • C, 0.83 • C, and 0.9 • C, respectively (Figure 13), which is attributed to the modification of ground materials, which could enhance the cooling effect of vegetation, with the low absorption of solar radiation and high ground surface moisture contents during the daytime. In fact, the combined-variable scenarios have a more significant cooling effect. For instance, Scenario A+E (increasing trees by 30%, increasing albedo by 0.2) provides higher reduction (0.96 ℃), when compared to Scenario A (0.32 ℃) and Scenario E (0.41 ℃) in LCZ-4 at 14:00. Additionally, for the value of two-variable doubles (compare Scenario A+E and Scenario B+F), temperature reduction is also enhanced during the day. Particularly, Scenario B+F has the most significant air temperature In addition, cooling method combinations can provide more advantageous cooling effects when applied in the appropriate LCZ. This makes a maximum average difference with magnitude up to 0.35 • C. High Ta reduction of up to 0.55 • C (in LCZ-4) can be observed when applying Scenario C+F owing to the relatively high reflection of hard paving surfaces, while LCZ-6 appears to have a decrease of 0.40 • C with Scenario C+F at 14:00. Accordingly, the combined implementation Scenario D+F is found to lead to the larger reduction of air temperature in LCZ-4 (0.85 • C), while it only causes 0.50 • C in LCZ-6. The results of Scenario A+E also vary in different LCZs, with a Ta decrease of 0.60 • C, 0.71 • C, and 0.78 • C for LCZ-4, LCZ-5, and LCZ-6, respectively. This is because the trees in the low-rise open areas (LCZ-6) are not well-shaded by buildings, reducing Ta more than the other two LCZs during the hottest period of the day.
Specifically, comparing the six simulated double-variable scenarios in LCZ-4 (open high-rise model), the rate of reductions in temperatures tend to decline during the day. Scenario A+C (an increase of tree quantities by 30% and green roof by 50%) has the lowest ∆Ta (0.13-0.43 • C), while Scenario B+F (an increase of tree quantities by 60% and albedo by 0.4) has the highest ∆Ta (0.26-0.96 • C) during the day. However, during 0:00-5:00, Scenario A+C, Scenario A+E, Scenario B+D, and Scenario B+F have slight Ta reductions compared to the Baseline scenario in LCZ-4.
Likewise, the results exhibit a similar trend in LCZ-5 and LCZ-6, where the magnitude of daytime Ta is cooler than the Baseline scenarios. Scenario A+E and Scenario B+F are found to have corresponding Ta reduction ranges of 0.25 to 0.73 • C and 0.30 to 0.82 • C in LCZ-5, respectively. Then again, for LCZ-6, the Ta reduced by 0.30-0.80 • C and 0.40-0.92 • C in Scenario A+E and Scenario B+F, which is attributed to the increased tree number and higher ground albedo.
Furthermore, it is noticeable that air temperature differences exist among all the scenarios. Comparing Scenario B+F (increasing trees by 60%, increasing albedo by 0.4) with Scenario D+F (green roof 100%, increasing albedo by 0.4) at 14:00, even though they both have albedo increased by 0.4, the cooling effect of Scenario B+F is higher than that in Scenario D+F during the day, and the difference is most significant in LCZ-6 (0.20-0.40 • C). On the one hand, it again demonstrates that tree quantities contribute more cooling than green roofs. On the other hand, it shows that the cooling effect mostly comes from adding trees in LCZ-6, because shading is found to be an effective measure to mitigate high Ta in the low building area, which lacks building shadows. However, when comparing Scenario B+D (increasing trees by 60%, green roof 100%) with Scenario D+F (green roof 100%, increasing albedo by 0.4), Scenario B+D has higher air temperature than that of Scenario D+F in LCZ-4; however, they are nearly equal in LCZ-5. In LCZ-6, Scenario B+D shows lower Ta than that of Scenario D+F, which means in LCZ-4, the heat dissipation potential is obtained through the modification of the ground material in the combination variable outweigh trees, while it seems the opposite in LCZ-6. The contribution rate of the two is equivalent for LCZ-5.
In general, all the modifications are effective in reducing the air temperature in LCZ-4, LCZ-5, and LCZ-6. This verifies the hypothesis that the urban air temperature conditions reduce as tree quantity, ground surface albedo values, and green roof area increase. The appropriate cooling method combination should also be applied when considering the maximum cooling effect. During the peak hour (14:00) the combined modification of trees and ground albedo created a more advantageous effect of air temperature reduction in LCZ-4-6. Furthermore, it can be noticed that the two-variable scenarios with the increase of trees produced a noticeable cooling effect and kept the area cooler during the peak hour in LCZ-5 and LCZ-6.

Discussion
The outputs of single-variable simulations show that tree quantities provide notable temperature reduction in LCZ-4, LCZ-5, and LCZ-6. The result agrees with some previous research results, which show that increasing tree coverage in open and bare urban areas can effectively mitigate thermal conditions significantly [77,78]. Comparing between Scenario A and Scenario B, it can be found that adding 60% more trees outweighs a 30% increase during the selected three hours. Moreover, while with the same vegetation coverage, LCZ-6 has the most significant decrease effect at 14:00. This may due to the fact that adding trees could provide more shadows to cool down the human-height air temperature of the Baseline scenario in LCZ-6, which has the highest air temperature among all three LCZs. Especially during the hottest hour, tree covers have significant impacts on cutting down solar radiation for LCZ-6.
Furthermore, the temperature reductions caused by adding trees are slightly greater than those for the high albedo modification cases. During the peak hour (14:00), the average temperature decreases 0.20-0.50 • C every 0.2 increase in ground albedo, which agrees with the previous research that the relatively low albedo cause more thermal energy storage in pavement and tends to generate a more severe UHI effect [79], while the higher albedo surfaces in these conditions are likely to reflect solar radiation quickly and store only a negligible amount [43].
As shown in Figure 12, there is a negative linear relationship between the number of trees and air temperature. Similar conclusions were also reached in some of the previous studies. For example, tree canopy coverage is negatively associated with air temperature in Phoenix and Changchun [17,80]; however, the cooling capacity of the same number of trees is quite different, due to the distinct arrangement [17,70], location [81], tree species [28,32], and climate background [17]. On the contrary, the study in Prague and Brno suggests that adding trees (a 30% increase) in the compact midrise development might even increase air temperature [82]. This may be explained by the fact that the cooling effect of urban greenery differs with the surrounding urban structure, which again reveals the importance of UHI study under the LCZ framework. Furthermore, higher ground albedo leads to a decrease of air temperature, which agrees with previous studies. When the albedo of road increased from 0.2 to 0.45, the daily-averaged urban heat island intensity of the urban area decreased by about 0.1 • C in Montreal [83]. Similarly, in Rome, a higher improvement of urban temperature is noticed during the daytime where temperature is decreased by 4 • C if the albedo of surfaces increases [84]. The model in Prague and Brno suggests that an albedo increased by 0.25 leads to a difference of −0.2 • C in daily average temperature, while an albedo increased by 0.5 leads to a daily average temperature difference of −0.5 • C in Prague and −0.4 • C in Brno [82], which also proves the negative linear relationship between the albedo and air temperature. Even though these results vary in spatial scales, the cooling capability of higher albedo cannot be ignored.
Among the three individual cooling measures, the green roofs have the least noticeable cooling effect on the thermal environment at outdoor pedestrian height, especially for open high-rise and open middle-rise residential context (LCZ-4 and LCZ-5). A previous study shows that the green roof strategy can effectively reduce the temperature of the outer surface of the roof and help reduce the building energy consumption [85]. However, the cooling impacts of the green roofs on sidewalks at pedestrian height depend on the building height [86,87]. Furthermore, no obvious linear relationship is observed between 50% and 100% of green roofs in the three LCZs. It is noticeable that the Ta drop of green roof on the pedestrian height in this study is quite insignificant compared to the very few relevant studies [87]. The reason may be attributed to the improper physical characteristics of the green roof setting in the simulation, such as leaf area index, plant height, and soil thickness [88]. Additionally, the low percentage of green roof of each site may also contribute to insignificant Ta reduction [89].
Comparing one-variable with two-variable scenarios, firstly, given the same percentage, two-variable scenarios are superior to one-variable scenarios in cooling air temperature. Secondly, sometimes single cooling methods with high value can also achieve the same effect as cooling method combinations. For instance, Scenario B (increasing trees by 60%) has higher Ta reduction (0.41-0.67 • C) than Scenario A+C (increasing trees by 30% and albedo by 0.2) (∆Ta 0.35-0.43 • C). Thirdly, even the same combination method has a different level of cooling in terms of different LCZ types. By providing an appropriate number of trees and high albedo materials, an optimum cooling level can be achieved, offering a significant reduction in the air temperature among all the simulated scenarios. In fact, a number of previous studies suggest that surface properties (i.e., building height, height-to width ratio, and sky view factor) are important factors of microclimate. As a new classification system, LCZ is a useful tool to quantify and generalize the variety of surface properties in urban areas and make comparisons between cities easily. From the discussions above, different LCZs have their favorable combinations of cooling methods. This finding agrees with some of the previous research. For instance, the study in Baghdad indicates that trees can cause a 4 K reduction of Ta in the medium density building area; however, there is nearly no influence in the high-density district [90]. Some studies also point out that trees are more effective in low urban density (high sky view factor) than in compact urban areas (low sky view factor) [91].

Limitations and Future Direction
Although the study simulates and compares 36 cooling scenarios with three current situations for LCZ-4, LCZ-5, and LCZ-6, there are still some practical limitations. Firstly, solar radiation is only generated by the assumed conditions in the model and the wind is constantly from one direction, which does not reveal the actual situation. Secondly, studies show that intensive green roofs have a better cooling effect than extensive green roofs [88,89,92]. However, the intensive green roof type is not considered as a variable in the simulation, which may lead to the result that there is a non-linear relationship between green roof and cooling effect. Further studies are needed in terms of considering leaf area index, plant height, and soil thickness as variables to explore the cooling effects of green roofs. The building's physical properties such as building layout, building surface fraction, and the number of buildings of the chosen sites may also lead to the limited area of green roof and unsatisfactory simulation results. Additionally, specific tree locations and arrangements are not taken into consideration in the research, which may have an influence on shading and ventilation [17], and the species is limited to Cinnamomum camphor, which is also one of the restricted settings. More studies are necessary to find the appropriate tree arrangements and species in terms of optimum cooling impact. Furthermore, each LCZ scenario was modelled in isolation without the surrounding environment; thus, the shading effects and ventilation between neighborhoods could not be taken into account. Lastly, the cooling methods may not be effective in another climate zone such as the tropical monsoon climate cities. Further discussions are needed in the rest of the LCZ classes and other cities of different climate zones.

Conclusions
The research aimed at exploring the cooling potential of tree quantity, ground albedo, green roofs and their combinations in LCZ-4, LCZ-5, and LCZ-6. Firstly, from the simulation results, it is demonstrated that all the cooling method scenarios have an influence on air temperature reduction, except for the extensive green roof index, which presents a negligible effect on air temperature at human height in all three LCZs regardless of its proportion. Even though 100% green roof modification affects pedestrian air temperature in low building areas (LCZ-6), the reduction of which is only 0.04 • C. In that case, the relationship between the cooling capacity and the green roof is not linear, whereas tree shade and higher albedo can provide effective local cooling.
Secondly, the contributory factor effects differ among different urban morphology types (LCZs). Increasing the quantity of trees by 60% leads to lower outdoor temperatures at pedestrian height for LCZ-5 and LCZ-6. However, due to the large but unshaded open space and inter-spacing between buildings, the effect of increasing trees by 30% or 60% is not as good as that in mid-rise and low-rise building areas, which makes some regions of LCZ-4 suffer from high air temperature. On the contrary, the effect of albedo modification on air temperature is likely to be valid in LCZ-4. Even though green roofs play the most insignificant role among three methods, they displays some Ta reduction in low buildings (LCZ-6). As they do not need to use valuable and highly competitive land at ground level, they can be used as a cooling method when the public space is limited in LCZ-6.
Thirdly, by applying combined cooling factors, a decrease in air temperature can be provoked in an open residential context (LCZ-4, LCZ-5 and LCZ-6), and especially with higher levels of tree quantity (60%) coupled with cool ground material (increased albedo of 0.4), the cooling effect is significantly enhanced. The variations in the estimated benefits among 36 scenarios reflect the importance of cooling analysis aimed at different open LCZs.
In conclusion, vegetated scenarios in LCZ-5 and LCZ-6 enhance cooling effect, especially in the 60% more trees scenarios due to shade density, while higher albedo scenarios in LCZ-4 lead to a significant reduction. Therefore, when considering one cooling strategy, this study recommends adding 60% more trees in LCZ-5 and LCZ-6, while increasing albedo (by 0.4) in LCZ-4. Additionally, Senario B+E (60% more trees and 0.4 higher ground albedo) is the most effective scenario in LCZ-4, LCZ-5, and LCZ-6 at pedestrian height. Thus, it is recommended that two major factors (trees and ground albedo) need to be considered in combination when planning and renewing open building areas, ultimately offering improvement in the thermal environment. Furthermore, it is not recommended to apply an extensive type of green roof in LCZ-4 and LCZ-5 when only considering decreasing air temperature at pedestrian height due to its negligible cooling effect. The findings obtained from this study should be extended to other LCZ types and climatic regions in future research.