The Impact of Tipuana tipu Species on Local Human Thermal Comfort Thresholds in Different Urban Canyon Cases in Mediterranean Climates : Lisbon , Portugal

Based upon the case of Lisbon, this article examined the in-situ effects of vegetation upon pedestrian thermal comfort levels. Focussing specifically upon the historic quarter that often witnesses the highest Tamb values and Urban Heat Island (UHI) intensities during the summer, the most common urban canyon cases (UCCs) were modelled, along with one of the most commonly used vegetative semi-deciduous species found in the city, Tipuana tipu. Based upon a reference point (RP) system, the assessments were undertaken through the use of a new version of the SkyHelios model, local obtained Grad values, and the modified physiologically equivalent temperature (mPET) index calculated through the human-biometeorological model RayMan. The study identified the in-situ thermo-physiological influences of Tipuana tipu during different periods of the year: (1) during the summer, which revealed considerable reductions of PET/mPET of up to 15.6 ◦C/11.6 ◦C during a very hot day (where daily maximum Tamb surpassed 35 ◦C); and (2) during the winter, which revealed the risks of oversharing as a result of the species keeping its foliage during the winter with reductions of PET/mPET of up to 2.7 ◦C/2.6 ◦C. Furthermore, the study utilised the climate tourism/transfer information scheme (CTIS) to categorise and facilitate the interpretation of the results.


Introduction
Currently, and within existing cities, particularly those concomitant to Mediterranean climates with hot-dry summers and classified with a Köppen Geiger (KG) of 'Csa', the importance of local human thermo-physiological thresholds is gaining new weight for local urban design and planning.Still, and often within southern Europe, many cities habitually present a noteworthy absence of climatological data/know-how that could prove useful for local decision making and design within the urban public realm [1].
As a result, and within the context of Lisbon, there have been a substantial amount of studies that have provided valuable contributions towards the comprehension of the overall bioclimatic conditions within its public realm [2,3], UHI intensities [4,5], and integration with planning policy [6,7].Moreover, and given the progress of the climate change adaptation agenda, studies specifically pertaining to impacts within Lisbon have also been undertaken [8][9][10].
Within the international arena, including within other cities located in other KG classifications, the influences of urban vegetation have already been debated within existing review studies [11][12][13][14].Within such studies, the benefits of vegetation upon the urban public realm were reviewed, including their capacity to regulate urban ambient temperature (T amb ), and UHI effects.In addition, the scientific community has also endeavoured to expand upon the known in-situ effects of vegetation, including their respective thermo-physiological influences upon pedestrians.Within this line of study, the particular influences of tree species have also been identified within different urban conditions [15][16][17][18][19][20].It is also worth noting that, within studies which identified bioclimatic conditions/opportunities within different default urban canyon cases (UCCs) (both symmetrical and asymmetrical), the thermo-physiological benefits from vegetation were also mentioned, namely by (i) Ali-Toudert and Mayer [21] in Ghardaia, Algeria (with KG of 'Bwh'); (ii) Ketterer and Matzarakis [22] in Stuttgart, Germany (with KG of 'Cfa'); Qaid and Ossen [23] in Putrajaya, Malaysia (with KG of 'Af '); (iii) Morakinyo and Lam [24] in Hong Kong (with KG of 'Cwa'); (iv) Algeciras, Consuegra [25] and Algeciras, Tablada [26] in Camagüey, Cuba (with KG of 'Aw'); and lastly, (v) Nouri, Costa [27] in Lisbon, Portugal (with KG of 'Csa').
Notwithstanding, and as suggested by Kong, Lau, et al. [28] within the scientific community, there are still a greater amount of studies which have focussed upon T amb variations rather than thermo-physiological impacts as a result of local vegetation.Concomitant with this perspective, this article is centred upon the importance of going beyond sole climatic characteristics, and aims to add to the existing studies which identify the significance of considering 'in-situ' thermo-physiological influences resultant of urban vegetation.
In line with this objective, this study (i) utilises a new version of the SkyHelios model [29,30] as a new means to address microclimatic characteristics (namely Wind (V) speed); and, (ii) examines and applies new modified thermo-physiological indices [31] to conduct more accurate evaluations of thermal comfort conditions during different periods of the year.Such methods are applied to the case of Lisbon to identify how one of the most commonly-found shading trees within the city (i.e., Tipuana tipu) can influence the in-situ bioclimate conditions within symmetrical urban UCCs which are typically located within the city's historical district.The results of the study indicate (1) how simulations can be united with field measurements to compare and adapt climatic data from the local meteorological station; (2) the in-situ thermo-physiological effects that a common shading tree can have upon pedestrians during the summer and winter periods in areas prone to higher thermal stimuli, particularly during the summer; and lastly, (3) how such outcomes can ease the transition and 'shared language' [1] between the fields of urban climatology and local urban design/planning.

Site
Located on the western coast of Portugal at 38 • 42 N and 9 • 8 W, Lisbon has a KG classification of 'Csa', which constitutes a Mediterranean climate with dry and hot summers [32].Accordingly, and as identified by studies conducted by Miranda [33] and Calheiros [34], the urban microclimatic conditions present annual periods in which outdoor thermal comfort thresholds can be strained as a result of numerous of climatic occurrences, namely (i) between 10 and 20 days where daily maximum T amb surpasses 35 • C; (ii) an occurrence of between 100 and 120 'summer days' where daily maximum T amb exceeds 25 • C; and additionally (iii) frequent occurrences of heat waves, where daily maximum T amb exceeds that of 32 • C for various periods of successive days.
As identified in the study conducted by Alcoforado, Andrade [1], both North and Northwest wind directions are the most common throughout the year, particularly during the summer.This being said, due to the proximity to the Tagus, Lopes [35] identified that estuarine breezes reach adjacent urban areas on 30% of the late mornings and early afternoons during the summer.This study focuses predominantly upon the historical quarter of 'Baixa Chiado', which due to its general morphological composition has been identified to frequently witness the highest UHI intensities [2] and highest temperatures during the summer [6].

Data
To retrieve the base climatic data for the study, meteorological recordings were attained from the World Meteorological Organisation weather station with the Index N • 08535, located within Lisbon with the latitude of 38 • 43 N, 9 • 9 W, and an altitude of 77 m.As with similar existing studies [1,19,25,31,[36][37][38][39], the extracted data was converted into a variation of the thermo-physiological index, physiologically equivalent temperature (PET) [40][41][42].The respective index was used due to the (i) feasibility of being calibrated upon easily obtainable microclimatic parameters; and (ii) base measuring unit being • C, which facilitates its comprehension by professionals such as urban planners and designers when approaching climatological facets.
Before considering the influences of UCCs and vegetation, an initial analysis was established to determine general diurnal thermal conditions during different climatological circumstances, i.e., during the summer and winter within the city (Table 1).Such an approach enabled a 'base understanding' of Lisbon's thermo-physiological conditions as presented by the weather station.

July
3rd July December 16th December In order to carry out this initial assessment, the RayMan Pro © model [43,44] was used to process the following retrieved parameters: T amb , relative humidity (RH), total cloud oktas (Okt), and V 10 .In addition to these climatological aspects, the calibration of the RayMan model was constructed upon the default standing 'standardised man'; i.e., a height of 1.75 m, weight of 75 kg, aged 35, with a clothing 0.9 clo, and an internal heat production of 80 w [40,41].In order to account for the deceleration effect of 'urban roughness' upon wind speeds obtained from the meteorological station [45], the values were adapted to permit the estimation of V values upon the gravity centre of the human body.As a result, V 10 values were interpreted to a height of 1.1 m (henceforth expressed as V 1.1 ) by using the formula presented by Kuttler [46] (Equation (1)).Given Lisbon's denser downtown district with frequent open spaces, and similar to the morphological layout/compositions examined within comparable bioclimatic studies in Barcelona [38] and within the historical district of Lisbon [27], the study applied the following calibrations to the formula: z 0 = 1.00 m, and α = 0.35.
where V h is the m/s at a height of h (10 m), α is an empirical exponent, depending upon urban surface roughness, and z 0 is the corresponding roughness length.At a later stage in the study, another variable tuning was considered when addressing the 'in-situ' effects of vegetation upon pedestrian comfort as a result of evapotranspiration.As identified by the early studies of Oke [47], McPherson [48], and Brown and Gillespie [49], T amb and RH are usually not meaningfully modified by landscape elements such as trees since the encircling atmosphere quickly dissipates any such in-situ oscillations.Such results were also obtained by later studies [19,[50][51][52][53] who identified such limitations resulting from evapotranspiration effects, particularly from single trees.Moreover, such outcomes were also obtained by a recent study conducted within one of the widest public spaces located in Lisbon's historical centre.Carried out in July 2015, and as shown in Figure 1 the limited effects of evapotranspiration beneath the crown were also verified.
Atmosphere 2018, 9, 12 4 of 28 where Vh is the m/s at a height of h (10 m), α is an empirical exponent, depending upon urban surface roughness, and z0 is the corresponding roughness length.At a later stage in the study, another variable tuning was considered when addressing the 'in-situ' effects of vegetation upon pedestrian comfort as a result of evapotranspiration.As identified by the early studies of Oke [47], McPherson [48], and Brown and Gillespie [49], Tamb and RH are usually not meaningfully modified by landscape elements such as trees since the encircling atmosphere quickly dissipates any such in-situ oscillations.Such results were also obtained by later studies [19,50,[51][52][53] who identified such limitations resulting from evapotranspiration effects, particularly from single trees.Moreover, such outcomes were also obtained by a recent study conducted within one of the widest public spaces located in Lisbon's historical centre.Carried out in July 2015, and as shown in Figure 1 the limited effects of evapotranspiration beneath the crown were also verified.As a result, the lower impact that evapotranspiration could generate upon comfort conditions at pedestrian height was not considered in this particularly study.Instead, the assessments were predominantly based upon evaluating the in-sit' impacts of a specific tree species beneath its crown, or in other words "the cool feeling people experience as a result of a reduction in radiation" [15].Although such information has grown in recent years (as exemplified by numerous studies mentioned in this study), the authors argue that there is still a need to build upon estimations of radiation modifications as a result of vegetative deflection as suggested by Brown and Cherkezoff [55].As a result, such outputs should continue to build upon the growing comprehension of vegetation influences upon microclimatic characteristics at the pedestrian level [56].
Therefore, supplementary data modifications/calibrations were applied to assess the specific thermo-physiological in-situ effects of Tipuana tipu upon pedestrian comfort levels.In addition to the preliminary climatic variables, global radiation (Grad) retrieved by the authors were used to supplement the radiation flux assessments.These values were obtained through field surveys through the use of the handheld apparatus KKMOON SM206, with an accuracy of ±10 and a resolution of 0.1 W/m 2 .Similar to the approach carried out by Nouri and Costa [54], such measurements were established to obtain hourly oscillations of Grad (i) in specific locations of different UCCs; and (ii) beneath vegetative canopies to identify the amount radiation that was diluted by the vegetative crown mass.Such an approach enabled solar attenuation radiation levels to be determined specifically beneath the tree crowns, and later introduced within the biometeorological model.As a result, the lower impact that evapotranspiration could generate upon comfort conditions at pedestrian height was not considered in this particularly study.Instead, the assessments were predominantly based upon evaluating the in-sit' impacts of a specific tree species beneath its crown, or in other words "the cool feeling people experience as a result of a reduction in radiation" [15].Although such information has grown in recent years (as exemplified by numerous studies mentioned in this study), the authors argue that there is still a need to build upon estimations of radiation modifications as a result of vegetative deflection as suggested by Brown and Cherkezoff [55].As a result, such outputs should continue to build upon the growing comprehension of vegetation influences upon microclimatic characteristics at the pedestrian level [56].
Therefore, supplementary data modifications/calibrations were applied to assess the specific thermo-physiological in-situ effects of Tipuana tipu upon pedestrian comfort levels.In addition to the preliminary climatic variables, global radiation (G rad ) retrieved by the authors were used to supplement the radiation flux assessments.These values were obtained through field surveys through the use of the handheld apparatus KKMOON SM206, with an accuracy of ±10 and a resolution of 0.1 W/m 2 .Similar to the approach carried out by Nouri and Costa [54], such measurements were established to obtain hourly oscillations of G rad (i) in specific locations of different UCCs; and (ii) beneath vegetative canopies to identify the amount radiation that was diluted by the vegetative crown mass.Such an approach enabled solar attenuation radiation levels to be determined specifically beneath the tree crowns, and later introduced within the biometeorological model.

Applied Methodology and Structure
The modelling, refinement, and representation of the UCC assessments were processed through a three-step approach.Firstly, a new version of the SkyHelios model (i)  [29] was used in order to assess the microclimatic conditions and implications upon the introduced parameters with respect to the urban morphologies.As presented by Fröhlich [30], the model has been developed to analyse the spatial dimension of a specific microclimate.The short runtime and wide range of support input formats, as well as the capability of dealing with different projected coordinate systems permit the applicability of the model to evaluate the comparison of different UCCs, tree compositions and their effect on thermal biometeorology.In processing terms, the new SkyHelios model essentially follows a similar approach to that of the Rayman model, but it has the additional feature of also considering the spatial dimension of the parameters.As a result, SkyHelios makes use of the graphic processor for computing the sky view factor (SVF) [29].Based upon the SVF, the radiation fluxes can be determined for any location within the respective model area, and summarised as mean radiant temperature (T mrt ).As defined by Oke [47], the SVF is the fraction of the visible sky seen from a specified point.Within the software module, the first operation is the rendering of a fisheye image for that stipulated location, and secondly, the SVF is determined by distinguishing transparent and coloured pixels within the generated image.Similar to real environments, the fisheye image is a half-sphere, and not all of the pixels should have the same influence upon SVF.Therefore, a dimensionless weighting factor represented as ωproj (Equation ( 2)) is utilised to consider the projections that adjusts the impact of a pixel by the sine of the zenith angle ϕ (•).
This results in a spherical SVF, whereby if a planner SVF is desired, another correction by ωplanar (again dimensionless) needs to be performed (Equation (3)).Such a modification increases the impact of objects close to the ground by the cosine of the azimuth angle (counted from the ground to the top).
Furthermore, a three-dimensional diagnostic wind model was integrated within the updated version in order to provide estimations of 'in-situ' V measurements, which as discussed can differ significantly from those presented by the meteorological station.Such a model was based upon the approach conducted by Röckle [57], but, and as discussed by Fröhlich [30], with updated parameterizations, namely (i) an improved upwind cavity as discussed in Bagal, Pardyjak [58]; and (ii) an improved description of street canyon vortices as described in Singh, Hansen [59].As a result, the SkyHelios model enables the identification of spatially-resolved V and associated direction.The wind field calculated by the given functions most likely contains a certain degree of divergence.Assuming incompressible air, such a divergence has to be minimized in order to get a valid wind field.In mathematical terms, this is performed by minimizing the functional for the scalar H (Equation ( 4)).
As shown in Equation (4) b h and b v are horizontal stability factors in s/m; u, v and w are the stream components in m/s, u 0 , v 0 and w 0 are representative of the initial stream components in m/s while dx, dy and dz are the grid spacing in metres.Due to this integration, the model is capable of estimating thermal indices Perceived Temperature (PT) [60], Universal Thermal Comfort Index (UTCI) [61], and PET that are spatially resolved at high resolutions (e.g., 1 × 1 m).
The second step was orientated at processing the results through the updated version of the RayMan model (ii) in order to interpret them into the modified physiologically equivalent temperature (mPET) as discussed in Chen and Matzarakis [31].As presented in their study, the predominant differences of the mPET index are the integrated thermoregulation model (based upon a multiple-segment model), and the clothing model which relays a more accurate analysis of the human bio-heat transfer mechanism.For this reason, the initial examination of the diurnal conditions both during the summer and winter conditions were thus presented in both PET and mPET; to firstly examine such differences identified by Chen and Matzarakis [31], and secondly, present more accurate evaluations of the human thermal comfort conditions examined later in the study.
The third step of the study was to ensure that the study outputs were communicated in a way to facilitate their comprehension by non-experts in the subjects linked to climatology.As a result, the obtained results were processed with the climate tourism/transfer information scheme (iii)  (CTIS) [62].Such means of communication have been growing within the scientific community in similar climatic studies [31,38,[63][64][65].

Identification of Bioclimatic Conditions During Summer/Winter Periods
As discussed, the first stage was to identify the general bioclimatic conditions (i.e., PET and mPET values) retrieved from the specified weather station.As a result of this exercise, not only was it possible to identify how physiological stress (PS) grades vary during different times of the year, but moreover to identify the days on which the simulations would be based.Based upon the configuration of the meteorological station, since the total Okt values were only recorded at 09:00, 12:00, and 15:00, it was necessary to approximate values for the remaining hours between 09:00 and 18:00.Such an approximation was undertaken by (i) the delineation of the mid-range values between the three daily recordings; and, in addition, (ii) taking into account the qualitative meteorological descriptions provided by the station.Such inputs were accompanied by the introduced T amb , RH, and V 1.1 diurnal recordings for the months of July and December, each representative of the hottest and coldest annual conditions for 2016 found in Lisbon (Table 1).Furthermore, and at this stage of the study, the applied SVF value was calibrated at 1.00 (or 100%) which would reproduce an assessment with total exposure to solar radiation.Once the PET and mPET values were processed, the comparative chart presented by Matzarakis, Mayer [42] was used in order to categorise the obtained temperature values into specific PS levels, shown in Table 2.
Table 2. Ranges of the thermal index physiologically equivalent temperature (PET) for different grades of physiological stress (PS) on human beings; internal heat production: 80 W, heat transfer resistance of the clothing: 0.9 clo according to [66] | Source: Adapted from, [42].Once the boundaries for each temperature range were stipulated, it was possible to input such margins into CTIS to facilitate the representation of the diurnal modifications of PS during the months of July and December.Such methods of representation have already been used in similar bioclimatic studies [31,38,[63][64][65] which also aimed at facilitating the comprehension of the obtained thermo-physiological results.

Defining Tipuana tipu Characteristics and Layout
The Tipuana tipu species originates from the region adjacent to the Tipuani River in Bolivia, and is also known for its lineages from both Brazil and Argentina.Given the right conditions, the species can reach considerable sizes, and presents a dark trunk, twisted branches, and a densely vegetated round crown.Resulting from the 'Csa' climate, the equally ornamental and rustic species are particularly compatible with Lisbon's climate.Consequently, they are one of the most commonly used deciduous species [67], with approximately 1900 examples within the city [68,69].In a few locations, a particular few have been listed as 'public interest' due to their considerable age and size as exemplified in 'Cais Sodré', 'São Bento' Plaza, and the 'Nove de Abril' Park, registering ages of 100, 80, and 128, respectively.Recognised as an effective urban shading tree specifically within Mediterranean climates [17,54,67,[69][70][71][72] it is also frequently mixed with other species within Lisbon.Such circumstances can be exemplified by a mixture with other identified 'shading trees' such as Jacaranda mimosifolia within the urban garden of 'Santos' and the public plaza of 'Rossio' within the historical quarter.
The example of the latter is illustrated in Figure 2A, which shows a lined sidewalk seating area beneath the vegetative crowns of a linear plantation of Tipuana tipu with a parallel row of Jacaranda mimosifolia on the edge of the sidewalk.Common to sidewalk plantation typologies as identified by Torre [70], the linear plantation configuration has been identified to effectively attenuate (i) radiation, especially when the tree lines are parallel/close to a building frontage [48,49,73,74]; and (ii) wind patterns that take place perpendicularly to the tree line [70,72,75].

Defining Tipuana tipu Characteristics and Layout
The Tipuana tipu species originates from the region adjacent to the Tipuani River in Bolivia, and is also known for its lineages from both Brazil and Argentina.Given the right conditions, the species can reach considerable sizes, and presents a dark trunk, twisted branches, and a densely vegetated round crown.Resulting from the 'Csa' climate, the equally ornamental and rustic species are particularly compatible with Lisbon's climate.Consequently, they are one of the most commonly used deciduous species [67], with approximately 1900 examples within the city [68,69].In a few locations, a particular few have been listed as 'public interest' due to their considerable age and size as exemplified in 'Cais Sodré', 'São Bento' Plaza, and the 'Nove de Abril' Park, registering ages of 100, 80, and 128, respectively.Recognised as an effective urban shading tree specifically within Mediterranean climates [17,54,67,[69][70][71] it is also frequently mixed with other species within Lisbon.Such circumstances can be exemplified by a mixture with other identified 'shading trees' such as Jacaranda mimosifolia within the urban garden of 'Santos' and the public plaza of 'Rossio' within the historical quarter.
The example of the latter is illustrated in Figure 2A, which shows a lined sidewalk seating area beneath the vegetative crowns of a linear plantation of Tipuana tipu with a parallel row of Jacaranda mimosifolia on the edge of the sidewalk.Common to sidewalk plantation typologies as identified by Torre [70], the linear plantation configuration has been identified to effectively attenuate (i) radiation, especially when the tree lines are parallel/close to a building frontage [48,49,73,74]; and (ii) wind patterns that take place perpendicularly to the tree line [70,72,75].Within urban pavements, the species often present smaller sizes [17,67,72,76] which is often due to inferior soil conditions and embedding methods.Illustrated in Figure 2B, and based upon the aforementioned studies, the dimension of the trees used within the study were calibrated with a height (h) of 12.5 m, vegetative crown radius (r) of 5.0/6.5 m, trunk length (l) of 3 m, and a trunk Within urban pavements, the species often present smaller sizes [17,67,72,76] which is often due to inferior soil conditions and embedding methods.Illustrated in Figure 2B, and based upon the aforementioned studies, the dimension of the trees used within the study were calibrated with a height (h) of 12.5 m, vegetative crown radius (r) of 5.0/6.5 m, trunk length (l) of 3 m, and a trunk diameter (d) of 0.55 m.Furthermore, since G rad values were measured beneath the crowns, parameters such as leaf area index and leaf angle distribution were not required to determine the vegetative density/transmissivity of the tree crown.
Beyond being one of the most commonly used 'shading' tree species, and unlike most deciduous species, it also has the particularity of maintaining its foliation both during the summer and winter period.As demonstrated in Figure 2C, and as identified by Viñas, Solanich [72] within the Iberian Peninsula, the only period in which the Tipuana tipu loses its leaves is during the early spring.Accordingly, and given the permanency of its vegetative crown during the winter, it was also possible to assess how the vegetative crown would influence 'in-situ' thermal conditions during colder conditions.As identified by early studies [15,49,55,77], it was possible to examine the consequences of reductions of radiation fluxes during the winter; a known essential microclimatic variable with both physiological and psychological connotations to counteract colder conditions during the winter in outdoor urban public spaces [78].

Application and Construction of the SkyHelios Simulations
Within the 'Obstacle' plugin associated with both the RayMan and the SkyHelios models, various morphological compositions were constructed.Based upon modifying one of the characteristics of the canyons, it was possible to obtain dissimilar height-to-width (HW) ratios which were the most commonly found within Lisbon's historical district.Accordingly, the canyon height was maintained at 20 m (equating roughly to a building of five stories), and the canyon widths fluctuated between 10, 20, 40, and 120 metres to obtain 'very high', 'high', 'medium', and 'very low' UCCs, respectively.Although focused upon a larger scale, a similar approach was applied by Norton, Coutts [79], who categorised broad priority requirements for green infrastructure within different urban canyons.Previously, an UCC with a description of 'low' was also initially considered for the study.Such a canyon presented a canyon width of 80 m, and a HW ratio of 0.25, yet based upon a previous study conducted by the authors, it was identified that it presented almost identical thermal conditions to the 0.17 ratio.As the latter presented greater dissimilarities to the higher 0.50 ratio, it was decided to focus upon the four identified UCCs within this specific study as presented in Table 3.In accordance with the proportional and rigid morphological composition succeeding the reconstruction of Lisbon's historical district following the great earthquake of 1755, all modelled canyons were symmetrically configured.Although such an approach is frequently 'presumed' when approaching H/W ratios, recent microscale bioclimatic studies [23,26] have also produced important results in asymmetrical canyons whilst also addressing urban thermal comfort conditions.
Through the use of the SkyHelios model, the HW 2.00 , HW 1.00 , HW 0.50 , and HW 0.17 were processed under two conditions: (i) without the presence of vegetation; and (ii) with the presence of vegetation.In addition, each canyon was aligned into a north-to-south orientation (NSO) and a west-to-east orientation (WEO) in order to assess the influence of the geo-referenced summer/winter sun path upon the two alignments.The layout of the NSO simulations with the presence of the Tipuana tipu is represented in Figure 3 indicates how, depending upon the UCCs, the tree layouts were structured within each undertaken assessment.Such an adjustment was based upon maintaining a similar amount of vegetative coverage ratio (VCR) throughout the width of each UCC.Within this study, such an indicative value was obtained by applying the straightforward formula as shown in Equation (5).
where W is the width of the aspect ratio, CS is the crown spread, n is the number of trees, and r is the radius (5, 6.5).
Atmosphere 2018, 9, 12 9 of 28 vegetative coverage ratio (VCR) throughout the width of each UCC.Within this study, such an indicative value was obtained by applying the straightforward formula as shown in Equation 5.
where W is the width of the aspect ratio, CS is the crown spread, n is the number of trees, and r is the radius (5, 6.5).Although the VCR was based upon determining the amount of vegetation against the width of the canyons, comparable ratios were also utilised within similar studies [17,56,80] who quantified the amount of vegetative biomass within an entire outdoor area.Based upon the discussed parameters in Figure 2B, the dimensions of the introduced trees were calibrated with a diameter of Although the VCR was based upon determining the amount of vegetation against the width of the canyons, comparable ratios were also utilised within similar studies [17,56,80] who quantified the amount of vegetative biomass within an entire outdoor area.Based upon the discussed parameters in Figure 2B, the dimensions of the introduced trees were calibrated with a diameter of 13 m, with the exception for HW 2.00 , which accommodated trees with a lower diameter of 10 m.Such a reduction was undertaken to account for the shorter width of the street, whilst respecting the smaller, yet feasible dimensions of the Tipuana tipu.
As presented within the study conducted by Lin, Tsai [81], the calculation of SVF within this study was based upon the classic single-point SVF within a specific location in the canyon to obtain a fisheye image at a calibrated height of 1.1 m.Similar to an approach applied by Algeciras, Consuegra [25] both the NSO and WEO accommodated three specific reference points (RPs) which were established to identify the diurnal oscillations of solar radiation.As shown in Table 4, the RPs were attributed a specific location within the different UCCs, and their respective coordinates within the SkyHelios model are presented and discussed in Table 5.Additionally, and even though canyon length was un-associated to the HW ratio, each was attributed a length of 200 m in order ensure that 'edges' of the default canyon would not meaningfully affect the obtained SVF values obtained by the SkyHelios.13 m, with the exception for HW2.00, which accommodated trees with a lower diameter of 10 m.Such a reduction was undertaken to account for the shorter width of the street, whilst respecting the smaller, yet feasible dimensions of the Tipuana tipu.
As presented within the study conducted by Lin, Tsai [81], the calculation of SVF within this study was based upon the classic single-point SVF within a specific location in the canyon to obtain a fisheye image at a calibrated height of 1.1 m.Similar to an approach applied by Algeciras, Consuegra [25] both the NSO and WEO accommodated three specific reference points (RPs) which were established to identify the diurnal oscillations of solar radiation.As shown in Table 4, the RPs were attributed a specific location within the different UCCs, and their respective coordinates within the SkyHelios model are presented and discussed in Table 5.Additionally, and even though canyon length was un-associated to the HW ratio, each was attributed a length of 200 m in order ensure that 'edges' of the default canyon would not meaningfully affect the obtained SVF values obtained by the SkyHelios.

2.00
Central and lateral RPs were varied little due to the limited width of the canyon

1.00
Central RPs were placed beneath tree crown, while the lateral RPs were situated amid the edge of the tree crown and the lateral façades of the canyon

0.50
Central RPs were placed between the edges of two tree crowns, and the lateral RPs placed between the edge of the tree crown and the lateral façades of the canyon 0.17 Central RPs were slightly offset to be directly beneath a tree crown, and lateral RPs which were also placed beneath the trees located adjacent to the façades of the canyon Based upon the calibration of a three-dimensional diagnostic wind model integrated within the program, the last of the eight-tree line was selected to permit the identification of the influence of preceding trees in UCCs (Figure 3).As a result, it was possible to determine how such obstacles within each canyon would influence V1.1 once it reached each of the designated RPs.Configured within the diagnostic tool within SkyHelios (Table 6), the initial V1.1 values obtained through Equation 1 could be recalibrated.Based upon the calibration of a three-dimensional diagnostic wind model integrated within the program, the last of the eight-tree line was selected to permit the identification of the influence of preceding trees in UCCs (Figure 3).As a result, it was possible to determine how such obstacles within each canyon would influence V 1.1 once it reached each of the designated RPs.Configured within the diagnostic tool within SkyHelios (Table 6), the initial V 1.1 values obtained through Equation ( 1) could be recalibrated.Table 6.Specific configurations of SkyHelios wind diagnostic tool.

0.50
Central RPs were placed between the edges of two tree crowns, and the lateral RPs placed between the edge of the tree crown and the lateral façades of the canyon 0.17 Central RPs were slightly offset to be directly beneath a tree crown, and lateral RPs which were also placed beneath the trees located adjacent to the façades of the canyon Based upon the calibration of a three-dimensional diagnostic wind model integrated within the program, the last of the eight-tree line was selected to permit the identification of the influence of preceding trees in UCCs (Figure 3).As a result, it was possible to determine how such obstacles within each canyon would influence V1.1 once it reached each of the designated RPs.Configured within the diagnostic tool within SkyHelios (Table 6), the initial V1.1 values obtained through Equation 1 could be recalibrated.

Disparities between July and December Bioclimatic Conditions
Demonstrated in Figure 4, the processed microclimatic parameters of T amb , RH, V 1.1 , and Okt obtained from the meteorological station were translated into PS levels for both July and December.As expected, between these two months, there was a large variation of diurnal thermal stress.
During July 2016, PS grades predominantly varied between 'slight heat stress' and 'extreme Heat stress'.In addition, and due to the thermo-physiological indices reaching values which considerably surpassed that of 41 • C, an 'extreme' PS grade was added ('extreme heat stress 2') to the CTIS calibration to effectively present results which reached 46 • C. Consequently, such an adjustment also relays back to the possibility for expanding the existing ranges of PS upon human beings as presented in Table 2.Such PS stress levels were obtained during days with particularly elevated T amb values.More precisely, the occurrence of 'extreme heat stress 2' took place under two extreme heat events: (1) during a sequential set of days where maximum diurnal T amb values surpass that of 32 • C-described as an 'urban heat wave' by Calheiros [34]; and (2) during a specific day where maximum diurnal T amb values surpass that of 35 • C-described as a 'very hot day' by Miranda [33].As revealed in Figure 4, July 2016 witnessed three urban heat waves and very hot days, whereby most extreme events were interlinked with one another, with the exception of the 3rd of July.During this specific day, recorded T amb values almost constantly remained at 35 • C with low V 1.1 values ranging between 1.0 m/s and 2.0 m/s between 12:00 and 16:00 (Table 1).
When considering the results for December 2016, it was evident that the PS thresholds varied differently than those observed for the summer period.More concretely, although only heat stress was observed during July, December revealed both types of thermal stress, particularly during the first half of the month.Such results imply that, whilst Mediterranean climates do witness rare occasions of considerable cold stress (as exemplified during the mornings of the 23rd and 27th), the majority of PS tends to vary between 'moderate cold stress' and 'no thermal stress'.Nonetheless, it was also noted that, particularly between the hours of 12:00 and 15:00, there were also periods surpassing 'slight heat stress', as exemplified by the 6th and 10th-12th of December.These outcomes indicate that Lisbon's thermo-physiological conditions during the winter fluctuate considerably less from thermal comfortable conditions than those obtained during the summer.Regardless it was noted the occurrence of days where, even between 12:00 and 15:00, PS grades still presented 'Moderate Cold Stress' as was the case during the 16th of December.Such conditions took place due to a combination of microclimatic characteristics, such as lower diurnal T amb values (i.e., oscillating between 8.9 • C and 12.4 • C), and higher V 1.1 values (i.e., oscillating between 2.0 m/s and 4.1 m/s) (Table 1).When considering the results for December 2016, it was evident that the PS thresholds varied differently than those observed for the summer period.More concretely, although only heat stress was observed during July, December revealed both types of thermal stress, particularly during the first half of the month.Such results imply that, whilst Mediterranean climates do witness rare

Comparisons between PET and mPET Indices
The inclusion of both thermo-physiological indices permitted (i) the comparison of PET results against those obtained by a newly adapted mPET index; and (ii) the scrutiny of whether the application of mPET in this study would also collaborate with the general the results obtained by Chen and Matzarakis [31] for western European climates.When cross-examining the outcomes of the two indices for July and December in Figure 4, two major distinctions were acknowledged.
Firstly, in the case of July, it was possible to verify that mPET values very rarely reached a PS of 'extreme heat stress'.In the sporadic case when mPET did reach such a PS threshold (such as 15:00 during the 3rd of July), it was associated to an extreme heat event (a very hot day).Yet, even under such conditions, and unlike the results obtained by the PET index, the adapted index was distant from 'extreme heat stress 2'.Concordant with these results, and as discussed within the original study, "almost no extreme hot events were given by the estimation of mPET and only moderate heat stress occurred [ . . .] on the contrary, PET [led to] the occurrence of extreme heat stress in Freiburg during the summer" ( [31], p. 7).
Secondly, and focusing upon the month of December, it was demonstrated that overall PS values were closer to the PS grade of 'no thermal stress' with the use of the mPET index, as opposed to the PS levels obtained by the PET index.Again, such outcomes were also obtained within the original study in Freiburg.Similar to the divergence between the indices for the 3rd of July, colder days such as the 16th of December revealed also lower stress levels through the estimations obtained by mPET.
For these reasons, and based upon the more accurate analysis of the human bio-heat transfer system when exposed to thermal stimuli, the mPET index was used throughout the rest of the study.Such a decision by far devalued the use of the PET index in such studies; rather, it proposed settings in which the preceding index may be complemented b: (i) new improvements to the original Munich energy balance model for individuals model parameters; and (ii) an adjustment of the existing ranges of PS upon human beings.

Day Selection for simulations
Based upon the outputs obtained for both July and December 2016, it was possible to identify the days in which to undertake the simulations within the SkyHelios model.Beyond providing the daily input variables for the particular days, this approach also enabled specific Summer/Winter sun paths to be calibrated into the model.As the intention was to examine the influence of the Tipuana tipu species during annual extreme conditions during the different times of year, both the hottest and coldest days were selected for the simulation of the UCCs.Based upon the identification of the monthly variations of diurnal PS, the 3rd of July and the 16th of December were chosen to represent extreme summer and winter conditions, respectively.

Canyon RP Outcomes
As identified by existing studies, thermal comfort within urban contexts is strongly affected by urban configurations, whereby morphological elements such as width, height and orientation of a specific canyon are imperative to evaluate precise microclimatic conditions [25,30,[82][83][84][85]. Additionally, and through SVF assessments, the impact of vegetation has also been identified as exemplified by Charalampopoulos, Tsiros [86], who determined that, particularly with low SVF values, the presence of vegetation presented positive influences upon thermal comfort conditions.Analogous results obtained by Shashua-Bar, Potchter [17] also indicated such an intrinsic relationship between the thermal effect of trees, with that of urban street geometry.
Undertaking a similar approach to the aforementioned studies, and processed by the SkyHelios model, Figures 5 and 6 present the obtained SVF values for the UCCs.Moreover, based upon the stipulation of the summer paths, it was possible to ascertain the amount of time that each RP was exposed to the sun between 09:00 and 18:00.Such variables were expressed as M Sum and M Win to portray the number of minutes each RP was cast in the sun, during the summer and winter, respectively.As expected, M Sum at each RP was higher given the absence of vegetation.M Win did not follow this tendency due to its lower winter sun path with the exception for 0.17 RP N which revealed a reduction in M Win of 270.
Atmosphere 2018, 9, 12 14 of 28 values, the presence of vegetation presented positive influences upon thermal comfort conditions.Analogous results obtained by Shashua-Bar, Potchter [17] also indicated such an intrinsic relationship between the thermal effect of trees, with that of urban street geometry.Undertaking a similar approach to the aforementioned studies, and processed by the SkyHelios model, Figure 5 and Figure 6 present the obtained SVF values for the UCCs.Moreover, based upon the stipulation of the summer paths, it was possible to ascertain the amount of time that each RP was exposed to the sun between 09:00 and 18:00.Such variables were expressed as MSum and MWin to portray the number of minutes each RP was cast in the sun, during the summer and winter, respectively.As expected, MSum at each RP was higher given the absence of vegetation.MWin did not follow this tendency due to its lower winter sun path with the exception for 0.17 RPN which revealed a reduction in MWin of 270.Between the two UCCs orientations, WEO generally resulted in higher durations of solar exposure.Yet, such durations varied between the RPs within the specific canyons.As exemplified by the WEO, considerable differences in MSum were identified, as exemplified by 2.00 RPS obtaining the MSum of 130, when 2.00 RPC2 and 2.00 RPN revealed higher MSum values of 500 and 480, respectively.In the case of HW0.50, which was the only UCCs in which the central RPs were not placed directly beneath a tree crown, 0.50 RPC1 and 0.50 RPC2 still revealed considerable reductions in MSum, particularly 0.50 RPC1.Nevertheless, the highest reductions of MSum were obtained in HW0.17, within all RPs, with a maximum reduction witnessed at 0.17 RPC1 where MSum was reduced by 470.

Canyon mPET Outcomes
Within this section, the results obtained from SkyHelios and subsequently processed by RayMan to obtain mPET values, and respective PS grades, were demonstrated for each of the four UCCs (Figure 7/Figure 8).In order to facilitate the comprehension of the results, the CTIS program was used to present the (i) differences obtained between the simulations with no vegetation (NVSim), Between the two UCCs orientations, WEO generally resulted in higher durations of solar exposure.Yet, such durations varied between the RPs within the specific canyons.As exemplified by the WEO, considerable differences in M Sum were identified, as exemplified by 2.00 RP S obtaining the M Sum of 130, when 2.00 RP C2 and 2.00 RP N revealed higher M Sum values of 500 and 480, respectively.In the case of HW 0.50 , which was the only UCCs in which the central RPs were not placed directly beneath a tree crown, 0.50 RP C1 and 0.50 RP C2 still revealed considerable reductions in M Sum , particularly 0.50 RP C1 .Nevertheless, the highest reductions of M Sum were obtained in HW 0.17 , within all RPs, with a maximum reduction witnessed at 0.17 RP C1 where M Sum was reduced by 470.

Canyon mPET Outcomes
Within this section, the results obtained from SkyHelios and subsequently processed by RayMan to obtain mPET values, and respective PS grades, were demonstrated for each of the four UCCs (Figure 7/Figure 8).In order to facilitate the comprehension of the results, the CTIS program was used to present the (i) differences obtained between the simulations with no vegetation (NV Sim ), against simulations with the presence of the Tipuana tipu species (TT Sim ); and, (ii) identified disparities from the results as initially obtained from the meteorological station.As anticipated, the UCCs with the least amount of thermal stress were HW 2.00 and HW 1.00 which generally presented both lower periods and quantities of thermo-physiological stress.Such was predominantly recognised within the NV Sim .More precisely, and with regards to the HW 2.00 canyon, it was established that during the:

•
3rd of July within the NSO, the presence of the tree crowns was able to reduce PS levels between 12:00-14:00 by one whole stress grade with a maximum mPET reduction of 6.8 • C (obtained at 14:00 in RPE).Such a decrease represented the capacity of the tree crown to decrease the vulnerability at pedestrian height during the few hours in which the canyon became exposed to direct solar radiation; • 3rd of July within the WEO, the TT Sim presented much more prolonged reductions of PS thresholds, especially in RP N and RP C2 which were cast in the shade of the vegetative crown between 10:00-15:00 and 10:00-16:00, respectively.During these hours, and in the case of RP C2 , the reduced mPET values ranged between 29.1 • C-37.7 • C, instead of 34.6 • C-42.9 • C as obtained in the NV Sim ; • 16th December for both the NSO and WEO, it was possible to identify that the absence/presence of vegetation led to no significant differences between the NV Sim and the TT Sim .In the case of the NSO, the very subtle variations of PS such as those presented at 12:00 and 16:00 were a result of V 1.1 values oscillating by 0.1 m/s.Within the WEO such distinctions were even less significant, yet as a result of generally higher V 1.1 values (i.e., between 0.2 m/s and 0.7 m/s), all RPs presented slightly colder PS levels.
Within the HW 1.00 canyon, it was determined that during the: • 3rd of July within the NSO, results were not excessively dissimilar from those obtained in HW 2.00 .However, within this UCC, PS levels increased both in stress intensity and duration.Both lateral RPs revealed periods of 'extreme heat stress' even in the TT Sim as exemplified between 12:00-13:00 for RP W , and at 13:00 for RP E .Both of these circumstances took place when the sun path crossed the gap between the crown and that of the western/eastern canyon facades.Nevertheless, the diurnal mPET results obtained in the NV Sim presented much higher PS grades when the canyon was directly exposed to radiation fluxes.Also obtained at RP E , a maximum reduction in mPET of 7.7 • C was identified when the PS in the NV Sim almost reached 'extreme heat stress 2'; • 3rd of July within the WEO, the results obtained in RP N show the influence that the vegetative shade could present within the RP until 15:00.At this specific time, both mPET values for the NV Sim and TT Sim reached 43.6 • C and 44.1 • C, respectively, thus leading to equally extreme PS grades, which gradually subsided by 18:00; • 16th December for both the NSO and WEO, and similar to the HW 2.00 canyon, there were very limited variations of thermal conditions between the NV Sim and TT Sim .Nevertheless, it was noted that the WEO presented much more stable PS grades between in RP N/C2/S in comparison to those within the NSO.The justification for this was attributed to the considerable variations of V 1.1 conditions in the entire canyon which were later presented in Figure 9.   Within the HW1.00 canyon, it was determined that during the: • 3rd of July within the NSO, results were not excessively dissimilar from those obtained in HW2.00.However, within this UCC, PS levels increased both in stress intensity and duration.Both lateral RPs revealed periods of 'extreme heat stress' even in the TTSim as exemplified between 12:00-13:00 for RPW, and at 13:00 for RPE.Both of these circumstances took place when the sun path crossed the gap between the crown and that of the western/eastern canyon facades.Nevertheless, the diurnal mPET results obtained in the NVSim presented much higher PS grades when the canyon was directly exposed to radiation fluxes.Also obtained at RPE, a maximum reduction in mPET of 7.7 °C was identified when the PS in the NVSim almost reached 'extreme heat stress 2'; • 3rd of July within the WEO, the results obtained in RPN show the influence that the vegetative shade could present within the RP until 15:00.At this specific time, both mPET values for the   Analogous to existing studies [1,25,30], the disparity of thermal conditions between those presented from the meteorological station and those obtained from the in-situ assessments were substantial, even in the case of the NVSim.Therefore, such outcomes continue to infer the risk of solely considering station values for microscale and studies/projects.On the other hand, these studies also present the methodical means in which to adapt such results into useful data that can be much more valuable for local/microscale thermal sensitive urban design and planning.
Regardless of the considerable differences, there were however associations between the results obtained from the station and those from the 'in-situ' assessments.For example, within most of the UCCs it was possible to identify the influence of the accentuated climatic circumstances as initially identified by the station, namely: (i) the elevated Tamb (35.9 °C) and simultaneously low V10 (2.0 m/s) When considering the results for the HW 0.50 and HW 0.17 , it was verified that overall PS levels increased, and the TT Sim presented even stronger attenuation results during the summer.Within the HW 0.50 canyon, during the

•
3rd of July within the NSO, due to the increased width of the canyon, the susceptibility to solar radiation was augmented in both the NV Sim and the TT Sim .Within the NV Sim , elevated PS levels increased in duration in comparison to HW 1.00 , especially during the morning in RP W , and the afternoon in RP C1 .Given the placement of the RP C1 between the edges of two tree crowns within the TT Sim , it witnessed similar PS grades between 12:00-14:00 to those obtained in the NV Sim ; • 3rd of July within the WEO, RP N and RP C2 presented similar occurrences within the TT Sim .
Whilst cast in the shade within the middle of the canyon, mPET values decreased by up to 6.3 • C as illustrated at 13:00 in RP C2 .When not cast in the shade, however, even the TT Sim presented PS thresholds reaching 'extreme heat stress 2' with mPET values reaching 45.2 • C. Such a value actually exceeded slightly the mPET value of 44.9 • C obtained in NV Sim as a result of a slightly higher V 1.1 value (+0.2 m/s); • 16th December for both the NSO and WEO, there still were no clear variations between PS grades between the NV Sim and the TT Sim .It was however noted that, between all RPs, the PS grades within RP C1 and RP E revealed higher cold stress as a result of higher diurnal V 1.1 values within these two locations (Figure 9).
Lastly, and within the HW 0.17 canyon, it was determined that during the:  Analogous to existing studies [1,25,30], the disparity of thermal conditions between those presented from the meteorological station and those obtained from the in-situ assessments were substantial, even in the case of the NV Sim .Therefore, such outcomes continue to infer the risk of solely considering station values for microscale and studies/projects.On the other hand, these studies also present the methodical means in which to adapt such results into useful data that can be much more valuable for local/microscale thermal sensitive urban design and planning.
Regardless of the considerable differences, there were however associations between the results obtained from the station and those from the 'in-situ' assessments.For example, within most of the UCCs it was possible to identify the influence of the accentuated climatic circumstances as initially identified by the station, namely: (i) the elevated T amb (35.9 • C) and simultaneously low V 10 (2.0 m/s) obtained at 15:00 during the 3rd of July 2016; and also, (ii) the low T amb (9.8 • C) and simultaneously high V 10 (8.0 m/s) obtained at 13:00 during the 16th of December 2016.

Wind Influences
As expected, the lower the UCC, the higher the variation of V 1.1 between the different RPs within the canyons.Such variations are presented in Figure 9, which summarises the NSO canyon results obtained through the use of the three-dimensional tool integrated within SkyHelios during the winter period.When considering the modifications of mPET/PS during December, it was possible to make direct associations between the thermal comfort conditions and the V 1.1 results.
When considering the NSO canyons in December, the variations of PS between the RPs (shown in Figures 7 and 8) can be directly linked to the variations of V 1.1 shown in Figure 9.More specifically, to the: (1) HW 2.00 canyon presenting almost identical values at each RP; (2) HW 1.00 canyon revealing the highest values in RP E , in-between values in RP C1 , and the lowest in RP W ; (3) HW 0.50 revealing equally higher values in RP E and RP C2 , with lowest values in RP W ; and lastly, (4) HW 0.17 again revealing greater values in RP E and RP C2 , and RP W once again revealing the lowest values amongst the three.Thus, it was possible to identify why almost all RP W PS levels were slightly closer to 'comfortable conditions' in comparison to the other two locations in the canyon.
Finally, when considering the vacillation of V 1.1 values between the NV Sim and the TT Sim , it was identified that the estimated effects upon wind currents at pedestrian height were fairly inconsequential.Such a result can be attributed to the output of the wind simulations being set at a height of 1.1 m, which fell well beneath the calibrated trunk height (l) of 3 m.Nevertheless, such a result calls for a future study into the effects of vegetation when they are configured closer to one another (i.e., in a cluster) rather than spaced at a considerable distance as calibrated in the TT Sim .Such a study would also certainly render very different results upon the radiation fluxes at pedestrian height as well.

Tmrt Influences
Contrasting the V 1.1 results, the NV Sim and the TT Sim demonstrated very clear deviations of T mrt in all canyons.As a result of the tree crowns, the casted shade within the RPs led to considerable reductions of radiation fluxes at the pedestrian level.
More specifically, Figure 10 exemplifies such variations of T mrt within the northern area of the different WEO canyons, whereby (i) even in the case of HW 2.00 (with a canyon width of only 10 m) the crown of the Tipuana tipu was able to reduce T mrt by 15.6 • C; (ii) in both HW 1.00 and HW 0.50 which presented a similar SVF value, similar T mrt reductions were of 13.1 • C and 13.3 • C, respectively; and lastly, (iii) the largest difference between the four locations was HW 0.17 which presented a dramatic T mrt decrease of 29.5 • C. Subsequently, and when referring to the obtained mPET differences between the NV Sim and the TT Sim , it was possible to identify the strong correlation with these T mrt results in each specific RP.When observing the Tmrt deviations during the winter between the two simulations, smaller, yet important differences were also acknowledged.As shown in Figure 11, due to the Tipuana tipu species maintaining its foliation during the winter, the vegetative crown also reduced radiation fluxes at pedestrian level during December.As a result, in some circumstances (particularly in lower UCCs) this led to the TTSim presenting lower comfort conditions.Within the HW0.17 canyon, and principally in RPN and RPC2, it was possible to verify direct results at 13:00 upon the obtained PS grades, even with small reductions of 5.5 °C and 4.6 °C, respectively.When observing the T mrt deviations during the winter between the two simulations, smaller, yet important differences were also acknowledged.As shown in Figure 11, due to the Tipuana tipu species maintaining its foliation during the winter, the vegetative crown also reduced radiation fluxes at pedestrian level during December.As a result, in some circumstances (particularly in lower UCCs) this led to the TT Sim presenting lower comfort conditions.Within the HW 0.17 canyon, and principally in RP N and RP C2 , it was possible to verify direct results at 13:00 upon the obtained PS grades, even with small reductions of 5.5 • C and 4.6 • C, respectively.When observing the Tmrt deviations during the winter between the two simulations, smaller, yet important differences were also acknowledged.As shown in Figure 11, due to the Tipuana tipu species maintaining its foliation during the winter, the vegetative crown also reduced radiation fluxes at pedestrian level during December.As a result, in some circumstances (particularly in lower UCCs) this led to the TTSim presenting lower comfort conditions.Within the HW0.17 canyon, and principally in RPN and RPC2, it was possible to verify direct results at 13:00 upon the obtained PS grades, even with small reductions of 5.5 °C and 4.6 °C, respectively.

Discussion
When considering the results obtained by this study, it became clear that the use of the updated SkyHelios version provided an effective tool to (i) downscale results obtained from Lisbon's meteorological station and render more accurate outputs for the microclimatic examinations as exemplified by this investigation; and (ii) estimate the in-sit' effects of the Tipuana tipu species by considering a multitude of climatic variables.Such an examination was complemented with the use of onsite measurements of G rad in order to attribute more precision to the results obtained by the simulations within each UCC.
Naturally, the results obtained from the model should be approached as estimations of what would take place within actual conditions, particularly with regards to the V 1.1 presented by the new integrated wind analysis tool.As a result, this raises the opportunity to undertake further field studies to further elaborate the results presented by this article.Nevertheless, it was noted that the output speeds during July within the lower UCCs were very similar to those obtained by the field measurements obtained by Nouri and Costa [54] during the summer of 2015 in Lisbon's historical quarter.In addition, the variation of different speeds within the different regions of the canyon were also analogous, whereby (i) the western regions of the NSO canyon also revealed the lowest wind speeds; and (ii) the central and eastern regions respectively revealed the highest values.Furthermore, and when considering the relationship of all outputs from the NV Sim and the TT Sim , it was possible to see clear associations between the different microclimatic variables, namely between the (i) quantity of M Sum /M Win values associated to each RP, and the general diurnal thermal conditions within that location; (ii) variations of thermal conditions suggested by the meteorological station, and those obtained within each assessed canyon/location; and, (iii) the impact of climatic variables such as T amb and V 10/1.1 upon the overall thermo-physiological projections.
The application of the new mPET index within this study also provided the opportunity to both compare and apply a newer PET version based upon an improved thermoregulation model.The results obtained in this study were concomitant with those revealed by Chen and Matzarakis [31], whereby the initial identification of bioclimatic conditions (PET and mPET variables) during the month of July and December did not reveal PS values exceeding that of 'Extreme Heat Stress'.Although later in the study the mPET results did in fact lead to higher PS grades within the NV Sim and the TT Sim , such exacerbations were attributed to the identification of dramatic 'in-situ' microclimatic occurrences such as when G rad values exceeded that of 900 W/m 2 , and other climatic variables such as T amb were particularly high (such as at 15:00 during the 3rd of July).Such extreme heat conditions, as a result, point to the need to improve the prediction and prevention of such events particularly for the context of Western Europe [87,88].
Within this study, numerous opportunities for future study were identified, namely (1) how mPET/PET entirely based upon variables obtained from the field surveys could differ and/or complement those obtained in this study, especially V values as identified in Nouri and Costa [54]; (2) the opportunity to consider how a group/cluster of trees would present different results, including possible oscillations of V between the NV Sim and the TT Sim [76]; (3) the expansion and of PS thresholds to situate mPET/PET results beyond existing levels (e.g., that of 41 • C) as discussed in Nouri, Lopes [10]; (3) following from the previous point, to moreover investigate how mPET outputs may need a general recalibration of existing PS boundaries as already indicated by Chen and Matzarakis [31]; and, lastly (iii) consider how other types of public space design measures, such as misting systems, can be incorporated within vegetative crowns to further augment thermal comfort through the effects of evaporative cooling [89,90].
As stipulated, the 3rd of July and the 16th of December were selected due to their particularly high estimated diurnal stress levels in order to evaluate the potential 'in-situ' attenuation effects from the Tipuana tipu.Both days revealed important results, whereby: (i) the summer simulation revealed considerable reductions of PET/mPET of up to 15.6 • C/11.6 • C during a very hot day where diurnal T amb surpassed that of 35 • C; and, (ii) the winter simulations revealed the risk of over shading within certain conditions, which in turn, led to a reduction of PET/mPET of up to 2.7 • C/2.6 • C. Such results relate directly to the significance of considering year-round implications of shading patterns within outdoor spaces, especially due to their identified dynamic and occasionally dramatic effect upon comfort conditions year round.

•
Today, in most cases, the urban fabric is already consolidated.As a result, the interdisciplinary practices of urban design and climatology must frequently seek resolutions to address thermal comfort factors in existing UCCs.Consequently, this results in the need for these two disciplines to symbiotically collaborate with one another.This originates the need for practices such public space design to (i) identify existing thermal risk factors; and just as importantly, (ii) produce creative solutions to generate better thermal environments for pedestrians in an era moreover prone to climate aggravations.The types of solutions are numerous, yet this study has focussed upon the use of urban vegetation to obtain such outcomes.

•
The study has sought to evaluate how an updated simulation tool and a newly updated thermo-physiological index can contribute to existing studies which recognise the importance of evaluating the in-situ influences of urban vegetation within different urban morphological compositions.

•
Focused upon the case of Lisbon, the study identified the in-situ influences of one of the most common shading trees in UCCs which are commonly found within the city's historical district; one which is known to particularly suffer from high T amb values and UHI effects during the summer.Notwithstanding, the study also examined the possible negative 'in-situ' thermo-physiological impacts that the Tipuana tipu species can have during the winter due to the permanency of its foliage during the colder periods of the year.As a result, the presence of vegetation in different UCCs was approached as an urban element which can both improve and decrease pedestrian thermal comfort during different periods of the year.

•
Beyond presenting bioclimatic results (which deliberates upon both singular climatic variables and thermo-physiological stress thresholds), the study suggests the importance of facilitating the transversal comprehension of the outcomes for non-climatic experts to easily interpret such outcomes discussed in the study.As an example, through the use of the CTIS software, the data from the meteorological station and the attained G rad were translated into easily interpretable figures which are constructed upon the common measuring unit of • C.Even when considering elements such as urban tourism (which is particularly elevated within Lisbon's historical centre during the summer), such a facilitation can also prove important when safeguarding other urban socio-economic aspects which are also intrinsically dependent upon a thermally comfortable public realm.

Figure 1 .
Figure 1.Comparison of monthly averages of Tamb and relative humidity (RH) field recordings obtained without the presence of a shading tree against those obtained directly beneath a vegetative crown in the central area of a low urban canyon case (UCC) in Lisbon's historical district | Source: Adapted from [54].

Figure 1 .
Figure 1.Comparison of monthly averages of T amb and relative humidity (RH) field recordings obtained without the presence of a shading tree against those obtained directly beneath a vegetative crown in the central area of a low urban canyon case (UCC) in Lisbon's historical district | Source: Adapted from [54].

Figure 2 .
Figure 2. Representation of the Tipuana tipu species | (A) Example of linear planation of species within the eastern lateral sidewalk in Rossio during the afternoon | (B) Dimensions of tree of which were later used in simulations | (C) Annual foliation and coloration of vegetative crown -adapted from Viñas, Solanich [72].

Figure 2 .
Figure 2. Representation of the Tipuana tipu species | (A) Example of linear planation of species within the eastern lateral sidewalk in Rossio during the afternoon | (B) Dimensions of tree of which were later used in simulations | (C) Annual foliation and coloration of vegetative crown-adapted from Viñas, Solanich [72].

Figure 3 .
Figure 3. North-to-south orientation (NSO) layout of the determined urban canyon cases (UCCs) and the central/lateral reference points (RPs) with specified crown spreads and vegetative coverage ratios (VCR); A-A cross-section in HW2.00, B-B cross-section in HW1.00, C-C cross-section in HW0.50, D-D cross-section in HW0.17

Figure 3 .
Figure 3. North-to-south orientation (NSO) layout of the determined urban canyon cases (UCCs) and the central/lateral reference points (RPs) with specified crown spreads and vegetative coverage ratios (VCR); A-A cross-section in HW 2.00 , B-B cross-section in HW 1.00 , C-C cross-section in HW 0.50 , D-D cross-section in HW 0.17

Figure 4 .
Figure 4. Variations of diurnal physiological stress (PS) based upon physiologically equivalent temperature (PET) and modified PET (mPET) at an hourly interval between 09:00 and 18:00 for July and December 2016 with identification of urban heat waves (UHW) and very hot days (VHD).

Figure 4 .
Figure 4. Variations of diurnal physiological stress (PS) based upon physiologically equivalent temperature (PET) and modified PET (mPET) at an hourly interval between 09:00 and 18:00 for July and December 2016 with identification of urban heat waves (UHW) and very hot days (VHD).

Figure 5 .
Figure 5. Identification of SVF and minutes cast in the sun during the summer/winter (M Sum /M Win ) between 09:00-18:00 for each RP for HW 2.00 and HW 1.00 .

Figure 7 .
Figure 7. Canyon variations of diurnal PS resultant of hourly mPET between canyon simulations with no vegetation (NV) vs. with the Tipuana tipu (TT) and those originally presented by the meteorological station for HW2.00 and HW1.00.

Figure 7 .
Figure 7. Canyon variations of diurnal PS resultant of hourly mPET between canyon simulations with no vegetation (NV) vs. with the Tipuana tipu (TT) and those originally presented by the meteorological station for HW 2.00 and HW 1.00 .

Figure 8 .
Figure 8. Canyon variations of diurnal PS resultant of hourly mPET between canyon simulations with NV vs. with the Tipuana tipu (TT) and those originally presented by the meteorological station for HW0.50 and HW0.17.

Figure 8 .
Figure 8. Canyon variations of diurnal PS resultant of hourly mPET between canyon simulations with NV vs. with the Tipuana tipu (TT) and those originally presented by the meteorological station for HW 0.50 and HW 0.17.

Figure 9 .
Figure 9.Diurnal wind speed (V1.1) variations between the RPs within all NSO canyons for the 16th of December.

Figure 9 .
Figure 9.Diurnal wind speed (V 1.1 ) variations between the RPs within all NSO canyons for the 16th of December.

Figure 10 .
Figure 10.Mean radiant temperature (Tmrt) differences during the summer at selected RPs between simulations with No Vegetation (NVSim) vs. with the Tipuana tipu (TTSim) during the 3rd of July.

Figure 10 .
Figure 10.Mean radiant temperature (T mrt ) differences during the summer at selected RPs between simulations with No Vegetation (NV Sim ) vs. with the Tipuana tipu (TT Sim ) during the 3rd of July.

Figure 10 .
Figure 10.Mean radiant temperature (Tmrt) differences during the summer at selected RPs between simulations with No Vegetation (NVSim) vs. with the Tipuana tipu (TTSim) during the 3rd of July.

Figure 11 .
Figure 11.Mean radiant temperature (Tmrt) differences during the summer at selected RPs between simulations with no vegetation (NVSim) vs. with the Tipuana tipu (TTSim) during the 16th of December.Figure 11.Mean radiant temperature (T mrt ) differences during the summer at selected RPs between simulations with no vegetation (NV Sim ) vs. with the Tipuana tipu (TT Sim ) during the 16th of December.

Figure 11 .
Figure 11.Mean radiant temperature (Tmrt) differences during the summer at selected RPs between simulations with no vegetation (NVSim) vs. with the Tipuana tipu (TTSim) during the 16th of December.Figure 11.Mean radiant temperature (T mrt ) differences during the summer at selected RPs between simulations with no vegetation (NV Sim ) vs. with the Tipuana tipu (TT Sim ) during the 16th of December.

Table 1 .
Average and diurnal (3rd of July and 16th of December) climatic data obtained from weather station (Index N • 08535; Okt: oktas).

Table 3 .
Description and categorisation of utilised urban canyon cases (UCCs) and their respective height, width, and height-to-width (HW) ratio.

Table 4 .
Stipulation of reference points (RPs) within each of the assessed urban canyon cases (UCCs) based upon single-point sky view factors (SVFs).

Table 4 .
Stipulation of reference points (RPs) within each of the assessed urban canyon cases (UCCs) based upon single-point sky view factors (SVFs).

Table 5 .
Input coordinates and description of reference points (RPs) within each assessed canyon, each with a total length of 200 m.Atmosphere 2018, 9, 12 11 of 29

Table 5 .
Input coordinates and description of reference points (RPs) within each assessed canyon, each with a total length of 200 m.

Table 6 .
Specific configurations of SkyHelios wind diagnostic tool.

Table 6 .
Specific configurations of SkyHelios wind diagnostic tool.
[1,54]Direction' 340°Accounting for the predominant wind patterns predominantly between North and Northwest within the historical quarter of the city[1,54]Vegetation offset in simulated canyons 50 m Accounting for possible interference of the wind vortices at the entrance/exit of the examined canyons
July, PS values often reached 'Extreme Heat Stress 2' within the NV Sim .In contrast, within all locations, the highest PS obtained by the TT Sim was just below the 'extreme heat stress' grade as revealed at 17:00 in RP C1 with a mPET value of 39.6 • C. When comparing the general mPET results between the NV Sim and the TT Sim , it was identified that the largest variation between PS thresholds took between 16:00-17:00.As summarised in Table7, the highest variations took place within the WEO, whereby a maximum reduction of mPET of 11.6 • C (equating to a PET of 15.6 • C) was observed at 16:00 in RP N ; • 16th of December, the obtained PS levels revealed, for the first time, noteworthy variations between the NV Sim and the TT Sim .Such was witnessed particularly within RP N and RP C2 between the hours of 12:00 and 15:00, where the TT SIM presented comparatively colder PS levels as a result of the shade cast by the tree crown.Between the two, RP N revealed the highest mPET reduction of 2.6 • C (equating to a PET of 2.7 • C) at 12:00.As result, such outcomes suggest the potential of vegetation to induce colder conditions due to obstructing radiation fluxes during the winter.

Table 7 .
Maximum reductions (K) in thermo-physiological indices between simulations with no vegetation (NV Sim ) vs. with the Tipuana tipu (TT Sim ) at each specific location within HW 0.17.