Characteristics of vertical turbulent velocities in the urban convective boundary layer

Turbulence measurements of the vertical velocity component were obtained by an instrumented aircraft under fair weather conditions in the St. Louis, Missouri, metropolitan area. Time series of vertical velocity fluctuations from horizontal flight segments made in the lower part of and near the middle of the convective boundary layer (CBL) over the urban area and surrounding region were subjected to various statistical and objective analyses. Higher order vertical velocity moments, and positive and negative velocity statistics, were computed. The horizontal dimensions of updrafts and downdrafts, and related properties of these turbulent eddies were derived by conditional sampling analysis. Emphasis is on a comparison of the results from urban and selected rural measurements from the lower part of the CBL.The vertical velocity probability density distribution for each segment was positively skewed and the mode was negative. The means and standard deviations of positive and negative velocity fluctuations were greater over the urban area. The urban vertical velocity variance was 50% greater than rural values, and power spectra revealed greater production of vertical turbulent energy in the urban area over a wide frequency range.The mean and maximum widths of downdrafts were generally larger than the corresponding values for updrafts. Differences between urban and rural eddy sizes were not statistically significant. The widths of the largest updraft and downdraft are comparable to the boundary-layer depth, Z _i, and the mean value of the ratio of spectral peak wavelength to Z _iwas about 1.3 and 1.1 for urban and rural areas, respectively. Convective similarity scaling parameters appeared to order both the urban and rural measurements.On assignment from the National Oceanic and Atmospheric Administration, U.S. Dept. of Commerce.     

Relating Urban Surface-layer Structure to Upwind Terrain for the Salford Experiment (Salfex)

Profiles of wind and turbulence over an urban area evolve with fetch in response to surface characteristics. Sodar measurements, taken on 22 April 2002 during the Salford Experiment in the UK (Salfex), are here related to upstream terrain. A logarithmic layer up to z = 65m was observed in all half-hour averaged profiles. Above this height the profile showed a different vertical gradient, suggesting a change in surface cover upstream. The drag coefficient varied by a factor of two over only a 20° direction change. Turbulence intensity (σ _x ) for each wind component (x) decreased with height, but the ratio suggested an underestimate of σ _u compared to previous results. Mean urban and suburban cover fraction within the source area for each height decreased sharply between z = 20 and 50m, increasing slightly above. The near-convergence of cover fractions thus occured for source areas of minimum length ≈ 2,200 m. In comparison, the mean length scale of heterogeneity L _P was calculated from surface cover data to be 1,284 m, and the corresponding mean blending height h _b was 175 m. Finally, the mean streamline angle, α, was negative and the magnitude decreased with height. An exponential fit to α for z ≤ 65m gave an e-folding height scale of 159 m. A simple relationship between this height scale and L _P was assumed, giving L _P ≈ 1,080 m, which is in reasonable agreement with the estimate from surface cover type. The results suggest that more emphasis is required on modelling and measuring surface-layer flow over heterogeneous urban canopies.     

Measurements Of Thermal Updraft Intensity Over Complex Terrain Using American White Pelicans And A Simple Boundary-Layer Forecast Model

An examination of boundary-layer meteorological and avian aerodynamic theories suggests that soaring birds can be used to measure the magnitude of vertical air motions within the boundary layer. These theories are applied to obtain mixed-layer normalized thermal updraft intensity over both flat and complex terrain from the climb rates of soaring American white pelicans and from diagnostic boundary-layer model-produced estimates of the boundary-layer depth z_i and the convective velocity scale w_*. Comparison of the flatland data with the profiles of normalized updraft velocity obtained from previous studies reveals that the pelican-derived measurements of thermal updraft intensity are in close agreement with those obtained using traditional research aircraft and large eddy simulation (LES) in the height range of 0.2 to 0.8 z_i. Given the success of this method, the profiles of thermal vertical velocity over the flatland and the nearby mountains are compared. This comparison shows that these profiles are statistically indistinguishable over this height range, indicating that the profile for thermal updraft intensity varies little over this sample of complex terrain. These observations support the findings of a recent LES study that explored the turbulent structure of the boundary layer using a range of terrain specifications. For terrain similar in scale to that encountered in this study, results of the LES suggest that the terrain caused less than an 11% variation in the standard deviation of vertical velocity.     

Influence of the boundary layer height on the global air–sea surface fluxes

Results from large-eddy simulations and field measurements have previously shown that the velocity field is influenced by the boundary layer height, z _i , during close to neutral, slightly unstable, atmospheric stratification. During such conditions the non-dimensional wind profile, φ _m , has been found to be a function of both z/L and z _i /L. At constant z/L, φ _m decreases with decreasing boundary layer height. Since φ _m is directly related to the parameterizations of the air–sea surface fluxes, these results will have an influence when calculating the surface fluxes in weather and climate models. The global impact of this was estimated using re-analysis data from 1979 to 2001 and bulk parameterizations. The results show that the sum of the global latent and sensible mean heat fluxes increase by 0.77 W m⁻² or about 1% and the mean surface stress increase by 1.4 mN m⁻² or 1.8% when including the effects of the boundary layer height in the parameterizations. However, some regions show a larger response. The greatest impact is found over the tropical oceans between 30°S and 30°N. In this region the boundary layer height influences the non-dimensional wind profile during extended periods of time. In the mid Indian Ocean this results in an increase of the mean annual heat fluxes by 2.0 W m⁻² and an increase of the mean annual surface stress by 2.6 mN m⁻².     

The Effects of Inner- and Outer-Layer Turbulence in a Convective Boundary Layer on the Near-Neutral Inertial Sublayer Over an Urban-Like Surface

A large-eddy simulation of the atmospheric boundary layer, large enough to contain both an urban surface layer and a convective mixed layer, was performed to investigate inner-layer and outer-layer scale motions. The objective was to determine the applicability of Monin–Obukhov similarity theory to inner-layer motions, to investigate the influence of outer-layer motions on surface-layer structure, as well as to assess the interaction of the two scales of motion. The urban surface roughness consisted of square-patterned cubic buildings of dimension H (40 m). A spatial filter was used to decompose the two scales in the inertial sublayer. The horizontal square filter of size 10H was effective in separating the inner-layer (surface-layer height ≈ 2 H) and outer-layer scales (boundary-layer height δ ≈ 30H), where the Reynolds stress contribution of the inner layer dominates in the logarithmic layer (depth 2H). Similarity coefficients for velocity fluctuations were successfully determined for inner-layer motions in the surface layer, proving the robustness of Monin–Obukhov similarity for surface-layer turbulence. The inner-layer structures exhibit streaky structures that have similar streamwise length but narrower spanwise width relative to the streamwise velocity fluctuation field, consistent with observations from an outdoor scale model. The outer-layer motions to some extent influence the location of ejections and sweeps through updraft and downdraft motions, respectively, thus, disturbing the homogeneity and similarity of inner-layer motions. Although the horizontal averages of the variances and covariance of motions reveal that the Reynolds stresses are dominated by inner-layer structures, the localized influence of the interaction of outer-layer horizontal and inner-layer vertical motions on the Reynolds stress is not insignificant.     

Velocity and Surface Shear Stress Distributions Behind a Rough-to-Smooth Surface Transition: A Simple New Model

A simple new model is proposed to predict the distribution of wind velocity and surface shear stress downwind of a rough-to-smooth surface transition. The wind velocity is estimated as a weighted average between two limiting logarithmic profiles: the first log law, which is recovered above the internal boundary-layer height, corresponds to the upwind velocity profile; the second log law is adjusted to the downwind aerodynamic roughness and local surface shear stress, and it is recovered near the surface, in the equilibrium sublayer. The proposed non-linear form of the weighting factor is equal to ln(z/z ₀₁)/ln(δ _i /z ₀₁), where z, δ _i and z ₀₁ are the elevation of the prediction location, the internal boundary-layer height at that downwind distance, and the upwind surface roughness, respectively. Unlike other simple analytical models, the new model does not rely on the assumption of a constant or linear distribution for the turbulent shear stress within the internal boundary layer. The performance of the new model is tested with wind-tunnel measurements and also with the field data of Bradley. Compared with other existing analytical models, the proposed model shows improved predictions of both surface shear stress and velocity distributions at different positions downwind of the transition.     

Formulation of the thermal internal boundary layer in a mesoscale model

Mesoscale models using a non-local K-scheme for parameterization of boundary-layer processes require an estimate of the planetary boundary layer (PBL) height z _i at all times. In this paper, two-dimensional sea-breeze experiments are carried out to evaluate three different formulations for the advective contribution in the z _i prognostic equation of Deardorff (1974).Poor representation of the thermal internal boundary layer in the sea breeze is obtained when z _i is advected by the wind at level z _i . However, significantly better results are produced if the mean PBL wind is used for the advecting velocity, or if z _i is determined simply by checking for the first sufficiently stable layer above the ground.A Lagrangian particle model is used to demonstrate the effect of each formulation on plume dispersion by the sea breeze.     

Concentration and gas/particle partitioning of polychlorinated biphenyls (PCBs) at an industrial site at Bursa, Turkey

Gas and particle phase concentrations of atmospheric polychlorinated biphenyls (PCBs) were measured at an urban/industrial site in the city of Bursa, Turkey. PCB concentration levels were presented between July 2004 and May 2005. Average particle and gas phase concentrations of individual PCB congeners ranged from 0.08 (PCB-183) to 6.86 (PCB-49) pg m^− 3 and from 0.01 (PCB-209) to 47.2 (PCB-33) pg m^− 3, respectively. The mean concentration of total (gas + particle) PCBs varied between 24.27 and 666.21 pg m^− 3 with an average of 287.27 ± 174.80 pg m^− 3. PCB concentrations at the sampling site were higher than the concentrations reported at non-urban sites. PCBs partitioned between gas and particle phases and the partitioning was examined according to different approaches such as logK_p–logP_L^o, logK_p–logK_OA and the Junge–Pankow model. In order to present possible interactions, a correlation matrix based on PCB congeners and meteorological parameters was constructed. Application of the Clausius–Clapeyron equation yielded a low slope value indicating possible emissions from local and regional sources originating mainly from urban/industrial areas, landfill and waste incineration plant. Then, likely dry deposition fluxes were estimated depending on reported dry deposition velocity and atmospheric concentration values.     

A numerical study of daily transitions in the convective boundary layer

We examine daily (morning–afternoon) transitions in the atmospheric boundary layer based on large-eddy simulations. Under consideration are the effects of the stratification at the top of the mixed layer and of the wind shear. The results describe the transitory behaviour of temperature and wind velocity, their second moments, the boundary-layer height Z _m (defined by the maximum of the potential temperature gradient) and its standard deviation σ _m , the mixed-layer height z _i (defined by the minimum of the potential temperature flux), entrainment velocity W _e, and the entrainment flux H _i . The entrainment flux and the entrainment velocity are found to lag slightly in time with respect to the surface temperature flux. The simulations imply that the atmospheric values of velocity variances, measured at various instants during the daytime, and normalized in terms of the actual convective scale w_*, are not expected to collapse to a single curve, but to produce a significant scatter of observational points. The measured values of the temperature variance, normalized in terms of the actual convective scale Θ_*, are expected to form a single curve in the mixed layer, and to exhibit a considerable scatter in the interfacial layer.     

The characteristics of turbulent velocity components in the surface layer under convective conditions

It is proposed that the ratios of the standard deviations of the horizontal velocity components to the friction velocity in the surface layer under convective conditions depend only onz _i /L wherez _i is the height of the lowest inversion andL is the Monin-Obukhov length. This hypothesis is tested by using observations from several data sets over uniform surfaces and appears to fit the data well. Empirical curves are fitted to the observations which have the property that at largez _i /-L, the standard deviations become proportional tow _*, the convective scaling velocity.Fluctuations of vertical velocity obtained from the same experiments scale withz/L, wherez is the height above the surface, in good agreement with Monin-Obukhov theory.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Performance of various similarity functions for nondimensional wind and temperature profiles in the surface layer in stable conditions

The linear functions for non-dimensional wind and temperature profiles are commonly used to describe the surface layer fluxes in atmospheric models. However, their applicability is limited to smaller values of the stability parameter z/L (where z is the height above ground and L is the Obukhov length) i.e. z/L < 1.0. These linear functions have been modified (Webb 1970, Quart. J. Roy. Meteor. Soc. 96, 67–90; Clarke 1970, Quart. J. Roy. Meteor. Soc. 96, 91–114; Hicks 1976, Quart. J. Roy. Meteor. Soc. 102, 535–551; Beljaars and Holtslag 1991, J. Appl. Meteorol. 30, 327–341; Cheng and Brutsaert 2005, Boundary-Layer Meteorol. 114, 519–538) over the years for calculating fluxes when z/L > 1.0 under strongly stable conditions. In view of this, the objective of the present study is to analyze the performance of these similarity functions to compute surface fluxes in stable conditions.The meteorological observations from the Cooperative Atmosphere-Surface Exchange Study (CASES-99) experiment are utilized for computing the surface fluxes in stable conditions. The computed fluxes are found to be reasonably close to those observed. The ratio of observed to computed fluxes reveals that the computed fluxes are close to the observations for all the similarity functions for z/L < 1.0 whereas the computed values show relatively a large scatter from observations for z/L > 1.0. The computed values of u and heat flux do not show significant differences from those observed at 99% confidence limit. The performance of all the similarity functions considered here is found to be comparable to each other in strongly stable conditions.     

Flux Footprints in the Convective Boundary Layer: Large-Eddy Simulation and Lagrangian Stochastic Modelling

We investigated the flux footprints of receptors at different heights in the convective boundary layer (CBL). The footprints were derived using a forward Lagrangian stochastic (LS) method coupled with the turbulent fields from a large-eddy simulation model. Crosswind-integrated flux footprints shown as a function of upstream distances and sensor heights in the CBL were derived and compared using two LS particle simulation methods: an instantaneous area release and a crosswind linear continuous release. We found that for almost all sensor heights in the CBL, a major positive flux footprint zone was located close to the sensor upstream, while a weak negative footprint zone was located further upstream, with the transition band in non-dimensional upwind distances −X between approximately 1.5 and 2.0. Two-dimensional (2D) flux footprints for a point sensor were also simulated. For a sensor height of 0.158 z _i, where z _i is the CBL depth, we found that a major positive flux footprint zone followed a weak negative zone in the upstream direction. Two even weaker positive zones were also present on either side of the footprint axis, where the latter was rotated slightly from the geostrophic wind direction. Using CBL scaling, the 2D footprint result was normalized to show the source areas and was applied to real parameters obtained using aircraft-based measurements. With a mean wind speed in the CBL of U = 5.1 m s⁻¹, convective velocity of w _* = 1.37 m s⁻¹, CBL depth of z _i = 1,000 m, and flight track height of 159 m above the surface, the total flux footprint contribution zone was estimated to range from about 0.1 to 4.5 km upstream, in the case where the wind was perpendicular to the flight track. When the wind was parallel to the flight track, the total footprint contribution zone covered approximately 0.5 km on one side and 0.8 km on the other side of the flight track.     

Strong Updraft at a Sea-Breeze Front and Associated Vertical Transport of Near-Surface Dense Aerosol Observed by Doppler Lidar and Ceilometer

To study the wind field within the atmospheric boundary layer over the Tokyo metropolitan area, Doppler lidar observations were made 45 km north of Sagami Bay and 30 km west of Tokyo Bay, from 14 May to 15 June 2008. Doppler lidar on 27 May 2008 observed the vertical and horizontal wind structure of a well-developed sea-breeze front (SBF) penetrating from Sagami Bay. At the SBF, a strong updraft (maximum w approximately equal to 5 m s⁻¹) was formed with a horizontal scale of about 500 m and vertical scale of 2 km. The spatial relationship between the strong updraft over the nose of the SBF and prefrontal thermal suggests that the strong updraft was triggered by interaction between the SBF and the thermal. After the updraft commenced, a collocated ceilometer observed an intense aerosol backscatter up to 2 km above ground level. The observational results suggest that the near-surface denser aerosols trapped in the head region of the SBF escaped from the nose of the SBF and were then vertically transported up to the mixing height by the strong updraft at the SBF. This implies that these phenomena occurred not continuously but intermittently. The interaction situations between the SBF and prefrontal thermal can affect the wind structure at the SBF and the regional air quality.     

An Improved Ekman Layer Approximation for Smooth Eddy Diffusivity Profiles

The Ekman boundary-layer model is extended analytically for a gradually varying eddy diffusivity K(z) ≥ 0, z ≥ 0. A solution for the Ekman layer is provided having similar structure to the constant-K case; that is, exponentially decaying sine functions for the two horizontal wind components. The analytical asymptotic solution compares well with its numerical counterpart for various K(z). The result can be useful in theoretical studies such as Ekman pumping, for efficient estimation of the Ekman layer profiles in various analyses with near-neutral stratifications, or for a rapid initialization of mesoscale models.     

Turbulence Characteristics of the Shear-Free Convective Boundary Layer Driven by Heterogeneous Surface Heating

Large-eddy simulations (LESs) are employed to investigate the turbulence characteristics in the shear-free convective boundary layer (CBL) driven by heterogeneous surface heating. The patterns of surface heating are arranged as a chessboard with two different surface heat fluxes in the neighbouring patches, and the heterogeneity scale Λ in four different cases is taken as 1.2, 2.5, 5.0 and 10.0 km, respectively. The results are compared with those for the homogeneous case. The impact of the heterogeneity scale on the domain-averaged CBL characteristics, such as the profiles of the potential temperature and the heat flux, is not significant. However, different turbulence characteristics are induced by different heterogeneous surface heating. The greatest turbulent kinetic energy (TKE) is produced in the case with the largest heterogeneity scale, whilst the TKE in the other heterogeneous cases is close to that for the homogeneous case. This result indicates that the TKE is not enhanced unless the scale of the heterogeneous surface heating is large enough. The potential temperature variance is enhanced more significantly by a larger surface heterogeneity scale. But this effect diminishes with increasing CBL height, which implies that the turbulent eddy structures are changed during the CBL development. Analyses show that there are two types of organized turbulent eddies: one relates to the thermal circulations induced by the heterogeneous surface heating, whilst the other identifies with the inherent turbulent eddies (large eddies) induced by the free convection. At the early stage of the CBL development, the dominant scale of the organized turbulent eddies is controlled by the scale of the surface heterogeneity. With time increasing, the original pattern breaks up, and the vertical velocity eventually displays horizontal structures similar to those for the homogeneous heating case. It is found that after this transition, the values of λ/z _i (λ is the dominant horizontal scale of the turbulent eddies, z _i is the boundary-layer height) ≈1.6, which is just the aspect ratio of large eddies in the CBL.     

Observational study of surface ozone at an urban site in East China

In this study, we present the observational data of near surface ozone and some meteorological parameters during 2004, at an urban site (36°42′ N, 117°08′ E, 34.5 m a.s.l.) of Jinan, China. Hourly ozone concentrations exceeding the standard value of China, 100 ppbv, were observed for 65 h (in 23 days) from April to October, and values exceeding US NAAQS (National Ambient Air Quality Standard) for 1 h ozone, 120 ppbv, were observed for 15 h (in 7 days) from late May to early July. Ozone formation presented the phenomenon of “weekend effect”, especially in summer. Monthly variation of ozone coincided with temperature except for July and August. The low ozone levels in July and August may be due to the short sunshine duration and much rainfall during this period. Among these meteorological parameters, daily averaged ozone shows a significant correlation with temperature (r = 0.66) in the year and with relative humidity (r = − 0.75) in summer. Throughout the year, high ozone concentrations were mainly associated with the wind from 180 to 247.5°, while high ozone concentration seemed to have no obvious correlation with a given wind direction in summer. An anomalous nocturnal high ozone episode during 23–25 May 2004 was investigated. Growth fractions of ozone during the nighttime episode were 62.2% and 71.1% for 23 and 24 May, respectively. Synoptic analysis shows that favorable synoptic condition had presumably elevated the background ozone level in this region. Backward trajectory analysis shows that the increase of ozone concentration and the relatively constant high ozone concentrations during the night of May 23 might originate from the transport of ozone rich air mass above boundary layer. Transport of ozone from Yangtze Delta and East Central China might be a significant process for the high ozone level during night May 24 at Jinan.     

A numerical study of the convective boundary layer

Computations of the buoyantly unstable Ekman layer are performed at low Reynolds number. The turbulent fields are obtained directly by solving the three-dimensional time-dependent Navier-Stokes equations (using the Boussinesq approximation to account for buoyancy effects), and no turbulence model is needed. Two levels of heating are considered, one quite vigorous, the other more moderate. Statistics for the vigorously heated case are found to agree reasonably well with laboratory, field, and large-eddy simulation results, when Deardorff's mixed-layer scaling is used. No indication of large-scale longitudinal roll cells is found in this convection-dominated flow, for which the inversion height to Obukhov length scale ratio –z _i /L _*=26. However, when heating is more moderate (so that –z _i /L _*=2), evidence of coherent rolls is present. About 10% of the total turbulent kinetic energy and turbulent heat flux, and 20% of the Reynolds shear stress, are estimated to be a direct consequence of the observed cells.     

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})$ .     

Local Imbalance of Turbulent Kinetic Energy in the Surface Layer

We utilize experimental data collected in 2002 over an open field in Hanford, Washington, USA, to investigate the turbulent kinetic energy (TKE) budget in the atmospheric surface layer. The von Kármán constant was determined from the near-neutral wind profiles to be 0.36 ± 0.02 rather than the classical value of 0.4. The TKE budget was normalized and all terms were parameterized as functions of a stability parameter z/L, where z is the distance from the ground and L is the Obukhov length. The shear production followed the Businger–Dyer relation for −2 < z/L < 1. Contrary to the traditional Monin–Obukhov similarity theory (MOST), the shear, buoyancy and dissipation terms were found to be imbalanced due to a non-zero vertical transport over all stabilities. Motivated by this local imbalance, modified parameterizations of the dissipation and the turbulent transport were attempted and generated good agreement with the experimental data. Assuming stationarity and horizontal homogeneity, the pressure transport was estimated from the residual of the TKE budget.     

Flux-gradient profiles for momentum and heat over an urban surface

Summary In this paper, we evaluate the applicability of flux-gradient relationships for momentum and heat for urban boundary layers within the Monin-Obukhov similarity (MOS) theory framework. Although the theory is widely used for smooth wall boundary layers, it is not known how well the theory works for urban layers. To address this problem, we measured the vertical profiles of wind velocity, air temperature, and fluxes of heat and momentum over a residential area and compared the results to theory. The measurements were done above an urban canopy whose mean height z_h is 7.3 m. 3-D sonic anemometers and fine wire thermocouples were installed at 4 heights in the region 1.5z_h < z < 4z_h. We found the following: (1) The non-dimensional horizontal wind speed has good agreement with the stratified logarithmic profile predicted using the semi-empirical Monin-Obukov similarity (MOS) function, when it was scaled by the surface friction velocity that is derived from the shear stress extrapolated to the roof-top level. (2) The scaled gradient of horizontal wind speed followed a conventional semi-empirical function for a flat surface at a level (z/z_h = 2.9), whereas, in the vicinity of the canopy height was larger than the commonly-used empirical relationship. (3) The potential temperature profile above the canopy shows dependency on the atmospheric stability and the scaled gradient of temperature is in good agreement with a conventional shear function for heat. In the case of heat, the dependency on height was not found. (4) The flux-gradient relationship for momentum and heat in the region 1.5z_h < z < 4z_h was rather similar to that for flat surfaces than that for vegetated canopies.