EPA/600/A-97/078
           SCIENCE FEATURES IN MODELS-3 COMMUNITY MULTISCALE AIR QUALITY SYSTEM

  Jason Ching1', Daewon Byun1', Jeff Young1', Francis S. Binkowski1', Jonathan Pleim1', Shawn Roselle1', James Godowitch1',
                                         William Benjey1*, and Gerald Gipson2
                           'Atmospheric Sciences Modeling Division, Air Resources Laboratory,
                                   National Oceanic and Atmospheric Administration,
                                          Research Triangle Park, NC 27711
               2 Human Exposure and Atmospheric Sciences Division, National Exposure Research Laboratory,
                         U.S. Environmental Protection Agency, Research Triangle Park, NC 27711
 1.  INTRODUCTION

         Air quality simulation models are important tools for
 use by regulatory, policy and research communities. The Clean
 Air Act provides a societal mandate to assess and to manage air
 pollution levels to protect human health and the environment.
 The U.S. Environmental Protection Agency (USEPA) has
 established National Ambient Air Quality Standards (NAAQS),
 requiring  the development  of effective emissions  control
 strategies for such pollutants as ozone, paniculate matter and
 nitrogen species. National and regional policies are needed for
 reducing and managing the amount and type of emissions that
 cause acid, nutrient and toxic pollutant deposition to ecosystems
 at risk and for enhancing the visual quality of the environment.
 Air quality models are  used to develop emission  control
 strategies that achieve these objectives. Control strategies must
 be  both  environmentally  protective  and   cost  effective.
 However, for  effectiveness,  one must recognize that  air
 pollution problems and strategies for their mitigation are very
 complex and the linkages between sources, meteorology and
 natural sources and landscapes are highly varied, complex and
 not very well understood.  The goal of developing cost-effective
 control  strategies  is  challenging  and  is best considered
 holistically. The effectiveness of any control  strategy is very
 limiting when air pollution issues are  handled in isolation.
Emissions from chemical, manufacturing, and other industrial
 activities,  power generation, transportation and waste treatment
activities  contribute  to  a variety of  air pollution issues,
 including ozone, particulate matter (PM), acid, nutrient and
toxic deposition, and  visibility in complex ways, and at' a
 variety of spatial and temporal scales. The residence times of
pollutants in  the atmosphere can extend to multi day, thus,
transport  consideration must be at  least regional in scale.
NAAQS requirements and other goals  for a cleaner environment
vary over a large range of time scales, from  peak hourly to
annual averages.
Corresponding author address: Jason K.S. Ching, MD-80.  'On
assignment to National Exposure Research Laboratory, U.S.
Environmental Protection Agency, Research Triangle Park, NC
27711; email: ching.jason@epamail.epa.gov.
         Development of air quality simulation models started
in the late seventies.  The  Urban Airshed  Model (UAM)
followed by the Regional Oxidant Model  (ROM) provided
Eulerian-based models for ozone, the former for urban and the
latter for regional scale.  Strategies for State Implementation
Plans (SIPS) used ROM to provide boundary conditions for
UAM simulations.  Attention to acid deposition issues was
addressed in the eighties with the development and evaluation
of regional acid deposition models such as the Regional Acid
Deposition Model (RADM), the Acid Deposition and Oxidant
Model (ADOM), and the Sulfur Transport and Emissions Model
(STEM). Other major modeling systems included the Regional
Lagrangian Modeling of Air Pollution model  (RELMAP),  a
Lagrangian  framework  system,  and  semi-empirical  and
statistical models.  Models of this period were designed to
address  specific air  pollution issues, such as ozone or acid
deposition.  Thus, flexibility to deal with other issues such as
particulate matter or toxics was very limited.  With the passage
of the Clean Air Act Amendments of 1990 (CAAA-90), a wide
range of additional issues was identified including visibility,
and fine and coarse  particles, as well as indirect exposure to
toxic pollutants such as  heavy  metals, semi-volatile organic
species, and  nutrient deposition to water bodies. These latter
issues will require multi-media models.
         To meet the challenges posed by the CAAA-90, the
USEPA embarked upon  the development of an advanced
modeling  framework,   Models-3,  designed for  holistic
environmental modeling utilizing state of science representation
of atmospheric processes in a high performance computing
environment. Descriptions of Models-3 can be found in Novak
et al. (1998) and Byun et al. (1998). The science components
in Models-3 are called the Community Multiscale Air Quality
(CMAQ) system. The Models-3/CMAQ system is designed as
a multi-pollutant, multi-scale Eulerian framework air quality
and atmospheric deposition modeling system. It contains state-
of-science parameterizations of atmospheric processes affecting
transport, transformation  and deposition of such pollutants as
ozone, particulate matter, airborne toxics, and  acidic  and
nutrient pollutant species.  With science in a continuing state of
advancement and review, it is an important design feature that
the Models-3 framework and CMAQ have the capability to
integrate and test future formulations in an efficient manner,
without needing to develop a completely new modeling system.
         In June 1998, the first release version (CMAQ-98)

-------
 will become available.  It will contain options representing
 different model descriptions of the major science processes.
 The science options available to the user include the gas phase
 chemistry mechanisms, RADM2 and CB-IV, a set of numerical
 solvers for the mechanisms, options for horizontal and vertical
 advection schemes, algorithms for fine and coarse particulate
 matter  predictions,  photolysis rates  and a  plume-in-grid
 approach. The purpose of this paper is to briefly list the key
 science processes in general and the specific modeling options
 included throughout  CMAQ-98, including the  emissions and
 meteorological processors and their interface processors. Also
 discussed are methods  to determine  and test the relative
 contribution of different atmospheric processes to the air quality
 predictions, important functions built into the CMAQ system.

 2. SCIENCE AND MODELING FEATURES OF CMAQ

          In this section, the science components and related
 features of the  CMAQ-98 system is briefly described. The
 multi-scale capabilities of Models-3/CMAQ are handled by a set
 of nested domains, each with successively finer resolution. The
 selection of grid resolution, domain size and area of interest is
 dependent on the application,  and on computational resources,
 A configuration of three levels of nests with grid resolutions of
 36-12-4 km will be  provided in CMAQ-98.  With such a
 configuration, the 36km grid resolution  resolves the regional
 scale concentration fields and provides boundary conditions for
 an  intermediate scale at 12  km resolutions to set up the
 transition to the 4km  scale to handle features on urban scales.

 2.1.  CMAO Preprocessors:

 2.1.1 Meteorology

         The  CMAQ's  generalized  coordinate  system
 accommodates  meteorological  data   produced  by  any
 comprehensive meteorology model.  Initial implementation uses
 the Fifth Generation  PSU/NCAR Mesoscale Model (MM5),
 (Grell  et  al.,  1993).  The  configuration for preliminary
 evaluation includes non-hydrostatic dynamics, four-dimensional
 data-assimilation (FDDA) of winds, temperature, and humidity,
the Kain-Fritsch convective parameterization scheme, and the
 high resolution  boundary layer. A set of four nested grids
 (108/36/12/4)km is used to achieve urban scale resolution which
 is consistent with continental scale dynamics.
         In addition to the standard community  version of
 MM5,  a more  sophisticated  land-surface model has been
 developed and coupled to the MM5 to improve simulation of
 surface fluxes, surface temperature and humidity, and planetary
 boundary-layer (PEL) development  (Pleim and Xiu, 1995).
The surface scheme includes explicit simulation of soil moisture
 and  vegetative  evapotranspiration.  The CMAQ  has been
 designed to take advantage of the additional surface and PBL
parameters  produced  by this scheme including bulk stomatal
 resistance which  is  then  used to  compute  chemical  dry
 deposition velocities.  The surface model is being added to the
 community MM5 and will be available as an option in future
releases from NCAR.
         CMAQ  uses  an interface processor to link the
meteorology and the chemistry-transport model (CTM) called
MCEP. As part of its functionality, MOT contains a generalized
coordinate capability. It also contains algorithms for treating dry
deposition,

2.1.2 Emissions

         CMAQ-98 will use an emission processing system
called Models-3 Emissions Processing and Projection System
(MEPPS).  The MEPPS  is a crucial part  of the Models-3
modeling framework.  The MEPPS will processes emission
inventory data, performs future projections (including control
scenarios run using the Models-3  study planner) as well as
pre-processing of data for use by the Models-3 CMAQ model.
The MEPPS' input Data PROcessor (MEDPRO) will be capable
of importing, performing QA  and  converting emission data,
including but  not limited  to  the  1985  National  Acid
Precipitation Assessment  Program Inventory,  1990 Interim
National Inventory, and  1990 National Emissions  Trends
Inventory. MEPPS estimates mobile sources using the Mobile
5a model, and biogenic emissions using the BEIS-2  model.
Point and county  level area source emissions data are taken
from inventories.  Emission data are tracked and reported by
source category code. Spatial gridding is  performed using
Arc/Info G1S software.  Meteorological data from MM5 via
MCIP  are used in modeling  mobile sources  and  biogenic
emissions. CMAQ-98 MEPPS will  have the capability  to
provide speciated emissions  needed for using CB-IV and
RADM2 chemistry mechanisms or user-defined modifications,
Future versions will be able to support other  mechanisms.

2.1.3 Boundary and Initial Conditions

CMAQ-98 will establish the boundary conditions for the nests
by first creating modeled fields at the coarse mode, either at the
108 or 36km resolution. Alternatively, future capabilities can
incorporate observed data fields such as from remote sensors
and satellite platforms. Initial conditions will  be obtained from
the end of spinup runs of several days.

2.2 The Chemistry-Transport Model (CTM).

         The science core of the CMAQ system is the CTM.
Science processors in  CTM will  include the  following
capabilities:

2.2,1 Chemistry.

CMAQ includes both the RADM2 and  the CB-IV gas phase
mechanisms.  In addition the CMAQ provides the capability to
edit the RADM2 and the CB-IV mechanisms or to import a
completely new chemical mechanism by means of a generalized
chemical mechanism reader.   This reader is external to the
CTM,  and can significantly simplify the task of altering the
chemical representation in a CTM. However, any change to the

-------
 fundamental representation of the organic species in these two
 mechanisms will correspondingly require substantial changes to
 MEPPS.  The CMAQ also accounts for the formation  of
 secondary aerosols and the reactions of pollutants in the
 aqueous phase. Secondary aerosol formation is parameterized
 on the basis of important gas-phase reaction rates and aqueous
 reactions  are  simulated by means of the aqueous chemical
 mechanism  incorporated in  RADM.  Both the CB-IV and
 RADM2 gas-phase mechanisms are linked to these processes,
 however,  providing  the capability  to simulate multi-phase
 interactions using either gas-phase mechanism. Finally, two
 chemistry solvers are available — the Sparse Matrix Vectorized
 Gear(SMVGEAR) algorithm developed by Jacobson and Turco
 (1994)  and  the  Quasi-Steady State Approximation (QSSA)
 method used in the Regional Oxidant Model.  SMVGEAR is
 generally recognized  as the more accurate of the two, but it is
 much slower than QSSA on non-vector computers.

 2.2.2 Cloud processes:

         Modeling clouds are essential in air quality modeling
 due to their critical role in atmospheric pollutant transport and
 chemistry processes. Clouds have both direct and  indirect
 effects on the  pollutant concentrations: they directly modify
 concentrations  via aqueous chemical reactions, vertical mixing,
 and  wet deposition removal processes, and indirectly affect
 concentrations  by altering radiative transmittances which affect
 photolysis rates  and  biogenic fluxes. CMAQ-98 will model
 deep and  shallow clouds using  the RADM-type (Walcek-
 Taylor) algorithms for 36 and 12 km; an explicit type cloud
 scheme  is  being tested for the 4 km simulations.

 2.2.3 Diffusion and Advection

         Options for computing  subgrid  vertical transport
 include eddy diffusion, and the Asymmetric Convective Model
 (ACM)  (Pleim and Chang,  1992), the  latter applicable to
 convective conditions.  Horizontal diffusion is limited to use of
 a constant eddy diffusion coefficient. Numerical methods differ
 in the handling of advection of concentration fields.  Several
methods are implemented in the CMAQ; these include the
 method by Smolarkewicz (1983), a scheme by Sort (1989), and
a piecewise parabolic method (PPM) (Collela and Woodward,
 1984). The Smolarkiewicz technique is mass conservative, and
it is based on the first order upstream method. The Bott scheme
CPU costs  are  comparable to  Smolarkiewicz's,  negative
concentration values are limited, and it is highly accurate and
mass conservative. The PPM subgrid distributions of advective
quantities are  represented by a parabola in each grid.  PPM
provides local  fit to data, is monotonic, and contains special
steepening procedures where sharp gradients exist.

2.2.4 Particle Modeling and Visibility:

         One  of the major advancements in  CMAQ is the
modeling of fine and coarse mode particles. The fine fraction
model   described in  Binkowski  and  Shankar  (1995)  is
 incorporated into CMAQ-98.   CMAQ will predict hourly
 gridded concentration of fine particle mass whose size is equal
 to  or  less  that 2.5 microns,  speciated to sulfate, nitrate,
 ammonium, organics and aerosol water. Secondary sulfate is
 produced by chemical reactions of hydroxyl radicals with sulfur
 dioxide producing sulfuric acid that either condenses to existing
 particles or nucleates to form new particles. Anthropogenic and
 biogenic hydrocarbon precursors  also  react with hydroxyl
 radicals, and with ozone and nitrates to produce condensible
 material. Aerosol water is equilibrated with relative humidity
 and the ammonium to sulfur molar ratios. The fine range is
 modeled with two log-normal size distributions;  the coarse
 fraction uses one log-normal distribution. Nuclei and Aitken
 particles are modeled in the smaller of the two fine particle
 modes. The second mode is called the accumulation mode.  At
 the point of overlap, particles from the nuclei mode are shifted
 to the accumulation mode when the Aitken-mode growth rate
 exceeds that of the accumulation mode. CMAQ model outputs
 also include number densities for both fine modes and for the
 coarse modes. CMAQ-98 will provide methods for computing
 particle nucleation, aerosol dry deposition, and cloud processes
 including the oxidation of Sulfur(IV) to Sulfur(VI) by hydrogen
 peroxide. Primary particles enter from MEPPS.
         The modeling of  aerosols in CMAQ  provides the
 capability to handle other environmental  issues including
 visibility and semi-volatile air toxics. Visibility (in Deciview
 units) is output for each grid and model time step. The visibility
 formulation  in  CMAQ-98 is based on integrating  the Mie
 scattering over the distribution of the predicted particles.  In
 another potential application, CMAQ can provide the basis for
 modeling the atmospheric transport and deposition of semi-
 volatile organic compounds (SVOC) with parameterizations for
their rates of condensation  to and/or volatilization from the
 modeled particles.

 2.2.5 Plume-in-Grid (PinG) modeling:
         CMAQ-98 will include algorithms to provide a more
 realistic treatment of the subgrid scale physical and chemical
 processes impacting pollutant species in plumes released from
 selected Major Elevated Point Source Emitters (MEPSEs). The
 key PinG components developed to treat the relevant processes
 at the proper spatial and temporal scales for pollutant plumes
 include  a Plume Dynamics  Model  (PDM) processor, which
simulates plume rise, vertical/horizontal growth, and plume
 position and a Lagrangian Reactive Plume Module (LRPM),
which simulates the relevant dynamic and chemical  reaction
 processes of subgrid plumes  until certain criteria  are met
triggering the handover of plume material to intercepted grid
cells.   PinG will  be used  for the 36km and the  12  km
simulations while at 4 km resolutions, PinG is not invoked as
the MEPSE emissions are directly released into the CTM 3-D
grid cells. A detailed description of the PinG features is given
 in Gillani et al. (1998).

2.2.6  Photolysis Modeling:

        The photochemistry of pollution is initiated by photo

-------
 dissociation of smog precursors driven by solar radiation. The
 amount of radiation available for photochemistry is clearly
 dependent on sun angle (time of day), season, latitude, and land
 surface characteristics, and  greatly affected by atmospheric
 scatterers and absorbers, as well as cloud fields. Additionally,
 photolytic rates are wavelength  and temperature  dependent,
 thus, making their calculation very complex, and requiring great
 accuracy. An  advanced photolysis rate model will be installed
 in CMAQ-98, with capabilities to predict temporally resolved
 3-D gridded photolyses rates using input data from the MM5-
 FDDA and other available sources. It will have flexibility in the
 specification  of wavelength bands, extraterrestrial irradiances,
 vertical profiles of ozone and aerosols, cloud distribution and
 adsorption cross section and quantum yield data.  The model
 will have the capability of computing photolysis rates for any
 chemical  mechanism  with a user specified absorption cross
 section and quantum yields data.  CMAQ-98 defaults have been
 set up for CB-IV and  RADM2 mechanisms. Studies  have
 shown great  sensitivities  to the modeled cloud  fields, and
 clearly photolysis rates will  be limited by the accuracy of the
 modeled cloud fields.

 2.3.  Process Analysis

         Process analysis involves examining the effects that
 different physical and chemical atmospheric processes have on
 pollutant concentrations.   It is accomplished  in a model by
 quantifying the contributions of individual processes to the
 overall change- in a species concentration.  This information
 explains how individual model predictions come about and
 reveals the relative importance of each process.  It can be used
 diagnostically to identify potential sources of error in either the
 model formulation  or the model input data.   Further, it is
 particularly useful for gaining understanding of the effects of
 making changes to the model  or to its  inputs.  CMAQ-98
 provides the capability to perform process analyses using two
 different pieces of information: Integrated Process Rates (IPRs)
 and Integrated Reaction Rates (IRRs).
         The IPRs are obtained during a model simulation by
 computing and saving the change in concentration of each
 species  caused by  physical processes including  advection,
 diffusion, emissions, etc. IPRs are also calculated for chemical
 reaction, aerosol production, and aqueous chemistry, but values
provide no information on the particulars occurring within each
 process (i.e.,  they  only give the net  effect of each process).
These IPRs can be output for each species at varying times to
show how the effects of each process vary both in time and in
 space.   Thus,  they are  particularly  useful for identifying
unexpectedly low or high process contributions which could be
 indicative of model errors.
         The IRR analysis  deals with  the  details of the
 chemical transformations.  CMAQ-98 currently provides the
 capability to  conduct IRR analyses  for gas-phase chemical
 species only,  but extensions to aqueous chemistry and aerosol
formation are planned. For gas-phase chemistry, the CTM has
 been instrumented to compute not only the concentration of
each species, but also to compute the integral of the individual
chemical reaction rates. The computation and output of the
IRRs are synchronized with the standard concentration field
calculations and outputs to permit one to study details of the
atmospheric transformations  and  their effects on  pollutant
concentrations. IRR analyses have been primarily applied to aid
in the understanding of ozone formation (Jeffries and Tonnesen,
1994).  For example, IRRs have been used to quantify such
chemical characteristics as OH and  NO chain lengths, the
amount of new radical production and termination, the amount
of ozone produced by and the yield of radical species from each
VOC species.  These quantities have typically been used to
understand the reasons for differences in model predictions
obtained with different chemical mechanisms.

2.4.  CMAQ Code Integration:

         The  codes  for CMAQ system integrated into the
Models-3  framework are in  FORTRAN.  The  integration
procedure was facilitated by following a small set of design,
coding  and  implementation  standards  that  include:  (1)
modularity,  which employs the concept of interchangeable
modules within process classes, allowing easy exchange of
time-splitting, science process solvers such as Bott and PPM;
(2) a standard subroutine interface at the module level. (Each
science module is conceived as operating on the gas and aerosol
species concentrations within the  domain grid space for  its
specified time period (model synchronization time step) in
sequence throughout the modeling scenario.); (3) restriction of
coding practices, which conceal data dependencies,  hinder
maintenance and foster hidden bugs, e.g., eliminating common
blocks, at least across modules; and (4) the Models-3 I/O-API
(ref 9), which  contains standardized file  I/O functions and a
modeler-friendly interface  built on  top of self-describing
netCDF (ref 10)  files that are portable across most,  if not  all
Unix platforms,
         The  CMAQ  codes have  been generalized  to
accommodate a range of different  computational coordinate
systems. The essential solution in making the code sufficiently
general is to relegate the problem to the MCIP which produces
data in the correct coordinate form for CMAQ. Once a user has
selected the modules and options,  building that version of a
CMAQ is made relatively easy.  The selected domain and
chemistry information is  encapsulated in parameter or data
statements that are automatically written to FORTRAN include
files, which in turn, are linked into  the code during the build
process to produce the CMAQ executable.  Execution of this
CMAQ requires the correct linkage with the meteorological and
emissions interface  processor data  to  ensure  the correct
generalized coordinate implementation.
3.
DISCUSSION
          A typical simulation of Models-3/CMAQ provides
hourly air quality fields for regional to urban scales for multi-
day episodes, typically of up to five days in duration.  The new
PM-fme  standard will  be based on annual  averages; thus,
utilization of CMAQ in this application will require aggregation

-------
 techniques. One such technique, initially developed for RADM
 wet deposition  applications,  was  recently  modified  and
 successfully applied to  fine paniculate matter by Eder and
 LeDuc (1996). This approach, which utilizes visibility data as
 a surrogate for  fine PM, will be applied to CMAQ on a
 continental  scale (i.e.,  contiguous United  States,  southern
 Canada and northern Mexico).  Future efforts will be needed to
 validate  this approach when a network of fine particulate
 samplers  is  deployed;   also,  aggregation  approaches  for
 mesoscale domains will need to be developed perhaps utilizing
 the method by Eder et al., (1994).
         The design basis for the CMAQ is its ability to adapt
 to science advances in atmospheric process modeling.  Thus
 there will be a continuing need for evaluating the performance
 and the veracity of the process modules as the CMAQ advances
 its science basis. The model evaluation activity for the CMAQ
 will be  staged  with the initial efforts to  show  relative
 performance against the  RADM, which itself has undergone
 extensive model evaluation efforts. Diagnostic evaluation will
 continue using-data bases from different regional studies such
 as the 1995 Southern Oxidant Study conducted in the vicinity
 of Nashville, TN and the NARSTO-NE study.
         CMAQ can be configured  for a wide range of
 applications from  science  studies  and  investigations to
 regulatory applications. An evaluated CMAQ can provide the
 benchmark from which more operational configurations can be
 referenced. It  is  anticipated  that as science  advances are
 provided to CMAQ, future configurations of a more operational
 nature can be periodically rebenched as appropriate.
         Finally, we encourage  full participation  by the
 scientific and modeling communities to engage in the growth
 and use of Models-3/CMAQ.  Models-3/CMAQ has flexibility
 for incorporating  scientific and modeling advances in its air
 quality science process  modules, for testing  of alternative
 science descriptions of processes, and for extending its current
 capability  to  handle   multimedia   environmental  issues.
 Additionally, the  community of users should  be  vigilant in
 performing evaluations and testing against improved databases
 and measurement technology to insure model comparability to
 the real world.
Comput.Phys. 54: 114-201.
(5) Eder, B.K., J.M. Davis and  P.  Bloomfeld,  1994: An
automatic classification scheme designed to better elucidate the
dependence of ozone on meteorology. J. Appl. Meteor, 33:
1182-1199.
(6) Eder, B.K. and S. LeDuc,  1996:  Can selected RADM
simulations be aggregated to estimate annual concentrations of
fine particulate matter?   Reprint from  the  llth Annual
International  Symposium on the Measurement of Toxic and
Related Pollutants, May 7-10,1996, RTF, NC 732-739
(7) Gillani, N.V., A. Biazar, Y. Wu, J. Godowitch, J. Ching, R.
Imhoff, 1998: The Plume-in-Grid treatment of major elevated
point source emissions in Models-3.  10th Joint Conference on
Applications  of Air Pollution Meteorology   with AWMA,
January 11-16, 1998, Phoenix, AZ,  Amer. Meteorol. Soc.,
Boston, MA.
(8)  Grell, G.A., J.  Dudhia, and D.R.  Stauffer, 1993:  A
description  of  the   Fifth-Generation  Penn   State/NCAR
Mesoscale Model(MM5). NCAR Technical Note, NCAR/TN-
398+IA
(9) http://www.iceis.mcnc.org/EDSS/ioapi/index.htm I/
(10) http://www.unidataucar.edu/packages/netcdf/
(11) Jacobson M. and R.P.  Turco, 1994:  SMVGEAR:  A
Sparse-Matrix, vectorized Gear code for atmospheric models.
Atmos. Environ. 28: 273-284.
(12) Jeffries, H. E. and S. Tonnesen, 1994.  "A comparison of
two photochemical reaction mechanisms using mass balance
and process analysis", Atm. Em, 28(18): 2991-3003.
(13) Novak, Joan, et al., 1998: Models-3: A unifying framework
for environmental modeling and assessments. Preprint Volume,
10th Joint AMS and AW&MA Conference on the Applications
of Air Pollution Meteorology, Phoenix, AZ, Jan 11-16,1998.
(14) Pleim J.E., and J.S.  Chang,  1992: A non-local closure
model  in the  convective boundary layer. Atm Environ., 26A:
965-981.
(15) Pleim, J.E., and A. Xiu, 1995: Development and testing of
a  surface flux and  planetary boundary  layer model  for
application in mesoscale models, J. Appl. Meteor., 34:16-32.
(16) Smolarkiewicz, P.K., 1983: A simple positive definition
advection scheme with small implicit diffusion. Man Wea Rev.,
Ill: 479-486.
4. REFERENCES:

(1)  Binkowski; F.S., and U.  Shankar,  1995: The Regional
Particulate Model: Part I. Models description and  preliminary
results. J.Geophys. Res., 100(D12): 26,191-26,209.
(2)  Byun, Daewon, et al., 1998: Description of the Models-3
Community  Multi-scale Air Quality(CMAQ)  Modeling
System.  Preprint Volume,  10th  Joint AMS  and AW&MA
Conference on the Applications of Air Pollution Meteorology,
Phoenix, AZ, Jan 11-16,  1998.
(3)  Bott, A., 1989:, A  positive  definite advection scheme
obtained by nonlinear renormalization of the advective fluxes,
Man. Wea.  Rev. 117: 1006-1015.
(4)  Colella, P., and P. L. Woodward, 1984:  The Piecewise
Parabolic Method (PPM) for gas-dynamical simulations,  J.
DISCLAIMER

         This paper has been reviewed in accordance with the
U.S.   Environmental   Protection   Agency's   peer  and
administrative review policies and approved for presentation
and publication.  Mention  of trade names or commercial
products does not constitute endorsement or recommendation
for use.

-------
                                 TECHNICAL  REPORT DATA
 1. REPORT NO.
   EPA/600/A-97/078
                             2,
 4. TITLE AND SUBTITLE

 Science  Features in Models-3  Community Multiscale
 Air Quality System
                                                             5.REPORT DATE
                                                             6.PERFORMING ORGANIZATION CODE
 7, AUTHOR(S)

 Jason Ching,  et al.
            8.PERFORMING ORGANIZATION REPORT NO.
 9. PERFORMING ORGANIZATION NAME AND ADDRESS

 Same as Block 12
                                                             10.PROGRAM ELEMENT NO.
                                                             11. CONTRACT/GRANT NO.
12.  SPONSORING AGENCY NAME AND ADDRESS

U.S.  Environmental  Protection Agency
Office of Research  and Development
National Exposure Research Laboratory
Research Triangle Park, NC 27711
            13.TYPE OF REPORT AND PERIOD COVERED

            Preprint, FY-98
            14. SPONSORING AGENCY CODS

            EPA/600/9
15. SUPPLEMENTARY NOTES
16. ABSTRACT

Models-3 is framework for environmental modeling utilizing state of science representation of
atmospheric processes in a high performance computing environment. The Community-based Multi
scale Mr Quality  (CMAQ) system is Models-3 air quality modeling system.  Models-3/CMAQ model
is designed as a multi-pollutant Eulerian grid regional to urban scale air quality and
atmospheric deposition modeling system.  It contains state-of-science parameterizations and
algorithms of the relevant and contributing atmospheric processes affecting transport,
transformation and deposition of sources contributing to such pollutant issues as ozone,
particulate matter, airborne toxics and deposition of acidic and nutrient pollutant species.
The release version of Models-3/CMAQ will feature science options available to the user
including the gas phase chemistry mechanisms, RADM2, and CB-IV a set of chemistry solvers,
options for horizontal and vertical advection schemes, algorithms for fine and coarse
particulate matter, photolysis rates and a plume-in-grid approach.   These current science
options are discussed in the context of the overall structure of Models-3/CMAQ model which
includes an emission and meteorological preprocessors linked to the CMAQ by interface
processors. Also discussed are the  methods to- determine and test the relative contribution
of different atmospheric processes to the air quality predictions, important functions built
into the CMAQ system.
17.
                                 KEY WORDS AND DOCUMENT ANALYSIS
                  DESCRIPTORS
b.IDENTIFIERS/ OPEN ENDED
TERMS
                                                                             C. COSATI
18.  DISTRIBUTION STATEMENT
RELEASE TO PUBLIC
19.  SECURITY CLASS (This
Report)

UNCLASSIFIED
                                                                             21.NO.  OF PAGES
                                                 20. SECURITY CLASS (This
                                                 Page)

                                                 UNCLASSIFIED
                            22. PRICE

-------