A puff pollutant dispersion model with wind shear and dynamic plume rise

A puff diffusion model, which includes wind shear and dynamic plume rise, is developed for numerical prediction of pollutant concentrations under unsteady and non-uniform flow conditions. The plume from a continuous source is treated as a series of puffs emitted successively from the source. Each puff is represented by a set of six tracer particles, which define the size, shape and location of the puff. Initially these particles are located at the surface of the source, on arbitrarily chosen orthogonal axes. The location of the particles is computed at each time step by taking into account advection, eddy diffusion, wind shear and entrainment of ambient air during plume rise. The concentration distribution of each puff is determined by fitting an ellipsoid to the cluster of the six particles and assuming a three-dimensional Gaussian distribution, with standard deviation equal to the half-lengths of the principal axes of the ellipsoid. The concentration at a point of interest is obtained by summing the contributions from nearby puffs. The effect of wind shear on the pollutant concentration is investigated by use of a typical wind shear encountered in the atmosphere. The results show that, at 600 m downstream from the source, the present model gives concentrations a factor of 2 higher and lower at one standard deviation below and above the plume center, respectively, than that of conventional models in which no wind shear is considered. The plume-rise formulation is calibrated against the observations compiled by Briggs and the model is used to predict the trajectory of a plume observed by Slawson and Csanady. Excellent agreement between the prediction and the observation can be achieved if an appropriate eddy diffusivity is chosen.     

Numerical and approximate solutions for plume rise

Numerical and approximate analytical solutions are compared for turbulent plume rise in a crosswind. The numerical solutions were calculated using the plume rise model of Hoult, Fay and Forney (1969, J. Air Pollut. Control Ass.19, 585–590), over a wide range of pertinent parameters. Some wind shear and elevated inversion effects are included. The numerical solutions are seen to agree with the approximate solutions over a fairly wide range of the parameters. For the conditions considered in the study, wind shear effects are seen to be quite small. A limited study was made of the penetration of elevated inversions by plumes. The results indicate the adequacy of a simple criterion proposed by Briggs (1969, AEC Critical Review Series, USAEC Division of Technical Information extension, Oak Ridge, Tennesse).     

A chemical model for urban plumes: Test for ozone and particulate sulfur formation in St. Louis urban plume

A chemical model with anthropogenic sources of nitrogen oxides and hydrocarbons is applied to simulate the chemical behaviour of pollutants in the St. Louis urban plume. It is suggested that a substantial increase in peroxy radical concentrations (HO₂,RO₂) in the polluted air mass outside its source region leads to an effective formation of secondary pollutants like ozone and sulfate particles. The model indicates characteristic time for ozone generation in the plume of a few hours. Maximum ozone mixing ratio of 115 ppb is predicted after 4 h transport time outside the source region. Conversion rates of SO₂ to H₂SO₄ through gas phase reactions with hydroxyl and peroxy radicals are estimated to be 1–5%h⁻¹. This leads to an approx 25% conversion of SO₂ to paniculate sulfur in the plume during the day. Agreement with measured ozone concentrations and flow rates of ozone and participate sulfur in the St. Louis plume on 18 July 1975 can be taken as strong indications that ozone and sulfate particle formation in the plume proceeds through the suggested mechanisms.     

Development and evaluation of the PRIME plume rise and building downwash model

A new Gaussian dispersion model, the Plume Rise Model Enhancements (PRIME), has been developed for plume rise and building downwash. PRIME considers the position of the stack relative to the building, streamline deflection near the building, and vertical wind speed shear and velocity deficit effects on plume rise. Within the wake created by a sharp-edged, rectangular building, PRIME explicitly calculates fields of turbulence intensity, wind speed, and streamline slope, which gradually decay to ambient values downwind of the building. The plume trajectory within these modified fields is estimated using a numerical plume rise model. A probability density function and an eddy diffusivity scheme are used for dispersion in the wake. A cavity module calculates the fraction of plume mass captured by and recirculated within the near wake. The captured plume is re-emitted to the far wake as a volume source and added to the uncaptured primary plume contribution to obtain the far wake concentrations. The modeling procedures currently recommended by the U.S. Environmental Protection Agency (EPA), using SCREEN and the Industrial Source Complex model (ISC), do not include these features. PRIME also avoids the discontinuities resulting from the different downwash modules within the current models and the reported overpredictions during light-wind speed, stable conditions. PRIME is intended for use in regulatory models. It was evaluated using data from a power plant measurement program, a tracer field study for a combustion turbine, and several wind-tunnel studies. PRIME performed as well as or better than ISC/SCREEN for nearly all of the comparisons.     

Effective stratification for plume rise and boundary layer growth computations

A method is proposed for evaluating the effective stratification compatible with well-known plume rise and entrainment growth formulations when the vertical distribution of density differs from idealized conditions. In addition consideration is given to reducing a measured non-constant profile of wind speed to one compatible with theory. A nomogram method based on the commonalities of the developments is offered for small volume data reduction.     

A model for the maximum credible hourly impact on any ground receptor from point sources with thermal plume rise

A pollutant dispersion model is developed, allowing fast evaluation of the maximum credible 1-h average concentration on any given ground-level receptor, along with the corresponding critical meteorological conditions (wind speed and stability class) for stacks with buoyant plumes in urban or rural areas. Site-specific meteorological data are not required, as the computed concentrations are maximized against all credible combinations of wind speed, stability class, and mixing height. The analysis is based on the dispersion relations of Pasquill-Gifford and Briggs for rural and urban settings, respectively, the buoyancy induced dispersion correlation of Pasquill, the wind profile exponent values suggested by Irwin, the buoyant plume rise relations of Briggs, as well as the Benkley and Schulman's model for the minimum mixing heights. The model is particularly suited for air pollution management studies, as it allows fast screening of the maximum impact on any selected receptor and evaluation of the ways to have this impact reduced. It is also suited for regulatory purposes, as it can be used to define the minimum stack size requirements for a given source as a function of the exit gas volume and temperature, the pollutant emission rates and their hourly concentration standards, as well as the source location relative to sensitive receptors.     

Modelling plume rise and dispersion from pool fires

This paper describes an investigation into the behaviour of smoke plumes from pool fires, and the subsequent generation of empirical models to predict plume rise and dispersion from such a combustion source. Synchronous video records of plumes were taken from a series of small-scale (0.06–0.25m²) outdoor methanol/toluene pool fire experiments, and used to produce sets of images from which plume dimensions could be derived. Three models were used as a basis for the multiple regression analysis of the data set, in order to produce new equations for improved prediction. Actual plume observations from a large (20.7 m×14.2 m) aviation fuel pool fire were also used to test the predictions. The two theoretically based models were found to give a better representation of plume rise and dispersion than the empirical model based on measurements of small-scale fires. It is concluded that theoretical models tested on small-scale fires (heat output ≈70 kW) can be used to predict plume behaviour from much larger combustion sources (heat output ≈70 MW) under near neutral atmospheric conditions.     

A hybrid plume model for local-scale atmospheric dispersion

This paper presents a local-scale dispersion model, based on atmospheric boundary layer scaling theory. In the vicinity of the source, Gaussian equations are used in both the horizontal and vertical directions. After a specified transition distance, gradient transfer theory is applied in the vertical direction, while the horizontal dispersion is still assumed to be Gaussian. The dispersion parameters and eddy diffusivity are modelled in a form, which facilitates the use of a meteorological pre-processor. We present a novel model of the vertical eddy diffusivity (K_z), which is a continuous function of height in various atmospheric scaling regions. The model also includes a treatment of the dry deposition of gases and particulate matter. The accuracy of the numerical model was analysed by comparing the model predictions with two analytical solutions; the numerical deviations from these solutions were less than 2% for the computational regime. The model has been tested against the Kincaid experimental field data. The agreement of the predictions and the data is good on the average, although the internal variation of the predictions versus data scatter plot is substantial.     

Validation of a Lagrangian model plume rise scheme using the Kincaid data set

Correct prediction of the initial rise of a plume due to momentum and buoyancy effects is an important factor in dispersion modelling. A new plume rise scheme, based upon conservation equations of mass, momentum and heat, for the Lagrangian model, NAME, is described. The conservation equations are consistent with the well-known analytical plume rise formulae for both momentum- and buoyancy-dominated plumes. The performance of the new scheme is assessed against data from the Kincaid field experiment. Results show that the new scheme adds value to the model and significantly outperforms the previous plume rise scheme. Using data from assessments of atmospheric dispersion models using the Kincaid data set, it is shown that NAME is comparable to other models over short ranges.     

A plume blight visibility model for regulatory use

The ERT visibility model has been developed over the past two years to estimate, principally for regulatory purposes, plume blight and visibility impairment resulting from point sources. The model combines the concentration fields calculated from single- or multiple-source dispersion models with calculations of the radiative transfer at four wavelengths in the visual range. It accounts for the scattering and/or absorption of the reactive pollutants NO₂, NO₃⁻, SO₄²⁻ and carbonaceous and non-carbonaceous TSP. It can consider any geometrically specified line of sight (e.g., to a specific vista from a point in a Class I area). The model gives the following information about plume blight visibility degradation: (1) visual range reduction; (2) plume contrast; (3) object contrast degradation; (4) plume discoloration, expressed as a blue/red intensity ratio and (5) object and sky discoloration, expressed as a ΔE value.     

Vertical emission profiles for Europe based on plume rise calculations

The vertical allocation of emissions has a major impact on results of Chemistry Transport Models. However, in Europe it is still common to use fixed vertical profiles based on rough estimates to determine the emission height of point sources. This publication introduces a set of new vertical profiles for the use in chemistry transport modeling that were created from hourly gridded emissions calculated by the SMOKE for Europe emission model. SMOKE uses plume rise calculations to determine effective emission heights. Out of more than 40 000 different vertical emission profiles 73 have been chosen by means of hierarchical cluster analysis. These profiles show large differences to those currently used in many emission models. Emissions from combustion processes are released in much lower altitudes while those from production processes are allocated to higher altitudes. The profiles have a high temporal and spatial variability which is not represented by currently used profiles.     

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.     

Evaluating the performance of the horizontal radial plume mapping technique for locating multiple plumes

This paper evaluated the feasibility of using the horizontal radial plume mapping (HRPM) technique to locate multiple emission sources via computational simulation. Seventy-two test maps, each having two Gaussian distributions, were generated in a two-dimensional domain. The HRPM technique with the non-negative least square (NNLS) algorithm was then applied to reconstruct the plumes, assuming a nine-beam scanning beam geometry. The NNLS algorithm successfully reconstructed the source locations of 68 of the 72 test maps. However, when one of the plumes was near the origin, the NNLS did not always identify the peak locations correctly. Furthermore, when the two plumes were spaced closely, the NNLS tended to reconstruct a wide plume covering both plumes instead of separating them due to the resolution limitation of the current nine-beam geometry. In the sensitivity analysis, five sets of random error (1%, 5%, 10%, 20%, and 30%) were added in the path-integrated concentration (PIC) from the 72 test maps, and thus, an additional 360 reconstructions were implemented. Robust results were obtained when the noise added was less than 20%. The results generally support the implementation of the NNLS algorithm in the HRPM technique as described in the U.S. Environmental Agency (EPA) Other Test Method 10 (OTM-10).

Implications: The methodology evaluated in this paper provides near-real-time estimates about the locations of multiple emission sources. The involved optical remote sensing instruments can monitor large spatial areas (e.g., landfills) in a cost-effective way.     

A fluctuating plume model for concentration fluctuations in a plant canopy

A one-particle Lagrangian model for continuous releases in the non-Gaussian inhomogeneous turbulence of a canopy layer is derived based on the fluctuating plume model of Franzese [2003. Lagrangian stochastic modeling of a fluctuating plume in the convective boundary layer. Atmos. Environ. 37, 1691–1701.]. The model equations are filtered by a time-dependent low-pass filter applied to the turbulent kinetic energy in order to obtain a fluctuating plume model able to simulate the vertical meandering of the cloud centroid through non-stationary Lagrangian equations. The model satisfies the well-mixed condition. The relative dispersion of particles and the concentration fluctuation statistics of a passive tracer inside a modeled vegetal canopy are studied. The probability density function of the concentration relative to the plume centroid is parameterized and the mean and variance fields of concentration are simulated and compared with wind-tunnel data and numerical simulations. A skewed, reflected probability density function for the vertical position of the plume centroid is considered.     

A new sampler for collecting separate dry and wet atmospheric depositions of trace organic chemicals

Studies conducted in Saskatchewan and elsewhere have demonstrated the atmospheric transport of agricultural pesticides and other organic contaminants and their deposition into aquatic ecosystems. To date these studies have focused on ambient concentrations in the atmosphere and in wet precipitation. To measure the dry deposition of organic chemicals, a new sampler was designed which uses a moving sheet of water to passively trap dry particles and gasses. The moving sheet of water drains into a reservoir and, during recirculation through the sampler, is passed through an XAD-2 resin column which adsorbs the trapped organic contaminants. All surfaces which contact the process water are stainless steel or Teflon. Chemicals collected can be related to airborne materials depositing into aquatic ecosystems. The sampler has received a United States patent (number 5,413,003 – 9 May 1996) with the Canadian patent pending.XAD-2 resin adsorption efficiencies for 10 or 50 μg fortifications of ten pesticides ranged from 76% for atrazine (2-chloro-4-ethylamino-6-isopropylamino-S-triazine) to 110% for triallate [S-(2,3,3-trichloro-2-phenyl)bis(1-methylethyl)carbamothioate], dicamba (2-methoxy-3,6-dichlorobenzoic acid) and toxaphene (chlorinated camphene mixture). Field testing using duplicate samplers showed good reproducibility and amounts trapped were consistent with those from high volume and bulk pan samplers located on the same site. Average atmospheric dry deposition rates of three chemicals, collected for 5 weeks in May and June, were: dicamba, 69 ng m^-2 da^-1; 2,4-D (2,4-dichlorophenoxyacetic acid), 276 ng m^-2 da^-1: and, γ-HCH (γ-1, 2, 3, 4, 5, 6-hexachlorocyclohexane), 327 ng m^-2 da^-1.