Ventilation and Air Quality in City Blocks Using Large-Eddy Simulation — Urban Planning Perspective

Buildings and vegetation alter the wind and pollutant transport in urban environments. This comparative study investigates the role of orientation and shape of perimeter blocks on the dispersion and ventilation of traffic-related air pollutants, and the street-level concentrations along a planned city boulevard. A large-eddy simulation (LES) model PALM is employed over a highly detailed representation of the urban domain including street trees and forested areas. Air pollutants are represented by massless and passive particles (non-reactive gases), which are released with traffic-related emission rates. High-resolution simulations for four different city-block-structures are conducted over a 8.2 km2 domain under two contrasting inflow conditions with neutral and stable atmospheric stratification corresponding the general and wintry meteorological conditions. Variation in building height together with multiple cross streets along the boulevard improves ventilation, resulting in 7–9% lower mean concentrations at pedestrian level. The impact of smaller scale variability in building shape was negligible. Street trees further complicate the flow and dispersion. Notwithstanding the surface roughness, atmospheric stability controls the concentration levels with higher values under stably stratified inflow. Little traffic emissions are transported to courtyards. The results provide urban planners direct information to reduce air pollution by proper structural layout of perimeter blocks.


Introduction
Decreased air quality is one of the major environmental challenges urban areas are facing today.An increasing number of people are exposed to high air pollution levels due to the ongoing intense urbanisation [1].Exposure to air pollution has several acute and chronic health effects including respiratory and cardiovascular diseases that further increase mortality (e.g., [2,3]).To illustrate, in 2016 around 4.1 million premature deaths worldwide [4], of which almost 10% in Europe [5], were linked to exposure to elevated concentrations of ambient air pollution.
In urban areas, road traffic accounts for a significant share of the local air pollutant emissions (e.g., [6]).If the prevailing wind conditions within the street canyons lead to inefficient transport and mixing of air, the traffic emissions can accumulate and have longer residence times near the ground level.Consequently, the highest pollution levels within the street canyons are commonly observed at pedestrian level [7].In order to reduce human exposure to air pollutants within urban areas, the mechanisms affecting the dispersion conditions at street level must be understood and thereby examined in greater detail.To this end, the study of ventilation becomes of critical importance as it addresses the capacity with which a densely built urban structure is capable of replacing the contaminated air with ambient fresh air.Here, ventilation is recognised as a transport process that improves local air quality and closely relates to the term breathability [8,9].The efficiency at which street canyon ventilation occurs depends on the complex interaction between the atmospheric boundary layer (ABL) flow and the local urban structures.Recognising that the structural layout of urban landscape plays a critical role in determining the local urban air quality, opens up the possibility for the urban planners to incorporate air quality considerations into structural urban design.However, this necessitates sufficiently detailed modelling solutions to address the problem of street canyon ventilation at the scale of individual structures, such as perimeter blocks which are ubiquitous in European cities.
The current capacity to examine street canyon ventilation relies mainly upon parametric models, wind tunnel simulations and field measurements [7].Parametric models are computationally inexpensive but have decreased accuracy and their applicability is limited to very simplified model set-ups.Reduced-scale wind tunnel simulations, in turn, can suffer from similarity constraints, and field measurements from uncontrollable and unrepeatable boundary conditions and inability to account for spatial variability [10].Recently, along with rapidly advancing computer power, the application of computational fluid dynamics (CFD) using either models based on the Reynolds-averaged Navier-Stokes equations (RANS) with parametrised turbulence or the large-eddy simulation (LES) method has increased.CFD models provide complete and full-scale three-dimensional flow and concentration fields, which is a major advantage in real and complex urban areas.Majority of the previous urban CFD studies have employed the RANS method, owing to its lower computational costs.Nonetheless, LES outperforms above a complex urban surface due to its ability to resolve instantaneous turbulence structures (e.g., [11][12][13]) also when taking street trees into account [14].
Most of the urban LES studies have considered an idealistic two-dimensional street canyon or a simplified urban topography without including the aerodynamic impact of street trees (e.g., [15][16][17]).This simplification can be dangerous as porous trees decelerate the flow and generate turbulence (e.g., [18]) as well as influence the canyon vortex [19] that is commonly observed in idealised street canyons [20].Thereby, trees can modify the vertical transport of air pollutants within street canyons [21][22][23][24].Thanks to the availability of detailed urban topography and land use datasets [25,26], the application of LES to real urban environments (e.g., [27][28][29][30][31]) and also to directly support urban planning [32] has become possible over the last decade.In these cases the computational domain has several requirements to meet.Firstly, the computational domain has to be large enough to capture all relevant turbulence scales (e.g., [33]) and to minimize uncertainties related to the domain boundary conditions [11].Furthermore, it is assumed that the whole ABL should be included vertically.Secondly, the grid spacing has to be small enough to explicitly resolve the turbulence scales containing most of the energy [31], which requires large computational resources.
The few ventilation studies conducted over a real urban surface [30,32,34] have applied a simple, indirect ventilation indicator, velocity ratio v r = v p /v ∞ , i.e., the ratio between the wind velocities at pedestrian level v p (height z = 2 m above ground level (a.g.l.)) and at the top of the model domain v ∞ , where the flow is no longer influenced by the urban surface, rather than having examined the pollutant concentration, and processes of ventilation and dispersion together in more detail.This study employs an LES model, coupled with a Lagrangian stochastic particle model, to perform high-resolution urban flow simulations in order to answer: How does the structural layout of densely arranged building blocks along a city boulevard influence the ventilation of the local traffic emissions?The accumulation and ventilation of traffic-related air pollutants are investigated in four virtual city-block-design alternatives that are immersed within a real complex urban environment.The question is approached by quantifying and comparing the ventilation efficiency within perimeter blocks using three different measures of air pollutant ventilation and dispersion that are more sophisticated than v r : concentration, vertical turbulent transport and dilution rate of Lagrangian particles that represent non-reactive gaseous air pollutants.The ventilation measures are further examined in parallel.Simulations are conducted applying two contrasting meteorological conditions for the inflow.

Model Description
The LES model employed is the Parallelized Large-Eddy Simulation Model (PALM) version 4.0 (revision 1904) for atmospheric and oceanic flows [35,36], which solves the three-dimensional fields of wind and scalar variables (e.g., potential temperature and scalar concentrations).PALM has been applied to various types of ABL studies, for example to study cloudy boundary layers (e.g., [37]), the aerodynamic impact of a plant canopy (e.g., [38]) and stable boundary layer [39].The performance of PALM over an urban-like surface has been validated against wind tunnel simulations, previous LES studies and field measurements [28][29][30]40].The aerodynamic impact of vegetation is taken into account in PALM by means of an embedded canopy model.Vegetation decelerates the flow due to the form and viscous drag forces, and the decelerating force depends on the wind velocity and plant area density (PAD, m 2 m −3 ).Despite the model has previously been developed for continuous vegetation, it has also been used for individual trees (e.g., [19,41]).For this study, the canopy model was revised to allow a heterogeneous distribution of tree canopy.
Technical specifications of PALM are represented in Appendix A. As a new feature in PALM, a full three-dimensional two-way self-nesting is utilised for the first time [42].In the nesting approach, a "child" computational domain with chosen grid spacing and dimensions is defined inside the "parent" computational domain.PALM is run in parallel in both domains with respective computational set-ups (Table 1) and the domains communicate with each other.Nesting enables to have both a large computational domain and high enough resolution in the main area of interest without making the simulation computationally too expensive.[35].

Bottom and solid walls
No-slip condition for the horizontal wind components u and v (i.e., u = v = w = 0 m s −1 ).For potential temperature θ, the vertical gradient ∂θ/∂z = 0 K m −1 .Monin-Obukhov similarity theory (MOST) is applied between any solid-wall boundary and the first grid level normal to the respective boundary surface.
Same as for the parent.

Top
Dirichlet condition, i.e., u = U g (geostrophic wind speed) and v = w = 0 m s −1 .θ is extrapolated using the initial gradient of θ from a precursor run.
Two-way nesting.Boundary conditions obtained from the parent domain.
Horizontal: Lateral Cyclic boundary conditions Two-way nesting.Boundary conditions obtained from the parent domain.
Horizontal: streamwise Non-cyclic.A time-dependent turbulent inflow is produced by a turbulence recycling method [43].Requires a precursor run that is carried out over a domain of the same vertical extent as the parent domain and 1/16 in area.
Two-way nesting.Boundary conditions obtained from the parent domain.
Ventilation of traffic-related air pollutants is studied by applying a Lagrangian stochastic particle model (LPM) [35].LPM allows studying pollutant accumulation inside urban structures as concentration values can become higher than at the source unlike in the commonly used Eulerian method (e.g., [31]) due to the maximum principle for scalar conservation laws (e.g., [44]).In LPM, Lagrangian particles are released inside assigned source volumes at selected moments in time, after which they are transported inside the computational domain by the flow field solved by PALM.In this study, Lagrangian particles are defined passive and massless, and hence they represent non-reactive gaseous air pollutants.As for reactive gases and aerosol particles, chemical reactions and dynamic processes manipulating the concentrations should be included.Furthermore, deposition of air parcels on surfaces is not considered here.LPM in PALM was revised for this study to allow for the horizontal heterogeneity and relative strengths of air pollutant sources, i.e., streets with traffic.Particles are released at a selected height with a constant release rate Q (s −1 m −2 ), whereas the horizontal locations and group numbers i (1 ≤ i ≤ 10) of different sources are read from an input file.In order to consider the difference in emission strengths, each source group i can be given a weighting factor w i and the weighted particle concentration pc weight is then calculated as where pc i is the particle concentration inside a grid box.Weighting method is applied in order to utilise the same, sufficiently high, particle release rate throughout the source area.This approach allows LPM to represent a nearly continuous emission source everywhere without the need to utilise an excessive number of particles there where emissions are high.

Modelling Area
The simulations are conducted in western Helsinki on a city boulevard that is planned to be built in the new City Plan of Helsinki for 2050.The boulevard will be framed by a densely built neighbourhood, and in this study, four different city-block-design versions are being investigated.The design versions V type , where type is either par (parallel), per (perpendicular), perHV (perpendicular with height variation) or parJJ (parallel with Jin-Jang shape variation) are visualised in Figure 1 and their specific characteristics are listed in Table 2.In all design versions, the width of the boulevard is 54 m and the total length is around 3.3 km (Figure 2).The average floor area is set to a constant value and the average number of floors is eight.The different city-block-designs are described at the level of detail which does not take into account balconies and bay windows, for example.Furthermore, all planned buildings along the boulevard have flat roofs.

V type Characteristics
V par Building blocks by the boulevard are oriented so that the longest side is parallel to the boulevard.Building heights are fixed to 30 m.
V per Building block by the boulevard are oriented so that the longest side is perpendicular to the boulevard.Building heights are fixed to 30 m.

V perHV
The orientation of the building blocks is similar to V per but the building height varies.The highest buildings are situated at the nodal points of the public transport, whereas the lowest buildings as well as open urban spaces are situated between the nodal points.Buildings on the eastern side of the boulevard are generally higher.
V parJJ A so-called "Jin-Jang" block model, in which the buildings are similar to those in V par but the base height is lower and tower-like structures set above the base.Thus the building shape and height are very irregular.

Figure 1.
The city-block-design alternatives.From left to right: V par (the longer sides of buildings parallel to the boulevard), V per (the longer sides perpendicular to the boulevard), V perHV (similar to V per but the building height varies) and V parJJ (similar to V par but with Jing-Jang building structure variation).The boulevard is marked in red and the major junction with a black dot.for southwesterly or easterly inflow conditions, respectively, is defined.The grid spacing in all directions (x, y, z) is 1.0 m in the child and 2.0 m in the parent domain.

Urban Surface Data
The information on the surface elevation and cover for each x,y-pixel is fed to the model as two-dimensional ASCII-formatted raster files, with separate files for the topography (Figure 2) and tree canopy (Figure 3).The topography data contains topography elevation Z (m) above sea level (a.s.l.) and solid flow obstacles that have a volume of at least one grid box, i.e., buildings but not vehicles.
Correspondingly, the tree canopy data contains the tree canopy height a.s.l.No overhanging structures are accepted in the current model version.The landform information is obtained from the archive of the National Land Survey of Finland with a grid resolution of 2 m and the mean forest height from the archive of the Natural Resources Institute of Finland with a grid resolution of 20 m.All the rest, including the existing buildings, planned city-blocks and their corresponding street network alternatives, modifications in the landform and height of the street trees to be planted, is provided by the City Planning Department of Helsinki.Topography and building information, and mean forest and street tree heights, are superimposed on a single raster map file, respectively, which is then pivoted according to the prevailing geostrophic wind direction in each simulation.Over the whole computational domain, Z is given in a vertical resolution of 2 m.The non-cyclic inflow condition applied in the simulation (see Table 1) requires a topography-free zone of around 1.4 km in width, starting from the inflow boundary, for turbulence recycling.In addition, buffer regions of 20 m and 80 m in width, where Z is smoothed towards 20 m and 35 m (see boundaries in Figure 2), are added at the outflow and remaining lateral boundaries, respectively.This is done in order to avoid computational instabilities and to satisfy the periodic boundary conditions at the lateral boundaries [45].

Tree Canopy Model
The tree canopy includes two rows of planted street trees in the middle of the boulevard and the surrounding forested areas (Figure 3).At each surface x,y-grid box with vegetation, a vertical profile of the plant area density (PAD) is defined (Figure 4).Street trees are modelled as 20-m-tall Tilia × vulgaris trees, for which the vertical PAD profile is constructed based on an experimental 5-year mean summertime value of the leaf area density LAD = 5.3 m 2 m −3 for Tilia × vulgaris trees in Viikki, Helsinki [46].In addition to leaf area, PAD considers also the trunk and stems.The trees are assumed to have a circular cone shape with the maximum PAD at the height of the lowest branches at 6.5 m and a surface value of 0.3 m 2 m −3 .Furthermore, the winter-time leaf-off period (see Section 2.3) is considered by decreasing PAD by 80% following previous studies (e.g., [47,48]).Here, the main interest is on the impact of the street trees along the boulevard.Thus, the PAD profiles of the trees outside the boulevard are constructed based on their height relative to the street trees.In all simulations, a drag coefficient C D = 0.2 is assigned following previous LES studies considering ABL flows individual trees (e.g., [41]) and forest canopies (e.g., [38,49]).Example profiles for smaller and taller trees outside of the boulevard plotted in green (dashed line with stars) and red (dashed line with triangles), respectively.Additionally, the winter time PAD profile for the street trees is given in black (solid line).

Lagrangian Stochastic Particle Model (LPM)
Equal to the surface cover information, the horizontal location and group of particle sources are given to LPM as a x,y-pixel raster file.The planned street network in the vicinity of the boulevard is defined as the particle source area (Figure 5).Hence, only the local traffic emissions are taken into account similar to previous numerical ventilation studies (e.g., [50]).This assumption is adequate enough since ventilation depends on the short-term temporal variation of the concentration of a substance, which is mainly governed by local sources.Particles are released at a constant rate of 0.25 m −2 s −1 within each 1 m × 1 m × 1 m source grid box.No particles are released below the street trees.The street network is divided into three particle groups (i = 1-3) according to the estimated traffic rates in year 2025 provided by the City Planning Department of Helsinki (Figure 5).Group 1 covers minor side streets with less than 3000 vehicles per day, group 2 medium streets with 3000 to 10,000 vehicles per day and group 3 the boulevard itself and a major crossroads with over 10,000 vehicles per day.To take into account the proportional difference of traffic rates, pc weight (see Equation ( 1)) is calculated by defining w i = 2i − 1 for each particle group.In order to remove the particles stuck in the computational domain and unnecessarily increasing the computational load, their maximum age is set to 800 s, after which particles disappear from the domain.This is a loosely calculated expected time for a particle to travel across the child domain with a speed of 2 m s −1 .
LPM is applied only inside the child domain (Figure 2).Particles are restrained from entering the parent domain by setting an absorption condition at all boundaries, except at the bottom where particles are reflected.Detailed information about LPM can be found in [51].

Computational Set-Up
The different boundary conditions used for the parent and child domains are listed in Table 1.The roughness length z 0 applied in the wall treatment using MOST is set to 0.05 m at all solid surfaces.A sufficiently developed turbulent inflow, with which the main runs are initialised, is created with a precursor runs.The characteristics of the inflow, including u, θ and Reynolds stress u w , are illustrated in Figure 6.The ABL depth is set in the precursor run by imposing a large vertical gradient to the potential temperature θ, i.e., a strong temperature inversion, from the ABL height up to the total height of the computational domain.This inversion prevents the ABL from increasing with time.In order to examine the influence of the meteorological conditions of the inflow on pollutant dispersion and ventilation, two contrasting conditions are employed in the simulations.The first runs are conducted applying "general" meteorological conditions in Helsinki with a neutrally stratified ABL of 200 m in height and a geostrophic wind of U g = 10 m s −1 from the southwest (WD = 225 • ) [52].The surface heat flux is set to zero to maintain the neutral atmospheric stability.Densely built urban areas tend to be unstably stratified in Helsinki (e.g., [53,54]), and thus applying a neutral stability is a conservative choice as the removal of air pollutants by turbulence is found to improve under unstable conditions [55,56].The second runs are performed applying "wintry" meteorological conditions with a moderately stable ABL of around 160 m in height and a geostrophic wind of U g = 8 m s −1 from the east (WD = 90 • ).This simulates conditions that usually lead to worst air quality events in Helsinki in winter if the Siberian high is prevailing.The atmospheric stability and vertical wind profile applied are follow those in Basu and Porté-Agel [57], and they are attained by initializing the precursor run with these profiles and applying a surface heat exchange rate of −0.006K m s −1 to maintain the atmospheric stability.Simulating a stably stratified ABL is a challenging task as turbulent eddies are small, and a relatively high grid resolution is needed in order to resolve most of the energy containing turbulence [39].Furthermore, two supportive simulations were conducted for all design versions: under a neutral stratification but with an easterly wind, and under stable stratification but with a southwesterly wind.Results of the simulations are shown in Appendix B, D and E.

Simulations and Data Output
Complete simulations are done in several parts.First, one precursor run per each inflow condition is carried out over one hour (3600 s) and the final state of the run is applied to initialize the main runs.The main run is carried out in batches of 55 min (batch 1) and 6 min (batch 2).The second batch starts from the final state of the first batch.The release of particles starts after 5 min from the start of batch 1 and is stopped at 56 min, i.e., one minute after the start of batch 2.
Data are output in three sequences over two different domains of 0.5 km 2 and 0.04 km 2 (Figure 2) after the transients formed in the initialisation of the main run have subdued.Data output 1 with a time interval of 5 s is collected over the large data output domain over a time span of 40 min starting at 15 min.The high frequency data output 2 starts after 50 min with an interval equal to the integration time step of the simulation of ca.0.07 s and 0.15 s under the neutrally and stably stratified inflow, respectively.Due to the high logging frequency, the data output 2 covers only the small data output domain.The data output 3 is collected over the large data output domain at a time interval of 5 s during the last 5 min of the simulation when the release of particles has been stopped.The vertical resolution is 1.0 m for all data output whereas the horizontal resolution is 1.0 m for the data output 2 and 2.0 m for the outputs 1 and 3.

Ventilation and Dispersion Measures
Ventilation and dispersion of pollutants in street canyons and courtyards is assessed by three different measures: number concentration pc (m −3 ), which from hereafter equals pc weight for the sake of simplicity, the vertical turbulent flux density F p (m −2 s −1 ) and dilution rate D (m −3 s −1 ) of particles.Various other measures have also been proposed in the literature, for instance the exchange velocity [58], the local purging flow rate and visitation frequency [59], the mean tracer age and age distribution [60], the net escape velocity [61], and the particle exchange rate which assumes horizontal homogeneity of p c [50].These measures are however applicable only for simplified urban areas.The velocity ratio [34], on the other hand, does not directly measure pollutant ventilation, and therefore is not used here.
F p is calculated as the covariance between the vertical wind velocity and the particle number concentration as follows: where w (t, x, y, z) and pc (t, x, y, z) are the instantaneous fluctuating vertical velocity and particle number concentration at point (x, y, z) at time t, and the overbar denotes the time average over a selected averaging period.Before calculating the covariance, linear de-trending is applied on both time series.Positive F p indicates upward transport, i.e., ventilation, and negative downward, i.e., re-entrainment of particles from air above.Hence, the higher F p the higher ventilation.F p has previously been observed to determine pollutant removal from street canyons in CFD studies over idealized street canyons (e.g., [15,17,[62][63][64]) as well as in wind tunnel simulations (e.g., [65]).By contrast, advective flux densities are shown to govern the redistribution of pollutants below the roof level [17].Therefore, to study pollutant ventilation, only the turbulent transport is examined here.F p is calculated at z = 20 m, which is the minimum roof height of all city-block-design versions and well above the surface.Furthermore, F p is calculated both from the high frequency data output 2 (F p,HF ) and the low frequency data output 1 (F p,LF ).As regards the low temporal frequency of 0.2 Hz of the data output 1, a proportion of the total F p may be missed and thus the exact values of F p,LF should be treated with caution.On the other hand, representativeness of F p,HF calculated only over the small data output domain is limited.Therefore, F p,HF is compared to F p,LF over the same domain in order to show the applicability of F p,LF in this comparative study.Particle dilution rate D is calculated after the particle source has been switched of at 56 min in batch 2 as where ∆t is the output time step.The higher D, the higher ventilation.According to the scalar conservation equation (e.g., Equation ( 5) in [35]), D equals to the sum of the advective and turbulent transport terms if no sources are present.D is calculated from the data output 3 inside an air volume below 20 m a.g.l.First, to minimize the dependence of D on the initial particle concentration, the initial total particle concentration values pc tot (t = 0) inside the analysis domain in each block-design version are normalized relative to that of design version V par .The analysis of D is based on two measures: volume average D(t) V over a volume of 0.015 km 3 (around 2.4 million grid cells) and time-column average D(x, y) t,z within each 2 m × 2 m × 20 m vertical column.The dilution of particles is observed to occur very rapidly as the total particle concentrations drop to half within 60-120 s.Hence, only the first 50 s of the data output 3 are selected for the analysis, which is also the averaging period for D(x, y) t,z .Sensitivity of D(x, y) t,z to the averaging period was tested by altering the averaging period by ±20 s (see Appendix F).

Results
From hereafter, all heights are given in heights above ground level (a.g.l.) unless otherwise specified.The major junction refers to the junction marked in Figure 1.Along with spatial visualisation, the different measures (pc, F p and D) are analysed separately over the boulevard (marked in Figure 1), other street canyons excluding the boulevard, courtyards, and surroundings, which is the surface area classified neither as a street canyon nor a courtyard, covering around 50% of the large data output domain.

Particle Concentration pc
The 40-min temporal mean and 90th percentile values of pc are calculated for two layers, z = 3-5 m and z = 9-11 m, in both inflow conditions.The mean values for both heights separated to the areas defined above are shown in Figure 7, and the horizontal distributions at z = 4 m are displayed in Figure 8.The mean distributions at z = 10 m and 90th percentiles at both heights follow closely the mean values at z = 4 m, which is hence considered to represent the horizontal variation of pc.
In general, the spatial variability of pc within the study domain is notable with higher values along the boulevard than in other areas.Under the general inflow conditions (Figure 8a), the concentration patterns in the southern and northern parts of the boulevard are distinctly different, as the mean wind direction relative to the boulevard changes from nearly parallel to oblique, respectively.The perpendicular wind component results in accumulation to the upwind side of the boulevard due to particle transport by canyon vortex circulation (e.g., [15]).This accumulation is however weak in V perHV .Even though the largest concentrations are commonly seen at the boulevard, some strong accumulation is formed along the side street west of the major junction.Accumulation here is particularly notable in V perHV with building height decreasing downwind, which is attributable to a typically weaker canyon vortex inside this type of step-down canyons (e.g., [65]).In the northern part, the hotspots of pc on the leeward side of the boulevard are nearly continuous along the boulevard in V par and V parJJ , whereas in V per and V perHV ventilation from the cross streets breaks the accumulation patterns.Under the wintry inflow conditions (Figure 8b), the mean wind direction is close to perpendicular to the boulevard, especially north of the major junction.A large forested area lies to the east of the boulevard, and hence its cross streets provide fresh, unpolluted air and decrease the mean concentrations along those streets and the boulevard.This is particularly emphasised in V per and V perHV which have more cross streets.In the southern part of the boulevard, the smallest particle accumulation on the leeward side of the boulevard is seen in V perHV .Yet, pc north of the major junction are not distinctly higher in V perHV than in V per despite building height decreasing downwind in V perHV in this wind direction, which is generally expected to weaken the canyon vortex and ventilation (e.g., [65]).
Notwithstanding the lower PAD values of street trees in the simulations under the wintry inflow conditions, pc at z = 4 m are approximately two-fold compared to the general inflow conditions, which demonstrates how strongly meteorology influences pollutant dispersion.Under both inflow conditions, pc values are of the same magnitude in V par and V parJJ .The courtyards are notably clean and pc remains low throughout all simulations.
The vertical dispersion of pc and mean flow structure inside the boulevard street canyon for two vertical cross Sections 1 and 2 marked in Figure 8 are shown in Figure 9a,b under the general, and Figure 9c,d under the wintry inflow conditions.In Figure 9a (cross Section 1), the wind is close to parallel to the canyon, i.e., towards the picture, and due to flow channelling and advective particle transport the concentrations remain low.In Figure 9b (cross Section 2), instead, the mean wind has a strong perpendicular component creating optimal conditions for a canyon vortex.However, the street trees block the flow and break the vortex in two.This explains the more consistent dispersion of pc inside the canyon and accumulation on both sides of the boulevard.Under the wintry inflow conditions, on the other hand, a uniform canyon vortex is formed, and it is emphasised in V par and V parJJ with longer street canyons along the boulevard.Yet the exact location varies between the city-block versions.Vertical maxima of pc under both inflow conditions are seen right below the maximum of the PAD profile at z = 6.5 m.In spite of the great spatial variability in pc between the block-design versions, the concentration characteristics of different block-design versions are similar at both 4 and 10 m levels in almost all cases, as shown in Figure 7. Above the boulevard, both the mean values and 90th percentiles are lowest in V perHV under both general and wintry inflow conditions, whereas the highest values are observed in V parJJ and V per under the general and wintry inflow conditions, respectively.In other street canyons, no systematic pattern is seen while in the courtyards and surroundings, V par has on average the lowest mean concentrations.

Turbulent Particle Flux F p
Due to data constrains, F p,HF was computed only for the small data output domain, which is visualised Appendix C. Furthermore, the lower resolution F p,LF was also calculated in order to compare the turbulent particle flux with the other metrics used in the study.To estimate possible uncertainties originating from the lower temporal resolution, Table 3 summarizes the mean values of both F p,HF and F p,LF in the same small data output domain separately for the boulevard, other street canyons, courtyards and surroundings.Overall, F p,HF and F p,LF show rather similar orders of magnitude between the design-versions except above the boulevard under the general inflow conditions and in the surroundings under the wintry inflow conditions.F p,LF generally overestimates flux density values, apart from above the courtyards and above the other street canyons under the general inflow conditions.However, a Student's one-sample t-test for the difference of the values of F p,HF and F p,LF shows that the difference is insignificant at 95% confidence level (p = 0.069).F p,HF and F p,LF values agree better under the wintry than the general inflow conditions, owing most likely to the nearly three times longer time step of the wintry simulations.The difference between F p,HF and F p,LF becomes even less significant (p = 0.18), if F p,HF under the general inflow conditions is calculated using only every other time step of the data output 2 when the data logging frequency is around 7 Hz.As a conclusion, F p,LF can be applied when comparing ventilation in different city-block-design versions.
Table 3. Horizontal mean of the high-frequency vertical turbulent particle flux density F p,HF (m −2 s −1 ) separately for the boulevard, other street canyons, courtyards and surroundings at z = 20 m under both general and wintry inflow conditions.Horizontal mean of the low-frequency vertical turbulent particle flux density F p,LF over the same domain is given in brackets.The horizontal distribution and mean values of F p,LF above the separate areas for all block-design versions under both inflow conditions are visualised in Figures 10 and 11.For all model runs, F p,LF appears to be on average positive indicating upward transport of particles and ventilation.This can be expected since the particle source is constantly maintaining pc at street level.In general, F p,LF inside the courtyards is close to zero or even negative, as the only sources are re-entrainment from above or advection from the building openings.F p,LF is slightly increased when the building height increases and correspondingly decreased when the building height decreases downwind in V perHV .Furthermore, F p,LF values are smaller under the wintry than the general inflow conditions.Under the general inflow conditions, F p,LF shows most positive values along the boulevard in V perHV , whereas in the other design versions areas of negative flux density (re-entrainment) also appear.Under the wintry inflow conditions, F p,LF is most positive in V per and V perHV .In V par , a vast area of negative flux density is observed south of the major junction.The same is observed in V parJJ as well but less pronounced.F p,LF is decreased on the leeward side of buildings and especially behind the corners of buildings and higher towers in V parJJ , which is reflected to particle concentrations.

Surroundings
Horizontal mean values in Figure 10 correspond well to the visual analysis based on Figure 11.As noted, values are higher under the general inflow conditions.F p,LF is mainly largest in V perHV and smallest in V par under both inflow conditions.Both V per and V parJJ show higher values than V par regardless of weather conditions.Furthermore, F p over the courtyards in V parJJ depends strongly on the inflow.

Particle Dilution Rate D
The time series of D(t) V under the general (a) and wintry (b) inflow conditions are illustrated in Figure 12.Values are given relative to D(t) V in V par .Overall, D(t) V shows high variability with time notwithstanding the spatial averaging of D over around 2.4 million grid points.Figure 12a for the general inflow conditions reveals no significant differences between V per , V perHV and V parJJ whereas they all show mostly smaller values than V par .Under the wintry inflow conditions (Figure 12b), D(t) V in V per and particularly in V perHV remain mostly above that in V par whereas D(t) V in V parJJ varies between above and below.
The horizontal distribution and mean values of D(x, y) t,z above the boulevard, other street canyons, courtyards and surroundings under both inflow conditions are displayed in Figures 13 and 14.The analysis area is the same as for pc and F p,LF .In general, D(x, y) t,z increases downwind under both general and wintry inflow conditions as particles in the upwind columns are transported quickly by the horizontal advection.Furthermore, values are twofold under the wintry compared to the general inflow conditions.This also shown in the supportive simulations in Appendix E. Especially V perHV shows efficient dilution to the surroundings.Sensitivity tests show how the averaging period over which D(x, y) t,z is calculated has minor impact on the results (Appendix F).
Horizontal mean values in Figure 14 quantify the differences in D(x, y) t,z between the design versions.Along the boulevard, dilution is most efficient in V par and V parJJ under the general inflow conditions, and in V perHV under the wintry inflow conditions.The efficient dilution to the surroundings in V perHV is seen also in the mean values, especially in the wintry inflow conditions.Moreover, V perHV shows mostly higher values than V per .Dilution inside the courtyards is weakest in V per .

Discussion
Pollutant dispersion is strongly governed by inflow conditions: wind parallel to a street canyon sweeps the particles away due to particle transport by horizontal advection, as also shown by Moon et al. [66], whereas perpendicular or oblique winds result in accumulation to the leeward side.Mean concentrations are two-fold under the wintry compared to the general inflow conditions, which stems from less efficient vertical mixing and pollutant transport in stable ABL (see also Appendix B).On the other hand, Figure 10 and Appendix D show that vertical transport is notably higher under a southwesterly wind despite the atmospheric stability.This is likely due to more efficient mechanical production of turbulence under a southwesterly wind as the inflow travels over a longer stretch of urban area before entering the area of interest, unlikely when the wind is from the east with mostly homogeneous vegetated areas.
Urban morphology is shown to have a significant role in the dispersion of particles, which is illustrated in Figures 7 and 8. Greatest concentrations along the boulevard and weakest vertical transport are observed when the longer sides of building blocks are parallel to the boulevard, which supports previous studies performed over idealised street canyons [63,67,68].Non-uniform building height in V perHV is shown to decrease pc accumulation on the leeward side of the boulevard street canyon, especially under a southwesterly wind when the building height mainly increases downwind [65,69,70].Apart from buildings, street trees can alter the flow field and canyon vortex, and thereby decrease pollutant ventilation.Therefore, a detailed representation of trees in urban flow and air quality simulations is of high importance.
Of all the design options, V perHV with building height variation results on average in the greatest ventilation rates and lowest street level concentrations above the street canyons, with around 7-9% lower concentrations at z = 4 m than the city-block-design alternative with highest concentrations inside street canyons.The result supports previous studies conducted over simplified urban topographies (e.g., [60,65,[71][72][73][74]).The building shape variability in V parJJ , instead, is mainly shown to increase F p,LF and improve ventilation in vertical, whereas the mean values of pc and D display equally efficient or even weaker ventilation than in V par .The result indicates that the aerodynamically rougher city-blocks in V parJJ do not notably improve the air quality near the pedestrian level.The irregularities in the building shape in V parJJ probably destroy the canyon vortex [75], which reduces advective particle transport.Further in-depth studies are needed to explain the last-mentioned result.
Applying a high-resolution LES model over a real urban topography of a vast extent provides a large amount of detailed information about the flow and concentration fields.However, the amount of output data is high and the complexity of the topography further complicates the analysis.The ventilation measures chosen to be applied in this study cover the pollutant dispersion mechanisms broadly and complement each other.Based on the results, the following generalisation can be deduced.The pedestrian level concentrations are lowest when F p are highest.Furthermore, under both inflow conditions the highest mean values of F p,LF are observed in the V perHV and lowest in V par , suggesting that F p is strongly determined by the urban morphology.Generally, D is correlated to pc, and the differences between their horizontal distributions stem from the flow changing the location of particles but not necessarily removing them from the study domain.A comparison of the horizontal patterns of D(x, y) t,z and F p,LF points out areas where the dilution in vertical is notable and vice versa.On the whole, however, the patterns are unlike, indicating that D measures mainly horizontal advection.This advective nature of D can also be concluded from the higher values when the flow is parallel to the longer sides of buildings as the flow is horizontally less blocked and disturbed by flow structures created by building corners.Last of all, pc as well as F p and D are smaller in courtyards and surroundings than inside street canyons indicating little particle transport from the streets.The vertical exchange above courtyards is shown to be low when the angle between the wind and principal courtyard axis is around 90 • , supporting results by Moonen et al. [76].
Owing to the large extent of the study domain, the analysis was confined to the temporal mean values of the analysis measures.Moreover, the number of different inflow conditions applied was limited to two due to the high computational expenses of the simulations.They were, however, carefully chosen to represent typical and worst-case conditions at the modelling site, but also for the results to apply as guidelines in all urban areas.Last of all, the study focused on the aerodynamic impacts and omitted, for instance, air pollutant chemistry and dynamics and anthropogenic heat sources.
Explicit rules can be difficult to give to optimize of flow and pollutant dispersion over complex urban terrains [71].To yield detailed information for a specific case, similar studies have to be conducted separately for each urban planning solution.However, the outcome of this study, related firstly to the most efficient dispersion in V perHV with variable building height and short wall facing the main road and secondly to the insignificant improvement of ventilation by building shape variation in V parJJ , can be used as a general guideline also for other areas.Such information is vital for urban planners in order to design dense urban areas where the level of pollution exposure is minimized.For decision making, a separate model could be used in which each of the ventilation measures and areas of interest would get a weighting factor.Based on the factors that are considered important, for instance good ventilation inside courtyards, one could define the most suitable city-block-design alternative.

Conclusions
This study examines how the structural arrangement and orientation of perimeter blocks affect the ventilation and dispersion of traffic-related emissions within their street canyons under two contrasting meteorological conditions of the inflow.The principal objective of this numerical study is to demonstrate means to generate information that enables future urban planners to improve the pedestrian level air quality within real urban areas.The study is comprised of simulations that feature four different virtual building-block arrangements that are proposed for a real city boulevard site.In addition to buildings, the representation of the urban area is highly realistic including surrounding landform as well as forested areas and street trees.The simulations are conducted using a LES model PALM with an embedded Lagrangian stochastic particle model.The Lagrangian particles representing massless and inert air pollutant tracers are released over streets with a source strength (i.e., weighted particle release rate) that is relative to the local traffic volumes.The study employs a full three-dimensional two-way self-nesting functionality, which is newly implemented to PALM.
The analysis is founded on three different measures: particle concentration pc, vertical turbulent particle transport F p and particle dilution rate D. Variation in building height is shown to enhance both F p and D, and thus decrease accumulation of pc at pedestrian level.Furthermore, short canyons by the boulevard are shown preferential for F p .However, building shape variability of smaller scale does not result in improved ventilation.D is governed by horizontal advection and is thus strongly determined by the horizontal wind direction.It is shown how the design version with the shortest wall parallel to the boulevard with variable building height creates on average the most optimal air quality at pedestrian level along the boulevard, but the spatial variability within the street is highly variable depending on the pollutant transport and dispersion.Courtyards remain clean throughout all simulations, implying that in general traffic-related pollutants are not easily transported to there.Despite the high roughness and complexity of the urban surface that leads to efficient mechanical production of turbulent, pollutant concentrations are clearly higher under stably stratified inflow conditions.This is the first LES study over a vast, urban area applying sophisticated measures to assess pollutant dispersion and ventilation.The numerical methods are novel and highly developed providing realistic estimations for the removal of nonreactive gaseous air pollutants from the pedestrian level.The results of this study provide unique information about the transport of traffic-related pollutants in this specific urban environment and the results can directly be applied by local urban planners, but naturally the findings support urban planning also in other cities.

Figure 2 .
Figure 2. The computational domain for southwesterly (WD = 225 • ) inflow conditions and V par where the longer sides of building blocks are parallel to the boulevard.The domain is separated into child (B) and parent (A) domains, where the white colour stands for the topography elevation Z = 0 m.The data output domains are marked with rectangles: the small domain with a solid black and the large domain with a black dashed line.The boulevard is marked in red.

Figure 3 .
Figure 3.A map of the tree canopy height Z canopy (m) for the whole computational domain.Orientation as in Figure 2.

Figure 4 .
Figure 4.The vertical profile of PAD (m 2 m −3 ) for the street along the boulevard (blue dashed line with dots).Example profiles for smaller and taller trees outside of the boulevard plotted in green (dashed line with stars) and red (dashed line with triangles), respectively.Additionally, the winter time PAD profile for the street trees is given in black (solid line).

Figure 5 .
Figure 5.The particle source area (in colours) in the block-design alternatives V par and V parJJ .Source areas are divided into groups 1, 2 and 3 based on the estimated mean traffic rates in year 2025 (see legends).Street surfaces below trees (white dots) are omitted as source areas.Orientation as in Figure 2.

Figure 6 .
Figure 6.The vertical profile of (a) streamwise velocity u (m s −1 ), (b) Reynolds stress u w (m 2 s −2 ), and (c) potential temperature θ (K) at the inflow boundary for the general (solid line) and wintry (dashed line) inflow conditions.The top x-axis for θ is for the wintry conditions.

Figure 7 .
Figure 7. 40-min horizontal mean particle concentrations pc (m −3 ) under the (a) general and (b) wintry inflow conditions separately above the boulevard, other street canyons, courtyards and surroundings at height z = 4 m (bars with solid lines) and z = 10 m (bar with dashed lines) for all runs.90th percentile values are given with errorbars.The mean difference to the values in V par is given in percentages (∆V par ).

Figure 8 .
Figure 8. 40-min temporal mean of particle concentration pc (m −3 ) at height z = 4 m under the (a) general and (b) wintry inflow conditions.Notice the orientation of the mean wind and different scales of pc.Cross sections in Figure 8 are marked with white lines.

Figure 9 .
Figure 9. 40-min mean particle concentration pc (m −3 ) across the boulevard.(Cross Sections 1 and 2, which are marked in Figure 8, are shown in (a,b) under the general inflow conditions, and in (c,d) under the wintry inflow conditions, respectively.Cross sections are viewed from the south.Lengths of the wind arrows are relative to the wind speeds normal to the cross section.Height in m a.g.l. are given in the left.Notice the different scales of pc.)

Figure 10 .
Figure 10.Horizontal mean of the low-frequency vertical turbulent particle flux density F p,LF (m −2 s −1 ) under the (a) general and (b) wintry inflow conditions separately for the boulevard, other street canyons, courtyards and surroundings at z = 20 m for all runs.The difference to the value in V par is given in percentages (∆V par ).

Figure 11 .
Figure 11.40-min mean low-frequency vertical turbulent particle flux density F p,LF (m −2 s −1 ) at z = 20 m under the (a) general and (b) wintry inflow conditions.Positive values indicate upward flux.

Figure 12 .
Figure 12.The volume averaged particle dilution rate D(t) V (m −3 s −1 ) between z = 1-20 m under the (a) general and (b) wintry inflow conditions.Results are represented relative to V par ( D(t) V,V par ).

Figure 13 .
Figure 13.Mean D(x, y) t,z (×10 −3 m −3 s −1 ) under the (a) general and (b) wintry inflow conditions separately for the boulevard, other street canyons, courtyards and surroundings between z = 1-20 m for all runs.D is calculated using data from the first 50 s after the particle source has been switched off.The difference to the value in V par is given in percentages (∆V par ).

Figure 14 .
Figure 14.The temporal mean particle dilution rate D(x, y) t,z (m −3 s −1 ) between z = 1-20 m for the first 60 s after the particle has been switched off: (a) general and (b) wintry inflow conditions.

Table 1 .
The boundary conditions of the model runs.Details of the conditions can be found in Maronga et al.

Table 2 .
The specific characteristics of different city-block-design alternatives.