Variability of physical meteorology in urban areas at different scales: implications for air quality

Air quality in cities is in ﬂ uenced not only by emissions and chemical transformations but also by the physical state of the atmosphere which varies both temporally and spatially. Increasingly, tall buildings (TB) are common features of the urban landscape, yet their impact on urban air ﬂ ow and dispersion is not well understood, and their e ﬀ ects are not appropriately captured in parameterisation schemes. Here, hardware models of areas within two global mega-cities (London and Beijing) are used to analyse the impact of TB on ﬂ ow and transport in isolated and cluster settings. Results show that TB generate strong updrafts and downdrafts that a ﬀ ect street-level ﬂ ow ﬁ elds. Velocity di ﬀ erences do not decay monotonically with distance from the TB, especially in the near-wake region where the ﬂ ow is characterised by recirculating winds and jets. Lateral distance from an isolated TB centreline is crucial, and ﬂ ow is still strongly impacted at longitudinal distances of several TB heights. Evaluation of a wake-ﬂ ow scheme (ADMS – Build) in the isolated TB case indicates important characteristics are not captured. There is better agreement for a slender, shorter TB than a taller non-cuboidal TB. Better prediction of ﬂ ow occurs horizontally further away and vertically further from the surface. TB clusters modify the shape of pollutant plumes. Strong updrafts generated by the overlapping wakes of TB clusters lift pollutants out of the canopy, causing a much deeper tracer plume in the lee of the cluster, and an elevated plume centreline with maximum concentrations around the TB mean height. Enhanced vertical spread of the pollutants in the near-wake of the cluster results in overall lower maximum concentrations, but higher concentrations above the mean TB height. These results have important implications for interpreting observations in areas with TB. Using real world ceilometer


Introduction
Air quality in cities is impacted not only by emissions and chemical transformations but also by the physical state of the atmosphere.In this paper we focus on the latter, particularly considering the impacts of tall buildings (TB).We consider both isolated tall buildings, which are common in Europe, and clusters of tall buildings, which occur more frequently in Asia and the Americas.
Critical in our consideration is the vertical spatial structure of the urban atmosphere.We identify three scales with the following horizontal and vertical dimensions [1][2][3] (Fig. 1): (1) Micro-scalebuildings to streets impact the urban canopy layer (UCL) and the roughness sub-layer (RSL).
(2) Local-scaleneighbourhoods (areas with common repeated characteristics, e.g.low density residential, high density commercial) inuence the inertial sub-layer (ISL) or constant ux layer above the RSL.
Fig. 1 Schematic of key urban climate scales.These vary with height of the roughness elements (H) and their spacing (W), the order of magnitude of the vertical extent and two proxy methods to determine the inertial sub-layer (ISL).UCLurban canopy layer, RSLroughness sub-layer, SLsurface layer, UBLurban boundary layer (note: not to scale).Inset (bottom) shows street-scale and tall-building flow regimes. 373) Meso-scalea city impacts the urban boundary layer (UBL), the depth of which changes with both time of day and season.
Flow, turbulence and dispersion vary rapidly in three dimensions based on the positioning of individual roughness elements (RE; e.g.buildings, trees) in the UCL-RSL. 4The UCL extends to the heights of the RE, but the vertical extent of the RSL varies with both the RE heights (mean, maximum, variability) and the spacing between them [5][6][7][8][9] (Fig. 1, purple equations).The UCL, RSL and ISL combine to make the surface layer (SL), or lowest $10% of the UBL (Fig. 1).
Observations and modelling within cities need to consider this structure of the urban atmosphere.Basic micrometeorological theory is developed for the ISL.For example, the logarithmic wind prole and the Monin-Obukhov similarity theory both apply in the ISL but not within the UCL-RSL. 2,10Thus, our understanding and the ability to predict the transport close to the surface (UCL + RSL) is poor.As bluff body RE (e.g.buildings) are not porous, unlike trees, they can channel ow within the UCL and inuence building ventilation. 11,12To enhance our theoretical understanding of processes and/or to undertake numerical predictions we typically simplify the urban form to capture what are thought to be the most relevant features. 4,13,14For example, buildings may be treated as having smooth walls, 15 rather than accounting for variations caused by the presence of windows, doors, balconies (etc.) or even non-cuboid shapes. 16Similarly, buildings may be treated as having at roofs that are all of the same height 13,17 rather than being sloped at different angles and/or of variable heights. 16iven issues of access, observations in the real world are oen undertaken in the UCL.These micro-scale measurements are challenging to interpret and generalise because small changes in sensor location can result in large changes in the observed ow and radiative characteristics because of the three-dimensional scaling well below the typical RE height. 2,15,18,19ndertaking measurements in the ISL, above the RSL (Fig. 1), becomes more representative (theoretically) because of the blended response from the surface so that individual micro-scale features no longer dominate, and the net transport is referred to as having one-dimensional scaling in extensive homogeneous areas.With sufficient homogeneous fetch, it is possible to use the eddy covariance technique to directly measure momentum, heat and mass uxes. 18Urban land surface models are ideally evaluated using observations in the ISL. 20,21The patchwork of neighbourhoods across a city create the UBL.Despite its importance to pollutant concentrations and hence air quality, there are few continuous observations of the UBL.
In this paper we explore the impacts of tall buildings on the urban atmosphere across these three scales (UCL + RSL, ISL, UBL) using observations undertaken in a wind tunnel, a water ume and in the real world.The advection and diffusion of pollutants in urban environments is governed by complex, building-induced wind elds (Fig. 1).Combined effects of channelled or recirculating ow in street canyons, ow branching at intersections and the turbulent vertical exchange at roof-level between the urban canopy and the boundary layer, contribute to the strong spatio-temporal variability of street-level concentrations and local human exposure levels. 22,23perational urban dispersion and air quality models need to parameterise building-induced dispersion features, as individual buildings are not resolved in the simulations.For example, the ADMS-Urban air quality model 24 combines a Gaussian plume dispersion approach with the OSPM street canyon model 25 to capture effects such as pollutant accumulation in recirculating ow or pollutant channelling through streets.The Lagrangian dispersion model QUIC-plume's 26 built-in ow-solver uses mass-conservation to give 3D building-resolved mean ow elds.SIRANE's 27 street-network approach parameterises dispersion through streets and intersections using bulk advection and vertical exchange velocities per street unit.
Typical operational settings rarely account for building effects down to the street-scale level.Urban applications of the widely used AERMOD 28 air quality model or the NAME 29 dispersion model, for example, use bulk parameterisations to reect changes in surface roughness and heat-ux characteristics over cities, but do not model canopy-layer ow effects.Both ADMS-Urban and AERMOD have built-in modules to represent wake-ow effects for individual (isolated) buildings (ADMS-Build, 30 AERMOD-PRIME 31 ).However, these modules are not routinely used in the context of urban air quality modelling.
Street-canyon ow parameterisations oen assume an equal-height canyon with the building height (H) to street width (W) ratio creating different ow regimes (0.35 # H/W, 0.35 < H/W < 0.65, H/W $ 0.65; isolated < wake interference < skimming). 23][34] Tall-building wake effects, i.e. modications of horizontal and vertical advection and diffusion characteristics on the leeward side of the buildings, can be important over large distances and hence can be crucial in urban air quality modelling on and across different spatial scales.However, little is known about the exact nature of ow and pollutant dispersion modications, especially in the case of clusters of tall buildings.This currently limits the capability of air-quality and emergency-response dispersion models to represent such processes adequately.
In this study, at the micro to local-scale we explore morphologically complex areas in London and Beijing using laboratory observations to assess the effect of isolated tall buildings and tall-building clusters on ow and dispersion elds, respectively.
At the local to city scale we examine the mixed-layer height (MLH) evolution.A simple empirical model, with parameters informed by London observations, 35,36 is evaluated using atmospheric prole measurements in Beijing and Paris.

Study areas
To explore the physical meteorological processes, we use real world observations and hardware models for parts of real cities.Naturally, the latter are simplications as they do not have all the real-world details.For example, the roofs and walls are at/smooth rather than sloped/crenulated, there are no anthropogenic heat emissions, nor radiative exchanges.The studies take place in three megacities: Beijing, London, and Paris.
In Beijing, the observations are taken in a residential-commercial area of the Haidian District between the third and fourth north ring roads.This area, typical of central Beijing, is the urban site (39 58 0 33 00 N, 116 22 0 41 00 E) of the Atmospheric Pollution and Human Health in a Chinese Mega-city (APHH-Beijing) 38 project.The Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences 325 m meteorological tower is located here. 38or the local-scale study, a hardware model is built of an area ($580 m diameter) north of the IAP tower with a cluster of 14 tall buildings surrounded by a low-rise canopy (Fig. 2b).The MLH is determined from an instrument located near the base of the IAP tower.
Similarly, in London the micro-to-local scale processes are studied in a small area within the London Borough of Southwark (Fig. 2a).The hardware model of a $700 m diameter region around St George's Circus (51.4987N, 0.1048 W) has 148 buildings that are mostly (49%) low-rise residential (building height H # 12 m) and mid-rise (12 to 24 m) office and commercial (39%).Only 17 of the buildings have H > 24 m (3, H > 32 m).
The MLH model coefficients used were derived 36 from observations taken in central London, where surface roughness is high. 10The sensor was operated at In Paris, only the larger scale MLH analyses are undertaken, using a sensor located at SIRTA (48.713N and 2.208 E) in a suburban location 20 km south of central Paris. 39

Wind-tunnel experiments of isolated tall buildings
An open-return boundary-layer wind-tunnel (test-section length: 20 m; cross-section: 3.5 m Â 1.5 m) was used for ow experiments in a hardware scale model of buildings for the area in south London (Fig. 2a, Section 2.1).The $700 m extent of the study domain has two isolated TB with heights of H TB ¼ H 81 ¼ 81 m (hereaer: T 81 ) and H 134 ¼ 134.3 m (T 134 ) surrounded by low-rise buildings (average H: $13.4 m) in a mixed residential-commercial neighbourhood (Fig. 2a).The hardware model is built at a scale of 1 : 200.T 134 with a triangular footprint is located at the southeastern edge of the domain (Fig. 2a) and the hexagonal T 81 is located close to the model centre.Hence, the ow at T 134 is also affected by the roughness transition between the approach-ow and model. 40The depth of the approach-ow internal boundary layer was $1 m, i.e. similar to the height of T 134 .
Point-wise velocity measurements 41,42 obtained using a two-component laser Doppler anemometry system, 43 are analysed.In each experiment, neutral atmospheric stability conditions (i.e.thermal ow effects are negligible) are used.The wind-tunnel free-stream velocity U ref is steady at 2 m s À1 AE1% in all cases.Experimental design, measurement techniques and uncertainties are available. 40ere we analyse two congurations of the model: (i) Core and (ii) No Tall.In the Core case the hardware model includes both the tall buildings and the low-rise building canopy (Fig. 2a).Whereas in the No Tall case only the low-rise canopy buildings (H < 32 m) are present (i.e.there are no TB).
Analysis of the wind-tunnel observations for these two conditions allows the inuence of the presence of the TB to be assessed.The buildings in this case have at roofs and smooth walls, but are close to the 'true' general footprint shape (i.e.not rectangular in cross-section).The street orientations and building arrangements are realistic, but other roughness elements (e.g.trees, bus shelters, lamp posts) are missing.Thus, the model can be regarded as a realistic characterisation of the actual 'chaotic' building forms found in a small area of London.The combination of large building scale (1 : 200), the tallest building (T 134 ) and wind tunnel dimensions allow proles throughout the UCL, but it is not possible to undertake measurements well beyond the RSL of the full model.

Water-ume experiments with a cluster of tall buildings
A closed-loop water ume (test-section length: 6.25 m; width: 1.2 m; constant water depth: 0.6 m) with particle image velocimetry (PIV) 44 and planar laserinduced orescence (PLIF) 45  The dispersion of a passive scalar tracer is determined using an aqueous solution of uorescent dye (Rhodamine 6G) released 5 mm upstream of the model domain from a 3 mm diameter ground source with a continuous source strength of 500 mm 3 s À1 .The tracer concentration at the source, C s , was 0.3 mg l À1 .The PLIF calibration procedure and general aspects of measurement uncertainty of the technique are discussed. 46n the $580 m extent model domain (Fig. 2b) there is a cluster of 14 tall buildings (maximum/minimum heights: H 72 ¼ 72 m/H 48 ¼ 48 m) surrounded by a low-rise canopy (typical building heights between 4.8 m and 14.4 m).The hardware model has a more detailed representation of the building cross-section than the London case (Fig. 2b cf.Fig. 2a), but still has at roofs.Some simplication is required because of the small scale of the models.The two scales used (1 : 2400 and 1 : 4800) 47 allow analysis of both the RSL and ISL characteristics, but limited detail within the UCL (Fig. 1).The depth of the approach-ow boundary layer is $85 mm in the model scale ($203 m at a scale of 1 : 2400; i.e. $3H 72 ).
To investigate the inuence of the Beijing TB cluster, multiple model congurations are used: (i) low-rise canopy without the tall-building cluster (Low), (ii) tall-building cluster with low-rise canopy (Tall + Low; Fig. 2b) and tall-building cluster in isolation without the low-rise buildings (Tall).

ADMS-Build wake model
The ADMS-Build wake model 30,48 parameterises the velocity-decit in the wake of an isolated, cuboidal-shaped building through a 3D extension of the 2D Counihan et al. 49 constant eddy-viscosity wake theory.ADMS-Build provides an analytical expression for the ow eld in the momentum decit region of the main building wake (see Hertwig et al. 40 eqn 3-5 and Fig. 12).
The beginning of the main wake is determined by modelling the extent of the recirculation (cavity) zone that develops on the lee of the building (see inset Fig. 1) as a function of height, crosswind width and along-wind building length.In this case, the inow to the model (i.e.ow upwind of the building) is based on wind tunnel ow conditions for the case without tall buildings (No Tall).Implementation of the model and further evaluation results are given elsewhere. 40

A simple urban boundary layer height model
The average MLH is estimated in Beijing and Paris using a simple parameterisation with coefficients derived in London. 36The modelled MLH varies both diurnally and seasonally with differences in cloud cover.The cloud conditions considered are (Table 1): clear nights followed by cumulus clouds during the day (clearCu); same conditions in both the day and nightcloudless (clear), cumulus (Cu) or stratus (St); or when the MLH is below the cloud base height (CBH).Cloud conditions cause variations in minimum MLH (z min ), the morning transition (MT) growth rate (GR MT ) and the nocturnal decay rate (DR noc ; Table 1).The MLH above ground level (agl) is approximated by calculations at six key times through the day (Table 1).These times vary with sunrise (SR) and sunset (SS) to identify the timing of minimum and maximum depth of the MLH, and the beginning and end of both the morning and the evening transition periods.
Automatic lidar and ceilometer (ALC) measurements of attenuated backscatter allow the mixed layer height to be determined from observations.Here Vaisala CL31 observations are used.A background correction 50 is undertaken prior to using the CABAM (characterising the atmospheric boundary layer based on ALC measurements) algorithm. 35The algorithm considers precipitation, cloud cover and cloud type, all obtained from the ALC measurements.
Instruments installed in Beijing, 38 very close to the TB cluster (Fig. 2b), and in suburban Paris 39 (Section 2.1; https://www.lmd.polytechnique.fr/sirta/mld/)are used to obtain the attenuated backscatter data.The APHH-Beijing 38 data available for analysis are from November 2016 to June 2017.Although observations occur in all four seasons, both summer and autumn are incomplete.For Paris, data are analysed for a year (2018). 51During the analysis, the diurnal variability of the MLH is assessed using the median and the interquartile range (IQR).The four seasons are dened by month: winter (DJF, December, January, February), spring (MAM), summer (JJA), and autumn (SON).

Tall-building wakes and model comparison in London
Prior analyses of tall-building wakes in the London domain 40 (Fig. 2a) found wind elds above the low-rise canopy to be impacted for downwind distances of at least 5H TB of the towers.As numerical weather prediction (NWP) models are used to drive dispersion and air quality models, their horizontal resolution is important.
With the resolution of next generation NWP models expected to be of the order of 100 m, 52 tall-building wakes will span several grid cells.We use a Cartesian coordinate system to present the results.The longitudinal axis (x direction; u velocity component) is aligned with the mean forcing wind direction (i.e.along-wind).The y axis (v velocity) is the lateral (transverse) direction.The z axis (w velocity) indicates the vertical component.Time-averaged velocities (e.g.U) are written in uppercase and are scaled by the free-stream velocity (U ref ).
The logarithmic law allows the wind prole to be modelled in the ISL. 53lthough not directly applicable within the RSL, 10 it is oen used above the UCL in air-quality models like ADMS-Urban, SIRANE or QUIC-plume.Analysis of measurements across the Core conguration (Fig. 2a) provides the ensemble mean vertical proles (z) of the longitudinal (along-wind) velocity (U/U ref ) at increasing radial distances (r i ) from the TB (Fig. 3).These results show the spatial variability of the ow in the presence of the tall buildings and the extent of wake effects with increasing distance and number of measurement points included in the ensemble average.Ensemble-averaged wind proles of the No Tall case (i.e. when T 81 and T 134 are removed) serve as a reference, as the proles follow the logarithmic law from above the height of the roughness length plus displacement height (z 0 + z d ).However, for both T 81 (Fig. 3a and c) and T 134 (Fig. 3b and d) the Core proles show very noticeable deviations from the logarithmic law for all radii.
The ow proles averaged over different radial distances only converge back to the reference state well above roof-level (z ¼ H 81 ) for the slender T 81 tower near the model centre.The taller and wider T 134 building has a much broader wake, 40 with lateral fanning enhanced by the shape and orientation of the tower relative to the south-easterly forcing wind direction (Fig. 3d).The differences from the No Tall reference are enhanced (cf.T 81 ) and remain large even when averaging over all available Core sites (note: no measurements are available above a height of H 134 ).
To explore further the impact of the T 134 wake (Fig. 3d) velocity differences (Core -No Tall) are analysed for both the mean along-wind velocity (DU/U ref ; Fig. 4a, c and e) and the mean vertical velocity (DW/U ref ; Fig. 4b, d and f) without the spatial averaging.The velocity differences are shown as a function of their straight-line distance from the centre of T 134 (d/H 134 ), with the lateral displacement of the sites (|y À y 134 |/ H 134 ) from the building centreline (y 134 ) indicated, at three heights: just above the low-rise canopy (z/H 134 ¼ 0.27; Fig. 4a and b), at close to half T 134 (z/H 134 ¼ 0.48; Fig. 4c and d) and near roof-level of T 134 (z/H 134 ¼ 0.91; Fig. 4e and f).
The strongest impact on the ow is measured directly behind the building, with a velocity decit DU of up to $1.2 times the reference velocity U ref (Fig. 4c) and residual updras (DW/U ref > 0) of up to $0.28 U ref (Fig. 4b) and downdras (DW/U ref < 0) of $0.22 U ref (Fig. 4f).The changes in the mean vertical ow over large distances are crucial, as in dispersion modelling it is typically assumed that on average W ¼ 0 m s À1 , i.e. the vertical pollutant exchange is taken to be governed by turbulent motions and not by the mean ow. 27This is clearly not the case in the presence of TB, which generate strong updras and downdras that affect street-level ow elds. 34The velocity differences (Fig. 4) do not decay monotonically with distance from the TB, especially in the near wake region where the ow is characterised by recirculating winds.Instead, the lateral distance of the measurement points from the building centreline is crucial.Despite being several building heights away (i.e.several 100 m), locations with small |y À y 134 |/H 134 are strongly impacted by the TB.
Some urban air-quality models (e.g.ADMS-Urban) have ow modules to account for the impact of isolated buildings on pollutant dispersion.Wake models (e.g.ADMS-Build 3 ) represent ow around cuboidal buildings of modest height to length (or width) ratios, without accounting for the impact of neighbouring low-rise buildings.We evaluate ADMS-Build (Section 2.4) at the London site using our measurements in the main wake (i.e.behind the recirculation zone of the near wake, starting $51 m downwind of T 81 ; $161 m downwind of T 134 ; 40 inset Fig. 1).Data are analysed for all available measurement heights, stratied by straight-line distance (d/H TB ) between the tall building and the comparison sites (Fig. 5).
Overall, the agreement between ADMS-Build and wind tunnel is better for the slender, shorter T 81 (Fig. 5a) than for the larger and wider T 134 (Fig. 5b).The percentage of locations where the model and measurement agree within a tolerance threshold of AE0.1 U/U ref for either T 81 or T 134 clearly depends on the respective distance from the tower.The ow at sites closer to the towers is less well predicted (cf.sites at greater distances).Model performance improves with height above the urban canopy (i.e.larger U/U ref conditions).This is partially explained by the neglecting of aerodynamic effects from low-rise buildings surrounding the towers within ADMS-Build.For some sites this causes an overprediction of velocities closer to the ground.The wake-ow close to T 134 (located at the edge of the wind-tunnel model; Fig. 2a) is also affected by the roughness transition between the upwind boundary-layer and the building model, which is not accounted for in ADMS-Build. 40Such changes in roughness characteristics (e.g.induced by the local building morphology of a neighbourhood) occur frequently in mega-cities (e.g.central business districts, residential quarters, urban parks).

Beijing: plume characteristics in a cluster of tall buildings (TB)
The effect of tall buildings on the dispersion of pollutants continuously emitted from an upwind ground source is analysed in the TB cluster in Beijing (Fig. 2b and  6).For reference, we compare the pollutant dispersion to an open area (i.e.without buildings, Fig. 6a) to identify changes in mean plume characteristics for dispersion through a low-rise building canopy (Low conguration; Fig. 6b) and through a high-rise building cluster surrounded by the low-rise buildings (Tall + Low, Fig. 2b, 6c).The time-averaged concentrations (C) are scaled by the concentration of the tracer solution at the release point (C s ).The elevation of the plume centreline is used to analyse the impact of the TB cluster.This is measured by the height of the maximum tracer concentration (C max /C s ) along the longitudinal direction.The vertical plume half-width is the distance over which C max /C s has dropped to half its local maximum value vertically.
Compared to the open-reference, the elevation of the plume centreline is noticeably enhanced in the presence of the low-rise canopy (Fig. 6a and b).An initial upward momentum of pollutants when released (i.e. at the source) is observed (Fig. 6a) causing C max to be initially away from the ground.However, in the presence of the building canopy, this is a secondary effect on the plume dispersion.The increased vertical transport and mixing of pollutants in the building canopy contribute to a faster decay of concentrations with increasing The vertical plume half-width is noticeably increased at the far end of the measurement domain ($8.5 building heights H 72 away from the source) in the Low conguration (cf. to the open case).Including the Tall-building cluster in the centre of the model domain leads to strong changes in the plume shape downwind of x/H 72 ¼ 2.5 (Fig. 6c).The strong updras generated in the overlapping wakes of the tall buildings li pollutants out of the canopy, resulting in a much deeper tracer plume on the leeward side of the cluster (larger vertical plume halfwidth) and an elevated plume centreline (maximum concentrations just below z ¼ H 72 ).Note that the rst building impacted by the plume is located at x/H 72 ¼ 3.25 (i.e. the TB at x/H 72 ¼ 2.75 is in front of the measurement plane, camera perspective).
The height of the plume envelope increases longitudinally (Fig. 6c; increasing from z/H 72 ¼ 1 to 2 over x/H 72 ¼ 3.25 to 8.5) and correlates with the region of maximum vertical wind shear (and maximum velocity variances; not shown).The enhanced mixing and vertical spread of the pollutants in the wake of the cluster results in overall lower maximum concentrations.Furthermore, on the windward face of the TB below the region of strong upward ow at roof-level, there is a region of enhanced downward motion (downwash; Fig. 1) and outow near the building's base resulting in a further lateral distribution of pollutants.A large portion of the pollutant plume behind the TB cluster is transported above the street canyons of the low-rise buildings.Hence, higher concentrations are observed above z/H 72 ¼ 1 compared to the open-terrain and Low cases at the same distances from the source.
To study the wake region behind the cluster over larger distances, measurements in a 1 : 4800 scale model of the Tall conguration (no low-rise buildings) Fig. 7 Vertical profiles of normalised concentrations (C/C max ) for the Tall case of the Beijing neighbourhood (q ¼ 0 ; arrow in inset of building geometry; scale 1 : 4800) showing measurements (circles) and Gaussian plume model fits (dashed lines) at two locations upwind of the building cluster (blue, cyan; modelled for a ground source release) and at three locations in the wake region (green, yellow, red; modelled for a release from an elevated source).The x/H 72 locations of the profiles are indicated in the inset showing the concentration plume (log 10 (C/C s ); contours as in Fig. 6).For clarity, profiles show only every third measurement point and measurements closest to the ground (with larger uncertainties) are removed.are analysed.Vertical proles of mean concentrations normalised by the local maximum concentration at various up-and down-wind distances to the tall building cluster show the change in plume shape (Fig. 7).While the region upwind of the buildings has maximum concentrations close to the ground and a rapid decay of concentrations with height, in the wake region the pollutants are distributed over a much larger vertical depth, with C max elevation growing with distance from the buildings.We compare the measured concentration proles with ts for a reected Gaussian plume model: 54 where z s is the height of the source and s z characterises the vertical plume spread (tted).
As the pollutants are released from a ground source (at x/H 72 ¼ 0), upwind of the buildings (i.e. over open terrain) the pollutant dispersion is modelled by setting z s ¼ 0 m.While the measured upwind proles decay less rapidly with height, overall their shape agrees with the Gaussian model (Fig. 7).Aer passing through the TB cluster, the plume shape in the wake is well approximated if a z s > 0 m is used, i.e. concentration proles in the wake of the tall buildings resemble those obtained as if an elevated source release had occurred.Here, the "virtual" source height z s /H 72 varies between 0.85 at x/H 72 ¼ 11 and 1.06 at x/H 72 ¼ 16 (Fig. 7).
Model development work supported by experimental evidence is needed for current (operational) air-quality and dispersion models to be able to capture such effects and make the models suitable for applications in cities with TB.

Mixed layer height (MLH) in Beijing and Paris
Here, the MLH is derived from observations of attenuated backscatter from aerosols and moisture.Hence, it is indicative of the depth of the atmospheric layer that the aerosols have mixed into, rather than the depth of the atmospheric layer that is undergoing active vertical motion (referred to as the mixing height). 55ayers of aerosols create gradients in the observed attenuated backscatter (here 10 m vertical and 15 s temporal resolution) and are used to identify the MLH. 35he observed diurnal maxima of MLH during the study years (Section 2.5) in Beijing and Paris are, for both cities, larger in spring and summer than in the autumn and winter (Fig. 8, row 3).With the warmer season peak daytime medians (thick lines) are similar to the 75 th percentile (thin upper lines) in autumn in both cities.
The simple MLH model (Table 1) applied in the two mega-cities uses parameters derived from six years of observations undertaken in central London. 36Synoptic conditions generally differ between the three cities, with a clear oceanic inuence in London, continental climate in Paris, and a mixture of orographic and maritime effects dominating Beijing. 38MLH growth is driven by the strength of the vertical motions and entrainment processes.The turbulent sensible heat (Q H ) as well as momentum exchange modify boundary layer dynamics within the synoptic context.The model accounts for differences in daylength (Section 2.5) as well as cloud cover and cloud type, but does not consider energy forcing differences associated with geographic location, surface cover, or human activities.As it effectively only has linear variations between the six key times within a day (Table 1, Fig. 8), the discussion of model performance should focus on these periods rather than other times.
On average the model (with London coefficients, Table 1) works surprisingly well in Beijing (Fig. 8, row 1) where parameterised diurnal patterns mostly fall within the inter-quartile range of the observations.For the London-Paris comparison (Fig. 8, row 2), the curves reveal an overprediction of MLH in winter and underprediction in summer.To evaluate the performance we compare the peak daytime values modelled with the corresponding observed MLH values (data point closest in time).Both the mean bias error (MBE) and the mean absolute error (MAE) are calculated.The overall MBE and MAE are greater for Paris (À230  1).The number of days (N) with data in each season varies.Observations (black, row 1 and 2, other colours in row 3) are Vaisala CL31 ceilometer attenuated backscatter data analysed with the CABAM algorithm (Section 2.5).Seasons are winter: DJF, spring: MAM, summer: JJA, autumn: SON.Note: y-axes are different for the two rows.

Faraday Discussions Paper
Faraday Discuss.This journal is © The Royal Society of Chemistry 2020 Given the lower latitudes (Section 2.1) of both cities (cf.London) the peak solar noon radiation on a clear day should be greater than in London.This could contribute to the MLH increment at the continental, suburban European site.The simple MLH parameterisation accounts for differences in cloud cover and cloud type, so that the greater probability of clear-sky days at the continental Paris site 57 does not explain the contrast to the London conditions.It is likely that the continental synoptic setting of Paris (e.g.absence of sea breeze effects) explains some of the larger growths rates and resulting greater MLH at this site during spring and summer (Fig. 8, row 3).A strong relation between peak MLH and peak incoming solar radiation has been found in Paris. 56Given clear skies in London are less frequent in summer than in Paris (12% versus 27%), 51 the London empirical data may not sample these synoptic conditions with the largest MLH as extensively as needed.
Potential increase in solar radiation might have less impact in the Asian megacity, given increased aerosols will reduce the shortwave receipt at the surface. 38,58,59Further, energy partitioning may restrict the growth of the mixed layer.In central London, the surface energy balance is dominated by Q H . [60][61][62] In Beijing, Q H is much greater in the summer at the IAP site than it is in suburban Meiyun. 63These two Beijing sites have very similar net all-wave radiation (Q*) but there are differences in outgoing longwave radiation characteristics, with the urban area having greater nocturnal emissions. 63The suburban site has larger latent heat uxes (Q E ) making the Bowen ratio (b ¼ Q H /Q E ) smaller.However, given the land cover characteristics, the urban IAP site has smaller summertime b than expected. 63This is thought to be associated with extensive external water use for irrigation of vegetation and for wetting roads to reduce airborne particles and cleaning. 63Thus, the summertime b values are much greater in central London than in urban Beijing which could mean that even if more energy were available from radiation, the MLH growth is still restricted.
In both Beijing and London, the MLH observations were taken near the centre of the mega-cities, whereas in Paris the site is in a suburban area some distance from the centre (Section 2.1).Urban density greatly affects anthropogenic heat ux (Q F ) emissions.This additional source of energy (i.e.][66][67] In colder climates the Q F emissions are oen greater in the winter than the summer 65,[68][69][70] because of building heating needs.Thus, winter Q F emissions provide a proportionally larger additional source of energy relative to Q* than in summer, given the shortwave radiation differences. 39As Q F is always a source of energy, unlike nocturnal Q*, it can have a large impact on the wintertime MLH.At night, this additional source of energy combined with storage heat ux (DQ S ) release and outgoing longwave radiation can contribute to Q H , thereby maintaining the MLH at a greater depth.Obviously, this is important as both winter and nocturnal periods are times when the MLH is lower (Fig. 8), hence any increase in boundary-layer volume will have proportionally larger impact on potential dilution/concentration of air pollutants and hence on local air quality.From machine learning analyses of the meteorological drivers of PM 1 variability at the Paris site, it is found that the impact of MLH up to about 900 m is particularly pronounced. 71uring winter, empirical coefficients (Table 1) derived from the central London observations appear to capture the MLH evolution at the central Beijing site well, while suburban Paris (Fig. 8) conditions are overestimated.This is likely explained by the lower Q F and DQ S at the suburban location providing less energy to maintain Q H and thereby a considerable MLH.MLH at SIRTA has been observed to be lower than in central Paris during winter and at night. 72verall, these results suggest the form of the empirical model is reasonable and further analyses of both ALC MLH data and surface energy balance partitioning with consideration of synoptic conditions would allow a simple, extremely fast, rst-order model to be developed with broader applicability.

Conclusions
Physical meteorology as well as chemical emissions and transformations need to be understood to interpret air quality changes and variations across mega-cities.Given that physical meteorological processes vary at micro, local and city scales (Fig. 1), awareness of these scales and their impact on atmospheric transport is important.Depending on the height of observations, physical and chemical characteristics of the atmosphere are inuenced more/less by individual roughness elements.Given the complexity of the urban canopy layer (UCL, Fig. 1), observations taken within this layer are very dependent on the exact measurement site.This creates a signicant challenge for those interested in air quality given that outdoor human exposure is most frequently at an order of 1-2 m agl.In mega-cities, there are increasing numbers of tall buildings (TB) in isolated and cluster settings.These expand the depth of both the UCL and the roughness sub-layer (RSL, Fig. 1), where atmospheric motions are more chaotic, and increase the inhomogeneity of the urban atmosphere, complicating the situation further.
Models that yield results rapidly enough to provide air quality forecasts, parameterise the urban physical meteorology.To date, these parameterisations have not considered TB adequately.Using hardware scale models that characterize real areas in two global mega-cities (London, Beijing) we explore the impact of isolated TB and a cluster of TB on ow and dispersion.The studies consider the UCL, RSL and ISL.Using real world observations in two global mega-cities (Beijing, Paris) we assess the UBL height modelled with parameters derived from a third mega-city (London).
Isolated TB wakes impact the ow above the background low-rise buildings for at least 5H TB downwind, with stronger impacts in regions close to the centre-line of the wake.Wind proles clearly still deviate from the logarithmic law at 4H TB radial distances from the TB.Flow proles only converge back to the reference state well above roof-level over the horizontal distances analysed in this study.A taller and wider TB, or a group of such buildings, has a much broader wake, varying with building shape and wind direction.
The strongest impact on the mean ow is evident directly behind the TB, based on the magnitude of the velocity decit and residual up-and down-dras.As dispersion models typically assume W ¼ 0 m s À1 with vertical pollutant exchange governed by turbulent motions and not by the mean ow, changes in mean vertical ow over distances of interest are crucial.
Velocity differences do not decay monotonically with distance from the TB, especially in the near-wake region where the ow is characterised by recirculation.Lateral distance from the TB centreline is crucial as strong effects persist close to the TB centreline even several H TB away (i.e.several 100 m).
Evaluation of ADMS-Build in the main TB wake considering straight-line distance (d/H TB ) found better agreement for a slender, shorter TB (cf.larger TB).Prediction of ow horizontally is better further away (cf.nearer); whereas the prediction is better higher above the background low-rise buildings while overprediction occurs closer to them given that they are unaccounted for by the model.
A TB cluster causes strong changes in the plume dispersion shape downwind.Strong updras generated in the overlapping wakes of the TB li pollutants out of the canopy.This results in a much deeper tracer plume in lee of the cluster (larger vertical plume half-width) and an elevated plume centreline (maximum concentrations just below H TB ).Enhanced mixing and vertical spread of the pollutants in the wake of the cluster result in overall lower maximum concentrations for a surface source, but higher concentrations above H TB (cf.open terrain and lowrise) at the same distance from the source.
A low-rise canopy (cf.open terrain) increases the vertical transport and mixing of pollutants in the building canopy.This contributes to a faster decay of concentrations with increasing distance from the release point.The along-wind decay is further enhanced by a stronger lateral plume spread from pollutant channelling and lateral plume branching in response to the local street topology.
The building cluster causes both vertical displacement of the maximum concentration and wider spread of pollutants.This vertical distribution of pollutants in the wake of the TB can be approximated using a reected Gaussian plume model for a release from an elevated source ("virtual" source height is $H TB but growing with distance from the buildings).
Observed diurnal maxima of mixed layer height (MLH) are greater in spring and summer (cf.autumn and winter), with median daytime peak in the warm seasons similar to the 75 th percentiles of the cooler seasons (Fig. 8, row 3).
A simple empirical MLH model (Table 1) with central London coefficients works surprisingly well in central Beijing, but less well in suburban Paris where values are oen over/under-predicted in winter/summer.Synoptic conditions, land cover characteristics, and human activities not accounted for by this simple parameterisation likely explain the discrepancies.Still, the results suggest the form of the model is reasonable, offering a simple, extremely fast, rst-order MLH calculation method.
Our results indicate a need for more observations at each of the scales considered in this paper.Given both the wind tunnel and water ume results are for neutral stability conditions, the inuence of TB (isolated and clusters) under non-neutral condition needs to be studied.At the microscale, large differences in surface temperatures occur. 15,16The variability of the urban form (e.g.TB height, This journal is © The Royal Society of Chemistry 2020 Faraday Discuss. number, low-rise building characteristics) and a wider range of upwind fetch conditions also need to be considered.Using the hardware model results, development of models to improve current (operational) air-quality and dispersion models should aim to enable them to capture the effects found here, thus making the models more suitable for applications in cities with TB.Future work should combine thermal inuences with TB observations (hardware, real world) for dispersion models so that there is improved capability to model vertical variability (wind, temperature, humidity) between buildings. 73he simple MLH model performance suggests further analyses of observations of both MLH data and surface energy balance partitioning across cities, with different synoptic and regional settings, would be benecial.Beyond improving our understanding of the variability of the MLH within cities, it would increase the range of data for establishing empirical coefficients.Thus, it is useful to consider these results in light of surface energy balance characteristics.Comparison of observed Bowen ratio (b) between central London 60,61 and suburban residential areas, 74,75 and similarly in Beijing, 63 indicates the latter is smaller.There is observational evidence that urban-rural contrasts in heat-ux partitioning explain spatial boundary layer depth variations 56,[76][77][78] when experiencing the same synoptic force.The timing of synoptic events, notably clear sky conditions, clearly needs to be considered to ensure the information within the empirical coefficients is appropriately disentangled.

Fig. 2
Fig. 2 Study domain with model buildings (lengths and heights in full-scale) for (a) London/Southwark and (b) Beijing/Haidian with 3D views in the horizontal plane (upper) and 2D views for vertical cross-sections along the east-west axis (lower).View from the south.Circles indicate the extent of the main model domains.Note that the model scales differ in (a) and (b).
systems allowed simultaneous velocity and scalar concentration measurements in the Beijing-Haidian hardware model (Fig. 2b, Section 2.1).A two-dimensional measurement plane is illuminated by a 100 mJ Nd:YAG pulsed laser operating at 2 Hz.Two CMOS cameras (4 MP resolution) are used for the PIV ow measurements and one 16 bit sCMOS camera (5.4 MP) for PLIF, both with a shutter speed of 1200 ms.The free-stream velocity U ref is kept at 0.45 m s À1 , assuring a constant water ow rate through the ume.No thermal effects are modelled.

Fig. 3
Fig. 3 Profiles of mean longitudinal (along-wind) velocity (U/U ref ) for the Core configuration of the London wind-tunnel model, which are ensemble averages for increasing distances from the tall building (radii r 1 to r 8 have increments of 0.5 H TB , colour-coded) plus the domain average across all sites for both Core and No Tall configurations for buildings (a) T 81 (wind direction q ¼ 0 ) and (b) T 134 (q ¼ 101 ).Logarithmic wind profile (log-law fit) using No Tall data with morphometrically 8,40 derived roughness length (z 0 ) and displacement height (z d ) is shown.Measurement locations (dots) and radii (arcs) for (c) T 81 and (d) T 134 , with the direct downwind transect in the building wakes (dashed lines) and wind direction of approaching flow (arrows) indicated.

Fig. 4
Fig. 4 Difference between Core and No Tall cases of the London model for the mean (a, c and e) longitudinal velocity (DU/U ref ) and (b, d and f) vertical velocity (DW/U ref ) as a function of straight-line distance (d/H 134 ) from the centre of building T 134 (south-east approach flow) at three heights z/H 134 (a, b) 0.27, (c, d) 0.48 and (e, f) 0.91.Lateral displacement of the measurement points (|y À y 134 |/H 134 , colour) from the centreline of T 134 (y 134 ; dashed line, Fig. 2d) indicated: along the T 134 centreline (>), off the centreline (B) and affected by the wake of building T 81 (O).

Fig. 5
Fig. 5 Mean longitudinal (along-wind) velocity (U/U ref ) in the London study area from the ADMS-Build wake model versus wind tunnel (WT) for profiles in the main wake behind (a) building T 81 (wind direction q ¼ 0 ) and (b) T 134 (q ¼ 101 ) with straight-line distance (d/H TB ) from the centre of the towers (colour).Inset maps show the (x,y) locations of the comparison sites in the main building wake; each site has data pairs at several heights.Ideal 1 : 1 (solid), AE0.1 U/U ref threshold (dotted) and linear regression fit (dashed).Percentage of the ADMS-Build data within AE0.1 U/U ref of the WT measurements are given for each distance range.
This journal is © The Royal Society of Chemistry 2020 Open Access Article.Published on 13 August 2020.Downloaded on 2/25/2021 2:34:42 PM. distance from the release point (cf.open reference).The along-wind decay is further enhanced by a stronger lateral plume spread from pollutant channelling and lateral plume branching in response to the local street topology (not measured in this experiment).

Fig. 6
Fig.6Contours of scaled tracer concentrations (log 10 (C/C s )) in the (x,z) plane in the surface layer for areas (a) without buildings, (b) with a Low configuration without the tallbuilding cluster, and (c) Tall + Low (i.e.tall building cluster with the surrounding low-rise buildings in the Beijing neighbourhood).For (b and c) the forcing wind direction is from the north (q ¼ 0 ; arrows in insets of building geometries).The heights of the maximum concentration (C max /C s ; solid lines) and of 0.5C max /C s (dashed lines) are indicated.In (c) grey lines indicate when values are less reliable, as tall buildings block the camera view.All spatial dimensions are scaled by the height of the tallest building in the cluster (H 72 ¼ 72 m).The tracer source is located at x/H 72 ¼ 0.
This journal is © The Royal Society of Chemistry 2020 Open Access Article.Published on 13 August 2020.Downloaded on 2/25/2021 2:34:42 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

Fig. 8
Fig. 8 Diurnal variations of observed (obs) and modelled mixed layer height (MLH) by season in, (row 1) Beijing and (row 2) suburban Paris, and (row 3) the observed MLH in Beijing and Paris showing median (thick line) and interquartile range (thin lines).Modelled MLH (colour) are stratified by cloud conditions (colour, Table 1).The number of days (N) with data in each season varies.Observations (black, row 1 and 2, other colours in row 3) are Vaisala CL31 ceilometer attenuated backscatter data analysed with the CABAM algorithm (Section 2.5).Seasons are winter: DJF, spring: MAM, summer: JJA, autumn: SON.Note: y-axes are different for the two rows.

Table 1
Simple mixed layer height (MLH) model equations for height and time (lower) with parameters derived for London 36 based on cloud cover state (upper).For the minimum height of the MLH (z min ) the value is the 75 th percentile, whereas the morning transition (MT) growth rate (GR MT ) and the nocturnal decay rate (DR noc ) are median values.The cloud classes are: clear nights followed by cumulus clouds during the day (clearCu); day and nightcloudless (clear), cumulus (Cu) or stratus (St); MLH is below the cloud base height (CBH) Faraday Discuss.This journal is © The Royal Society of Chemistry 2020 Open Access Article.Published on 13 August 2020.Downloaded on 2/25/2021 2:34:42 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.
Open Access Article.Published on 13 August 2020.Downloaded on 2/25/2021 2:34:42 PM.This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.m/431 m) than for Beijing (À136 m/306 m).The MBE by season indicates that peak values for suburban Paris are overestimated in winter by $100 m and underestimated in summer by $440 m.However, as the modelled MLH peak is dependent on the growth from the minimum it is important to assess this.The modelled z min (Table 1) in Beijing and in Paris are generally reasonable with the MBE and MAE below AE40 m and 100 m for Beijing, and below AE80 m and 150 m for Paris, across all seasons.Only in winter, are the Paris modelled values too high (MBE ¼ 142 m, Fig. 8).This suggests the growth rates derived from London characteristics are not large enough to capture MLH variability in Paris.Previous estimates of GR MT at SIRTA from ALC data were stratied by season (cf.cloud, Table 1).These varied from (mean AE standard deviation): 56 DJF: 149 AE 119, SON: 196 AE 137, MAM: 225 AE 138, JJA: 247 AE 142 m h À1 .These values are generally consistent with Table 1 but support the expectation of a larger peak MLH in Paris.