Analysis of Urban Greening Scenarios for Improving Outdoor Thermal Comfort in Neighbourhoods of Lecce (Southern Italy)

: This study analyses the interactions and impacts between multiple factors i.e., urban greening, building layout, and meteorological conditions that characterise the urban microclimate and thermal comfort in the urban environment. The focus was on two neighbourhoods of Lecce city (southern Italy) characterised through ﬁeld campaigns and modelling simulations on a typical hot summer day. Field campaigns were performed to collect greening, building geometry, and microclimate data, which were employed in numerical simulations of several greening scenarios using the Computational Fluid Dynamics-based and microclimate model ENVI-met. Results show that, on a typical summer day, trees may lead to an average daily decrease of air temperature by up to 1.00 ◦ C and an improvement of thermal comfort in terms of Mean Radiant Temperature (MRT) by up to 5.53 ◦ C and Predicted Mean Vote (PMV) by up to 0.53. This decrease is more evident when the urban greening (in terms of green surfaces and trees) is increased by 1266 m 2 in the ﬁrst neighbourhood and 1988 m 2 in the second one, with respect to the current scenario, proving that shading effect mainly contributes to improving the urban microclimate during daytime. On the contrary, the trapping effect of heat, stored by the surfaces during the day and released during the evening, induces an increase of the spatially averaged MRT by up to 2 ◦ C during the evenings and a slight deterioration of thermal comfort, but only locally where the concentration of high LAD trees is higher. This study contributes to a better understanding of the ecosystem services provided by greening with regard to microclimate and thermal comfort within an urban environment for several hours of the day. It adds knowledge about the role of green areas in a Mediterranean city, an important hot spot of climate change, and thus it can be a guide for important urban regeneration plans.


Introduction
The growing impact of global warming is causing cities to encounter increasingly intense phenomena such as the Urban Heat Island (UHI). Among other things, UHI leads to negative consequences on human wellbeing, one of which is the thermal discomfort. Outdoor human thermal comfort in cities not only depends on air temperature but on the combined effect of air temperature, wind speed, air humidity, and radiation. The outdoor thermal comfort evaluation is usually performed combining micro-meteorological measurements or modelling and a guided user questionnaire survey on thermal perception and preference, as recently briefly summarized by Gatto et al. [1]. Recently, the improvement of thermal comfort has become one of the cornerstones of urban design and planning; in this perspective, urban greening, including forests, road trees, trees in parks, gardens, and trees in abandoned corners [2], plays an important role by providing a multitude of benefits to urban populations in reducing air, water, noise pollution, and the urban heat island, working as natural air conditioners, preventing solar radiation from heating the buildings and surfaces, cooling the air by evapotranspiration and reducing wind speed [3,4].
Urban greening alters the microclimate by means of several processes. The absorbed solar radiation energy induces an evaporation of water from vegetation to atmosphere (evapotranspiration processes) and then the atmosphere is cooled; in addition, the shade of vegetation protects other urban surfaces from solar radiation and prevents the increase of air and surface temperatures [3,5,6]. All these processes are usually modelled in Computational Fluid Dynamics (CFD) simulations [7]. The influence of urban greening on air quality can be positive or negative depending on different factors such as urban morphology, green features, location of emissions, meteorology, etc. (e.g., [8][9][10][11][12]). Therefore, careful planning of urban green spaces is absolutely essential, and the use of numerical modelling, allowing the simulation of different scenarios and the combination of several factors, can help to understand the right way to plan and manage urban greenery, so as to enhance its positive effects [13].
This study aims to investigate the effects of urban greening on air temperature and thermal comfort in urban and suburban neighbourhoods of Lecce city (Italy). The CFD microclimate model ENVI-met is employed because of its capacity to combine the simulation of processes induced by urban greening with atmospheric processes occurring in the upper atmospheric layer of built environments [14]. Different scenarios are evaluated by removing/adding trees in critical areas tailored to decrease high air temperatures characterizing the study areas, improve the thermal comfort, and provide some useful planting guidelines for urban planners. The thermal comfort is evaluated in terms of the mean radiant temperature (MRT) and the predicted mean vote (PMV)-two indices widely employed in research literature. A brief overview of indices employed for the evaluation of thermal comfort is given in Yang et al. [13].
The novelty of this study is that it compares the effects of tree species, characterised by different crown structure, height, and Leaf Area Density (LAD) in two different neighbourhoods of the same city. This provides a better understanding of the factors that influence the greening performance in terms of provision of ecosystem services and where improvements can be made as part of an urban regeneration plan aimed at mitigating climate change.

Materials and Methodology
This study was conducted through field measurements aimed at characterising the geometry of the urban environment, greening, and meteorology (see Section 2.2). In particular, the experimental measurement campaign (see Section 2.5), which led to the data collection of air temperature (T air ), was performed to validate the ENVI-met model simulations (see Section 2.4).

Study Area: An Overview on Lecce City
The neighbourhoods analysed here are located in Lecce (Figure 1), a medium-size city located on the central part of the Salentine Peninsula, in south-eastern Apulia Region. With about 94,000 inhabitants (http://demo.istat.it/popres/index.php?anno=2020&lingua= ita, accessed on 20 April 2021), Lecce is the main town of the second province in the region by population and one of the most important cities of Apulia. The geometry and morphology of the city are typical of the Mediterranean, with two-to-three-storey buildings and narrow street canyons. Lecce is built on flat land, at about 40-50 m above sea level and it is particularly interesting to study because it is part Mediterranean region, which has recently been described as a climate change "hot spot", where warming is particularly pronounced in summer and changes in precipitation patterns have been documented during the 20th century [15,16]. Based on the Köppen-Geiger classification, the most widely used quantitative classification of world climates, Lecce belongs to the Warm Mediterranean Climate 'Csa' class, with temperate-hot climate and specifically with average temperature of the coldest month between 6 • C and 10 • C, annual average between 14.5 • C and 17 • C, and 4 months with average > 20 • C. buildings and narrow street canyons. Lecce is built on flat land, at about 40-50 m above sea level and it is particularly interesting to study because it is part Mediterranean region, which has recently been described as a climate change "hot spot", where warming is particularly pronounced in summer and changes in precipitation patterns have been documented during the 20th century [15,16]. Based on the Köppen-Geiger classification, the most widely used quantitative classification of world climates, Lecce belongs to the Warm Mediterranean Climate 'Csa' class, with temperate-hot climate and specifically with average temperature of the coldest month between 6 °C and 10 °C, annual average between 14.5 °C and 17 °C, and 4 months with average >20 °C. In addition to the studies of climatic aspects on a larger scale, Lecce has been the subject of smaller-scale studies that have investigated the effect of urban morphology and greening on the urban microclimate in several neighbourhoods ( Figure 1). Maggiotto et al. [17] investigated the air temperature in four neighbourhoods using field measurements and modelling simulations. Observations showed that the city is characterized by higher air temperatures compared to surrounding areas, and air temperature has a large variability within the city and is very sensitive to urban morphology. Di Sabatino et al. [18] investigated the effects of Tilia cordata Mill. on local meteorology along a road not far from the city centre and specifically on wind flow and turbulence, by also evaluating building façades and ground temperatures with infrared images. The analysis showed that the presence of trees in street canyons may trap heat closer to the ground and this effect is accentuated during the night hours due to the trapping effect. Moreover, a reduced wind speed was observed. Martinelli and Matzarakis [19] investigated its courtyard thermal comfort in terms of Physiologically Equivalent Temperature (PET) using the RayMan model [20] and using micrometeorological data from a 30-year period 1980-2009. Results showed Lecce has the lowest PET excursion and steadier conditions throughout the year, presumably because of the stabilizing influence of the sea. Specifically, results showed 108 days with median PET < 13 °C and 85 days with median PET > 29 °C. Recently, Gatto et al. [1] evaluated the effect of urban greening on thermal comfort (in terms of Mean Radiant Temperature MRT and Predicted Mean Vote PMV) following a methodology based on field measurements and simulations with ENVI-met (the same as that followed in this In addition to the studies of climatic aspects on a larger scale, Lecce has been the subject of smaller-scale studies that have investigated the effect of urban morphology and greening on the urban microclimate in several neighbourhoods ( Figure 1). Maggiotto et al. [17] investigated the air temperature in four neighbourhoods using field measurements and modelling simulations. Observations showed that the city is characterized by higher air temperatures compared to surrounding areas, and air temperature has a large variability within the city and is very sensitive to urban morphology. Di Sabatino et al. [18] investigated the effects of Tilia cordata Mill. on local meteorology along a road not far from the city centre and specifically on wind flow and turbulence, by also evaluating building façades and ground temperatures with infrared images. The analysis showed that the presence of trees in street canyons may trap heat closer to the ground and this effect is accentuated during the night hours due to the trapping effect. Moreover, a reduced wind speed was observed. Martinelli and Matzarakis [19] investigated its courtyard thermal comfort in terms of Physiologically Equivalent Temperature (PET) using the RayMan model [20] and using micrometeorological data from a 30-year period 1980-2009. Results showed Lecce has the lowest PET excursion and steadier conditions throughout the year, presumably because of the stabilizing influence of the sea. Specifically, results showed 108 days with median PET < 13 • C and 85 days with median PET > 29 • C. Recently, Gatto et al. [1] evaluated the effect of urban greening on thermal comfort (in terms of Mean Radiant Temperature MRT and Predicted Mean Vote PMV) following a methodology based on field measurements and simulations with ENVI-met (the same as that followed in this study). Results showed that, on a typical summer day, the greening has a positive effect on thermal comfort by decreasing PMV of 2 points and MRT by 7 • C and that these variations are also related to urban morphology and greening type. Buccolieri et al. [21] investigated the effect of greening in the "Santa Rosa" neighbourhood (also studied here) on the dispersion of nitrogen oxides (NOx) and on CO 2 storage. A further neighbourhood, "167", is under investigation and will be the subject of a future publication.

Description of the Neighbourhoods
The first neighbourhood, referred to as "Santa Rosa", is located in the suburban northern part of the city; the second one, referred to as "Libertini", is located in the city centre ( Figure 1). Santa Rosa is a residential neighbourhood, sparser than Libertini, with several sensitive spots (a church, a primary school, and a small square with a playground) and university buildings. Libertini is a dense area that includes part of the Old Town, a limited traffic road usually overcrowded with pedestrians, and a traffic-busy street canyon (a) (Figure 2).  . Results showed that, on a typical summer day, the greening has a positive effect on thermal comfort by decreasing PMV of 2 points and MRT by 7 °C and that these variations are also related to urban morphology and greening type. Buccolieri et al. [21] investigated the effect of greening in the "Santa Rosa" neighbourhood (also studied here) on the dispersion of nitrogen oxides (NOx) and on CO2 storage. A further neighbourhood, "167", is under investigation and will be the subject of a future publication.

Description of the Neighbourhoods
The first neighbourhood, referred to as "Santa Rosa", is located in the suburban northern part of the city; the second one, referred to as "Libertini", is located in the city centre ( Figure 1). Santa Rosa is a residential neighbourhood, sparser than Libertini, with several sensitive spots (a church, a primary school, and a small square with a playground) and university buildings. Libertini is a dense area that includes part of the Old Town, a limited traffic road usually overcrowded with pedestrians, and a traffic-busy street canyon (a) (Figure 2). To characterize the building heights and the main greening features, field measurements were performed in both neighbourhoods. Specifically, buildings heights were measured using a laser rangefinder (LEICA DISTOTMpro4, Swiss Technology by Leica Geosystems, St. Gallen, Switzerland; www.leicadisto.it, accessed on 4 April 2021). The measurement operation took place in indirect mode. This procedure is used when the distance to be measured is difficult to reach; in this case, the DISTO calculates the desired length based on three auxiliary measurements according to the Pythagorean Theorem.
Concerning the greening features, the Leaf Area Index (LAI, m 2 m -2 ) of tree crowns was estimated from measurements of the photo-synthetically active radiation (PAR) acquired by an Accu-PAR LP80 ceptometer (https://www.metergroup.com/environment/products/accupar-lp-80-leaf-area-index/, accessed on 4 April 2021). All measurements were taken parallel to the ground. Five replicas were done just near the crown (where the sensor measured unobstructed PAR) and at its base (where it was supposed that the LAI is maximum). The leaf area density (LAD, m 2 m -3 ) was thus estimated by dividing LAI with the depth of the tree crown.
Regarding the Santa Rosa neighbourhood, Buccolieri et al. [21] reported part of the experimental campaign, which is summarised here and compared with the Libertini one.
The main characteristics of urban morphology including greening features of the two neighbourhoods here modelled are as summarized below: To characterize the building heights and the main greening features, field measurements were performed in both neighbourhoods. Specifically, buildings heights were measured using a laser rangefinder (LEICA DISTOTMpro4, Swiss Technology by Leica Geosystems, St. Gallen, Switzerland; www.leicadisto.it, accessed on 4 April 2021). The measurement operation took place in indirect mode. This procedure is used when the distance to be measured is difficult to reach; in this case, the DISTO calculates the desired length based on three auxiliary measurements according to the Pythagorean Theorem.
Concerning the greening features, the Leaf Area Index (LAI, m 2 m -2 ) of tree crowns was estimated from measurements of the photo-synthetically active radiation (PAR) acquired by an Accu-PAR LP80 ceptometer (https://www.metergroup.com/environment/products/ accupar-lp-80-leaf-area-index/, accessed on 4 April 2021). All measurements were taken parallel to the ground. Five replicas were done just near the crown (where the sensor measured unobstructed PAR) and at its base (where it was supposed that the LAI is maximum). The leaf area density (LAD, m 2 m -3 ) was thus estimated by dividing LAI with the depth of the tree crown.
Regarding the Santa Rosa neighbourhood, Buccolieri et al. [21] reported part of the experimental campaign, which is summarised here and compared with the Libertini one.
The main characteristics of urban morphology including greening features of the two neighbourhoods here modelled are as summarized below: • Santa Rosa has a dimension of 300 m × 300 m, with two street canyons (highlighted in red in Figure 2) with an aspect ratio H/W (where H is the mean height of the buildings and W the mean width of the street) equal to 0.6. The neighbourhood is characterized by buildings heights ranging from 8 m to 26 m with an average height of 13 m. The  Figure 2) characterised mainly by unirrigated turfs.

Scenarios Investigated
Three scenarios were investigated for each (

Model Description and Simulation Set-Up
The investigated scenarios were simulated by means of ENVI-met (http://www.envimet.info/doku.php?id=kb:review, accessed on 13 May 2021), a prognostic non-hydrostatic model able to analyse the surface-plant-air interactions. It solves Reynolds-averaged Navier-Stokes equations using a 1.5 order turbulence closure k-ε model and is composed of a 3D main model with a one dimensional (1D) atmospheric boundary layer (ABL) model, which permits to extend the analysis from the ground surface up to 2500 m. Evapotranspiration processes are simulated during the entire daily cycle on the basis of biological and thermo-physical properties of plants and soils.
The simulation set-up, summarised in Table 1, is the same as adopted and validated in Gatto et al. [1] and Buccolieri et al. [21]. This choice is aimed at making the studies comparable and to further contribute to in-depth guidelines for policy makers. Hourly air temperature T air and relative humidity were forced at the model boundary to drive the simulation

Model Description and Simulation Set-Up
The investigated scenarios were simulated by means of ENVI-met (http://www.envimet.info/doku.php?id=kb:review, accessed on 13 May 2021), a prognostic non-hydrostatic model able to analyse the surface-plant-air interactions. It solves Reynolds-averaged Navier-Stokes equations using a 1.5 order turbulence closure k-ε model and is composed of a 3D main model with a one dimensional (1D) atmospheric boundary layer (ABL) model, which permits to extend the analysis from the ground surface up to 2500 m. Evapotranspiration processes are simulated during the entire daily cycle on the basis of biological and thermo-physical properties of plants and soils.
The simulation set-up, summarised in Table 1, is the same as adopted and validated in Gatto et al. [1] and Buccolieri et al. [21]. This choice is aimed at making the studies comparable and to further contribute to in-depth guidelines for policy makers. Hourly air temperature Tair and relative humidity were forced at the model boundary to drive the simulation with meteorological input obtained from a 10 m high meteorological station (ARPA-Puglia, https://www.arpa.puglia.it/pagina2839_meteo.html, accessed on 5 May 2021) located about 1 km from Libertini and 3 km from Santa Rosa. The day 12 June 2018 was selected as representative of the 24-h period of the meteorological campaign (June 2018), since the temperature profile was reasonably correlated with the average temperature profile. Wind velocity and wind direction conformed to the other cases studied.

Micrometeorological Measurements
To evaluate the simulation performance in both the neighbourhoods investigated, modelling results were compared against T air measurements recorded continuously in June 2018 using Vaisala HMP45C thermo-hygrometers ( Figure 4) at Point 1 shown in Figure 5. The sensors were mounted on the banister of a balcony at the first floor of a building, 5 m above ground level to be representative of the urban canyon conditions. This was considered a reasonable compromise overcoming some of the inherent difficulties in positioning instruments in urban street canyons such as security from vandalism and limited access permission. The sensors were connected to a Campbell CR3000 data logger, which stored 10 min averaged values.

Micrometeorological Measurements
To evaluate the simulation performance in both the neighbourhoods investigated, modelling results were compared against Tair measurements recorded continuously in June 2018 using Vaisala HMP45C thermo-hygrometers ( Figure 4) at Point 1 shown in Figure 5. The sensors were mounted on the banister of a balcony at the first floor of a building, 5 m above ground level to be representative of the urban canyon conditions. This was considered a reasonable compromise overcoming some of the inherent difficulties in positioning instruments in urban street canyons such as security from vandalism and limited access permission. The sensors were connected to a Campbell CR3000 data logger, which stored 10 min averaged values.   Figure 4).

Calculated and Analyzed Parameters
Microclimate and thermal comfort analysis were carried out at z = 1.4 m, taken here as pedestrian level [22], extracting daily contours of Tair at strategically important points within the areas ( Figure 5) and calculating MRT, which can be defined as the radiation  Figure 4).

Calculated and Analyzed Parameters
Microclimate and thermal comfort analysis were carried out at z = 1.4 m, taken here as pedestrian level [22], extracting daily contours of T air at strategically important points within the areas ( Figure 5) and calculating MRT, which can be defined as the radiation heat absorbed by humans [23], and PMV, which takes into consideration both micrometeorological parameters and personal factors such as heat resistance of clothing and human activity. The choice of the indices investigated was also partly motivated by the need to compare the results with those already obtained in Gatto et al. [1].
MRT was calculated by simulating the short and long-wave radiant flux densities from the 3-dimensional environment and it was investigated here because it is considered to have a major impact on the estimation of many thermo-physiological indices describing human thermal comfort, including PMV, PET, Standard Effect Temperature (SET), and Universal Thermal Comfort Index (UTCI) [14]. PMV was simulated on the basis of meteorological parameters of T air , Rh, wind velocity (v), and MRT. As a recent literature review on ENVImet [24] confirms, PMV is now an index adopted worldwide and used in various studies because it considered comprehensively T air , RH, WS, MRT, and personal factors such as clothing heat resistance and human activity. Outdoor PMV ranges from −4 (very cold) to 4 (very hot), with 0 being neutral. For outdoor conditions, these values differ from the indoor theoretical limits (−3; +3) (ASHRAE Standard 55-2010 Thermal environmental conditions for human occupancy), and summer heat stress situations or extreme cold conditions can produce PMV values much higher than +4 (+8 and more) and much lower than −4 (−8 and less), as can be inferred from the PMV maps shown in Figure 8 of Gatto et al. [1]. Specifically, spatial maps for 3 h (12:00, 16:00 and 20:00), considered interesting by analysing the daily T air profile and also investigated in other studies [24], were extracted. The intention is, in fact, to investigate the spatial distribution of variables in different parts of the day and two warmer hours (where heat stress is more significant in summer and greater the impact of vegetation) and a less warm hour like sunset (where the impact of vegetation is less) were chosen. To better discriminate the species-specific thermal performance and compare it with data already available in the literature, the Heat Reduction Potential (HRP) index was calculated in specific areas ( Figure 2) at three selected hours. Specifically, the HRP equation was used in Morakinyo et al. [25] to evaluate the thermal performance of several tree forms in specific urban morphologies (street canyons). It was defined as the potential time-dependent human thermal comfort improvement in a vegetated area against the reference area and, therefore, represents the cooling potential of trees relative to the heat-stress. Morakinyo et al. [25] assessed thermal comfort using the PET thermal comfort index because its performance had already been evaluated in a subtropical hot-humid climate. In this study, the equation was modified using the PMV index, both to test its effectiveness and to compare the results with a previous study [1]. It is therefore estimated as: where: HRP PMV = Heat Reduction Potential in terms of thermal comfort (%) PMV SAt = PMV in the selected area with trees PMV SAwt = PMV in selected area without trees PMV AAWT = PMV in the whole area without trees

Model Validation
In Gatto et al. [1], the performance of ENVI-met with the chosen setup was validated. To further assess the model performance, several standard parameters were calculated for T air , namely the mean, the standard deviation (st.dev.), the higher value, the normalized mean square error NMSE, the fractional bias FB, and the correlation coefficient R. According to COST Action 732 [26], recommended criteria are: NMSE ≤ 1.5; −0.3 ≤ FB ≤ 0.3. Results of the statistical analysis are presented in Table 2 and show a quite satisfactory model performance. The comparison shows that ENVI-met was able to reproduce the T air behaviour (see Figure 5 of Gatto et al. [1]), as can be also seen from the statistical analysis (Table 2), with a lower overestimation (especially in Santa Rosa). The limitation of this model evaluation is that measurements in only one point are available in each neighbourhood. However, the recent review by Tsoka et al. [14] summarized the main studies performed with ENVI-met concluding that most studies have evaluated the model performance during summer rather than in winter, with T air the most widely evaluated variable with the existing evidence suggesting the model can generally accurately capture the diurnal profile of this parameter.  Figure 6 shows the temporal profiles of T air calculated at Points 1, 2, 3, and 4 in Santa Rosa and 1, 2, and 3 in Libertini (see Figure 5 for the position of the points) during the chosen summer day. In general, at each point, both in Santa Rosa and Libertini, T air is lower in the presence of trees with respect to the no-greening scenarios. The general profile trend of all the scenarios investigated is similar: a gradual increase of T air occurs from 6:00 to 13:00 local time, followed by a small plateau until 16:00, and then a decrease until 5:00 (the first hours are not shown in the graphs).

Temporal Profiles of Air Temperature
Main results are as summarized below: • Santa Rosa: Comparing the no greening scenario with the current one, the T air average daily decrease is about 0.55 • C at all points, with a maximum decrease of 0.92 • C recorded at 15:00 at P1. At P2, the maximum decrease is 0.60 • C at 16:00; at P3 it is 0.68 • C at 18:00; and at P4 it is 0.73 • C at 18.00. Looking at the modified scenario, the T air average daily decrease is more variable between the points; it is about 1.00 • C at P1 and P4, 0.68 • C at P2, and 0.55 • C at P3. On average, it is higher than the current one by about 0.25 • C. The maximum decrease of 1.10 • C is recorded at 15:00 at P1. At P2, the maximum decrease is 0.84 • C at 20:00; at P3 is 0.71 • C at 18:00; and at P4 it is 1.18 • C at 16:00. So, in general, it is possible to observe a maximum decrease in the presence of trees in the central hottest hours of the day; • Libertini: Comparing the no-greening scenario with the current one, the T air average daily decrease is about 0.25 • C in all points and it is more evident at P2 (street canyon) where the maximum decrease of 0.57 • C is recorded at 10:00. At P1 and P3, the maximum decrease is 0.15 • C at 16:00. As for Santa Rosa, looking at the modified scenario, the T air average daily decrease is more variable between the points; it is about 0.74 • C at P2, 0.61 • C at P1, and 0.34 • C at P3. On average, it is higher than the current one by about 0. Main results are as summarized below:

Spatial Distribution of Mean Radiant Temperature and Predictive Mean Vote
Here the analysis is extended to the investigation of the spatial distribution of MRT ( Figure 7) and PMV (Figure 8) at 12:00 (noon), 16:00 (mid-afternoon), and 20:00 (evening) to better understand how tree morphology and location, combined with different urban geometries, influence the thermal comfort during the day in the whole area.    In general, the trees cooling effect between different scenarios is greater in daytime than that at evening-time as already observed from the temporal variation of T air at different points. The difference in cooling effect between different scenarios decreases after sunset. The main results are summarized below: • on observing the spatial averages during the hot hours of the day (12:00, 16:00), the presence of trees induces a decrease of MRT in both Santa Rosa and Libertini. This decrease is most pronounced at 16:00, where it possible to observe a spatially averaged (average) decrease of 3.63 • C and 4.35 • C in the current and modified scenario in Santa Rosa, respectively, and 2.16 • C and 5.53 • C in Libertini. It is therefore also observable that there is a clear improvement in the modified scenario compared to the current one equal to 0.72 • C in Santa Rosa and 3.37 • C in Libertini; • on observing the spatial averages at 20:00, the decrease of MRT in the presence of trees is less marked, and nil in the current scenario in Libertini. Specifically, the spatial average is strongly influenced by what occurs in areas where the greening is most present. In fact, in the presence of trees, there is a localised increase in the street canyons in Santa Rosa, which is more evident in the street canyon "b" ( 1 • C) in both current and modified scenarios. In Libertini, it is possible to see that this increase in the street canyon "a" in the current scenario ( 0.60 • C) is less evident in the modified one, where there instead appears to be an increase in the investigated square and where new trees with high LAD were added ( 2 • C).
Looking at the spatial distributions of PMV ( Figure 8) at 12:00, 16:00, and 20:00, it is notable that they are distributed in the range of 3.62-5.38, 4.02-5.78 and 2.03-2.27 points, respectively. According to PMV value scale (Section 2.6), during daylight hours a condition of extreme heat stress was experienced with values greater than 3 points while, with values from 2 to 3 points, during and after sunset, a moderate heat stress was experienced. In line with what was observed in the MRT analysis, thermal comfort, expressed in terms of PMV, improves in the presence of trees in both Santa Rosa and Libertini in all 3 h investigated.
The main results are summarized as follows: • there does not appear to be much difference between 12:00 and 16:00, although in this case, even an improvement in the modified scenario compared to the current one is noteworthy. In both Santa Rosa and Libertini, the maximum decrease is at 16:00 and is 0.35 and 0.53 respectively. The PMV decrease is strongest in the street canyon "b" in Santa Rosa and it appears localized where the trees are in Libertini; • the improvement in thermal comfort at 20:00 is visible over the whole area but the rate of decrease of PMV is reduced compared to daylight hours. In Santa Rosa the decrease is almost unchanged between the current scenario and the modified one and is about 0.10, while in Libertini it is almost zero in the current scenario and equal to 0.10 in the modified one.

Discussion
To better understand and discuss the local effects of trees, the HRP index was calculated in the selected street canyons and squares ( Figure 2) and results are summarized in Table 3.

Day-Time Hours of a Typical Summer Day
The presence of trees induced a positive effect leading to a T air average daily decrease of about 0.40 • C and 1.00 • C in the current and modified scenario in Santa Rosa, respectively, and 0.57 • C and 0.74 • C in the current and modified scenario in Libertini. This improvement was also evident when looking at MRT and PMV, with a maximum decrease of MRT values equal to 4.35 • C and 5.53 • C and a maximum decrease of PMV values equal to 0.35 and 0.53 points at 16:00 in the modified scenario in Santa Rosa and Libertini, respectively. Such an improvement in the modified scenario is due to the increase of green areas replacing the paved area and to the increase of number of trees.
Beyond the species, the canopy density, expressed in terms of LAD, and the crown geometry, are very important parameters determining the cooling effect of trees, as already found in scientific literature [27][28][29]. The tree's shade reduces the incoming short-wave radiation and the mean radiant temperature, which shows a direct relationship with thermal comfort [30,31]. By combining the effect of tree species, geometry and arrangement, and urban geometry, it can be concluded that: • HRP index in the street canyons "a" of Santa Rosa and Libertini are characterised by similar aspect ratios and by the presence of a double lateral row of Quercus ilex L. with similar crown geometry (5 m high and wide crown) and LAD equal to 0.71 m 2 /m 3 in Santa Rosa and 0.40 m 2 /m 3 in Libertini. The greening induces a greater thermal comfort in Santa Rosa, with a maximum percentage of improvement at 16:00 equal to 15% in the current scenario and to 20% in the modified one, compared to 6% and 11% in the current and modified scenario, respectively, of Libertini. Libertini is denser and located in the city centre where the air recirculation is expected to be worse than at the edge of the city (where Santa Rosa is located); the presence of trees in already dense areas was shown to worsen the street canyon ventilation in some circumstances [7,18], and thus, in this case, the positive effect of trees in improving thermal comfort is counterbalanced by the heat trapping effect; shading effect in street canyons and squares: The effect is enhanced in denser geometries. For the same reason, in general, the effect of greening is found to be larger in street canyons rather than in the squares; further, in Santa Rosa Square, the thermal comfort improvement (in terms of both MRT and PMV) is smaller than in street canyons due to the presence of Pinus halepensis Mill., which has a smaller shade projection effect than that provided by Quercus sp. [32,33]; • ventilation: the wind direction is 210 • , which is a parallel to the street canyon "a" in Santa Rosa and oblique to the street canyon "a" in Libertini. This promotes better heat dispersion in Santa Rosa and thus explains the better HRP compared with Libertini neighbourhood; • evapotranspiration effect: T air decrease and thermal comfort improvement in the presence of trees is also due to the evapotranspiration effect, which helped cool the air, as already observed in other studies [5]. This decrease is consistent with results obtained in other Italian cities (Teramo, Rome, Palermo, Avola, Perugia), where, in the presence of greening, a decrease in the range between 0.5 • C and 3.0 • C was observed [34][35][36][37][38]; • replacing the asphalt surface with grass contributes to reduced MRT of about 0.50 • C.

Evening-Time Hours of a Typical Summer Day
During the evening, the decrease of the cooling effect by trees, visible from the HRP index and the spatial distribution of MRT and PMV, could be due to the fact that trees, after the stomata closure due to the lack of radiation, have almost no cooling effect from transpiration [28].
The local increase of MRT and PMV where trees are more concentrated and have a greater LAD (e.g., the street canyon "b" in Santa Rosa and the square in Libertini) could be due to the trapping of heat stored by the surfaces during the day and released during the evening. This increase could be positive especially during cold winters and, therefore, it should be additionally investigated in other urban geometries and different meteorological conditions. Similar results were found by Di Sabatino et al. [18] in other streets of Lecce. The notable increase in street canyon "b" can be explained by the additional street canyon effect which, due to the already reduced ventilation with the surrounding atmosphere, enhanced the trapping of heat below the crown; ventilation here is also reduced due to the wind blowing perpendicular.

Conclusions
This paper contributes to a better understanding of the ecosystem services related to microclimate and thermal comfort provided by urban greening (and specifically trees). This is an important target for developing climate change mitigation strategies in increasingly fragile and vulnerable contexts such as cities.
Specifically, this work deals with the study of the effects of trees on air temperature and thermal comfort in two neighbourhoods of Lecce (Italy), the first located in a suburban area and the second in the city centre.
Based on the findings of the study, it is suggested that policy makers and urban planners give preference to selecting trees with a large and dense crown in open areas (such as squares), because the cooling effect is determined more by the tree cover (in terms of LAD and number of specimens), dimension, and shape of tree crown rather than the tree species. Further, the lowering of LAD decreases the shading effect by tree crowns in the daytime, but enhances the nocturnal cooling of the space beneath the crowns. In the case of street canyons, it is always preferable to conduct micrometeorological analyses since ventilation is critical to heat dispersion (parallel wind leads to a better thermal comfort).
The introduction of greening in an urban regeneration plan may be strategic in achieving the Sustainable Development Goals (SDGs), foreseen in Agenda 2030, as suggested by the present study, but it must be based on ad-hoc evaluation because the effects of urban greening depend on multiple factors.