Field Assessment of Neighboring Building and Tree Shading Effects on the 3D Radiant Environment and Human Thermal Comfort in SummerwithinUrban Settlements inNortheast China

Shading is one of the most effective strategies to mitigate urban local-scale heat stress during summer. (erefore, this study investigates the effects of shading caused by buildings and trees via exhaustive field measurement research on urban outdoor 3D radiant environment and human thermal comfort. We analyzed the characteristics of micrometeorology and human thermal comfort at shaded areas, and compared the difference between building and tree shading effects as well as that between shaded and sunlit sites. (e results demonstrate that mean radiant temperature Tmrt (mean reduction values of 28.1°C for tree shading and 28.8°C for building shading) decreased considerably more than air temperature Ta (mean reduction values of 1.9°C for tree shading and 1.2°C for building shading) owing to shading; furthermore, the reduction effect of shading on UTCI synthesized the variation in the above two parameters. Within the shaded areas, short-wave radiant components (mean standardized values of 0.104 for tree shading and 0.087 for building shading) decreased considerably more than long-wave radiant components (mean standardized values of 0.848 for tree shading and 0.851 for building shading) owing to shading; the proportion of long-wave radiant flux densities absorbed by the reference standing person was high, leading to a relatively high long-wave mean radiant temperature, and R between long-wave mean radiant temperature and air temperature exceeded 0.8. Moreover, the directional sky view factor (SVF) was utilized in this study, and it showed significant positive correlation with short-wave radiant flux densities, but no statistically evident correlation with long-wave radiant flux densities. Meanwhile, Tmrt was most relevant with SVFS⟶ with R of 0.9756. Furthermore, UTCI rose two categories at the sunlit areas compared with that at the shaded areas. In contrast, Ta and Tmrt played the first positive role in UTCI at shaded and sunlit areas, respectively.


Introduction
Urbanization is accelerating worldwide and is affecting urban living and transforming human society in various manners [1]. Among them, increasing urban severe heat in summer is a problem for many countries, which has strongly aggravated human thermal stress within urban spaces [2]. e outdoor thermal environment is significantly affected by the design of the built and vegetation environment [3], and it affects the thermal comfort experienced by people; furthermore, people's perception for the thermal environment as a result influences their usage of outdoor spaces [4,5]. As outdoor leisure activities can improve people's health and vitality [6,7], thermally comfortable and environmentally healthy outdoor spaces in densely urban areas are of significance for citizens.
For the reasons outlined above, urban planning requires quantified information to improve outdoor space utilization for a better thermal environment and thermal comfort. Several studies indicated that shading can considerably contribute to a higher quality urban living in summer, particularly on a small scale [8,9]. So far, related studies have been conducted in a few countries and regions with various climates, such as Hong Kong, China [10][11][12]; Guangzhou, China [13]; Tempe, USA [14]; Freiburg, Germany [15]; Russi, Italy [16]; and Malaysia [17]. Local shading to block the solar radiation within urban areas can be accomplished by (1) optimized design of buildings, streets, and open spaces [18][19][20], (2) natural shading, such as trees [12,13,[21][22][23][24], or (3) man-made devices and structures, such as awnings or sunshades and tunnels [14,25]. Some studies have analyzed the shading effects of these methods on the outdoor local thermal environment. For instance, height-to-width ratio (H/W) is often used as an index to assess the shading level of simple building objects or street canyon geometry [20,26,27]. As for the tree shading method, researchers have applied leaf area index (LAI) to describe the occlusion of tree canopies to solar radiation [28]. Sky view factor is another important factor to weigh the shading degree of complex urban structures [11,29,30]. Shading is utilized to obstruct the incident solar radiation directly into urban spaces and further influences the surface temperature of the surrounding surfaces such as building exterior surface, underlying surface, and tree surface [11,23,28,31]. Additionally, near-surface air temperature and flow are affected because of convective heat transfer [15]. As a result, the radiant environment, air temperature, and humanbiometeorology vary with shading circumstances, with outdoor local-scale energy reduction and with space cooling [28]. Studies have demonstrated that, under normal circumstances, a higher H/W, lower SVF, or higher LAI can create a more obstructed shaded area [11] and thus decrease the heat stress. Urban outdoor shading can also effectively reduce the severe heat that people receive, thus considerably influencing people's thermal sensation and comfort. To evaluate people's perception of the surrounding thermal environment based on the energy balance between human body and outdoor environment, several thermophysiological indices are employed, such as physiological equivalent temperature (PET) [32], outdoor standard effective temperature (OUT_SET * ) [33], and Universal ermal Climate Index (UTCI) [34]. According to a study by Lee et al. [15], building and tree shade decreases PET by 13.1°C and 15.7°C, respectively, compared with the sunlit areas nearby. In summary, shading can be an effective method to ameliorate outdoor heat stress and improve thermal comfort during summer.
ere are a few limitations in the existing studies: (1) they focused on the shading effects on traditional outdoor micrometeorological parameters such as air temperature and seldom concentrated on the urban outdoor radiant environment; (2) few have analyzed the shading effects considering radiant flux densities in detail, for example, directional radiant components; and (3) most existing studies have focused on a specific shading structure in a particular region and season. As a result, the comparisons between the effects of diverse shading forms, for example, shading by neighboring buildings and trees, on the radiant environment and human comfort were generally neglected. e above subjects have not been widely investigated.
In this study, exhaustive measurements and an investigation that targeted different practical effects of building and tree shading forms on the outdoor thermal and 3D radiant environment and human thermal comfort during summer were conducted in a severely cold region. e research involved the subdivision and proportion of the composition of radiant flux densities that related to shortand long-wave radiation; the reduction influence of the two shading methods on detailed directional short-and longwave radiant flux densities from 3D environment and thus the roles of short-and long-wave radiant flux densities in T mrt ; the radiant environment features and different wavelength radiant flux densities functions at the shaded areas; and finally, the different shading effects on human thermal comfort. e research also used novel applications of the long-wave mean radiant temperature (LT mrt ) and directional SVF [11] to emphasize the importance of long-wave radiant flux densities under shaded conditions and 3D directional radiant characteristics, respectively. Moreover, the research subject of this study is in a severely cold region. Unlike studies that were conducted in warm or hot areas, the localscale ambient temperature in severely cold regions is not as high as in other climates, but the heat stress is still high during summer because of the high radiant flux densities outdoors. is indicates that the outdoor radiant environment may affect human thermal comfort differently in severely cold regions than in other regions. erefore, research subjects with the regional features of this study provide a necessary supplement to expand the research locations and results for related fields of study.
e main objectives of this study are as follows: (1) to explore the shading effects on micrometeorological parameters, particularly on the detailed directional 3D radiant components, and on human thermal comfort; (2) to investigate the radiant characteristics under shaded conditions; (3) to propose and apply directional SVF and evaluate its relationship with 3D radiant flux densities; and (4) to probe human thermal comfort differences and the impact of micrometeorological parameters on thermal comfort under shaded and sunlit conditions. is study offers a practical research method and quantified data targeting summer conditions for developing general and specific strategies of rational urban planning.  [35]. Summertime here is ephemeral, from early July to late August; however, the solar radiation is relatively high, and air temperature may exceed 30°C [36]. In this study, shaded areas obstructed by buildings (SHA_Buil) and trees (SHA_Tree) and a sunlit background meteorological station (SUN_0), in a university campus in the city center, were selected. Figure 1 shows the location of Harbin and the distribution of each research site in the university campus. In detail, SHA_Buil was occluded by 45 m height nine-story (on the east, south, and west side of the research spot) and 20 m height four-story (on the north side of the research spot) buildings, with light gray elastic coating exterior surfaces. SHA_Tree was under tall elm trees with a mean height of 10 m. Underlying surface material under both shaded sites was seepage brick, with light gray and red colored at SHA_Buil and light red colored at SHA_Tree. In addition, the playground where SUN_0 was on was covered by imitation turf and plastic track. SHA_Buil and SHA_Tree were selected strictly according to the height and degree of the enclosure that blocked most of the direct solar radiation from all directions. Quantitatively, the sky view factor with a lens facing upward was closed to or less than 0.3, and the lateral sky view factors were close to or less than 0.15. Moreover, they also met the requirements for denizens' outdoor activities. SUN_0 was on a playground that was sufficiently broad to represent an open space that received high radiation.

Methodology
Measurements at all the research sites were conducted from 7:30 to 18:00 within a week, from July 28 to August 3 in 2018, simultaneously. e meteorological conditions of short-term sample days should meet the characteristics of typical seasonal weather conditions [14,27,37,38]. For this study, the days for measurement were selected under clearsky conditions with relatively intense solar radiation in July and August, which are the hottest months of the year. Moreover, according to the official historical summer weather database over the years for Harbin (refer to the data for July and August from 2002 to 2018), the daytime ranges of air temperature and global radiant flux densities were 19.0-36.0°C and 14.0-1031.0 W/m 2 , respectively [39]. In this study, during the same period with the official weather data, the ranges of the measured air temperature and global radiant flux densities were 26.7-35.5°C and 20-1025 W/m 2 , respectively. is indicates that the meteorological conditions during the sample measurement period in this study were representative of the norm.

Field Measurement Scheme.
At SHA_Buil, SHA_Tree, and SUN_0, air temperature (T a ), relative humidity (RH), wind speed (v), and the directional 3D short-and long-wave radiant flux densities reaching the reference standing person (K i and L i , respectively) were recorded every 1 min at 1.1 m above ground level using three sets of tailor-made humanbiometeorological measurement units. Among these, 3D K i and L i were from the vertical (downward: ↓ and upward: ↑) and horizontal directions (easterly: E⟶, westerly: W⟶, southerly: S⟶, and northerly: N⟶). Surface temperature T s was recorded every 1 min using button-type temperature recorders that were fixed to the measurement positions using industrial adhesive. In detail, the specific locations of the T s recorders were as follows: five were used at SHA_Buil, among which one was pasted on the ground under the measurement unit, and four were fixed to the four surrounding building walls at 1.1 m above ground level right facing the measurement unit; another five recorders were used at SHA_Tree, and one was also pasted on the ground  Advances in Meteorology under the measurement unit; the other four were fixed to the four tree trunks facing the measurement unit with a distance of 1.1 m to the ground; for SUN_0, the surrounding environment was much more open, so only one recorder was fixed to the ground under the measurement unit. Generally, 4410 data points were obtained for each parameter measured above. Moreover, directional sky view factor SVF i values were obtained through fisheye photos taken adjacent to each radiometer and calculated by RayMan software [40,41]. In detail, SVF was subdivided into five directions to represent the urban structures: lens facing upward for SVF↓ and facing lateral directions for SVF E⟶ , SVF W⟶ , SVF S⟶ , and SVF N⟶ . Table 1 lists the main technical parameters of the sensors used in this study.
Because RH depends on T a to a certain degree, it makes more sense under a human-biometeorological perspective to use the water vapor pressure (VP) [30]. erefore, in the following content, we converted RH to VP by multiplying RH and saturated water vapor pressure corresponding to various T a . To describe and compare the data under different shaded or sunlit conditions, the seven-day mean data and standard deviation (SD) were used as statistical indices to display the overall distribution, centralization tendency, and dispersion degree of the local-scale micrometeorological variables. Figure 2 shows the mean and SD distributions of T a , VP, global radiant flux densities G (equal to K↓ in this study), and v from all the research sites during the measurement period. e figure illustrates that there was a relatively large difference in micrometeorological conditions among the shaded and sunlit areas, as well as between different shading forms. e values and discrete degrees of the four types of micrometeorological parameters at the shaded sites were significantly lower than those at the sunlit background meteorological station, particularly for G, which indicated that both building and tree shading can effectively reduce outdoor radiation. Figure 3 shows the overall distribution, centralization tendency, and dispersion degree of the surface temperature T s of the exterior walls, tree trunks, and underlying grounds. As can be seen, T s at the shaded areas was evidently lower than that at the sunlit area. Mean ground surface temperature at SHA_Buil and SHA_Tree was 10.9°C and 11.4°C lower than that at SUN_0, respectively. Comparing the two shaded areas, T s of surfaces facing the same direction at near SHA_Buil was higher than that at near SHA_Tree, with 1.1°C mean facade temperature difference and 0.5°C ground surface temperature difference; furthermore, it fluctuated more strongly at SHA_Buil than at SHA_Tree comparatively.

Mean Radiant
Temperature. According to the radiation budget between the radiant environment and the reference standing person, the "six-directional method" was utilized to calculate T mrt in this study [43], as shown in the following equation: where K abs and L abs are the total of short-and long-wave radiant flux densities absorbed by the standing humanbiometeorological reference person [44], respectively (both in W/m 2 ). Parameter α l is the long-wave absorption coefficient, σ is the Stefan-Boltzmann constant (5.67 * 10 −8 W/ (m 2 ·K 4 )), and the unit of T mrt is°C. K abs and L abs are calculated using the following equations, respectively: where α k is the absorption coefficient of human body for short-wave radiant flux densities K i (W/m 2 ) and α l is the absorption coefficient of human body for long-wave radiant flux densities L i (W/m 2 ). e standard value of the two physical quantities above is as follows: α k � 0.70 and α l � 0.97. W i are the angular factors of the reference standing person for K i and L i and are set to 0.06 for the two vertical directions as well as 0.22 for the four horizontal directions [43]. e total short-or long-wave radiant flux densities absorbed by the reference standing person can be classified into the vertical and horizontal ones (K ver,abs , K hor,abs , L ver,abs , and L hor,abs ) by adding the corresponding components.
e main purposes of using T mrt in this study are summarized as follows: (1) to comprehensively quantify the short-and long-wave radiation from the 3D environment and to evaluate the specific shading effects on the radiant environment; (2) to determine the thermal comfort index UTCI and to assess its role in human thermal comfort under different shaded and sunlit conditions.
In this study, long-wave mean radiant temperature (LT mrt ) was inducted to represent the impact of long-wave radiant flux densities at the shaded areas by substituting L abs for K abs + L abs , as shown in the following equation [11]:

Standardization Processing.
In this study, the data were presented by means of difference related to the micrometeorological parameters and thermal comfort index such as T a , T mrt , and UTCI, and of ratio for the radiant flux densities (short-and long-wave radiant flux densities), instead of absolute values. e impact of different shading forms on the micrometeorological parameters, thermal comfort index, and radiant flux densities was then compared. Difference of T a , T mrt , or UTCI was calculated as follows: where ΔT represents the T a , T mrt , or UTCI difference between the building or tree shading and the sunlit background areas (ΔT a , ΔT mrt , or ΔUTCI); T SHA_Buil and T SHA_Tree are T a , T mrt , or UTCI at the building and tree shading areas, respectively; and T SUN_0 is T a , T mrt , or UTCI at the sunlit background meteorological station. With regard to the radiant flux densities, ratios were utilized, as shown below: where rK and rL are the standardized short-and long-wave radiant ratios of the shaded to the sunlit background areas, respectively; K SHA_Buil and K SHA_Tree are the short-wave radiant flux densities at the building and tree shading areas, respectively; L SHA_Buil and L SHA_Tree are the long-wave radiant flux densities at the building and tree shading areas, respectively; and K↓ SUN_0 and L↑ SUN_0 are the reference short-and long-wave radiant flux densities at the sunlit background meteorological station, respectively (downward ↓ short-wave radiant flux densities and upward ↑ long-wave radiant flux densities).

Overall Shading Effects by Neighboring Buildings and
Trees. Based on the abovementioned calculation processing, the shading effects of neighboring buildings and trees on the micrometeorological parameters including T a , the shortand long-wave radiant components, and thus T mrt , and the human thermal comfort evaluation index UTCI are analyzed in this section.
3.1.1. Shading Effects on T a . 29 July 2018 was selected as the typical weather day. Figure 4(a) shows the temporal variation in T a at the shaded and sunlit sites. Different T a patterns and magnitudes appeared among the three areas that peak T a was observed at noon at SUN_0 and SHA_Buil, but T a at SHA_Tree peaked in the morning, then fluctuated at noon and early afternoon, and dropped in the late afternoon and dusk. By applying equation (5), the effect of shading by neighboring buildings and trees on T a can be analyzed, as presented in Figure 4(b). Tree shading brought about a stronger decrease in T a especially in the afternoon, with a mean reduction value of 1.9°C and a peak value of −2.8°C at about 12:10; while building shading caused T a to decrease relatively smoothly with a mean reduction value of 1.2°C and a peak value of −2.0°C at 12:10. e main reason of shading's decreasing effect on T a is that shading can occlude a considerable part of incident solar radiation. is weakened the radiant exchange and decreased the surface temperature of the surrounding solid surfaces (see Figure 3), thus lowering the convective heat transfer between the exterior surfaces and local air. Moreover, the evapotranspiration cooling effect of the vegetation [21][22][23]45], as well as the more enclosed environment (involving more occlusion to incident solar radiation), around SHA_Tree reduced T a more than SHA_Buil.

Shading Effects on 3D Radiant
Environment. e diurnal temporal variations of short-and long-wave radiant components reaching the reference standing person are plotted in Figures 5(a) and 5(b). It can be seen that the shortwave radiant component's variation ranges remained steady at a lower level between 27 and 107 W/m 2 at SHA_Buil regardless of the solar position. In contrast, at SHA_Tree, short-wave radiation fluctuated during the pre-and postperiods of high solar altitude angles from the incident directions (↓, S, and E in the morning and W in the afternoon), respectively, because of solar radiation transmission through the leaves. e highest values occurred at SUN_0, especially for K↓ and K S⟶ with the peak value reaching 1025 W/m 2 when the solar altitude angle was the highest during the day. Besides, K E⟶ and K W⟶ were symmetrically reversed at approximately the highest solar altitude angle, and they reached their peak values in the morning and afternoon, respectively. Figure 5(b) presents the changing profiles of the long-wave radiant components. It can be seen that the patterns of the long-wave radiant components and the relative magnitudes among the research sites were similar to those of the local-scale T a , indicating that the temporal variation in T a could affect the long-wave radiant components. ere might exist a certain correlation between longwave radiation and T a (it will be further analyzed in Section 3.2.2). In addition, directional differences among the components at the two shaded areas were smaller than those at SUN_0. Comparing the two shaded areas, the long-wave radiant flux densities were lower at SHA_Tree than at SHA_Buil. As mentioned in Section 3.1.1, T a could be lowered near greenery. erefore, it is well-founded that the long-wave radiation could also be lower, as T a was, because of the effect of the evapotranspiration of the trees around SHA_Tree.
We also used equations (6) and (7) to compare the different summer shading effects on outdoor 3D short-and long-wave radiant flux densities reaching the reference standing person on July 29, as shown in Figures 5(c) and 5(d). For rK i in Figure 5(c), before 17:00, the reduction effect of building shading on short-wave radiation had much weaker fluctuations than that of tree shading. In particular, rK at SHA_Buil remained below 0.15 until late afternoon; tree shading led to a higher standardized short-wave radiation in the short-wave incident directions from ↓, S, and E in the morning and from W in the early afternoon, which was mainly caused by the transmission of solar radiation under the tree shading condition. After 17:00, shading effects on reducing the short-wave radiation began to decrease at both shaded areas.
As delivered in Figure 5(d), the directional standardized long-wave radiation under both shaded conditions presented different patterns compared with standardized short-wave radiation that they reduced from the beginning of the measurement, reached their lowest values close to 0.75 at noon, and increased till the measurement ended. By contrast, the standardized long-wave radiation variations from different directions were larger at SHA_Buil than at SHA_Tree. A special circumstance appeared from 9:30 to 10:50 when the downward standardized long-wave radiation had a mutation, which corresponded to the time at which rK showed an abrupt change at SHA_Tree because of the radiant transmission effects. Figure 6(a), the temporal trends of T mrt were primarily similar to those of the short-wave radiant flux densities from the main incident directions (downward and southerly), which was caused by the significant correlation between T mrt and K↓ and K S⟶ . It has been proved by some studies [30,35]. Moreover, when taking T a in the foregoing content into account, T mrt was much higher than T a with a mean difference (T mrt -T a ) of 13.5°C at SHA_Buil, 15.5°C at SHA_Tree, and 40.7°C at SUN_0 during the daytime. Figure 6(b) shows the shading effects by neighboring buildings and trees on T mrt on July 29 using equation (5). As can be seen in the figure, building shade reduced T mrt more significantly than tree shade did before 11:00; however, after this time, the two shading forms showed similar reduction effects on T mrt . Meanwhile, shading could significantly reduce T mrt with mean values of 28.8°C for building shading and 28.1°C for tree shading 6 Advances in Meteorology      Advances in Meteorology compared with the sunlit background meteorological station. According to the results discussed in Section 3.1.2, summer shading had a stronger impact on short-wave radiant components than on long-wave radiant components. erefore, T mrt reduction by shading was mainly caused by the decrease in short-wave radiant flux densities. Furthermore, T mrt reduction by shading was stronger than the T a reduction (Figure 4), revealing that shading can lead to a more significant effect on outdoor radiant environment than on air temperature.  Fiala multinode model [34]. It is appropriate for thermal assessments on any scale as well as in all climates and seasons [34,46]. In this study, we calculated UTCI by inputting T a , T mrt -T a , and VP at 1.1 m above ground level and v at 10 m above ground level (extrapolated from 1.1 m to 10 m above ground level by equation (8) [47]) for the shaded and sunlit areas ( [46], http://www.utci.org) based on the adaptive clothing model [48] and reference human activity level of walking at a speed of 4 km/h [46]. To quantify the shading impact on UTCI, the differences between the building or tree shading and the sunlit background meteorological station on July 29 were calculated using equation (5), as shown in Figure 7. Overall, tree shading brought about a stronger decrease in UTCI than building shading did and that shading decreased UTCI with mean values of 9.2°C for tree shading condition and 7.8°C for building shading condition. Moreover, both shading conditions reduced UTCI most strongly at noon. An opposite pattern appeared from 9:30 to 10:50 where ∆UTCI was higher at SHA_Tree than at SHA_Buil, the period of which was corresponding to the radiant flux densities' sudden increase at SHA_Tree (see Section 3.1.2). e shading reduction effect on UTCI depended on its comprehensive impact on the micrometeorological parameters, especially on T a and T mrt . As a result, ∆UTCI integrated the characteristics of ∆T a along with ∆T mrt :

Shading
where v 10 is the wind speed at a height of 10 m above ground level (m/s), v 1.1 is the wind speed measured by the sensors at 1.1 m above ground level in this study (m/s), z is the distance from the ground (10 m in this case), α is the mean speed exponent set to 0.33 in city-center areas, and z' is the height of the sensors installed above ground level (1.1 m in this study).
To sum up, both building and tree shading had considerable reduction effects on the radiant components and thus on T mrt , which was greater than that on T a . eir influence on UTCI consequently integrated the characteristics of ∆T a along with ∆T mrt . Meanwhile, the differences between the two types of shading forms were nonnegligible.

Radiant Characteristics under Shaded Conditions.
As the radiant environment is an important factor in the summer shading function, detailed analyses considering directional short-and long-wave radiant components are conducted in this section. In particular, the different influences and contributions of short-and long-wave radiation are compared.     10 Advances in Meteorology based on the seven-day measurement data. From Figure 8(a), we can see that the variation ranges of K i,abs under the tree canopies were higher than those under the neighboring building obstacles because of the greater shortwave radiant flux densities transmission through the leaves under the tree shading condition. However, values of L i,abs were higher at SHA_Buil than those at SHA_Tree, as presented in Figure 8(b) is was mainly because there existed massed solid surfaces at the building shading area that received more long-wave radiant flux densities. Moreover, directional differences among the two vertical directions and the four horizontal directions were distinct in both L i,abs and K i,abs components; the average values of short-or long-wave radiant flux densities absorbed by the reference standing person from a single horizontal direction were three to four times those of a single vertical direction. is also resulted in the mean K hor,abs or L hor,abs being six to eight times greater than that of K ver,abs or L ver,abs . ese results mainly attributed to the larger emission intensity from the horizontal directions and the higher absorption for the horizontal short-or long-wave radiant flux densities by the reference standing person of 0.22 versus for the vertical ones of 0.06. In urban outdoor settlements, shaded conditions do not mean the complete absence of short-wave radiation during the daytime. ere still exist diffuse and a small amount of direct short-wave radiation. Some studies have proved that the proportion of long-wave radiant flux densities absorbed by the human body was relatively high in the outdoor environments [14,30,35]. is study further calculated and compared the proportion of K abs and L abs in K abs + L abs under different shaded and sunlit conditions. Figure 9 presents K abs /(K abs + L abs ) and L abs /(K abs + L abs ) at the shaded areas in this study. K abs /(K abs + L abs ) did not exceed 10% under the shaded condition by neighboring buildings, which implied that L abs /(K abs + L abs ) reached more than 90%. e proportion range of K abs in K abs + L abs under the tree shading condition was larger with a peak value of 23%, indicating a minimum L abs /(K abs + L abs ) value of 77%. By comparison, L abs /(K abs + L abs ) at the sunlit area was much lower, with a minimum value even below 65% and a maximum value below 88%. erefore, shading not only reduced the air temperature, 3D radiant flux densities, and T mrt , but also caused a higher portion of L abs in K abs + L abs under shaded conditions at the same time. Similar studies were carried out in summer by Middel and Krayenhoff [14] and Lee et al. [30]. e former distinguished the decomposed contributions of directional short-and long-wave radiation on T mrt at a series of sites during one day. ey concluded that there was a higher proportion of lateral and hence total long-wave radiation at shaded areas than at sunlit areas. e latter considered clear-sky summer days from 2007 to 2010 and did not differentiate between shaded and sunlit conditions but also pointed out that the

Advances in Meteorology 11
proportion of K abs did not exceed 40% with L abs proportion not falling below 60%. Because summer shading can greatly increase L abs / (K abs + L abs ), we applied LT mrt by substituting L abs for K abs + L abs under shaded conditions (see equation (4)). LT mrt was the degenerated T mrt considering only long-wave radiant flux densities absorbed by the reference standing person [11]. It also reacted to the mean surface temperature of surrounding objects within shaded areas where the longwave radiant flux densities played a more essential role. Figure 10 shows the distributions of LT mrt at the two shaded research sites. As can be seen, LT mrt at SHA_Buil fluctuated more strongly (with a peak of 41.8°C and a mean value of 39.6°C) than that at SHA_Tree (with a peak of 40.4°C and a mean value of 38.8°C).
is again confirmed the results shown in Figure 8(b) that the L i,abs values and variation range were higher at SHA_Buil than those at SHA_Tree.

Relationship between LT mrt and T a .
e outdoor radiant environment is the decisive factor for urban micrometeorological and human-biometeorological conditions. First, as shown in the foregoing section, shading can reduce T a slightly. As the atmosphere is almost transparent to solar radiation (short-wave), it is very weak for the atmosphere to increase air temperature by receiving solar radiation directly. In general, convective heat transfer between the local atmosphere and the solid surfaces (including underlying surfaces, neighboring building exterior walls, and vegetation surfaces) is the principal factor causing T a to change, in which the surface temperature of the solid surfaces plays a direct role. e surface temperature varies with the incident solar radiation, which affects the radiation exchange. By taking the surface temperature at the same position (the underlying ground surface) of the three research sites for instance (Figure 3), the mean T s at the shaded research sites was approximately 11.1°C lower than that at the sunlit sites, thus reducing T a to a certain extent. As a result, shading not only decreased radiation intuitively but also decreased T a indirectly.
LT mrt reflects the long-wave radiant flux densities emitted from the surrounding solid surfaces. erefore, the surface temperature has a higher explanatory power for the outdoor long-wave radiant flux densities emissions under shaded conditions [11]. As proved in the previous paragraph, T a was significantly influenced by the surface temperature. us, we analyzed the relationship between these two affected factors influenced by surface temperature of the shaded areas.
First of all, the time during the measurement was divided into four periods: morning (7:30-10:00), noon (10:00-13: 00), afternoon (13:00-16:00), and dusk (16:00-18:00) according to the solar altitude angle. As shown in Figure 11, the solar altitude angle range was 6-60°d uring the research period. 30°and 50°were selected as the degree nodes to classify, and their nearest sharp-time points were set as the division bases. en, the analysis of variance was conducted at SHA_Buil and SHA_Tree taking LT mrt , T a , and time into account. Finally, the estimated marginal mean results at the two shaded areas are obtained, as shown in Figure 12.
It can be seen that LT mrt varied with T a and elapsed time. On the one hand, the time segments for LT mrt and T a to reach their highest values were different at the two shaded sites. In detail, the highest LT mrt with the corresponding T a appeared in the afternoon and at dusk at SHA_Buil; however, LT mrt and T a were higher at noon and in the afternoon than at other periods at SHA_Tree. It can be interpreted that, as time lapsed, the heat accumulation increased till afternoon and dusk at SHA_Buil, where there were numerous structures with strong heat storage capacity. However, because of the transmission at SHA_Tree, radiation from various directions from the external environment was allowed to enter the research space to heat solid surfaces such as trunks. As is known, radiation is strongest around midday, which directly causes LT mrt and the corresponding T a during noon and afternoon being higher than those in the morning and at dusk.
On the other hand, from the trend in the scatter plots between the estimated marginal mean LT mrt and T a , linear increase in the estimated marginal mean LT mrt with T a was found at both SHA_Buil and SHA_Tree. erefore, Table 2 summarizes their relationship that is quantified by the coefficient of determination R 2 (LT mrt and T a at the two shaded areas satisfied the normal distribution), based on the data throughout the measurement period. Evidently, the correlation coefficients were intensely significant at the significance level sig � 0.000 regardless of SHA_Buil, SHA_Tree, and the total values. erefore, LT mrt had a significant positive correlation with T a , in which the correlation was stronger at SHA_Buil. Meanwhile, the difference in the analyses of variance and regression results between the building and tree shading was attributed to the urban microclimatic factors owing to the neighborhood: the organization, disposition, and material of the surrounding environment. To conclude, T a can offer a basic statistical comprehension to a local-scale LT mrt at shaded areas.
From the results in this section, we discover that summer shading also led to variations in the short-and long-wave proportions. LT mrt varied with T a and elapsed time under different shaded conditions. Yet how the neighboring building and tree obstacles, that are ascribed to urban morphology, block radiation and quantitatively influence the radiant components should be investigated further, as presented in the following section.

Quantified Directional Urban Morphology and Radiant
Environment.
e shading effects on urban local-scale microclimate are mainly caused by the blocking of direct short-wave radiation. As mentioned in Section 3.1.2, shading primarily influenced short-wave radiant flux densities from the incident directions of the direct radiation to the greatest extent. Furthermore, shading of the direction-dependent short-wave radiant flux densities had consequences for other radiant flux densities. Ultimately, different shading modes mainly affect the short-and long-wave radiant flux densities in the spatial distribution and range. eir temporal variations primarily change with the solar altitude and orientation. Moreover, the urban morphology reflected by SVF of the research site plays an essential role in the microclimate. us, directional SVF was proposed to examine the dependence of microclimate on the urban structure in terms of different shading methods in this study [11]. Table 3 summarizes the SVF i from five directions at the research sites. From the SVF images at SHA_Tree, we can infer that it offered a higher probability of short-wave radiant flux densities entering the space from more tilt angles. To conclude, although the directional SVF i values at SHA_Tree were lower than those at SHA_Buil in the corresponding direction, the short-wave radiant flux densities reaching the tree shading site still maintained a higher level and a stronger fluctuation (see Figure 8(a)).
To analyze the detailed relationship between the radiant components and SVF i , five directional short-wave (longwave) averaged radiant data at each research site were calculated over the seven measurement days and were associated with the five corresponding directional SVF i values [11,14]. erefore, fifteen data points in total were linearly fitted for the directional short-or long-wave and the corresponding SVF i with the same direction, as presented in Table 2: Linear fitting analysis between estimated marginal mean LT mrt and T a at the building and tree shading areas and the total.  Figure 13. Evidently, directional SVF i caused the short-wave radiant flux densities to linearly increase and the long-wave radiant components to decrease. However, no clear dependence of long-wave radiant flux densities on SVF i was found because of the exceedingly low R 2 . erefore, SVF had a higher explanatory power for the short-wave radiant flux densities, but this was not evident for the long-wave radiant flux densities. We further investigated the impact of the directional SVF i on T mrt . e nondirectional T mrt of the three research sites was averaged over the seven measurement days [14,15,20,30,49] and was linearly fitted with the five directional SVF i and the mean SVF, with three data points around one fitting line, as shown in Figure 14. Table 4 summarizes R 2 values, indicating that T mrt was strongly correlated with the mean SVF, in which T mrt was most relevant to SVF S⟶ . Figure 14 simultaneously shows that the slope values were all positive, which once again proved that T mrt was mainly governed by short-wave radiant flux densities, that shading in terms of a lower SVF brought out a lower K abs + L abs even if the long-wave radiant flux densities were increasing.
Similarly, Middel and Krayenhoff also demonstrated that the lateral short-wave radiation correlated more strongly with SVF↓ (with R 2 of 0.78 at noon and 0.70 in the afternoon) than the lateral long-wave radiation did (with R 2 of 0.02 at noon and 0.08 in the afternoon). Besides, T mrt was found to have R 2 values of 0.66 at noon and 0.64 in the afternoon related to SVF↓ during the daytime [14]. In brief, it has turned out that directional SVF i values were meaningful for the study of the effects of shading on outdoor directional radiant components and the comprehensive radiant parameter T mrt .

Human ermal Comfort (Distributions and Influencing Micrometeorological Parameters).
e effects of shading on micrometeorological parameters such as radiant components and air temperature have been discussed in the previous sections. e outdoor shaded or sunlit conditions would affect the denizens' thermal experience and thermal comfort. erefore, further analyses of the impact of different shading forms on human thermal comfort will be explored in this section.

UTCI Distribution.
Jin et al. conducted the outdoor meteorological measurements and thermal sensation and comfort questionnaire survey also in Harbin and determined the ranges of UTCI that represent people's feeling of comfort based on the field research [50]. In this study, we utilized the UTCI values categorized in terms of thermal stress reclassified by them (Table 5). Seven-day UTCI values were separated into different thermal stress categories, as shown in Figure 15. At the shaded areas, UTCI was primarily distributed over "strong heat stress" and "moderate heat stress," in which the proportion of "moderate heat stress" at SHA_Tree was 22% higher than that at SHA_Buil. By contrast, UTCI at the sunlit areas was raised by two categories under similar background weather conditions that the thermal stress was categorized as "extreme heat stress," "very strong heat stress," and "strong heat stress." It indicates that the outdoor heat stress during summer daytime was at a high level, which also reflected the necessity of summer shading to improve human thermal comfort outdoors.

Impact of Micrometeorological Parameters on Human
ermal Comfort. Urban design is interested in humanbiometeorological results, whose application contributes to the citizens' thermal comfort, health, outdoor activities, work efficiency, and well-being [51]. Based on all the sevenday field measurement data at shaded and sunlit areas, the partial correlation analysis results in Table 6 demonstrated that UTCI was positively determined by T a , T mrt , and VP but negatively depended on v on summer days. It can be seen that both T a and T mrt were highly correlated with UTCI, with the partial correlation coefficient r values exceeding 0.96. Nevertheless, subtle differences appeared when distinguishing shaded and sunlit sites, that UTCI firstly depended on T a at the shaded areas, followed by T mrt : but conversely, T mrt played the most important role in UTCI at the sunlit areas where the radiant flux densities were significantly high and denizens experienced strong heat stress.
e reason can be explained as follows. e radiation exchange, quantified by T mrt , is the key meteorological variable in summer. At shaded areas in the daytime, there is usually a strong correlation between T mrt and T a [52]. Meanwhile, there is less incident short-wave radiation entering these enclosed spaces. is is the reason why we found that T a was the first positive variable for the human thermophysiological   Figure 14: Linear fitting between seven-day mean T mrt and the directional SVF as well as the mean SVF. perception to the outdoor environment. In other words, T a had the role of a proxy variable for the radiation exchange under shaded condition. While at the sunlit area, T mrt was still the dominant factor for UTCI. e above outcomes were also illustrated by Lee et al., who indicated that the correlation between human thermal comfort and T mrt was slightly stronger than with T a on clearsky summer days [15,30]. In comparison, Middel and Krayenhoff conducted day-night measurements in 22 sunlit or shaded sites on the hottest day in Tempe and found that UTCI was most sensitive to T a , whereas it was about half as sensitive to T mrt as T a [14]. e differences of the results between their research and this study are mainly because (1) their study considered not only the daytime but also the evening and night after sunset when solar radiation passed off and T a was dominant and (2) they focused on all the sunlit areas together with the shaded sites, so the relationship between UTCI and each individual micrometeorological variable under shaded and sunlit conditions was not identified. e relatively low partial correlation coefficient values between UTCI and VP as well as v indicated that these two microclimatic parameters at the shaded and sunlit areas played a more subsidiary part in the perception of outdoor heat by denizens, which was also proved by Middel and Krayenhoff using sensitivity analysis [14].
So far, human thermal comfort within urban settlements affected by neighboring building and tree shading has been taken into account by implementing the "human-oriented" concept in urban planning and design objectives.

Conclusions
is study investigates the summer shading effects caused by neighboring buildings and trees and analyzes the difference between shaded and sunlit conditions as well as that between different shading forms in a typical city in the severely cold region. Summer shading can significantly reduce outdoor heat stress and offer the diversity in thermal environment and spaces for denizens. Based on the analyses, we obtained the following conclusions: (1) Shading led to a stronger reduction on T mrt than on T a , and its effect on UTCI synthesized the variation in the above two parameters that building shading decreased T mrt , T a , and UTCI with mean values of 28.8, 1.2, and 7.8°C and tree shading of 28.1, 1.9, and 9.2°C, respectively, compared with the background meteorological site with sunlit condition. (2) Within the shaded areas, short-wave radiant components decreased considerably more than longwave radiant components owing to shading; the proportion of L abs in K abs + L abs was high, and it led to a relatively high LT mrt , which had an R 2 with T a exceeding 0.8, and reached their highest values in the afternoon and at dusk at SHA_Buil and at noon and in the afternoon at SHA_Tree. (3) Directional SVF i exhibited a significantly positive correlation with short-wave radiant flux densities; however, the relationship between SVF i and 3D long-wave radiant flux densities was not statistically evident. Meanwhile, T mrt was most relevant with SVF S⟶ with a R 2 of 0.9756. (4) UTCI at the sunlit areas rose two categories compared with that at the shaded areas under the similar background weather conditions; T a and T mrt played the first positive part in UTCI under shaded and sunlit conditions, respectively. (5) Future research will aim to conduct long-term measurements considering urban morphology variations. and to obtain more novel insights.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.