Abstract

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 are necessary to assess outdoor human bioclimate. We observed in front of a building façade in Berlin, Germany, which is half-greened while the other part is bare. was reduced (mean 2 K) 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 (MD = −1.93 K) and downward short-wave radiation (MD = 14.39 W/m2) were sufficiently simulated in contrast to upward short-wave and long-wave radiation. Finally, we compare the simulated reduction of 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.

1. 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 () 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].

Faç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 [79] or serving as habitat for wildlife [6]. The knowledge about the effectiveness of faç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 façade greening regarding surface and air temperature depend on solar irradiance, vegetation properties, and the particular greening system [1113]. Also, the mechanisms of cooling (e.g., through shadowing or transpiration) vary between different plant species [13]. The effects of faç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]. 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 façade greening on 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 façade greening on human bioclimate. Microclimate models, such as ENVI-met, RayMan, or SOLWEIG, are often used in different disciplines. Particularly over the last few years, the number of studies applying ENVI-met 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 façade greening. A comparison of 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 façade greening and we are not aware of other models that are able to simulate façade greening with regard to outdoor human bioclimate. ENVI-met, RayMan, and SOLWEIG include plants, but specific vegetation types, such as façade greening, are not explicitly declared to be supported. Nevertheless, all of the models should be able to consider some aspects of faç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 faç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 façade greening. ENVI-met moreover simulates evapotranspiration. Thus, several effects of façade greening are expected to be reproduced, such as changes in SVF, increased water-vapor fluxes, decreased long-wave radiation, or modified short-wave radiation.

In order to study the effects of façade greening on outdoor human bioclimate, we will firstly use observational data to answer the following question: (1) How large is the reduction of in front of a greened façade compared to a bare one at a study site? Afterwards, we will apply the models ENVI-met, RayMan, and SOLWEIG to the same site in order to (2) evaluate the general performance of the models in simulating 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 in front of the façade greening. In a wider sense the last question addresses if microclimate models are able to represent specific types of urban vegetation such as façade greening.

2. Material and Methods

2.1. Study Site

In order to observe and simulate the effects of façade greening on , we chose a south-southwest oriented building facade, which is half-covered with Parthenocissus tricuspidata while the other part is bare. The faç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 façade without technical climbing support. They rooted in a raised bed on the west oriented side of the building (unsealed area about 6 m2). 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 façade greening, and 0.18 for the ground surface. The mean depth of the vegetation layer is  m (Table 1). The average leaf area density (LAD) of the façade greening is 1.85 m2/m3 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.).

2.2. Observation

We measured on 7, 23, 2013 (00–23 CET), in front of the building faç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 ENVI-met, such as soil and indoor temperature (Figure 2).

The mobile measurements quantified 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 was calculated from mean radiant flux density () with the Stefan-Boltzmann law (2). The net radiometers were oriented along the faç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/m2), = long-wave radiation (W/m2), = 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/m2 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.

2.3. Models and Simulation

For the same day (7, 23, 2013, 00–23 CET), we applied the microclimate models: ENVI-met Version 3.1 Beta 5, ENVI-met Version 4.0 Preview [18, 38], SOLWEIG 2013a [24, 25], and RayManPro [22, 23]. The faç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 faç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 ENVI-met 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 on-site observations averaged between the greened and the bare sites in front of the façade for all variables with the exceptions of wind speed, wind direction, and global radiation. To initialize ENVI-met, daily mean values were applied.

2.3.1. 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 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 faç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.

2.3.2. SOLWEIG

SOLWEIG simulates radiation fluxes and based on digital surface models [24, 25]. In contrast to ENVI-met and RayMan, SOLWEIG calculates shadow patterns and consequently 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 faç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.

2.3.3. RayMan

Matzarakis et al. [22] developed the model RayMan (RayMan Pro 2.2) to calculate 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 faç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.

2.4. Analyses

The analyses of the simulation results refer to two points in front of the faç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 ()). 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.).

3. Results

3.1. Measuring the Effects of Façade Greening on

decreased by 2.13 K (1000–1600 CET, 7, 23, 2014) in front of the greened façade compared to the bare one (Figure 4). The components for calculating varied slightly between both sites and not only the short-wave radiation reflected and the long-wave radiation emitted from the faç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.

3.2. Evaluation of ENVI-met, RayMan, and SOLWEIG
3.2.1. Mean Radiant Temperature

The models simulated reasonably well in pattern and amplitude compared to the observations (Figures 6(a) and 6(b)). Experiments with SOLWEIG showed the lowest normalized standard deviation (Figure 6(b)) and high (Table 5). Despite high absolute errors in RMSD and MAD, experiments with ENVI-met produced the lowest MD (Table 5).

3.2.2. Long- and Short-Wave Radiation

Short-wave downward radiation in front of the greened and the bare site was simulated with high 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, SOLWEIG (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)).

3.2.3. 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).

3.3. Comparison of the Observed and Simulated Effects of Faç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 façade greening, however, was only simulated by S0 and S1, although too small. The experiments with ENVI-met showed a higher 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 the greened site was reduced in the experiments with ENVI-met and SOLWEIG (Figure 10(a)), but to a smaller extent (Figure 10(b)).

4. Discussion

4.1. Measuring the Effects of Façade Greening on

The observed differences in 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 façade greening in a tropical environment. They only detected a reduction of when the façade greening was opposed to direct sunlight and not shaded [15]. In addition, the reduction depended on the distance to the faç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 façade of 1.2 m. Furthermore, the effects of faç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 façade greening was only attached to one façade wall. Hence, the impact on modified long- and short-wave radiation from the greened façade accounts to only 22% in calculating (1). Therefore, the impact on was limited, even though emitted long-wave radiation was clearly reduced in front of the greened façade. The albedo also interferes with because the light bare wall with its higher albedo reflected more short-wave radiation than the darker greened façade.

Furthermore, the differences in cannot be reduced to the effects of faç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 faç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 faç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 long-wave radiation, the greened site received less energy from all directions than the bare one. Consequently, the cooling effect of 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 faç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 façade greening for specific sites and climates.

4.2. Evaluation of ENVI-met, RayMan, and SOLWEIG
4.2.1. Mean Radiant Temperature

Simulated 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 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 differed hardly between the computed and predefined meteorological input data. Consequently, calculated global radiation by RayMan and SOLWEIG was sufficient for simulating , 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 as well.

In SOLWEIG, the transmissivity is of importance for calculating beneath vegetation [25]. We applied an averaged value for the whole model domain because we measured transmissivity only at two points of the faç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 in this study probably as a consequence of too low short-wave downward radiation. On the contrary, Krüger et al. [20] examined an overestimation for various cloudiness conditions (RMSD 14.93 K, MAD 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 [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 during the day and an overestimation during the night. In this study, however, the overestimation of short-wave radiation led to an overestimation of . 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 closest to the observation both in this study and in former ones (Table 9). ENVI-met 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.

4.2.2. 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 . The small deviation in short-wave downward radiation was a main driver for the sufficient simulation of .

Short-wave upward radiation has less impact on due to the generally smaller energy flux density and low in the calculation of . Thus, the overestimation in all experiments was less decisive for simulating . ENVI-met overestimated the duration of intense short-wave 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 long-wave downward radiation compensate each other and thus the total impact of the inaccuracies on 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 (ENVI-met 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 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 as we detected only three other studies (Table 10), which leads to a high uncertainty regarding the reliability in simulating these variables.

4.2.3. 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, ENVI-met 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 .

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).

4.3. Comparison of the Observed and Simulated Effects of Façade Greening

The applied models were not able to reproduce the observed effects of façade greening. SOLWEIG was the only model that simulated at least a small reduction of . The reductions cannot be related to transpirative effects of the plants because such effects are not considered in SOLWEIG. Thus, the decrease in was produced by radiation changes. Simulated long-wave radiation emitted from the wall was clearly reduced in front of the façade greening (S0, S1).

ENVI-met, which is able to simulate transpiration, reproduced higher in front of the greened site compared to the bare site. Emitted long-wave radiation was reduced in front of the greened façade, but the effect was not permeated to . Increasing the spin-up time in ENVI-met might amplify the generation of effects of façade greening. Specifically in ENVI-met and RayMan, the inabilities in simulating a reduction of in front of the greened façade could not be easily traced back because they do not return all relevant components for .

To conclude, the models showed a limited applicability to façade greening in this study. The small alterations of long- and short-wave radiation due to façade greening were indistinct. Thus, specific types of urban green cannot be included 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.

5. Conclusions

The effect of façade greening on outdoor human bioclimate was limited in this case study because only a small reduction in front of the façade greening was detected. Hence, façade greening has only a minor effect in reducing outdoor heat stress. With a façade greening attached to more than one façade in a street canyon or court yard the effect on , however, might be enlarged.

The general ability of ENVI-met, RayMan, and SOLWEIG to simulate 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 . In the specific case of faç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.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors wish to thank their colleagues for the great support during the measurement campaign and the enriching discussions. The study is part of the Research Unit 1736 “Urban Climate and Heat Stress in Mid-Latitude Cities in View of Climate Change (UCaHS)” (http://www.UCaHS.org) funded by the Deutsche Forschungsgemeinschaft (DFG) under the codes SCHE 750/8-1, SCHE 750/9-1, and WE 1125/30-1.