Strategic guidelines for street canyon geometry to achieve sustainable street air quality

This paper is concerned with the motion of air within the urban street canyon and is directed towards a deeper understanding of pollutant dispersion with respect to various simple canyon geometries and source positions. Taking into account the present days typical urban configurations, three principal flow regimes “isolated roughness flow”, “skimming flow” and “wake interference flow” (Boundary Layer Climates, 2nd edition, Methuen, London) and their corresponding pollutant dispersion characteristics are studied for various canopies aspect ratios, namely relative height (h₂/h₁), canyon height to width ratio (h/w) and canyon length to height ratio (l/h). A field-size canyon has been analyzed through numerical simulations using the standard k-ε turbulence closure model. It is found that the pollutant transport and diffusion is strongly dependent upon the type of flow regime inside the canyon and exchange between canyon and the above roof air. Some rules of thumbs have been established to get urban canyon geometries for efficient dispersion of pollutants.     

Numerical modeling on air quality in an urban environment with changes of the aspect ratio and wind direction

Due to heavy traffic emissions within an urban environment, air quality during the last decade becomes worse year by year and hazard to public health. In the present work, numerical modeling of flow and dispersion of gaseous emissions from vehicle exhaust in a street canyon were investigated under changes of the aspect ratio and wind direction. The three-dimensional flow and dispersion of gaseous pollutants were modeled using a computational fluid dynamics (CFD) model which was numerically solved using Reynolds-averaged Navier–Stokes (RANS) equations. The diffusion flow field in the atmospheric boundary layer within the street canyon was studied for different aspect ratios (W/H?=?1/2, 3/4, and 1) and wind directions (θ?=?90°, 112.5°, 135°, and 157.5°). The numerical models were validated against wind tunnel results to optimize the turbulence model. The numerical results agreed well with the wind tunnel results. The simulation demonstrated that the minimum concentration at the human respiration height within the street canyon was on the windward side for aspect ratios W/H?=?1/2 and 1 and wind directions θ?=?112.5°, 135°, and 157.5°. The pollutant concentration level decreases as the wind direction and aspect ratio increase. The wind velocity and turbulence intensity increase as the aspect ratio and wind direction increase.     

Modeling reactive pollutant dispersion in an urban street canyon

Reactive pollutant dispersion in an urban street canyon with a street aspect ratio of one is numerically investigated using a computational fluid dynamics (CFD) model. The CFD model developed is a Reynolds-averaged Navier–Stokes equations (RANS) model with the renormalization group (RNG) k–ε turbulence model and includes transport equations for NO, NO₂, and O₃ with simple photochemistry. An area emission source of NO and NO₂ is considered in the presence of background O₃ and street bottom heating (ΔT=5 °C) with an ambient wind perpendicular to the along-canyon direction. A primary vortex is formed in the street canyon and the line connecting the centers of cross-sectional vortices meanders over time and in the canyon space. The cross-canyon-averaged temperature and reactive pollutant concentrations oscillate with a period of about 15 min. The averaged temperature is found to be in phase with NO and NO₂ concentrations but out of phase with O₃ concentration. The photostationary state defect is small in the street canyon except for near the roof level and the upper downwind region of the canyon and its local minimum is observed near the center of the primary vortex. The budget analysis of NO (NO₂) concentration shows that the magnitude of the advection or turbulent diffusion term is much larger (larger) than that of the chemical reaction term and that the advection term is largely balanced by the turbulent diffusion term. On the other hand, the budget analysis of O₃ concentration shows that the magnitude of the chemical reaction term is comparable to that of the advection or turbulent diffusion term. The inhomogeneous temperature distribution itself affects O₃ concentration to some extent due to the temperature-dependent photolysis rate and reaction rate constant.     

On the correlation of air and pollutant exchange for street canyons in combined wind-buoyancy-driven flow

The ventilation and pollutant transport in a two-dimensional (2D) street canyon of building-height-to-street-width (aspect) ratio h/b = 1 under different unstable stratifications were examined. To characterize the combined wind-buoyancy-driven flow and pollutant transport at different Richardson number Ri, a computational fluid dynamics (CFD) model based on the Reynolds-averaged Navier–Stokes (RANS) equations with the Renormalization Group (RNG) k ? ε turbulence model was adopted. Unlike the isothermal condition, a secondary recirculation is initiated at the ground-level windward corner of the street canyon once the unstable stratification is switched on (Ri < 0). It traps the ground-level pollutant leading to elevated pollutant concentration there. As Ri further decreases, the enlarging secondary recirculation enables direct pollutant removal from its core to the shear layer that offsets the ground-level pollutant accumulation. The ventilation and pollutant removal performance under different unstable stratifications are compared by the air (ACH) and pollutant (PCH) exchange rates, and pollutant retention time (τ). Both the mean and turbulent components of ACH are found to increase with decreasing Ri, suggesting that unstable stratification promotes ventilation in street canyons. Moreover, the CFD results agree well with our theoretical model that ACH² varies linearly with Ri. Turbulent transport originally dominates the pollutant removal under isothermal condition. However, progressive domination of pollutant removal by mean wind can be observed with decreasing stability (decreasing Ri from 0 to ?10.6). The critical value is estimated to be Ri = ?8, below which mean wind is the major pollutant removal carrier. Reduction in τ is also observed with decreasing Ri. Hence, in unstable stratification, pollutant resides shorter time in the street canyon compared with its isothermal counterpart, and the ventilation and pollutant removal are more favorable.     

Description of pollutant dispersion in an urban street canyon using a two-dimensional lattice model

The pollutant dispersion in a street canyon has been described in this work by using an isothermal two-dimensional lattice model coupled to the Smagorinsky sub-grid scale model. The influence of the ratio between the height of the upstream and downstream canyon walls, as well as the gap distance between them on the flow pattern, was analyzed considering the situations of ‘open country’ or isolated street canyon and ‘urban roughness’ in which the influence of an urban fabric was considered. The model determined the trajectories of a large number of passive tracer particles released in the computational domain, making it easy to visualize the flow regimes established in each case. The results agreed with the observations reported from the experiments showing a strong influence on the flow inside the canyon exerted by the upstream landscape configuration.     

Three-dimensional modeling of air flow and pollutant dispersion in an urban street canyon with thermal effects

Effects of excess ground and building temperatures on airflow and dispersion of pollutants in an urban street canyon with an aspect ratio of 0.8 and a length-to-width ratio of 3 were investigated numerically. Three-dimensional governing equations of mass, momentum, energy, and species were modeled using the RNG k-epsilon turbulence model and Boussinesq approximation, which were solved using the finite volume method. Vehicle emissions were estimated from the measured traffic flow rates and modeled as banded line sources, with a street length and bandwidths equal to typical vehicle widths. Both measurements and simulations reveal that pollutant concentrations typically follow the traffic flow rate; they decline as the height increases and are higher on the leeward side than on the windward side. Three-dimensional simulations reveal that the vortex line, joining the centers of cross-sectional vortexes of the street canyon, meanders between street buildings and shifts toward the windward side when heating strength is increased. Thermal boundary layers are very thin. Entrainment of outside air increases, and pollutant concentration decreases with increasing heating condition. Also, traffic-produced turbulence enhances the turbulent kinetic energy and the mixing of temperature and admixtures in the canyon. Factors affecting the inaccuracy of the simulations are addressed.     

Impact of building facades and ground heating on wind flow and pollutant transport in street canyons

This paper investigates the impacts of building facades and ground heating on the wind flow and pollutant transport in street canyons using the computational fluid dynamic (CFD) technique. Street canyons of H/W (H representing the building height and W the street width) varied from 0.1 to 2, which covered the basic flow regimes of skimming flow (H/W=1 or 2), wake interference flow (H/W=0.5), and isolated roughness flow (H/W=0.1), were examined in a series of sensitivity tests. Heating that occurred on different surfaces, including ground surface and building façades, posed considerable effects on the street canyon wind flow and pollutant transport compared with those under isothermal conditions. The CFD results showed that the mechanically induced wind flow and pollutant transport were complicated by the buoyancy under temperature stratification. Individual street canyons of different H/W and surface-heating scenarios exhibited their unique wind flow structure and pollutant transport behaviors. Two counter-rotating vortices were calculated in the street canyons of H/W=1, in which the zone of higher pollutant concentration under isothermal conditions was switched from the leeward side to the windward side. In the street canyon of H/W=2, the recirculating wind pattern was perturbed by surface heating that led to the development of either one primary vortex or three closely coupled vortices. Because of the complicated wind structure, the zones of higher pollutant concentration located either on the leeward or windward ground level were subjected to the surface-heating scenarios. Only two vortices were developed inside the street canyon of H/W=0.5. The large primary vortex, centered inside the street canyon, extended above the roof level of the street canyon. Meanwhile, a small secondary vortex was found at the ground-level windward corner whose size results as a function of surface-heating configurations. Finally, in the street canyon of H/W=0.1, an isolated clockwise-rotating vortex was developed beside the leeward building while the wind in the windward side blew in the prevailing wind direction. As a result, air pollutant emitted at the street centerline was unlikely to be carried into the leeward vortex. Instead, it was dispersed rapidly on the windward side before being removed from the street canyon.     

Analysis of air quality within a street canyon using statistical and dispersion modelling techniques

The dispersion model, ADMS-Urban, alongside the statistical modelling technique, generalized additive modelling, have been used to predict hourly NO_x and nitrogen dioxide (NO₂) concentrations at a busy street canyon location and the results compared with measurements. Generalized additive models (GAMs) were constructed for NO₂ and NO_x concentrations using input data required to run ADMS-Urban. Bivariate polar plots have been produced from the wind flow (speed and direction) and pollution data (measured and predicted concentrations) to provide further information regarding the complex wind-pollutant interactions in an urban street canyon. The predictions made with the GAMs show excellent agreement with measured concentrations at this location, reproducing both the magnitude of NO_x and NO₂ concentrations and also the wind speed-wind direction dependence of pollutant sources within the canyon. However, the predictions made with ADMS-Urban under-estimated the measured NO_x by 11% and NO₂ by 21% and there are clear differences in the bivariate polar plots. Several sensitivity tests were carried out with ADMS-Urban in an attempt to produce predictions in closer agreement to those measured at Gillygate. Increasing the primary NO₂ fraction in ADMS-Urban (from 10% to 20%) had a considerable effect on the predictions made with this model, increasing NO₂ predictions by ∼20%. However, the bivariate plots still showed major differences to those of the measurements. This work illustrates that generalized additive modelling is a useful tool for investigating complex wind-pollutant interactions within a street canyon.     

Modelling dispersion of traffic pollution in a deep street canyon: Application of CFD and operational models

In this study, numerical modelling of the flow and concentration fields has been undertaken for a deep street canyon in Naples (Italy), having aspect ratio (i.e. ratio of the building height H to the street width W) H/W = 5.7. Two different modelling techniques have been employed: computational fluid dynamics (CFD) and operational dispersion modelling. The CFD simulations have been carried out by using the RNG k–? turbulence model included in the commercial suite FLUENT, while operational modelling has been conducted by means of the WinOSPM model. Concentration fields obtained from model simulations have been compared with experimental data of CO concentrations measured at two vertical locations within the canyon. The CFD results are in good agreement with the experimental data, while poor agreement is observed for the WinOSPM results. This is because WinOSPM was originally developed and tested for street canyons with aspect ratio H/W ≌ 1. Large discrepancies in wind profiles simulated within the canyon are observed between CFD and OSPM models. Therefore, a modification of the wind profile within the canyon is introduced in WinOSPM for extending its applicability to deeper canyons, leading to an improved agreement between modelled and experimental data. Further development of the operational dispersion model is required in order to reproduce the distinct air circulation patterns within deep street canyons.     

Air ventilation impacts of the “wall effect” resulting from the alignment of high-rise buildings

The objective of this study is to investigate the air ventilation impacts of the so called “wall effect” caused by the alignment of high-rise buildings in complex building clusters. The research method employs the numerical algorithm of computational fluid dynamics (CFD – FLUENT) to simulate the steady-state wind field in a typical Hong Kong urban setting and investigate pollutant dispersion inside the street canyon utilizing a pollutant transport model. The model settings of validation study were accomplished by comparing the simulation wind field around a single building block to wind tunnel data. The results revealed that our model simulation is fairly close to the wind tunnel measurements. In this paper, a typical dense building distribution in Hong Kong with 2 incident wind directions (0° and 22.5°) is studied. Two performance indicators are used to quantify the air ventilation impacts, namely the velocity ratio (VR) and the retention time (T_r) of pollutants at the street level. The results indicated that the velocity ratio at 2 m above ground was reduced 40% and retention time of pollutants increased 80% inside the street canyon when high-rise buildings with 4 times height of the street canyon were aligned as a “wall” upstream. While this reduction of air ventilation was anticipated, the magnitude is significant and this result clearly has important implications for building and urban planning.     

A multi-approach monitoring of particulate matter, metals and PAHs in an urban street canyon

For the first time until now, the results from a prediction model (Atmospheric Dispersion Modelling System (ADMS)-Road) of pollutant dispersion in a street canyon were compared to the results obtained from biomonitors. In particular, the instrumental monitoring of particulate matter (PM10) and the biomonitoring of 14 polycyclic aromatic hydrocarbons (PAHs) and 11 metals by Quercus ilex leaves and Hypnum cupressiforme moss bags, acting as long- and short-term accumulators, respectively, were carried out. For both PAHs and metals, similar bioaccumulation trends were observed, with higher concentrations in biomonitors exposed at the leeward canyon side, affected by primary air vortex. The major pollutant accumulation at the leeward side was also predicted by the ADMS-Road model, on the basis of the prevailing wind direction that determines different exposure of the street canyon sides to pollutants emitted by vehicular traffic. A clear vertical (3, 6 and 9 m) distribution gradient of pollutants was not observed, so that both the model and biomonitoring results suggested that local air turbulences in the street canyon could contribute to uniform pollutant distribution at different heights.     

Simulation of wind-driven dispersion of fire pollutants in a street canyon using FDS

Air quality in urban areas attracts great attention due to increasing pollutant emissions and their negative effects on human health and environment. Numerous studies, such as those by Mouilleau and Champassith (J Loss Prevent Proc 22(3): 316–323, 2009), Xie et al. (J Hydrodyn 21(1): 108–117, 2009), and Yassin (Environ Sci Pollut Res 20(6): 3975–3988, 2013) focus on the air pollutant dispersion with no buoyancy effect or weak buoyancy effect. A few studies, such as those by Hu et al. (J Hazard Mater 166(1): 394–406, 2009; J Hazard Mater 192(3): 940–948, 2011; J Civ Eng Manag (2013)) focus on the fire-induced dispersion of pollutants with heat buoyancy release rate in the range from 0.5 to 20 MW. However, the air pollution source might very often be concentrated and intensive, as a consequence of the hazardous materials fire. Namely, transportation of fuel through urban areas occurs regularly, because it is often impossible to find alternative supply routes. It is accompanied with the risk of fire accident occurrences. Accident prevention strategies require analysis of the worst scenarios in which fire products jeopardize the exposed population and environment. The aim of this article is to analyze the impact of wind flow on air pollution and human vulnerability to fire products in a street canyon. For simulation of the gasoline tanker truck fire as a result of a multivehicle accident, computational fluid dynamics large eddy simulation method has been used. Numerical results show that the fire products flow vertically upward, without touching the walls of the buildings in the absence of wind. However, when the wind velocity reaches the critical value, the products touch the walls of the buildings on both sides of the street canyon. The concentrations of carbon monoxide and soot decrease, whereas carbon dioxide concentration increases with the rise of height above the street canyon ground level. The longitudinal concentration of the pollutants inside the street increases with the rise of the wind velocity at the roof level of the street canyon.     

Assessing the representativeness of monitoring data from an urban intersection site in central London,UK

The wind flow field around urban street-building configurations has an important influence on the microscale pollutant dispersion from road traffic, affecting overall dilution and creating localised spatial variations of pollutant concentration. As a result, the “representativeness” of air quality measurements made at different urban monitoring sites can be strongly dependent on the interaction of the local wind flow field with the street-building geometry surrounding the monitor. The present study is an initial attempt to develop a method for appraising the significance of air quality measurements from urban monitoring sites, using a general application computational fluid dynamics (CFD) code to simulate small-scale flow and dispersion patterns around real urban building configurations. The main focus of the work was to evaluate routine CO monitoring data collected by Westminster City Council at an intersection of street canyons at Marylebone Road, Central London. Many monitors in the UK are purposely situated at urban canyon intersections, which are thought to be local “hot spots” of pollutant emissions, however very limited information exists in the literature on the flow and dispersion patterns associated with them. With the use of simple CFD simulations and the analysis of available monitoring data, it was possible to gain insights into the effect of wind direction on the small-scale dispersion patterns at the chosen intersection, and how that can influence the data captured by a monitor. It was found that a change in wind direction could result in an increase or decrease of monitored CO concentration of up to 80%, for a given level of traffic emissions and meteorological conditions. Understanding and de-coupling the local effect of wind direction from monitoring data using the methods presented in this work could prove a useful new tool for urban monitoring data interpretation.     

City breathability and its link to pollutant concentration distribution within urban-like geometries

This paper is devoted to the study of pollutant concentration distribution within urban-like geometries. By applying efficiency concepts originally developed for indoor environments, the term ventilation is used as a measure of city “breathability”. It can be applied to analyse pollutant removal within a city in operational contexts. This implies the evaluation of the bulk flow balance over the city and of the mean age of air. The influence of building packing density on flow and pollutant removal is, therefore, evaluated using those quantities. Idealized cities of regular cubical buildings were created with packing density ranging from 6.25% to 69% to represent configurations from urban sprawl to compact cities. The relative simplicity of these arrangements allowed us to apply the Computational Fluid Dynamics (CFD) flow and dispersion simulations using the standard k–? turbulence model. Results show that city breathability within the urban canopy layer is strongly dependent from the building packing density. At the lower packing densities, the city responds to the wind as an agglomeration of obstacles, at larger densities (from about 44%) the city itself responds as a single obstacle. With the exception of the lowest packing density, airflow enters the array through lateral sides and leaves throughout the street top and flow out downstream. The air entering through lateral sides increases with increasing packing density.At the street top of the windward side of compact building configurations, a large upward flow is observed. This vertical transport reduces over short distance to turn into a downward flow further downstream of the building array. These findings suggest a practical way of identifying city breathability. Even though the application of these results to real scenarios require further analyses the paper illustrates a practical framework to be adopted in the assessment of the optimum neighbourhood building layout to minimize pollution levels.     

Measurements and analyses of nitrogen oxides and ozone in the yard and on the roof of a street-canyon in Suzhou

The concentrations of air pollutants such as nitrogen oxides and ozone characterised by very fast chemical reactions can significantly vary within urban street-canyon due to the short distances between sources and receptor. With the primary objective to analyse this issue, NO, NO₂, NO_x, O₃, BTX, and wind flow field were continuously measured for 1 week at two heights (a street-level yard and a 25-m-high rooftop) in an urban canyon in Suzhou (China). The yard ozone concentrations were found to be up to six times lower than on the roof. Different frequency distributions (FD), dynamical and chemical processes of the pollutant variations from yard to roof are discussed to explain the findings. The predominant factors for the dissimilar pollutant vertical diffusion at the two measurement locations were associated to dissimilar fluid-dynamic and heterogeneous removal effects that likely induced dissimilar ozone chemical processes relative to NO_x and BTX precursors.     

A note on average vertical profiles of vehicular pollutant concentrations in urban street canyons

Recent observations of air pollutant concentrations measured within and above street canyons were used to study the average vertical profiles of vehicular pollutant concentrations in the urban environment. The idea of an exponential vertical concentration distribution, exp( −Bz^q), resulted from a near ground-level source diffusing over flat terrain, was tentatively extended to the urban street canyons, where the empirical parameters B and q are generally dependent on the atmospheric stability and the aerodynamic characteristics of the canyon.     

A concentration correction scheme for Lagrangian particle model and its application in street canyon air dispersion modelling

Pollutant dispersion in street canyons with various configurations was simulated by discharging a large number of particles into the computation domain after developing a time-dependent wind field. Trajectory of the released particles was predicted using a Lagrangian particle model developed in an earlier study. A concentration correction scheme, based on the concept of “visibility”, was adopted for the Lagrangian particle model to correct the calculated pollutant concentration field in street canyons. The corrected concentrations compared favourably with those from wind tunnel experiments and a linear relationship between the computed concentrations and wind tunnel data were found. The developed model was then applied to four simulations to test for the suitability of the correction scheme and to study pollutant distribution in street canyons with different configurations. For those cases with obstacles presence in the computation domain, the correction scheme gives more reasonable results compared with the one without using it. Different flow regimes are observed in the street canyons, which depend on building configurations. A counter-clockwise rotating vortex may appear in a two-building case with wind flow from left to right, causing lower pollutant concentration at the leeward side of upstream building and higher concentration at the windward side of downstream building. On the other hand, a stable clockwise rotating vortex is formed in the street canyon with multiple identical buildings, resulting in poor natural ventilation in the street canyon. Moreover, particles emitted in the downstream canyon formed by buildings with large height-to-width ratios will be transported to upstream canyons.     

Pollutant dispersion in urban street canopies

The effects of building configurations on pollutant dispersion around street canopies were studied numerically. The dispersion of pollutants emitted from ground sources was simulated by continuously discharging large number of particles into the computation domain. The mean wind velocities at each time-step were firstly computed by solving the time-dependent incompressible Navier–Stokes equations, while the fluctuated velocities were determined using a statistical procedure. The trajectories of the discharged particles were obtained from a Lagrangian particle model. Three categories of numerical simulation were conducted to study the effect of different canopy geometries on the pollutant dispersion. The computed wind field data were consistent with the wind field characteristics described in the previous wind tunnel studies. A counter-clockwise vortex was found resulting in high pollutant concentration at the windward side of the downstream building of the street canopy and low pollutant concentration at the leeward side of the upstream building. The increase in height of the urban roughness buildings would facilitate the pollutant dispersion in urban street canopy under certain building configurations. Two or more vortices stacked vertically in a street canopy were found when height of the upstream and downstream buildings of a street canopy was increased, preventing pollutants from escaping out of the canopy.     

Flow and dispersion in an urban cubical cavity

Flow and dispersion in an urban cubical cavity are numerically investigated using a Reynolds-averaged Navier–Stokes equations (RANS) model with the renormalization group (RNG) k–? turbulence closure model. The urban cubical cavity is surrounded by flank walls that are parallel to the streamwise direction, called end-walls, as well as upstream and downstream walls. A primary vortex and secondary vortices including end-wall vortices are formed in the cavity. Because of the end-wall drag effect, the averaged mean-flow kinetic energy in the cavity is smaller than that in an urban street canyon that is open in the along-canyon direction. A trajectory analysis shows that the end-wall vortices cause fluid particles to move in the spanwise direction, indicating that flow in the cavity is essentially three-dimensional. The iso-surfaces of the Okubo–Weiss criterion capture cavity vortices well. The pollutant concentration is high near the bottom of the upstream side in the case of continuous pollutant emission, whereas it is high near the center of the primary vortex in the case of instantaneous pollutant emission. To get some insight into the degree of pollutant escape from the cavity according to various meteorological factors, extensive numerical experiments with different ambient wind speeds and directions, inflow turbulence intensities, and cavity-bottom heating intensities are performed. For each experiment, we calculate the time constant, which is defined as the time taken for the pollutant concentration to decrease to e^?1 of its initial value. The time constant decreases substantially with increasing ambient wind speed, and tends to decrease with increasing inflow turbulence intensity and cavity-bottom heating intensity. The time constant increases as the ambient wind direction becomes oblique. High ambient wind speed is found to be the most crucial factor for ventilating the cavity, thus improving air quality in an urban cubical cavity.     

Impacts of Upstream Building Width and Upwind Building Arrangements on Airflow and Pollutant Dispersion in a Street Canyon

A two-dimensional numerical model for evaluating the wind flow and pollutant dispersion within a street canyon was first developed using the FLUENT code, which was then validated against a wind tunnel experiment. Then, the effects of the upstream building width and upwind building arrangement on the airflow and pollutant dispersion inside an isolated street canyon were investigated numerically. The numerical results revealed that: (1) the in-canyon vortex center shifts downwards as the upstream building width increases; (2) the recirculation zone covers the entire upstream building roof for the cases when W/H = 0.5, 1.0, 1.5, and 2.0 (W is the upstream building width and H is the building height), whereas the flow reattaches the upstream building roof for the cases when W/H = 2.5 and 3.0; (3) when the upstream building width is shorter than the critical width W_C (= 2H), an increase in the upstream building width leads to an increase in the pollution level on the leeward wall of the canyon and a decrease in the roof-level concentrations at the upstream building; (4) when the upstream building width is longer than the critical width, the roof-level concentrations at the upstream building are negligibly small and the pollution level on the leeward wall of the canyon is almost unaffected by a further increase in the upstream building width; (5) when the buildings are placed upwind of the canyon, the flow attaches the upstream building roof and, therefore, almost none of the pollutants are distributed on the upstream building roof; and (6) the pollution levels inside the canyon and on the downstream building roof increase significantly with the number of upwind buildings.