Evaluating the Effects of Façade Greening on Human Bioclimate in a Complex Urban Environment

The evaluation of the effectiveness of countermeasures for a reduction of urban heat stress, such as façade greening, is challenging due to lacking transferability of results from one location to another. Furthermore, complex variables such as the mean radiant temperature (Tmrt) are necessary to assess outdoor human bioclimate. We observed Tmrt in front of a building façade in Berlin, Germany, which is half-greened while the other part is bare. Tmrt was reduced (mean 2K) in front of the greened compared to the bare façade. To overcome observational shortcomings, we applied the microscale models ENVI-met, RayMan, and SOLWEIG. We evaluated these models based on observations. Our results show that Tmrt (MD = −1.93 K) and downward short-wave radiation (MD = 14.39W/m) were sufficiently simulated in contrast to upward short-wave and long-wave radiation. Finally, we compare the simulated reduction of Tmrt with the observed one in front of the façade greening, showing that the models were not able to simulate the effects of façade greening with the applied settings. Our results reveal that façade greening contributes only slightly to a reduction of heat stress in front of building façades.


Introduction
Heat stress risk in cities threatens human health [1] and effective countermeasures are not clearly identified.Even though many measures to reduce outdoor urban heat stress are proposed from various disciplines and at different spatial scales [2], their effectiveness is still disputed.Assessing the effectiveness of countermeasures is challenging for many reasons.For example, the transferability of results from one location to another is limited due to the complexity of the urban system and various climates [3].Moreover, the impact of a countermeasure on outdoor human bioclimate cannot be described sufficiently by simple climate elements, such as surface or air temperature.The mean radiant temperature ( mrt ) is an important variable for the assessment of human bioclimate in urban environments because it includes long-wave and short-wave radiation that reaches the human body [4,5].
Fac ¸ade greening is a promising countermeasure to reduce urban heat.It can be attached to a large area in cities [6] and features several cobenefits, such as insulating buildings [7][8][9] or serving as habitat for wildlife [6].The knowledge about the effectiveness of fac ¸ade greening to reduce outdoor heat stress risks, however, is still incomplete due to limits in transferring of results to other sites [10].The cooling effects of fac ¸ade greening regarding surface and air temperature depend on solar irradiance, vegetation properties, and the particular greening system [11][12][13].Also, the mechanisms of cooling (e.g., through shadowing or transpiration) vary between different plant species [13].The effects of fac ¸ade greening on air temperature were found to be small to negligible [7,11,14].Wall temperatures decreased more strongly [14,15] as well as the emitted long-wave radiation [16]. mrt was reduced by 2 K to 13 K in the tropics [15], but the effects for other locations are unknown.Studies that quantify the influence of fac ¸ade greening on  mrt are rare [15] because so far architecture and engineering disciplines are dominant with a focus on air and surface temperature and its influence on the building [10].
Microclimate models can be applied to expand the knowledge of the effects of fac ¸ade greening on human bioclimate.Microclimate models, such as ENVI-met, RayMan, or 2 Advances in Meteorology SOLWEIG, are often used in different disciplines.Particularly over the last few years, the number of studies applying ENVImet and RayMan has grown rapidly.This is an encouraging trend as it may present an increasing awareness of the topic of human bioclimate in cities and may accelerate the implementation of countermeasures.Nevertheless, the reliability and the uncertainty of the results from studies solely based on simulations may vary appreciably.The deviations between observations and simulations have sparsely been evaluated in a comprehensive way.Thus, the inaccuracies in simulating countermeasures or real case situations are insufficiently known even apart from the special case of fac ¸ade greening.A comparison of  mrt simulated by ENVI-met, RayMan, and SOLWEIG under the same conditions was so far performed by Chen et al. [17].Other variables, such as short-wave and long-wave radiation, did not undergo an intercomparison [36].
So far the models have not been applied to simulate fac ¸ade greening and we are not aware of other models that are able to simulate fac ¸ade greening with regard to outdoor human bioclimate.ENVI-met, RayMan, and SOLWEIG include plants, but specific vegetation types, such as fac ¸ade greening, are not explicitly declared to be supported.Nevertheless, all of the models should be able to consider some aspects of fac ¸ade greening.RayMan supports different trees with a specific emissivity and albedo as well as changes due to vegetation in the sky view factor (SVF).Thus, albedo and emissivity are modified due to fac ¸ade greening.SOLWEIG considers vegetation in SVF and reduction of short-wave radiation through transmissivity and sets sunlit leaf temperatures to air temperature.Hence, especially long-wave radiation is expected to alternate in front of a building with fac ¸ade greening.ENVI-met moreover simulates evapotranspiration.Thus, several effects of fac ¸ade greening are expected to be reproduced, such as changes in SVF, increased water-vapor fluxes, decreased long-wave radiation, or modified shortwave radiation.
In order to study the effects of fac ¸ade greening on outdoor human bioclimate, we will firstly use observational data to answer the following question: (1) How large is the reduction of  mrt in front of a greened fac ¸ade compared to a bare one at a study site?Afterwards, we will apply the models ENVImet, RayMan, and SOLWEIG to the same site in order to (2) evaluate the general performance of the models in simulating  mrt and other variables relevant to assess human bioclimate.Thereby, we will contribute to the intercomparison and evaluation of these models in a complex urban environment.Finally, we investigate (3) if the models are able to simulate the observed alteration of  mrt in front of the fac ¸ade greening.In a wider sense the last question addresses if microclimate models are able to represent specific types of urban vegetation such as fac ¸ade greening.

Study Site.
In order to observe and simulate the effects of fac ¸ade greening on  mrt , we chose a south-southwest oriented building facade, which is half-covered with Parthenocissus tricuspidata while the other part is bare.The fac ¸ade is located at the campus of the Berlin Institute of Technology, Germany, Hardenbergstaße 38 (52.5112N, 13.3241E, 31-35 m above sea level), in a compact midrise zone (Local Climate Zone 2) [37].Parthenocissus tricuspidata covers the building wall from the bottom to the top (9 m) and nearly half of the building width (11 m).The plants are attaching themselves to the fac ¸ade without technical climbing support.They rooted in a raised bed on the west oriented side of the building (unsealed area about 6 m 2 ).The substrate consists of loamy sand.The plants are irrigated in irregular intervals, but we irrigated the cavity before the measurement campaign to guarantee sufficient water supply.The albedo (calculated from observations on 7, 23, 2013, 1000-1600 Central European Time, CET) is 0.36 in front of the bare building wall, 0.28 in front of the fac ¸ade greening, and 0.18 for the ground surface.The mean depth of the vegetation layer is 0.27 ± 0.08 m (Table 1).The average leaf area density (LAD) of the fac ¸ade greening is 1.85 m 2 /m 3 according to our measurements at one vertical transect of 1 m × 9 m.On 8, 19, 2013, we harvested all leaves at this vertical transect and measured the depth of the vegetation layer and the size of all leaves with an area meter (Model 3100, LI-COR, Inc.).

Observation.
We measured on 7, 23, 2013 (00-23 CET), in front of the building fac ¸ade.This day was characterized by clear sky conditions and high air temperatures (Figure 1).The measurements comprised mobile integral radiation measurements as well as wind speed, relative humidity, air temperature, and several other variables to initialize ENVImet, such as soil and indoor temperature (Figure 2).
The mobile measurements quantified  mrt using integral radiation measurements based on three ventilated net radiometers (CNR4, Kipp & Zonen Corp., accuracy ±10% for daily totals) and the calculation of angular factors with the formula by Thorsson et al. [4] ((1) and ( 2)).Equation (1) summarizes the long-and short-wave radiations from the bottom, the top, and the four cardinal points.The variables were weighted according to the angular factors in order to estimate the shape of a standing person.Afterward  mrt was calculated from mean radiant flux density ( str ) with the Stefan-Boltzmann law (2).The net radiometers were oriented along the fac ¸ade, which means towards north-northeast (NNE) and not towards the main cardinal points: where   = angular factors for weighting the radiation fluxes (0.06 for up-and downward and 0.22 for the cardinal points),   = short-wave radiation (W/m 2 ),   = long-wave radiation (W/m 2 ),   = absorption coefficient of short-wave radiation (0.07),   = the emissivity of the human body (0.97), and  = the Stefan-Boltzmann constant (5.67⋅10 −8 W/m 2 K −4 ).The mobile station also included a ventilated temperature and relative humidity sensor (CS215, Campbell Scientific Inc., accuracy ±0.4 ∘ C for +5 to +40 ∘ C and ±4% for 0 to 100%).The station was set up at a distance of 1.2 m to the building.The measurement height was 1.1 m.We moved the station every 15 min between the greened and the bare site and aggregated the observed data afterwards from a minutely to an hourly resolution for each site (Figure 2).
Additional measurements especially for initializing ENVI-met comprised soil temperatures at a depth of 0.3 m, 0.2 m, and 0.1 m (analog mercury-thermometers), soil humidity (TDR probes) with a resolution of 10 min, and indoor air temperatures averaged over four rooms (Testo 174H, Testo AG, accuracy ±0.5 ∘ C for −20 to +70 ∘ C) with a 5 min resolution.Furthermore, short-wave transmissivity through the leaves was measured at two heights, each with two pyranometers, one behind and the other one in front of the greening (SP 110, Apogee).Mean transmissivity was 0.28.Off-site observations of global radiation and wind speed were carried out above roof level at the department of ecology (52.4572N, 13.3158E) with a distance of about 6 km to the study site.
The fac ¸ade greening was introduced into the models as a flat plant in front of the building with leaves from ground to top.Furthermore, we specified plant parameters (albedo, transmissivity, LAD, etc.) to represent the fac ¸ade greening as later described for the different models.
In order to analyze the influence of predefined meteorological data from observations and to compare it to the calculations by the models, we performed seven experiments (Table 2).The current version of ENVI-met 3 (E0a) does not support forcing of air temperature and relative humidity and thus additionally the unreleased version 4.0 Preview was applied (E0 unforced and E1 forced).The new version has already been used in other publications [33,34,39].SOLWEIG and RayMan require predefined values of air temperature and relative humidity.Both models can calculate short-wave downward radiation internally, which is the maximum short-wave downward radiation that can be derived at the specific date and location without clouds.These estimated values are used as global radiation to create experiments (R0 and S0) without measurement data of global radiation.
Building data for the input files was generated from the 3D city model of Berlin, which is available in the CityGML data format (http://www.citygml.org).We derived the spatial distribution and the height of vegetation from the Urban and Environmental Information System by the Senate Department for Urban Development and the Environment, Berlin.Terrain is neglected in all simulations because ENVImet 3 does not support it.Furthermore, we assume its effects to be negligible due to the flat terrain at the study site.Meteorological data for the simulation were derived from onsite observations averaged between the greened and the bare sites in front of the fac ¸ade for all variables with the exceptions of wind speed, wind direction, and global radiation.To initialize ENVI-met, daily mean values were applied.

ENVI-met. ENVI-met is a CFD model for simulating
surface-plant-air interactions from micro-to local scale with a focus on urban environments [18,38,40].It requires intensive computational resources of a personal computer as each of the experiments lasted over 10 days.ENVI-met is, in contrast to RayMan and SOLWEIG, able to simulate the physiological processes of vegetation and to describe vegetation not only as a porous obstacle.
The domain in ENVI-met had a spatial resolution of 1 m × 1 m and consisted of 110 × 80 grid points after removing 30 nesting grid points at each border.The domain size in ENVI-met was limited.Thus, the model domain is smaller than in SOLWEIG and RayMan.We used telescoping grids (factor 15%) starting at 2 m to include the highest building (43 m) because vertical grids were limited to 30 layers.In ENVI-met the wind field and turbulent heat flux were simulated as well.Therefore, the model domain needed to be

Integral radiation measurement
Indoor measurements rotated by 30 ∘ to avoid unrealistic roughness at the edges of the building walls (Figure 3).With the high resolution vegetation data, we introduced new vegetation classes based on the default plant database of ENVI-met to account for vegetation heights with a resolution of 1 m.The fac ¸ade greening was also introduced as a new class with the measured LAD values (Table 1).We interpolated the 9 measured levels to 10 relative height levels to adapt to the structure of ENVI-met's plant database (see Table 3).According to on-site experiences, soils beneath vegetation were classified as loamy sand, streets as asphalt, areas beneath buildings as pavement, and nonstreet areas as light brick roads.1.9 2.5 2.1 1.9 1.9 1.9 1.8 1.7 1.5 1.3 0.1

SOLWEIG. SOLWEIG simulates radiation fluxes and
mrt based on digital surface models [24,25].In contrast to ENVI-met and RayMan, SOLWEIG calculates shadow patterns and consequently  mrt for the middle of an hour; thus meteorological input was averaged accordingly.Vegetation and building data were also compiled with a resolution of 1 m.SOLWEIG uses trunk zones to calculate shadow patterns of vegetation, which were about 30% of the vegetation height at the study site.The fac ¸ade greening was represented as a flat plant without a trunk characterized by the observed values of albedo and transmissivity.In order to generate a SOLWEIG run (S0) without measurement data of global radiation, SOLWEIG1D was used.SOLWEIG1D calculates the maximum global radiation for the geographical location, which SOLWEIG2013a is not able to do.

RayMan.
Matzarakis et al. [22] developed the model RayMan (RayMan Pro 2.2) to calculate  mrt and biometeorological indices in complex environments.The building and vegetation data were created manually with the supplied editor on the basis of the digital surface models.Due to this approach, buildings had a flat roof in RayMan.In contrast to SOLWEIG and ENVI-met, RayMan calculates fluxes only for one point of interest.Thus, two obstacle files were used with different center points.Similar to SOLWEIG, RayMan can calculate global radiation itself.Consequently, we did experiments with calculated (R0) and measured global radiation (R1).The fac ¸ade greening was represented as series of plants in front of the building with a diameter of 1 m, a trunk length of 0 m, and an albedo of 0.28.

Analyses.
The analyses of the simulation results refer to two points in front of the fac ¸ade in each model domain, the greened site and the bare site (Figure 3).For the model evaluation (Section 3.2), time series plots are only presented for the greened site, but the statistical analysis includes both sites.
The models that did not calculate a specific climate element (e.g., air temperature) were excluded for the section.Long-wave radiation from RayMan and ENVI-met was calculated from surface temperatures using the Stefan-Boltzmann law (with an emissivity of 0.95).
The selected statistical measures for the model evaluation are a combination of generally recommended ones [41,42] and the ones mostly found in former studies (e.g., coefficient of determination ( 2 )).Root-mean-square deviation (RMSD), mean deviation (MD), and mean absolute deviation (MAD) were calculated as described by Schlünzen and Sokhi [42].For all statistical analyses we used IDL 8.2.2 (2007-2012 Exelis Visual Information Solutions, Inc.) and the library "Coyote" (1996-2014 Fanning Software Consulting, Inc.).

Measuring the Effects of Fac ¸ade
Greening on   . mrt decreased by 2.13 K (1000-1600 CET, 7, 23, 2014) in front of the greened fac ¸ade compared to the bare one (Figure 4).The components for calculating  mrt varied slightly between both sites and not only the short-wave radiation reflected and the long-wave radiation emitted from the fac ¸ade from NNE (Figure 5).Thus, the sums of short-wave and long-wave radiation downward, upward, and from the cardinal points differed between the bare and the greened site (Table 4).Nevertheless, long-wave radiation emitted from the greened site is clearly lower compared to the bare one.5).Despite high absolute errors in RMSD and MAD, experiments with ENVI-met produced the lowest MD (Table 5).

Long-and Short-Wave
Radiation.Short-wave downward radiation in front of the greened and the bare site was simulated with high  2 and low MD.RMSD showed larger deviation, especially in experiments with ENVI-met.Larger deviation also occurred in the morning and in the evening, which are the transition times between direct sun and shadow (Figure 7(a), Table 6).Simulated short-wave upward radiation differed much more from observations than short-wave downward radiation (Figure 7(b), Table 6).ENVI-met and SOLWEIG are able to compute short-wave upward radiation in contrast to RayMan.ENVI-met produced the amplitude more precisely than SOLWEIG, but SOLWEIG simulated the duration of intensive short-wave downward radiation closer to the observations (Figure 7(a)).
Concerning the long-wave downward radiation, SOL-WEIG (S0 and S1) was the closest to the observations (Figures 7(c) and 7(d)).Long-wave upward radiation was underestimated during the night in all models.During the day, the differences between the simulations and observations decreased (Figure 7(c)).

Air
Temperature and Specific Humidity.Simulations with RayMan and SOLWEIG are based on predefined air temperatures and specific humidity; only ENVI-met computes these variables itself.ENVI-met captured the diurnal cycle of air temperature, but specific humidity was reproduced with low deviations from observations only in the forced run (E1) (Figures 8(a) and 8(b)).E1 showed lower deviations than E0a and E0 in simulating air temperatures especially in the early morning (Figure 8(a), Table 7).E1 was the closest to the observation for specific humidity as well (Figure 8(b) and Table 7).The unforced experiments E0a and E0 overestimated specific humidity by 1.5 g/kg (RMSD).

Comparison of the Observed and Simulated Effects of
Fac ¸ade Greening.The observed differences between the bare and the greened site were not well reproduced by the models (Figure 9 and Table 8).The peaks in the transition times were simulated in all models to different extent.The cooling effect of fac ¸ade greening, however, was only simulated by S0 and S1, although too small.The experiments with ENVI-met showed a higher  mrt in front of the greened site than in front of the bare site, while R0 and R1 produced no differences between the sites (Table 8).Emitted long-wave radiation from

Measuring the Effects of Fac ¸ade
Greening on   .The observed differences in  mrt of 2 K were distinct but lower than the measurement accuracy of about 4 K. Tan et al. [15] detected larger differences of 2 to 13 K between sites with and without fac ¸ade greening in a tropical environment.They only detected a reduction of  mrt when the fac ¸ade greening was opposed to direct sunlight and not shaded [15].In addition, the reduction depended on the distance to the fac ¸ade and was mainly limited to 1 m.Thus, the small observed differences in this study can be partly explained with the larger distance of the measurement to the fac ¸ade of 1.2 m.Furthermore, the effects of fac ¸ade greening on air temperature were found to be small (1.5 K, Djedjig et al. [11], 1 K, Berry et al. [16], 0.25 K, Gross [14], and negligible, Pérez et al. [7]).In this case study fac ¸ade greening was only attached to one fac ¸ade wall.Hence, the impact on modified long-and short-wave radiation from the greened fac ¸ade accounts to only 22% in calculating  mrt (1).Therefore, the impact on  mrt was limited, even though emitted long-wave radiation was clearly reduced in front of the greened fac ¸ade.The albedo also interferes with  mrt because the light bare wall with its higher albedo reflected more short-wave radiation than the darker greened fac ¸ade.Furthermore, the differences in  mrt cannot be reduced to the effects of fac ¸ade greening exclusively.Boundary conditions were different between both sites regarding times of shadowing (Table 4).The opposing building structure is higher in front of the bare site than in front of the greened site (Figure 3), which leads to different shadowing times.In the evenings, at around 1700 CET, the greened fac ¸ade was completely shadowed while the bare site still received direct sunlight for another half an hour.In the morning, at around 0900 CET, when the sunlight appeared at the fac ¸ade, some differences existed as well.Furthermore, the opposing building might also have an influence regarding the reflection of short-wave radiation.With respect to short-wave and longwave radiation, the greened site received less energy from all directions than the bare one.Consequently, the cooling effect of  mrt was overestimated in the observations.This demonstrates shortcomings in the case study as both sites were not identical regarding the radiation characteristics, which introduces uncertainties when interpreting deviations between the sites.The problem of comparability and usability of case studies to observe the effects of fac ¸ade greening occurs often as criticized by Hunter et al. [10].Specifically in complex urban environments optimal study sites are difficult to find.Elaborated assessments (e.g., integral radiation measurements) can reveal such problems in contrast to simple single-variable measurements (e.g., air temperature).Our findings further support the need of reliable model tools in combination with observations to study the effects of fac ¸ade greening for specific sites and climates.

Mean Radiant
Temperature.Simulated  mrt deviated about 7 K (mean RMSD) from the observation in this and about MD 2-8 K and 2-15 K RMSD in former studies (Table 9).This uncertainty range is reasonable compared to the average accuracy in measurements of about 4 K.The uncertainty, however, varies largely depending on the selected model and study design (Table 9).Krüger et al. [20] concluded for RayMan that the uncertainty in calculating  mrt depends on different factors, such as morphology, meteorological conditions, and surface properties.For example, correlation coefficients were much higher in this study compared to the investigations of Chen et al. [17].Such deviations are caused by model specific parameterizations.
Other reasons for deviations in all applied models were discrepancies in building or plant shapes.All experiments produced larger deviations during the transition between shadow and direct radiation, which was also reported by Thorsson et al. [4] and Matzarakis et al. [23].Such deviations are caused by inaccuracies in building shapes due to the spatial resolution of 1 m or by errors in the basic data of building shapes.
Simulated  mrt differed hardly between the computed and predefined meteorological input data.Consequently, calculated global radiation by RayMan and SOLWEIG was sufficient for simulating  mrt , and observations of global radiations did not enhance the model performance considerably, at least for this nearly cloud-free day.In ENVI-met, the forcing of air temperature and relative humidity showed only minor effects on simulating  mrt as well.
In SOLWEIG, the transmissivity is of importance for calculating  mrt beneath vegetation [25].We applied an averaged value for the whole model domain because we measured transmissivity only at two points of the fac ¸ade greening and SOLWEIG does not allow any spatial variations.Nevertheless, the range of deviations in this study is similar  to findings by Lindberg and Grimmond [25], who detected a RMSD of 3.1 K and a MAD of 2.74 K.
RayMan underestimated  mrt in this study probably as a consequence of too low short-wave downward radiation.On Figure 10: Simulated long-wave radiation (a) emitted from the greened site (wall or leaves) and (b) differences between bare and greened fac ¸ade in simulations and observation for the experiments with ENVI-met 3 (E0a), unforced (E0) and forced ENVI-met 4 (E1), and SOLWEIG without (S0) and with global radiation from observation (S1) (7,23,2014).
12.88 K).Additionally, they detected a higher agreement with observations on clear sunny days compared to cloudy ones [20].During sunny days the accuracy of SVF has a high influence on the uncertainty in simulating  mrt [20].Hence, inaccuracies in the building and vegetation shape had a large impact in this study.The deviation in RayMan varies remarkably in this study and other ones [4,19,22] (Table 9).
The uncertainty also differs strongly between different studies with ENVI-met.Ali-Toudert [43] detected an underestimation of  mrt during the day and an overestimation during the night.In this study, however, the overestimation of short-wave radiation led to an overestimation of  mrt .Yang et al. [44] also detected an overestimation in ENVI-met due to an overestimation of short-wave downward radiation (without statistical values).
Overall, SOLWEIG simulated  mrt closest to the observation both in this study and in former ones (Table 9).ENVImet showed the largest deviations in this study.Furthermore, shortcomings in the evaluation of all three models become visible regarding the amount of studies and in comparability due to different statistics to evaluate models uncertainty.

Long-and Short-Wave
Radiation.Short-wave downward radiation showed less deviations compared to other radiation components.ENVI-met did not capture the amplitude of global radiation well, partly because only water vapor alternates short-wave radiation [44].Additionally, Ali-Toudert [43] detected problems in simulating global radiation with ENVI-met and consequently a radiation adjustment factor was introduced to ENVI-met.Here, short-wave downward radiation in ENVI-met was reduced in advance by 25% to match the amplitude of observations (Table 2).Additional deviations in short-wave downward radiation were produced by the resolution of grid points and by inaccuracies of the input data of buildings and vegetation as discussed for  mrt .The small deviation in short-wave downward radiation was a main driver for the sufficient simulation of  mrt .
Short-wave upward radiation has less impact on  mrt due to the generally smaller energy flux density and low   in the calculation of  mrt .Thus, the overestimation in all experiments was less decisive for simulating  mrt .ENVI-met overestimated the duration of intense shortwave upward radiation because shadowing effects of the surface are neglected.SOLWEIG captured the duration but overestimated the amplitude of short-wave upward radiation due to a higher albedo.SOLWEIG cannot consider different albedo values of ground surfaces and walls.In future, this problem might be fixed because Lindberg and Grimmond [25] proposed for the next SOLWEIG version the possibility to specify albedo and emissivity spatially.
Long-wave downward radiation was underestimated.The overestimations of short-wave upward radiation and longwave downward radiation compensate each other and thus the total impact of the inaccuracies on  mrt was reduced.ENVI-met 3 (E0a) failed to simulate the diurnal variation of long-wave downward radiation as it uses averaged emissivity values and surface temperatures to calculate it.
Long-wave upward is underestimated during the night because the models do not consider heat storage of the buildings sufficiently (ENVI-met 4) or not at all (ENVImet 3, SOLWEIG, and RayMan).Hwang et al. [19] and Lin et al. [21] discussed these shortcomings in ENVI-met 3 and detected the fact that daytime surface temperature was overestimated and nighttime surface temperature was underestimated.ENVI-met 4 should be able to consider these effects in principle but failed in this case study.Possibly, the heat storage was insufficiently parameterized or the spin-up time was too short.Follow-up studies might acquire a more suitable description of heat storage in the walls.During the night, the underestimation of emitted long-wave radiation also led to undervalued  mrt in all experiments.
Overall, the individual terms of the radiation balance were, with the exception of short-wave downward radiation, not plausibly simulated by all models.Short-wave downward radiation may also show much higher inaccuracies under cloudy conditions as Krüger et al. [20] detected for RayMan.Furthermore, the lack of model evaluation is more apparent for the individual terms of the radiation balance than for  mrt as we detected only three other studies (Table 10), which leads to a high uncertainty regarding the reliability in simulating these variables.

Air
Temperature and Specific Humidity.ENVI-met reproduced air temperature close to the observations with a RMSD between 0.9 and 1.6 K and MD between 0.0 and 0.4 K in this case study and with a mean RMSD of 1.9 K and MD −0.65 K in former studies (Table 11).During the day, ENVImet tended to underestimate air temperature but overestimated it during the night [21,43].Berkovic et al. [45], who compared simulations and observations qualitatively, found larger differences during the evenings for air temperatures and relative humidity (about 3 K, 15%).Forcing in E1 clearly decreased the uncertainty compared to E0 especially during the day but showed only minor effects on  mrt .
The uncertainty in simulating air temperature with ENVI-met varies between different studies remarkably (Table 11).A high agreement was achieved by Skelhorn et al. [35] and Srivanit and Hokao [26] after optimization and calibration efforts of different parameters and with a longer spin-up time.Air temperature has been evaluated in several studies.Thus, the simulation seems to be reliable, with an uncertainty of about 1.88 K (RMSD).
Specific humidity was evaluated by Chen et al. [17] with a RMSD of 0.82 g/kg, which is higher than in this study (RMSD = 0.35 g/kg).Forcing (E1) improved the simulation of specific humidity compared to the unforced runs (E0a, E0).

Comparison of the Observed and Simulated Effects of
Fac ¸ade Greening.The applied models were not able to reproduce the observed effects of fac ¸ade greening.SOLWEIG was the only model that simulated at least a small reduction of  mrt .The reductions cannot be related to transpirative effects of the plants because such effects are not considered in SOLWEIG.Thus, the decrease in  mrt was produced by radiation changes.Simulated long-wave radiation emitted from the wall was clearly reduced in front of the fac ¸ade greening (S0, S1).ENVI-met, which is able to simulate transpiration, reproduced higher  mrt in front of the greened site compared to the bare site.Emitted long-wave radiation was reduced in front of the greened fac ¸ade, but the effect was not permeated to  mrt .Increasing the spin-up time in ENVI-met might amplify the generation of effects of fac ¸ade greening.Specifically in ENVImet and RayMan, the inabilities in simulating a reduction of  mrt in front of the greened fac ¸ade could not be easily traced back because they do not return all relevant components for  mrt .
To conclude, the models showed a limited applicability to fac ¸ade greening in this study.The small alterations of longand short-wave radiation due to fac ¸ade greening were indistinct.Thus, specific types of urban green cannot be included Advances in Meteorology in these models without further adaptations.Restrictions must be made, however, regarding the complex real-case study site.Moreover, modification of model parameters and settings (e.g., spin-up time) especially for ENVI-met could lead to other results but would go beyond the scope of this study.

Conclusions
The effect of fac ¸ade greening on outdoor human bioclimate was limited in this case study because only a small reduction  mrt in front of the fac ¸ade greening was detected.Hence, fac ¸ade greening has only a minor effect in reducing outdoor heat stress.With a fac ¸ade greening attached to more than one fac ¸ade in a street canyon or court yard the effect on  mrt , however, might be enlarged.
The general ability of ENVI-met, RayMan, and SOL-WEIG to simulate  mrt was reasonable as expected for well-established models.Nevertheless, the deviations from observations vary largely between different studies.Additionally, the deviations from observations for other variables (specific humidity, long-wave downward or short-wave upward radiation) were higher and might impede the models' ability in assessing heat stress.When considering the large differences in complexity and computational time, the good performance of the simple SOLWEIG and RayMan models contrary to the elaborate ENVI-met model is encouraging.ENVI-met, however, offers more opportunities for various issues, such as studies of plant-air interactions or effects of changes in albedo of individual surfaces.It also provides more options for tuning and modifications by the users, which were not completely exhausted in this study.Moreover, we recognized a lack of model evaluations regarding the amount of evaluation studies and the considered variables.An explicit statement of model uncertainties for interpreting the results should be included in every study and not only in rare exceptional cases.
The applied models are helpful for assessing human bioclimate in general due to the acceptable uncertainty in simulating  mrt .In the specific case of fac ¸ade greening in a complex urban environment, however, their usability is limited in the current set-up.Generic studies or simple environments combined with modified parameterizations might improve the usability.Specific types of vegetation besides trees should not be incorporated in these models without modifications and extensive evaluation.Therefore, new simulation tools or advancements in existing models are desirable to complement observational case studies.The combination of biometeorological microclimate models and observations is helpful in order to complement benefits of each method.More effort in bridging the gaps between case studies and large-scale applications of countermeasures is needed to detect an effective countermeasure against heat stress risks in cities.

2 )Figure 1 :
Figure 1: Air temperature in front of the fac ¸ade and global radiation during the measurement campaign (7, 23, 2013).

Figure 2 :
Figure 2: Measurement arrangement and study site in Berlin, Germany.

Figure 3 :
Figure 3: Model domains of the different models.In the domain of SOLWEIG, dark colors represent taller buildings and vegetation compared to pale colors.Red circles indicate the greened and the bare site of observation and analyses.

Figure 4 :
Figure 4: Mean radiant temperature observed at the bare and the greened site in front of the building on 7, 23, 2013.

Figure 5 :
Figure 5: Short-wave radiation at the greened (a) and the bare (b) site and long-wave radiation at the greened (c) and the bare (d) site on 7, 23, 2013.All components used to calculate mean radiant temperature are represented, which are up-and downward radiation (up, down) and radiation from the cardinal points (west-northwest (WNW), east-southeast (ESE), south-southwest (SSW), and north-northwest (NNE) (direction from the fac ¸ades)).

Figure 7 :
Figure 7: Comparison of simulated and observed (a) short-wave downward radiation, (b) short-wave upward radiation, (c) long-wave downward radiation, and (d) long-wave upward radiation on 7, 23, 2013, for the greened site in front of the fac ¸ade simulated with ENVI-met 3 (E0a), unforced (E0) and forced ENVI-met 4 (E1), RayMan without (R0) and with global radiation from observation (R1), and SOLWEIG without (S0) and with global radiation from observation (S1).The gray areas indicate the accuracy range of the observations.

Figure 8 : 10 Δ
Figure 8: Comparison of simulated and observed (a) air temperatures and (b) specific humidity of the experiments with ENVI-met 3 (E0a) and unforced and forced ENVI-met 4 (E0 and E1) on 7, 23, 2013.The gray area indicates the accuracy range of the observations.

Table 1 :
Leaf area density (LAD) and depth of the fac ¸ade greening based on a vertical transect of (1 m × 9 m).

Table 2 :
Overview of the input parameters and model settings in the seven experiments."X" represents input of times series of meteorological data.

Table 3 :
Measured leaf area density (LAD) profile of fac ¸ade greening (bold) and assumed parameter for the plant database in ENVI-met.Plant Type Albedo Height (m) Root zone (m) Leaf area density (LAD) at level Root area density, all

Table 4 :
Sum (00-23 CET) of long-and short-wave radiation at the greened and the bare site up-and downward (up, down) and from the cardinal points (west-northwest (WNW), east-southeast (ESE), south-southwest (SSW), and north-northwest (NNE) (direction from the fac ¸ades)) on 7, 23, 2013.  is the angular factor used in calculation of mean radiant temperature (0.06 for up-and downward radiation and 0.22 for radiation from the cardinal points).

Table 6 :
Overview of the performance in simulating downward (a) and upward (b) short-wave and long-wave radiation on 7, 23, 2013.For further explanations refer to Table5.

Table 7 :
Overview of ENVI-met's performance in simulating air temperatures and specific humidity on 7, 23, 2013.For further explanations refer to Table5.

Table 9 :
Overview of deviations in mean radiant temperature between simulations and observations in other studies concerning root-mean-square deviation (RMSD), mean deviation (MD), mean absolute deviation (MAD), and coefficient of determination ( 2 ).Standard deviations (±) are stated, if more than one analysis is performed per study.

Table 10 :
Overview of deviations in short-wave downward radiation between simulation and observation in other studies.For further explanation refer to Table9.

Table 11 :
Overview of deviations in air temperature and specific humidity between simulation and observation in other studies.For further explanations refer to Table9.