The Cooling Effect of Large-Scale Urban Parks on Surrounding Area Thermal Comfort

This empirical study investigates large urban park cooling effects on the thermal comfort of occupants in the vicinity of the main central park, located in Madrid, Spain. Data were gathered during hot summer days, using mobile observations and a questionnaire. The results showed that the cooling effect of this urban park of 125 ha area at a distance of 150 m could reduce air temperatures by an average of 0.63 ◦C and 1.28 ◦C for distances of 380 m and 665 meters from the park. Moreover, the degree of the physiological equivalent temperature (PET) index at a distance of 150 meters from the park is on average 2 ◦C PET and 2.3 ◦C PET less compared to distances of 380 m and 665 m, respectively. Considering the distance from the park, the correlation between occupant perceived thermal comfort (PTC) and PET is inverse. That is, augmenting the distance from the park increases PET, while the extent of PTC reduces accordingly. The correlation between these two factors at the nearest and furthest distances from the park is meaningful (p-value < 0.05). The results also showed that large-scale urban parks generally play a significant part in creating a cognitive state of high-perceived thermal comfort spaces for residents.


Introduction
Due to climate change and the growing urbanization, the heat in cities is rising rapidly [1,2].Today, heat has adversely affected urban life across the world [3], including urban areas of the Mediterranean climate [4,5].The increase in temperature in urban areas, especially densely populated areas has given rise to the phenomenon of urban heat islands (UHI) [6,7], which can threaten the health and comfort of citizens [8,9].The green urban spaces have been researched by numerous studies as an adaptive strategy to reduce the effect of urban heat and improve the health of citizens, by considering thermal comfort [10,11] as well as their socio-economic role [12].Various studies have been conducted on the effects of various types of green infrastructures aiming at reducing urban heat and thermal discomfort [13,14], which include different scales and forms such as small local parks [15], large urban parks [16], urban gardens [17], green roofs [18,19], green walls [20,21] and street trees [22,23].This study emphasizes the cooling effect of large urban parks.The studies conducted on large urban parks have shown that such parks have a significant effect on air temperature reduction and UHI [24][25][26], which is especially noticeable during the summer [27,28].However, it is noteworthy that the temperature drop does not only take place within the park but also surrounding areas [29].The cooling effect of urban parks, referred to as a park cooling island (PCI) [30], is known to affect the surrounding environment depending on the area of the green space and quality of green coverage [31,32].Two relevant indices to measure the cooling effect of urban parks are the cooling effect distance (CED) and cooling effect intensity (CEI) [33], which have been extensively used in various scales and climates [34][35][36].
Studies on the cooling effect of urban parks in different regions have shown that large scale urban parks with areas over 10 ha can have an average of 1 to 2 • C effect on surrounding areas that are an average of 350 meters away [33].In general, air temperature reductions in urban parks are typically up to 0.5-4 • C and may even cause up to 5-7 • C reduction [37].Despite the cooling effect of urban parks on their surroundings, such an outcome is not merely dependent on the traits of the green space [38].The morphology of the surrounding area of parks, the sky view factor (SVF), the spatial configuration of the location, and the covered area, also affect the perception of the cooling effect of urban parks and thermal comfort [37,39].
Numerous studies have been carried out on the cooling effect of urban parks.However, research pertaining to the extent of this effect on thermal comfort perception from psychological and physical perspectives is limited.This gap necessitates examining the cooling effect of urban parks on the aforementioned parameters at different intervals from the park to determine factors that may be applied in sustainable urban development.

Park Cooling Effect and Thermal Comfort Perception
A set of thermal indices was investigated in two experimental (effective temp (ET), resultant temp (RT), humid operative temp (HOP), operative temp (OP) and wind chill index (WCI)) and analytical (index of thermal stress (ITS), heat stress index (HIS), effective temp (ET*), stand.effective temp (SET*), out.stand.eff.temp (OUT_SET), predicted mean vote (PMV), perceived temp (PT) and physiological equivalent temperature (PET)) groups [40][41][42][43].The basis of these analytical indices was the energy balance (produced and wasted human energy).The main issue of thermal indices pertains to the average thermal comfort assessment and climate conditions of each area [44].The result of various studies on validating other indices shows that examples such as SET and PET have a high correlation (89%) with thermal comfort in open spaces [45].Most studies in recent years have utilized SET, PMV and PET to predict comfort levels in open spaces [46][47][48].
PET enables the comparison of the full effect of thermal conditions about the outside environment with individual experience [49,50].PET is one of the recommended indices in urban and regional urban planning around the world, used to predict thermal changes of urban or regional clusters [51].This index has shown a significant correlation to thermal comfort in various climatic conditions in open urban spaces as validation [52].One of the prominent influencing factors on PET condition is the T mrt (mean radiant temperature) climate variable.The indices mentioned above provide a single image of a set of individual and climate variables and enable the comparison of comfort conditions in various environments (due to global factors) [53].
In studies conducted to investigate the thermal comfort of urban parks, the PET index is commonly used.In a study conducted on a warm sunny day in Beijing, China (21 August, 2 pm) [54], it was shown that the effect of the Yuan Dynasty Relics green space with 102 ha area, in close approximately to the area of Retiro Park, reduces the PET by an average of 2 • C and a maximum of 15.6 • C. Another study conducted at Zhongshan Park in Shanghai city center with an area of 21.42 ha, showed that the PCI led to thermal comfort during Shanghai's winter as well as summer and had a PET of 15-29 • C [55].Additionally, a study conducted at the Central Park of Cairo with an area of approximately 26.01 ha [56] showed that parks had a significant effect in enhancing thermal comfort during summer, which, based on this effect, entailed a PET value of 22-30 • C throughout the day.
While the cooling effect of urban parks has generally been recognized as a strategy for mitigating urban heat and providing thermal comfort [57,58], the extent of perceived thermal comfort stemming from the cooling effect of urban parks is difficult to predict, necessitating investigations of people's attitudinal and bodily experiences concerning physical thermal conditions at various intervals in the surrounding areas of these parks [59,60].The degree of perceived thermal comfort by an occupant, which is caused by the cooling effect of urban parks, is dependent on factors such as individual behavioral and psychological traits, in addition to distance and cooling effect intensity [61][62][63].Individual demographic traits include age, gender, and physical characteristics such as height and weight of occupants [64,65].Psychological traits include the individual's experience of being in the environment and their thermal expectations and tolerance [66,67].Behavioral aspects include the extent of being covered (in terms of clothes) and the type of activities conducted by the individual [68].
A common method to determine the extent of thermal comfort experienced by occupants is the use of surveys [69].Essentially, by using such survey methods as open, semi-open, or multiple-choice questions [70], individuals can be asked about their thermal comfort experiences [71].In order to measure the human thermal comfort perception level, cognitive mapping can be utilized in combination with questionnaires.Employing this method and asking residents to identify places where they feel more comfortable from a thermal point of view provided a more comprehensible image of the residents' thermal comfort level for understanding [72,73].In recent years, utilizing questionnaires and asking direct questions about the perceived comfort of citizens alongside micro-climatic perceptions have yielded valuable results [74].
Many studies have been conducted on the cooling effect of urban parks, particularly their CEI and CED effects.However, less research has been carried out on the thermal comfort of this effect at different distances from the park.Additionally, few studies on the cooling effect of urban parks on thermal comfort have mainly investigated through physiological indices.As both psychological and physical perspectives together provide the complete concept of thermal comfort, it, therefore, necessitates examining the extent of thermal comfort created by urban green spaces from either perspective.
The purpose of this study is to investigate the cooling effect of the large urban park and its effect on the thermal comfort of occupants in the areas around the park.This study focuses on the perceived thermal comfort (PTC) and physiological equivalent temperature (PET) of locations, at a defined distance from the park and similar in terms of urban aspects and influential factors such as floor coverings, enclosures, street canyon and vegetation, whilst responding to the following research questions:

•
What is the extent of the large urban park cooling effect on thermal comfort based on PET? • What is the extent of the large urban park cooling effect on occupants' PTC?

•
What is the relationship between the measured PTC and the measured occupant PET?

Location and Selection of the Sites
This research was located in Madrid (40 • 25 08" N; 3 • 41 31" O), the capital of Spain with a population of 3,223,334 and population density of 5265.91 /km 2 , and a hot-summer Mediterranean climate (Csa) according to the Köppen classification [75].The average annual temperature in downtown Madrid is 19.9 • C during the day and 10.1 • C at night.The warmest month of the year is July, with an average temperature of 32.Madrid from 26 July 2015 represents the effect of Retiro Park in counteracting to the effect of UHI in the northern, eastern and southern areas of the park [78].
In order to accurately assess the cooling effect of Retiro Park on its surrounding areas, it was necessary to select places on a micro-scale with common features that are located at various distances from the park.The northern area of Retiro Park, due to its well-organized form, has more regular urban-type compared to other areas (east and south).Based on the temperature zones of the updated thermal map of Madrid's urban heat islands in 2015 [79], three street crossings (intersections) located in the Salamanca and Recoletos neighborhoods in the north of Retiro Park along the Lagasca Street, were chosen.These sites were located at 150 meters, 380 meters and 665 meters from the park at Lagasca intersection with the Conde de Aranda, Jorge Juan and Hermosilla streets, giving zones of different temperature ranges (yellow, orange and red; Figure 1).These types of selected local intersections represent most of Retiro Park's northern area intersections (Figure 2).Based on this map [77,79], each of the three intersections A, B and C had a temperature difference of about 0. Essentially, the UHI map of Madrid from 26 July 2015 represents the effect of Retiro Park in counteracting to the effect of UHI in the northern, eastern and southern areas of the park [78].
In order to accurately assess the cooling effect of Retiro Park on its surrounding areas, it was necessary to select places on a micro-scale with common features that are located at various distances from the park.The northern area of Retiro Park, due to its well-organized form, has more regular urban-type compared to other areas (east and south).Based on the temperature zones of the updated thermal map of Madrid's urban heat islands in 2015 [79], three street crossings (intersections) located in the Salamanca and Recoletos neighborhoods in the north of Retiro Park along the Lagasca Street, were chosen.These sites were located at 150 meters, 380 meters and 665 meters from the park at Lagasca intersection with the Conde de Aranda, Jorge Juan and Hermosilla streets, giving zones of different temperature ranges (yellow, orange and red; Figure 1).These types of selected local intersections represent most of Retiro Park's northern area intersections (Figure 2).Based on this map [77,79], each of the three intersections A, B and C had a temperature difference of about 0.8 °C.All three assessed intersections were selected based on a common geometric configuration.The considered geometric configuration included an enclosure, sky view factor (SVF), and street height and width (Figure 3, Table 1).The considered material structure included the ground surface coverage material, pavement material and building facade material, which were similar in all three regions; the ground surface material in all three regions was asphalt; the pavement material was gray tiles and the building material is mostly red-colored bricks and bright cement (Figure 3).The presence of green spaces in the street canyon was highly significant in terms of micro-climatic issues [80].In All three assessed intersections were selected based on a common geometric configuration.The considered geometric configuration included an enclosure, sky view factor (SVF), and street height and width (Figure 3, Table 1).The considered material structure included the ground surface coverage material, pavement material and building facade material, which were similar in all three regions; the ground surface material in all three regions was asphalt; the pavement material was gray tiles and the building material is mostly red-colored bricks and bright cement (Figure 3).The presence of green spaces in the street canyon was highly significant in terms of micro-climatic issues [80].In this regard, vegetation was also considered as well as the aforementioned parameters (Table 1), so the selected regions were also similar in terms of this parameter (Figure 3).

Microclimate Parameter Measurements
Micro-climatic measurements were conducted during six days in hot summer, on 22 June, 10 and 24 of July, 10 and 24 August and 10 September 2018.Data extraction started from 22 June and was spaced at roughly every two weeks, targeting sunny and calm days (no clouds) up to the end of summer (10 September).Ta and relative humidity (RH) were collected during the six days, were measured by mobile microclimate stations (HOBO MX2301A Temperature/RH Data Logger, manufactured by Onset Computer Corporation, MA, USA) with precision Ta; ± 0.2 °C and RH; ± 2.5% between 10:40 and 12:10 CEST (Central European Summer Time).The data was collected at oneminute intervals, and the weather data collection duration at each location was 10 minutes (Table 2).There were two priorities in selecting the intersections.Firstly, the points were selected according to Madrid's UHI map (Figure 1) where the distances were of different temperature ranges, and secondly, locations were selected that shared the greatest similarity in terms of physical and structural traits.In this vein, three sites A, B and C were selected to extract data, in accordance with the goals of this study.

Microclimate Parameter Measurements
Micro-climatic measurements were conducted during six days in hot summer, on 22 June, 10 and 24 of July, 10 and 24 August and 10 September 2018.Data extraction started from 22 June and was spaced at roughly every two weeks, targeting sunny and calm days (no clouds) up to the end of summer (10 September).Ta and relative humidity (RH) were collected during the six days, were measured by mobile microclimate stations (HOBO MX2301A Temperature/RH Data Logger, manufactured by Onset Computer Corporation, MA, USA) with precision T a ; ± 0.2 • C and RH; ± 2.5% between 10:40 and 12:10 CEST (Central European Summer Time).The data was collected at one-minute intervals, and the weather data collection duration at each location was 10 minutes (Table 2).The devices were placed 1.5 meters from the ground and covered by a radiation shield.The wind speed (WS) was determined by a Proster Digital Anemometer MS6252a placed at the same distance.The type of ground coverings, walls and types of and distance to vegetation that affect the temperature conditions were checked according to field studies; high-resolution images were taken by a Canon Eos 600D, 5184 pixels × 3456 pixels digital camera.Arial images were taken from Google Earth, 2018-2019.Fisheye images were taken at three intersections using a fisheye (Sigma 8 mm circular) lens.Additionally, the weather data of Retiro Park were collected via the AEMET (Agencia Estatal de Meteorologia) fixed station [76] located inside the park.These parameters were previously identified to show the effect of urban parks on thermal comfort [54].

The Questionnaire Survey
A random semi-structured survey was conducted on 145 pedestrians (nA: 49, nB: 45 and nC: 51) at intersections where microclimatic data were collected (workplace and residential).The interviewees comprised different gender and age groups (with the exception of children), and were active at different levels.Table 3 shows detailed information on the number and characteristics of individuals, on different days and intersections.
Data extraction from the surveys took place during the period of determining micro-climatic parameters at three intersections A, B and C (more and less than 10 minutes).Questions were divided into three sections.
The second section comprised four questions about the degree of thermal comfort perception during the interview.The questions were short and designed as five options (very high = 5 to very low = 1).Research questions concerning the level of individual thermal comfort included: (Q1) How much thermal comfort do you feel? (not too hot, not too cold); (Q2) How warm do you feel? (Q3) What is the extent of thermal comfort perceived through Retiro Park?(Q4) How much heat can you tolerate in this location?
The third section included open personal questions such as gender, height, weight, level of activity and type of clothing.Essentially, the questionnaire questions were designed so that in a short period (10 minutes), pedestrians in the neighborhood of different age and gender groups could be interviewed to extract information about their level of perceived thermal comfort (Figure 4).The total score from the questions is presented as an indicator of the perceived thermal comfort (PTC).However, the inquiry of PTC was ambiguous for occupants, so this question was addressed by asking how comfortable do you feel in terms of temperature (not too hot not too cold)?For the second inquiry, the question that was also directly asked was how much heat do you feel?Therefore, data relevant to this question was considered in reverse form to determine the perceived thermal comfort index (very high = 1 to very low = 5).It should be noted that in order to accurately derive the research data, prior to the start of summer on 22 May, 18 experimental micro-climatic data and questionnaires (six questionnaires at each intersection) were compiled (10 minutes at each intersection) in order to address the issues and queries of questionnaire gaps and to identify bio-climatic notions amidst the main questionnaire compilation stage.

Physiological Parameters
In order to analyze thermal comfort and derive PET, the relevant parameters should be calculated.The calculation of parameters for PET assessment includes clo, level of activity, SVF and T mrt .In this study, Ray Man 1.2 [81,82] software was used to derive SVF, T mrt and PET.The calculation of clo, which pertains to the extent of being covered for an individual, was based on its computational indices [83].Given that each clothing item has its own index, the clo for each person was based on the collected data.The greatest extent of clothing coverage for participants was a T-shirt with trousers (34%), shirt with trousers (25.4%) and the least was shorts with a t-shirt (22.6%), which were of 61 4).
pedestrians in the neighborhood of different age and gender groups could be interviewed to extract information about their level of perceived thermal comfort (Figure 4).The total score from the questions is presented as an indicator of the perceived thermal comfort (PTC).However, the inquiry of PTC was ambiguous for occupants, so this question was addressed by asking how comfortable do you feel in terms of temperature (not too hot not too cold)?For the second inquiry, the question that was also directly asked was how much heat do you feel?Therefore, data relevant to this question was considered in reverse form to determine the perceived thermal comfort index (very high = 1 to very low = 5).It should be noted that in order to accurately derive the research data, prior to the start of summer on 22 May, 18 experimental micro-climatic data and questionnaires (six questionnaires at each intersection) were compiled (10 minutes at each intersection) in order to address the issues and queries of questionnaire gaps and to identify bio-climatic notions amidst the main questionnaire compilation stage.

Physiological Parameters
In order to analyze thermal comfort and derive PET, the relevant parameters should be calculated.The calculation of parameters for PET assessment includes clo, level of activity, SVF and Tmrt.In this study, Ray Man 1.2 [81,82] software was used to derive SVF, Tmrt and PET.The calculation of clo, which pertains to the extent of being covered for an individual, was based on its computational indices [83].Given that each clothing item has its own index, the clo for each person was based on the collected data.The greatest extent of clothing coverage for participants was a T-shirt with trousers (34%), shirt with trousers (25.4%) and the least was shorts with a t-shirt (22.6 %), which were of 61 clo, 65 clo and 40 clo, respectively.Determination of activity levels was based on individual activity in the environment.Every activity had an indicator, walking 115 W/m2, sitting 60 W/m2, standing 70 W/m2, and fast walking 220 W/m2 [84] (Table 4).SVF calculation was conducted by importing fisheye images in the RayMan software, and 3D simulation of spatial geometry and environment vegetation in the Obstacle section of the software.The SVF index is one of the effective parameters in deriving the thermal comfort of the environment [85].
The PET index was considered a vital indicator in studying thermal comfort and has been used in numerous studies as a reliable indicator in deriving the thermal comfort of the external environment [46].T mrt and PET were derived based on extracted micro-climatic data (Table 5), and SVF data derived via the Ray Man software (Table 1).

Psychological Parameters
SPSS software was utilized, according to the statistical data obtained from the questionnaire to examine the significance of the data in three selected intersections, and also to examine the significance of the relationships between various indices with PTC [86,87].In order to examine cognitive maps more accurately, AramMMA software was used [88].In this program, all the cognitive maps were overlapped, and then according to the color spectrum, it determined which locations had the most point of reference (Figure 5).Each color represents the number of pointing a location (for example red means one time or purple means 12 times).In order to convert the qualitative data of cognitive maps into quantitative data, inside the AramMMA software, the Retiro area was rated, so that 100 percent of the score was dedicated to the maps that noted parks, while zero percent was considered for those that Energies 2019, 12, 3904 9 of 21 did not mention to the park.By using the cognitive map analyses, the importance of park cooling effect on resident psychological thermal comfort in different parts of the Retiro Park was represented.3.4.2.Psychological Parameters SPSS software was utilized, according to the statistical data obtained from the questionnaire to examine the significance of the data in three selected intersections, and also to examine the significance of the relationships between various indices with PTC [86,87].In order to examine cognitive maps more accurately, AramMMA software was used [88].In this program, all the cognitive maps were overlapped, and then according to the color spectrum, it determined which locations had the most point of reference (Figure 5).Each color represents the number of pointing a location (for example red means one time or purple means 12 times).In order to convert the qualitative data of cognitive maps into quantitative data, inside the AramMMA software, the Retiro area was rated, so that 100 percent of the score was dedicated to the maps that noted parks, while zero percent was considered for those that did not mention to the park.By using the cognitive map analyses, the importance of park cooling effect on resident psychological thermal comfort in different parts of the Retiro Park was represented.

Results
Data was gathered six days in hot summer, on 22 June, 10 and 24 July, 10 and 24 August and 10 September at three different intersections at different distances from the northern area of Retiro Park (Table 2).The selected data collection time interval was generally between 10:40 and 12:10 CEST (the exact time of each intersection is given in Table 2).This time was chosen since it was between the maximum and minimum temperatures in the park, so the temperature data extracted at points A, B and C were close to the average park temperature (mid-temperature taken by the AEMET weather station [76] on the data extraction days; Table 6).Based on Figure 6, the temperature range for points A, B and C (gray rectangle) was between the maximum and minimum extracted park temperature (linear range).It is noteworthy that among the three areas, the average temperature of point A was closer to the average temperature of the park (mid-temperature; Table 6) compared to the other points.AEMET data showed that the lowest and highest temperatures pertaining to the data extraction days in the Retiro Park were related to the 4:50-6:00 and 13:40-14:40 CEST timeframes.

Results
Data was gathered during six days in hot summer, on 22 June, 10 and 24 July, 10 and 24 August and 10 September at three different intersections at different distances from the northern area of Retiro Park (Table 2).The selected data collection time interval was generally between 10:40 and 12:10 CEST (the exact time of each intersection is given in Table 2).This time was chosen since it was between the maximum and minimum temperatures in the park, so the temperature data extracted at points A, B and C were close to the average park temperature (mid-temperature taken by the AEMET weather station [76] on the data extraction days; Table 6).Based on Figure 6, the temperature range for points A, B and C (gray rectangle) was between the maximum and minimum extracted park temperature (linear range).It is noteworthy that among the three areas, the average temperature of point A was closer to the average temperature of the park (mid-temperature; Table 6) compared to the other points.AEMET data showed that the lowest and highest temperatures pertaining to the data extraction days in the Retiro Park were related to the 4:50-6:00 and 13:40-14:40 CEST timeframes.Gray rectangles: Temperature at three intersections (bottom side: T a A, top Side: T a C).Linear range: The temperature range of the park is based on AEMET (bottom line: T a mid of park, top line: T a Max of the park. The mean T a for points A, B and C are 25.99 • C, 26.62 • C and 27.27 • C, respectively.The T a difference over the measurement days between distance of 150 m (A) and 380 m (B) from the park was about 0.63 • C, between B and C (distance of 665 m) was about 0.65 • C and for the difference between A and C, the temperature difference was about 1.28 • C. According to the results, as the distance from the park increases, the temperature and its difference with areas closer to the park will be greater.

Assessment of the Park's Cooling Effecting on Thermal Comfort Indices
In this section, PET data and questionnaire results were used to assess the cooling effect of Retiro Park on thermal comfort.PET data are derived based on the indices standard and the perceived thermal effect is obtained based on the average questionnaire and cognitive maps data.By using both these data, the extent of thermal comfort originating from the Retiro Park cooling effect at distances A, B and C could be deduced as physiological and psychological points of view.

Comparison of the Park Cooling Effect Significance on PET at Different Distances from the Park
In a more accurate analysis using SPSS software, a one-way ANOVA was used to analyze the relationship between distances A, B and C from the park and the relevant PET differences.Upon every variance analysis and in the case of significant mean difference (the significant level of a p-value less than 0.05), Tukey test, which is a series of post-hoc tests, were used to accurately determine which of the average of a variable has a significant difference [89].Essentially, by using this test, the thermal comfort relationship was deduced based on the PET index, which was obtained using the Ray Man software on distances A, B and C.
The results of the ANOVA test (Table 7) show that concerning the PET index, there was a lower significant error rate (0.05) and 95% confidence level of a significant difference between average PET data in each of the A, B and C points (p-value = 0.029).Therefore, Tukey's post-hoc tests were used to determine whether there is a significant difference between each of the A, B and C intersections.The results of this studied grouping test are presented in Table 8.According to this table, it can be stated that there was a significant relationship between the mean PET index at points A and C but the mean PET index between point B with points A and C did not entail a significant difference.The mean PET for A, B and C was 29.3 • C, 31.3 • C and 31.6 • C, respectively.Essentially, it can be stated that the PET index difference among intersections A and C was explicit and significant (p-value = 0.036).The mean of this index in A had a difference of 2.3 • C with C, while the difference between A with B was about 2 • C, which was not considered significant in terms of the Tukey test (p-value = 0.09).In order to compare the significance of the PTC at different distances from the park, the mean dataset extracted from the questionnaire for each section was derived as the PTC.The questions answered by citizens via the questionnaire were about how they felt about the ambient temperature and the effect of Retiro Park on such feelings in the form of four questions (Table 3).Based on the results of the ANOVA statistical test (Table 9) and a lower significance error level (p-value < 0.05), there was 95% confidence pertaining to a significant difference between mean thermal comfort index (PTC) in each of the A, B and C regions (p-value = 0.032).Thus, Tukey's post-hoc test was used to determine the significant difference between each of the intersections.Since the rating of questions was from very low to very high in the form of numbers 1 to 5, Table 10 shows that the average of the total questions during the 6 days for A, B and C was approximately 3.3., 2.9 and 2.8, such that area A was of the highest share indicating more thermal comfort experienced by people in this area.In interpreting the table above, it can be stated that there was a significant difference between the mean PTC index at points A and C (p-value < 0.05), but the mean PTC at point B with points A and C had no significant difference.Based on this relationship, it is also clear that with an increase in distance from the green area, the perceived comfort experienced by citizens was reduced; thus there was an inverse correlation between increased distance and PTC.The meaningfulness of the cognitive maps data with the distance from the park was evaluated using the ANOVA test.The analysis of this section was performed by using a AramMMA software for analyzing cognitive maps.In this analysis, the maps pointing to the Retiro park were distinguished from maps that did not mention the park.In fact, 100 percent of the score was dedicated to the maps noting parks, while zero percent was considered for those that did not mention to the park.The result of the ANOVA test (Table 11) illustrated that there was an acceptable agreement (95%) between the correlation of cognitive maps and the distance from the park (p-value < 0.05).Owing to the agreement, in each of the A, B and C regions the Tukey test was used separately, and the results were as follows.As illustrated (Table 12), 87.7%, 60% and 50.9% of the respondents in the A, B and C areas, respectively, referred to the park.As expected, the highest level of cognitive mapping from the park is in the vicinity of the park area (150 to 380 meters).Nevertheless, in the C region, with 665 meters distance from the park, more than half of the respondents claimed that the site was comfortable thermally thus preferring to spend more time there.The difference between A and C was also meaningful (p-value = 0.000).The good agreement of this relationship depicts the impact of the Retiro park (as a place with a high level of thermal comfort) on the resident's mental images and psychological dimensions.For a closer assessment of the effect of Retiro Park on thermal comfort, the PET index and PTC by citizens were compared on the basis of the total data extracted from the questionnaire.Mainly, this PET analysis, which is a standard index for perceived thermal comfort, was compared with the personal opinion of people about their perceptions concerning the thermal comfort of the environment.SPSS software was used to assess the correlation test as well as other statistical tests.
Based on the results of the Pearson correlation test (Table 13), the relationship between the PET index with thermal comfort in the selected points of A and C was significant (p-value < 0.05) and was of inverse correlation (−0.404 and −0.379 ) .Due to the medium correlation coefficient (medium = −0.3 to −0.5) [90][91][92], the relationship between the PET index and the PTC in area A (p-value = 0.004) was more significant that of area C (p-value = 0.006).Moreover, the significance of PET index relationship with perceived thermal comfort in area B (p-value = 0.061) was not accepted (p-value > 0.05).Regarding the results, by increasing the distance from the park at intersections B and C, despite physical and structural similarities to the intersection the A, there was an increase in air temperature (T a ), and consequently, the extent of PET and T mrt also changed accordingly.
The data pertaining to the perceived thermal comfort (PTC) acquired from the questionnaire was rated from 1 to 5 (1 = very low to 5 = very high) and at points A, B and C during six days were on average approximately 3.3, 2.9 and 2.8, respectively (Table 10).At point A, residents experienced higher thermal comfort compared to the other points.Essentially, the level of PTC at this location was medium to high according to citizens' opinions (medium = 3), whereas for B and C, the average perceived thermal comfort for the data extraction days was lower than medium.Furthermore, the results show that the extent of PTC at intersection A was more significant in all the days (Figure 7).It is noteworthy that based on statistical data, the perceived thermal comfort was correlated to PET and was of inverse correlation, such that its graph behavior was inversely correlated with the PET graph (Figures 7 and 8).For example, on 10th August, where the PET value was at a minimum at point A (PET 24.3 • C), the relevant perceived thermal comfort was at a maximum (3.6) or on 10th July, at point C, the PET was 36.3 • C, which was higher compared to other points and days which was inversed in the perceived thermal comfort graph and the least value of this index was at point C with a value of 2.4.In addition to the similar behaviors of PET, T mrt , T a and PTC parameters at each intersection during the six days of assessment (Figures 7 and 8), the credibility of the data and calculations conducted in this study were approved.Based on the results, it can be said that Retiro Park played an essential role in providing thermal comfort to the citizens during the summer days in downtown Madrid, which was physically and psychologically debatable.As mentioned, Madrid has a Mediterranean climate that is characterized by hot summers.Studies in the Mediterranean climate have shown how urban green spaces can reduce the impact of urban heat [93,94].A study carried out in Greece, Athens [95], found that the urban park on the western margin of Athens during the hot summer days could reduce daytime air temperature between 0.2 °C and 2.6 °C.In a similar study conducted in Lisbon, another Mediterranean city during the six days of summer 2007 and 2008 [96], the results showed that a 0.24 ha urban park was able to reduce the air temperature inside the park to 6.9 °C compared to its surrounding air temperatures.These studies and the studies mentioned in the first part of the article indicate that large-scale parks affect their surroundings through cooling down their environment.This is significant for Retiro Park due to its variety of vegetation and vegetation density compared to previous studies.The cooling effect of the park can cool over a distance of 350 m in its northern part, with a dense and regular texture, between 0.06 and 1.28 °C compared to the 655 m range near Heat Island.
As noted earlier, little is known about the cooling effect of urban parks on thermal comfort, and most of the studies have focused on the CED and CEI levels.However, in the Mediterranean areas, a study was carried out in Israel investigating the thermal comfort of urban parks using the PET index.Based on the results, it can be said that Retiro Park played an essential role in providing thermal comfort to the citizens during the summer days in downtown Madrid, which was physically and psychologically debatable.As mentioned, Madrid has a Mediterranean climate that is characterized by hot summers.Studies in the Mediterranean climate have shown how urban green spaces can reduce the impact of urban heat [93,94].A study carried out in Greece, Athens [95], found that the urban park on the western margin of Athens during the hot summer days could reduce daytime air temperature between 0.2 • C and 2.6 • C. In a similar study conducted in Lisbon, another Mediterranean city during the six days of summer 2007 and 2008 [96], the results showed that a 0.24 ha urban park was able to reduce the air temperature inside the park to 6.9 • C compared to its surrounding air temperatures.These studies and the studies mentioned in the first part of the article indicate that large-scale parks affect their surroundings through cooling down their environment.This is significant for Retiro Park due to its variety of vegetation and vegetation density compared to previous studies.The cooling effect of the park can cool over a distance of 350 m in its northern part, with a dense and regular texture, between 0.06 and 1.28 • C compared to the 655 m range near Heat Island.
As noted earlier, little is about the cooling effect of urban parks on thermal comfort, and most of the studies have focused on the CED and CEI levels.However, in the Mediterranean areas, a study was carried out in Israel investigating the thermal comfort of urban parks using the PET index.
In this study conducted by Cohen, Potchter and Matzarakis [97], a total of 10 urban parks of various areas (0.2 to 0.36 ha) were assessed in Tel Aviv.Although the investigated parks were smaller in terms of area compared to Retiro Park, the results illustrated that parks with richer vegetation density had greater cooling effects and thermal comfort, and could reduce temperatures by up to 3.8 • C whilst bringing PET to 18 • C during hot summer.
This study proved how important this is in areas with hot summers.However, in this study and other studies on the thermal comfort of large urban parks, less attention was paid to the thermal comfort from a psychological point of view.In this study, in addition to assessing CED and CEI of a park and thermal comfort from the physiological point of view (PET), mental thermal comfort through questionnaires (PTC) and cognitive maps were also investigated.Based on the results, the significance of both PTC and PET indices was confirmed.The two indices were inversely correlated, and when the cooling effect was reduced by the distance from the park, the PET rate increased and the PTC level decreased, which indicates the impact of the cooling effect of the large urban park on the thermal comfort from both psychological and physical perspectives.
In this regard, the role of large urban parks was tangible in providing thermal comfort for citizens and it is necessary to pay more attention to various levels of urban planning and sustainable development.

Conclusions
This study examined the potential of large urban parks in providing thermal comfort for citizens living within the park perimeter.The results extracted amidst six hot summer days in Madrid show that large urban parks exhibit a cooling effect.Considering the significance of the mean T a difference between the distance of 150 m and 665 m from the park (p-value < 0/05), as well as the lower temperature of about 1.28 • C pertaining to the distance closest to the park compared to distances further away from the park (under equal conditions), in addition to the insignificant mean T a difference at the distance of 380 m compared to 150 m and 665 m (p-value > 0/05), we found that the cooling effect of the large urban park at distances close to the park (under 380 meters) had a significant role in temperature reduction.
Accordingly, the level of PET would increase as the distance from the park increased, and residents would perceive less thermal comfort compared to distances closer to the park.The degree of PET index at a distance of 150 meters from the park was on average 2 • C PET and 2.3 • C PET less compared to distances of 380 meters and 665 meters respectively.The PTC of citizens was acquired based on the average obtained data from the questionnaire, which showed that people in the vicinity of the park experienced more thermal comfort.For the other two regions of the park, which were more distant, less thermal comfort was experienced (less than average).
The results of the cognitive maps analyses demonstrated that large parks played a significant role in thermal comfort improvement affecting people's mental map.For people, such parks are a space where they feel more comfortable, thereby spending more time to enjoy the desirable temperature.Although the resultant mental maps were closer to residents in the vicinity of parks, the results demonstrated that at long distances to the park (665 meters) those locations still had a psychological dimension in offering thermal comfort (more than 50%).
PET and PTC as the main variables of this research were inversely correlated such that when the distance from the park was increased, the PET was increased and thermal comfort was decreased.The correlation between the two indices was significant in the nearest and furthest distance from the park (p-value < 0.05) and the highest correlation coefficient of the two indices pertains to the distance closest to the park (150 meters).Essentially, this indicates the high level of perceived thermal comfort from the citizens' point of view as well as the PET compared to distant locations.

Figure 1 .
Figure 1.Madrid's urban heat island (UHI) map on 26 July 2015 [77,79] and the selected northern region of the Retiro Park.Intersection A is in the yellow zone, intersection B is in the orange zone and intersection C is in the red zone, close to one of the UHIs of Madrid.

Figure 1 .
Figure 1.Madrid's urban heat island (UHI) map on 26 July 2015 [77,79] and the selected northern region of the Retiro Park.Intersection A is in the yellow zone, intersection B is in the orange zone and intersection C is in the red zone, close to one of the UHIs of Madrid.Energies 2019, 12, x FOR PEER REVIEW 5 of 23

Figure 2 .
Figure 2. The northern region map of Retiro Park and the location of the selected sites.Point A is located 150 meters away from the park; point B is located 380 meters away from the park and point C is located 665 m away from the park.

Figure 2 .
Figure 2. The northern region map of Retiro Park and the location of the selected sites.Point A is located 150 meters away from the park; point B is located 380 meters away from the park and point C is located 665 m away from the park.

Energies 2019 , 23 Figure 3 .
Figure 3. Street view of the three intersections (A, B and C) as well as fisheye shots.

Figure 3 .
Figure 3. Street view of the three intersections (A-C) as well as fisheye shots.
clo, 65 clo and 40 clo, respectively.Determination of activity levels was based on individual activity in the environment.Every activity had an indicator, walking 115 W/m 2 , sitting 60 W/m 2 , standing 70 W/m 2 , and fast walking 220 W/m 2 [84] (Table

Figure 4 .
Figure 4.The average frequency of questionnaire data (Questions 1-4) at three intersections.Intersection A is gray, intersection B is orange and intersection C is blue.

Figure 4 .
Figure 4.The average frequency of questionnaire data (Questions 1-4) at three intersections.Intersection A is gray, intersection B is orange and intersection C is blue.

Figure 5 .
Figure 5. Cognitive map analyses using the AramMMA software [88] during six hot summer days of 2018.

Figure 5 .
Figure 5. Cognitive map analyses using the AramMMA software [88] during six hot summer days of 2018.

Table 6 .
Mean values for air temperature (T a ), in the three intersection (10:40-12:10 CET) and values for air temperature (T a ), relative humidity (RH) and wind velocity (W) in the Retiro Park on all the measurement days.

Figure 6 .Table 6 .
Figure 6.The diagram shows the temperature range of the three intersections A, B and C measured and the temperature range of Rétiro Park during the six days of hot summer 2018 in Madrid.Gray rectangles: Temperature range at three intersections (bottom side: Ta A, top Side: Ta C).Linear range: The temperature range of the park is based on AEMET (bottom line: Ta mid of park, top line: Ta Max of the park.Table 6.Mean values for air temperature (Ta), in the three intersection (10:40-12:10 CET) and values for air temperature (Ta), relative humidity (RH) and wind velocity (W) in the Retiro Park on all the measurement days.

Figure 6 .
Figure 6.The diagram shows the temperature range of the three intersections A, B and C measured and the temperature range of Rétiro Park during the six days of hot summer 2018 in Madrid.

4. 1 . 2 .
Comparison of the Significance of the Questionnaire Dataset PTC at Different Distances from the Park

4. 1 . 3 .
Comparing the Cognitive Maps Significance at Different Distances from the Park

Energies 2019 , 23 Figure 7 .
Figure 7.The average data for questionnaire data titled PTC at each three intersections A, B and C during six days.

Figure 7 .
Figure 7.The average data for questionnaire data titled PTC at each three intersections A, B and C during six days.

Figure 8 .
Figure 8.The relationship among mean PET (gray), mean Ta (blue) and Tmrt (orange), at each three intersections A, B and C during six days.

Figure 8 .
Figure 8.The relationship among mean PET (gray), mean T a (blue) and T mrt (orange), at each three intersections A, B and C during six days.
[77]during the day.Then August, June and September are the hottest months with average daily temperatures of 31.3 • C, 38.2 • C and 26.4 • C, respectively[76].Retiro Park has an area of approximately 125 ha.It is the largest and oldest park in the center of Madrid with rich plant diversity.One of the prominent issues of this study was the method of selecting distances from Retiro Park.Although studies have not been conducted on the average cooling effect of Retiro Park, an updated UHI map of Madrid was presented in 2015 by Núñez Peiró et al.[77], where the effect of Retiro Park was made evident based on the temperature color spectrum.Essentially, the UHI map of 8 • C.

Table 1 .
The geometric configuration of the intersections and street trees properties in three investigated intersections.Part Street

Table 3 .
The proportional percentages questionnaire data in the three investigated intersections on all the measured days.

Table 4 .
Clothing and activity level of responders in the three investigated intersections.

Table 4 .
Clothing and activity level of responders in the three investigated intersections.

Table 5 .
Values of the T mrt ( • C) and average values for physiological equivalent temperature (PET; • C) in part A, B and C.

Table 7 .
ANOVA PET analyses in the three investigated selected points.

Table 8 .
Tukey PET analyses the three investigated intersections (A, B and C).

Table 9 .
ANOVA perceived thermal comfort (PTC) analyses in the three investigated selected points.

Table 10 .
Tukey PTC analyses in the three investigated intersections (A, B and C). (I

Table 11 .
ANOVA maps analyses in the three investigated selected points.

Table 12 .
Tukey cognitive maps analyses in the three investigated intersections (A, B and C).

Table 13 .
Pearson correlation analyses between PET and PTC.