Improvement Strategies for Microclimate and Thermal Comfort for Urban Squares: A Case of a Cold Climate Area in China

: Urban squares are an important part of a city’s overall spatial environment. However, many urban squares lack rational designs, causing the thermal environment to deteriorate. To ensure sustainable urban development, urban square microclimates should be improved. Given that, this study investigates the effects of three coverages of three landscape elements of urban squares through modeling and simulation using the ENVI-met model validated by ﬁeld measurements. The correlation between physiological equivalent temperature (PET) and different amounts of landscape elements is investigated using Spearman analysis. This study presents a case study of a typical urban square in a cold climate area. Design strategies in the area are proposed. The results show that the microclimate and thermal comfort of the urban square can be improved by expanding water bodies, modest increasing buildings and optimizing vegetation. Vegetation is the most important landscape element affecting thermal comfort in the urban square. The PET can be reduced by about 1.5 ◦ C by increasing the vegetation cover from 40% to 70%. However, the degree of microclimate regulation by vegetation is disturbed by water bodies and buildings (| ρ | ≥ 0.5). Therefore, to achieve a more comfortable thermal environment, a combination of landscape elements should be considered.


Introduction
Large-scale climate projections show a strong trend of increasing temperatures over the next few decades, including all-day temperature increases, longer warm periods and more intense heatwaves [1].As temperatures rise in the future, extreme weather events will become increasingly common [2].Urban populations will be at greater risk from the impact of the urban heat island (UHI) effect, including low-quality outdoor environments [3], atmospheric pollution [4], risks to human health [5] and increased energy consumption in buildings [6].At present, much progress has been made in designing strategies to reduce the impact of indoor comfort on occupants [7].However, due to the complexity of outdoor environments compared to indoor environments, the thermal comfort of open spaces has not been adequately studied [8].As part of cities, the mitigation of UHI effects in urban open spaces is necessary to improve the outdoor thermal comfort of residents.Outdoor thermal comfort affects a person's physical stress and emotional calmness, which are important for the user's choice of space.The most important climatic factor affecting outdoor thermal comfort is the temperature and solar radiation in the microclimate [9][10][11][12][13].Therefore, improving the urban microclimate to create a healthy outdoor thermal environment is an important topic for urban planning.Some human thermal comfort indicators based on microclimate parameters have been widely used for thermal comfort assessment in different outdoor climate regions, such as Standard Effective Temperature (SET*) [14] and physiological equivalent temperature (PET) [15].Urban squares are an important part of urban open space and the main venue for outdoor activities.In recent years, with the development of urbanization, many large squares have been built in China and other countries, but these spaces need improvement [16].Since people are exposed to the outdoors for long periods when using urban squares, improving the thermal comfort of urban squares has become a top priority for urban designers.
Previous research has shown that urban landscape configurations can effectively alter outdoor thermal environments at the microscale to mitigate the UHI effect [17].It is evident from various research that vegetation and shading facilities can be used as an effective method to improve thermal comfort in urban squares.Research in Germany showed that maximizing the surface area of tree shade could provide a 5.2% cooling effect on a local urban square at 15:00 [18].Another example of vegetation is in the United States, in California, where increasing vegetation was the most effective daytime strategy for extreme and average days.The maximum cooling temperature of vegetation reached −3.5 • C on extreme local days [19].In Argentina, the concentration of 60 percent of the forest land around closed centers is the most effective program for the restoration of memorial squares [20].There is also a growing amount of literature that suggests that the height of buildings plays an important role in the quality of the thermal environment.Changing building height has been reported to result in significant changes in surrounding thermal comfort [21,22].The tallest buildings can provide more comfortable thermal conditions when they are in the center of the layout [23].Similarly, the thermal comfort of an area can be improved by increasing building coverage [24].However, over-increasing building coverage can lead to the deterioration of thermal comfort [22].Finally, although few studies have discussed the design of water bodies in urban plazas, in general, water bodies are an essential element in urban areas.A study in Iran showed that evaporative cooling of water is an effective strategy to improve microclimatic conditions in urban squares in temperate and arid regions [25].From a case in Italy, it was seen that expanding water surfaces could reduce the temperature of a city square in the center of Rome by 0.9 • C. Locally, the expansion of water surfaces improved thermal comfort, which was second only to the improvement of vegetation [12].
In general, previous studies have provided a lot of new insight into the variation of thermal comfort in urban squares, such as exploring the influence of the combination type of trees on microclimate [26][27][28], exploring the influence of the shape and size of water bodies on urban climate [29], the influence of building orientation on the thermal comfort of urban open spaces [30], etc.However, these studies have mainly focused on the improvements brought by optimal planning layouts, and there is a lack of quantitative studies on landscape elements.Moreover, papers on this topic have generally focused on a single landscape element and still lack a systematic understanding of how multiple factors contribute to thermal comfort.This paper uses multi-factor case studies (including buildings, water bodies and vegetation) and quantitative research methods to make a comparative study of various scenarios, thus filling the gap in comprehensive research in this field.The simulation software ENVI-met was also used.Experimental data such as air temperature and relative humidity were collected by actual measurements.Boundary conditions were determined, and numerical models were calibrated.
Numerous studies have confirmed that the same measure has different effects in different zones [31,32].Therefore, it is necessary to consider the local climate, the actual situation and the user's perception in practice.Studies on the outdoor thermal environment of squares in China have been conducted in severely cold climate areas [33], hot summer and cold winter climate areas [30] and hot summer and warm winter climate areas [34].However, few studies have been conducted in cold climate areas due to the lack of educational resources in northern China [35].Therefore, it is necessary to supplement the thermal comfort studies of urban squares in this area.
This study is based on climatic conditions in cold climates and has three main research objectives.The first objective is to assess and rank the impact of three configurations of three landscape elements on microclimate parameters in urban squares; the second objective is to assess the correlation between the impact of building coverage, water coverage and vegetation coverage on thermal comfort and to explore its internal mechanisms; and the third objective is to develop and test the applicability of various thermal environment improvement design strategies to provide a reference for urban designers.

Research Area
Xi'an is located between 107 • 40 -109 • 49 E and 33 • 39 -34 • 45 N. It is geographically located in the center of China.Among the five climate areas in China, Xi'an is a typical city in a cold climate area.A cold climate area is defined as an area with an average temperature of the coldest month being −10-0 • C and the average daily temperature being ≤5 • C for 90-145 days [36].For a long time, urban planning and construction in cold climate areas have focused on keeping warm in winter, and little attention has been paid to the heat in summer.The historical meteorological data of Xi'an are shown in Table 1.Xi'an has a hot summer with a significant UHI effect, which is not conducive to the daily life and outdoor activities of residents.The Big Wild Goose Pagoda Square was built in 589 A.D. After more than 1000 years of rising and falling and expansion, it is an iconic square in Xi'an and was listed among the World Heritage Sites in 2013 [37].This study was conducted in the Big Wild Goose Pagoda's North Square (Figure 1).The Square is made up of a water fountain, a cultural square, a garden landscape, a cultural promenade and tourism commercial facilities, with the tower serving as the center axis.It measures 480 m east to west and 350 m north to south and covers 168,000 m 2 .
Buildings 2022, 12, 944 3 of 20 second objective is to assess the correlation between the impact of building coverage, water coverage and vegetation coverage on thermal comfort and to explore its internal mechanisms; and the third objective is to develop and test the applicability of various thermal environment improvement design strategies to provide a reference for urban designers.

Research Area
Xi'an is located between 107°40′-109°49′ E and 33°39′-34°45′ N. It is geographically located in the center of China.Among the five climate areas in China, Xi'an is a typical city in a cold climate area.A cold climate area is defined as an area with an average temperature of the coldest month being −10-0 °C and the average daily temperature being ≤5 °C for 90-145 days [36].For a long time, urban planning and construction in cold climate areas have focused on keeping warm in winter, and little attention has been paid to the heat in summer.The historical meteorological data of Xi'an are shown in Table 1.Xi'an has a hot summer with a significant UHI effect, which is not conducive to the daily life and outdoor activities of residents.The Big Wild Goose Pagoda Square was built in 589 A.D. After more than 1000 years of rising and falling and expansion, it is an iconic square in Xi'an and was listed among the World Heritage Sites in 2013 [37].This study was conducted in the Big Wild Goose Pagoda's North Square (Figure 1).The Square is made up of a water fountain, a cultural square, a garden landscape, a cultural promenade and tourism commercial facilities, with the tower serving as the center axis.It measures 480 m east to west and 350 m north to south and covers 168,000 m 2 .Xi'an has four distinct seasons, with the hottest month of July and the coldest month of January throughout the year.In this study, field measurements were performed in the Big Wild Goose Pagoda's North Square from 27-29 July 2020.Four measurement sites were set up at the north-south road, entrance square, east-west road and green space for this study to investigate the various influencing elements (Figure 1).Site 2 was used as a reference site because it is the closest to the boundary without the influence of shade all day.All instruments were ISO 7726 compliant and were mounted on tripods in the square at a height of 1.4 m above ground.Detailed information on the instruments is in Table 2.In addition to the measured climatic parameters, data from local weather stations were consulted.These observations were utilized to determine the microclimate conditions in the study area and served as the foundation for the square simulation model.Figure 2 shows the framework for this study.Xi'an has four distinct seasons, with the hottest month of July and the coldest month of January throughout the year.In this study, field measurements were performed in the Big Wild Goose Pagoda's North Square from 27-29 July 2020.Four measurement sites were set up at the north-south road, entrance square, east-west road and green space for this study to investigate the various influencing elements (Figure 1).Site 2 was used as a reference site because it is the closest to the boundary without the influence of shade all day.All instruments were ISO 7726 compliant and were mounted on tripods in the square at a height of 1.4 m above ground.Detailed information on the instruments is in Table 2.In addition to the measured climatic parameters, data from local weather stations were consulted.These observations were utilized to determine the microclimate conditions in the study area and served as the foundation for the square simulation model.Figure 2 shows the framework for this study.

Model Setup
ENVI-met V4.4 (ENVI-met) was used for simulations in this study.A complete ENVI-met simulation model has four subsystems (atmosphere, soil, vegetation and buildings), and the model database provides a large number of environmental parameter options [38].The software includes a tool called "Biomet" which allows the math of various thermal comfort indexes.The ENVI-met model is one of the most widely used dynamic simulation tools for microclimate analysis.Researchers use it not only to investigate current microclimate conditions but also to compare and evaluate the effects of various strategies on UHI effects [39].

Model Setup
ENVI-met V4.4 (ENVI-met) was used for simulations in this study.A complete ENVImet simulation model has four subsystems (atmosphere, soil, vegetation and buildings), and the model database provides a large number of environmental parameter options [38].The software includes a tool called "Biomet" which allows the math of various thermal comfort indexes.The ENVI-met model is one of the most widely used dynamic simulation tools for microclimate analysis.Researchers use it not only to investigate current microclimate conditions but also to compare and evaluate the effects of various strategies on UHI effects [39].
The ENVI-met model of the study area was constructed based on field surveys and aerial photographs.The horizontal modeling grid cells were 4 × 4 m, and the height was 2 m.This model was used to create a simulation model with 83 × 99 horizontal cells and Buildings 2022, 12, 944 5 of 21 30 vertical cells.To avoid boundary effects, 5 layered grids were placed on either side of the simulation model.The materials of each element in the square were chosen based on their current state.It should be noticed that the water body in the study area includes timed fountains, which, in this study, were simplified to shallow pools because the fountains only run for 20 min during the measuring period.An actual measured uniform water depth of 0.5 m in the ENVI-met model was used because water levels greater than 0.45 m did not affect the simulation results [40].
Previous studies have emphasized that changing boundary settings, material properties and experimental protocols can improve the correlation between simulations and measurements [41,42].In this study, the weather station was located on the outskirts of Xi'an city, far from the measurement site in the city center.To reduce the error between the simulation and measurements, the input meteorological parameters of the model needed to be modified, and the modified data were used as the original boundary condition of modeling.Upon completion of the data collection, we compared the measurements from the reference station with the data from the weather station.As shown in Figure 3a,b, the measured temperature per hour at the weather station was 1-4.5 • C below the measured data, and the RH was 3-30% percent higher than the measured data.Therefore, in this study, modified weather station data that increased the AT of the weather station by 7% and decreased the RH by 20% were used as the input data for the model.Table 3 shows the microclimate simulation settings and meteorological parameters that were entered by the model in ENVI-met.
The ENVI-met model of the study area was constructed based on field surveys and aerial photographs.The horizontal modeling grid cells were 4 × 4 m, and the height was 2 m.This model was used to create a simulation model with 83 × 99 horizontal cells and 30 vertical cells.To avoid boundary effects, 5 layered grids were placed on either side of the simulation model.The materials of each element in the square were chosen based on their current state.It should be noticed that the water body in the study area includes timed fountains, which, in this study, were simplified to shallow pools because the fountains only run for 20 min during the measuring period.An actual measured uniform water depth of 0.5 m in the ENVI-met model was used because water levels greater than 0.45 m did not affect the simulation results [40].
Previous studies have emphasized that changing boundary settings, material properties and experimental protocols can improve the correlation between simulations and measurements [41,42].In this study, the weather station was located on the outskirts of Xi'an city, far from the measurement site in the city center.To reduce the error between the simulation and measurements, the input meteorological parameters of the model needed to be modified, and the modified data were used as the original boundary condition of modeling.Upon completion of the data collection, we compared the measurements from the reference station with the data from the weather station.As shown in Figure 3a,b, the measured temperature per hour at the weather station was 1-4.5 °C below the measured data, and the RH was 3-30% percent higher than the measured data.Therefore, in this study, modified weather station data that increased the AT of the weather station by 7% and decreased the RH by 20% were used as the input data for the model.Table 3 shows the microclimate simulation settings and meteorological parameters that were entered by the model in ENVI-met.

Redesign Model
The current coverage of water bodies, vegetation and buildings in the study area is 10%, 20% and 15%, respectively.Previous studies have concluded that people prefer to be exposed to green spaces with water [43].According to the characteristics of the cultural theme square, more than 80% of these spaces are designed with water bodies.The average proportion of water area ranges from 5% to 25%, with a minimum of about 5% and a maximum of about 30% [44].On this basis, the water coverage in the model was redesigned to 5% and 30% to compare with the current coverage of the study area (10%).These cases are symbolized by W5%, W10% and W30%, as shown in Figure 4a.Greenspace has been identified as an important factor influencing the landscape preference of the square [45].According to the Chinese urban green space classification standard, the percentage of the green area is classified as less than 35%, 35% to 65% and greater than 65% [46].Based on this criterion, the model was redesigned to 40% and 70% for comparison with the status quo (20%).These cases are symbolized by T20%, T40% and T70%, as shown in Figure 4b.Previous studies have proven that building thermal comfort for the whole day outdoors is mainly affected by building height, and building distance and layout form have negligible effects on the average outdoor thermal comfort of spaces [47,48].Buildings in the study area are currently predominantly two-story, with some single-story structures and a building coverage of 15%.Considering the influence of building coverage and height on the thermal environment of the site, the building coverage was redesigned to 5% and 30% to compare with the current situation (15%).These cases are symbolized by B5%, B15% and B30%, as illustrated in Figure 4c.

Redesign Model
The current coverage of water bodies, vegetation and buildings in the study area is 10%, 20% and 15%, respectively.Previous studies have concluded that people prefer to be exposed to green spaces with water [43].According to the characteristics of the cultural theme square, more than 80% of these spaces are designed with water bodies.The average proportion of water area ranges from 5% to 25%, with a minimum of about 5% and a maximum of about 30% [44].On this basis, the water coverage in the model was redesigned to 5% and 30% to compare with the current coverage of the study area (10%).These cases are symbolized by W5%, W10% and W30%, as shown in Figure 4a.Greenspace has been identified as an important factor influencing the landscape preference of the square [45].According to the Chinese urban green space classification standard, the percentage of the green area is classified as less than 35%, 35% to 65% and greater than 65% [46].Based on this criterion, the model was redesigned to 40% and 70% for comparison with the status quo (20%).These cases are symbolized by T20%, T40% and T70%, as shown in Figure 4b.Previous studies have proven that building thermal comfort for the whole day outdoors is mainly affected by building height, and building distance and layout form have negligible effects on the average outdoor thermal comfort of spaces [47,48].Buildings in the study area are currently predominantly two-story, with some single-story structures and a building coverage of 15%.Considering the influence of building coverage and height on the thermal environment of the site, the building coverage was redesigned to 5% and 30% to compare with the current situation (15%).These cases are symbolized by B5%, B15% and B30%, as illustrated in Figure 4c.

Thermal Comfort Index
In the 55th ASHRAE Standard, human thermal comfort is clearly defined as a psychological state of satisfaction with the thermal environment [49].The relationship between physiological factors and outdoor thermal comfort begins with heat exchange between the body and its surroundings.In normal times, people can dissipate heat by evaporation, convection and heat transfer to maintain heat balance with their surroundings [50].This thermoregulatory activity can change the physiological parameters of the body accordingly.Therefore, given the correlation between physiological parameters and thermal comfort, although human thermal comfort is a subjective feeling, the relevant physiological parameters can be used as an objective evaluation indicator of human thermal comfort [51].Dozens of thermal comfort evaluation criteria can be used for reference by researchers, including many thermal comfort evaluation criteria based on the steady-state heat transfer model and dynamic heat transfer model [52,53].Over the last two decades, thermal comfort research has moved toward more complex dynamic conditions [54].Thermal comfort prediction methods can be divided into three categories: (1) thermal equilibrium methods (mainly laboratory studies); (2) adaptive regression models (mainly field studies); and (3) adaptive heat balance [52].
Several common metrics such as Standard Effective Temperature (SET*) [14], Predicted Average Vote (PMV), General Thermal Climate Index (UTCI) [55] and Physiological Equivalent Temperature (PET) [15] have been developed to evaluate thermal comfort in different types of outdoor environments [56,57].Based on the existing literature review, PET and UTCI are the main evaluation indices that can fit cold climate areas [58].Physiologically equivalent temperature (PET) has been used more frequently in outdoor thermal comfort studies in recent years because it is the only one that can be applied to a variety of climate areas with temperatures ranging from −50 to 50 • C [15].Many studies have demonstrated that the range of thermal comfort varies with area [40,[59][60][61].Since both Tianjin and Xi'an belong to a cold climate area, the PET adapted to the local climate of Tianjin, as modified by D Lai et al., which was used as a criterion for evaluating the thermal environment in this study (Table 4) [62].

Simulation Data Validation
In this study, a statistical analysis of the simulation and measurement data was performed using the Statistical Package for the Social Sciences (SPSS).As previous studies have proposed, AT and relative humidity RH have been used to test ENVI-met models as important microclimate parameters [12,63].Therefore, the ENVI-met model was validated using AT and RH measurements from the four sites.The results are considered more accurate when the root-mean-square deviation (RMSE) values approach zero and when R 2 values approach one.Figure 5 shows a comparison between the simulated model and the field data.It was calculated that the RMSE in AT was 0.92 • C and that the RMSE in RH was 1.061.In the ENVI-met model of simple force, the measurement deviated from the simulation, which may be attributed to the static cloud and wind condition assumptions.In addition, inadequate measurements of air flow [64], wall and roof material properties and the use of default modeling were also contributing factors.Nevertheless, the RSME values that were obtained in the current study are in the range of those summarized by Salata et al. in a recent review of ENVI-met simulations (RMSE in AT is less than 4.30 • C, and RMSE in RH is less than 10.2) [65].In addition, as shown in Figure 5a,b, the R 2 in AT and RH were 0.79 and 0.86, indicating a high correlation between the measured and simulated values.The above data analysis shows that ENVI-met can be applied to the microclimate simulation of Xi'an urban squares.
values approach one.Figure 5 shows a comparison between the simulated model and the field data.It was calculated that the RMSE in AT was 0.92 °C and that the RMSE in RH was 1.061.In the ENVI-met model of simple force, the measurement deviated from the simulation, which may be attributed to the static cloud and wind condition assumptions.In addition, inadequate measurements of air flow [64], wall and roof material properties and the use of default modeling were also contributing factors.Nevertheless, the RSME values that were obtained in the current study are in the range of those summarized by Salata et al. in a recent review of ENVI-met simulations (RMSE in AT is less than 4.30 °C, and RMSE in RH is less than 10.2) [65].In addition, as shown in Figure 5a,b, the R 2 in AT and RH were 0.79 and 0.86, indicating a high correlation between the measured and simulated values.The above data analysis shows that ENVI-met can be applied to the microclimate simulation of Xi'an urban squares.

PET Distribution in the Status Quo
The PET values were calculated for the research area by Biomet.As previous studies have proposed, the worst thermal comfort conditions always occur between 14:00 and 16:00 [30,66,67], both in real measurements and in ENVI-met modeling.Figure 6 shows the PET values for the north square of the Big Wild Goose Pagoda from 9:00 to 18:00.Since the distribution of PETs within each hour does not satisfy a normal distribution, the mean of this data set is more representative than the median.The line of the means shows that the mean PET values in the study area were the highest at 14:00, reaching "very hot" (PET > 46 • C) in summer.Therefore, this section focuses on the microclimate and thermal comfort of the square at 12:00.
the distribution of PETs within each hour does not satisfy a normal distribution, th of this data set is more representative than the median.The line of the means show the mean PET values in the study area were the highest at 14:00, reaching "very hot > 46 °C) in summer.Therefore, this section focuses on the microclimate and therma fort of the square at 12:00.

Microclimate Conditions of Redesigned Cases
The ENVI-met model was tested in Section 3.1, and it was shown to provide co ing simulation results.Given this, simulations can be used as an alternative to r reconstructions.ENVI-met was used to simulate and test the effects of microclima thermal comfort modifications in the redesigned study area.According to the rede scheme (Figure 4), the correlated ENVI-met model was built and simulated to obt spatial and temporal distributions of microclimate data, including AT, RH and W PET values were calculated for the four measurement sites by Biomet.

The Impact of Differing Water Coverage on the Microclimate Parameters
The distribution maps (Figure 7) of microclimate parameters (AT, RH and WS that increasing the coverage of the water body effectively changes the AT above an the water, and the temperature distribution within the square decreases with d from the water body for all three cases.This is because the specific heat capacity o is greater than that of the ground.Comparing the three cases, we can see that adding bodies at different locations in the square also has different effects, with the dow areas of the square being more influenced by the water bodies than the upwind a is also clear that there are spatial differences in the effects of increasing the water a AT.For example, the cooling effect was greater in the not-under-tree space than under-tree space.The difference in RH due to increasing water coverage was mor ous than the difference in AT.The difference in RH between the water body a ground stems from their difference in evaporation capacity.
With respect to wind speed, as shown in Figure 7, the difference in wind spe tween the three cases was less than 0.06 m/s.The increase in water area changed t face roughness, resulting in a slight change in wind speed.The difference in tempe between land and water resulted in different wind speed trends.In general, an in in water body coverage led to a decrease in AT, an increase in RH and a slight cha

Microclimate Conditions of Redesigned Cases
The ENVI-met model was tested in Section 3.1, and it was shown to provide convincing simulation results.Given this, simulations can be used as an alternative to realistic reconstructions.ENVI-met was used to simulate and test the effects of microclimate and thermal comfort modifications in the redesigned study area.According to the redesigned scheme (Figure 4), the correlated ENVI-met model was built and simulated to obtain the spatial and temporal distributions of microclimate data, including AT, RH and WS.The PET values were calculated for the four measurement sites by Biomet.

The Impact of Differing Water Coverage on the Microclimate Parameters
The distribution maps (Figure 7) of microclimate parameters (AT, RH and WS) show that increasing the coverage of the water body effectively changes the AT above and near the water, and the temperature distribution within the square decreases with distance from the water body for all three cases.This is because the specific heat capacity of water is greater than that of the ground.Comparing the three cases, we can see that adding water bodies at different locations in the square also has different effects, with the downwind areas of the square being more influenced by the water bodies than the upwind areas.It is also clear that there are spatial differences in the effects of increasing the water area on AT.For example, the cooling effect was greater in the not-under-tree space than in the under-tree space.The difference in RH due to increasing water coverage was more obvious than the difference in AT.The difference in RH between the water body and the ground stems from their difference in evaporation capacity.With respect to wind speed, as shown in Figure 7, the difference in wind speed between the three cases was less than 0.06 m/s.The increase in water area changed the surface roughness, resulting in a slight change in wind speed.The difference in temperature between land and water resulted in different wind speed trends.In general, an increase in water body coverage led to a decrease in AT, an increase in RH and a slight change in WS.Therefore, an appropriate increase in water area can provide an effective cooling effect in the design of urban squares.

The Impact of Differing Vegetation Coverage on the Microclimate Parameters
Previous studies have demonstrated that the increased canopy density of vegetation attenuates the effect of solar radiation and decreases AT [68].In this study, the distribution of microclimate was investigated for vegetation coverage of 10%, 40% and 70%.The results show that the relationship between vegetation coverage and AT in the square followed this pattern (Figure 8).Moreover, the effect of vegetation coverage on AT had spatial differences.For instance, because buildings block part of the solar radiation, the cooling effect of vegetation is reduced under building shade [69,70].In addition, vegetation has different cooling effects around buildings with different orientations.As shown in Figure 8, the cooling effect of vegetation around the east-west oriented buildings was better than that of the north-south oriented buildings.With respect to relative humidity, vegetation converts liquid water into water vapor by transpiration to reduce the temperature and thermal comfort of leaves and surrounding air.The effect of RH by changes in vegetation coverage was pronounced in the summer.With respect to wind speed, the temperature difference caused airflow, resulting in a weak change in wind speed in the area where trees increased or decreased.WS was reduced in increased arborization because the dense tree canopy blocked the wind flow.Since this variation difference is less than 0.15 m/s, it can be considered negligible.In general, because the wind environment is quiet all year in cold climate areas, changing the vegetation coverage of the square has no significant effect on the WS in summer.

The Impact of Differing Building Heights and Densities on the Microclimate Parameters
Xi'an has a long history as the capital of thirteen ancient dynasties.The planning of the modern city of Xi'an still retains the style of the ancient capital city, following the layout of horizontal and vertical streets and symmetrical buildings.The vertical layout of the building complex is mostly north-south oriented, with the main building in the middle.The ancillary buildings are oriented east-west and are located on either side.The overall height of the complex is high in the north and low in the south.The layout of the study area also follows this style.It has been demonstrated that different urban spatial patterns lead to multiple reflections of solar radiation between buildings, which affects the flow of wind in the city and ultimately the dissipation of heat within the city [42].In case B30%, the accumulated solar radiation was absorbed, and the heat that was released into the surrounding space was greater than it was in the other two cases as the buildings increased overall, resulting in a higher AT in this case than in the other two cases.For the overall square, the buildings caused changes in the thermal environment mainly by affecting solar radiation.In this study, although buildings have a local impact on the thermal environment through their processes of absorbing, releasing and blocking heat, the AT of the whole square was not significantly affected by changing the building height and density (Figure 9).
With respect to relative humidity, a comparison of the three cases reveals that the changes in building height and density caused only slight changes in the RH distribution within the square.With respect to wind speed, past studies have shown that increasing building height and density decreases the overall wind speed.In this study, B30% was the case with the largest building coverage and height.In the summer, the wind speed distribution of B30% was lower than it was in the other two cases.The increase in buildings increased the wind shadow area in the vicinity, which led to a decrease in WS in the wind shadow area.

Thermal Comfort of Redesigned Cases
Since sites 2 and 4 are more representative, the PET values were calculated and compared between 9:00 and 18:00 for sites 2 and 4 after the square redesign (Figure 10). Figure 6 shows that the thermal comfort of the study area in summer is classified as "very hot" after 12:00.The PET values in the redesigned square did not decrease significantly, but the human comfort was reduced from "very hot" to "hot".This indicates that the redesign of the square had a moderating effect on thermal comfort.For instance, as shown in Figure 10, rising vegetation coverage usually improves thermal comfort in summer.

Internal Mechanics
The simulation results show that the thermal comfort regulation effect of a single landscape element can be influenced by other landscape elements (Figures 7-9), For example, the cooling effect of trees was lower near buildings than it was in open spaces.Therefore, increasing the understanding of the correlation between landscape elements can help to improve their thermal comfort enhancement efficiency.A correlation analysis is an analysis of the degree of correlation between two variables.There are three ways to calculate the correlation analysis, which are Pearson's correlation coefficient, Spearman's correlation coefficient and Kendall's tau-b correlation coefficient.In this study, the correlation between the effects of vegetation, water bodies and buildings on PET was analyzed.It was necessary to test whether each variable conformed to a normal distribution to select a suitable analysis method.Since the median deviation from the mean was large and most of the variables did not conform to a normal distribution, Spearman correlation analysis (Spearman's ρ) was used.In statistics, Spearman's rank correlation coefficient is named for Charles Spearman.The rank score can be introduced in a correlation analysis when the data do not satisfy the normal distribution, or when some quantities cannot be expressed in the data [71].More specifically, the two variables are first sorted to obtain the corresponding rank numbers, and the rank numbers are used to replace the original data and are then brought into the Pearson correlation coefficient formula to obtain the Spearman correlation coefficient:

Internal Mechanics
The simulation results show that the thermal comfort regulation effect of a single landscape element can be influenced by other landscape elements (Figures 7-9), For example, the cooling effect of trees was lower near buildings than it was in open spaces.Therefore, increasing the understanding of the correlation between landscape elements can help to improve their thermal comfort enhancement efficiency.A correlation analysis is an analysis of the degree of correlation between two variables.There are three ways to calculate the correlation analysis, which are Pearson's correlation coefficient, Spearman's correlation coefficient and Kendall's tau-b correlation coefficient.In this study, the correlation between the effects of vegetation, water bodies and buildings on PET was analyzed.It was necessary to test whether each variable conformed to a normal distribution to select a suitable analysis method.Since the median deviation from the mean was large and most of the variables did not conform to a normal distribution, Spearman correlation analysis (Spearman's ρ) was used.In statistics, Spearman's rank correlation coefficient is named for Charles Spearman.The rank score can be introduced in a correlation analysis when the data do not satisfy the normal distribution, or when some quantities cannot be expressed in the data [71].More specifically, the two variables are first sorted to obtain the corresponding rank numbers, and the rank numbers are used to replace the original data and are then brought into the Pearson correlation coefficient formula to obtain the Spearman correlation coefficient: Data were assigned a corresponding rank based on their average descending position in the overall data.The PET values of site 2 from 10:00 to 18:00 in the summer are taken as an example (Table 5): By subtracting the corresponding elements of the two observed variables to obtain a difference, the above equation can also be transformed into: where is the difference between the two ranks of each observation, and n is the number of observations.Spearman's ρ is a parentless measure of two-variable dependence that uses a monotonic equation to evaluate the correlation between them.If the data do not contain any repeated values and the two variables are exactly monotonically connected, the Spearman's ρ is either +1 or −1.Spearman's ρ uses the rank order magnitude of the two variables for analysis and is a nonparametric statistical method [72].It applies to continuous variables that do not meet the requirements of a normal distribution of Pearson correlation coefficients.It can also be used to measure the correlation between ordered categorical variables.
In this study, nine sets of cases were analyzed to clarify the internal mechanics between the three elements.The results are shown in Figure 11.There was a correlation between PET values between changes in the quantity of water bodies and buildings.The PET values were between the quantity of vegetation and water bodies, and quantity of vegetation and buildings were positively correlated.The PET values between the quantity of vegetation and water bodies had a lesser correlation, and all other elements' relationships were positive.Landscape architects are often confused by different possible choices in the design process.This finding can provide a reference for which landscape factor should be adjusted as a priority in the process of optimizing thermal comfort in completed projects.We hope that our study can provide design guidance for urban plaza landscape configurations under similar climatic conditions. of vegetation and water bodies had a lesser correlation, and all other elements' relationships were positive.Landscape architects are often confused by different possible choices in the design process.This finding can provide a reference for which landscape factor should be adjusted as a priority in the process of optimizing thermal comfort in completed projects.We hope that our study can provide design guidance for urban plaza landscape configurations under similar climatic conditions.

Discussion
A theoretical implication that can be drawn from the third section is that the design of squares in cold climate areas should consider the coverage of landscape elements such as water bodies, vegetation and buildings.For water bodies, the rise in evaporation and ventilation due to the increase in water bodies can improve thermal comfort.It is possible to reduce the thermal comfort of a square from near "very hot" to "hot" or "warm".For vegetation, increasing the vegetation coverage usually increases the thermal comfort of the square.Therefore, the vegetation proportion design needs to be optimized to achieve the desired thermal perception.For buildings, the shading provided by buildings increases thermal comfort to some extent in the summer.For example, as shown in Figure 12, the PET value of case B30% is about 0.82 • C lower than that of cases B5% and B15%.The importance of landscape elements impacting summer thermal comfort is ranked as follows: vegetation, water bodies and buildings.
This study also explores the role of interactions among the three landscape elements.Specifically, there is a correlation between the effects of water bodies, vegetation and buildings on PET, and thus, a linear model should not be used when building the mathematical model.Water bodies and buildings do not influence each other to a high degree on the site's PET, so designers do not have to consider the effect of the combination of these two elements when designing.In contrast, the influence brought by vegetation can be influenced by both water bodies and buildings, which prompts designers to give special consideration to the location of vegetation when designing a square.In the summer, there is a slight interaction between vegetation and water bodies, which is influenced by buildings at a higher level.At the same time, there is also an interaction between water bodies and buildings, which suggests that designers can achieve better results when choosing the right combination according to the interactions of several factors in the design.
the square.Therefore, the vegetation proportion design needs to be optimized to achieve the desired thermal perception.For buildings, the shading provided by buildings increases thermal comfort to some extent in the summer.For example, as shown in Figure 12, the PET value of case B30% is about 0.82 °C lower than that of cases B5% and B15%.The importance of elements impacting summer thermal comfort is ranked as follows: vegetation, water bodies and buildings.This study also explores the role of interactions among the three landscape elements.Specifically, there is a correlation between the effects of water bodies, vegetation and buildings on PET, and thus, a linear model should not be used when building the mathematical model.Water bodies and buildings do not influence each other to a high degree on the site's PET, so designers do not have to consider the effect of the combination of these two elements when designing.In contrast, the influence brought by vegetation can

Conclusions
This study is a pioneering study that has successfully quantified the impact of multifactor landscape elements (vegetation, buildings and water bodies) on thermal comfort in summer urban squares.Our results successfully quantify and compare the effect of landscape element change ratios on thermal comfort improvement in urban squares in cold climate areas.The main findings of the study are as follows: (1) Modified boundary conditions allow ENVI-met simulations to achieve convincing simulation accuracy in urban squares in cold climate areas (R 2 > 0.78, RMSE < 2).(2) Increasing vegetation coverage is the most effective strategy to improve microclimate and thermal comfort, followed by increasing water coverage and modest increasing building height and density.The ranking of important landscape elements that affect thermal comfort in summer is vegetation, water bodies and buildings.(3) There is a correlation between different types of landscape elements for thermal comfort enhancement (|ρ| ≥ 0.5), so they should be considered together.This method of arrangement can be used by designers in the pre-planning and project renovation stages.
It should be mentioned that this study focuses on only three main landscape elements, whereas in reality, the types of elements and their interactions with microclimate conditions are much more complex, including elements such as the airflow around the perimeter of the building enclosing the square.Despite these limitations, this study has certainly increased

Figure 2 .
Figure 2. The framework of this study.

Figure 3 .
Figure 3.The comparison of the reference point data and weather station data: (a) AT; (b) Rh.

Figure 3 .
Figure 3.The comparison of the reference point data and weather station data: (a) AT; (b) Rh.

Figure 4 .
Figure 4. Redesigned scenarios of the (a) water bodies, (b) trees and (c) buildings of the square.

Figure 5 .
Figure 5. Linearity between the simulated data and the measured data: (a) AT; (b) RH.

Figure 5 .
Figure 5. Linearity between the simulated data and the measured data: (a) AT; (b) RH.

Figure 6 .
Figure 6.PET variations at different times.

Figure 6 .
Figure 6.PET variations at different times.

Figure 7 .
Figure 7. Distribution maps of microclimate parameters at 14:00 with redesigned water body coverage.

Figure 8 .
Figure 8. Distribution maps of microclimate parameters at 14:00 with redesigned vegetation coverage.

Figure 9 .
Figure 9. Distribution maps of microclimate parameters at 14:00 with redesigned building heights and densities.

Figure 9 .Figure 10 .
Figure 9. Distribution maps of microclimate parameters at 14:00 with redesigned building heights and densities.

Figure 10 .
Figure 10.Time course of PET values at sites 2 and 4 in the redesigned cases: (a) water bodies; (b) vegetation; and (c) buildings.

Figure 11 .
Figure 11.The Spearman's rank correlation coefficient of landscape elements based on PET.Figure 11.The Spearman's rank correlation coefficient of landscape elements based on PET.

Figure 11 .
Figure 11.The Spearman's rank correlation coefficient of landscape elements based on PET.Figure 11.The Spearman's rank correlation coefficient of landscape elements based on PET.

Table 1 .
Historical weather in Xi'an (from May 2020 to May 2021).

Table 1 .
Historical weather in Xi 'an (from May 2020 to May 2021).

Table 2 .
The specifications of the instruments that were used.

Table 2 .
The specifications of the instruments that were used.
Figure 2. The framework of this study.

Table 3 .
Microclimate simulation settings and meteorological parameters of the model input in ENVI-met.

Table 3 .
Microclimate simulation settings and meteorological parameters of the model input in ENVI-met.

Table 4 .
Relationship between thermal sensation and PET for Tianjin residents in cold climate areas.