Impact of Megacity Shanghai on the Urban Heat-Island Effects over the Downstream City Kunshan

The impact of upstream urbanization on the enhanced urban heat-island (UHI) effects between Shanghai and Kunshan is investigated by analyzing seven years of surface observations and results from mesoscale model simulations. The observational analysis indicates that, under easterly and westerly winds, the temperature difference between Shanghai and Kunshan increases with wind speed when the wind speed \(<\) 5 m s \(^{-1}\) . The Weather Research and Forecasting (WRF) numerical model, coupled with a one-layer urban canopy model (UCM), is used to examine the UHI structure and upstream effects by replacing the urban surface of Shanghai and/or Kunshan with cropland. The WRF/UCM modelling system is capable of reproducing the surface temperature and wind field reasonably well. The simulated urban canopy wind speed is a better representation of the near-surface wind speed than is the 10-m wind speed at the centre of Shanghai. Without the urban landscape of Shanghai, the surface air temperature over downstream Kunshan would decrease by 0.2–0.4  \(^{\circ }\) C in the afternoon and 0.4–0.6  \(^{\circ }\) C in the evening. In the simulation with the urban landscape of Shanghai, a shallow cold layer is found above the UHI, with a minimum temperature of about \(-0.2\) to \(-\) 0.5  \(^{\circ }\) C during the afternoon hours. Strong horizontal divergence is found in this cold layer. The easterly breeze over Shanghai is strengthened at the surface by strong UHI effects, but weakened at upper levels. With the appearance of the urban landscape specific humidity decreases by 0.5–1 g kg \(^{-1}\) within the urban area because of the waterproof property of an urban surface. On the other hand, the upper-level specific humidity is increased because of water vapour transferred by the strong upward vertical motions.     

Transfer Coefficients of Momentum, Heat and Water Vapour in the Atmospheric Surface Layer of a Large Freshwater Lake

In studies of lake–atmosphere interactions, the fluxes of momentum, water vapour and sensible heat are often parametrized as being proportional to the differences in wind, humidity and air temperature between the water surface and a reference height above the surface. Here, the proportionality via transfer coefficients in these relationships was investigated with the eddy-covariance method at three sites within an eddy-covariance mesonet across Lake Taihu, China. The results indicate that the transfer coefficients decreased with increasing wind speed for weak winds and approached constant values for strong winds. The presence of submerged macrophytes reduced the momentum transfer (drag) coefficient significantly. At the two sites free of submerged macrophytes, the 10-m drag coefficients under neutral stability were 1.8 $(\pm \,0.4) \times \,10^{-3}$ ( ± 0.4 ) × 10 ? 3 and $1.7\,(\pm \,0.3) \times \,10^{-3 }$ 1.7 ( ± 0.3 ) × 10 ? 3 at the wind speed of $9\,\text{ m } \text{ s }^{-1}$ 9 m s ? 1 , which are 38 and 34 % greater than the prediction by the Garratt model for the marine environment.     

Streamwise Versus Spanwise Spacing of Obstacle Arrays: Parametrization of the Effects on Drag and Turbulence

A Reynolds-averaged Navier–Stokes model is used to investigate the evolution of the sectional drag coefficient and turbulent length scales with the layouts of aligned arrays of cubes. Results show that the sectional drag coefficient is determined by the non-dimensional streamwise distance (sheltering parameter), and the non-dimensional spanwise distance (channelling parameter) between obstacles. This is different than previous approaches that consider only plan area density $(\lambda _\mathrm{p})$ . On the other hand, turbulent length scales behave similarly to the staggered case (e. g. they are function of $\lambda _\mathrm{p}$ only). Analytical formulae are proposed for the length scales and for the sectional drag coefficient as a function of sheltering and channelling parameters, and implemented in a column model. This approach demonstrates good skill in the prediction of vertical profiles of the spatially-averaged horizontal wind speed.     

A New Aerodynamic Parametrization for Real Urban Surfaces

This study conducted large-eddy simulations (LES) of fully developed turbulent flow within and above explicitly resolved buildings in Tokyo and Nagoya, Japan. The more than 100 LES results, each covering a 1,000 $\times $ 1,000 m $^{2}$ area with 2-m resolution, provide a database of the horizontally-averaged turbulent statistics and surface drag corresponding to various urban morphologies. The vertical profiles of horizontally-averaged wind velocity mostly follow a logarithmic law even for districts with high-rise buildings, allowing estimates of aerodynamic parameters such as displacement height and roughness length using the von Karman constant $=$ 0.4. As an alternative derivation of the aerodynamic parameters, a regression of roughness length and variable Karman constant was also attempted, using a displacement height physically determined as the central height of drag action. Although both the regression methods worked, the former gives larger (smaller) values of displacement height (roughness length) by 20–25 % than the latter. The LES database clearly illustrates the essential difference in bulk flow properties between real urban surfaces and simplified arrays. The vertical profiles of horizontally-averaged momentum flux were influenced by the maximum building height and the standard deviation of building height, as well as conventional geometric parameters such as the average building height, frontal area index, and plane area index. On the basis of these investigations, a new aerodynamic parametrization of roughness length and displacement height in terms of the five geometric parameters described above was empirically proposed. The new parametrizations work well for both real urban morphologies and simplified model geometries.     

Surface Heterogeneity Effects on Regional-Scale Fluxes in the Stable Boundary Layer: Aerodynamic Roughness Length Transitions

The effects of abrupt streamwise transitions of the aerodynamic roughness length ( $z_\mathrm{o}$ z o ) on the stable atmospheric boundary layer are evaluated using a series of large-eddy simulations based on the first Global Energy and Water Cycle Experiment Atmospheric Boundary Layer intercomparison study (GABLS1). Four $z_\mathrm{o}$ z o values spanning three orders of magnitude are used to create all possible binary distributions with each arranged into patches of characteristic length scales equal to roughly one-half, one, and two times the equivalent homogeneous boundary-layer height. The impact of the heterogeneity on mean profiles of wind speed and temperature, on surface fluxes of heat and momentum, and on internal boundary-layer dynamics are considered. It is found that $z_\mathrm{o}$ z o transitions do not significantly alter the functional relationship between the average surface fluxes and the mean profiles of wind speed and potential temperature. Although this suggests that bulk similarity theory is applicable for modelling the stable boundary layer over $z_\mathrm{o}$ z o heterogeneity, effective surface parameters must still be specified. Existing models that solve for effective roughness lengths of momentum and heat are evaluated and compared to values derived from the simulation data. The existing models are unable to accurately reproduce both the values of the effective aerodynamic roughness lengths and their trends as functions of patch length scale and stability. A new model for the effective aerodynamic roughness length is developed to exploit the benefits of the other models tested. It accurately accounts for the effects of the heterogeneity and stratification on the blending height and effective aerodynamic roughness length. The new model provides improved average surface fluxes when used with bulk similarity.     

Large-Eddy Simulation of Mesoscale Circulations Forced by Inhomogeneous Urban Heat Island

The large-eddy simulation mode of the Weather Research and Forecasting model is employed to simulate the planetary boundary-layer characteristics and mesoscale circulations forced by an ideal urban heat island (UHI). In our simulations, the horizontal heterogeneity of the UHI intensity distribution in urban areas is considered and idealized as a cosine function. Results indicate that the UHI heating rate and the UHI intensity heterogeneity affect directly the spatial distribution of the wind field; a stronger UHI intensity produces a maximum horizontal wind speed closer to the urban centre. The strong advection of warm air from the urban area to the rural area in the upper part of the planetary boundary-layer causes a more stable atmospheric stratification over both the urban and rural areas. The mesoscale sensible heat flux caused by the UHI circulation increases with UHI intensity but vanishes when the background wind speed is sufficiently high $(>$ 3.0  $\mathrm{{m\,s}}^{-1})$ .     

Numerical Simulations of Laminar Air–Water Flow of a Non-linear Progressive Wave at Low Wind Speed

A numerical simulation for two-dimensional laminar air–water flow of a non-linear progressive water wave with large steepness is performed when the background wind speed varies from zero to the wave phase speed. It is revealed that in the water the difference between the analytical solution of potential flow and numerical solution of viscous flow is very small, indicating that both solutions of the potential flow and viscous flow describe the water wave very accurately. In the air the solutions of potential and viscous flows are very different due to the effects of viscosity. The velocity distribution in the airflow is strongly influenced by the background wind speed and it is found that three wind speeds, $U=0$ , $U=u_m$ (the maximum orbital velocity of a water wave), and $U=c$ (the wave phase speed), are important in distinguishing different features of the flow patterns.     

Determination Of The Surface Drag Coefficient

This study examines the dependence of the surface drag coefficienton stability, wind speed, mesoscale modulation of the turbulent flux and method of calculation of the drag coefficient. Data sets over grassland, sparse grass, heather and two forest sites are analyzed. For significantly unstable conditions, the drag coefficient does not depend systematically on z/L but decreases with wind speed for fixed intervals of z/L, where L is the Obukhov length. Even though the drag coefficient for weak wind conditions is sensitive to the exact method of calculation and choice of averaging time, the decrease of the drag coefficient with wind speed occurs for all of the calculation methods. A classification of flux calculation methods is constructed, which unifies the most common previous approaches.The roughness length corresponding to the usual Monin–Obukhovstability functions decreases with increasing wind speed. This dependence on wind speed cannot be eliminated by adjusting the stability functions. If physical, the decrease of the roughness length with increasing wind speed might be due to the decreasing role of viscous effectsand streamlining of the vegetation, although these effects cannot be isolated from existing atmospheric data.For weak winds, both the mean flow and the stress vector often meander significantly in response to mesoscale motions. The relationship between meandering of the stress and wind vectors is examined. For weak winds, the drag coefficient can be sensitive to the method of calculation, partly due to meandering of the stress vector.     

Storm Surges in New York During Hurricane Sandy in 2012: A Verification of the Wind-Stress Tide Relation

In October 2012 Hurricane Sandy devastated New York City and its vicinity caused mainly by the storm surge, which is the water height above normal astronomical tide level. The meteorological conditions were as follows: minimum central pressure, 962 hPa, highest sustained wind speed 27.1 m s $^{-1}$ ? 1 and maximum gust 37.8 m s $^{-1}$ ? 1 . The peak storm surge was at 3.9 m and the peak storm tide at 4.4 m (which is referenced above mean lower low water). The wind-stress tide relation shows that $S=K\,V^{2}$ S = K V 2 , where $S$ S is the storm surge, $V$ V is the wind speed and $K$ K is the coefficient. It is found that with $S$ S in units of m, and $V$ V in  m s $^{-1}$ ? 1 , $K = 0.0051$ K = 0.0051 with $R^{2}= 0.91$ R 2 = 0.91 ( $R$ R is the correlation coefficient) indicating that 91 % of the total variation of the storm surge can be explained by variations in the wind stress, which is proportional to $V^{2}$ V 2 . Similar results were obtained during Hurricane Irene in 2011, which also affected the New York area. Therefore, this simple wind stress-tide relation should be useful in coastal engineering, urban planning, and emergency management.     

Equilibrium Atmospheric Boundary-Layer Flow: Computational Fluid Dynamics Simulation with Balanced Forces

Forcing relationships in steady, neutrally stratified atmospheric boundary-layer (ABL) flow are thoroughly analyzed. The ABL flow can be viewed as balanced between a forcing and a drag term. The drag term results from turbulent stress divergence, and above the ABL, both the drag and the forcing terms vanish. In computational wind engineering applications, the ABL flow is simulated not by directly specifying a forcing term in the ABL but by specifying boundary conditions for the simulation domain. Usually, these include the inflow boundary and the top boundary conditions. This ‘boundary-driven’ ABL flow is dynamically different from its real counterpart, and this is the major reason that the simulated boundary-driven ABL flow does not maintain horizontal homogeneity. Here, first a dynamical approach is proposed to develop a neutrally stratified equilibrium ABL flow. Computational fluid dynamics (CFD) software (Fluent 6.3) with the standard \(k\) – \(\varepsilon \) turbulence model is employed, and by applying a driving force profile, steady equilibrium ABL flows are simulated by the model. Profiles of wind speed and turbulent kinetic energy (TKE) derived using this approach are reasonable in comparison with the conventional logarithmic law and with observational data respectively. Secondly, the equilibrium ABL profiles apply as inflow conditions to simulate the boundary-driven ABL flow. Simulated properties between the inlet and the outlet sections across a fetch of 10 km are compared. Although profiles of wind speed, TKE, and its dissipation rate are consistently satisfactory under higher wind conditions, a deviation of TKE and its dissipation rate between the inlet and outlet are apparent (7–8 %) under lower wind-speed conditions (2 m s \(^{-1}\) at 10 m). Furthermore, the simulated surface stress systematically decreases in the downwind direction. A redistribution of the pressure field is also found in the simulation domain, which provides a different driving pattern from the realistic case in the ABL.     

The Effect of Wind-Turbine Wakes on Summertime US Midwest Atmospheric Wind Profiles as Observed with Ground-Based Doppler Lidar

We examine the influence of a modern multi-megawatt wind turbine on wind and turbulence profiles three rotor diameters ( $D$ D ) downwind of the turbine. Light detection and ranging (lidar) wind-profile observations were collected during summer 2011 in an operating wind farm in central Iowa at 20-m vertical intervals from 40 to 220 m above the surface. After a calibration period during which two lidars were operated next to each other, one lidar was located approximately $2D$ 2 D directly south of a wind turbine; the other lidar was moved approximately $3D$ 3 D north of the same wind turbine. Data from the two lidars during southerly flow conditions enabled the simultaneous capture of inflow and wake conditions. The inflow wind and turbulence profiles exhibit strong variability with atmospheric stability: daytime profiles are well-mixed with little shear and strong turbulence, while nighttime profiles exhibit minimal turbulence and considerable shear across the rotor disk region and above. Consistent with the observations available from other studies and with wind-tunnel and large-eddy simulation studies, measurable reductions in wake wind-speeds occur at heights spanning the wind turbine rotor (43–117 m), and turbulent quantities increase in the wake. In generalizing these results as a function of inflow wind speed, we find the wind-speed deficit in the wake is largest at hub height or just above, and the maximum deficit occurs when wind speeds are below the rated speed for the turbine. Similarly, the maximum enhancement of turbulence kinetic energy and turbulence intensity occurs at hub height, although observations at the top of the rotor disk do not allow assessment of turbulence in that region. The wind shear below turbine hub height (quantified here with the power-law coefficient) is found to be a useful parameter to identify whether a downwind lidar observes turbine wake or free-flow conditions. These field observations provide data for validating turbine-wake models and wind-tunnel observations, and for guiding assessments of the impacts of wakes on surface turbulent fluxes or surface temperatures downwind of turbines.     

Intercomparison of Methods for the Simultaneous Estimation of Zero-Plane Displacement and Aerodynamic Roughness Length from Single-Level Eddy-Covariance Data

We applied three approaches to estimate the zero-plane displacement $d$ through the aerodynamic measurement height $z$ (with $z = z_{m}- d$ and $z_{m}$ being the measurement height above the surface), and the aerodynamic roughness length $z_{0}$ , from single-level eddy covariance data. Two approaches (one iterative and one regression-based) were based on the universal function in the logarithmic wind profile and yielded an inherently simultaneous estimation of both $d$ and $z_{0}$ . The third approach was based on flux–variance similarity, where estimation of $d$ and consecutive estimation of $z_{0}$ are independent steps. Each approach was further divided into two methods differing either with respect to the solution technique (profile approaches) or with respect to the variable (variance of vertical wind and temperature, respectively). All methods were applied to measurements above a large, growing wheat field where a uniform canopy height and its frequent monitoring provided plausibility limits for the resulting estimates of time-variant $d$ and $z_{0}$ . After applying, for each approach, a specific data filtering that accounted for the range of conditions (e.g. stability) for which it is valid, five of the six methods were able to describe the temporal changes of roughness parameters associated with crop growth and harvest, and four of them agreed on $d$ to within 0.3 m most of the time. Application of the same methods to measurements with a more heterogeneous footprint consisting of fully-grown sugarbeet and a varying contribution of adjacent harvested fields exhibited a plausible dependence of the roughness parameters on the sugarbeet fraction. It also revealed that the methods producing the largest outliers can differ between site conditions and stability. We therefore conclude that when determining $d$ for canopies with unknown properties from single-level measurements, as is increasingly done, it is important to compare the results of a number of methods rather than rely on a single one. An ensemble average or median of the results, possibly after elimination of methods that produce outliers, can help to yield more robust estimates. The estimates of $z_{0}$ were almost exclusively physically plausible, although $d$ was considered unknown and estimated simultaneously with the methods and results described above.     

Thermal Stratification Effects on Flow Over a Generic Urban Canopy

The influence of local surface heating and cooling on flow over urban-like roughness is investigated using large-eddy simulations. By adjusting the incoming or outgoing heat flux from the ground surface, various degrees of local thermal stratification, represented by a Richardson number \((Ri_\tau )\) , were attained. Drag and heat transfer coefficients, turbulence structure, integral length scales, and the strength of quadrant events that contribute to momentum and heat fluxes were obtained and are compared with locally stable, neutral and unstable flows. With increasing \(Ri_\tau \) , or equivalently as the flow characteristics change from local thermal instability to stability, a gradual decline in the drag and heat transfer coefficients is observed. These values are found to be fairly independent of the type of thermal boundary condition (constant heat flux or constant temperature) and domain size. The maps of anisotropy invariants showed that for the values of \(Ri_\tau \) considered, turbulence structures are almost the same in shape for neutral and unstable cases but differ slightly from those in the stable case. The degree of anisotropy is found to decrease as \(Ri_\tau \) increases from \(-2\) to 2.5. Compared to the neutral case, the integral length scales are shortened in the streamwise and vertical direction by ground cooling, but enhanced in the vertical direction with ground heating. Quadrant analysis showed that an increase in floor heating increases the strength of ejections above the canopy. However, the contributions of updrafts or downdrafts to the heat flux are found not to be significantly influenced by the type of local thermal stratification for the values of \(Ri_\tau \) considered. From the octant analysis, the transport mechanisms of momentum and heat above the canopy are found to be very similar in both locally unstable and stable flows.     

A Scale-Adaptive Approach for Spatially-Varying Urban Morphology Characterization in Boundary Layer Parametrization Using Multi-Resolution Analysis

Urban morphology characterization is crucial for the parametrization of boundary-layer development over urban areas. One complexity in such a characterization is the three-dimensional variation of the urban canopies and textures, which are customarily reduced to and represented by one-dimensional varying parametrization such as the aerodynamic roughness length $z_{0}$ and zero-plane displacement $d$ . The scope of the paper is to provide novel means for a scale-adaptive spatially-varying parametrization of the boundary layer by addressing this 3-D variation. Specifically, the 3-D variation of urban geometries often poses questions in the multi-scale modelling of air pollution dispersion and other climate or weather-related modelling applications that have not been addressed yet, such as: (a) how we represent urban attributes (parameters) appropriately for the multi-scale nature and multi-resolution basis of weather numerical models, (b) how we quantify the uniqueness of an urban database in the context of modelling urban effects in large-scale weather numerical models, and (c) how we derive the impact and influence of a particular building in pre-specified sub-domain areas of the urban database. We illustrate how multi-resolution analysis (MRA) addresses and answers the afore-mentioned questions by taking as an example the Central Business District of Oklahoma City. The selection of MRA is motivated by its capacity for multi-scale sampling; in the MRA the “urban” signal depicting a city is decomposed into an approximation, a representation at a higher scale, and a detail, the part removed at lower scales to yield the approximation. Different levels of approximations were deduced for the building height $\bar{{H}}$ and planar packing density $\lambda _\mathrm{p}$ . A spatially-varying characterization with a scale-adaptive capacity is obtained for the boundary-layer parameters (aerodynamic roughness length $z_{0}$ and zero-plane displacement $d$ ) using the MRA-deduced results for the building height and the planar packing density with a morphometric model; an attribute that is shown to be of great advantage to multi-scale and multi-resolution numerical weather prediction models.     

Impacts of Mixing Processes in Nocturnal Atmospheric Boundary Layer on Urban Ozone Concentrations

A number of open questions remain regarding the role of low-level jets (LLJs) and nocturnal mixing processes in the buildup of tropospheric ozone. The prevalence of southerly winds and LLJs in the U.S. Southern Great Plains during summer makes this region an ideal site for investigating the structure of the nocturnal boundary layer and its impacts on urban air quality. Ozone $(\mathrm{O}_{3})$ and nitrogen oxide concentrations measured at regulatory monitoring sites in the Oklahoma City (OKC) area and simulations with the Weather Research and Forecasting with Chemistry (WRF/Chem) model were analyzed to show how the nocturnal LLJ moderates boundary-layer mixing processes and air quality. Datasets collected during the Joint Urban 2003 campaign, which took place in July 2003 in OKC, provided detailed information about nocturnal boundary-layer structure and dynamics. In general, ${\mathrm{O}_{3}}$ time series show the expected behavior that urban ${\mathrm{O}_{3}}$ concentrations decrease at night due to nitrogen oxide titration reactions, but elevated ${\mathrm{O}_{3}}$ concentrations and secondary ${\mathrm{O}_{3}}$ peaks are also seen quite frequently after sunset. LLJs developed on most nights during the study period and were associated with strong vertical wind shear, which affected the boundary-layer stability and structure. Near-surface ${\mathrm{O}_{3}}$ concentrations are higher during less stable nights when active mixing persists throughout the night. The WRF/Chem model results agree well with the observations and further demonstrate the role of LLJs in moderating nocturnal mixing processes and air quality. The highest nocturnal ${\mathrm{O}_{3}}$ concentrations are linked to a strong LLJ that promotes both nocturnal long-range transport and persistent downward mixing of ${\mathrm{O}_{3}}$ from the residual layer to the surface.     

Similarity Scaling Over a Steep Alpine Slope

In this study, we investigate the validity of similarity scaling over a steep mountain slope (30–41 $^\circ $ ). The results are based on eddy-covariance data collected during the Slope Experiment near La Fouly (SELF-2010); a field campaign conducted in a narrow valley of the Swiss Alps during summer 2010. The turbulent fluxes of heat and momentum are found to vary significantly with height in the first few metres above the inclined surface. These variations exceed by an order of magnitude the well-accepted maximum 10 % required for the applicability of Monin–Obukhov similarity theory in the surface layer. This could be due to a surface layer that is too thin to be detected or to the presence of advective fluxes. It is shown that local scaling can be a useful tool in these cases when surface-layer theory breaks down. Under convective conditions and after removing the effects of self-correlation, the normalized standard deviations of slope-normal wind velocity, temperature and humidity scale relatively well with $z/\varLambda $ , where $z$ is the measurement height and $\varLambda (z)$ the local Obukhov length. However, the horizontal velocity fluctuations are not correlated with $z/\varLambda $ under all stability regimes. The non-dimensional gradients of wind velocity and temperature are also investigated. For those, the local scaling appears inappropriate, particularly at night when shallow drainage flows prevail and lead to negative wind-speed gradients close to the surface.     

Large-Eddy Simulation of Heat Transfer from a Single Cube Mounted on a Very Rough Wall

The local thermal effects in the wake of a single cube with a strong heated rear face, representing a large building in an urban area, are studied using large-eddy simulations (LES) for various degrees of heating, which are characterized by the local Richardson number, $Ri$ . New wall models are implemented for momentum and temperature and comparison of the flow and thermal fields with the wind-tunnel data of Richards et al. (J Wind Eng Ind Aerodyn 94, 621–636, 2006) shows fair agreement. Buoyancy effects are quite evident at low $Ri$ and a significant increase in the turbulence levels is observed for such flows. Apart from the comparisons with experiments, further analysis included the estimation of the thermal boundary-layer thickness and heat transfer coefficient for all $Ri$ . For sufficiently strong heating, the heat transfer coefficient at the leeward face is found to be higher than the roof surface. This suggests that, beyond a certain $Ri$ value, buoyancy forces from the former surface dominate the strong streamwise convection of the latter. Quadrant analysis along the shear layer behind the cube showed that the strength of sweeps that contribute to momentum flux is considerably enhanced by heating. The contribution of different quadrants to the heat flux is found to be very different to that of the momentum flux for lower $Ri$ .     

On the Mechanism of Air Pollutant Removal in Two-Dimensional Idealized Street Canyons: A Large-Eddy Simulation Approach

Flow resistance, ventilation, and pollutant removal for idealized two-dimensional (2D) street canyons of different building-height to street-width (aspect) ratios $AR$ are examined using the friction factor $f$ , air exchange rate (ACH), and pollutant exchange rate (PCH), respectively, calculated by large-eddy simulation (LES). The flows are basically classified into three characteristic regimes, namely isolated roughness, wake interference, and skimming flow, as functions of the aspect ratios. The LES results are validated by various experimental and numerical datasets available in the literature. The friction factor increases with decreasing aspect ratio and reaches a peak at $AR = 0.1$ in the isolated roughness regime and decreases thereafter. As with the friction factor, the ACH increases with decreasing aspect ratio in the wake interference and skimming flow regimes, signifying the improved aged air removal for a wider street canyon. The PCH exhibits a behaviour different from its ACH counterpart in the range of aspect ratios tested. Pollutants are most effectively removed from the street canyon with $AR = 0.5$ . However, a minimum of PCH is found nearby at $AR = 0.3$ , at which the pollutant removal is sharply weakened. Besides, the ACH and PCH are partitioned into the mean and turbulent components to compare their relative contributions. In line with our earlier Reynolds-averaged Navier–Stokes calculations (Liu et al., Atmos Environ 45:4763–4769, 2011), the current LES shows that the turbulent components contribute more to both ACH and PCH, consistently demonstrating the importance of atmospheric turbulence in the ventilation and pollutant removal for urban areas.     

Observations of the Temperature and Humidity Structure Parameter Over Heterogeneous Terrain by Airborne Measurements During the LITFASS-2003 Campaign

We used numerical simulations to investigate the general relationship between urban morphology and the intensity of wind gusts in built-up areas at the pedestrian level. The simulated urban boundary layer developed over a 19.2 km (length) \(\times \) 4.8 km (width) \(\times \) 1.0 km (height) simulation domain, with 2-m resolution in all directions, to explicitly resolve the detailed shapes of buildings and the flow at the pedestrian level. This complex computation was accomplished using the lattice Boltzmann method and by implementing a large-eddy simulation model. To generalize the results, a new parameter that expresses the intensity of gusts (the gust index, \({\tilde{U}}_{ max})\) was defined as the local maximum wind speed divided by the freestream velocity. In addition, this parameter was decomposed into the mean wind-speed ratio, \({\tilde{U}} \) and turbulent gust ratio, \({\tilde{U}}^{{\prime }}\) to evaluate the qualities of gusts. These parameters were useful for quantitatively comparing the gust intensities within urban canopies at different locations or even among different experiments. In addition, the entire horizontal domain was subdivided into homogeneous square patches, in which both the simulated gust parameters and the morphological characteristics of building geometries were averaged. This procedure masked the detailed structure of individual buildings but retained the bulk characteristics of the urban morphology. At the pedestrian level, the gust index decreased with increasing building cover. Compared to \({\tilde{U}} \), the quantity \({\tilde{U}}^{{\prime }}\) notably contributed to the index throughout the range of plan area index \((\lambda _p)\) values. The dependences of all normalized wind-speed ratios transiently changed at \(\lambda _p =~0.28\). In cases where \(\lambda _p < 0.28, {\tilde{U}} \) decreased with increasing \(\lambda _p \), although \({\tilde{U}}^{{\prime }}\) was almost constant. In cases where \(\lambda _p > 0.28, {\tilde{U}}\) was almost constant and \({\tilde{U}}^{{\prime }}\) decreased with increasing \(\lambda _p \). This was explained by the change in flow regimes within the building canyon. At a higher elevation above the canopy layer, \(\lambda _p \) becomes less relevant to normalized wind-speed ratios, and instead the aerodynamic roughness length became important.