Next Article in Journal
Psychometric Investigation of the Cultural Intelligence Scale Using the Rasch Measurement Model in South Korea
Next Article in Special Issue
Evaluation of Results of City Sustainable Transformation Projects in the Fields of Mobility and Energy Efficiency with Real Application in a District in Valladolid (Spain)
Previous Article in Journal
Evaluating the Implementation of the “Build-Back-Better” Concept for Critical Infrastructure Systems: Lessons from Saint-Martin’s Island Following Hurricane Irma
Previous Article in Special Issue
Dynamic and Heterogeneous Demand for Urban Green Space by Urban Residents: Evidence from the Cities in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sustainable Urban Greening and Cooling Strategies for Thermal Comfort at Pedestrian Level

1
Department of Civil Engineering and Architecture, University of Catania, 95125 Catania, Italy
2
Department of Electric, Electronics and Computer Engineering, University of Catania, 95125 Catania, Italy
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(6), 3138; https://doi.org/10.3390/su13063138
Submission received: 18 February 2021 / Revised: 9 March 2021 / Accepted: 10 March 2021 / Published: 12 March 2021

Abstract

:
The increase of the urban warming phenomenon all over the world is gaining increasing attention from scientists as well as planners and policymakers due to its adverse effects on energy consumption, health, wellbeing, and air pollution. The protection of urban areas from the outdoor warming phenomenon is one of the challenges that policy and governments have to tackle as soon as possible and in the best possible way. Among the urban heat island mitigation techniques, cool materials and urban greening are identified as the most effective solutions in reducing the urban warming phenomenon. The effects produced by the adoption of cool materials and urban forestation on the urban microclimate were investigated through a computational fluid-dynamic (CFD) model. The CFD model was calibrated and validated thanks to experimental surveys within the Catania University campus area. The urban microclimate thermal comfort analysis and assessment were carried out with the Klima–Michel Model (KMM) and Munich Energy Balance Model for Individuals (MEMI). In particular, three scenarios were performed: cool, low, and high levels of urban greening. The cool scenario, although it produces air temperature at around 1.00 °C, determines the worst condition of outdoor thermal comfort, especially at the pedestrian level. On the contrary, a high level of urban greening, obtained by the extensive green roofs together with an urban forestation, guarantees the wellbeing of pedestrians, showing more convenient values of PMV and PET.

Graphical Abstract

1. Introduction

The increase of urban population all over the world is leading to environmental, energy, and health implications of urban warming [1]. A recent study has shown that heatwaves, which are estimated to become more frequent and last longer in a hot climate, interact with urban warming to produce rather high heat stress for urban residents [2].
Furthermore, the features of the urbanized environments exacerbate the effects of global climate change [3]. Urban areas tend to have air temperatures higher than the rural surroundings as a result of gradual ground changes that include replacing the vegetation with buildings and roads [4].
The increase in temperatures in urban spaces is known as urban heat island (UHI) [5]. This phenomenon is due to several factors such as the specific composition of urban areas in terms of people and activities, morphology of neighborhoods, geometries, urban canyons [6], as well as surfaces with low albedo [7]. On a clear summer afternoon, the air temperature in a typical city could come up 2.5 °C higher than in the surrounding rural areas [8]. These factors affect the urban climate, energy uses, and livability of cities [9]. Urban areas are often characterized by higher vehicular traffic and less vegetation than their surrounding environment [10].
The radiative properties of the surface used in outdoor spaces have a key role, as they influence the air temperature of the lowest layers of the urban canopy layer [11]. Road pavements, claddings of facades, and roofs of the buildings are fundamental in the energy balance of the urban areas. The materials absorb solar radiation and store and release it by means of infrared radiation and convection heat transfer to their surrounding area during the evening and night-time [12]. Urban surfaces are characterized by values of solar absorption higher than rural surfaces [13]. Thereby, they have a superficial temperature lower than that of the building surfaces. Building materials emit an amount of energy higher than rural surfaces, and it leads to persistent warming in urban areas [14].
Increased urban temperatures produce the rising of the energy demand for building cooling, which in turn increases the greenhouse gas emissions and the outdoor thermal discomfort [15,16].
For evaluating outdoor thermal comfort in the urban contexts, several models are available in the scientific literature. Fanger’s indexes based on a steady-state energy balance model are implemented in many international standards [17,18]. Currently, the original equation of Fanger indexes is modified by Jendritzky and Nubler, 1981 [19,20], to evaluate the thermal comfort in external environments. They have produced the Klima–Michel Model (KMM) by adding the complex radiation balance in the outdoor environment and considering basic thermoregulatory processes through the human body, such as the constriction or dilation of peripheral blood vessels and the rate of physiological sweat [21].
The “Munich Energy-Balance Model for Individuals” (MEMI) [22,23] is the basis for the calculation of another index, called Physiological Equivalent Temperature (PET), which collected all the meteo-climatic factors in a single temperature index, namely the equivalent temperature. In KMM and MEMI models, the thermoregulatory processes are evaluated considering the human activity and climatic conditions.
Nowadays, mitigation strategies such as cool and green surfaces can enhance the outdoor thermal comfort in the urban areas, since both strategies reduce the sensible heat reaching the building surfaces, but the processes for reducing UHI are different [24]. Cool materials increase the reflectiveness of incoming solar radiation in urban areas reducing the net radiation [25]. Green surfaces increase the latent heat flux through the evapotranspiration process, while, thanks to the shading effect, reducing the incoming net radiation [26]. Cool materials applied on road pavements, with reflectivity 0.60, prove to be able to reduce the surface temperature by 5–6 °C compared to the conventional asphalt road [11]. In the Greek city of Florina, the application of cool pavements and cool roofs led to a reduction in the maximum air temperature by 1.39 °C and a decrease in soil temperature by 3.52 °C [27]. In Acharnes (Greece), reflective materials on the roofs allowed a reduction in surface temperature on the roofs between 1.51 °C to 10 °C and a decrease of 17% in energy consumptions [28]. The analysis of several mitigation strategies such as green roofs, green park, cool roofs, and cool pavements in an old densely populated neighborhood of Avola (Italy) highlighted that cool pavements reduce the environmental air temperature by over 2.0 °C [29]. Nevertheless, cool pavements can induce thermal discomfort to pedestrians and residents because of the high reflected solar radiation [30].
The adoption of vegetation in urban areas produces significant benefits by lowering air temperature and reducing air conditioning in the building [26,31]. It was found that medium trees and shrubs can cool down cities by from about 0.3 °C to 1.0 °C [4]. The results of CFD simulations on the effectiveness of several different mitigation strategies on the microclimate in Vienna showed that, at the pedestrian level, the planting trees allow achieving a decrease in air temperature of about 0.7 °C. A slight cooling in the air temperature, about 0.3 °C, was found when applying green roofs on all buildings of the investigated area. Moreover, the combined adoption of trees and green roofs can lead to a lower air temperature by around 1.0 °C [32].
However, few studies analyze critically the mitigation strategies that involve cool materials and urban forestation from the point of view of thermal comfort at pedestrian level. The purpose of this paper is to provide a comprehensive assessment of the effectiveness of three different scenarios based on the adoption of cool materials and low and high levels of urban greening.
These analyses have been developed through computational thermo-fluid-dynamic (CFD) simulations, calibrated with experimental measurements within the investigated area. Different parameters such as air temperature, wind speed, relative humidity, mean radiant temperature, predicted mean vote, and physiological equivalent temperature indexes were calculated and analyzed for each scenario. Finally, a comparative analysis was carried out to highlight the weaknesses and strengths of the different strategies.
The paper is organized as follows:
  • Section 2 describes the methodological approach, the tools and software used for CFD simulations, and analytic models for calculating comfort indicators. The equipment used for a ground survey is described too;
  • Section 3 provides detailed descriptions of the case study;
  • Section 4 discusses the methods used for calibrating and validating the model through statistical analysis approach;
  • Section 5 provides the results of computational simulations (CFD) for Cool and Green scenarios;
  • Section 6 reports scenarios’ results and some critical issues through comparative analysis.

2. Materials and Methods

2.1. CFD Simulations

Computational fluid-dynamic simulations have been carried out through the ENVI-met software [33], which is a three dimensional numerical prognostic calculation code based on a holistic model. This tool is able to recreate the microclimate and thermo-physical behavior of the urban microscale. Three-dimensional modeling of the urban area is based on the discretization of simulation domain in finite volumes. Thereby, the model is composed of many three-dimensional cells with a typical horizontal resolution from 0.5 to 10 m.
The software implements theoretical physical models of nonlinear equations, including three main prognostic variables, that allow the analysis of the interactions between environment and atmosphere on a local basis [34,35]. The main variables are average air flow, temperature and humidity, and turbulence. The analytic model is based on Navier–Stokes equations for wind flow, atmospheric flow turbulence equations, energy and momentum equation, and boundary condition parameters [34]. ENVI-met software implements the energy and mass balance equations related to:
  • Shortwave and longwave radiation fluxes attributable to shading, reflection, and re-radiation from building systems and the vegetation;
  • Transpiration, evaporation, and sensible and latent heat flux from the vegetation, including full simulation of all plant physical parameters [36];
  • Water and heat transfer inside the soil system.
The calculation procedure implies the following steps:
  • Geometry modeling of simulation domain;
  • Generation of the mesh of the calculation grid;
  • Implementation of buildings, obstacles geometry, trees, shrubs, and plants;
  • Definition of the boundary conditions and entering climate data;
  • Definition of the radiative and thermo-physical properties of buildings and urban materials;
  • Calculation of the field of temperature, humidity, wind speed, and comfort indicators;
  • Extraction of the simulation outcomes in 2D and 3D maps and in specific points at different height.
The results of the simulations are collected in several outputs folders that can be extracted as plots or text files.
In this study, outdoor air temperature (To), wind speed (w), relative humidity (RH), soil temperature (Tsoil), and mean radiant temperature (MRT) were considered useful parameters for comparing different scenarios. The evaluation of outdoor thermal comfort at pedestrian level, 1.00 m above ground level, in the different scenarios is carried out using the Klima–Michel Model (KMM) for calculating PMV index and Munich Energy-Balance Model for Individuals (MEMI) based on the Physiological Equivalent Temperature (PET).

2.2. Comfort Indexes

2.2.1. PMV

Predicted mean vote (PMV) is a well-known index defined by ISO 7730 [17]. It was initially intended for indoor use. In 1990, Jendritzky [20] introduced the Klima–Michel model (KMM) in the PMV in order to take into account the outdoor conditions. It is an adaptation of the Fanger model [17], with a short wave radiation model, for the calculation of the mean radiant temperature. This model assumes a typical activity and clothing of a man, 35 years old, 1.75 m tall, 75 kg weight, walking at 4 km/h (2.3 met).
Therefore, the values of the adapted PMV index [19,20] for the outdoor environment can overcome the international standard graduation used for indoor spaces. Indeed, while the PMV scale is defined between −3 (very cold) and +3 (very hot), where 0 is the thermal neutral (comfort) condition, the PMV for outdoor environment, may reach values above +4 or below −4 [37].

2.2.2. PET

The “MEMI” model is used for calculating the physiological equivalent temperature (PET). It is based on the energy balance of the human body (Equation (1)) and some parameters of the Gagge two-node model [38]. The most important advantage of the model is the possibility of calculating the physiological sweat rate, as a function of the average skin temperature (Tsk) and the core temperature (Tcr). The evaluation of both the heat flowing from the parts of the body surface covered by clothes or not is also allowed.
According to the Gagge model, the human body is divided into the outer region “skin” and internal body “core”. Therefore, the equation of the thermal balance between the human body and the outdoor environment is given by:
M + W + R + C + Ed + Ere + Esw = S = Scr + Ssk
where S—the heat stored through body heating or cooling; Scr—the heat stored in the core; Ssk—the heat stored in the skin; M—metabolic rate; W—Physical work output; R—Net radiative heat flux; C—convective heat flux; E—Latent heat flux exchanged through evaporation of moisture diffused through the skin (Ed), respiratory (Ere), and evaporation sweat (Esw).
The terms reported in Equation (1) are in (W/m2) and depend on the meteorological parameters [38], heat resistance of clothing, and human activity.
In particular, Esw depends on sweat frequency that is a function of the core temperature (Tcr), outdoor environmental conditions, and human activity. As the mean clothing temperature, mean skin temperature and sweat rate depend also on the climatic conditions; a substantial difference emerges with respect to Fanger’s theory.
The equation is solved under typical indoor settings (mean radiant temperature equal to air temperature MRT = Ta, water vapor pressure (pa) of 12 hPa, approximately equivalent to a relative humidity of 50% at Ta = 20 °C, and wind speed (w) of 0.1 m/s. In these conditions, the air temperature obtained solving the system is the PET.
PET is defined as the air temperature at which, in a typical indoor setting, the heat balance of the human body (work metabolism 80 W of light activity, added to basic metabolism; heat resistance of clothing 0.9 clo) is maintained with core and skin temperatures equal to those under the actual environmental conditions. The relationship between PMV and PET is shown in Figure 1.

2.3. Field Survey

The weather data implemented in ENVI-met are as follows: dry bulb temperature, relative humidity, global, direct and diffuse solar radiation, wind speed, and direction taken from a permanent weather station. The uncertainty, measurements range, and characteristics of the sensors are reported below:
  • Air temperature and relative humidity sensors: (uncertainty: 0.10 °C and 0.1%; measurements range: −50–100 °C and 0–100%);
  • Radiometer: (spectral response: 300–3000 nm; operative temperature −40 °C/+80 °C; uncertainty: ±4 W/m2);
  • Anemometer: (measurement range: 0 ÷ 50 m/s; threshold: 0.36 m/s; uncertainty 1%; below 3 m/s and 1.5% above 3 m/s; resolution: 0.06 m/s).
For the validation of the CFD simulations, the measurements were carried out with a portable microclimate station Heat Shield Master ELR615M by LSI Lastem.
Figure 2 shows all the components and sensors, such as a globe-thermometer, a thermo-hygrometric probe, and a rotor anemometer. Table 1 reports the technical characteristics of each probe.
The values of air temperature (To), globe temperature (Tg), wind speed (w), and relative humidity (RH) were recorded every minute. The instruments were tested and calibrated before the field survey, while the measurement process complied with the Standard ISO 7726 [40].

2.4. Statistical Analysis

Table 2 reports the statistical indexes for quantifying the discrepancies between the observed (mi) and simulated (si) data.
The mean absolute error (MAE) is the average of all absolute errors and expresses the difference between actual and predicted values. It represents the vertical or horizontal distance between each point and the equality line [41].
The root mean square error (RMSE) is instead the standard deviation of the residuals (i.e., the differences between observed and predicted values). Residuals are a measure of how far from the regression line data points are. It is worth noticing that few large errors can lead to a great RMSE value [42,43,44].
The Pearson coefficient (r) is given by the ratio between the covariance of the two variables and the product of their mean square deviations. It measures the intensity of the correlation between two random variables or two quantitative statistical characters.
The coefficient of determination (R2) is defined as the proportion of the variance in the dependent variable (si) that is predictable from the independent variable (mj). It indicates how simulated values are close to the regression line of the measured values. The index R2 ranges between 0 and 1; the upper value means that the simulated values match the measured ones perfectly [45].

3. The Case Study

The investigated area is a part of the University Campus of the metropolitan city of Catania along the Mediterranean coast in Southern Italy (latitude 37°30′ North and longitude 15°04′ East). This area was selected for its peculiar urban texture, high population density, as well as urban density.
According to the international Köppen–Geiger climate classification, Catania is characterized by warm and humid summer and moderately cool, wet winter. In the summer, the average outdoor temperature ranges from 23 °C to 35 °C, with peaks of 39 °C when hot winds blow from North Africa. Figure 3 shows the investigated area.
The area has a surface of around 33,600 m2 and includes five buildings 5.00 to 15.00 m high that host classrooms and Department offices. The squared buildings have huge sizes (35 × 35 m2), as shown in Figure 3b. Figure 4 shows the different typology of the land surface of the investigated area at current state.
Currently, about 8496 m2 are the building roofs and about 12,160 m2 is a paved area (pavements, roads), while the presence of vegetation is estimated at around 6610 m2 of the total investigated area.
The measurements were taken with a LSI-Lastem weather station (equipped with the air temperature sensor, black globe thermometric sensor, thermo-hygrometric sensor, and anemometer) located on the roof of the building of the Electric, Electronics, and Computer Engineering Department.

3.1. Baseline Scenario

3.1.1. Model Simulation

The simulation model presents a domain with a number of grids equal to 50 × 42 × 30 featured by a mesh of 4.0 × 4.0 × 1.0 m. The total domain dimension is 200.0 × 168.0 × 100.0 m. The height along the z-axis of the domain is constant until 15.0 m, which represents the size of the tallest building; after that, it increases with a telescopic factor rate of 10%. Along all borders of the domain, six cells of nesting grids are set. The time step of the simulation was constant and equal to two seconds to avoid numerical instability and convergence issues.
Figure 5 displays the 3D Model of the investigated area in the baseline scenario where green color represents the vegetated area (trees or shrubs), black color represents the paved area (road, parking, square), dark gray represents photovoltaic panels, light gray represents pedestrians’ roads, and light brown represents the unpaved area (ground, soil).
In the domain, three receptors, namely P1, P2, and P3, have been selected for the purpose of validating the CFD model (see Figure 6).
In the same way, two receptors (S1 and S2) have been selected to investigate locally the variation of microclimate and thermal comfort under the three proposed different scenarios (see Figure 5).
Two sequences of simulations were carried out: from 6.00 a.m. on July 28 to 6.00 a.m. on July 31 and from 6.00 a.m. on August 22 to 6.00 a.m. on August 25. The simulations started at sunrise when the energy balance is almost equal to zero.
Finally, the parameters for calculating PMV and PET were set up as follows: age of users, 24 years old, medium height of 1.70 m, average weight of 70 kg, metabolic rate equal to 1.50 Met, and clothing thermal insulation (clo) equal to 0.50.

3.1.2. Materials and Thermo-Physical Properties

Materials and thermo-physical properties have been modeled according to international literature sources [12]. The outer surfaces of buildings, roads, pedestrian pavements, and other outdoor spaces are made up of “standard” materials that absorb remarkable solar radiation during daytime contributing to the release of high heat flux.
Globally, eight different kinds of urban surfaces have been identified and defined. Aiming to simplify the model, for each of them, a single value of albedo and emittance based on the material predominance was assigned.
Table 3 reports the values of albedo and solar emittance adopted for the model.
It is worth highlighting that PV systems are installed on many roofs. Consequently, they have been modeled too. Indeed, they have been considered as conventional roofs with a reduction in reflectivity of around 0.625. The latter was calculated considering the conversion rate of the solar radiation into electricity (15%), the PV panel reflection (10%), and the portion of absorbed energy that comes back to the external environment by convection and radiation from the upper and lower surface of the PV modules (12.5%).
The vegetation occupies about 20% of the total land area in the model. Indeed, most trees and shrubs are placed to the north-west of the selected area, whereas the extensive vegetation and bare ground are facing south-east. The trees are predominantly conifers, but some deciduous plants are present, while flowerbeds are occupied by some small grass and evergreen plants. The plants identified on the area and simulated are conic: small trunk, sparse, medium 15 m (conifer); tree: 10.0 m height, dense foliage, leafless (conifer); and Populus Alba: 7.0 m height (deciduous) and base grass, average density, 0.50 m tall. All the plants are characterized by the albedo estimated in the range 0.25–0.40.
Thermo-physical and optical properties of plants as well as the leaf area density (LAD) are reported in Table 4.

3.2. On Site Weather Measurements

Figure 6 shows the actual locations where the portable weather stations have been placed during the field measurement and surveys. Table 5 reports the locations and time periods of the survey.
The portable station was settled at 1.2 m on the roof of the building n. 13 (P1 point) from 10.00 a.m. on July 29 to 6.00 p.m. on July 30 and in P3 point in close proximity to the stationary weather station from 12.00 a.m. on August 22 to 12.00 p.m. on August 24 (Figure 6a). The mobile station was installed on a tripod of 1.2 m height and placed close to the building n. 15 (P2 point) from 10.00 a.m. on July 30 to 6.00 p.m. on July 30.

4. Calibration and Validation

The simulation model was calibrated using hourly data of weather station reported in Figure 7.
The reliability of the ENVI-met model was evaluated against the statistical indicators introduced in Section 2.4 using the data recorded by the portable weather station in P1, P2, and P3. Table 6 reports the calculated statistical indicators.
The analysis of MAE and RMSE values for the air temperature and relative humidity reveal a not perfect correlation at the points P1 and P2, while at P3 the discrepancies are quite low. For the air temperatures, MAE at P3 is less of 1.0 °C, while anywhere else it exceeds 1.0 °C, and it is even above 2.0 °C at P2.
Nevertheless, a good correlation with experimental measurements is confirmed in terms of air temperature by the high values of the correlation coefficient “r“ and of the determination coefficient R2, which are always higher than 0.90. Although the statistical indicators do not achieve optimal values, the profiles of air temperature and relative humidity are very close and similar in all points.
Indeed, as an example, the comparison between the measured and simulated values at P1, displayed in Figure 8a,b, highlight a good correlation for air temperature (To) and relative humidity (RH).
Nevertheless, during the night, the temperatures simulated are around 1.0 °C lower compared with the measured ones. This could be due to an overestimation of the low long-wave radiation exchange with the sky because of low sky temperatures or high values of emittance of the outer layer of the building’s surface. Despite this, the model does not seem to be affected by the fraction of anthropogenic heat that actually contributes to climate microscale.
During the hours of maximum solar radiation, the profile of simulated air temperature is lower than about 1.0–2.0 °C with respect to that measured trend. A similar result was also found in other studies [45]. For example, Mohamed et al. [46] have evaluated the microclimate of an urban area located within the Islamic Quarter of Cairo (Egypt). They have compared ENVI-met and measured data, discovering a difference in peak value of the temperature of 2.0 °C between 11.00 a.m. and 3.00 p.m. on a summer day.
Gusson and Duarte performed on-site microclimate measurements for the calibration of ENVI-met models of two densely populated districts in São Paulo in Brasil (Brasilândia and Bela Vista). A grid with a mesh of 3 × 3 × 3 m3 was adopted for modeling both districts in ENVI-met. Despite good correlation coefficients, they found out a difference in peak value of temperatures of 3.2 °C between simulated and recorded data at around 12 p.m. [47].

5. Proposed Scenarios

Three scenarios were proposed and evaluated to enhance outdoor environmental microclimate conditions:
  • Cool scenario (Cs), in which a reflective paint is applied on the pedestrian pavements, roads, and building roofs;
  • Green Roof scenario (GRs), in which extensive green roofs are installed with an increase of surface cover at around 25%;
  • Green Saturation scenario (GSs), in which any available area is covered with vegetation and green roofs. In this way, the green surface cover increases by around 44%.
More in detail, the Cool scenario (Cs) consists in the replacement of conventional materials with cool paints for pedestrian roads, pavements roads, and the application of cool paint on the roofs of all buildings. The cool materials are supposed to have high albedo and high emissivity to reflect the solar radiation. Commercial products were adopted: their thermal and optical features are reported in Table 7 [48].
The total surface treated with cool materials is 20,650 m2, which represents 61% of the investigated area. The Green Roof scenario (GRs) is based on the installation of 8500 m2 extensive green roofs on the top of all buildings. The features of the substrate and vegetation layer of the green roof are reported in Table 8 and Table 9, respectively.
Data available in the scientific literature [49] were considered for the thermo-physical properties of substrate and vegetation layers. The Green Saturation scenario (GSs), in addition to the green roofs, considers the planting of new trees, such as Populus alba and pine, with an average height of 10 m, and new grass until it reaches the saturation of bare soil. Pines and Populus alba have been adopted, because they are endemic species and have physical characteristics that are suitable for urban greening. The surface of vegetation is 14,800 m2 in the green saturation scenario.

6. Results and Discussion

6.1. Microclimate and Thermal Comfort Daily Variation

This section reports the results of thermal comfort simulation for the Baseline, Cool, and Green scenarios. An accurate analysis has highlighted that Green Roofs (GRs) do not affect the microclimate at pedestrian level. For this reason, Figure 9a–d reports only the profile of air temperature, mean radiant temperature, predicted mean vote and physiological Equivalent Temperature for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios at the S1 and S2 receptors.
Under the Baseline scenario (Bs), the air temperature ranges from 22.8 °C to 31.0 °C for both S1 and S2 receptors. At noon, the Cool scenario (Cs), compared to the baseline scenario (Bs), gives place to a decrease of the air temperature by 0.5 °C at the S2 receptor and by 0.8 °C at the S1 receptor.
The Green Saturation scenario (GSs) shows the best performance at the receptors S2, which is close to green spaces, with a reduction in the air temperature of about 0.8 °C. Nevertheless, a decrease in the air temperature is also achieved at the receptor S1, which is quite far from green areas.
It is to outline that during the daytime the Cool scenario (Cs) shows values of MRT higher than in the green saturation (GSs) and baseline (Bs) scenarios; this happens for the zone closest to the building facades due to highest albedo of these surfaces. An increase of 4.0 °C is observed at S2 receptor and of 3.4 °C at S1 receptor.
As to the PMV, a value of 2.5 is exceeded for all scenarios for both receptors from 8:00 a.m. to 4.00 p.m. The only exception is observed at the receptor S2 for the Green Saturation scenario (GSs). A similar result is observed for the PET index, which is higher than 35.0 °C, “very strong heat stress”, excepting at the receptor S2 for the Green Saturation scenario (GSs).
However, during the hour of maximum solar radiation, the highest values of PMV and PET, according to the value of MRT, have been observed for the Cool scenario. It has to be highlighted that the trend of all investigated parameters show a rapid drop of the curve at around 2.00 p.m. on August 23 due to a decrease of direct solar radiation caused by clouds. In any case, the best results are achieved by the GSs, for which the PET ranges from 19.9 °C to 36.4 °C.
The cumulative frequency distribution was also calculated for the following parameters: To, MRT, PMV, and PET during the selected period from 12.00 p.m. on August 22 to 12.00 p.m. on August 24.
Actually, Figure 10a–d reports only the cumulative frequency distribution of air temperature and mean radiant temperature of the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios in the S1 and S2 receptors.
It can be seen that in the Cool scenario (Cs), the air temperature To exceeds 28 °C for 10% of the time. The situation is thus better than in the Baseline scenario (Bs), where To exceeds 28 °C for about 30% of the time. That may be explained on the basis of the previous remarks.
From Figure 10b, it is evident that the three scenarios are better differentiated in terms of MRT, especially for MRT higher than 33 °C. As to the parameters more strictly related to thermal comfort (PMV end PET), it can be seen from Figure 10c that all investigated scenarios reach the “Warm” level. Actually, the “warm level” is exceeded 25% of time in the (Cs) and 20% in the (Bs). The GSs exceeds for 10% of time the “Warm” category in the receptor S1. The effect of the Green Saturation scenario (GSs) is remarkable in receptor S2.
It worth highlighting that, as far as the PMV is concerned, the worst results were obtained in the (Cs), while the best performance was achieved by the GSs. This trend is even more evident in the cumulative frequency distribution of PET.

6.2. Spatial Analysis

6.2.1. Air Temperature, Wind Speed, and Relative Humidity

Figure 11 collects a horizontal view (x-y) of air temperature, wind speed, and relative humidity for the three scenarios at 1.00 m height above the ground level at 12:00 am on August 23.
In this section, the scenarios are analyzed according to a spatial point of view, considering the following parameters: air temperature (To), wind speed (w), and relative humidity (RH). In the Baseline scenario (Bs), the air temperature of 30.0 °C is reached in most of the areas analyzed.
The highest values of air temperature occur on waterproof surfaces, such as pedestrian pavement and roads, while the minimum value of temperature, of about 28.4 °C, occurs on unpaved surfaces. Under the Cool scenario (Cs), temperatures lower than 29.0 °C have been found in all the investigated areas.
Above soil surfaces, the temperatures range between 28.5 to 28.8 °C and over the roads and pedestrian pavements, the temperature is about 28.2 °C. Globally, cool pavements and roofs allow to reduce the peaks of about 1.0 °C.
For the Green Saturation scenario (GSs), the maximum external air temperature is 29.2 °C, whereas the minimum is 28.0 °C. Consequently, this latest scenario is the most effective, especially regarding the magnitude of minimum values of air temperature.
However, it is worth noticing that the cooling effect of vegetation is not limited to spaces where the plants were added, but their effects are spread to all the surrounding areas.
As regards the field of wind, it is possible to observe that the wind comes from the east, and its speed ranges from 0 to 2.25 m/s. As predicted, the highest wind speed is observed in open areas, whereas areas of stagnation appear behind the buildings or in the narrow spaces between the buildings. Moreover, in areas where vegetation is most present, substantial disturbances of wind speed occur as a consequence of the highest roughness of the soil surface.
The maps of relative humidity, for the three scenarios, indicate that the areas where vegetation is present are characterized by the highest relative humidity (up to 49%). In the Green Saturation scenario (GSs), the highest values of relative humidity were found. A maximum value of 51.6% was recorded in areas where the vegetation is very dense.

6.2.2. Soil Temperature, Mean Radiant Temperature, and Thermal Comfort Indexes

In this section, the scenarios are analyzed by a spatial point of view considering the following parameters: soil temperature (Tsoil), mean radiant temperature (MRT), predicted mean vote (PMV), and physiological equivalent temperature (PET) (Figure 12).
For the Baseline scenario (Bs), the soil temperature reaches the maximum values of 44.2 °C in paved areas, while just 27.0 °C is observed in the vegetated or shaded areas. The MRT has an average value of 52.0 °C in the whole area, which drops to 43.0 °C in shaded areas. The PMV index has a maximum value of 3.03 and a minimum value of 1.99.
The areas occupied by pedestrian pavements and roads have PMV values of about 3.00 and can be classified as hot areas, whereas the green areas with PMV values of 2.00 can be classified as warm areas. The PET index has a maximum value of 47.0 °C in paved areas and a minimum of 30.5 °C in the vegetation areas. Therefore, the Baseline scenario can be classified as hot and characterized by rather strong heat stress. Thus, for most of the investigated areas, the environment is very far from comfortable conditions, and high degrees of discomfort are perceived by the people.
The Cool scenario (Cs) shows values of soil temperature between 25.1 °C and 38.8 °C. On the surface of cool pavements, one can observe a reduction at around 5.0 °C compared to the maximum value of the Baseline scenario (Bs) and a reduction by 2.0 °C in the minimum value. It has to be highlighted that the mean radiant temperature reaches values of 55.0 °C almost in the whole area with an increment of about 3.0 °C compared to the Baseline scenario. This quite unexpected outcome depends on the remarkable amount of solar radiation reflected from the cool pavements and able to cause an increase in the temperature of the surrounding surfaces. The PMV ranges from a minimum value of 2.08 to a maximum of 3.50. In the Cool scenario, both the minimum and maximum values of PMV are higher than in the Baseline case. The PET index has a maximum value of 47.6 °C and a minimum value of 31.4 °C, with an increment of about 1.6 °C and 1.0 °C with respect to the Baseline scenario. It is worth to be highlighted that in the areas between neighboring buildings and in the spaces between the cool pavements and buildings is apparently a strong increase of PET. These areas suffer because of high reflected radiation and relative humidity.
Although the surface temperature of the pavements is reduced, the multiple heat reflections lead to an increase of the human thermal stress, as deductible by observing the magnitude of the maximum and minimum values of PMV and PET, compared to the baseline scenario.
Looking at the map of the reflected radiation, it is worth underlining that the cool pavements reflect in the average 200 W/m2 against 170 and 130 W/m2 of the Baseline (Bs) and Green Saturation (GSs) scenarios, respectively. The previous outcome confirms that the application of cool paint over a paved area contributes to the deterioration of the outdoor thermal comfort due to the increase of the reflected solar irradiance.
Figure 13 reports the 2D maps of reflected solar radiation (Gr) for all investigated scenarios on a horizontal plane at pedestrian level, at 12.00 a.m., on August 23.
The Green Saturation scenario (GSs) offers the best behavior in terms of MRT, PMV, and PET. As a matter of fact, MRT ranges between 34.4 °C and 52.3 °C. In addition, PMV has a minimum value of 1.4 and a maximum value of 3.0. Finally, PET assumes values between 28.0 °C and 41.1 °C.
To further characterize the various scenarios, it is useful to introduce the indicator IPET defined by the variation of PET between the Baseline (Bs) and the proposed scenarios (Ps,i), as follows:
I PET , i ( % ) = PET Ps , i PET Bs PET Bs
The i index is referred to the different proposed scenarios, i.e., Ps = Cs, Ps = GSs.
Figure 14 depicts the IPET for the Cool and Green Saturation scenarios at z = 1.00 m at noon on August 23. The positive sign means an increase of PET, while the negative sign indicated a reduction.
As it is possible to see, the Cool scenario shows an increase in the PET in all points of the investigated area. Figure 14 shows that the highest increase of the PET (10–15%) occurs in the internal courtyards of the buildings, mostly due to the multiple reflections. On the contrary, the Green Saturation scenario shows a reduction in the PET values in all investigated areas.
The highest reduction (−30% to −20%) is reached in the areas with trees and shrubs. It is worth pointing out that a decrease in PET of about 5% is achieved even in areas not provided with vegetation.
Of course, the beneficial effect of the Green scenario has to be attributed to the evapotranspiration of plants, which leads to a cooling effect in the surroundings. Moreover, the shading effect of trees and shrubs allows the creation of small areas where the warming sensation is attenuated.
Generally, under the Green Saturation scenario, it is possible to achieve acceptable outdoor thermal comfort with the exception of very limited areas represented by the surfaces between neighboring buildings and spaces around the buildings.

6.2.3. Spatial Cumulative Frequency Distribution

Figure 15 shows the cumulative frequency distribution of To, MRT, PMV, and PET for Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios with reference to the whole investigated area at height of 1.0 m above the ground level at 12.00 p.m. on August 23.
Overall, the graphs show that in spite of the fact that the Cool scenario (Cs) exhibits the best performance in terms of To, it proves to be the worst case in terms of MRT, because of the high reflectivity assumed by the surfaces treated with cooling paints.
The consequences can be seen in the PMV and PET, considering that both PMV and PET are more sensitive to MRT than to To.
This is a clear hint about the adverse effect that reflective paints may cause on the thermal sensation in the outdoor environment. In conclusion, compared to the baseline scenario, the best performance thus pertains to the GSs.

6.3. Spatial Averages

In order to provide a synthetic comparison of the three proposed scenarios, the average values of the thermal parameters and the comfort Indicators were calculated. The spatial average values on the whole area relating to air temperature, mean radiant temperature, relative humidity, predicted mean vote, and physiological equivalent temperature is as follows:
L ¯ = i = 1 n L n n
Here, L ¯ is the average spatial value of any previously mentioned (L) parameter extended to all cells (n) of the domain.
The average values of the investigated parameters for each scenario are compared and reported in Figure 16. In particular, Figure 16a displays a comparison in terms of outdoor air temperature, mean radiant temperature, and relative humidity. Figure 16b reports the comparison in terms of the mean values of PMV and PET according to the thermal sensation scale and heat stress categories. All the microclimate parameters and thermal comfort indicators refer to 1.00 m above the ground level and at 12.00 a.m. on August 23. Finally, PMV and PET normalized to their maximum values for the three analyzed scenarios are reported in Figure 16c.
Again, overall, Figure 16 shows that the best results are achieved by the Green Saturation scenario (GSs), which has the lowest of value To, MRT, PMV, and PET. These results are attributable to the synergic effect of shading and evapotranspiration processes. On the other hand, the Cool scenario (Cs) has the worst performance even with respect to the Baseline scenario (Bs), for all the mentioned parameters, with the exception of To.
It is to outline, that although the surfaces treated with green materials are less by 17% than those treated with cool materials, they allow achieving the greatest benefits for the urban microclimate and thermal comfort at the pedestrian level.

7. Conclusions

This paper aims to identify the effectiveness of urban greening and cool materials as strategies for the enhancement of urban micro-climate and well-being conditions of pedestrians. To this purpose, a micro-scale analysis, based on a numerical model and experimental field measurements, was carried out in the University Campus of the metropolitan city of Catania. The influence of materials and vegetation was analyzed to assess pedestrian comfort in different scenarios. The outcomes of the simulations reveal that the green roof does not affect the microclimate at the pedestrian level, while the application of cool materials can reduce the air temperature by almost 1.0 °C in the peak value compared to the other scenarios, while it is the Green Saturation scenario that ensured the best performance in terms of MRT, PMV, and PET.
Other meaningful results were also provided by the cumulative frequency curves. They allowed the estimation of how long a given condition of To, MRT, PMV, or PET is holding, and showed that the three scenarios are more sensitive to MRT than to To. Finally, the evaluation of PET and PMV provided comparable results. All the gathered information made it easy to assess the effectiveness of the various mitigation measures.
In conclusion, the methodology used here may be useful to urban planners and policymakers searching for the optimal urban greening and cooling strategies in the urban environment.

Author Contributions

Conceptualization: A.G. and F.N.; methodology: F.N., M.D., A.G., and L.M.; software: M.D. and F.N.; validation: M.D., A.G., and F.N.; formal analysis: L.M. and F.N.; investigation: M.D., A.G., and F.N.; resources: A.G., L.M., and F.N.; data curation: M.D., A.G., and F.N.; writing—original draft preparation: M.D., F.N., L.M., and A.G.; writing—review and editing: L.M., F.N., M.D., and A.G.; visualization: L.M. and A.G.; supervision: L.M., F.N., and A.G.; project administration: M.D., A.G., and F.N.; funding acquisition: A.G. and F.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Heatwave Plan for England–Making the Case: The Impact of Heat on Health–Now and in the Future. Department of Health, 1–21, 2015. Available online: https://www.england.nhs.uk/2015/06/2015-heatwave-plan (accessed on 11 March 2021).
  2. Li, D.; Boud-Zeid, E. Synergistic interactions between urban heat islands and heat waves: The impact in cities is larger than the sum of its parts. J. Appl. Meteorol. Climatol. 2013, 52, 2051–2064. [Google Scholar] [CrossRef] [Green Version]
  3. Climate Change Impacts and Vulnerabilty European. An Indicator Based Report. European Environment Agency (EEA): Copenhagen, 2012. Available online: https://energee-watch.eu/wp-content/up-loads/2014/05/Climate-change-impacts-and-vulnerability-in-Europe-2012.pdf (accessed on 11 March 2021).
  4. Akbari, H.; Levinson, R.; Miller, W.; Berdahl, W.P. Energy Saving Potentials and Air Quality Benefits of Urban Heat Island Mitigation, Lawrence Berkeley National Laboratory. 2005. Available online: http://escholarship.org/u-c/item/20j676c9 (accessed on 11 March 2021).
  5. Akbari, H.; Cartalis, C.; Kolokotsa, D.; Muscio, A.; Pisello, A.L.; Rossi, F.; Santamouris, M.; Synnefa, A.; Wong, N.; Zinzi, M. Local climate change and urban heat island mitigation techniques–the state of the art. J. Civ. Eng. Manag. 2016, 22, 1–16. [Google Scholar] [CrossRef] [Green Version]
  6. Alexandri, E.; Jones, O. Temperature decreases in an urban canyon due to green walls and green roofs in diverse climates. Build. Environ. 2008, 43, 480–493. [Google Scholar] [CrossRef]
  7. Akbari, H.; Konopacki, S. Calculating energy-saving potentials of heat-island reduction strategies. Energy Policy 2005, 33, 721–756. [Google Scholar] [CrossRef]
  8. Tzavali, A.; Paravantis, J.P.; Mihalakakou, G.; Fotiadi, A.; Stigka, E. Urban heat island intensity: A literature review. Fresenius Environ. Bull. 2015, 24, 4537–4554. [Google Scholar]
  9. Santamouris, M.; Synnefa, A.; Karlessi, T. Using advanced cool materials in the urban built environment to mitigate heat islands and improve thermal comfort conditions. Sol. Energy 2011, 85, 3085–3102. [Google Scholar] [CrossRef]
  10. Takebayashi, H.; Moriyama, M. Surface heat budget on green roof and high reflection roof for mitigation of urban heat island. Build. Environ. 2007, 42, 2971–2979. [Google Scholar] [CrossRef] [Green Version]
  11. Dimoudi, A.; Zoras, S.; Kantzioura, A.; Stogiannou, X.; Kosmopoulos, P.; Pallas, C. Use of cool materials and other bioclimatic interventions in outdoor places in order to mitigate the urban heat island in a medium size city in Greece. Sustain. Cities Soc. 2014, 13, 89–96. [Google Scholar] [CrossRef]
  12. Akbari, H.; Matthews, H.D. Global cooling updates: Reflective roofs and pavements. Energy Build. 2012, 55, 2–6. [Google Scholar] [CrossRef]
  13. Dandou, A.; Santamouris, M.; Synnefa, A.; Soulakellis, N.; Tombrou, M. On the Use of Cool Materials as a Heat Island Mitigation Strategy. J. Appl. Meteorol. Climatol. 2008, 47, 2846–2856. [Google Scholar]
  14. Oke, T.R. The energetic basis of the Urban Heat Island. Q. J. R. Meteorol. Soc. 1982, 108, 1–24. [Google Scholar] [CrossRef]
  15. Salvati, A.; Roura, H.C.; Cecere, C. Assessing the urban heat island and its energy impact on residential buildings in Mediterranean climate: Barcelona case study. Energy Build. 2017, 146, 38–54. [Google Scholar] [CrossRef] [Green Version]
  16. Mihalakakou, G.; Santamouris, M.; Papanikolaou, N.; Cartalis, C.; Tsangrassoulis, A. Simulation of the Urban Heat Island Phenomenon in Mediterranean Climates. Pure Appl. Geophys. 2004, 161, 429–451. [Google Scholar] [CrossRef]
  17. Fanger, P.O. Thermal Comfort. Analysis and Applications in Environmental Engineering; McGrawHill: New York, NY, USA, 1970. [Google Scholar]
  18. ISO 7730. Ergonomics of the Thermal Environment–Analytical Determination and Interpretation of Thermal Comfort Using Calculation of the PMV and PPD Indices and Local Thermal Comfort Criteria; International Standards Organization: Geneva, Switzerland, 2005. [Google Scholar]
  19. Honyo, T. Thermal Comfort in Outdoor Environment. Glob. Environ. Res. 2009, 13, 43–47. [Google Scholar]
  20. Jendritzky, G.; Nübler, W. A model analysing the urban thermal environment in physiologically significant terms. Meteorol. Atmos. Phys. 1981, 29, 313–326. [Google Scholar] [CrossRef]
  21. Höppe, P. Heat balance modelling. Experientia 1993, 49, 741–746. [Google Scholar] [CrossRef] [PubMed]
  22. Höppe, P. Die Energiebilanz des Menschen. Wiss. Mitt. Meteor. Inst. Univ. Miinchen 1984, 49, 1–12. [Google Scholar]
  23. Höppe, P. The physiological equivalent temperature—A universal index for the biometeorological assessment of the thermal environment. Int. J. Biometeorol. 1999, 43, 71–75. [Google Scholar] [CrossRef] [PubMed]
  24. Li, D.; Boud-Zeid, E.; Oppenheimer, M. The effectiveness of cool and green roof as urban heat mitigation strategies. Environ. Res. Lett. 2014, 9, 055002. [Google Scholar] [CrossRef]
  25. Santamouris, M. Cooling the cities—A review of reflective and green roof mitigation technologies to fight heat island and improve comfort in urban environments. Sol. Energy 2014, 103, 682–703. [Google Scholar] [CrossRef]
  26. Taha, H.; Akbari, H.; Rosenfeld, A. Heat island and oasis effects of vegetative canopies: Micro-meteorological field-measurements. Theor. Appl. Climatol. 1991, 44, 123–138. [Google Scholar] [CrossRef] [Green Version]
  27. Zoras, S.; Tsermentselis, A.; Kosmopoulos, P.; Dimoudi, A. Evaluation of the application of cool materials in urban spaces: A case study in the center of Florina. Sustain. Cities Soc. 2014, 13, 223–229. [Google Scholar] [CrossRef] [Green Version]
  28. Kolokotsa, D.D.; Giannariakis, G.; Gobakis, K.; Giannarakis, G.; Synnefa, A.; Santamouris, M. Cool roofs and cool pavements application in Acharnes, Greece. Sustain. Cities Soc. 2018, 37, 466–474. [Google Scholar] [CrossRef]
  29. Gagliano, A.; Nocera, F.; Detommaso, M.; Evola, G. Thermal behavior of an extensive green roof: Numerical simulations and experimental investigations. Int. J. Heat Technol. 2006, 34, S226–S234. [Google Scholar] [CrossRef]
  30. Santamouris, M. Urban warming and mitigation: Actual status, impacts and challenges. In Urban Climate Mitigation Techniques, 1st ed.; Santamouris, M., Kolokotsa, D.D., Eds.; Taylor & Francis Group: London, UK, 2016; pp. 1–26. [Google Scholar] [CrossRef]
  31. Susca, T.; Gaffin, S.R.; Dell’osso, G.R. Positive effects of vegetation: Urban heat island and green roofs. Environ. Pollut. 2011, 159, 2119–2226. [Google Scholar] [CrossRef] [PubMed]
  32. Maleki, A.; Mahdavi, A. Evaluation of urban heat islands mitigation strategies using Three dimensional urban micro-climate model ENVI-met. Asian J. Civ. Eng. 2016, 17, 357–371. [Google Scholar]
  33. ENVI-met V4.5. Urban Environment through Holistic Microclimate Modelling. Available online: https://www.envi-met.com/last (accessed on 13 February 2020).
  34. Bruse, M.; Fleer, H. Simulating surface–plant–air interactions inside urban environments with a three dimensional numerical model. Environ. Model. Softw. 1998, 13, 373–384. [Google Scholar] [CrossRef]
  35. Bruse, M. ENVI-met 3.0: Updated Model Overview. 2004. Available online: http://www.envi-met.net/documents/papers/overview30.pdf (accessed on 11 March 2021).
  36. Jacobs, C.M.J.; Van den Hurk, B.J.J.M.; de Bruin, H.A.R. Stomata behaviour and photosynthetic rate of unstressed grapevines in semi-arid conditions. Agric. For. Meteorol. 1996, 80, 111–134. [Google Scholar] [CrossRef]
  37. Fang, Z.; Feng, X.; Lin, Z. Investigation of PMV Model for Evaluation of the Outdoor Thermal Comfort. Procedia Eng. 2017, 205, 2457–2462. [Google Scholar] [CrossRef]
  38. Gagge, A.P.; Stolwijk, A.J.A.; Nishi, Y. An effective temperature scale based on a simple model of human physiological regulatory response. ASHRAE Trans. 1971, 77, 247–257. [Google Scholar]
  39. Giuffrida, S.; Gagliano, F.; Nocera, F.; Trovato, M.R. Landscape assessment and economic accounting in wind farm programming: Two cases in sicily. Land 2018, 7, 120. [Google Scholar] [CrossRef] [Green Version]
  40. ISO 7726. Ergonomics of the Thermal Environment–Instruments for Measuring Physical Quantities; International Organization for Standardization: Geneva, Switzerland, 1998. [Google Scholar]
  41. Despotovic, M.; Nedic, V.; Despotovic, D.; Cvetanovic, S. Review and statistical analysis of different global solar radiation sunshine models. Renew. Sustain. Energy Rev. 2015, 52, 1869–1880. [Google Scholar] [CrossRef]
  42. Ramos Ruiz, G.; Fernández Bandera, C. Validation of calibrated energy models: Common errors. Energies 2017, 10, 1587. [Google Scholar] [CrossRef] [Green Version]
  43. Trovato, M.R.; Nocera, F.; Giuffrida, S. Life-cycle assessment and monetary measurements for the carbon footprint reduction of public buildings. Sustainability 2020, 12, 3460. [Google Scholar] [CrossRef] [Green Version]
  44. Gagliano, A.; Patania, F.; Capizzi, A.; Nocera, F.; Galesi, A. A proposed methodology for estimating the performance of small wind turbines in urban areas. In Sustainability in Energy and Buildings, Smart Innovation, Systems and Technologies; M’Sirdi, N., Namaane, A., Howlett, R.J., Jain, L.C., Eds.; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar] [CrossRef]
  45. Salvati, A.; Kolokotroni, M. Microclimate Data for Building Energy Modelling: Study on ENVI-Met Forcing Data. In Proceedings of the 16th International Building Performance Simulation Association IBPSA Conference, Rome, Italy, 2–4 September 2019. [Google Scholar]
  46. Mohamed Elnabawi, H.; Hamza, N.; Dudek, S. Numerical modelling evaluation for the microclimate of an outdoor urban form in Cairo, Egypt. Hous. Build. Natl. Res. Center J. 2015, 11, 246–251. [Google Scholar] [CrossRef] [Green Version]
  47. Gusson, C.S.; Duarte, D.H.S. Effects of built density and urban morphology on urban microclimate-calibration of the model ENVI-met V4 for the subtropical Sao Paulo, Brazil. Procedia Eng. 2016, 169, 2–10. [Google Scholar] [CrossRef]
  48. CRRC Rated Products Directory. Available online: Coolroofs.org/products/search.php (accessed on 12 June 2020).
  49. Sailor, D.J. A green roof model for building energy simulation programs. Energy Build. 2008, 40, 1466–1478. [Google Scholar] [CrossRef]
Figure 1. The relationship values for PMV and PET model according to [39].
Figure 1. The relationship values for PMV and PET model according to [39].
Sustainability 13 03138 g001
Figure 2. Microclimate portable station: (a) LSI Lastem BVA304 Tripod/Mast Assembly 1.2 m; (b) Heat Shield Master ELR615M; (c) Globe thermometer ELR615M probe: Pt100; Thermo-Hygrometric sensors; (d) Rotor Anemometer DNA 202.
Figure 2. Microclimate portable station: (a) LSI Lastem BVA304 Tripod/Mast Assembly 1.2 m; (b) Heat Shield Master ELR615M; (c) Globe thermometer ELR615M probe: Pt100; Thermo-Hygrometric sensors; (d) Rotor Anemometer DNA 202.
Sustainability 13 03138 g002
Figure 3. University Campus in the neighborhood of S. Sofia: (a) 3D views; (b) 2D view.
Figure 3. University Campus in the neighborhood of S. Sofia: (a) 3D views; (b) 2D view.
Sustainability 13 03138 g003
Figure 4. Land use and surface cover (%).
Figure 4. Land use and surface cover (%).
Sustainability 13 03138 g004
Figure 5. Baseline scenario (BS): 3D view of the model.
Figure 5. Baseline scenario (BS): 3D view of the model.
Sustainability 13 03138 g005
Figure 6. Position of the receptors P1, P2, and P3: (a) real geometry; (b) ENVI-met model.
Figure 6. Position of the receptors P1, P2, and P3: (a) real geometry; (b) ENVI-met model.
Sustainability 13 03138 g006
Figure 7. Hourly values of outdoor air temperature (To), horizontal solar radiation (Hh), relative humidity (RH), wind speed (w): (a) 29–30 July; (b) 22–24 August.
Figure 7. Hourly values of outdoor air temperature (To), horizontal solar radiation (Hh), relative humidity (RH), wind speed (w): (a) 29–30 July; (b) 22–24 August.
Sustainability 13 03138 g007
Figure 8. Correlation between measured and simulated data at P1 point: (a) air temperature (To); (b) relative humidity (RH), on and July 29 and 30.
Figure 8. Correlation between measured and simulated data at P1 point: (a) air temperature (To); (b) relative humidity (RH), on and July 29 and 30.
Sustainability 13 03138 g008
Figure 9. Outdoor temperature (To) (a) and mean radiant temperature (MRT) (b) for the Baseline (Bs), Cool (Cs), and Green Saturation scenarios, Table 1 and S2. PMV (c) and PET (d) for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios at the receptor S1 and S2.
Figure 9. Outdoor temperature (To) (a) and mean radiant temperature (MRT) (b) for the Baseline (Bs), Cool (Cs), and Green Saturation scenarios, Table 1 and S2. PMV (c) and PET (d) for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios at the receptor S1 and S2.
Sustainability 13 03138 g009
Figure 10. Cumulative frequency distribution of outdoor temperature (To) (a) and mean radiant temperature (MRT) (b) for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios in the receptors S1 and S2. Cumulative frequency distribution of predicted mean vote (PMV) (c) and physiological equivalent temperature (d) for Table 1 and S2.
Figure 10. Cumulative frequency distribution of outdoor temperature (To) (a) and mean radiant temperature (MRT) (b) for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios in the receptors S1 and S2. Cumulative frequency distribution of predicted mean vote (PMV) (c) and physiological equivalent temperature (d) for Table 1 and S2.
Sustainability 13 03138 g010
Figure 11. Air temperature, wind speed, and relative humidity on horizontal plan (x-y) at z = 1.00 m and at 12.00 a.m. on August 23 for the Baseline (Bs), Cool (Cs), and Green Staturation (GSs) scenarios.
Figure 11. Air temperature, wind speed, and relative humidity on horizontal plan (x-y) at z = 1.00 m and at 12.00 a.m. on August 23 for the Baseline (Bs), Cool (Cs), and Green Staturation (GSs) scenarios.
Sustainability 13 03138 g011
Figure 12. Soil temperature, mean radiant temperature, predicted mean vote, and physiological equivalent temperature for Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios on a horizontal plan (x-y) at z = 1.00 m, at 12.00 a.m., on August 23.
Figure 12. Soil temperature, mean radiant temperature, predicted mean vote, and physiological equivalent temperature for Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios on a horizontal plan (x-y) at z = 1.00 m, at 12.00 a.m., on August 23.
Sustainability 13 03138 g012
Figure 13. Reflected solar radiation (Gr) on a horizontal plan (x-y) at z = 1.00 m, at 12.00 a.m., on August 23: (a) Baseline scenario (Bs); (b) Cool scenario (Cs); (c) Green Saturation scenario (GSs).
Figure 13. Reflected solar radiation (Gr) on a horizontal plan (x-y) at z = 1.00 m, at 12.00 a.m., on August 23: (a) Baseline scenario (Bs); (b) Cool scenario (Cs); (c) Green Saturation scenario (GSs).
Sustainability 13 03138 g013
Figure 14. Normalized percentage variation of PET (IPET) (a) Cool scenario (Cs); (b) Green Saturation scenario (GSs) on a horizontal plan (x-y) at z = 1.00 m, at 12.00 a.m., on August 23.
Figure 14. Normalized percentage variation of PET (IPET) (a) Cool scenario (Cs); (b) Green Saturation scenario (GSs) on a horizontal plan (x-y) at z = 1.00 m, at 12.00 a.m., on August 23.
Sustainability 13 03138 g014
Figure 15. Cumulative frequency distribution referred to the whole surface of investigated area at 12.00 p.m. on August 23 at height of 1.0 m above the ground level for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios: (a) outdoor temperature (To); (b) mean radiant temperature (MRT); (c) PMV; (d) PET.
Figure 15. Cumulative frequency distribution referred to the whole surface of investigated area at 12.00 p.m. on August 23 at height of 1.0 m above the ground level for the Baseline (Bs), Cool (Cs), and Green Saturation (GSs) scenarios: (a) outdoor temperature (To); (b) mean radiant temperature (MRT); (c) PMV; (d) PET.
Sustainability 13 03138 g015
Figure 16. The average values of the investigated parameters for the Baseline (Bs), Cool (Cs) and Green Saturation (GSs) Scheme 1.00 m, at 12.00 a.m., on August 23. (a) To, MRT, and RH; (b) PMV and PET with threshold values of thermal sensation and stress category; (c) PMV and PET normalized to their maximum values.
Figure 16. The average values of the investigated parameters for the Baseline (Bs), Cool (Cs) and Green Saturation (GSs) Scheme 1.00 m, at 12.00 a.m., on August 23. (a) To, MRT, and RH; (b) PMV and PET with threshold values of thermal sensation and stress category; (c) PMV and PET normalized to their maximum values.
Sustainability 13 03138 g016
Table 1. Technical features of probe.
Table 1. Technical features of probe.
Climatic VariablesType of ProbesAccuracyRange
Air temperatureThermometer
Silicon band gap
±0.3 °C, (@10, 40 °C)
±0.8 °C (@60 °C)
−20 °C ÷ 60 °C
Relative HumidityHygrometric Sensors1.8% (10 ÷ 90%)0 ÷ 100%
Wind speedDNA205 Rotor Anemometer2.5%0 ÷ 0.75 m/s
Globe temperatureELR615M globe thermometer probe (Ø = 150 mm), Sensor type: 1/3 DIN-A Pt100± 0.3 °C (@25 °C)−20 ÷ 125 °C
Table 2. Statistical indices used for validation.
Table 2. Statistical indices used for validation.
IndexNameFormulaU.M.
MAEMean absolute error MAE = i = 1 n | m i s i | n data-dependent
RMSERoot mean square error RMSE = i = 1 n ( m i s i ) 2 n data-dependent
rPearson correlation coefficient r = i = 1 n ( m i m ¯ ) ( s i s ¯ ) i = 1 n ( m i m ¯ ) 2 · i = 1 n ( s i s ¯ ) 2 -
R2Coefficient of determination R 2 = 1 i = 1 n ( m i s i ) 2 i = 1 n ( m i m ¯ ) 2 -
Table 3. Thermal and optical properties of existing material.
Table 3. Thermal and optical properties of existing material.
Urban ComponentsAlbedo (ρ)Solar Emittance (ε)
Roads, parking, square (asphalt)0.200.85
Pedestrian roads0.300.85
Ground, soil0.200.95
Concrete building roofs0.300.90
Clay brick building roofs0.300.90
Building walls0.200.90
Photovoltaic panels0.6250.90
Shelter parking0.200.90
Table 4. Thermo-physical and optical properties of plants included in the model.
Table 4. Thermo-physical and optical properties of plants included in the model.
TypeSpeciesHeight (m)LAD (-)ρ (-)ε (-)
Conic, small trunk, denseconifer15.02.30.120.10
Tree dense foliage leaflessconifer10.02.50.250.12
Populus albadeciduous7.02.00.400.20
Base grass, average density-0.300.300.250.15
Table 5. Locations and investigated periods selected for measurements survey.
Table 5. Locations and investigated periods selected for measurements survey.
NoLocationsTime PeriodsLocationsHeight a.g.l.
(m)
1P110.00 a.m. on July 29–6.00 p.m. on July 30Roof Terrace of building. n. 135.4
2P210.00 a.m.–6.00 pm. on July 30Pavements7.0
3P312.00 a.m. on August 22–12.00 a.m. on August 24Roof Terrace of building. n. 135.4
Table 6. Validation indices of the ENVI-met model.
Table 6. Validation indices of the ENVI-met model.
VariablesStatistical IndicesUnit Locations
P1 (h = 5.41 m)
29th–30th July
P2 (h = 1.20 m)
30th July
P3 (h = 1.20 m)
22th–24th August
ToMAE°C1.332.410.82
RMSE°C1.602.500.69
r-0.960.950.98
R2-0.930.910.94
RHMAE%4.764.953.04
RMSE%5.496.022.15
r-0.940.770.97
R2-0.880.600.95
Table 7. Surface cover (%) treated with cool materials.
Table 7. Surface cover (%) treated with cool materials.
Urban ComponentsSurface (m2)Surfaces Ratio (%) Albedo (ρ) Emittance (ε)
Cool pavements12,16036% 0.800.90
Cool roofs850025%0.800.90
Total20,65061%--
Table 8. Features of the substrate layer.
Table 8. Features of the substrate layer.
HeightH0.15m
Thermal conductivityλ1.00W·m−1·K−1
Absorptanceα0.60-
Albedoρ0.30-
Emittanceε0.90-
Moisture contentΘsat0.50m−3·m−3
Table 9. Features of the vegetation layer.
Table 9. Features of the vegetation layer.
HeightH0.30m
Leaf Area IndexLAI1.50m2·m2
Absorptanceα0.60-
Albedoρ0.25-
Emittanceε0.95-
Transmissivityτ0.15-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Detommaso, M.; Gagliano, A.; Marletta, L.; Nocera, F. Sustainable Urban Greening and Cooling Strategies for Thermal Comfort at Pedestrian Level. Sustainability 2021, 13, 3138. https://doi.org/10.3390/su13063138

AMA Style

Detommaso M, Gagliano A, Marletta L, Nocera F. Sustainable Urban Greening and Cooling Strategies for Thermal Comfort at Pedestrian Level. Sustainability. 2021; 13(6):3138. https://doi.org/10.3390/su13063138

Chicago/Turabian Style

Detommaso, Maurizio, Antonio Gagliano, Luigi Marletta, and Francesco Nocera. 2021. "Sustainable Urban Greening and Cooling Strategies for Thermal Comfort at Pedestrian Level" Sustainability 13, no. 6: 3138. https://doi.org/10.3390/su13063138

APA Style

Detommaso, M., Gagliano, A., Marletta, L., & Nocera, F. (2021). Sustainable Urban Greening and Cooling Strategies for Thermal Comfort at Pedestrian Level. Sustainability, 13(6), 3138. https://doi.org/10.3390/su13063138

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop