The Multiple-Scale Nature of Urban Heat Island and Its Footprint on Air Quality in Real Urban Environment

: The complex interaction between the Urban Heat Island (UHI), local circulation, and air quality requires new methods of analysis. To this end, this study investigates the multiple scale nature of the UHI and its relationship with ﬂow and pollutant dispersion in urban street canyons with and without the presence of vegetation. Two ﬁeld experimental campaigns, one in summer and one in winter, were carefully designed in two parallel urban street canyons in the city of Bologna (44 ◦ 29 (cid:48) N, 11 ◦ 20 (cid:48) E; Italy) characterized by a similar orientation with respect to the impinging background ﬂow but with a di ﬀ erent aspect ratio and a di ﬀ erent presence of vegetation. In addition to standard meteorological variables, the dataset collected included high-resolution ﬂow data at three levels and concentration data of several pollutants. The UHI has been evaluated by combining surface temperature of building facades and ground surfaces acquired during two intensive thermographic campaigns with air temperature from several stations in order to verify the presence of intra-city neighborhood scale UHIs additional to the more classical urban–rural temperature di ﬀ erences. The presence of trees together with the di ﬀ erent morphologies was shown to mitigate the UHI intensity of around 40% by comparing its value in the center of the city free of vegetation and the residential area. To capture the multiple-scale nature of UHI development, a simple relationship for the UHI convergence velocity, used as a surrogate for UHI strength, is proposed and used to establish the relationship with pollutant concentrations. The reliability of the proposed relationship has been veriﬁed using a Computational Fluid Dynamics (CFD) approach. The existence of a robust relationship between UHI strength and pollutant concentration may indicate that the positive e ﬀ ect of mitigation solutions in improving urban thermal comfort likely will also positively impact on air pollution. These results may be useful for a quick assessment of the pollutant accumulation potential in urban street canyons.


Introduction
More than half of the world's population live in urban areas, and especially in highly-dense cities, a proportion that is expected to increase to 68% by 2050 [1]. Long term projections indicate that urbanization together with the overall growth of the world's population could add further 2.5 billion people to urban areas by 2050, with percentages reaching 90% in Asia and Africa [1].
Increasing urbanization and associated change in land surface physical properties such as roughness, thermal inertia, and albedo [2] and land surface processes can deeply affect regional and local scale meteorology and air quality, through alteration of the surface-atmosphere coupling As discussed above, several methods of investigations have been used in the literature to study the UHI, the UPI, and their inter-linkages with urban flow and dispersion. As far as the numerical approaches are concerned, methods rely on both atmospheric models such as WRF (Weather Research and Forecasting), e.g., [33][34][35] and Computational Fluid Dynamics (CFD) approaches, e.g., [36][37][38]. While the first category is able to capture the large-scale interaction with the local features by solving the complete set of equations that capture the full relationship between surface and a rotating atmosphere, the latter uses generic thermo-fluid dynamics equations having the advantage of characterizing individual obstacles effects on flow and dispersion. Specifically, CFD approaches have been used to evaluate the role of canyon related features, such as the effect of aspect ratio and wind/ventilation characteristics on the urban thermal microclimate, UHI, heat, and mass transport, e.g., [39][40][41][42][43][44]. Within the present study, CFD simulations are complementary to the data analyses and used only for the purpose of assessing the robustness of a simple relationship between temperature gradients and flow and concentration field. To our knowledge this is the first study of this kind relating the multiple scale-nature of the UHI with the flow dynamics, providing a relationship between UHI and UPI in urban street canyons, and using CFD simulations to strengthen the validity of such a relationship. Following the above, this paper investigates the link between UHI and pollutant concentration in urban street canyons, assessed with and without vegetation. To evaluate how the two effects are linked, two experimental field campaigns were carried out in the city of Bologna ( Figure 1) in Italy (44 • 29 N, 11 • 20 E) within the "iSCAPE" ("Improving the Smart Control of Air Pollution in Europe") H2020 research and innovation project. Specifically, we try to respond to the following questions: • Could the UHI phenomenon be a multiple-scale phenomenon? What is the impact of the UHI at different scales on the flow dynamics? Are there any differences in behavior? • As UHI and UPI are commonly linked, is it possible to determine a relationship linking these two aspects at a very local scale? How the associated circulation would develop and how can relate to local pollutant concentrations?
Atmosphere 2020, 11, x FOR PEER REVIEW 3 of 32 As discussed above, several methods of investigations have been used in the literature to study the UHI, the UPI, and their inter-linkages with urban flow and dispersion. As far as the numerical approaches are concerned, methods rely on both atmospheric models such as WRF (Weather Research and Forecasting), e.g., [33][34][35] and Computational Fluid Dynamics (CFD) approaches, e.g., [36][37][38]. While the first category is able to capture the large-scale interaction with the local features by solving the complete set of equations that capture the full relationship between surface and a rotating atmosphere, the latter uses generic thermo-fluid dynamics equations having the advantage of characterizing individual obstacles effects on flow and dispersion. Specifically, CFD approaches have been used to evaluate the role of canyon related features, such as the effect of aspect ratio and wind/ventilation characteristics on the urban thermal microclimate, UHI, heat, and mass transport, e.g., [39][40][41][42][43][44]. Within the present study, CFD simulations are complementary to the data analyses and used only for the purpose of assessing the robustness of a simple relationship between temperature gradients and flow and concentration field. To our knowledge this is the first study of this kind relating the multiple scale-nature of the UHI with the flow dynamics, providing a relationship between UHI and UPI in urban street canyons, and using CFD simulations to strengthen the validity of such a relationship. Following the above, this paper investigates the link between UHI and pollutant concentration in urban street canyons, assessed with and without vegetation. To evaluate how the two effects are linked, two experimental field campaigns were carried out in the city of Bologna ( Figure 1) in Italy (44°29′ N, 11°20′ E) within the "iSCAPE" ("Improving the Smart Control of Air Pollution in Europe") H2020 research and innovation project. Specifically, we try to respond to the following questions: • Could the UHI phenomenon be a multiple-scale phenomenon? What is the impact of the UHI at different scales on the flow dynamics? Are there any differences in behavior? • As UHI and UPI are commonly linked, is it possible to determine a relationship linking these two aspects at a very local scale? How the associated circulation would develop and how can relate to local pollutant concentrations?  With the intent to respond to the previous questions, the paper is structured as follows. After the Introduction section, Section 2 describes the two field campaigns, detailing the characteristics of the two urban street canyons where the two campaigns were carried out and the instrumentation therein deployed. Section 3.1 focuses on the evaluation of the multiscale nature of the UHI effect at urban and neighborhood levels. Section 3.2 then derives a relationship between UHI and pollutant concentrations in urban street canyons, verified by means of high-resolution CFD simulations adequately setup. Finally, Section 4 draws the main conclusions of this work.

Rationale for the Study and Site Description
Because of the inherent complexity and multiple interactions between physical processes at different spatial and temporal scales, the study of the urban environment requires a multiapproach investigation. Indeed, flow dynamics and turbulence structures are well known to depend on spatial and temporal scales within the city [45,46], time of the day and insolation [47], and even traffic regimes [48,49]. Atmospheric processes such as ventilation, diffusion or flow circulation can affect or be affected by vegetation and urban heat island effects. With the aim to disentangle the contribution of UHI and pollutant concentration, in this study the analysis of field data collected during two experimental field campaigns was complemented with numerical simulations.
In particular, two intensive experimental field campaigns were carried out respectively during summer (10 August 2017-24 September 2017) and winter (16 January 2018-14 February 2018) to measure atmospheric flows and turbulence, air temperature and air quality within two different street canyons embedded in the city of Bologna. The city is located at the southern end of the Po Valley at the foothill of the Apennines chain, a region considered a hotspot for climate change and air pollution [50]. The city morphology can be considered as typical for European cities [51,52], characterized by a densely built historical center surrounded by residential suburbs and a small industrial area.
The center is characterized by a dense network of narrow street canyons occasionally interrupted by small parks or squares, while the residential suburb streets retain the canyon-like shape but diminishing the building packaging. Two parallel street canyons running north to south, Marconi St. and Laura Bassi St., were chosen as experimental sites for the campaigns, being the first representative of the historical area, and the second of the suburbs ( Figure 2). Laura Bassi St. is a 700 m long street surrounded by both small private houses (2-3 floors) and higher buildings (4-5 floors), accounting for a mean building height H of 17 m. The street is composed of two lanes surrounded by an almost regular deployment of deciduous trees, for a mean street width W of 25 m. Marconi St., about 600 m long, is surrounded by buildings with at least 4-5 floors up to 9-10 floors (H = 33 m). The street is composed of four lanes, two for private and two reserved for public transport (W = 20 m). Along the carriageway, there are porticos, covered pathways with sideways. No vegetative element is planted along this road, except for the last 50 m approaching one street end, where there is a single lane of deciduous trees placed on one side of the road. Top views of both canyons are pictured in Figure 2. The majority of the vegetative elements in Laura Bassi St. neighborhood belongs to the Platanus Acerifolia Mill family, a hybrid of Platanus Orientalis, and Platanus Occidentalis, widely spread in European urban habitats. The average height of the branch-free trunk is about 5 m and the crown extends up to 15 m in height. Along the street, tree trunks are regularly spaced at approximately 8 m distance, allowing different crowns superimposition.

Instrumentation Setup
As previously described, the main scope of this work is to investigate the analysis of the UHI at different spatial scales, and to address the interaction between UHI and UPI effects. As such, and according to the indications of [53] describing the need for urban-specific measurement techniques to provide guidance to policy and decision makers for the design of sustainable cities, we utilized observations from a network of meteorological and air quality monitoring stations covering the various spatial scales involved in the investigations (see Figure 3).

Instrumentation Setup
As previously described, the main scope of this work is to investigate the analysis of the UHI at different spatial scales, and to address the interaction between UHI and UPI effects. As such, and according to the indications of [53] describing the need for urban-specific measurement techniques to provide guidance to policy and decision makers for the design of sustainable cities, we utilized observations from a network of meteorological and air quality monitoring stations covering the various spatial scales involved in the investigations (see Figure 3).
Atmosphere 2020, 11, x FOR PEER REVIEW 6 of 32 As previously described in Section 2.1, the main investigation sites are the street canyons of Marconi St. and Laura Bassi St., where a full set of instrumentations was deployed to measure both meteorological and air quality quantities. Supporting measurements of ancillary variables are retrieved from the local environmental protection agency ARPAE (Agenzia Regionale per la Prevenzione, l'Ambiente e l'Energia dell'Emilia-Romagna) network of monitoring stations inside the  The main sites of the campaigns were extensively instrumented to cover a large variety of the local-scale dynamics and thermodynamics of vegetated and non-vegetated urban canopies, and to explore the links with air quality. Three sites were deployed at three different heights within and above each canyon, named ground, mid-canyon, and rooftop levels. The instrumentation deployed in each canyon is specular, but the installation heights are different due to the difference in the identified buildings where the instruments were deployed. In Marconi St., ground level (site MA1, 44 • Figure 4).
Two mobile laboratories, i.e., vans equipped for air quality and meteorological measurements, were deployed by ARPAE at MA1 and LB1 locations. Mobile laboratories were equipped for continuous measurements of air quality pollutants such as nitrogen dioxides (NO x -NO 2 -NO), carbon monoxide (CO), sulphur dioxide (SO 2 ), ozone (O 3 ), and particulate matter (PM 10 and PM 2.5 ), together with cup anemometer and wind vane for wind speed and direction measurements and thermohygrometers for temperature and relative humidity measurements, all at a sampling rate of 1 min except for particulate matter (1 day sampling rate). High-frequency meteorological instrumentation was installed at each level in and above the canyons to enhance the temporal resolution of the meteorological data and to enable the investigation of turbulent processes. The 3D-wind field was sampled at a 20 Hz with GILL Windmaster sonic anemometers (Gill Instruments Limited, Hampshire, UK), while air temperature, relative humidity and pressure were gathered at sampling rate of 1 Hz respectively by HCS2S3 Rotronic thermohygrometers (Rotronic Instruments Ltd., Crawley, UK) and Vaisala PTB110 barometers (Vaisala, Helsinki, Finland). Additional measurements of the energy balance between short-wave and long-wave far-infrared radiation versus surface-reflected short-wave and outgoing long-wave radiation were performed at sampling rate of 1 min at both MA3 and LB3 locations with CNR4 radiometers (Kipp & Zonen B.V., Delft, The Netherlands). A summary of the experimental setup at each location site is presented in Table 1. summer campaign to the 15 m of the winter one according to site availability (the locations of each site is displayed in Figure 4). Two mobile laboratories, i.e., vans equipped for air quality and meteorological measurements, were deployed by ARPAE at MA1 and LB1 locations. Mobile laboratories were equipped for continuous measurements of air quality pollutants such as nitrogen dioxides (NOx-NO2-NO), carbon monoxide (CO), sulphur dioxide (SO2), ozone (O3), and particulate matter (PM10 and PM2.5), together with cup anemometer and wind vane for wind speed and direction measurements and thermohygrometers for temperature and relative humidity measurements, all at a sampling rate of 1 min except for particulate matter (1 day sampling rate). High-frequency meteorological instrumentation was installed at each level in and above the canyons to enhance the temporal resolution of the meteorological data and to enable the investigation of turbulent processes. The 3Dwind field was sampled at a 20 Hz with GILL Windmaster sonic anemometers (Gill Instruments Limited, Hampshire, UK), while air temperature, relative humidity and pressure were gathered at sampling rate of 1 Hz respectively by HCS2S3 Rotronic thermohygrometers (Rotronic Instruments Ltd., Crawley, UK) and Vaisala PTB110 barometers (Vaisala, Helsinki, Finland). Additional measurements of the energy balance between short-wave and long-wave far-infrared radiation versus surface-reflected short-wave and outgoing long-wave radiation were performed at sampling rate of 1 min at both MA3 and LB3 locations with CNR4 radiometers (Kipp & Zonen B.V., Delft, The Netherlands). A summary of the experimental setup at each location site is presented in Table 1.        sideline to the urban ring-road next to a small vegetated area. Within the experimental campaigns, two intensive operational periods (one per campaign) were devoted to collecting extra measurements to investigate the UHI phenomenon. Specifically, the temperature distribution of building façades and ground surfaces in the two streets was analyzed collecting simultaneously images in the two canyons with two high-performance thermal cameras FLIR T620 ThermalCAMs (FLIR Systems Inc., Wilsonville, OR, USA) with uncooled microbolometer, 640 × 480 pixels resolution and an image acquisition frequency of 50-60 Hz, selecting two appropriate periods characterized by clear sky and calm winds, i.e., during synoptic conditions particularly favorable to the onset of UHI, i.e., characterized by weak synoptic forcing [54] (22-23 August 2017 during the summer campaign, and 8-9 February 2018 during the winter campaign). Specifically, we selected days for which the weather forecasts indicated wind speed below 5 m s −1 at 700 hPa and characterized by clear-sky conditions (see Appendix B). During the experiments, images were collected simultaneously by standing operators during a 24 h acquisition with regular intervals of 2 h (total of 12 acquisitions in 24 h at each site). This choice allowed the collection of images at 12:00 (close to the maximum surface temperature), 14:00, 16:00 (close to the maximum air temperature), 18:00, 20:00, 22:00 (close to maximum UHI intensity), 24:00, 02:00, 04:00, 06:00, 08:00, 10:00, and 12:00 (UTC + 2). Note that during the summer campaign a technical issue occurred to the thermal camera in Laura Bassi St., so measurements in that canyon from midnight to the end of the experiment were performed with a delay of 1 h with respect to Marconi St. Analyzed buildings were selected on the basis of the homogeneity of construction material and the absence of obstacles (balconies, eave, etc.), metal or glass. Several shots or portions of the same building façade at several different heights were taken in order to maintain a similar resolution for all images. Finally, measurements of ground surfaces were also collected. The precise positioning of the selected buildings is provided in Figure 4. In particular, in Marconi St. the buildings 1, 2, and 8 are located on the east side, while buildings number 4, 5, and 6 are located on the west side; finally, the positions 3 and 7 indicate the sites where the temperature of the asphalt surface was measured. Similarly, in Laura Bassi St., buildings 5, 7, and 8 are located on the east side, while buildings number 1, 2, and 3 are located on the west side, and finally numbers 4 and 6 indicate the positions of the measurements of the ground surface.

Urban Heat Island
To address the urban heat island development at neighborhood and city scale levels we used the data collected within the two intensive thermographic campaigns performed in the two street canyon areas described previously. Figure 5 presents and compares the diurnal evolution of the air temperature measured within and above the two street canyons of Marconi and Laura Bassi (sites 1 and 2 respectively) and at the rural meteorological station of Mezzolara (MZ in Figure 5 Figure 5 shows very clearly the UHI phenomenon, i.e., the effect of overheating of cities with respect to the countryside during the summer season. The effect is reduced at Laura Bassi St. with respect to Marconi St., both due to the presence of vegetation as well as being located further from the city center in a residential area in the outskirts of Bologna where the presence of small houses and apartments favor heat exchange and urban thermal comfort with respect to the highly packaged buildings typical of the historical city center where Marconi St. lies. This also indicates the positive role of trees in improving thermal comfort and reducing the urban heat island phenomenon. In addition, the plots also show that during the night both street canyons are isothermal, i.e., there is no temperature difference between the building façades on the two sides of the streets nor with the temperature of the air in the street canyon, while during the day, depending on the position of the Sun, one side is hotter than the other one: this effect is visible at both street canyons but is more evident at Marconi St. During winter (Figure 5b), similar UHI effects are observed at both sites, with  Figure 5 shows very clearly the UHI phenomenon, i.e., the effect of overheating of cities with respect to the countryside during the summer season. The effect is reduced at Laura Bassi St. with respect to Marconi St., both due to the presence of vegetation as well as being located further from the city center in a residential area in the outskirts of Bologna where the presence of small houses and apartments favor heat exchange and urban thermal comfort with respect to the highly packaged buildings typical of the historical city center where Marconi St. lies. This also indicates the positive role of trees in improving thermal comfort and reducing the urban heat island phenomenon. In addition, the plots also show that during the night both street canyons are isothermal, i.e., there is no temperature difference between the building façades on the two sides of the streets nor with the temperature of the air in the street canyon, while during the day, depending on the position of the Sun, one side is hotter than the other one: this effect is visible at both street canyons but is more evident at Marconi St.
During winter (Figure 5b), similar UHI effects are observed at both sites, with urban-rural differences of about 6 • C but with a nighttime increase in both air and buildings' temperatures observed during the night of 8 February 2018. This effect was partially attributed to the presence of a cloud cover limiting the nighttime radiative cooling, as indicated during night both by local observations and by thermal composite satellite images. In addition, the nighttime increase can be tentatively attributed to the on and off switching of residential heating: in this case, the difference between the patterns observed in the two street canyons might be due to the different kinds of buildings in the two street canyons, where Marconi St. is characterized by block of flats with central heating whereas Laura Bassi St. is instead characterized by small residential houses with independent heating. Considering the similar urban-rural temperature differences in the two seasons and the reduced differences at neighborhood scales observed in the winter season, in addition to the limited impact of vegetation during this season, on the following we chose to focus just on the summer observations.

City (Urban-Rural) and Neighborhood (Canyon-Canyon) Scales
The UHI is generally attributed to the differential heat release in the atmosphere between urban and rural environment, causing a horizontal temperature difference between them. However, the heterogeneous morphology of a city allows the development of horizontal temperature differences within the urban environment (for example between two neighborhoods characterized by different morphology and vegetation amounts). In the current study, we will address the classical urban-rural temperature difference as a city-scale UHI, while the intra-city (canyon-canyon) temperature difference will be interpreted as neighborhood-scale UHI. Figure 6a shows the evolution of the normalized UHI intensities ∆T + satisfying this classification. By "normalized" we intend the measured value of the temperature difference divided by its maximum in the analyzed period. We use the normalized values to enhance the shape of the quantity's evolution, highlighting the analogous evolution of the ∆T + signals notwithstanding the different spatial scales involved in the computation. It is also worth noticing that while an inversion in the UHI intensity is observed computing the temperature difference with above-canopy urban temperatures, this behavior almost disappears with the in-canyons and between-canyons temperatures, highlighting again the high impact of the morphology and heat release from urban surfaces (buildings and street).
Atmosphere 2020, 11, x FOR PEER REVIEW 13 of 32 urban-rural differences of about 6 °C but with a nighttime increase in both air and buildings' temperatures observed during the night of 8 February 2018. This effect was partially attributed to the presence of a cloud cover limiting the nighttime radiative cooling, as indicated during night both by local observations and by thermal composite satellite images. In addition, the nighttime increase can be tentatively attributed to the on and off switching of residential heating: in this case, the difference between the patterns observed in the two street canyons might be due to the different kinds of buildings in the two street canyons, where Marconi St. is characterized by block of flats with central heating whereas Laura Bassi St. is instead characterized by small residential houses with independent heating. Considering the similar urban-rural temperature differences in the two seasons and the reduced differences at neighborhood scales observed in the winter season, in addition to the limited impact of vegetation during this season, on the following we chose to focus just on the summer observations.

City (Urban-Rural) and Neighborhood (Canyon-Canyon) Scales
The UHI is generally attributed to the differential heat release in the atmosphere between urban and rural environment, causing a horizontal temperature difference between them. However, the heterogeneous morphology of a city allows the development of horizontal temperature differences within the urban environment (for example between two neighborhoods characterized by different morphology and vegetation amounts). In the current study, we will address the classical urban-rural temperature difference as a city-scale UHI, while the intra-city (canyon-canyon) temperature difference will be interpreted as neighborhood-scale UHI. Figure 6a shows the evolution of the normalized UHI intensities ∆ satisfying this classification. By "normalized" we intend the measured value of the temperature difference divided by its maximum in the analyzed period. We use the normalized values to enhance the shape of the quantity's evolution, highlighting the analogous evolution of the ∆ signals notwithstanding the different spatial scales involved in the computation. It is also worth noticing that while an inversion in the UHI intensity is observed computing the temperature difference with above-canopy urban temperatures, this behavior almost disappears with the in-canyons and between-canyons temperatures, highlighting again the high impact of the morphology and heat release from urban surfaces (buildings and street). Refs. [55][56][57] have derived and proven a simple relation between the urban wind speed and the UHI intensity directly from the temperature equation Refs. [55][56][57] have derived and proven a simple relation between the urban wind speed and the UHI intensity directly from the temperature equation The advection term (first term of the right-hand side of the equation, hereafter RHS) is expressed with its scale values, where U is the urban wind speed, T u is the urban temperature, T r is the rural temperature and L = √ 4A/π = 13,400 m is the city diameter. This term is responsible for the dissipation of the UHI due to the temperature advection and so can be interpreted as a cooling factor. The last term of the RHS of Equation (1)  H 0 , and sustains the UHI slowing the cooling. As this last term is not always measured, ref. [56] parametrized its contribution in a simple attenuation factor K of the cooling term. Therefore, integrating Equation (1) considering only the first term of the RHS multiplied by K and using ∆T = T u − T r , we get where ∆T 0 is the initial temperature difference before the onset of the cooling, i.e., is the maximum temperature difference. The factor K ranges between 0 and 1. When K = 1 we are in the limiting case of .
H 0 = 0, while for K < 1 the diffusive heat flux retarding the night cooling process is taken into account. Figure 6b shows the dependency of the city and neighborhood-scales UHI intensity on velocity, with the solid lines retrieved using Equation (2). ∆T 0 and K are extrapolated at U = 0.5 ms −1 (activation threshold for cup and vane anemometers) fitting Equation (2) on the measured data for a fixed time t f = 6 h identified as the time of cooling between sunset and the nocturnal ∆T equilibrium (see Table 2). Table 2. Fit characteristics and likelihood parameters (coefficient of determination R 2 and root mean square error RMSE) between data in Figure 6b and Equation (2) Figure 6b clearly shows the UHI-intensity dependency on the urban wind speed, as suggested by Equation (2). Moreover, it highlights the similar behavior of the city and neighborhood-scale UHI. The different degree of dependency observed between the two scales of UHI can attributed to the smaller initial (and maximum) temperature difference ∆T 0 measured at the neighborhood scale. Indeed, this evidence is quite expected as the temperature gradients within the same city should be constrained by the temperature diffusion within the urban environment. Conversely, the largely different land use of the countryside is expected to enlarge the temperature difference with the city, explaining the larger initial temperature observed in Figure 6b. The values of K also highlight the larger contribution of the heat flux divergence in the canyon-rural temperature cooling process, as the canopy environment tends to preserve the diurnal temperature more easily than the atmosphere above (leading to a smaller contribution of the heat flux divergence in the urban-rural temperature cooling process) due to the vicinity of multiple heat sources. The presence of vegetation in Laura Bassi St. does not impact on the heat fluxes contribution, but together with the larger aspect ratio of the canyon (with respect to Marconi St.), contributes to reduce the UHI intensity. This difference between the business-center canyon of Marconi St. and the vegetated-residential canyon of Laura Bassi St. is certified by the presence of a small UHI phenomenon between the two neighborhoods, due to the different cooling efficacy. The impact of the heat flux divergence is instead small as already observed for the city-scale UHI involving the canyons.

Canyon (Facade-Facade) Scale: The In-Canyon Thermal Circulation
As by the mere definition of temperature difference between two locations, a very local UHI effect can be observed also within an urban canopy, where a horizontal temperature gradient can establish between the opposite sides of the canyon (building façade to building façade). More precisely, this temperature gradient acts as the primary forcing of the in-canyon thermal circulation. As for the large scale UHI phenomenon previously discussed, this local phenomenon establishes a relationship between the in-canyon wind speed and the building façade temperature difference (evaluated using the east and west canyon-side building façades shown in Figure 5). Figure 7a shows the similar behavior of these two quantities, highlighting that a direct proportionality relates the temperature difference and the wind speed (while it was inverse at larger scales). While at larger scales the wind speed was a cooling factor, within the canopy the temperature difference forces the circulation. The following quadratic fit suits this in-canyon dependency between the building façade temperature difference from the thermographic data and the in-canyon wind speed measured at the ground level of each canyon as displayed in Figure 7b. Fit parameters are reported in Table 3. Despite the clear dependency of the wind speed on the temperature difference, and hereby the importance of the thermal forcing in driving the in-canyon circulation, the presence of a non-negligible ambient wind speed U a measured above the canopy (as observed at the site SI and AS, Figure 6b) can be responsible to drive an additional inertial component within the canyon. In this scenario, the measured wind speed would be the sum of a thermally driven U r and an inertially driven U I component. While the inertial component is typically a fraction of the wind speed above the canopy, the thermal one requires a different parametrization. As by the mere definition of temperature difference between two locations, a very local UHI effect can be observed also within an urban canopy, where a horizontal temperature gradient can establish between the opposite sides of the canyon (building façade to building façade). More precisely, this temperature gradient acts as the primary forcing of the in-canyon thermal circulation. As for the large scale UHI phenomenon previously discussed, this local phenomenon establishes a relationship between the in-canyon wind speed and the building façade temperature difference (evaluated using the east and west canyon-side building façades shown in Figure 5). Figure 7a shows the similar behavior of these two quantities, highlighting that a direct proportionality relates the temperature difference and the wind speed (while it was inverse at larger scales). While at larger scales the wind speed was a cooling factor, within the canopy the temperature difference forces the circulation. The following quadratic fit suits this in-canyon dependency between the building façade temperature difference from the thermographic data and the in-canyon wind speed measured at the ground level of each canyon as displayed in Figure 7b. Fit parameters are reported in Table 3. Despite the clear dependency of the wind speed on the temperature difference, and hereby the importance of the thermal forcing in driving the in-canyon circulation, the presence of a non-negligible ambient wind speed measured above the canopy (as observed at the site SI and AS, Figure 6b) can be responsible to drive an additional inertial component within the canyon. In this scenario, the measured wind speed would be the sum of a thermally driven and an inertially driven component. While the inertial component is typically a fraction of the wind speed above the canopy, the thermal one requires a different parametrization.  Table 3. Fit characteristics and likelihood parameters (coefficient of determination and root mean square error (RMSE)) between data in Figure 7b Table 3. Fit characteristics and likelihood parameters (coefficient of determination R 2 and root mean square error (RMSE)) between data in Figure 7b and Equation (3)  As the thermally driven component must be a quadratic function of the temperature difference within the canyon (we must satisfy Equation (3) and the inertial component cannot by definition), we adopt the UHI convergence velocity [58,59] defined as U r = C(gαL|∆T|) 1/2 (4) where g is the gravitational acceleration, α is the thermal expansion coefficient, L is a length scale set equal to the canyon width, and C = 0.08 is a constant computed by [58] from laboratory experiments and the similarity theory. Specifically, C is given by the independency of U r of the Reynolds number Re = UW/ν (where U in-canyon ground-level wind speed, W is the mean canyon width, and ν is the kinematic viscosity of the air) evaluated comparing Equation (4) with Re, as done in Figure 8a. In our case, scaling the measured wind speed and evaluating its independency on the Reynolds number, we obtain different values of C, named C 1 , equal to 0.13 and 0.32 for Laura Bassi St. and Marconi St., respectively. This discrepancy in the constant between the laboratory and the current investigation is associated with the inertial component of the wind speed, not accounted by Equation (4). To verify this hypothesis, we compute the total parametrized velocity U s as where U a is the wind speed measured in Silvani Street. Equation (5) provides a good approximation of the in-canyon measured velocity, as shown in Figure 8a. The clear difference between the two canyons unveils the total absence of a diurnal evolution of the wind speed within Laura Bassi St., while in Marconi St. the behavior approaches the ambient wind dynamics. It is suggested that the presence of trees is such that the wind in the surface layer of the canopy (below the tree crown) is forced to remain small during the whole day, causing a decrease in the circulation intensity. Moreover, the presence of trees is a larger reduction factor for the inertial wind-speed component than the disadvantageous aspect ratio of Marconi St.
Atmosphere 2020, 11, x FOR PEER REVIEW 16 of 32 As the thermally driven component must be a quadratic function of the temperature difference within the canyon (we must satisfy Equation (3) and the inertial component cannot by definition), we adopt the UHI convergence velocity [58,59] defined as where is the gravitational acceleration, is the thermal expansion coefficient, is a length scale set equal to the canyon width, and = 0.08 is a constant computed by [58] from laboratory experiments and the similarity theory. Specifically, is given by the independency of of the Reynolds number = / (where in-canyon ground-level wind speed, is the mean canyon width, and is the kinematic viscosity of the air) evaluated comparing Equation (4) with , as done in Figure 8a. In our case, scaling the measured wind speed and evaluating its independency on the Reynolds number, we obtain different values of , named , equal to 0.13 and 0.32 for Laura Bassi St. and Marconi St., respectively. This discrepancy in the constant between the laboratory and the current investigation is associated with the inertial component of the wind speed, not accounted by Equation (4). To verify this hypothesis, we compute the total parametrized velocity as = 0.08( |∆ |) / + (5) where is the wind speed measured in Silvani Street. Equation (5) provides a good approximation of the in-canyon measured velocity, as shown in Figure 8a. The clear difference between the two canyons unveils the total absence of a diurnal evolution of the wind speed within Laura Bassi St., while in Marconi St. the behavior approaches the ambient wind dynamics. It is suggested that the presence of trees is such that the wind in the surface layer of the canopy (below the tree crown) is forced to remain small during the whole day, causing a decrease in the circulation intensity. Moreover, the presence of trees is a larger reduction factor for the inertial wind-speed component than the disadvantageous aspect ratio of Marconi St. The different behavior of UHIs depending on the horizontal scale suggests a complex interaction between different areas of the same city, where different internal boundary layers (IBLs) may develop and interact with each other. Specifically, [60] proposed a conceptual model where the presence of "overlapping IBLs" is such that nearer the urban canopy there are spatially localized patches of flow in local equilibrium with the surface. This model suggested that atmospheric phenomena may be influenced by being confined in the same IBL or developed between different IBLs. Homogeneous conditions (as an isothermal atmosphere within the canopy, see Figure 5) defines a single IBL, which may develop gradients with the neighboring IBLs, leading to spatial heterogeneities. For this reason, the city and neighborhood scale UHIs are expected to behave similarly as they involve multiple IBLs (in the case of the city-scale UHI, the whole urban and rural boundary layers). Conversely, the façade- The different behavior of UHIs depending on the horizontal scale suggests a complex interaction between different areas of the same city, where different internal boundary layers (IBLs) may develop and interact with each other. Specifically, [60] proposed a conceptual model where the presence of "overlapping IBLs" is such that nearer the urban canopy there are spatially localized patches of flow in local equilibrium with the surface. This model suggested that atmospheric phenomena may be influenced by being confined in the same IBL or developed between different IBLs. Homogeneous conditions (as an isothermal atmosphere within the canopy, see Figure 5) defines a single IBL, which may develop gradients with the neighboring IBLs, leading to spatial heterogeneities. For this reason, the city and neighborhood scale UHIs are expected to behave similarly as they involve multiple IBLs (in the case of the city-scale UHI, the whole urban and rural boundary layers). Conversely, the façade-to-façade temperature difference occurs within a single IBL and, as such, it presents a different and unique interaction with the local flow field.

Canyon-Scale UHI and UPI Interaction
It has been documented that UHI and UPI, two major problems of the urban environment which are becoming more serious with rapid urbanization, can interact with each other, through their effect on the intensity of turbulence mixing and on radiation transfer processes [16]. This makes it necessary to study both effects concurrently. To establish a link between the thermal circulation and the pollutant concentrations measured in the Marconi St., we derived a parameterization for both quantities. In particular, the UHI convergence velocity U r toward an urban area was parameterized following the previously introduced Equation (4) using the value of 0.08 for the empirical constant C = 0.08 previously derived by [58] in laboratory experiments analyzing the horizontal radial convergence flow velocity in the presence of heated elements in a fluid. Implementing this parameterization for a street canyon, we considered the characteristic length scale of the UHI equal to the mean street canyon width, while the temperature difference causing the thermal circulation was considered that between the west and east building façades.
In the case of pollutant concentrations, we derived a normalization considering all the processes impacting on pollutant concentrations in a canyon, clearly depending first of all on the emission source rate but impacted by other atmospheric local processes. To derive a normalization for pollutant concentrations in an urban street canyon, we selected CO as a passive pollutant emitted by vehicle exhausts. As previously presented in [61], we then identified the key influencing factors on pollutant concentrations in an urban street canyon: (1) the source strength; (2) the canyon geometry, related to the dispersion volume and affecting the mass exchange between the canyon and the upper atmosphere; (3) the background wind velocity impinging on the canyon. A detailed explanation of the algorithm will be presented in [62], but briefly in analogy with [63], measured CO concentrations were normalized with a reference concentration CO* defined as CO * = Q e HWU b (6) where Q e is the pollutant source rate [g s −1 ], U b is the bulk velocity measured above the canyon, H is the mean building height, and W is the mean street canyon width. In particular, bulk velocity was derived from the measurements of the sonic anemometers placed on the rooftops in the two street canyons, while the mean building height and street canyon width were estimated from detailed digital maps of the two areas. An algorithm was instead developed to estimate the pollutant source rate. To this end, the traffic counts available as number of vehicles travelling in Marconi St. with a 5 min time resolution from inducive loops technology were classified upon fleet composition, i.e., fuel type (gasoline, diesel, LPG -Liquid Petroleum Gas-, CNG -Compressed Natural Gas-), vehicle type (motorcycle, light vehicles, heavy vehicles, buses) and EURO technology. The number of buses passing through Marconi St. was calculated from the bus time schedules in Bologna (Regional transport company TPER (Trasporto Passeggeri Emilia Romagna) [64]). The local fleet composition was then extracted from the regional inventory of circulating vehicles (Italian Car Club company (ACI, Automobile Club d'Italia) [65]) and was used to disaggregate the difference between the total traffic counts and the number of buses derived from bus schedules into traffic counts per vehicle type. Pollutant emission rates [g/km] were then finally estimated using the joint EMEP/EEA (European Monitoring and Evaluation Programme/European Environment Agency) air pollutant emission inventory guidebook for each vehicle category [66]. Finally, the pollutant source rate Q e was estimated from pollutant emissions using a representative vehicle urban speed of 19 km h −1 [66].  Figure 9 shows very clearly the existence of a relationship between the UHI velocity and the pattern of normalized concentrations in Marconi St., whereas in Laura Bassi the thermal circulation is not directly related with the pollutant distribution, which is driven by other mechanisms owing to the presence of trees and to the different morphology of the neighborhood. It should be noted that although the good relationship between Ur and the normalized concentration may appear counter intuitive, it is a consequence of the fact that the presence of a temperature gradient between two opposite façades may lead to the establishment of an opposite vortex [67], which confines the primary vortex closest to the ground and diminish the vertical exchange. The establishment of this new vortex depends on the strength of the temperature gradient. If not so strong and of the order of 10 °C as in our case, we will see that the effect of the temperature gradient will be to deform the vortex and finally diminish its ability to remove pollutants from the near surface. Simulations performed in the following section will be used to prove this argument. The rise in normalized concentration observed at Laura Bassi St. at 3:00 (local time) is due to the normalization procedure that uses a reference concentration CO* (Equation (6)). Further, while traffic pollutant concentrations (CO, but also NOx) show generally lower nighttime values and a first peak during the morning rush hours, the normalization procedure involves calculating the ratio with the reference CO* estimated with the pollutant source rate and the bulk velocity (the other values in the equation are linked to geometric factors and as such are constant values). Average conditions at 3:00 were characterized by high bulk velocities and low pollutant source rates (as connected with the reduced transit of vehicles during nighttime), altogether resulting in a very low CO* value and therefore a sharp increase of the normalized CO/CO* concentration. The scatterplot (Figure 10) of normalized CO concentrations vs. UHI velocity in Marconi St. confirms the presence of a linear relationship between the two quantities,  Figure 9 shows very clearly the existence of a relationship between the UHI velocity and the pattern of normalized concentrations in Marconi St., whereas in Laura Bassi the thermal circulation is not directly related with the pollutant distribution, which is driven by other mechanisms owing to the presence of trees and to the different morphology of the neighborhood. It should be noted that although the good relationship between Ur and the normalized concentration may appear counter intuitive, it is a consequence of the fact that the presence of a temperature gradient between two opposite façades may lead to the establishment of an opposite vortex [67], which confines the primary vortex closest to the ground and diminish the vertical exchange. The establishment of this new vortex depends on the strength of the temperature gradient. If not so strong and of the order of 10 • C as in our case, we will see that the effect of the temperature gradient will be to deform the vortex and finally diminish its ability to remove pollutants from the near surface. Simulations performed in the following section will be used to prove this argument. The rise in normalized concentration observed at Laura Bassi St. at 3:00 (local time) is due to the normalization procedure that uses a reference concentration CO* (Equation (6)). Further, while traffic pollutant concentrations (CO, but also NO x ) show generally lower nighttime values and a first peak during the morning rush hours, the normalization procedure involves calculating the ratio with the reference CO* estimated with the pollutant source rate and the bulk velocity (the other values in the equation are linked to geometric factors and as such are constant values). Average conditions at 3:00 were characterized by high bulk velocities and low pollutant source rates (as connected with the reduced transit of vehicles during nighttime), altogether resulting in a very low CO* value and therefore a sharp increase of the normalized CO/CO* concentration. The scatterplot (Figure 10) of normalized CO concentrations vs. UHI velocity in Marconi St. confirms the presence of a linear relationship between the two quantities, suggesting that the thermal circulation inside the street canyon may increase the pollutant concentrations at pedestrian level, a result in agreement with previous modeling studies by [68,69] and with the nighttime negative correlations between UPI and UHI obtained by [16]. Previous literature [70] showed the reverse effect, i.e., that haze pollution may contribute further to the UHI, confirming again the presence of an interaction between air pollution and the UHI phenomenon. The reduced magnitude of normalized concentrations in Marconi St. during the summer thermographic campaign is due to the fact that our normalization procedure successfully removes all other factors driving pollutant concentrations in the street canyon, i.e., the geometry of the canyon, the pollutant emissions and the bulk wind speed above the canyon. Therefore, our results show the existence of a relation between the thermal forcing exerted by UHI and pollutant concentrations in a street canyon. In the following section, the scaling for heat and concentration and the existence of a linear relationship between concentrations and the thermal circulation is verified with CFD simulations adequately set up.
Atmosphere 2020, 11, x FOR PEER REVIEW 19 of 32 confirming again the presence of an interaction between air pollution and the UHI phenomenon. The reduced magnitude of normalized concentrations in Marconi St. during the summer thermographic campaign is due to the fact that our normalization procedure successfully removes all other factors driving pollutant concentrations in the street canyon, i.e., the geometry of the canyon, the pollutant emissions and the bulk wind speed above the canyon. Therefore, our results show the existence of a relation between the thermal forcing exerted by UHI and pollutant concentrations in a street canyon. In the following section, the scaling for heat and concentration and the existence of a linear relationship between concentrations and the thermal circulation is verified with CFD simulations adequately set up.

Verification with CFD Simulations
In order to examine the results presented in the previous section, high resolution numerical CFD (Computational Fluid Dynamics) were carried out on an idealized street canyon with the same geometry as Marconi St. and considering the effect of temperature difference between building facades on the two sides of the street canyon.
The CFD code CD-adapco STAR-CCM+ 19.2 has been employed to solve the steady-state Reynolds-averaged Navier-Stokes (RANS) equations with a realizable k-ε turbulence model, similar to [71]. The buoyancy thermal effects have been considered in this work, as generated by the temperature differences between the building facades and air temperature. The effects related to radiation and material properties have been included in the temperature boundary conditions, that in our case were known from measurements. A detailed description of the model is available in the Appendix A.

Verification with CFD Simulations
In order to examine the results presented in the previous section, high resolution numerical CFD (Computational Fluid Dynamics) were carried out on an idealized street canyon with the same geometry as Marconi St. and considering the effect of temperature difference between building facades on the two sides of the street canyon.
The CFD code CD-adapco STAR-CCM+ 19.2 has been employed to solve the steady-state Reynolds-averaged Navier-Stokes (RANS) equations with a realizable k-ε turbulence model, similar to [71]. The buoyancy thermal effects have been considered in this work, as generated by the temperature differences between the building facades and air temperature. The effects related to radiation and material properties have been included in the temperature boundary conditions, that in our case were known from measurements. A detailed description of the model is available in the Appendix A.
The inflow wind profile has been calculated solving a two equations system using data from the two previously cited meteorological stations (Asinelli and Silvani), in order to obtain the friction velocity u * and the roughness length z 0 for each wind condition. The displacement height d has been set equal to 1/3 of the street canyon medium height, i.e., 11 m for Marconi St. This approach highlights the relationship between wind direction and urban geometry.
The temperature boundary conditions were instead set using the street canyon façade temperatures measured during the summer thermographic campaigns, thoroughly described in the previous section. In particular, the differences between the following two cases are considered: 22/08/2018 12:00 (UTC+2) 2.
22/08/2018 16:00 (UTC+2) The velocity and temperature boundary conditions for these two cases are reported in Table 4. The two cases shown in Table 4 are particularly interesting since the wind velocity magnitude is about the same, the wind direction is perpendicular to the canyon (coming from east) in both cases, but the temperature difference presents the opposite sign. Therefore, these two cases represent two opposite scenarios in terms of thermal boundary conditions. The numerical simulations are aimed to support the analysis and the results about the interaction between the UHI and the UPI obtained at street-canyon from the observations collected in the experimental campaign in Marconi St., presented in the previous section. For this reason, since the simulations refer to Marconi street where no tree is present, the role of trees is not considered in the setup of CFD simulations. Figures 11 and 12 show the velocity vectors and CO concentration obtained in three street canyons with different aspect ratios, with case 1 and case 2 boundary conditions. Considering a canyon with AR = 1 and comparing the two cases with opposite building temperature differences (Figures 11 and 12a), we observe that the location of the vortex center slightly moves lower on the right. This explains the reduced CO concentration differences obtained in the two scenarios, for both the leeward side and the windward side.
Considering a canyon with AR equal to 1.65 and comparing the results obtained in the two cases with opposite temperature difference (Figures 11 and 12b), we observe that the buoyancy-induced flow breaks the strong recirculating vortex. The wind flow structure within the canyon is very peculiar, still resulting in only one primary vortex but located far away from the ground. Indeed, an upward buoyancy flux opposes the downward airflow along the wall. As a result, the CO pollutant plume increases significantly its shape, above all in lower areas and especially along the windward side: overall, the pollutant dispersion is greatly reduced.  Considering a canyon with AR = 1 and comparing the two cases with opposite building temperature differences (Figures 11 and 12a), we observe that the location of the vortex center slightly moves lower on the right. This explains the reduced CO concentration differences obtained in the two scenarios, for both the leeward side and the windward side.
Considering a canyon with AR equal to 1.65 and comparing the results obtained in the two cases with opposite temperature difference (Figures 11 and 12b), we observe that the buoyancy-induced flow breaks the strong recirculating vortex. The wind flow structure within the canyon is very peculiar, still resulting in only one primary vortex but located far away from the ground. Indeed, an upward buoyancy flux opposes the downward airflow along the wall. As a result, the CO pollutant plume increases significantly its shape, above all in lower areas and especially along the windward side: overall, the pollutant dispersion is greatly reduced.
Finally, considering a canyon with AR equal to 2 and observing the results obtained in the two cases with opposite temperature difference (Figures 11 and 12c), we can remark that the combination of thermal effects and buildings' narrow shape creates a multiple vortex structure. When the windward wall and the ground are warmer than the air and the leeward-side, a weaker ground-level secondary recirculating vortex is easily recognizable near the windward side, right below the primary one. Indeed, the upward buoyancy flux strongly contrasts with the downward airflow along the wall, dividing the wind structure into two counter-rotating vortices, and in particular a clockwise top vortex and a reverse secondary vortex. In this scenario, the pollutant dispersion is enhanced with respect to the preceding scenario considering AR equal to 1.65, but still significantly reduced and worrying if compared to a regular canyon as in the first scenario, especially close to the ground along pedestrians' walkways.
By considering the whole diurnal cycle for the cases AR equal to 1 and 2, respectively, a similar relation between Ur and CO/CO* as that derived in the previous section for AR equal to 1.65 and illustrated by Figures 9 and 10 is obtained. Figure 13 shows the pattern of normalized CO Finally, considering a canyon with AR equal to 2 and observing the results obtained in the two cases with opposite temperature difference (Figures 11 and 12c), we can remark that the combination of thermal effects and buildings' narrow shape creates a multiple vortex structure. When the windward wall and the ground are warmer than the air and the leeward-side, a weaker ground-level secondary recirculating vortex is easily recognizable near the windward side, right below the primary one. Indeed, the upward buoyancy flux strongly contrasts with the downward airflow along the wall, dividing the wind structure into two counter-rotating vortices, and in particular a clockwise top vortex and a reverse secondary vortex. In this scenario, the pollutant dispersion is enhanced with respect to the preceding scenario considering AR equal to 1.65, but still significantly reduced and worrying if compared to a regular canyon as in the first scenario, especially close to the ground along pedestrians' walkways.
By considering the whole diurnal cycle for the cases AR equal to 1 and 2, respectively, a similar relation between U r and CO/CO* as that derived in the previous section for AR equal to 1.65 and illustrated by Figures 9 and 10 is obtained. Figure 13 shows the pattern of normalized CO concentrations obtained in the CFD simulations, as compared to the observed normalized concentrations. Also shown are the plots for the vortices for five times. As such, Equation (4) provides a simple model to describe how the temperature difference between building facades within the street canyon can drive the airflow dynamics and the shape of the vortices within the canyon and therefore the pollutant concentration.
concentrations obtained in the CFD simulations, as compared to the observed normalized concentrations. Also shown are the plots for the vortices for five times. As such, Equation (4) provides a simple model to describe how the temperature difference between building facades within the street canyon can drive the airflow dynamics and the shape of the vortices within the canyon and therefore the pollutant concentration.

Conclusions
In this paper, we have investigated the multiple scale nature of the UHI phenomenon and its impact on flow dynamics, and have derived a relationship between UHI and pollutant concentrations in a real urban environment. To this aim, two intensive experimental field campaigns were carried out in two urban street canyons in the city of Bologna characterized by a different presence of vegetation, i.e., one being a tree-free street canyon and one characterized by the presence of a tree line along both sides of the road.
The comparison of the air temperatures measured in the city within two intensive thermographic campaigns with those measured in the countryside revealed a UHI effect of roughly 6-7 °C; however, the comparison between the temperatures measured in the two canyons indicated on average 2 °C higher temperatures in the tree-free street canyon with respect to the vegetated canyon located in a residential area away from the city center, as an indication of the presence of local UHIs inside the urban texture owing to the different morphologies and presence of vegetation of the neighborhoods. Different UHI scales have been identified and their impact on the flow velocity is quantified. City (urban-rural) and neighborhood (canyon-canyon) scales share the same exponential decrease of the UHI intensity as a function of wind speed, confirming the existence of an internal UHI induced by morphology differences and vegetation presence between different areas of the city. A third local UHI effect (induced by opposite façade temperature difference within the same canyon) is instead identified as direct forcing of a local thermal circulation, quantified as the UHI convergence velocity. The different behaviors of UHIs depending on the horizontal scale is hypothetically attributed to the development of overlapping IBLs within the urban one. As soon as the UHI intensity is associated to different IBLs, their behavior will be similar to the urban-rural UHI where the "IBLs"

Conclusions
In this paper, we have investigated the multiple scale nature of the UHI phenomenon and its impact on flow dynamics, and have derived a relationship between UHI and pollutant concentrations in a real urban environment. To this aim, two intensive experimental field campaigns were carried out in two urban street canyons in the city of Bologna characterized by a different presence of vegetation, i.e., one being a tree-free street canyon and one characterized by the presence of a tree line along both sides of the road.
The comparison of the air temperatures measured in the city within two intensive thermographic campaigns with those measured in the countryside revealed a UHI effect of roughly 6-7 • C; however, the comparison between the temperatures measured in the two canyons indicated on average 2 • C higher temperatures in the tree-free street canyon with respect to the vegetated canyon located in a residential area away from the city center, as an indication of the presence of local UHIs inside the urban texture owing to the different morphologies and presence of vegetation of the neighborhoods. Different UHI scales have been identified and their impact on the flow velocity is quantified. City (urban-rural) and neighborhood (canyon-canyon) scales share the same exponential decrease of the UHI intensity as a function of wind speed, confirming the existence of an internal UHI induced by morphology differences and vegetation presence between different areas of the city. A third local UHI effect (induced by opposite façade temperature difference within the same canyon) is instead identified as direct forcing of a local thermal circulation, quantified as the UHI convergence velocity. The different behaviors of UHIs depending on the horizontal scale is hypothetically attributed to the development of overlapping IBLs within the urban one. As soon as the UHI intensity is associated to different IBLs, their behavior will be similar to the urban-rural UHI where the "IBLs" are the urban and rural. Conversely, within the same IBL, temperature differences will most likely be a forcing for local thermal circulations.
Further, in order to examine the concurrent UHI and UPI phenomena, we derived the UHI convergence velocity resulting from the temperature difference between the west and east building façades in the non-vegetated canyon and a normalization for pollutant concentrations taking into account all key influencing factors, i.e., the emission source rate, canyon geometry, and background wind velocity. The analysis indicated the presence of a relationship between UHI and normalized concentrations, suggesting indeed an interaction between UHI and UPI phenomena, and in particular that the thermal circulation inside street canyons may be responsible of increasing the pollutant concentrations at pedestrian level.
The scaling adopted for the UHI convergence velocity was then verified by high resolution CFD simulations in an idealized street canyon considering two cases characterized by similar wind direction and intensity but opposite building temperature differences. CFD simulations demonstrated that the temperature difference between the building façades on the two sides of the canyons drives the airflow dynamics determining the position and shape of the vortices, definitely governing the pollutant concentrations in the street canyon.
Overall, the results of this work indicate that the positive role of mitigation solutions in reducing UHI effect may further extend to air pollution and therefore may be seen as an effective solution in mitigating two of the more pressing modern issues with important consequences on human health, therefore potentially greatly enhancing the sustainability of modern cities.
containing a box with elements with a medium size and another box surrounding the area and the street canyons, containing the finest elements.
Atmosphere 2020, 11, x FOR PEER REVIEW 26 of 32 containing a box with elements with a medium size and another box surrounding the area and the street canyons, containing the finest elements. Figure A1. Example of mesh zones (a) and grid refinement near walls (b).
A geometry with a canyon having the same dimensions of canyon A has been used as a base geometry for the first part of the sensitivity tests. Table A1 shows the dimensions of the box containing the finest elements, the size of the elements, and the total number of N_cells in the computational domain. The numerical results have been shown to verify and to support the analysis and the results about the interaction between UHI and UPI at street canyon scale, obtained from the observations collected in the experimental field campaign previously presented. Since the CFD simulations refer to a street canyon where no tree is present (Marconi St.), the setup does not consider the effect of trees.  Figure A2 shows the normalized root mean square deviation obtained from a vertical line in the middle of the canyon, defined as where x is the variable used for the comparison (velocity, CO concentration, or temperature), g1 and g2 are the grid number of two refinements, p is the point along a probe line where the comparison is made, and Np is the number of points in the probe line. A geometry with a canyon having the same dimensions of canyon A has been used as a base geometry for the first part of the sensitivity tests. Table A1 shows the dimensions of the box containing the finest elements, the size of the elements, and the total number of N_cells in the computational domain. The numerical results have been shown to verify and to support the analysis and the results about the interaction between UHI and UPI at street canyon scale, obtained from the observations collected in the experimental field campaign previously presented. Since the CFD simulations refer to a street canyon where no tree is present (Marconi St.), the setup does not consider the effect of trees. Five refinements have been compared, in the range N_cells = [0.250-6.181] millions of elements. Figure A2 shows the normalized root mean square deviation obtained from a vertical line in the middle of the canyon, defined as where x is the variable used for the comparison (velocity, CO concentration, or temperature), g 1 and g 2 are the grid number of two refinements, p is the point along a probe line where the comparison is made, and Np is the number of points in the probe line.
Atmosphere 2020, 11, x FOR PEER REVIEW 27 of 32 Figure A2. Convergence tests results. Figure A2 shows that the convergence is achieved for velocity, CO concentration, and temperature results.
A series of validation test has been performed by comparing the numerical results with experiments performed in Marconi street. Table A1 shows the comparison between experimental measurements and numerical results.
The deviation between the numerical results and measurements is then 1.4 °C, while the norm is 2.5 °C and the standard deviation is 0.7 °C.

Appendix B
This appendix provides a brief overview on the criteria we adopted to identify the atmospheric conditions suitable for the two intensive thermographic campaigns. After the identification of suitable periods from the regional synoptic charts provided from regional weather forecast service of the Regional Environmental Protection Agency (ARPAE), we analyzed the data recorded from the atmospheric soundings launched at San Pietro in Capofiume, located 20 km to the north-north-east of Bologna, to verify the absence of strong synoptic forcing and ensure the minimum criteria for the onset of the UHI with minimum synoptic perturbations. As from the soundings reported below, we were able to verify and check the fulfilment of our criteria in two consecutive days during both the summer and winter experimental campaign: wind speed at 700 hPa lower than 5 m s −1 and clear-sky conditions. The first sounding was meant to identify the day, while the second one to check the persistence of the required conditions within the experiment. Figures A3 and A4 show the soundings for the summer and winter UHI experiment, respectively.  Figure A2 shows that the convergence is achieved for velocity, CO concentration, and temperature results.
A series of validation test has been performed by comparing the numerical results with experiments performed in Marconi street. Table A1 shows the comparison between experimental measurements and numerical results.
The l ∞ deviation between the numerical results and measurements is then 1.4 • C, while the l 2 norm is 2.5 • C and the standard deviation is 0.7 • C.

Appendix B
This appendix provides a brief overview on the criteria we adopted to identify the atmospheric conditions suitable for the two intensive thermographic campaigns. After the identification of suitable periods from the regional synoptic charts provided from regional weather forecast service of the Regional Environmental Protection Agency (ARPAE), we analyzed the data recorded from the atmospheric soundings launched at San Pietro in Capofiume, located 20 km to the north-north-east of Bologna, to verify the absence of strong synoptic forcing and ensure the minimum criteria for the onset of the UHI with minimum synoptic perturbations. As from the soundings reported below, we were able to verify and check the fulfilment of our criteria in two consecutive days during both the summer and winter experimental campaign: wind speed at 700 hPa lower than 5 m s −1 and clear-sky conditions. The first sounding was meant to identify the day, while the second one to check the persistence of the required conditions within the experiment. Figures A3 and A4 show the soundings for the summer and winter UHI experiment, respectively.  Within both periods, clear-sky conditions are observed through the absence of convective available potential energy and small amounts of condense and gaseous water, ensuring the absence of both convective and stratus clouds as well as precipitation. Only during the second night of the winter period, the convective available potential energy was not zero but small enough to ensure that no convective events would have developed during the following day. Wind speed in the lower troposphere (below 700 hPa) was always weak, reaching the 5 m s −1 only during the first night of the summer period.
Therefore, both periods satisfied the criteria of weak synoptic forcing suitable for the best characterization of the UHI within the summer and winter intensive thermographic campaigns.  Within both periods, clear-sky conditions are observed through the absence of convective available potential energy and small amounts of condense and gaseous water, ensuring the absence of both convective and stratus clouds as well as precipitation. Only during the second night of the winter period, the convective available potential energy was not zero but small enough to ensure that no convective events would have developed during the following day. Wind speed in the lower troposphere (below 700 hPa) was always weak, reaching the 5 m s −1 only during the first night of the summer period.
Therefore, both periods satisfied the criteria of weak synoptic forcing suitable for the best characterization of the UHI within the summer and winter intensive thermographic campaigns. Within both periods, clear-sky conditions are observed through the absence of convective available potential energy and small amounts of condense and gaseous water, ensuring the absence of both convective and stratus clouds as well as precipitation. Only during the second night of the winter period, the convective available potential energy was not zero but small enough to ensure that no convective events would have developed during the following day. Wind speed in the lower troposphere (below 700 hPa) was always weak, reaching the 5 m s −1 only during the first night of the summer period.
Therefore, both periods satisfied the criteria of weak synoptic forcing suitable for the best characterization of the UHI within the summer and winter intensive thermographic campaigns.