Simulating Dispersion in the Evening-Transition Boundary Layer

We investigate dispersion in the evening-transition boundary layer using large-eddy simulation (LES). In the LES, a particle model traces pollutant paths using a combination of the resolved flow velocities and a random displacement model to represent subgrid-scale motions. The LES is forced with both a sudden switch-off of the surface heat flux and also a more gradual observed evolution. The LES shows ‘lofting’ of plumes from near-surface releases in the pre-transition convective boundary layer; it also shows the subsequent ‘trapping’ of releases in the post-transition near-surface stable boundary layer and residual layer above. Given the paucity of observations for pollution dispersion in evening transitions, the LES proves a useful reference. We then use the LES to test and improve a one-dimensional Lagrangian Stochastic Model (LSM) such as is often used in practical dispersion studies. The LSM used here includes both time-varying and skewed turbulence statistics. It is forced with the vertical velocity variance, skewness and dissipation from the LES for particle releases at various heights and times in the evening transition. The LSM plume spreads are significantly larger than those from the LES in the post-transition stable boundary-layer trapping regime. The forcing from the LES was thus insufficient to constrain the plume evolution, and inclusion of the significant stratification effects was required. In the so-called modified LSM, a correction to the vertical velocity variance was included to represent the effect of stable stratification and the consequent presence of wave-like motions. The modified LSM shows improved trapping of particles in the post-transition stable boundary layer.     

Effect of Vertical Wind Shear on Concentration Fluctuation Statistics in a Point Source Plume

Measurements of concentration fluctuation intensity, intermittency factor, and integral time scale were made in a water channel for a plume dispersing in a well-developed, rough surface, neutrally stable, boundary layer, and in grid-generated turbulence with no mean velocity shear. The water-channel simulations apply to full-scale atmospheric plumes with very short averaging times, on the order of 1–4 min, because plume meandering was suppressed by the water-channel side walls. High spatial and temporal resolution vertical and crosswind profiles of fluctuations in the plume were obtained using a linescan camera laser-induced dye tracer fluorescence technique. A semi-empirical algebraic mean velocity shear history model was developed to predict these concentration statistics. This shear history concentration fluctuation model requires only a minimal set of parameters to be known: atmospheric stability, surface roughness, vertical velocity profile, and vertical and crosswind plume spreads. The universal shear history parameter used was the mean velocity shear normalized by surface friction velocity, plume travel time, and local mean wind speed. The reference height at which this non-dimensional shear history was calculated was important, because both the source and the receptor positions influence the history of particles passing through the receptor position.     

On the dispersion parameters of plumes from tall stacks in a shoreline environment

Experimental data for buoyant plumes released from high sources into layers having little ambient turbulence show that plume dispersion parameters vary in a manner similar to that during initial plume rise. This is consistent with general plume rise theory. Dispersion of plumes from tall stacks in a shoreline environment where a thermal internal boundary layer is formed often demonstrates this behaviour.     

An Auxiliary Tool to Determine the Height of the Boundary Layer

Results from radiosoundings, performed both over land and over sea, show that the ascent rate of a radiosounding balloon, the vertical velocity of the balloon, can be used to determine the height of the boundary layer. In many cases the balloon has a higher ascent rate in the boundary layer and a lower, less variable, ascent rate above. The decrease in ascending velocity appears as a jump at the top of the boundary layer. Two examples of potential temperature profiles for unstable stratification and one profile for stable conditions are shown with the corresponding ascent rates. A comparison between the boundary-layer height determined from potential temperature profiles and from ascent rates is presented for a larger dataset. The different ascent rates of the balloon in the boundary layer and above can be explained by a decrease in drag on the balloon in combination with a lowering of the critical Reynolds number in the boundary layer caused by turbulence. Hence, by simply logging the time from release of a radiosonde, it is possible to obtain additional information that can be used to estimate the height of both the unstable and stable boundary layers.     

Impact of the Nocturnal Low-Level Jet and Orographic Waves on Turbulent Motions and Energy Fluxes in the Lower Atmospheric Boundary Layer

The nocturnal low-level jet (LLJ) and orographic (gravity) waves play an important role in the generation of turbulence and pollutant dispersion and can affect the energy production by wind turbines. Additionally, gravity waves have an influence on the local mixing and turbulence within the surface layer and the vertical flux of mass into the lower atmosphere. On 25 September 2017, during a field campaign, a persistent easterly LLJ and gravity waves were observed simultaneously in a coastal area in the north of France. We explore the variability of the wind speed, turbulent eddies, and turbulence kinetic energy in the time–frequency and space domain using an ultrasonic anemometer and a scanning wind lidar. The results reveal a significant enhancement of the turbulence-kinetic-energy dissipation (by?50%) due to gravity waves in the LLJ shear layer (below the jet core) during the period of wave propagation. Large magnitudes of zonal and vertical components of the shear stress (approximately 0.4 and 1.5 m² s^?2, respectively) are found during that period. Large eddies (scales of 110 to 280 m) matching the high-wind-speed regime are found to propagate the momentum downwards, which enhances the mass transport from the LLJ shear layer to the roughness layer. Furthermore, these large-scale eddies are associated with the crests while comparatively small-scale eddies are associated with the troughs of the gravity wave.

     

Turbulence and diffusion over inhomogeneous terrain

This paper provides an overview of some aspects of atmospheric boundary-layer dispersion processes over homogeneous and complex terrain. Special emphasis is placed on a discussion of the boundarylayer scaling regimes over homogeneous terrain and the characteristics of the dispersion processes associated with each of these regimes. The paper points out that vertical concentration profiles usually deviate substantially from a Gaussian distribution. The mean flow and turbulence over a low hill is dealt with, and in the inner layer the turbulence levels are increased due to the mean flow speed-up. In the outer layer the turbulence is modified by the rapid distortion effect. In a middle layer the turbulence is reduced due to the effect of a hill-induced streamline curvature. The paper concludes that the flow perturbations introduced by large-scale hills and valleys invalidate the use of simple approximations for describing atmospheric dispersion processes, and that it is necessary to utilize the full set of equations of motion.     

大气对流边界层中的涡漩结构

大气边界层中存在尺度从几百米到几十公里的大涡漩运动。它们在边界层中动量、热量、水汽等垂直输送中起重要作用。作者从边界层中对流和上部稳定层中波动相互作用的观点，发展得出大涡结构的对流波动理论。根据此理论，大涡的波谱构成主要由上、下层大气中风向、风速、层结以及两层之间的温度跃变等因素决定。本文根据卫星云图和天气资料分析了一次冷空气爆发流经暖洋面上形成云街、对流单体以及它们之间的相互演化的过程，并用对流波动理论，依据各阶段的大气条件计算出它们的波数构成，并得出了垂直速度、辐合带、界面扰动的分布，解释了云街、对流单体的形成、结构及相互转化的原因     

The Influence Of Vertical Wind Direction Shear On Dispersion In The Convective Boundary Layer, And Its Incorporation In Coastal Fumigation Models

The mean concentration distributionwithin a plume released from a point source in the atmosphericboundary layer can be greatly influenced by the systematic turningof wind with height (i.e. vertical wind direction shear). Such aninfluence includes a deflection of the plume centroid, with anassociated shearing of the vertical plume cross-section, and anenhancement of dispersion, in the horizontal plane. Wind directionshear is normally not accounted for in coastal fumigation models,although dispersion observations with shear acting as acontrolling parameter are not uncommon. A three-dimensionalLagrangian stochastic model is used to investigate the influenceof uniform wind direction shear on the diffusion of a point-sourceplume within the horizontally homogeneous convective boundarylayer, with the source located at the top of the boundary layer.Parameterisations are developed for the plume deflection andenhanced dispersion due to shear within the framework of aprobability density function (PDF) approach, and compared with theLagrangian model results. These parameterisations are thenincorporated into two applied coastal fumigation models: a PDFmodel, and a commonly used model that assumes uniform andinstantaneous mixing in the vertical direction. The PDF modelrepresents the vertical mixing process more realistically. A moreefficient version of the PDF model, which assumes a well-mixedconcentration distribution in the vertical at large times, isapplied to simulate sulfur dioxide data from the Kwinana CoastalFumigation Study. A comparison between the model results and thedata show that the model performs much better when the wind-sheareffects are included.     

Scaling of the Vertical Spreading of a Plume of a Passive Tracer in a Rough-Wall Neutral Boundary Layer

The influence of surface roughness on the dispersion of a passive scalar in a rough wall turbulent boundary layer has been studied using wind-tunnel experiments. The surface roughness was varied using different sizes of roughness elements, and different spacings between the elements. Vertical profiles of average concentration were measured at different distances downwind of the source, and the vertical spread of the plume was computed by fitting a double Gaussian profile to the data. An estimate of the integral length scale is derived from the turbulence characteristics of the boundary layer and is then used to scale the measured values of plume spread. This scaling reduces the variability in the data, confirming the validity of the model for the Lagrangian integral time scale, but does not remove it entirely. The scaled plume spreading shows significant differences from predictions of theoretical models both in the near and in the far field. In the region immediately downwind of the source this is due to the influence of the wake of the injector for which we have developed a simple model. In the far field we explain that the differences are mainly due to the absence of large-scale motions. Finally, further downwind of the source the scaled values of plume spread fall into two distinct groups. It is suggested that the difference between the two groups may be related to the lack of dynamical similarity between the boundary-layer flows for varying surface roughness or to biased estimates of the plume spread.     

边界层急流型重力波─-飞机颠簸的一种形成机制

采用线性化的Boussinesq流体边界层绝热流动方程,比较一维边界层急流型重力波的垂直运动量级,讨论边界层急流型重力波中的湍流发展,认为边界层急流型重力波是造成边界层飞机颠簸的一种机制。     

边界层急流型重力波——飞机颠簸的一种形成机制

采用线性化的Boussinesq流体边界层绝热流动方程，比较一维边界层急流型重力波的垂直运动量级，讨论边界层急流型重力波中的湍流发展，认为边界层急流型重力波是造成边界层飞机颠簸的一种机制。     

Effect of wind vertical shear on diffusion characteristics in the mesoscale range

Using the turbulent statistical form of the vertical vortex diffusion coefficient K, in the planetary boundary layer (PBL) and Ekman spiral wind profile, the three-dimensional diffusion equation is solved by the numeri-cal method. The influences of vertical shear of both wind direction and wind speed on pollution trajectory and horizontal diffusion parameters σy are numerically analysed. The expressions of both pollution trajectory and σy, including the factor of wind shear, are obtained. The results show that the vertical shear of wind is important among all factors affecting the mesoscale dispersion. Specifically, from neutral to stable atmospheric conditions, vertical shear of wind makes greater contribution to σy than turbulence, thus it is the most important factor. In this paper, we have compared horizontal dispersion pattern with both Pasquill's dispersion pattern considering wind direction shear, and experimental data collected at 9 different sites rang-ing from 10 to 100 km, and the results show that our dispersion pattern is closer to the experimental values than Pasquill's results, and his correction to shear of wind direction is too large under the stable conditions.     

Large-Eddy Simulation of Katabatic Flows

A large-eddy simulation model with rotated coordinates and an open boundary is used to simulate the characteristics of katabatic flows over simple terrain. Experiments examine the effects of cross winds on the development of the slope-flow boundary layer for a steep (20°) slope and the role of drainage winds in preventing turbulence collapse on a gentle slope (1°). For the steep flow cases, comparisons between model average boundary-layer velocity, temperature deficit, and turbulence kinetic energy budget terms and tower observations show reasonable agreement. Results for different cross slope winds show that as the cross slope winds increase, the slope flow deepens faster and behaves more like a weakly stratified, sheared boundary layer. Analysis of the momentum budget shows that near the surface the flow is maintained by a balance between downslope buoyancy forcing and vertical turbulence flux from surface drag. Above the downslope jet, the turbulence vertical momentum flux reverses sign and acceleration of the flow by buoyancy is controlled by horizontal advection of slower moving ambient air. The turbulence budget is dominated by a balance between shear production and eddy dissipation, however, buoyancy and pressure transport both are significant in reducing the strength of turbulence above the jet. Results from the gentle slope case show that even a slight terrain variation can lead to significant drainage winds. Comparison of the gentle slope case with a flat terrain simulation indicates that drainage winds can effectively prevent the formation of very stable boundary layers, at least near the top of sloping terrain.     

A Large-Eddy Simulation and Lagrangian Stochastic Study of Heavy Particle Dispersion in the Convective Boundary Layer

Large-eddy simulation and Lagrangian stochastic dispersion models were used to study heavy particle dispersion in the convective boundary layer (CBL). The effects of various geostrophic winds, particle diameters, and subgrid-scale (SGS) turbulence were investigated. Results showed an obvious depression in the vertical dispersion of heavy particles in the CBL and major vertical stratification in the distribution of particle concentrations, relative to the passive dispersion. Stronger geostrophic winds tended to increase the dispersion of heavy particles in the lower CBL. The SGS turbulence, particularly near the surface, markedly influenced the dispersion of heavy particles in the CBL. For reference, simulations using passive particles were also conducted; these simulation results agreed well with results from previous convective tank experiments and numerical simulations.     

Turbulent transport from an arctic lead: A large-eddy simulation

The upward transfer of heat from ocean to atmosphere is examined for an Arctic lead, a break in the Arctic ice which allows contact between the cold atmosphere and the relatively warm ocean. We employ a large-eddy model to compute explicitly the three-dimensional turbulent response of the atmosphere to a lead of 200 m width. The surface heat flux creates a turbulent plume of individual quasi-random eddies, not a continuous updraft, which penetrate into the stable atmosphere and transport heat upward.Maximum updraft velocities and turbulence occur downwind of the lead rather than over the lead itself, because the development time of an individual thermal eddy is longer than its transit time across the lead. The affected vertical region, while shallow over the lead itself, grows to a height of 65m at 600 m downwind of the lead; beyond that, the depth of the turbulent region decreases as the eddies weaken. The maximum vertical turbulent heat flux occurs at the downwind edge of the lead, beyond which a relative maximum extends upward into the plume. Negative surface heat flux immediately downwind of the lead creates a growing stable layer, but above that internal boundary layer the turbulent heat flux is still positive. Updraft maxima are typically 28 cm/s, but compensating downdrafts result in time-averaged vertical velocities of less than 1 cm/s in the plume. Conditional sampling separates the updraft and downdraft contributions. Formulas for the horizontal eddy development distance and for the vertical plume penetration height are presented. The relative importance of mean and turbulent transport is compared for both vertical and horizontal heat transfer: turbulence dominates the vertical heat transport whereas mean advection dominates the horizontal transport, these offsetting transports producing a quasi-stationary state.     

VORTEX MOTIONS IN THE ATMOSPHERIC CONVECTIVE BOUNDARY LAYER

Large vortices with scales ranging from hundreds meters to tens of kilometers are generallyfound in the atmospheric convective boundary layer(CBL).These vortices play important roles in the vertical transport of momentum,heat,water vaporand other tracers in the boundary layer.On the basis of the view of interaction between theconvection in CBL and the gravity waves in the upper stable layer the authors developed aconvection-wave theory on the formation of large vortices.According to the theory thewavenumber spectrum of the large vortices mainly depends on the atmospheric conditions in bothof the upper and lower layers,such as wind speed,wind direction shear,stratification as well astemperature jump.In the present paper satellite image and weather data in a case of cold air outbreak over warmocean are analyzed to study every stage of the convective processes,such as cloud street,convective cell as well as their transformation.According to the theory the wavenumbercompositions for cloud street and convective cell are calculated,respectively,on the basis of theatmospheric conditions at every stage.The distributions of vertical motions,convergent band anddisturbed interface are obtained and compared with the cloud patterns in the convective processes.Thus the study seems to offer a likely explanation for the origin of large vortices in CBL.     

Radiative Processes in the Stable Boundary Layer: Part II. The Development of the Nocturnal Boundary Layer

The interaction between radiation and turbulence in the stable boundary layer over land is explored using an idealized model, with a focus on the surface layer after the evening transition. It is shown that finer vertical resolution is required in transitional boundary layers than in developed ones. In very light winds radiative cooling determines the temperature profile, even if similarity functions without a critical Richardson number are used; standard surface similarity theory applied over thick layers then yields poor forecasts of near-surface air temperatures. These points are illustrated with field data. Simulations of the developing nocturnal boundary layer are used to explore the wider role of radiation. Comparatively, radiation is less significant within the developed stable boundary layer than during the transition; although, as previous studies have found, it remains important towards the top of the stable layer and in the residual layer. Near the ground, reducing the surface emissivity below one is found to yield modest relative radiative warming rather than intense cooling, which reduces the potential importance of radiation in the developed surface layer. The profile of the radiative heating rate may be strongly dependent on other processes, leading to quite varied behaviour.     

Pollutant Plume Dispersion in the Atmospheric Boundary Layer over Idealized Urban Roughness

The Gaussian model of plume dispersion is commonly used for pollutant concentration estimates. However, its major parameters, dispersion coefficients, barely account for terrain configuration and surface roughness. Large-scale roughness elements (e.g. buildings in urban areas) can substantially modify the ground features together with the pollutant transport in the atmospheric boundary layer over urban roughness (also known as the urban boundary layer, UBL). This study is thus conceived to investigate how urban roughness affects the flow structure and vertical dispersion coefficient in the UBL. Large-eddy simulation (LES) is carried out to examine the plume dispersion from a ground-level pollutant (area) source over idealized street canyons for cross flows in neutral stratification. A range of building-height-to-street-width (aspect) ratios, covering the regimes of skimming flow, wake interference, and isolated roughness, is employed to control the surface roughness. Apart from the widely used aerodynamic resistance or roughness function, the friction factor is another suitable parameter that measures the drag imposed by urban roughness quantitatively. Previous results from laboratory experiments and mathematical modelling also support the aforementioned approach for both two- and three-dimensional roughness elements. Comparing the UBL plume behaviour, the LES results show that the pollutant dispersion strongly depends on the friction factor. Empirical studies reveal that the vertical dispersion coefficient increases with increasing friction factor in the skimming flow regime (lower resistance) but is more uniform in the regimes of wake interference and isolated roughness (higher resistance). Hence, it is proposed that the friction factor and flow regimes could be adopted concurrently for pollutant concentration estimate in the UBL over urban street canyons of different roughness.     

Forward-in-time and Backward-in-time Dispersion in the Convective Boundary Layer: the Concentration Footprint

The turbulence field obtained using a large-eddy simulation model is used to simulate particle dispersion in the convective boundary layer with both forward-in-time and backward-in-time modes. A Lagrangian stochastic model is used to treat subgrid-scale turbulence. Results of forward dispersion match both laboratory experiments and previous numerical studies for different release heights in the convective boundary layer. Results obtained from backward dispersion show obvious asymmetry when directly compared to results from forward dispersion. However, a direct comparison of forward and backward dispersion has no apparent physical meaning and might be misleading. Results of backward dispersion can be interpreted as three-dimensional or generalized concentration footprints, which indicate that sources in the entire boundary layer, not only sources at the surface, may influence a concentration measurement at a point. Footprints at four source heights in the convective boundary layer corresponding to four receptors are derived using forward and backward dispersion methods. The agreement among footprints derived with forward and backward methods illustrates the equivalence between both approaches. The paper shows explicitly that Lagrangian simulations can yield identical footprints using forward and backward methods in horizontally homogeneous turbulence.