Current status of air quality models for particulate matter.

The current status of the mathematical modeling of atmospheric particulate matter (PM) is reviewed in this paper. Simulating PM requires treating various processes, including the formation of condensable species, the gas/ particle partitioning of condensable compounds, and in some cases, the evolution of the particle size distribution. The algorithms available to simulate these processes are reviewed and discussed. Eleven 3-dimensional (3-D) Eulerian air quality models for PM are reviewed in terms of their formulation and past applications. Results of past performance evaluations of 3-D Eulerian PM models are presented. Currently, 24-hr average PM2.5 concentrations appear to be predicted within 50% for urban-scale domains. However, there are compensating errors among individual particulate species. The lowest errors tend to be associated with SO4(2-), while NO3-, black carbon (BC), and organic carbon (OC) typically show larger errors due to uncertainties in emissions inventories and the prediction of the secondary OC fraction. Further improvements and performance evaluations are recommended.     

Different splitting techniques with application to air pollution models

Splitting techniques are commonly used when large-scale models, which appear in different fields of science and engineering, are treated numerically. Four types of splitting procedures are defined and discussed. The problem of the choice of a splitting procedure is investigated. Several numerical tests, by which the influence of the splitting errors on the accuracy of the results is studied, are given. It is shown that the splitting errors decrease linearly when (1) the splitting procedure is of first order and (2) the splitting errors are dominant. Three examples for splitting procedures used in all large-scale air pollution models are presented. Numerical results obtained by a particular air pollution model, Unified Danish Eulerian Model (UNI-DEM), are given and analysed.     

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.     

Comparison of a new Eulerian model with a modified Lagrangian approach for particle distribution and deposition indoors

Understanding of aerosol dispersion characteristics has many scientific and engineering applications. It is recognized that Eulerian or Lagrangian approach has its own merits and limitations. A new Eulerian model has been developed and it adopts a simplified drift–flux methodology in which external forces can be incorporated straightforwardly. A new near-wall treatment is applied to take into account the anisotropic turbulence for the modified Lagrangian model. In the present work, we present and compare both Eulerian and Lagrangian models to simulate particle dispersion in a small chamber. Results reveal that the standard k–ε Lagrangian model over-predicts particle deposition compared to the present turbulence-corrected Lagrangian approach. Prediction by the Eulerian model agrees well with the modified Lagrangian model.     

Modeling atmospheric mercury deposition in the vicinity of power plants

Two mathematical models of the atmospheric fate and transport of mercury (Hg), an Eulerian grid-based model and a Gaussian plume model, are used to calculate the atmospheric deposition of Hg in the vicinity (i.e., within 50 km) of five coal-fired power plants. The former is applied using two different horizontal resolutions: coarse (84 km) and fine (16.7 km). More than 96% of the power plant Hg emissions are calculated with the plume model to be transported beyond 50 km from the plants. The grid-based model predicts a lower fraction to be transported beyond 50 km: >91% with a coarse resolution and >95% with a fine resolution. The contribution of the power plant emissions to total Hg deposition within a radius of 50 km from the plants is calculated to be <8% with the plume model, <14% with the Eulerian model with a coarse resolution, and <10% with the Eulerian model with a fine resolution. The Eulerian grid-based model predicts greater local impacts than the plume model because of artificially enhanced vertical dispersion; the former predicts about twice as much Hg deposition as the latter when the area considered is commensurate with the resolution of the grid-based model. If one compares the local impacts for an area that is significantly less than the grid-based model resolution, then the grid-based model may predict lower local deposition than the plume model, because two compensating errors affect the results obtained with the grid-based model: initial dilution of the power plant emissions within one or more grid cells and enhanced vertical mixing to the ground.     

Modeling Atmospheric Mercury Deposition in the Vicinity of Power Plants

Abstract

Two mathematical models of the atmospheric fate and transport of mercury (Hg), an Eulerian grid–based model and a Gaussian plume model, are used to calculate the atmospheric deposition of Hg in the vicinity (i.e., within 50 km) of five coal–fired power plants. The former is applied using two different horizontal resolutions: coarse (84 km) and fine (16.7 km). More than 96% of the power plant Hg emissions are calculated with the plume model to be transported beyond 50 km from the plants. The grid–based model predicts a lower fraction to be transported beyond 50 km: >91% with a coarse resolution and >95% with a fine resolution. The contribution of the power plant emissions to total Hg deposition within a radius of 50 km from the plants is calculated to be <8% with the plume model, <14% with the Eulerian model with a coarse resolution, and <10% with the Eulerian model with a fine resolution. The Eulerian grid–based model predicts greater local impacts than the plume model because of artificially enhanced vertical dispersion; the former predicts about twice as much Hg deposition as the latter when the area considered is commensurate with the resolution of the grid–based model. If one compares the local impacts for an area that is significantly less than the grid–based model resolution, then the grid–based model may predict lower local deposition than the plume model, because two compensating errors affect the results obtained with the grid–based model: initial dilution of the power plant emissions within one or more grid cells and enhanced vertical mixing to the ground.     

Using prognostic model-generated meteorological output in the AERMOD dispersion model: an illustrative application in Philadelphia, PA

In this study, we introduce the prospect of using prognostic model-generated meteorological output as input to steady-state dispersion models by identifying possible advantages and disadvantages and by presenting a comparative analysis. Because output from prognostic meteorological models is now routinely available and is used for Eulerian and Lagrangian air quality modeling applications, we explore the possibility of using such data in lieu of traditional National Weather Service (NWS) data for dispersion models. We apply these data in an urban application where comparisons can be made between the two meteorological input data types. Using the U.S. Environment Protection Agency's American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model (AERMOD) air quality dispersion model, hourly and annual average concentrations of benzene are estimated for the Philadelphia, PA, area using both hourly MM5 model-generated meteorological output and meteorological data taken from the NWS site at the Philadelphia International Airport. Our intent is to stimulate a discussion of the relevant issues and inspire future work that examines many of the questions raised in this paper.     

Air quality modeling of hazardous pollutants: current status and future directions

The current requirements and status of air quality modeling of hazardous pollutants are reviewed. Many applications require the ability to predict the local impacts from industrial sources or large roadways as needed for community health characterization and evaluating environmental justice concerns. Such local-scale modeling assessments can be performed by using Gaussian dispersion models. However, these models have a limited ability to handle chemical transformations. A new generation of Eulerian grid-based models is now capable of comprehensively treating transport and chemical transformations of air toxics. However, they typically have coarse spatial resolution, and their computational requirements increase dramatically with finer spatial resolution. The authors present and discuss possible advanced approaches that can combine the grid-based models with local-scale information.     

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.     

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.     

Interregional Exchanges of Air Pollution: Model Types and Applications

Regional scale air quality simulation models covering spatial scales of thousands of kilometers are finding increasing applications in studies of acid deposition and other air pollution problems. The purpose of this paper is to familiarize the nonexpert with the characteristics of the major types of interregional air quality models currently in use: Eulerian grid, statistical trajectory, and Lagrangian trajectory. The basic features, advantages, and disadvantages of each of these modeling approaches are summarized, as are the important limitations and problems associated with interregional modeling in general. Typical applications are illustrated using examples from the use of a representative Lagrangian trajectory model, ENAMAP, over the eastern North American area.     

Chemical Leasing business models and corporate social responsibility

Chemical Leasing is a service-oriented business model that shifts the focus from increasing sales volume of chemicals towards a value-added approach. Recent pilot projects have shown the economic benefits of introducing Chemical Leasing business models in a broad range of sectors. A decade after its introduction, the promotion of Chemical Leasing is still predominantly done by the public sector and international organizations. We show in this paper that awareness-raising activities to disseminate information on this innovative business model mainly focus on the economic benefits. We argue that selling Chemical Leasing business models solely on the grounds of economic and ecological considerations falls short of branding it as a corporate social responsibility initiative, which, for this paper, is defined as a stakeholder-oriented concept that extends beyond the organization’s boundaries and is driven by an ethical understanding of the organization’s responsibility for the impact of its business activities. For the analysis of Chemical Leasing business models, we introduce two case studies from the water purification and metal degreasing fields, focusing on employees and local communities as two specific stakeholder groups of the company introducing Chemical Leasing. The paper seeks to demonstrate that Chemical Leasing business models can be branded as a corporate social responsibility initiative by outlining the vast potential of Chemical Leasing to improve occupational health and safety and to strengthen the ability of companies to protect the environment from the adverse effects of the chemicals they apply.     

Numerical validation of a Eulerian hydrochemical code using a 1D multisolute mass transport system involving heterogeneous kinetically controlled reactions

It is demonstrated that at steady state, the 1D thermo-kinetic hydrochemical Eulerian mass balance equations in pure advective mode are indeed identical to the governing mass balance equations of a single reaction path (or geochemical) code in open system mode. Thus, both calculated reaction paths should be theoretically identical whatever the chemical complexity of the water-rock system (i.e., multicomponent, multireaction zones kinetically and equilibrium-controlled). We propose to use this property to numerically test the thermo-kinetic hydrochemical Eulerian codes and we employ it to verify the algorithm of the 1D finite difference code KIRMAT. Compared to the other methods to perform such numerical tests (i.e., comparisons with analytical, semi-analytical solutions, between two Eulerian hydrochemical codes), the advantage of this new method is the absence of constraints on the chemical complexity of the modelled water-rock systems. Moreover, the same thermo-kinetic databases and geochemical functions can be easily and mechanically used in both calculations, when the numerical reference comes from the Eulerian code with no transport terms (u and D = 0) and modify to be consistent with the definition of the open system mode in geochemical modelling. The ability of KIRMAT to treat multicomponent pure advective transport, subjected to several kinetically equilibrium-controlled dissolution and precipitation reactions, and to track their boundaries has been successfully verified with the property of interest. The required numerical validation of the reference calculations is bypassed in developing the Eulerian code from an already checked single reaction path code. A forward time-upstream weighting scheme (a mixing cell scheme) is used in this study. An appropriate choice of grid spacing allows to calculate within the grid size uncertainty the correct mineral reaction zone boundaries, despite the presence of numerical dispersion. Its correction enables us to improve the convergence and to extend the numerical test to mixed advective-dispersive mass transport. However, the skewness factor involves numerical oscillations that prevent to compute different grid spacing. The use of a different chemically controlled time step constraint in both calculations induces some inconsistencies into the validation tests. This numerical validation method may be applied as well as to check a thermo-kinetic hydrochemical finite element based code, from a 1D heterogeneous systems, and 2D-3D systems provided that they are designed so as to be 1D equivalent. A one-step algorithm and the use of a numerical reference coming from the Eulerian code to be tested ensure the potential success (accuracy) of the numerical validation method.     

A Lagrangian approach for modelling turbulent transport and chemistry

A Lagrangian two-particle model for the relative diffusion and mixing of two reactive species is proposed. The model has been tested for consistency in simple geometrical configurations and has been compared to experimental data, Eulerian first- and second-order closure models as well as to other Lagrangian models. The model includes the covariance between the two species and therefore gives better predictions of both diffusion and chemistry than first-order closure models or the one-particle model proposed by Chock and Winkler (1994a,b, Journal of Geophysical Research, 99 D1, 1019–1031, 99 D1, 1033–1041). The model is a generalisation of the model presented by Komori et al. (1991, Journal of Fluid Mechanics 228, 629–659). The results of the calculations, although preliminary in character, indicate that the proposed algorithm is robust and efficient, and yields satisfactory results in turbulence with slow to moderate chemistry.     

Current Status of Air Quality Models for Particulate Matter

ABSTRACT

The current status of the mathematical modeling of atmospheric particulate matter (PM) is reviewed in this paper. Simulating PM requires treating various processes, including the formation of condensable species, the gas/ particle partitioning of condensable compounds, and in some cases, the evolution of the particle size distribution. The algorithms available to simulate these processes are reviewed and discussed. Eleven 3-dimensional (3-D) Eulerian air quality models for PM are reviewed in terms of their formulation and past applications. Results of past performance evaluations of 3-D Eulerian PM models are presented. Currently, 24-hr average PM_2.5 concentrations appear to be predicted within 50% for urban-scale domains. However, there are compensating errors among individual particulate species. The lowest errors tend to be associated with SO₄ ^2-, while NO₃ ^-, black carbon (BC), and organic carbon (OC) typically show larger errors due to uncertainties in emissions inventories and the prediction of the secondary OC fraction. Further improvements and performance evaluations are recommended.     

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₄²⁻.     

湖泊富营养化模型的研究进展

湖泊的富营养化是全球普遍关注的环境问题之一.湖泊的富营养化模型是防治、修复和治理湖泊富营养化的重要决策工具.按研究的侧重点不同,将湖泊富营养化模型分为简单回归模型、水质模型、生态模型和生态-水动力水质模型,并分别回顾了四类模型的研究进展.最后指出湖泊富营养化模型的发展趋势,强调不确定理论、3S技术、耦合模型是今后湖泊富营养化模型研究的重点,应在此基础上建立通用的模拟、预测、评价和优化模型,为湖泊富营养化管理提供科学依据.     

湖泊富营养化模型的研究进展

湖泊的富营养化是全球普遍关注的环境问题之一.湖泊的富营养化模型是防治、修复和治理湖泊富营养化的重要决策工具.按研究的侧重点不同,将湖泊富营养化模型分为简单回归模型、水质模型、生态模型和生态-水动力水质模型,并分别回顾了四类模型的研究进展.最后指出湖泊富营养化模型的发展趋势,强调不确定理论、3S技术、耦合模型是今后湖泊富营养化模型研究的重点,应在此基础上建立通用的模拟、预测、评价和优化模型,为湖泊富营养化管理提供科学依据.     

Mesoscale modeling of wintertime particulate matter episodes in eastern Washington, USA

The Spokane, Washington area is classified as a non-attainment area for the 24-h PM₁₀ standard due to a history of high particulate matter concentrations. A Eulerian regional air quality model (CALMET/CALGRID) has been used to characterize the emission, transport and dispersion of PM₁₀ and PM_2.5 in Spokane. Observations from a residential site (Rockwood, RW) and an industrial site (Crown Zellerbach, CZ), spanning July 1994–August 1996 were used to evaluate the current emission inventory. Two major tasks were devised to conduct the objectives of this investigation. First, a simple and efficient urban dispersion model (WYNDValley) was used to simulate important episodes characterized by the highest PM₁₀ and PM_2.5 concentrations. The selected episodes included four days with wet conditions for which no roads would have been emitting and seven days with dry conditions for which roads would emit. In the second step, a single road-emitting event was selected from the previous predicted results for further analysis using the Eulerian regional air quality model to examine the emission inventory. The urban and regional models predicted the observed concentration distributions reasonably well for the source emissions inventoried in Spokane. The mass concentrations of PM₁₀ were well predicted for the roads emitting case examined by both models indicating that the emission inventory based primarily upon area sources including roads is reasonably well characterized, at least at the RW site. The area sources around CZ are less well characterized, so that the PM₁₀ concentrations are underpredicted at CZ. The models appear unable to reach an equilibrium mass balance status at the beginning of the simulation, and the urban model seems unable to properly resolve the nocturnal boundary layer.     

The MICS-Asia study: model intercomparison of long-range transport and sulfur deposition in East Asia

An intercomparison study involving eight long-range transport models for sulfur deposition in East Asia has been initiated. The participating models included Eulerian and Lagrangian frameworks, with a wide variety of vertical resolutions and numerical approaches. Results from this study, in which models used common data sets for emissions, meteorology, and dry, wet and chemical conversion rates, are reported and discussed. Model results for sulfur dioxide and sulfate concentrations, wet deposition amounts, for the period January and May 1993, are compared with observed quantities at 18 surface sites in East Asia. At many sites the ensemble of models is found to have high skill in predicting observed quantities. At other sites all models show poor predictive capabilities. Source–receptor relationships estimated by the models are also compared. The models show a high degree of consistency in identifying the main source–receptor relationships, as well as in the relative contributions of wet/dry pathways for removal. But at some locations estimated deposition amounts can vary by a factor or 5. The influence of model structure and parameters on model performance is discussed. The main factors determining the deposition fields are the emissions and underlying meteorological fields. Model structure in terms of vertical resolution is found to be more important than the parameterizations used for chemical conversion and removal, as these processes are highly coupled and often work in compensating directions.