Evaluating the Cooling Potential of Urban Green Spaces to Tackle Urban Climate Change in Lisbon

: The increase and optimization of urban vegetation has been considered an e ﬀ ective mitigation measure of an urban heat island (UHI), with positive e ﬀ ects on human thermal comfort. In this study, the cooling potential of all green spaces in Lisbon was estimated. For that, several mobile measurements of air temperature data were made in a single park (Gulbenkian’s Garden). These measurements were used for the interpolation of air temperature. Furthermore, urban biomass was estimated using remote sensing products, namely Landsat satellite images. Ultimately, a linear regression model was built from the relation between vegetation density and air temperature. Results regarding the estimation of biomass (AGB) in the city of Lisbon were higher in winter than in summer. The urban green spaces cooling potential model showed that for every decrease of 1 ◦ C in air temperature between a measuring point and a reference station we need to increase the area covered by vegetation by 50 m 2 (planar measure). This methodology can be applied in other urban areas for the quantiﬁcation of the cooling e ﬀ ect provided by vegetation in order to improve urban climate thermal conditions and human well-being and, consequently, to mitigate some consequences of future climate change.


Introduction
The fast process of urbanization has caused drastic changes in soil use and, especially in urban areas. The most pronounced modification is the increase of temperatures within cities, comparatively to the surrounding countryside. One of the most effective mitigation measures is the increase and optimization of urban green spaces. In fact, vegetation has been considered a key element for the quality of life of urban population, especially considering climate change scenarios induced by anthropogenic activities, with several consequences in health and human thermal comfort.
Green spaces provide other ecosystem services in addition to urban climate improvements on bioclimatic comfort: maintenance and conservation of biodiversity, air quality (absorption of atmospheric pollutants and particles; reduction of CO 2 levels; release of oxygen; fixing of dust) and water quality, protection from strong winds, reduction of noise levels, hydrological benefits (decrease in surface runoff and in the frequency of floods through interception of precipitation and increase of soil water retention), reduction in risk of soil erosion and sociocultural (recreation areas, aesthetic appreciation of landscape) and economic benefits (for example, valuation of real estate) [1][2][3][4][5]. However, in this study, the focus is on urban vegetation thermal benefits, namely the reduction of air and surface temperature through two physical processes: on one side, the shadow effect that consists in the interception of solar radiation, reflecting it in larger quantities and reducing the absorption and storage of energy by urban surfaces and, consequently, the atmospheric heating [5][6][7][8]. On the other In view of these limitations, the definition of specific recommendations for planning, management and optimization of urban green spaces is an obvious requirement to mitigate UHI effects.
Therefore, the purpose of this study is to quantify the cooling effect of urban green spaces in the urban canopy layer (UCL) in Lisbon and to estimate the amount of AGB needed to reduce 1 • C in air temperature and to fight climate change in the near future and reduce urban vulnerability.

Study Area
The city of Lisbon is located on the west coast of Portugal, about 30 km east of Atlantic Ocean, on the right bank of Tagus river (Figure 1a). In view of these limitations, the definition of specific recommendations for planning, management and optimization of urban green spaces is an obvious requirement to mitigate UHI effects.
Therefore, the purpose of this study is to quantify the cooling effect of urban green spaces in the urban canopy layer (UCL) in Lisbon and to estimate the amount of AGB needed to reduce 1 °C in air temperature and to fight climate change in the near future and reduce urban vulnerability.

Study Area
The city of Lisbon is located on the west coast of Portugal, about 30 km east of Atlantic Ocean, on the right bank of Tagus river (Figure 1a). Regarding the topographic features of the region, maximum altitudes range between 160 m, on the urbanized area, and 226 m, on the Monsanto hill (west part of the city). This hill is a vast forest park (990ha) created in 1934 with the primary function of leisure of the Lisbon population and constitutes a major obstacle to the penetration of maritime air [28].
The city has a temperate climate, with a moderate winter and the hot season is the driest period of the year, since precipitation concentrates between October and April. According to Alcoforado and Andrade [2,35], this type of climate (Köppen Csa) is modified by the rough relief and the proximity to the Tagus estuary and the Atlantic Ocean. These features, along with the latitudinal position, give the city a certain thermal amenity.
Currently, the available space for new vegetated areas is very limited in Lisbon, even though this scenario was very different in the past decades. The abundance and spatial distribution of green spaces has suffered notable alterations until now and it is far from being homogeneous throughout the county. According to Franco & McDonald [36], this city has over 120 local parks and gardens plus 13 urban forests with denser vegetation, from where stands out the Monsanto hill, which cover about 10% of the city's total area. In 2006, public urban green spaces in Lisbon occupied 1303.6 ha, whereas the area covered by tree canopy and dense shrub spots corresponded to 1558 ha [37]. Mendes et al. in 2015 [38] estimated that 16.8% of the city is covered by trees.
Focusing only on trees, there are over 600 000 specimens in Lisbon belonging to about 200 different species, indigenous and, especially, exotic [37,39]. According to Soares [3], in 2011 there were about 41,000 trees in public spaces (streets and avenues), having increased in number since then. Furthermore, of the 1500 km of roads in the city, 500 km (33%) have trees. Soares [3] prepared a record of the main street trees in Lisbon and concluded that the majority are deciduous, along with 42% of trees in public gardens.
Gulbenkian's garden (Figure 1b), the green space selected as sample in this study, has an area of 8.5 ha and is composed by several buildings (2.5 ha), while vegetated area covers the rest of the area [40]. Regarding vegetation, 34.4% corresponds to shrubs, 30.8% to trees, 29.3% to herbaceous species and only 5.5% are climbing plants. Perennial species account for 74% of total vegetation, while only 25.3% corresponds to deciduous plants.

Acquisition and Treatment of Climatic Data
In order to obtain climatic data for the analysis of the thermal behavior of urban green spaces and to interpolate air temperatures, several mobile measurements were made within and around Gulbenkian's garden ( Figure 1b). Four measurement routes (north, south, west and east) were prepared outside the green space and only one on the inside, covering distinct thermal conditions (shading vs in the sun; points with very dense vegetation vs open and wide spaces covered only with lawns; points near water bodies). All these routes took less than an hour to complete, so that the data recorded derive from different types of urban morphology, urban materials, street orientation topographic position and, mostly, the presence or absence of green spaces, and not from normal daily variation of meteorological variables [41]. A thermocouple (type k) was used to measure air temperature. All the measurements were taken, as much as possible, on days with clear sky (cloud cover equal or less than 4/8), since these conditions promote great thermal differences between urban green spaces and the surrounding built area. The data correction procedure was based on the estimation of the difference between each air temperature and relative humidity value registered in every measuring point and the same data recorded, at the same time, on a reference meteorological station (IGOT-University of Lisbon- Figure 1a).
The obtained results correspond to air temperature and relative humidity local anomalies (∆) and were used for the interpolation of Gulbenkian's garden thermal behavior ( Figure 2 and Table 1).
Temperature Anomalies (ºC)  A day with the least cloud cover was selected (20 May 2018) and air temperatures registered inside and in the surrounding area were spatialized in a GIS (Geographic Information System) environment, using ArcGIS software (v. 10.5.1): a shapefile with several points referring to all measuring locations was generated and air temperature data already corrected were introduced. For the interpolation the geostatistical method of kriging was chosen because it is based on the principle that points relatively close in space have a similar behavior than points that are away from each other and, for that, less correlated in space.  A day with the least cloud cover was selected (20 May 2018) and air temperatures registered inside and in the surrounding area were spatialized in a GIS (Geographic Information System) environment, using ArcGIS software (v. 10.5.1): a shapefile with several points referring to all measuring locations was generated and air temperature data already corrected were introduced. For the interpolation the geostatistical method of kriging was chosen because it is based on the principle that points relatively close in space have a similar behavior than points that are away from each other and, for that, less correlated in space.

Estimation of Urban Biomass in Lisbon
For the determination of AGB, data from remote sensing products (Landsat 8 images) were used. This satellite was launched in February of 2013, in a polar and sun synchronous orbit at an altitude of about 705 km. This platform is composed by two sensors, OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor). Concerning the reasons of its choice, Landsat 8 images have been widely used for the determination of biomass at local and regional scales and this satellite has the most complete data series, especially available to the public [22,[42][43][44][45]. Taken this into account, two Lansdat 8 images with the least cloud cover possible over Lisbon, representing summer (17 July 2017) and winter (5 February 2016), were chosen and the determination of green volume consisted in an extensively search of models based in remote sensing products that were already tested in areas with some climatic proximity to this city, with relatively simple calculations and using vegetation indexes like NDVI, previously tested and validated for this type of applications. Only three models were selected since they were the most feasible, simple and quick to test: Filella et al. [46], Pereira et al. [47] and Chang & Shoshany, 2016 [48]. [47] proposed a model for the estimation of biomass in Aire and Candeeiros Hill Natural Park (Portugal), using satellite images from Landsat 5. On the other side, [46] applied an empirical model to bushes on Northeast Spain (typical Mediterranean climate). Lastly, [48] presented a model for the estimation of Mediterranean shrub biomass from the fusion of Sentinel-1 and Sentinel-2 images. The equations presented in each study are synthetized in Table 2 and the results of their application were, then, used in a simple linear regression model to estimate the urban green spaces cooling potential of Lisbon. For the application of these models, NDVI (Equation (1)) was calculated from the following formula: where, NIR corresponds to near-infrared band and; Red represents the red band.
The NDVI values below 0.15 were excluded in both summer and winter situations since according to [49] this threshold is the value representative of green areas in Lisbon. The two outputs of NDVI are presented in Figure 3. For the application of all three biomass models, it should be noted that the measure units differ, reason why all results were converted to a common measure unit (kg/m 2 ).  (Figure 3b). In both maps, Monsanto Hill stands out as the "lung of the city". However, in a more detailed analysis, some of the areas with higher NDVI values (> 0,6) in winter, for example on the airport surroundings, correspond to lawns that aren't irrigated and, for that, present lower values in summer season (on most cases lower than 0.41). Furthermore, other areas (especially streets, for  Figure 3b). In both maps, Monsanto Hill stands out as the "lung of the city". However, in a more detailed analysis, some of the areas with higher NDVI values (>0.6) in winter, for example on the airport surroundings, correspond to lawns that aren't irrigated and, for that, present lower values in summer season (on most cases lower than 0.41). Furthermore, other areas (especially streets, for example Avenida da Liberdade) covered by trees present lower NDVI values in winter (less than 0.41) and slightly higher in summer.

Estimation of Cooling Potential of Green Spaces
In order to evaluate the influence of vegetation on the behavior of air temperatures, several attempts to build a simple linear regression model were made: • On a first attempt, the correlation between AGB and the temperature differences at Gulbenkian's Garden was tested and turned out not statistically significant for both seasons (R = 0.6).

•
On a second and final attempt, the linear regression model incorporated a map with the density of vegetation produced from green volume estimations of the city instead of AGB. This parameter was calculated for both winter and summer ( Figure 4) using the Kernel density, which estimates the magnitude of green mass per area.
Two different linear regression models were generated (condition of vegetation in summer and in winter), having previously been verified the correlations between the variables. example Avenida da Liberdade) covered by trees present lower NDVI values in winter (less than 0.41) and slightly higher in summer.

Estimation of Cooling Potential of Green Spaces
In order to evaluate the influence of vegetation on the behavior of air temperatures, several attempts to build a simple linear regression model were made: • On a first attempt, the correlation between AGB and the temperature differences at Gulbenkian's Garden was tested and turned out not statistically significant for both seasons (R=0,6).

•
On a second and final attempt, the linear regression model incorporated a map with the density of vegetation produced from green volume estimations of the city instead of AGB. This parameter was calculated for both winter and summer (Figure 4) using the Kernel density, which estimates the magnitude of green mass per area.
Two different linear regression models were generated (condition of vegetation in summer and in winter), having previously been verified the correlations between the variables. .

Biomass in Lisbon
The amount of AGB estimated in Lisbon is higher in winter than in summer ( Figures 5 and 6), as for the NDVI.

Biomass in Lisbon
The amount of AGB estimated in Lisbon is higher in winter than in summer ( Figures 5 and 6), as for the NDVI.  For a better interpretation of the obtained results, Table 3 presents some of the characteristics of AGB for the whole city, green spaces and street trees (source: City Hall of Lisbon).  For a better interpretation of the obtained results, Table 3 presents some of the characteristics of AGB for the whole city, green spaces and street trees (source: City Hall of Lisbon).  For a better interpretation of the obtained results, Table 3 presents some of the characteristics of AGB for the whole city, green spaces and street trees (source: City Hall of Lisbon). It is possible to verify that the average biomass for the whole city and the average biomass of green spaces are higher in winter than in summer in all models tested (positive differences reach 0.3 kg/m 2 using the model as described in [48]). However, the total green volume estimated is considerably higher in summer than in winter (positive differences reach 16.5 ton using the model described in [47]), even though the summer season in Lisbon is usually dry (Csa) and only foliage and irrigated green areas are more active. Furthermore, if one considers only the biomass of street trees in the city, the results are similar for all models in both winter and summer.

The Cooling Potential of Green Spaces in Lisbon
The correlation coefficients referring to the cooling potential model in Gulbenkian's Garden for both summer (−0.74) and winter (−0.83) are negative and strong, especially during winter (time of the year when the amount of vegetation in the garden is slightly lesser than in summer). However, since air temperatures are higher in summer and UHI is more pronounced in Lisbon in this season [50], with negative effects of comfort and human health, only the regression model for this season (Figure 7 and Table 4) is presented in this paper. It is possible to verify that the average biomass for the whole city and the average biomass of green spaces are higher in winter than in summer in all models tested (positive differences reach 0.3 kg/m 2 using the model as described in [48]). However, the total green volume estimated is considerably higher in summer than in winter (positive differences reach 16.5 ton using the model described in [47]), even though the summer season in Lisbon is usually dry (Csa) and only foliage and irrigated green areas are more active. Furthermore, if one considers only the biomass of street trees in the city, the results are similar for all models in both winter and summer.

The Cooling Potential of Green Spaces in Lisbon
The correlation coefficients referring to the cooling potential model in Gulbenkian's Garden for both summer (−0.74) and winter (−0.83) are negative and strong, especially during winter (time of the year when the amount of vegetation in the garden is slightly lesser than in summer). However, since air temperatures are higher in summer and UHI is more pronounced in Lisbon in this season [50], with negative effects of comfort and human health, only the regression model for this season ( Figure  7 and Table 4) is presented in this paper.  Density of vegetation nuclei (m 2 ) Air temperature anomalies (ºC) Figure 7. Air temperature anomalies vs density of vegetation nuclei (summer) in Gulbenkian's Garden. Note: thermal anomalies correspond to temperature differences between each value measured and reference station (IGOT). As one might see, the distribution of variables is close to the line in summer, especially in higher temperature differences and, consequently, lower densities of vegetation. With the slope of the line is possible to calculate the amount of vegetation (in the plane and not in volume) required for every increase or decrease of 1 • C in air temperature. Since the slope is negative, for every decrease of 0.02 • C in air temperature we need an increase of 1 m 2 of vegetation, in both models. Thus, the greater the density of vegetation, the lower the air temperature. Besides, the correlation (R) between air temperature and vegetation is negative and very strong, that is, 83% of variation of air temperature anomalies in winter and 74% in summer is due to vegetation. Still, the determination coefficient (R 2 ) shows that more than half (about 70% in winter and 55% in summer) of air temperature variation inside and in the surroundings of Gulbenkian's Garden is explained by the presence of vegetation.
Using this equation, an extrapolation of the cooling potential of green spaces was made for the whole city (Figure 8a). The distribution of areas with high and low potential has some similarities with the behavior of UHI in this city (Figure 8b). In fact, the areas with the fewest amount of vegetation and, for that, the lowest cooling potential, mainly the riverside front, correspond to the areas with the highest intensities of UHI (> 1.5 °C) and vice-versa. Only the city airport (NE) presents a different behavior because in summer the surrounding area is covered by dry lawns, so vegetation cooling potential is low but UHI is also weak because of its location on the city fringe.

Discussion
In general, all models of biomass in Lisbon estimated higher values in winter than in summer. However, according to Soares [3], street trees in the city are mostly deciduous, along with about 42% The distribution of areas with high and low potential has some similarities with the behavior of UHI in this city (Figure 8b). In fact, the areas with the fewest amount of vegetation and, for that, the lowest cooling potential, mainly the riverside front, correspond to the areas with the highest intensities of UHI (>1.5 • C) and vice-versa. Only the city airport (NE) presents a different behavior because in summer the surrounding area is covered by dry lawns, so vegetation cooling potential is low but UHI is also weak because of its location on the city fringe.

Discussion
In general, all models of biomass in Lisbon estimated higher values in winter than in summer. However, according to Soares [3], street trees in the city are mostly deciduous, along with about 42% of trees in public gardens. In order to better understand the reasons that explain this result, it's important to analyze the characteristics of each biomass model tested and the vegetation characteristics of the city. First, all authors used field samples, either to generate the model or to calibrate it. Regarding the satellite images, only Filella et al. [46] didn't use this type of data, while Pereira et al. [47] used Landsat 5 images and, lastly, Chang & Shoshany [48] used two different types of satellite images, Santinel-1 and Sentinel-2, generating a final fusion model. Regarding the area features where these investigations were made, all studies were applied in areas with Mediterranean climates, one of them tested in Portugal.
Concerning the models tested outside the country, on one side, Filella et al. [46], whose model generated the least amount of biomass for Lisbon in both seasons, estimated biomass of Mediterranean forest on Garraf's Natural Park in Barcelona. More recently, Chang and Shoshany [48] performed the mapping of Mediterranean forests in climatic transition zones that represent several Mediterranean environments. In these areas, the annual total precipitation varies between 200 and 600 mm, being, then, in most cases, lesser than annual average precipitation recorded in Lisbon. Nevertheless, they generated a robust model with a correlation coefficient of 0.86.
Concerning the studies made in Portugal, Pereira et al. [47] estimated the amount of biomass in Aire and Candeeiros' Hill Natural Park, with an area of 37,000 ha, altitudes ranging between 50 and 680 m, annual air temperatures of 14.8 • C and precipitation (846 mm) higher than the annual average in Lisbon. They concluded that the best green mass estimation (R = 0.76) can be obtained using NDVI.
Furthermore, all three models were applied to the hole vegetation in this city. By analyzing the values of total green volume and biomass of green spaces, the unexpectedly higher amounts in winter and lower in summer are influenced by extensive areas of smaller species/lawn and some of them aren't irrigated frequently. For that, they drier in summer (biomass in lower in this season) than in winter. Still, the estimation of biomass was based only on one case/satellite image per season. In this case, the satellite image used for winter is from the end of the season/ beginning of spring, so the amount of biomass is relatively higher than in peak winter (December or January). At the same time, the green volume of only street trees is theoretically similar in summer and in winter, even though in this latter season fall foliage occurs, so it's expected a fewer amount of biomass. For these reasons, it is important to estimate the volume of vegetation per stratum or specie.
In order to establish a parallel with other investigations that have estimated biomass in urban areas, Raciti et al. [30] produced a high-resolution map of tree biomass in Boston (EUA). During the estimation process, the authors excluded NDVI values inferior to 0.1 considering that above this threshold the amount of vegetation is sufficiently significant to be included in the analysis. They estimated an average biomass of 1.3 kg/m 2 for the city and this result is similar to the average of Filella et al.;  models for both seasons. According to the same authors, the mapping of urban vegetation might give relevant information to urban planning, especially concerning UHI mitigation and thermal comfort improvement.
Likewise, Goh et al. [31] presented a method for the estimation of biomass in large metropolitan areas from SPOT 5 satellite images. For that, they used as study area all green spaces managed in Singapura, corresponding to more than 42,000 ha, having concluded that total biomass of these areas (excluding trees on streets and natural reserve areas) is 4 kg/m 2 -much higher than average biomass estimated for Lisbon in every model (summer and winter). Besides this, studies that have estimated total biomass of an urban area are scarce, mainly not making distinction between species.
The indicator of vegetation (density of vegetation nuclei) introduced in the model of estimation of green spaces cooling potential proved to be useful and productive. In fact, many of the studies that sought to estimate the influence of vegetation on air temperature reduction used NDVI to represent it. For example, Feyisa et al. [34] studied 21 green spaces in Addis Abba (Ethiopia) and sought to evaluate the relation between vegetation characteristics and air temperature recorded, having found a positive relation between PCI and NDVI. Therefore, they concluded that the cooling effect is determined mainly by the specie type, tree cover, dimension and shape of green spaces. In the case of generated models, a major part of air temperature variation inside and in the surroundings of Gulbenkian's Garden is due to the amount of vegetation represented by the Kernel density. In Lisbon, it is clear that the larger the green volume, the lower air temperature and the bigger the Park Cooling Effect. For these reasons, the choice of vegetation type to integrate in a green space is very important since different species generate different amounts of biomass. Indeed, the characteristics of vegetation have proved determinant factors on cooling effect of green spaces. Monteiro et al. [51] highlighted the area occupied by tree species and by lawn, while Potchter et al. [16] verified that green spaces composed by wide, dense and tall tree covers provide a bigger cooling effect during the moments of the day with higher temperatures, unlike green spaces composed by lawn and fewer trees, which were hotter than the surrounding area, especially during the day [52].
Regarding the linear regression models that were generated, it was verified that the use of a statistical method adjusted to biomass generate more satisfactory results than the volume of vegetation (for a vegetation increase of 1 m 2 air temperature will reduce about 0.02 • C). In fact, Kernel density allowed a statistical transformation of biomass spatialization and generated a vegetation gradient which is better correlated with air temperature. Thus, it's important to investigate deeper the relation between these two variables in order to apply this methodology to other urban areas.
The green spaces cooling potential map generated for Lisbon showed that the areas with the highest intensities of UHI (>1.5 • C) had the least amount of biomass and, therefore, the lowest cooling potential. In fact, these areas correspond mainly to the city historic center and the riverside front, with high density of buildings and the lowest amount of available space for increasing vegetation. As one might expect, here the climatic modifications are more pronounced, and this explain the highest densities of UHI and the lowest wind velocity [53].
Several studies have explored the relation between urban green space coverage and some characteristics of urban spaces like city area and population size in order to see if there is enough vegetation for the effective improvement of urban microclimate. Fuller and Gaston [54] evaluated this relationship across 386 European cities and concluded that South European and Mediterranean cities have lowest percentages of green space coverage and per capita green space provision (m 2 per person) compared to North European cities, even though future climate change scenarios anticipate an increase of temperatures in the former regions.

Conclusions
The present study estimated that the amount of AGB in Lisbon is slightly higher in the winter, contrary to what could be expected since the leaf fall occurs at the time of year in which higher volumes were found. Such results can be explained by the fact that biomass was estimated combining all layers of vegetation (arboreal, shrub and herbaceous), so the extensive areas occupied by lawn which are very dry in summer contributed to the apparent reduction on vegetation volume in this season. In a future work we recommended that the estimation of biomass should consider vegetation layers or species.
It was found that an area of 50 m 2 covered by vegetation can cool the air by 1 • C. So, this area can be considered the minimal size of a green space in order for it to make a significant contribution to the improvement of urban microclimate. In fact, the improvement of urban climate is one of the ecosystem services provided by vegetation which contributes to the sustainability of cities. With an increase of vegetation, it was concluded that the urban green spaces cooling potential would be higher, with positive feedbacks on human thermal comfort, contributing for healthier and more resilient cities. Therefore, the authors recommend that this measure should be taken into account when planning new green areas in the city.
Considering that the areas with the highest intensities of UHI (>1.5 • C) correspond to the ones with relatively low cooling potential, we acknowledge that the green space area currently existent in Lisbon is not enough to promote sufficient cooling, especially during summer.
In addition, it is relevant to enumerate the limitations of this work. Regarding AGB, firstly, the models used weren't built considering the topographic and climatic features of the city of Lisbon, so adjustments will be necessary in order to better fit the study area, since the amount and type of vegetation vary geographically. Besides this, it's relevant to collect field samples in the future that would allow the validation of results and the calibration of models. Still, the estimation of AGB from models based on satellite images and other remote sensing products entails numerous limitations associated to the characteristics of satellite/platform, spatial and temporal resolution, etc. In this case, the spatial resolution of Landsat 8 Images (30 m on the visible spectrum) is too coarse for the effective estimation of the volume of vegetation. Lastly, the estimation of the urban green spaces cooling potential for the whole city is based on an extrapolation of a model built with data from only one green space.
This work presents a new methodology for the quantification of the amount of vegetation in urban areas considering their height and volume (aboveground biomass), which was applied in a few cities and is still in the experimentation phase due to its challenges related to the heterogeneity of soil use and fragmentation of green areas. This parameter can be used as a model to quantify the effective reduction on air temperature provided by vegetation.
Considering the present work can be a starting point for the estimation of the cooling potential of urban green spaces in other cities, it would be interesting to hereafter evaluate the quality of ABG estimations, with field data and LIDAR images (height of vegetation), and the quality of the cooling potential, using a mesoscale meteorological network installed in Lisbon This network could also be useful for a further analysis of the green spaces cooling potential and the validation of the cooling potential map generated in this study in order to build a set of recommendations for planning, management and optimization of urban green spaces: What would be a good ratio between the built up area and the green spaces? Is there an optimal spatial configuration of green spaces (e.g., alignments, etc.), including ensuring their spatial contiguity?; Is there a certain structure of vegetation needed for a significant contribution to the improvement of urban microclimate? Are different categories (trees, shrubs, grassland) more efficient than others?