Study of Urban Greenery Models to Prevent Overheating of Parked Vehicles in P + R Facilities in Ljubljana, Slovenia

Parking in park-and-ride (P + R) facilities on the outskirts of a city reduces the traffic inside the cities and follows the principles of sustainable mobility. However, large paved (asphalt) surfaces create urban heat islands (UHI). This causes the temperature to rise in vehicles during full-day parking, which has a negative effect on comfort level and driving performance. This study was conceptualized as two-stage research. The first (preliminary) stage dealt with the measurement and analysis of temperature data at two nearby open parking lots in the city of Ljubljana, one of which was the main research spatial area, P + R Barje (L1), and Trnovo parking (L2), which was used for comparison in the first-stage research. In the preliminary research, we underlined the problem of long-term parking in parking areas exposed to heatwaves (HW). The second stage involved the studying of greening schemes in P + R facilities, which would allow for optimal shading during parking. Advanced 3D and 4D models using ENVI-met and LEONARDO software were developed, which assessed the parking surfaces and the areas of optimum outdoor comfort. Shading by greenery was adopted in this paper, as the strategy aimed at improving the conditions by modelling different variants of greening the parking lot L1.


Introduction
Urban areas are becoming the main target of climate change [1]. The effects of global warming and climate change are expected to be exaggerated in cities, particularly with respect to extreme heatwave events [2][3][4][5]. The phenomenon of urban climate must be comprehended in the local and regional contexts, in which all dynamics and impacts are undertaken in relation with the expansion of urban sprawls. Both the diminishing of green urban areas and the increase of artificial surfaces are aggravating the impact of heatwaves (HW). The greatest HW temperature increases are expected in Central European cities, whereas cities such as Ljubljana, Prague, and Zagreb will experience a significant rise in HW [4]. A number of studies, including observational, modelling, or both, focusing on the micro to the meso-scale, have investigated the effectiveness of using different mitigation strategies to reduce the urban heat island (UHI) in various cities. Strategies aimed at increasing urban vegetation in urban areas are referred to as green infrastructure (GI) strategies [6][7][8][9][10].
Urban areas are characterized by a large proportion of artificial surfaces, such as concrete and asphalt, which absorb and store more heat than natural vegetation, leading to the UHI effect [11]. Measures to heat were significant and proportionally greater in the second half-hour for subjects <40 years and for speeds below 60 km/h (i.e., in city traffic). The latter finding suggests that heat may have increased arousal, and there was some indication of a redistribution of attention away from the most peripheral signals at the higher temperature.
Urban greenery can help to regulate urban microclimates. Green areas are actually the most efficient element that help to reduce the UHI effect [28,29]. Adding more urban trees, parks, gardens, wetlands, and green roofs within urban areas is generally referred to as the implementation of GI strategies. Trees mitigate air overheating, offer shade, and increase the feeling of humidity with their breathing leaves. Vegetation is the necessary element of parking areas, with trees and green areas in between the parking places, as only this can prevent the overheating of concrete, asphalt, and metal surfaces to intolerable temperatures. Importantly, the specificity of P + R sites is the following: the car is parked at the P + R site for several hours, while we return to it from a hard day's work from overheated spaces. In any case, the city is a much friendlier urban place if high-quality green mobility is taken care of. Increasing the vegetation land cover could considerably reduce surface temperatures at parking lots [20,30]; trees of a height of 5-10 m help to control the overheating of surfaces [18]. Gillner et al. [31] recommend decreasing the thermal load in urban areas for future tree planting by choosing species with high cooling potential. The location of trees at a parking lot has an influence on solar exposure reduction. Bajsanski et al. [18] developed and demonstrated an algorithm that optimizes the location of trees, aiming to provide the maximum overshadowing of parking lots.
The main purpose of this paper was to underline the problem of the rate of temperature rise in parked vehicles in large parking areas exposed to sunshine and propose solutions with an optimal introduction of vegetation as protection against the sun. The research of the relation between built and open (void) spaces refers to the effects of surface and atmospheric UHI. The case study of Ljubljana in Slovenia offers an opportunity to explore spatial and microclimate qualities referring to sustainable mobility systems. Besides the pioneering exploration of the aforementioned conditions at open P + R lots in the city of Ljubljana, the novelty of the presented study is the research methodology, which represents a combination of several types of methods, from ex ante problem observation, to experimental work in situ, to computer simulations and modeling, to an ex post multicriteria comparative analysis of the obtained results.

Case Study
The city of Ljubljana is the largest urban area in Slovenia; 14.0% of the total Slovenian population-i.e., 289,830 inhabitants-reside in the city of Ljubljana [32]. The natural conditions in the Municipality of Ljubljana (MOL) carry environmental risks, especially regarding floods and earthquakes. The weather in Ljubljana is conditioned by the city's geographical position in a vast basin extending from Slovenia's pre-Alpine to karst areas. In the study of 571 European cities [4], Ljubljana is one of 28 cities in the top 20% changes in both flooding and HW maximum temperature indices.
Ljubljana was awarded the 2016 European Green Capital title [33], which gives it the designation of a sustainable city. An important criterion for this flattering title was the progress made by the city in terms of greening schemes, sustainable mobility, and energy performance. In the application portfolio, the following information regarding green areas was provided: there are 542 m 2 of public green areas per capita, more than 46% of the urban area is covered by native forest, 75% of all surfaces are green, and the Path of Remembrance and Comradeship (34 km) is the longest city tree line with 7000 trees [33].
The first P + R facility in Ljubljana was built on the Dolgi Most location in 2001 [17]. Currently, there are five active P + R city sites, and P + R Barje represents the largest one. It is positioned three kilometers away from the city center and equipped with elements of sustainable mobility. The Transport Policy, adopted by MOL in 2012 and included in the Integrated Traffic Strategy of Ljubljana Urban Region (LUR) in 2018, puts emphasis on the restriction of personal motor traffic and the change of Sustainability 2020, 12, 5160 4 of 17 travel habits "with improving the distribution of mobility by one third of journeys by public transport, one third by bicycle and on foot and one third by personal vehicle" [34]. The sixth Ljubljana's P + R facility will be completed during 2020 in Stanežiče [35].
In relation to the meaning of P + Rs for urban planning, Žaucer et al. emphasize the following: "the multifunctionality of the space, the mixing of uses and the appropriate interweaving of functions enable the P + R node to be alive, pleasant and safe regardless of the part of the day, time of year or day of the week. A successful P + R hub is an urban area integrated into a settlement. It is a public space, which is also a place to meet, socialize and create mutual ties in the local community" [17]. Until today, however, Ljubljana has not developed a model for P + R areas, and their evolution rather follows other spatial interventions in the city. There is no clear position on whether multifunctional P + Rs are more appropriate, nor does there exist guidelines for the multifunctional combining of the importance of locations. On the other hand, the significance of existing P + Rs in the city can be proved by analyzing the public traffic network. The largest P + R Barje is directly linked to the highway and the city ring, thus, covering the arrival to Ljubljana from the Western to the Southwest regions of Slovenia, as well as the LUR area tied to this part. This represents more than one third of the territory of Slovenia. Furthermore, P + R Barje is positioned beside road traffic connections with Italy and Croatia.
Most of the publicly available parking areas in Ljubljana are surrounded, at least on two sides, with green areas, which indirectly beneficially reduce overheating. The construction of additional P + R parking areas on the city outskirts should reduce the occupancy of parking areas in the center and the pressures in peak periods, as estimated by MOL. The P + R occupancy rate in Ljubljana depends on several factors, such as the trends of using the highway system, the changing of the local transport, the price of city parking areas, access to desired locations, etc. P + Rs are locations of static activities: i.e., cars are parked for 8-10 h nowadays on the average workday, plus the outbound and inbound travel to the site. Any car parked here absorbs the quality of the urban climate through external surfaces of sheet metal. At the same time, the degree of overheating is returned inside and directly affects human well-being.
For this research case study, P + R Barje in Ljubljana ( Figure 1) was selected (area of 208 ha, constructed in 2015). A particularity of this parking area is the soft Barje ground, which necessitated the installation of a concrete foundation slab with a surface area of 13,000 m 2 , supported by 1117 piles. All the necessary installations were added: sewage, public lighting, optical connection for Urbanomats, and a Bicikelj self-service bike rental point. of travel habits "with improving the distribution of mobility by one third of journeys by public transport, one third by bicycle and on foot and one third by personal vehicle" [34]. The sixth Ljubljana's P + R facility will be completed during 2020 in Stanežiče [35].
In relation to the meaning of P + Rs for urban planning, Žaucer et al. emphasize the following: "the multifunctionality of the space, the mixing of uses and the appropriate interweaving of functions enable the P + R node to be alive, pleasant and safe regardless of the part of the day, time of year or day of the week. A successful P + R hub is an urban area integrated into a settlement. It is a public space, which is also a place to meet, socialize and create mutual ties in the local community" [17]. Until today, however, Ljubljana has not developed a model for P + R areas, and their evolution rather follows other spatial interventions in the city. There is no clear position on whether multifunctional P + Rs are more appropriate, nor does there exist guidelines for the multifunctional combining of the importance of locations. On the other hand, the significance of existing P + Rs in the city can be proved by analyzing the public traffic network. The largest P + R Barje is directly linked to the highway and the city ring, thus, covering the arrival to Ljubljana from the Western to the Southwest regions of Slovenia, as well as the LUR area tied to this part. This represents more than one third of the territory of Slovenia. Furthermore, P + R Barje is positioned beside road traffic connections with Italy and Croatia.
Most of the publicly available parking areas in Ljubljana are surrounded, at least on two sides, with green areas, which indirectly beneficially reduce overheating. The construction of additional P + R parking areas on the city outskirts should reduce the occupancy of parking areas in the center and the pressures in peak periods, as estimated by MOL. The P + R occupancy rate in Ljubljana depends on several factors, such as the trends of using the highway system, the changing of the local transport, the price of city parking areas, access to desired locations, etc. P + Rs are locations of static activities: i.e., cars are parked for 8-10 h nowadays on the average workday, plus the outbound and inbound travel to the site. Any car parked here absorbs the quality of the urban climate through external surfaces of sheet metal. At the same time, the degree of overheating is returned inside and directly affects human well-being.
For this research case study, P + R Barje in Ljubljana ( Figure 1) was selected (area of 208 ha, constructed in 2015). A particularity of this parking area is the soft Barje ground, which necessitated the installation of a concrete foundation slab with a surface area of 13,000 m 2 , supported by 1117 piles. All the necessary installations were added: sewage, public lighting, optical connection for Urbanomats, and a Bicikelj self-service bike rental point.

Methods and Tools
The materials and tools for this study consisted of surveys of key precedents, detailed mapping, and modelling using ENVI-met. The study implemented morphological and typological research methodologies, field observations, and on-site analyses, with the data used as follows: • Geomorphological data by the Surveying and Mapping Authority of the Republic of Slovenia were used [36]. Analyses of geomorphological (topography, reliefs, terrain type, and climate conditions), and mapping of landscape characteristics were conducted. Spatial data were mapped using GIS ESRI ArcGIS software (ArcMap version 10.3.1) for visual representations and analyses. • Thermal satellite imagery was done based on Landsat 8 data obtained from the satellite thermal infrared sensor (TIRS) [37]. The Landsat 8 satellite (launched as the Landsat Data Continuity Mission-LDCM on 11 February, 2013) mission continues the acquisition of high-quality data that meet both NASA and USGS scientific and operational requirements for observing land use and land change [37]. We used Landsat 8 images for the production of the thermal map. The images were processed using Esri ArcGis Software, which was transformed in a way to calculate the temperature for each pixel separately by converting the raw bands into top of atmosphere radiance (TOAr). The temperature calculations were done with the Raster calculator tool, which calculated the temperature values for each pixel separately. For the visual representation, the initial 30 × 30 m grid was resampled using the ArcMap Resample tool to a 1 × 1 m grid. A bilinear resampling technique was used.

•
Remote sensing images: Google satellite images of P + R Barje (L1) and Trnovo parking (L2) [38] ( Figure 2) and orthophoto images from the remote sensing database, sheet orthophotos DOF050 2018 [36] were used. Remote sensing data were used for a manual line-art drawing of the spatial situation in Autodesk AutoCAD 2019. Furthermore, we elaborated a spatial representation of the P + R Barje parking area. Sheets of DOF050 images covering the area were combined and clipped to the study area using ESRI ArcGIS 10.6 and the ArcToolbox clip tool.

•
The ENVI-met [39] database of materials and 3D vegetation was used in ENVI-met SCIENCE in ENVI-met SPACES for the manual drawing of the cell map in the size of 100 × 100 × 20 pixels. We mapped the existing vegetation and surfaces using input cartographic data.

Preliminary Research
The study on ambient temperatures inside vehicles in large parking areas exposed to sun was done in the preliminary research step to identify the relationship between UHI areas and open parking spaces. The preliminary research included the study on temperature measurements and

Preliminary Research
The study on ambient temperatures inside vehicles in large parking areas exposed to sun was done in the preliminary research step to identify the relationship between UHI areas and open parking spaces. The preliminary research included the study on temperature measurements and thermal comfort. When checking the temperature conditions, we compared two parking areas that were 1300 m apart. The parking areas had the same orientation and other boundary conditions, but they differed in terms of their greening schemes. Parking area L1 contained trees that had no canopy and did not provide protection against sun, so the vehicles were exposed to direct solar radiation all day long; L2 had dense trees and the vehicles were in the shade throughout the day (Figure 3). During the parking, the outside and vehicle cabin ambient temperatures were taken at the height of the driver's head. During the night, the vehicles were parked in a garage (sheltered area) at a distance of 2200 and 3500 m, respectively, from the selected parking areas (a 6-and 8-min drive, respectively), for about 6 h each night.
To focus on the worst-case scenario, the preliminary study went underway in summer months, which is when the temperatures in vehicles rise the most [26]. The temperature was measured upon arrival to the parking area and before the departure. Days with a different cloud coverage during parking were selected. Only air temperatures were measured, although the comfort level in a vehicle is also influenced by the relative air humidity and air speed. Under such high thermal stress, air temperature has the highest influence on a person's comfort level and driving performance.

Elaboration of Greening Scheme Models at P + R Sites
This study focused on P + R Barje (L1) (Figures 2 and 3), analyzing the various scenarios of introducing vegetation as a protection against sun. The scenarios for developing the greening schemes were investigated using microclimate simulations by applying a numerical model at the micro level. The ENVI-met software allowed for the simulation of the progression of various thermodynamic parameters in a 3D space (cell grid, where each cell is 2 × 2 × 2 m) between morphology-microclimate and vegetation. The programmatic model considered the physical considerations based on the principles and laws of fluid mechanics and thermodynamics used to calculate three-dimensional wind fields, turbulence, air temperature and humidity, radiation flows, and pollutant dispersion [40,41]. The central part of the P + R Barje parking area is in the South of Ljubljana, at latitude 46 and longitude 14.5. The model was presented at a 3D level used to simulate interactions between the surface area, vegetation, and atmosphere in urban areas in a rectangular Two cars with similar characteristics were used in the study: Peugeot 208, white-colored, with dark grey textile seats. In both cars, the windows were closed at all times. The vehicles were parked in both parking areas at the same time: (1) outdoor parking area L1 without trees, where the vehicle was parked in direct sunlight, with no shade, and (2) parking area with trees L2, where the vehicle was parked in the shade of trees all day long ( Figure 3). The orientation of parked vehicles was the same in both parking areas. In line with the objective of this paper, warm days with various temperatures and cloud coverage were selected: i.e., when temperatures inside vehicles rise significantly. The windows in the vehicles were closed during the parking. Other studies also have found that cracking windows by a few centimeters does not contribute to continuous cooling, as temperatures stay high inside vehicles after several hours of parking under direct sunlight [22].
During the parking, the outside and vehicle cabin ambient temperatures were taken at the height of the driver's head. During the night, the vehicles were parked in a garage (sheltered area) at a distance of 2200 and 3500 m, respectively, from the selected parking areas (a 6-and 8-min drive, respectively), for about 6 h each night.
To focus on the worst-case scenario, the preliminary study went underway in summer months, which is when the temperatures in vehicles rise the most [26]. The temperature was measured upon arrival to the parking area and before the departure. Days with a different cloud coverage during parking were selected. Only air temperatures were measured, although the comfort level in a vehicle is also influenced by the relative air humidity and air speed. Under such high thermal stress, air temperature has the highest influence on a person's comfort level and driving performance.

Elaboration of Greening Scheme Models at P + R Sites
This study focused on P + R Barje (L1) (Figures 2 and 3), analyzing the various scenarios of introducing vegetation as a protection against sun. The scenarios for developing the greening schemes were investigated using microclimate simulations by applying a numerical model at the micro level. The ENVI-met software allowed for the simulation of the progression of various thermodynamic parameters in a 3D space (cell grid, where each cell is 2 × 2 × 2 m) between morphology-microclimate and vegetation. The programmatic model considered the physical considerations based on the principles and laws of fluid mechanics and thermodynamics used to calculate three-dimensional wind fields, turbulence, air temperature and humidity, radiation flows, and pollutant dispersion [40,41]. The central part of the P + R Barje parking area is in the South of Ljubljana, at latitude 46 and longitude 14.5. The model was presented at a 3D level used to simulate interactions between the surface area, vegetation, and atmosphere in urban areas in a rectangular grid with a variable accuracy and 3D distances in x, y, and z directions. The parking area extended in a 150 × 200 m area, which was transferred into a digital model of 100 × 100 × 20 cells. This resolution allowed for the analysis of the interactions between the various buildings, surface, and plants for various scenarios in a 24 h period.
Parking spaces are organized in the West-East direction with green areas, with trees parallel to the parking spaces. The parking area is located in an open plain area surrounded by a low lawn in the South and Western sides, asphalt surfaces in the North, and a green area and trees in the Eastern side. Based on microclimate parameters, structure, and the composition of the ground surface, ENVI-met [39] was used to simulate alternative scenarios of vegetation and the organization of parking spaces.
The digital surveying groundwork data for the parking areas were obtained from the [36]. For the simulation model, relevant microclimate data from the nearest weather station, Ljubljana Bežigrad [42], were used in order to study how urban microclimate develops under the current conditions and how vegetation improves the thermal comfort at the P + R Barje parking area. We used the graphic interface ENVI-met SPACES for the manual plotting of the site as cells, along with the value parameters containing information on plants, land cover, and materials. All cartographic data first had to be collected and used for input into the ENVI-met graphic editor, where the basic databases were used in the first stage as generic values of the initial simulations (for the following elements: soil and surface, vegetation, buildings, and greenings). The maps with localized 3D data were then transferred into the model, where new input data on the atmosphere for modelling the context were added, which was the second stage of simulations. All the input data used were taken from the day when the temperatures were measured: i.e., 29 June 2018. The first phase was intended for the accurate modelling of the static context, while the second phase included some variables connected to surface and atmospheric temperature, such as solar radiation duration, wind direction, humidity, etc. The variant solutions were limited to the distribution of trees and green walls of a height of 3 m. Green canopies were not included in the simulation because of the software limitation of spatial representations with cells and the available size of the model with 100 × 100 × 20 cells.

Preliminary Research Outcome
In June, 2018, Ljubljana had seven hot days with temperatures exceeding 30 • C and 20 days with temperatures of 25 • C and above [43]; in July, 2018, there were 10 days with temperatures above 30 • C and 24 warm days above 25 • C [44]; and there were 18 days with temperatures above 30 • C and 27 days with temperatures 25 • C and above in August, 2018 [45]. Together, with the position of open parking lots L1 and L2 in the dynamic UHI map of Ljubljana (Figure 2), the conducted measurements confirmed that the temperature in the vehicle parked at the sun-exposed lot reached values that caused • the difference between outdoor morning temperature at locations L1 (strongly sun-exposed P + R Barje) and L2 (well-landscaped Trnovo parking), measured at 8.00 a.m., was small (from 0.4 • C min. to 1.4 • C max.); • the temperature in the cabin of the sun-exposed vehicle at L1 was significantly higher than the temperature in the vehicle in the shade at L2 (both were measured at 2.00 p.m. during the whole examination period); • on July 7 at 2.00 p.m., the temperature in the cabin of the vehicle parked at the exposed lot L1 was high even in partly cloudy weather; • on 29 June, 3 July, and 4 July at 2.00 p.m., the temperature in the cabin of the car parked at L1 reached a dangerous level despite the more favorable external conditions (cloudy weather and milder outside air temperature) ( Table 1); • at 2.00 p.m. on the excessively hot day of 7 August, the temperature inside the car parked at L1 was extremely high, thus, endangering the health and well-being of vehicle users. The greenery on the parking lot played a key role in preventing the overheating of the vehicle cabin. The temperature in the cabin of the vehicles parked in greened parking lot L2, measured at 2.00 p.m., was always significantly lower than the temperature in lot L1 without shading. To avoid the overheating of a vehicle cabin, it seems vital to provide shading during whole parking time. Such a prerequisite was adopted and implemented in the second stage of the research, which was modelling the greenery in lot L1 to improve temperature-related conditions.

Modelling the Greenery at P + R Barje (L1)
The central part of the study dealt with the design and comparison of different variants of greenery distribution at L1. By producing advanced 3D and 4D numerical models with the software tools ENVI-met and LEONARDO (Table 2), we evaluated the parking lot against the areas of optimal outdoor comfort and minimum environmental damage caused by the designed solutions. Table 2. Input data of the study area in ENVI-met 1 .

Input Data
Latitude and longitude  To control the relevance of shading and the temperature, we selected the date of 29 June, 2018 as the first measurement day at both L1 and L2 locations. Being close to the annual solstice, the angle of incidence of solar radiation was high on this day. The decrease of the angle in the following days was considered more favorable for shading by greenery. The initial analysis provided an insight into modelling possibilities and helped to hypothesize alternative planting systems, the position of green surfaces, and the use of shading bodies. A comparative analysis of designed greening models allowed us to determine the optimal shading of the parking lot and, consequently, to reduce the maximum surface and air temperature of the parking area at the height of 1.2 m.
In the outcome of the modelling process were four solutions for greening ( Figure 4). In three of these, the existing functional organization of the parking lot was preserved (Models A, C, and D), while in Model B, the layout was changed but the number of parking places remained unchanged.
In terms of greening elements, the modelling included trees with different layouts and vertical green walls with different orientations.
The main tasks in defining four diverse models were to ensure that the position of the inserted greening elements did not impact the capacity of the parking area, and to examine how their shadow changed throughout the day. Therefore, the adjustment of the position of greening elements did not impact the organizational quality of the parking lot, but it did cause changes in the air temperature.
By comparing the used greening elements-trees vs. green walls-we found that a tree provided better conditions for creating shade, while a green wall allowed a more effective temperature decrease in the parking lot, at the height of 1.2 m.
In terms of greening elements, the modelling included trees with different layouts and vertical green walls with different orientations.
The main tasks in defining four diverse models were to ensure that the position of the inserted greening elements did not impact the capacity of the parking area, and to examine how their shadow changed throughout the day. Therefore, the adjustment of the position of greening elements did not impact the organizational quality of the parking lot, but it did cause changes in the air temperature.  The designed models (Figure 4, Tables 3 and 4) featured the following: • Model S stands for the current situation at site L1. We included trees (44) that were there in 2018 and are still there now (May 2020). Model S2 did not include the existing large trees outside L1; however, the model detected them (there was, however, no direct impact on the shading, depiction Model S3, while there was an impact on temperatures in Model S4).

•
Model A presents the variant of tree densification in direction N-S (between the individual lines of parking). Trees in the E-W direction were added. The parking design affected the reorganization of the individual parking areas, which was still five parking spaces per unit. In Model A5 (in Table 5), we see that the shading was distributed evenly, meaning that there were several quality conditions of shading the parking areas. According to Model A3, we found that higher temperatures were evenly distributed across the parking area. This was the result of an increased number of trees and airflow stabilization. In this case, the number of trees did not bring about a reduction in temperatures.

•
Model B also presents the variant of tree densification in the N-S direction, but with additional trees along a zigzagging line, with one tree per five parking spaces, switching to the East and West orientation. The design of the parking was not affected by this; the parking spaces remained organized in groups of five. The variant had an interesting composition, but in the horizontal direction, it created major islands without trees (compared to Model A), causing less favorable shading with trees and a poorer distribution of quality parking spaces.
• Model C involved green walls between the parking lines in a N-S direction. The design of the parking remained the same, while a passage between the groups of five parking spaces must be provided. The model confirmed the thesis that green walls favorably affect air temperatures; however, in the sense of creating well-shaded parking areas (Model C5), it was less effective than Model A5 and B5. Interestingly, with this positioning, we managed to avoid all the worst shading positions. This variant also significantly affected the lowering of temperatures as, compared to Model A and B, the temperatures were, on average, lower by 2 • C.

•
Model D was an organization variant of the parking area, with walls installed in the E-W direction. The parking area design was not changed with this; the discontinuity of the green wall between the groups of five parking spaces was not necessary. The model negated the vision of the significance of green walls, as shades in Model D5 were positioned the worst compared to all the previous models. The importance of the orientation of the greening element related to the rate of sunshine (sun hours) was confirmed. Air temperatures were less favorable in Model D3 than in Model C and Model C3 (which was related to the element of the green wall in N-S orientation).  Model X ( Figure 5). The assessment of the results included the significance of a different kind of designing P + Rs by offsetting the position of roads at an angle of 5°-10° and the inclusion of larger green areas in-between ( Figure 5), providing better tree growth and the possibility of tree densification (comparison in Table 5, Figure 6).

Discussion
Further evaluation and optimization of the developed models can be done on the basis of the "aesthetics," "amenity," and "shading benefits" that were explained by Nieuwenhuijsen [46] in the context of urban and transport planning for livable and healthy cities as elements that "have also been shown to provide significant economic benefits, while benefits in terms of water regulation, carbon reduction and air quality are usually more modest" [47].
To assess the visual quality of the derived variants, trees and living walls on the parking lot and trees as a boundary condition should all be taken into consideration. Two main assessment principles can be established: the level of the visual perception of the parking lot from afar, and the level of the visual perception of the belonging parking system-i.e., the density of greening elements. Following these criteria, it can be concluded that densely green Model A and B (Figure 4; Models A1 and B1) interrupt the visual continuity, but at the same time offer a pleasant feeling. Models C and D provide for the inclusion of the green wall element, which is of interest in terms of composition, but does not improve shading. Model C is approaching the values of shading Models A and B (Figure 4-Models A5, B5, and C5); however, in terms of area perception, this implies restricted visual contacts. According to all the conditions investigated, Model D is the worst variant, even though the position of elements in the E-W direction allows for a visual connection when approaching the P + R (the access road is in the N-S direction). Even though the green wall is interesting, it is, in fact, less suitable in terms of maintenance. For discussion reasons, we added Model X, which offers tree densifications, creating the perception of green/forest islands. Model X does not follow the existing parking area

Discussion
Further evaluation and optimization of the developed models can be done on the basis of the "aesthetics," "amenity," and "shading benefits" that were explained by Nieuwenhuijsen [46] in the context of urban and transport planning for livable and healthy cities as elements that "have also been shown to provide significant economic benefits, while benefits in terms of water regulation, carbon reduction and air quality are usually more modest" [47].
To assess the visual quality of the derived variants, trees and living walls on the parking lot and trees as a boundary condition should all be taken into consideration. Two main assessment principles can be established: the level of the visual perception of the parking lot from afar, and the level of the visual perception of the belonging parking system-i.e., the density of greening elements. Following these criteria, it can be concluded that densely green Model A and B (Figure 4; Models A1 and B1) interrupt the visual continuity, but at the same time offer a pleasant feeling. Models C and D provide for the inclusion of the green wall element, which is of interest in terms of composition, but does not improve shading. Model C is approaching the values of shading Models A and B ( Figure 4-Models A5, B5, and C5); however, in terms of area perception, this implies restricted visual contacts. According to all the conditions investigated, Model D is the worst variant, even though the position of elements in the E-W direction allows for a visual connection when approaching the P + R (the access road is in the N-S direction). Even though the green wall is interesting, it is, in fact, less suitable in terms of maintenance. For discussion reasons, we added Model X, which offers tree densifications, creating the perception of green/forest islands. Model X does not follow the existing parking area design (the starting condition in the study was the existing design); nevertheless, it is important in the sense of comparing green areas and the number of trees. The tree concentration is suitable; however, the (non-functional) driving dynamics is put under question. Shading brings an added aesthetic value to a P + R facility, as it positively affects the perception of space under extreme heat. Models A and B offer the most shaded areas; from the point of view of a visual continuity, it would be recommended to remove the trees representing the boundary conditions of the site. Given the number of trees and considering all the criteria, Model B is the most rational one and gives an added value to the artistic organization of the green element.
In this study, we addressed the thermal comfort level related to the rise in temperatures inside vehicles parked at a P + R facility and the threats to health and drivers' psychophysical performance. Moreover, it is important to understand the interplay of the results given in Figure 4 with temperatures and shading. The sense of ventilation, airiness, tree density, etc., are the elements that affect the comfort level at a P + R facility. It is necessary to explain that all these elements decrease the temperatures on the ground (as indicated by the models); however, trees are three times more beneficial than other green elements. The aforementioned elements and model depictions confirm the positive roles of Models A and B, but with restrictions referring to tree density and the visual significance of the P + R role, which is not a green area, per se.
Adequate green elements introduced in an open parking area not only provide shade and prevent the overheating of parked vehicles, but also contribute to the reduction of the outdoor air temperature due to evaporation. The temperature values of the areas shaded in vegetation were perceived as the most comfortable in Model C (Figure 4, Model C5). Nevertheless, given the progression of shading, we find that the green wall as a shading element present in this, as well as in Model D, does not adequately support the P + R organization. Furthermore, shading by a green wall does not affect a person in the way that a cool tree canopy would (it is more like shading provided by a building). The most favorable natural shading element would be a pergola, which would cover the whole parking area and, thus, provide permanent shading. However, it is necessary to note that this decreases the significance of the elements presented in Sections 3.1 and 3.2 of this paper.
Overall, given the primary function of an open parking area, this study demonstrated the need for optimization between the capacity and organization of a parking lot, on the one hand, and the type and spatial distribution of green elements that provide shading, on the other. In the existing parking area P + R Barje, the trees were identified as the most important element in providing shading in the variant when the shading is provided between parking spaces.

Conclusions
This study revealed the interrelations between the greenery in open parking areas, the temperature inside parked vehicles, and the local UHI effect. This study was conceptualized as two-stage research. The first (preliminary) stage dealt with the measurement and analysis of temperature data at two nearby open parking lots in the city of Ljubljana, one of which was the main research spatial area (P + R Barje, L1), and the other (Trnovo parking, L2) was used for comparison. The second phase of the study was dedicated to modelling the greenery at the principal location L1. Accordingly, the presentation of findings in this paper referred separately to each of the two stages. Progressive result generation allows for adjustments and, thus, enables a wider application of the designed research methodology.
Like most European capitals, Ljubljana is susceptible to the manifestation of UHI [15]. While urban morphology initially influences the occurrence and intensity of UHI, ultimate consequences are felt by the users of the space. As demonstrated in this paper on the example of two nearby open parking locations in Ljubljana that were similar in all parameters except for the density of the greenery, the UHI effect is more intensive on the lot characterized by sealed artificial surface and scarce greenery (location L1), in spite of adjacent natural green areas. Likewise, the temperature in the cabin of the vehicle parked at location L1 was significantly higher in comparison with location L2, which had sufficient greenery to provide constant shade to parked vehicles.
Parking the vehicles in the open space during the warmer part of the year can easily result in the overheating of the cabin, especially when the duration of parking time is longer. The measurements showed that parking on an open surface without solar protection raises the temperature in the cabin to such an extent that it becomes dangerous for the health of the driver and other passengers in the vehicle, but also for traffic safety. For illustration, on the warmest day of measurement period, when outdoor air temperature was 35.5 • C, the temperature of the air in the vehicle cabin rose to 59.1 • C. High air temperatures in vehicles were measured also during moderately hot sunny days, which was also found by other authors [22,23].
The key measure to prevent overheating in the cabin of the vehicle parked in an open space is shading. Though the shade can be produced by natural (e.g., trees) and man-made objects (e.g., the overhang) [21], greenery offers far more benefits: inter alia, the reduction of the aforementioned UHI effect [6-10]. Therefore, shading by greenery was adopted in this paper as the strategy aimed at improving the conditions by modelling different variants of greening parking lot L1.
Another important input to modelling was the fact that the vehicle must be parked in shade at all times, without (even short-term) exposure to direct sunlight. The desired state currently represents a challenge to P + R lots in many European cities. To that end, a transitional practical solution could concern the categorization of individual parking places according to the estimated length of parking time. In research terms, the continuation of the presented study could deal with the length of time during which the indoor air temperature of a vehicle parked in an open space rises above the recommended thresholds, examined against various shading-related conditions. An alternative research direction that would help to strengthen the validity of the obtained results refers to the in-depth analysis of drivers' behavior at open P + R surfaces in terms of selecting a parking place, and the practice of precooling the vehicle cabin before leaving the parking area.