Large-eddy simulation of pollution dispersion in an urban street canyon—Part I: comparison with field data

Large-eddy simulations (LESs) are applied to the problem of pollution dispersion within the urban canopy layer, specifically street canyons. The objective is to study the turbulence structure and hence the physical dispersion mechanisms of pollutants. LESs are implemented by incorporating the dynamic sub-grid scale stress model into the commercial computational fluids dynamics code CFX. To gain confidence in the approach, simulations are performed for a canyon-like geometry (roof garden) for which experimental measurements were also made. The experimental campaign consisted of using sonic anemometers to measure mean flow and turbulence intensities at a high sample rate of 60 Hz. Good agreement between simulations and experimental data are obtained. Real geometric features, such as non-uniform wall heights, result in a very much three-dimensional flow distribution. Comparisons with the k–ε model show that LESs are able to predict more accurately the turbulence statistics of the flow.     

A field study of factors influencing the concentrations of a traffic-related pollutant in the vicinity of a complex urban junction

The paper describes a field study focused on the dispersion of a traffic-related pollutant within an area close to a busy intersection between two street canyons in Central London. Simultaneous measurements of airflow, traffic flow and carbon monoxide concentrations ([CO]) are used to explore the causes of spatial variability in [CO] over a full range of background wind directions. Depending on the roof-top wind direction, evidence of both flow channelling and recirculation regimes were identified from data collected within the main canyon and the intersection. However, at the intersection, the merging of channelled flows from the canyons increased the flow complexity and turbulence intensity. These features, coupled with the close proximity of nearby queuing traffic in several directions, led to the highest overall time-average measured [CO] occurring at the intersection. Within the main street canyon, the data supported the presence of a helical flow regime for oblique roof-top flows, leading to increased [CO] on the canyon leeward side. Predominant wind directions led to some locations having significantly higher diurnal average [CO] due to being mostly on the canyon leeward side during the study period. For all locations, small changes in the background wind direction could cause large changes in the in-street mean wind angle and local turbulence intensity, implying that dispersion mechanisms would be highly sensitive to small changes in above roof flows. During peak traffic flow periods, concentrations within parallel side streets were approximately four times lower than within the main canyon and intersection which has implications for controlling personal exposure. Overall, the results illustrate that pollutant concentrations can be highly spatially variable over even short distances within complex urban geometries, and that synoptic wind patterns, traffic queue location and building topologies all play a role in determining where pollutant hot spots occur.     

A numerical experiment of roadside diffusion under traffic-produced flow and turbulence

Roadside air pollution due to heavy traffic is one of the unsettled issues in the atmospheric environment in urban areas. As a practical application of a Computational Fluid Dynamics (CFD) model, a coupled mesoscale-CFD model was applied to the Ikegamicho area of Kawasaki City, Japan. For this study, the effects of traffic-produced flow and turbulence (TPFT) on the dispersion of the pollutants near the heavy traffic road were mainly investigated in an actual urban area. First, a series of preliminary CFD calculations was conducted for a road tunnel field experiment to obtain a fitting parameter for the traffic-produced flow. The calculation was then performed for 24 h in December 2005 around Ikegamicho, and the results were compared with the data at a roadside monitoring post in the area, located 10 m from the boundary of the ground road. In general, the effect of traffic-produced flow and turbulence was limited at the downstream side of the roads. The maximum concentration of NO_x was reduced and smoothed out along the traffic flow by the traffic-produced flow and turbulence on the road. The effects of traffic-produced turbulence on the dispersion of pollutants were greater than those of traffic-produced flow; however, the effects of traffic-produced flow were not negligible. The concentration of pollutants was not particularly dependent on the turbulent Schmidt number because most of the emission sources were introduced as volume sources in the present calculations, and the effect caused by differences in the material diffusion coefficient was not particularly significant at the outside of the road.     

Application of the k–ε turbulence model to the high Reynolds number skimming flow field of an urban street canyon

A two-dimensional, steady, k–ε turbulence model was used to investigate the high Reynolds number skimming flow field of an urban street canyon. We describe the critical canyon width-to-height ratios that distinguish a cascade of vortex patterns that form in an urban street canyon. Details of the flow field are reported that includes the structure of the mean flow field, turbulent kinetic energy, turbulent length scale, turbulent eddy viscosity, and Reynolds stress for three typical different aspect ratios, W/H, of a street canyon. The consequences of vortex layering on vertical transport are explored.     

Spatio-temporal covariation of urban particle number concentration and ambient noise

Mobile measurements of ambient noise and particle number concentrations were carried out within an urban residential area in Essen, Germany, during summer 2008. A busy major road with a traffic intensity of about 44,000 vehicles per day was situated within the study area. The spatio-temporal distribution of noise and particles was closely coupled to road traffic on the major road. Total particle number concentrations in proximity to the main road were on average between 25,000 cm^?3 and 35,000 cm^?3 while sound levels reached 70–78 dB(A). These estimates were more than double-fold (factor 2.4) in comparison to the urban residential background. At a 50 m distance off the road particle number concentrations were decaying to about 50% of the initial value. The measurements were characterised by close spatial correlation between total particle number concentration and ambient noise with correlation coefficients of up to r = 0.74. However, during one measurement day coupling between both quantities was weak due to higher turbulent mixing within the canopy layer and a change in ambient wind directions. Enhanced dilution of particle emission from road traffic by turbulent mixing and ‘decoupling’ from the influence of road traffic are believed to be responsible.     

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.     

On the escape of pollutants from urban street canyons

Pollutant transport from urban street canyons is numerically investigated using a two-dimensional flow and dispersion model. The ambient wind blows perpendicular to the street and passive pollutants are released at the street level. Results from the control experiment with a street aspect ratio of 1 show that at the roof level of the street canyon, the vertical turbulent flux of pollutants is upward everywhere and the vertical flux of pollutants by mean flow is upward or downward. The horizontally integrated vertical flux of pollutants by mean flow at the roof level of the street canyon is downward and its magnitude is much smaller than that by turbulent process. These results indicate that pollutants escape from the street canyon mainly by turbulent process and that the net effect of mean flow is to make some escaped pollutants reenter the street canyon. Further experiments with different inflow turbulence intensities, inflow wind speeds, and street aspect ratio confirm the findings from the control experiment. In the case of two isolated buildings, the horizontally integrated vertical flux of pollutants by mean flow is upward due to flow separation but the other main results are the same as those from the control experiment.     

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.     

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.     

Carbon dioxide fluxes over an urban park area

From September 2006 to October 2007 turbulent fluxes of carbon dioxide were measured at an urban tower station (26 m above ground level, z/z_h = 1.73) in Essen, Germany, using the eddy covariance technique. The site was located at the border between a public park area (70 ha) in the south–west of the station and suburban/urban residential as well as light commercial areas in the north and east of the tower. Depending on the land-use two different sectors (park and urban) were identified showing distinct differences in the temporal evolution of the surface-atmosphere exchange of CO₂. While urban fluxes appear to be governed by anthropogenic emissions from domestic heating and traffic (average flux 9.3 μmol m^?2 s^?1), the exchange of CO₂ was steered by biological processes when the park contributed to the flux footprint. The diurnal course during the vegetation period exhibited negative daytime fluxes up to ?10 μmol m^?2 s^?1 on average in summer. Nevertheless, with a mean of 0.8 μmol m^?2 s^?1 park sector fluxes were slightly positive, thus no net carbon uptake by the surface occurred throughout the year.In order to sum the transport of CO₂ a gap-filling procedure was performed by means of artificial neural network generalisation. Using additional meteorological inputs the daily exchange of CO₂ was reproduced using radial basis function networks (RBF). The resulting yearly sum of 6031 g m^?2 a^?1 indicates the entire study site to be a considerable source of CO₂.     

孤立与非孤立城市街道峡谷内污染物扩散

通过求解二维不可压N-S方程、k-ε方程及污染物对流扩散方程,模拟了孤立街道峡谷与非孤立街道峡谷内的流场及交通污染物浓度场.计算结果与风洞试验结果总体趋势一致.非孤立街道峡谷内污染物壁面浓度要大于孤立街道峡谷内的壁面浓度.通过计算街道峡谷建筑屋顶高度处的垂直方向污染物通量,说明了湍流扩散是污染物扩散出街道峡谷的主要原因,其污染物通量总为正,而平均流通量可以为负.非孤立街道峡谷由于平均流流动和湍流流动的总扩散通量减少,造成污染物在街道峡谷内集聚,从而理论上解释了非孤立街道峡谷与孤立街道峡谷污染扩散的差别.     

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.     

Assessment of natural components of PM10 at UK urban and rural sites

Proposals from the European Commission have raised the possibility that Member States may be able to subtract the concentrations of natural components of airborne particulate matter from measured concentrations when evaluating compliance with EU Limit Values. By applying the pragmatic mass closure model [Harrison et al., 2003. A pragmatic mass closure model for airborne particulate matter at urban background and roadside sites. Atmospheric Environment 37, 4927–4933] to chemical composition data for PM₁₀, it has been possible to estimate the concentrations of natural sea salt, strongly bound water and secondary organic carbon (which is assumed wholly biogenic) to the measured mass of PM₁₀. Because of the difficulty in distinguishing between natural and anthropogenic crustal dusts, the contribution of natural windblown dust and soil has not been accounted for. When the natural components are estimated for two urban and one rural site in the UK, the long-term mean PM₁₀ concentration is reduced by between 5.2 and 7.3 μg m⁻³. The number of exceedences of the 50 μg m⁻³ 24-h limit value falls dramatically from 54 to 21 (from a total of 291 days) at an urban street canyon site, 7 to 3 (n=292 days) at an urban background site and from 8 to 0 (n=241 days) at a rural site when using gravimetric PM₁₀ concentrations. The calculations have also been performed using PM₁₀ concentrations measured by TEOM increased by a factor of 1.3 as recommended by the European Commission as an interim means of estimating gravimetric equivalency, and the number of exceedences of the 24-h limit value fell from 92 to 47 (from a total of 291 days) at the urban street canyon site, from 11 to 3 (n=292 days) at the urban background site and from 6 to 3 (n=241) at the rural site. Clearly, therefore, application of this proposed measure would make a very major difference to the likelihood of compliance or otherwise with the 24-h limit value for PM₁₀.     

On the influence of the urban roughness sublayer on turbulence and dispersion

The concept of the urban roughness sublayer is discussed and this lowest atmospheric layer over a rough surface is shown to have a non-negligible vertical extension over typical urban surfaces. The existing knowledge on the turbulence and flow structure within an urban roughness sublayer is reviewed, focusing on the height dependence of turbulent fluxes and a scaling approach for turbulence statistics, such as velocity variances, in the above-roof part of the roughness sublayer. Finally, the implication of this turbulence and flow structure upon dispersion characteristics is investigated. The most prominent difference of explicitly taking into account the roughness sublayer in a dispersion simulation (as compared to assuming a `constant flux layer') is a clearly enhanced ground level concentration far downwind from the source. For the example of a tracer release experiment over a (sub) urban surface (Copenhagen) it is shown that introducing the roughness sublayer clearly improves the model performance.     

Computational Fluid Dynamics (CFD) Simulations on the Effect of Rough Surface on Atmospheric Turbulence Flow Above Hilly Terrain Shapes

The behavioral distribution of the atmospheric turbulence flow over the terrain with changes in a rough surface has become one of the most important topics of air pollution research, among such other topics as transportation and dispersion pollutants. In this study, a computational model on atmospheric turbulence flow over a terrain hill shaped with rough surface was investigated under neutral atmospheric conditions. The flow was assumed to be 2D and modeled using computational fluid dynamics (CFD) models, which were numerically solved using Reynolds-averaged Navier-Stokes equations. Rough surface conditions were modeled using a number of windbreak fences regularly spaced on the hill. The mean velocity and turbulent structures such as turbulence intensity and turbulent kinetic energy were investigated in the upwind and downwind regions over the hill, and the numerical models were validated against the wind-tunnel results to optimize the turbulence model. The computational results agreed well with the results obtained from the wind tunnel experiments. The computational results indicate that the mean velocity was observed to increase dramatically around the crest of the upwind slope of the hill. A thick internal boundary layer was observed with a fence on the crest and downwind region of the hill. The reversed flow and recirculation zone were formed in the wake region behind the hill. It was thus determined that turbulent kinetic energy decreases as the mean velocity increases.     

Observation of SO2 dry deposition velocity at a high elevation flux tower over an evergreen broadleaf forest in Central Taiwan

A 60-m flux tower was built on a 2100 m mountain for the measurement of the air pollutant concentration and the evaluation of dry deposition velocity in Central Taiwan. The tower was constructed in an evergreen broadleaf forest, which is the dominant species of forest in the world. Multiple-level SO₂ concentrations and meteorological variables at the site were measured from February to April 2008. The results showed that the mean dry deposition velocities of SO₂ were 0.61 cm s^?1 during daytime and 0.27 cm s^?1 during nighttime. From the comparison of the monthly data, a tendency was observed that the dry deposition velocity increases with LAI and solar radiation. Furthermore, it was observed that the deposition velocity was larger over wet canopy than over dry canopy, and that higher deposition velocities in the wet season were mainly caused by non-stomatal uptake of wet canopy. Over wet canopy, the mean dry deposition velocities of SO₂ were estimated to be 0.83 cm s^?1 during daytime and 0.47 cm s^?1 during nighttime; and 0.44 cm s^?1 during daytime and 0.19 cm s^?1 during nighttime over dry canopy. There is good agreement between the results of this study and those in other studies and the predictions of Zhang et al. (2003a). The medians (geometric means) of derived r_c during daytime are 233 (266) m s^?1 over dry canopy and 147 (146) m s^?1 over wet canopy. It was found that solar radiation is the critical important meteorological variable determining stomatal resistance during daytime. For non-stomatal resistance, clear dependencies were observed on the friction velocity and relative humidity.     

Characterization of PM2.5-bound polycyclic aromatic hydrocarbons in Atlanta

Twenty-four hour PM_2.5 samples from a rural site, an urban site, and a suburban site (next to a major highway) in the metropolitan Atlanta area in December 2003 and June 2004 were analyzed for 19 polycyclic aromatic hydrocarbons (PAH). Extraction of the air samples was conducted using an accelerated solvent extraction method followed by isotope dilution gas chromatography/mass spectrometry determination. Distinct seasonal variations were observed in total PAH concentration (i.e. significantly higher concentrations in December than in June). Mean concentrations for total particulate PAHs in December were 3.16, 4.13, and 3.40 ng m^?3 for the urban, suburban and rural sites, respectively, compared with 0.60, 0.74, and 0.24 ng m^?3 in June. Overall, the suburban site, which is impacted by a nearby major highway, had higher PAH concentration than did the urban site. Total PAH concentrations were found to be well correlated with PM_2.5, organic carbon (OC), and elemental carbon (EC) in both months (r² = 0.36–0.78, p < 0.05), although the slopes from the two months were different. PAHs represented on average 0.006% of total PM_2.5 mass and 0.017% of OC in June, compared with 0.033% of total PM_2.5 and 0.14% of OC in December. Total PAH concentrations were also correlated with potassium ion (r² = 0.39, p = 0.014) in December, but not in June, suggesting that in winter biomass burning can potentially be an important source for particulate PAH. Retene was found at a higher median air concentration at the rural site than at the urban and suburban sites—unlike the rest of the PAHs, which were found at lower levels at the rural site. Retene also had a larger seasonal difference and had the weakest correlation with the rest of the PAHs measured, suggesting that retene, in particular, might be associated with biomass burning.     

Chemical transport models

Background, aim, and scope  Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. Materials and methods  The upward mixing rate forming the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. Results  This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. Discussion  In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants calculated with chemical and air-quality models. In addition, this scheme was also compared with a commonly used, local, eddy-diffusivity scheme. Simulated concentrations of NO₂ by the COM scheme and new parameterization of the in-canopy resistance are closer to the observations when compared to those obtained from using the local eddy-diffusivity scheme. Conclusions  Concentrations calculated with the COM scheme and new parameterization of in-canopy resistance, are in general higher and closer to the observations than those obtained by the local, eddy-diffusivity scheme (on the order of 15–22%). Recommendations and perspectives  To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO₂) were compared for the years 1999 and 2002. The comparison was made for the entire domain used in simulations performed by the chemical European Monitoring and Evaluation Program Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.     

Seasonal ozone behavior along an elevation gradient in the Colorado Front Range Mountains

Ambient surface ozone was monitored for one year at a series of seven sites along an elevation gradient from 1600 m to 3500 m above sea level (ASL) in Boulder County, Colorado. Spatial variability of ozone, quantified as the root mean squared deviation of hourly ozone per kilometer horizontal separation, decreased with elevation and distance from local sources, validating the assumption that (except at the City of Boulder (BO) site) the results of the study are representative of the Colorado Front Range. The northern hemisphere (NH) tropospheric spring ozone peak was clearly apparent in late April and early May and affected ozone at all elevations. Ozone consistently increased with elevation during winter, with a mean monthly rate of 1.5 ppbv per 100 m elevation. In summer, this monotonic increase in ozone with elevation was not observed; instead mean monthly ozone increased in two steps, by ～15 ppbv between 1610 m and 1940 m ASL and by ～10 ppbv between 3350 m and 3530 m ASL to a maximum of 60 ppbv. The amplitude of the diurnal ozone cycle decreased with increasing elevation. Average summertime diurnal swings in ozone concentration had a magnitude of 29 ppbv at 1610 m ASL, and 7–16 ppbv at the mid-elevation sites. In winter a diurnal cycle was observed only at the BO site, ozone concentrations at the remaining six locations changed on a multi-day timescale, indicating regional background behavior as the primary factor for wintertime ozone. Even the highest elevation site was influenced by transported urban air pollution in summer, indicated by the average 5 ppbv diurnal increase in ozone. Ozone exposure at the mid- to high-elevation sites in many instances approached and exceeded the 8-h National Ambient Air Quality Standard of 75 ppbv. The elevated ozone levels along this transect were interpreted to be caused by the confounding effects of the high elevation of these sites, increased ozone in long-range transported air, and anthropogenic ozone production in air transported from the nearby urban and suburban areas east of the Colorado Front Range Mountains.     

Large-eddy simulation for flow and dispersion in urban streets

Large-eddy simulations (LES) with our recently developed inflow approach (Xie, Z.-T., Castro, I.P., 2008a. Efficient generation of inflow conditions for large-eddy simulation of street-scale flows. Flow Turbul. Combust., vol. 81(3), pp. 449–470.) have been used for flow and dispersion within a genuine city area – the DAPPLE site, located at the intersection of Marylebone Rd and Gloucester Pl in Central London. Numerical results up to second-order statistics are reported for a computational domain of 1.2 km (streamwise) × 0.8 km (lateral) × 0.2 km (in full scale), with a resolution down to approximately one meter in space and one second in time. They are in reasonable agreement with the experimental data. Such a comprehensive urban geometry is often, as here, composed of staggered, aligned, square arrays of blocks with non-uniform height and non-uniform base, street canyons and intersections. Both the integrative and local effect of flow and dispersion to these geometrical patterns were investigated. For example, it was found that the peaks of spatially averaged u_rms, v_rms, w_rms and <u′w′> occurred neither at the mean height nor at the maximum height, but at the height of large and tall buildings. It was also found that the mean and fluctuating concentrations in the near-source field is highly dependent on the source location and the local geometry pattern, whereas in the far field (e.g. >0.1 km) they are not. In summary, it is demonstrated that full-scale resolution of around one meter is sufficient to yield accurate prediction of the flow and mean dispersion characteristics and to provide reasonable estimation of concentration fluctuations.