A simple and fast model to compute concentration moments in a convective boundary layer

Recently, a modified meandering plume model for concentration fluctuations in a convective boundary layer has been formulated (Atmos. Environ. 34 (2000) 3599). This model is based on a hybrid Eulerian–Lagrangian approach and it accounts for the skewed and inhomogeneous turbulence characteristics of the convective flow. Using the same hypotheses, but eliminating the need for Lagrangian particle model, we propose a generalised approach, that only requires the knowledge of mean concentration field. The proposed model is independent from the method used to obtain the mean concentration field. The evaluation of the concentration fluctuation field needs a computational time of only few seconds on a standard PC. Therefore, the model is suitable for practical applications.     

Non-linear response of wet deposition to emissions reduction: A model study

An Eulerian atmospheric model with complex chemistry (Acidic Deposition and Oxidant Model) and a Lagrangian model with linear chemistry (Ontario Ministry of the Environment Trajectory Model) were used to simulate the wet SO₄²⁻ deposition pattern over eastern North America for 16 days during April 1981.The two model results agree reasonably well with each other when the 16 day average values are compared. They also show reasonable agreement with observed data. Having established the ability of the models to predict deposition patterns for 1981 emissions, reduction scenarios with 50% SO_x and 50% SO_x and NO_x of the 1981 emissions were studied through the Eulerian model. Near the heavy emissions area, the reduction in SO₄²⁻ wet deposition is only about 30–40%. In this respect the linear Lagrangian model departs significantly from the Eulerian model. This non-linearity in response is attributed to the role of oxidants in controlling the conversion of SO₂ to SO₄²⁻.     

Prediction of atmospheric dispersion of pollutants in an airport environment

In this article we discuss the development of a methodology to predict atmospheric turbulent dispersion of pollutants generated from air traffic in an airport. It is based on the Lagrangian stochastic model (LSM), developed by Das and Durbin [2005. A Lagrangian stochastic model for dispersion in stratified turbulence, Physics of Fluids 17, 025109]. The approach is via the backward trajectory formulation of the model. The sources and receptors in an airport type problem are modeled as spheres and procedures have been derived for concentration calculation by both forward and backward trajectory methods. Some tests are performed to highlight certain features of the method. The turbulence statistics that are required as input are provided in terms of similarity profiles. The airport domain is partitioned to make the required search algorithms efficient. Pollutant concentration profiles are calculated over a range of meteorological data.     

An efficient Lagrangian stochastic model of vertical dispersion in the convective boundary layer

We consider the one-dimensional case of vertical dispersion in the convective boundary layer (CBL) assuming that the turbulence field is stationary and horizontally homogeneous. The dispersion process is simulated by following Lagrangian trajectories of many independent tracer particles in the turbulent flow field, leading to a prediction of the mean concentration. The particle acceleration is determined using a stochastic differential equation, assuming that the joint evolution of the particle velocity and position is a Markov process. The equation consists of a deterministic term and a random term. While the formulation is standard, attention has been focused in recent years on various ways of calculating the deterministic term using the well-mixed condition incorporating the Fokker–Planck equation. Here we propose a simple parameterisation for the deterministic acceleration term by approximating it as a quadratic function of velocity. Such a function is shown to represent well the acceleration under moderate velocity skewness conditions observed in the CBL. The coefficients in the quadratic form are determined in terms of given turbulence statistics by directly integrating the Fokker–Planck equation. An advantage of this approach is that, unlike in existing Lagrangian stochastic models for the CBL, the use of the turbulence statistics up to the fourth order can be made without assuming any predefined form for the probability distribution function (PDF) of the velocity. The main strength of the model, however, lies in its simplicity and computational efficiency. The dispersion results obtained from the new model are compared with existing laboratory data as well as with those obtained from a more complex Lagrangian model in which the deterministic acceleration term is based on a bi-Gaussian velocity PDF. The comparison shows that the new model performs well.     

Validation of the operational emergency response model at the Swedish Meteorological and Hydrological Institute using data from etex and the chernobyl accident

The Eulerian atmospheric tracer transport model MATCH (Multiscale Atmospheric Transport and Chemistry model) has been extended with a Lagrangian particle model treating the initial dispersion of pollutants from point sources. The model has been implemented at the Swedish Meteorological and Hydrological Institute in an emergency response system for nuclear accidents and can be activated on short notice to provide forecast concentration and deposition fields.The model has been used to simulate the transport of the inert tracer released during the ETEX experiment and the transport and deposition of ¹³⁷Cs from the Chernobyl accident. Visual inspection of the results as well as statistical analysis shows that the extent, time of arrival and duration of the tracer cloud, is in good agreement with the observations for both cases, with a tendency towards over-prediction for the first ETEX release. For the Chernobyl case the simulated deposition pattern over Scandinavia and over Europe as a whole agrees with observations when observed precipitation is used in the simulation. When model calculated precipitation is used, the quality of the simulation is reduced significantly and the model fails to predict major features of the observed deposition field.     

Testing the importance of accurate meteorological input fields and parameterizations in atmospheric transport modelling using dream - validation against ETEX-1

A tracer model, the DREAM, which is based on a combination of a near-range Lagrangian model and a long-range Eulerian model, has been developed. The meteorological meso-scale model, MM5V1, is implemented as a meteorological driver for the tracer model. The model system is used for studying transport and dispersion of air pollutants caused by a single but strong source as, e.g. an accidental release from a nuclear power plant. The model system including the coupling of the Lagrangian model with the Eulerian model are described. Various simple and comprehensive parameterizations of the mixing height, the vertical dispersion, and different meterological input data have been implemented in the combined tracer model, and the model results have been validated against measurements from the ETEX-1 release. Several different statistical parameters have been used to estimate the differences between the parameterizations and meterological input data in order to find the best performing solution.     

Evaluation of a Gaussian and a Lagrangian model against a roadside data set,with emphasis on low wind speed conditions

The evaluation of the high percentiles of concentration distributions is required by most national air quality guidelines, as well as the EU directives. However, it is problematic to compute such high percentiles in stable, low wind speed or calm conditions. This study utilizes the results of a previous measurement campaign near a major road at Elimäki in southern Finland in 1995, a campaign specifically designed for model evaluation purposes. In this study, numerical simulations were performed with a Gaussian finite line source dispersion model CAR-FMI and a Lagrangian dispersion model GRAL, and model predictions were compared with the field measurements. In comparison with corresponding results presented previously in the literature, the agreement of measured and predicted data sets was good for both models considered, as measured using various statistical parameters. For instance, considering all NO_x data (N=587), the so-called index of agreement values varied from 0.76 to 0.87 and from 0.81 to 1.00 for the CAR-FMI and GRAL models, respectively. The CAR-FMI model tends to slightly overestimate the NO_x concentrations (fractional bias FB=+14%), while the GRAL model has a tendency to underestimate NO_x concentrations (FB=−16%). The GRAL model provides special treatment to account for enhanced horizontal dispersion in low wind speed conditions; while such adjustments have not been included in the CAR-FMI model. This type of Lagrangian model therefore predicts lower concentrations, in conditions of low wind speeds and stable stratification, in comparison with a standard Lagrangian model. In low wind speed conditions the meandering of the flow can be quite significant, leading to enhanced horizontal dispersion. We also analyzed the difference between the model predictions and measured data in terms of the wind speed and direction. The performance of the CAR-FMI model deteriorated as the wind direction approached a direction parallel to the road, and for the lowest wind speeds. However, the performance of the GRAL model varied less with wind speed and direction; the model simulated better the cases of low wind speed and those with the wind nearly parallel to the road.     

Study of expiratory droplet dispersion and transport using a new Eulerian modeling approach

Understanding of droplet nuclei dispersion and transport characteristics can provide more engineering strategies to control transmission of airborne diseases. Droplet dispersion in a room under the conventional well-mixed and displacement ventilation is simulated. Two droplet nuclei sizes, 0.01 and 10 μm, are selected as they represent very fine and coarse droplets. The flow field is modeled using k–ε RNG model. A new Eulerian drift-flux methodology is employed to model droplet phase. Under the conventional ventilation scheme, both fine and coarse droplets are homogeneously dispersed within approximately 50 s. Droplet nuclei exhibit distinctive dispersion behavior, particularly for low airflow microenvironment. After 270 s of droplet emission, gravitational settling influences the dispersion for 10 μm droplets, and concentration gradient can still be observed for displacement ventilation.     

Coastal and synoptic recirculation affecting air pollutants dispersion: A numerical study

This study examines the spatial distribution of potential recirculation over the East Mediterranean Sea, and the combined effect of synoptic and meso-scale recirculations on plume dispersion in the region. For this purpose, three case studies are performed by the RAMS–HYPACT modeling system, each for a different synoptic scale flow pattern. Both a quantitative measure of the recirculation potential at each grid cell and particle dispersion are calculated. Although the recirculation index is an Eulerian quantity for the wind field and plume dispersion is a manifestation of the Lagrangian behavior of the wind, good correlation is found between the two.Several locations are identified as having high recirculation potential, including southern Cyprus, the coasts of Israel and Lebanon, the eastern slopes of the Judean Mountains and the Haifa Bay in particular. In the latter location, high recirculation potential could be explained by strong interaction between the land–sea surfaces, curvature of the bay and proximity of the Carmel ridge. It is shown that the synoptic and meso-scale recirculations may, under certain conditions, act together and at the same time in determining particle distribution. Under weak synoptic scale flows, particles are recirculated over the entire East Mediterranean Sea basin, returning onshore after a period of 2–3 days to join freshly emitted particles. At the same time, near-shore land–sea breeze effects cause particles to recirculate on smaller time scales of less then one day, sometimes passing as much as three times over the same airshed. A single elevated emission source is shown to have the potential to impair air quality at a coastal strip as long as 100–200 km upon returning onshore.     

Examination of a subgrid-scale parameterization for the transport of pollutants in a nonprecipitating cumulus cloud ensemble

The mass flux based subgrid-scale parameterization technique of Gidel (1983, J. geophys. Res.88, 6587–6599) is re-examined for use in Eulerian long-range transport models. Specifically, the parameterization scheme is incorporated into the STEM-II Eulerian transport/transformation/removal model and the model is used to investigate pollutant transport in a nonprecipitating cumulus cloud ensemble. The effects of entrainment, detrainment, evaporation and the transport by subsidence, updrafts and turbulent diffusion are included in the analysis. Presented simulation results indicate that the parameterization is able to treat the rapid vertical transport by cloud updrafts, enables the calculation of reaction rates based on subgrid-scale concentrations, and is readily adopted by Eulerian models.     

Continuous four-dimensional source attribution for the Berlin area during two days in July 1994. Part I: the new Euler–Lagrange-model system LaMM5

The multiple nested three-dimensional mesoscale Eulerian grid point model MM5 was directly coupled with a Lagrangian particle trajectory model in order to perform a four-dimensional source attribution for the area of Berlin based on the import probability density (IPD) distribution of the according receptor box. Within the resulting model system LaMM5 introduced here, the IPD distributions are not based on backward trajectories, which lack the recognition of the turbulent environment, but on the forward integration of a huge amount (order of Million per day) of particle releases according to an emission scenario which is approximately continuous in space and time. Hence the receptor import record yields an accordingly continuous IPD distribution. Much attention has been paid on spatial and temporal resolution at the interface between both model parts (online-coupling) and the interface itself has been extended by the turbulent quantities resulting from the higher-order turbulence closure of the Eulerian model part. LaMM5 is applied on an episode with high photochemical activity across Berlin at two consecutive days (25th and 26th of July in 1994) with varying meteorological conditions leading to an accordingly different source attribution. The main results are:

• •The decay of the Berlin IPD with increasing source-receptor distance and time appears in an exponential manner if only sources out of a constant level (z=25 m) are regarded.
• •Heterogeneous wind fields in time and space enhance the contributions (emissions) of nearby sources to the total import of the receptor in contrast to stationary wind fields which increase the scope of the IPD distribution in upstream direction.

There are further results from several additional sensitivity studies presented in a companion paper B (Part II).     

Erratum to “Continuous four-dimensional source attribution for the Berlin area during two days in July 1994. Part I: The new Euler–Lagrange-model system LaMM5”: [Atmospheric Environment 35 (32) (2001) 5497]

The multiple nested three-dimensional (3D) mesoscale Eulerian grid point model MM5 was directly coupled with a Lagrangian particle trajectory model in order to perform a 4D source attribution for the area of Berlin based on the import probability density (IPD) distribution of the according receptor box. Within the resulting model system LaMM5 introduced here, the IPD distributions are not based on backward trajectories, which lack the recognition of the turbulent environment, but on the forward integration of a huge amount (order of million per day) of particle releases according to an emission scenario which is approximately continuous in space and time. Hence, the receptor import yields an accordingly continuous IPD distribution. Much attention has been paid on spatial and temporal resolution at the interface between both model parts (online-coupling) and the interface itself has been extended by the turbulent quantities resulting from the higher-order turbulence closure of the Eulerian model part. LaMM5 is applied on an episode with high photochemical activity across Berlin at two consecutive days (25 and 26 July 1994) with varying meteorological conditions leading to an accordingly different source attribution. The main results are:

• •The decay of the Berlin IPD with increasing source-receptor distance and time appears in an exponential manner if only sources out of a constant level (z=25 m) are regarded.
• •Heterogeneous wind fields in time and space enhance the contributions (emissions) of nearby sources to the total import of the receptor in contrast to stationary wind fields which increase the scope of the IPD distribution in upstream direction.

There are further results from several additional sensitivity studies presented in a companion paper B (Part II).     

Numerical modeling of dry deposition coupled to 22 photochemical reactions

Three Eulerian models for the dry deposition of photochemically reactive species were formulated and evaluated: a K-theory model with independent transport and deposition of each species, a K-theory model coupled with 22 gas-phase reactions, and a second-order flux-budget model coupled with 22 reactions. Operator splitting was used to separately solve the reaction and dispersion terms in the models including photochemistry. In an evaluation of numerical method performance, the Adams–Moulton method with a pseudo-steady-state approximation for the free radicals was found to be consistent with, but more computationally efficient than Gear’s method for the solution of the reaction terms. The sensitivity of profiles of vertical concentration and flux to the Eulerian model formulation varied according to the species’ Damköhler number. Although a K-theory model is adequate for weakly depositing species with Damköhler numbers <10^-3, a second-order flux-budget model is required for species with Damköhler numbers that exceed unity, such as nitrogen oxides, free radicals, and those sensitive to net flux production by chemical reactions. Selection of an appropriate model formulation for the entire system depended on the most reactive species. Simple K-theory models may not accurately predict dry deposition fluxes in the urban surface layer.     

Application of MM5 and CAMx4 to local scale dispersion of particulate matter for the city of Christchurch,New Zealand

A numerical model, Mesoscale Model version 5 (MM5), is used in conjunction with a three-dimensional Eulerian/Lagrangian dispersion model (CAMx4) to model PM₁₀ dispersion for a period of 48 h for the city of Christchurch, New Zealand. In a typical winter, Christchurch usually experiences severe degradation in air quality. The formation of a nocturnal temperature inversion layer during stagnant synoptic conditions, and the emissions of particulate matter (PM₁₀) mainly from solid fuel home heating appliances (the ‘Domestic’ factor) leads to severe smog episodes on about 30 nights each winter. The modelling results from the highest resolution computational grid are compared with observed meteorology and air pollution dispersion for winter 2000, when the Christchurch Air Pollution Study (CAPS2000) was underway. The numerical modelling system is able to simulate surface-layer meteorology and PM₁₀ spatial distribution with a good level of skill, with the Index of Agreement and Pearson's correlation coefficient greater than 0.8 for PM₁₀.     

A trajectory plume model for simulating air pollution transients

An air quality simulation model that is simple, yet capable of accurately estimating concentrations under unsteady meteorological conditions, has been developed. This trajectory plume model uses the Gaussian plume equation, but has an applicability that is approximately as wide as the Lagrangian puff model. The plume axis is represented by a series of straight-line plume segments. The performance of this model was evaluated by comparing it with other diffusion models. A comparison between simulation results using the present model and those using integrated puff and Eulerian diffusion models for three different metropolitan areas (one in Japan and two in the U.S.) has indicated that a simple trajectory plume model performs as well as the two other more complex models in simulating pollutant dispersion under complicated meteorological conditions such as those which occur during the transition period from a sea breeze to a land breeze.     

CFD modelling of small particle dispersion: The influence of the turbulence kinetic energy in the atmospheric boundary layer

When considering the modelling of small particle dispersion in the lower part of the Atmospheric Boundary Layer (ABL) using Reynolds Averaged Navier Stokes simulations, the particle paths depend on the velocity profile and on the turbulence kinetic energy, from which the fluctuating velocity components are derived to predict turbulent dispersion. It is therefore important to correctly reproduce the ABL, both for the velocity profile and the turbulence kinetic energy profile.For RANS simulations with the standard k–ε model, Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the k–ε turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46–47, 145–153.) proposed a set of boundary conditions which result in horizontally homogeneous profiles. The drawback of this method is that it assumes a constant profile of turbulence kinetic energy, which is not always consistent with field or wind tunnel measurements. Therefore, a method was developed which allows the modelling of a horizontally homogeneous turbulence kinetic energy profile that is varying with height.By comparing simulations performed with the proposed method to simulations performed with the boundary conditions described by Richards and Hoxey (1993. Appropriate boundary conditions for computational wind engineering models using the k–ε turbulence model. Journal of Wind Engineering and Industrial Aerodynamics 46–47, 145–153.), the influence of the turbulence kinetic energy on the dispersion of small particles over flat terrain is quantified.     

Removal of lead and zinc ions from aqueous solution using Amasya zeolites from Turkey

A 3-D Eulerian-Lagrangian approach to moving vehicles is presented that takes into account the traffic-induced flow rate and turbulence. The method is applied to pollutant dispersion in an individual street canyon and a system of two street canyons forming a perpendicular intersection. The approach is based on computational fluid dynamics (CFD) calculations using a Eulerian approach for continuous phase and a Lagrangian approach for moving vehicles. The wind speed was assigned values of 4, 7 and 12 m/s. One-way and two-way traffic with different traffic rates per lane is considered. In the case of the intersection, a longitudinal wind direction was assumed. Predictions show differences in the pollutant dispersion in the case of one-way and two-way traffic.     

Monte Carlo simulation of nitrogen oxides dispersion from a vehicular exhaust plume and its sensitivity studies

The pollutant dispersion behavior from the vehicular exhaust plume has a direct impact on human health, particularly to the drivers, bicyclists, motorcyclists, pedestrians, people working nearby and vehicle passengers. A two-dimensional pollutant dispersion numerical model was developed based on the joint-scalar probability density function (PDF) approach coupled with a k–ε turbulence model to simulate the initial dispersion process of nitrogen oxides, temperature and flow velocity distributions from a vehicular exhaust plume. A Monte Carlo algorithm was used to solve the PDF transport equations in order to obtain the dispersion distribution of nitrogen oxides concentration. The model was then validated by a series of sensitivity experimental studies in order to assess the effects of vehicular exhaust tailpipe velocities, wind speeds and chemistry on the initial dispersion of NO and NO₂ mass concentrations from the vehicular exhaust plume. The results show that the mass concentrations of nitrogen oxides decrease along the centerline of the vehicular exhaust plume in the downstream distance. The dispersion process can be enhanced when the vehicular exhaust tailpipe velocity is much larger than the wind speed. The oxidation reaction of NO plays an important role when the wind speed is large and the vehicular exhaust exit velocity is small, which leads to chemical reduction of NO, and the formation and accumulation of NO₂ in the exhaust plume. It is also found that the effect of vehicular exhaust-induced turbulence in the vicinity of the exhaust tailpipe exit is more dominant than the effect of wind turbulence, while the wind turbulence gradually shows a significant role for the dispersion of nitrogen oxides along with the development of exhaust plume. The range of dispersion of nitrogen oxides in the radial direction is increased along with the development of vehicular exhaust plume.     

Eulerian—Lagrangian relationships in monte carlo simulations of turbulent diffusion

A Monte Carlo technique for simulating turbulent diffusion is developed based on the Eulerian space-time velocity autocorrelation function. This method is an improvement over the usual Lagrangian approach since the Eulerian statistical properties of the turbulence are more easily measured. The analysis leads to some simple analytic relations between Eulerian and Lagrangian integral scales.