Influence of Green Spaces on Outdoors Thermal Comfort — Structured Experiment in a Mediterranean Climate

As a growing part of the global population lives in cities, green spaces are an essential asset for improving quality of life. This study aims to address the role of green spaces in providing favorable thermal comfort conditions for the use of outdoor spaces. The research methodology consisted of a structured experiment where a defined set of individuals from different age groups was exposed to differentiated microclimate conditions. Four nearby locations were considered, ranging from a stone-paved surface without shade to high tree canopy coverage over grass. This experiment took place in three different days in summer and early autumn conditions, with a total of 432 questionnaires. Results show a wide range of thermal sensations found during this experiment, while more favorable thermal sensations were found in shaded locations. To investigate the role of prevailing meteorological and personal conditions on thermal sensations, multinomial logistic regression analysis was applied. Results show the influence of air temperature, global radiation, wind speed, and interviewees’ gender. As meteorological variables were influenced by the diverse contexts found within a close distance inside the studied green space, results from this structured experiment suggest the need for micrometeorological diversity in the local context as a means to promote greater adaptive opportunities for green spaces users.


Introduction
About three-quarters of the European population lives in cities [1] facing various environmental stresses including air pollution, noise, and uncomfortable thermal conditions.Urban greening has been widely recognized as a key factor to mitigate the adverse effects of urbanization [2][3][4][5].However, urban greening may also be associated with multiple negative effects, such as allergies [6], reduced passive pollutant dispersion [7], and undesirable shade [8].Both the positive and negative effects of green spaces are related to the urban context and the quality of urban design.Amongst all recognizable benefits from urban greening, the improvement of thermal conditions has been addressed by different studies demonstrating the positive influence of green spaces at multiple levels, including the regional [9], urban (such as the urban heat island) [10][11][12], and microscales [13][14][15][16][17][18][19][20][21][22][23][24].
Green spaces can also provide ideal locations for social interaction and healthy living by different age groups [25][26][27] and as several studies [28][29][30][31][32] suggest, there is a strong relation between microclimatic conditions and the patterns of use of open spaces.Survey studies also show that most users tend to perceive the importance of green spaces in improving the thermal environment [33][34][35], particularly by providing shade during the warmer periods.By understanding human relations with outdoor spaces, relevant suggestions can be provided to urban designers.
Trees influence the thermal environment by providing shade, intercepting radiation with their canopy and reducing surface temperature, followed by heat transference through convection from warmer areas.Additional cooling is provided by the conversion of radiation into latent heat and evapotranspiration [36].Grass surfaces can also provide a cooler environment, especially when wet [37], although as mentioned by Doick et al. [38], grass surfaces are less effective than trees in cooling the environment.
Urban environments are often complex and offer diverse comfort conditions within short distances [39], as is the case of the differences between green spaces and its surroundings [40], in what is known as the "Park Cool Island" (PCI) effect [22,40,41].
The research trying to address outdoor thermal comfort is complex, as it is driven not only by meteorological conditions but also by the physiological and psychological attributes of individuals, which are not static, as psychological adaptation plays an essential role [42,43].Consciously or unconsciously, people tend to adapt to local circumstances by modifying their clothing and activity patterns [44][45][46][47], looking to balance the meteorological conditions by improving their personal thermal sensations.
Moreover, people tend to adapt their behavior by visiting different locations and performing diverse activities, choosing spatial attributes that adjust to their thermal and psychological preferences at a given moment [48,49].They often show preferences for spaces that include vegetation [42,50] and particularly trees [20,30].This latter option may be the consequence of the cooling effects of trees during the summer in warm climate conditions [31,51].
Several authors [46,49] have tried to establish the relationship between meteorological circumstances and human response through behavior and perception studies.When compared with indoor spaces, outdoor environments offer a particular context, with a differentiated contribution from direct solar radiation as well as higher and more variable air velocities [52].
Due to their subjective nature, perceptions are traditionally addressed by surveys (questionnaires), combined with the assessment of both meteorological site conditions and physiological attributes of the participants [53].Under this assumption, many studies take place in outdoor spaces (e.g., [29,46,48,[54][55][56]), and particularly in parks (e.g., [32,48,49,57]).The common characteristic of most of such studies is the fact that they allow users to act on their own will before interviews.These circumstances lead to variability in study conditions, including both personal and physiological parameters such as clothing and metabolism, respectively, and some psychologically-driven factors, such as perceived control, time exposure, and environmental stimulation, due to multiple adaptation processes (physiological or psychological) [49], thus increasing the chances of finding persons under comfortable conditions.
A different approach to the study of outdoor thermal comfort consists of developing structured experiments, trying to strengthen the control over both physiological and psychological adaptation factors, addressing selected locations, and designing experiments which restrict clothing choices, time exposure, and activity levels, while reducing the variation in such psychological effects as perceived control and environmental stimulation.These structured experiments have rarely been used to address outdoor thermal comfort (with the exception of the studies by Givoni et al. [58]), but are closer to the traditional approach to indoor thermal comfort studies, as is the case of several studies developed by Fanger [59] using pre-established clothing insulation (clo) and metabolic rates (met), and controlled thermal conditions.
Although some general assumptions can be made regarding the relevance of green spaces and vegetation in improving outdoor thermal comfort, there is the need for further studies regarding specific climate contexts [60].Despite the recent attempt to determine a common Mediterranean Outdoor Comfort Index by Golasi et al. [46,61] and the existence of several studies regarding personal thermal sensations in Mediterranean climate locations [48,55,56,[62][63][64], there is still the need for wider recognition of specific thermal comfort conditions, as is the case of the specific warm-summer Mediterranean climate, category Csb under the Köppen-Geiger classification [65].
The aim of this study was to evaluate the role of green spaces in providing favorable thermal comfort conditions for the use of outdoor spaces, with the development of a structured experiment in selected green spaces locations, under summer and early autumn conditions.

Study Area
The study area is located in the city of Bragança (41 • 44 N, 6 • 48 W, 690 m above a.m.s.l.) in the Northeast of Portugal.The structured experiment took place in a garden inside a University Campus (Instituto Politécnico de Bragança).
The climate of Bragança is characterized by an annual precipitation that averages 742 mm, with an irregular distribution through the year.Winter rainfall represents about 60% of the annual total.The annual average maximum and minimum temperatures are 17.6 • C and 6.7 • C, respectively.The monthly climograph of Bragança is shown in Figure 1.The highest monthly average air temperature was reached in July with 21.1 • C (average maximum temperature, 28.3 • C), and the coldest in January with 4.5 • C (average minimum temperature, 0.5 • C).The average relative humidity (RH) during the year ranges from 40% in July to 80% in January.Meteorological data shows that summers in Bragança are warm and dry and winters are cold and wet, fitting the climatic classification of Köppen-Geiger under category Csb [65].

Research Methodology
This study took place in four different locations within a close distance of no more than 50 m apart from each other, under contrasting conditions, with the objective of evaluating the influence of microclimate variables such as global radiation (St), relative humidity (RH), wind speed (V), air temperature (T a ), and mean radiant temperature (T mrt ) on thermal sensation.The locations had the following conditions (Figures 2 and 3b  A total number of 12 persons (six men and six women), grouped into four age groups (20-30, 30-40, 40-50, and 50-60 years) participated in this study, wearing jeans and a white t-shirt (approximately 0.5 clo [66,67]).Each participant filed in a questionnaire assessing individual thermal perceptions.Each age group, with three participants, stayed seated, facing the sun and resting (metabolic rate of 55 to 70 W/m 2 [67]) in each location for over 15 minutes (10 minutes adjusting to local conditions and five minutes filling in the questionnaire), moving to the next location in a rotation scheme (Figure 3a).After approximately 60 minutes, every group had gone through the four locations.This procedure took place in three periods during the day, starting in the morning at 09:00 h, again after lunch at 14:00 h, and in the mid-afternoon at 17:00 h.This experiment was repeated in three sunny days during the summer and early fall, in the months of June, July, and October, with 432 questionnaires, 36 for each of the surveyed individuals.This research combined the application of on-site questionnaires and meteorological data collection.Questionnaires were used to assess subjects' thermal perceptions taking into consideration the recommendations provided by ISO 10551 [68], including the evaluation of the personal thermal sensation (or personal thermal state) in a symmetrical seven-degree two-pole scale.Other  A total number of 12 persons (six men and six women), grouped into four age groups (20-30, 30-40, 40-50, and 50-60 years) participated in this study, wearing jeans and a white t-shirt (approximately 0.5 clo [66,67]).Each participant filed in a questionnaire assessing individual thermal perceptions.Each age group, with three participants, stayed seated, facing the sun and resting (metabolic rate of 55 to 70 W/m 2 [67]) in each location for over 15 minutes (10 minutes adjusting to local conditions and five minutes filling in the questionnaire), moving to the next location in a rotation scheme (Figure 3a).After approximately 60 minutes, every group had gone through the four locations.This procedure took place in three periods during the day, starting in the morning at 09:00 h, again after lunch at 14:00 h, and in the mid-afternoon at 17:00 h.This experiment was repeated in three sunny days during the summer and early fall, in the months of June, July, and October, with 432 questionnaires, 36 for each of the surveyed individuals.This research combined the application of on-site questionnaires and meteorological data collection.Questionnaires were used to assess subjects' thermal perceptions taking into consideration the recommendations provided by ISO 10551 [68], including the evaluation of the personal thermal sensation (or personal thermal state) in a symmetrical seven-degree two-pole scale.Other A total number of 12 persons (six men and six women), grouped into four age groups (20-30, 30-40, 40-50, and 50-60 years) participated in this study, wearing jeans and a white t-shirt (approximately 0.5 clo [66,67]).Each participant filed in a questionnaire assessing individual thermal perceptions.Each age group, with three participants, stayed seated, facing the sun and resting (metabolic rate of 55 to 70 W/m 2 [67]) in each location for over 15 min (10 min adjusting to local conditions and five minutes filling in the questionnaire), moving to the next location in a rotation scheme (Figure 3a).After approximately 60 min, every group had gone through the four locations.This procedure took place in three periods during the day, starting in the morning at 09:00 h, again after lunch at 14:00 h, and in the mid-afternoon at 17:00 h.This experiment was repeated in three sunny days during the summer and early fall, in the months of June, July, and October, with 432 questionnaires, 36 for each of the surveyed individuals.
This research combined the application of on-site questionnaires and meteorological data collection.Questionnaires were used to assess subjects' thermal perceptions taking into consideration the recommendations provided by ISO 10551 [68], including the evaluation of the personal thermal sensation (or personal thermal state) in a symmetrical seven-degree two-pole scale.Other information included respondents' perspective on existent conditions, such as thermal preference and the perception on wind, radiation and temperature.Measurements of outdoor meteorological data provided simultaneous information on local conditions.

Meteorological Measurements
The monitoring of meteorological variables took place simultaneously with the users' perceptions surveys.Multiple sensors were used, according to the specifications in Table 1, to measure meteorological variables (air temperature, relative humidity, wind speed and direction, global solar radiation, and globe temperature) near the survey locations (approximately 120 cm over the ground).Data were collected and averaged using the CR10X datalogger (Campbell Scientific, Logan, UT, USA).The sampling interval was one minute for all the meteorological variables.The mean radiant temperature (T mrt ) was calculated using the globe temperature and wind speed according to Equation [69]: where T g is the globe temperature ( • C), T a is the air temperature ( • C), V is the wind speed (m•s −1 ), ε is the globe emissivity (0.95), and D is the globe diameter (mm).The empirically-derived parameter (1.1 × 10 8 ) and the wind exponent (V 0.6 ) together represent the globe's mean convection coefficient (1.1 × 10 8 V 0.6 ).The micrometeorological conditions in the four locations were monitored using multiple instrumentation, evaluating simultaneous conditions.The exceptions were: transmitted radiation, measured under the tree shadow with a tube solarimeter (Delta-T Devices) placed over 0.2 m above ground; and wind speed and direction, measured only in two locations: behind a windscreen upstream from a wind barrier (TENT) and in a exposed wind location (SUN) representative of the prevailing wind on all other locations.Data were analyzed using SPSS statistics 18 [70], including both generic graphic interpretation and multinomial logistic regression analysis.Multinomial logistic regression was applied because it allows for the evaluation of the statistical relation between a dependent categorical variable and multiple independent continuous variables (covariates) or categorical factors (design variables) [71].SPSS uses the maximum likelihood method to calculate the logistic coefficients for this analysis.

Meteorological Variables in Different Experimental Locations
The different locations had contrasting meteorological conditions during the study: similar distributions in air temperature (Ta), with higher variations in both TENT and WINDSHELT locations (Figure 4a).Relative humidity (RH) distribution had slightly higher values in the TREE location, despite the similarities among the different locations (Figure 4b).Major differences were found regarding global radiation (St), and shaded locations (TREE and TENT), as expected, had much lower values when compared with sun-exposed locations (WINDSHELT and SUN) (Figure 4c).Mean radiant temperature results show an increase in the intensity of this variable when moving from more vegetated locations (TREE) to sunny paved locations (SUN) (Figure 4d).Finally, wind speed values were lower behind a wind shelter (WINDSHELT) when compared with the reference measurement for other locations (Figure 4e).vegetated locations (TREE) to sunny paved locations (SUN) (Figure 4d).Finally, wind speed values were lower behind a wind shelter (WINDSHELT) when compared with the reference measurement for other locations (Figure 4e).

Thermal Comfort on Different Locations
As the participants in the survey were asked to express their personal thermal sensations in a scale ranging from cold (−3) to hot (+3), results show a wide range (Figure 5a), with a large number of answers in the warm side of the scale when compared with cooler personal thermal sensations.Regarding gender, women presented a higher number of responses in the lower personal thermal sensations when compared to men (Figure 5b).When looking at the distribution within the four study locations (Figure 5c), results show a higher number of answers indicating a neutral (0) sensation or slight discomfort (−1, 1) in places under shade, while sun-exposed locations had a higher number of answers at the warm end of the personal thermal sensations scale.

Thermal Comfort on Different Locations
As the participants in the survey were asked to express their personal thermal sensations in a scale ranging from cold (−3) to hot (+3), results show a wide range (Figure 5a), with a large number of answers in the warm side of the scale when compared with cooler personal thermal sensations.Regarding gender, women presented a higher number of responses in the lower personal thermal sensations when compared to men (Figure 5b).When looking at the distribution within the four study locations (Figure 5c), results show a higher number of answers indicating a neutral (0) sensation or slight discomfort (−1, 1) in places under shade, while sun-exposed locations had a higher number of answers at the warm end of the personal thermal sensations scale.

Influence of Microclimatic and Personal Variables on Personal Thermal Sensations
To address the relation between personal thermal sensations and the meteorological conditions on the different locations, a multinomial logistic regression is applied.The data shows that only females felt Cold (−3), to avoid a quasi-complete separation bias, the corresponding two records were removed from the dataset.Table 2 presents the likelihood ratio test for the best results considering the variables and factors studied in this research (model fitting sig.<0.001; pseudo R-squared measures: Cox and Snell= 0.706; Nagelkerke = 0.734; McFadden = 0.374).These results show an adjustment (p <0.001) for some of the meteorological variables (covariates), namely air temperature (Ta), wind (V), and global radiation (St).From the personal factors tested, including age groups, only gender (GEN), considering a binary condition (0-male and 1-female), had an adjustment in the multinomial logistic regression model.3 shows the coefficients (B) for independent variables in the logistic regression equations along with other statistical parameters, considering a neutral thermal sensation state (0) as the reference level for the analysis.Results from the significance test (sig.<0.05) show that, with the exception of the cold condition (−3), for each personal thermal sensation there are at least two variables that can explain the differences between such condition and the reference thermal sensation (Neutral).The odds ratio (Exp(B)) is the factor by which the odds of being in the predicted level of the dependent variable are multiplied when the independent variable increases one unit.No effect is expected when the odds ratio equal 1.0, below 1.0 indicates a negative effect, while above 1.0 indicates a positive effect [71].

Influence of Microclimatic and Personal Variables on Personal Thermal Sensations
To address the relation between personal thermal sensations and the meteorological conditions on the different locations, a multinomial logistic regression is applied.The data shows that only females felt Cold (−3), to avoid a quasi-complete separation bias, the corresponding two records were removed from the dataset.Table 2 presents the likelihood ratio test for the best results considering the variables and factors studied in this research (model fitting sig.<0.001; pseudo R-squared measures: Cox and Snell= 0.706; Nagelkerke = 0.734; McFadden = 0.374).These results show an adjustment (p <0.001) for some of the meteorological variables (covariates), namely air temperature (Ta), wind (V), and global radiation (St).From the personal factors tested, including age groups, only gender (GEN), considering a binary condition (0-male and 1-female), had an adjustment in the multinomial logistic regression model.Table 3 shows the coefficients (B) for independent variables in the logistic regression equations along with other statistical parameters, considering a neutral thermal sensation state (0) as the reference level for the analysis.Results from the significance test (sig.<0.05) show that, with the exception of the cold condition (−3), for each personal thermal sensation there are at least two variables that can explain the differences between such condition and the reference thermal sensation (Neutral).The odds ratio (Exp(B)) is the factor by which the odds of being in the predicted level of the dependent variable are multiplied when the independent variable increases one unit.No effect is expected when the odds ratio equal 1.0, below 1.0 indicates a negative effect, while above 1.0 indicates a positive effect [71].
-When in a hot condition (+3), there is a similar pattern, as the increase by one unit in temperature (1.408) and radiation (1.016) means an increase by 40.8% and 1.6%, respectively, in the odds of moving from a neutral thermal sensation (0) and into this condition (+3).Conversely, the increase by one unit in wind (0.066) means a decrease by 93.4%, in the odds of moving from a condition of neutrality (0) towards the hot condition (+3).Regarding gender, results show that there is a strong relation between gender and negative thermal sensations.Thus, for females, the odds of feeling cool (−2) are 9.53 times higher than the odds for male to have that same thermal sensation, under the same conditions.A similar situation was found for the slightly cool condition (−1), as the odds of being in this condition are 5.9 times larger for females than the odds for males to have that same thermal sensation.

Discussion
Results show that under the conditions prevailing during the structured experiment there is a consistent relationship between different meteorological conditions and personal thermal sensations, which in turn correspond to the effects generated by local attributes of the case study locations.Such effects have already been studied in many studies, such as the general park cooling effect provided by tree shade, including reduced temperature [20,22,72,73], radiation [19,31,37,50], and wind speed [74][75][76].Fewer studies address strategies to improve negative thermal sensations (−3 to −1) often felt during the cooler seasons (such as promoting solar radiation intake) [32,56] to promote a temperature increase or wind speed reduction [56], providing similar interpretations to the results presented in this study.
Differences in personal thermal sensation between males and females have already been identified in multiple studies as described by Wang et al. [77], suggesting that females are more sensible to cooler environments, thus tending to express lower (or cooler) thermal sensations.
Results from this structured experiment came to suggest the need for micrometeorological diversity in the local context, as a means to promote greater adaptive opportunities for green spaces users.Such an interpretation has already been presented by Hirashima et al. [49], who reached this conclusion looking at different behaviors from user of diverse green areas in Belo Horizonte, Brazil, during summer and winter conditions.Tseliu et al. [56] also addressed the need for contrasting conditions taking into consideration the differences between summer and winter conditions in Athens (Greece), identifying different needs for thermal adaptation in those opposing seasons.
As this study could not cover all factors influencing microclimate conditions, additional aspects such as individual adaptation to seasonal conditions and the effects of diverse vegetative treatment [56], as well as different levels of tree shade intensity [31], should be considered in future studies, as they may influence the results in this kind of research.

Conclusions
In this study, questionnaire surveys were conducted aiming at addressing the influence of local micrometeorological conditions on thermal comfort, looking at the potential effect of small distance changes inside green spaces and trying to understand their potential for ameliorating thermal conditions during warm sunny days.
Under this research conditions, results on personal thermal sensations show that there were more frequent neutral (0) to slight discomfort (−1, 1) conditions in shaded areas, particularly over grass and under tree cover when compared with sun-exposed locations.Conversely, respondents in sun-exposed locations presented higher levels of dissatisfaction concerning the warmer conditions (+1 to +3).However, during the prevailing meteorological conditions, all locations offered opportunities for thermal comfort or slight dissatisfaction throughout the experiment, thus showing the potential for the study locations to provide comfort while users perform metabolic activities (met) and wear light clothing (clo).
The use of multinomial logistic regression demonstrated the influence of microclimate conditions on personal thermal sensations.Changes in local conditions such a global radiation and wind speed, which could be found within close distance in this study, had an impact on personal thermal sensations.Results suggest the relevance of conditions provided by green spaces (tree shadow and wind penetration) on thermal comfort, but also the need for diversity in green space design, as unshaded areas offer opportunities for increasing solar radiation intake, improving the chances for thermal comfort under cooler temperature conditions, while the shaded areas provide the opposite effect in warmer temperature conditions.
The results from this study only have local validity and may not be applicable to other climate conditions and cultural contexts.This study may help to inform urban and green space designers on the importance of contrasting microclimate conditions within public open spaces, helping to inform the development of adequate design and maintenance practices that may improve thermal comfort and promote the use of outdoor spaces.

Figure 1 .
Figure 1.Monthly climograph of Bragança city.The gray band is the minimum-maximum temperature range ( • C); the vertical bars represent the monthly rainfall (mm).
): (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under direct solar radiation (SUN); (C) over grass under the shade of a tree (Cupressus lusitanica) (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Figure 2 .
Figure 2. Location of the city of Bragança (Portugal) and of the structured experiment locations: (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under direct solar radiation (SUN); (C) over grass under the shade of a tree (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Figure 3 .
Figure 3. (a) Schematic representation of the structured experiment; (b) Study locations: (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under solar direct radiation (SUN); (C) over grass under the shade of a tree (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Figure 2 .
Figure 2. Location of the city of Bragança (Portugal) and of the structured experiment locations: (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under direct solar radiation (SUN); (C) over grass under the shade of a tree (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Figure 2 .
Figure 2. Location of the city of Bragança (Portugal) and of the structured experiment locations: (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under direct solar radiation (SUN); (C) over grass under the shade of a tree (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Figure 3 .
Figure 3. (a) Schematic representation of the structured experiment; (b) Study locations: (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under solar direct radiation (SUN); (C) over grass under the shade of a tree (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Figure 3 .
Figure 3. (a) Schematic representation of the structured experiment; (b) Study locations: (A) over grass, under solar direct radiation and near a wind shelter (WINDSHELT); (B) over stone-paved ground and under solar direct radiation (SUN); (C) over grass under the shade of a tree (TREE); and (D) over grass under artificial shade provided by a tent (TENT).

Table 1 .
Measured meteorological variables and instruments.

Table 2 .
Influence of each variable on thermal sensation (likelihood ratio tests).GEN: gender; Ta: air temperature; V: wind; St: global radiation.

Table 2 .
Influence of each variable on thermal sensation (likelihood ratio tests).GEN: gender; Ta: air temperature; V: wind; St: global radiation.

Table 4 .
Exponent of the correlation coefficient for the variables with influence on personal thermal sensations, expressed in odd ratios.