&EPA
United Slifcs
Envirwiirwilfll Piotoctnn
Agflncy
COMPARISON OF REGULATORY
DESIGN CONCENTRATIONS
AERMOD
vs
ISCST3, CTDMPLUS, ISC-PRIME

-------
                                         EPA-454/R-03-002
                                              June 2003
COMPARISON OF REGULATORY
    DESIGN CONCENTRATIONS
                   AERMOD
                        vs
       ISCST3, CTDMPLUS, ISC-PRIME

                     Staff Report
              U.S. Environmental Protection Agency
            Office of Air Quality Planning and Standards
            Emissions Monitoring and Analysis Division
              Research Triangle Park, North Carolina

-------
                                   Disclaimer
This report has been reviewed by the Office of Air Qualtiy Planning and Standards, U.S.
Environmental Protection Agency, and has been approved for publication. Mention of trade
names or commercial products does not constitute endorsement or recommendation for use.
                                Acknowledgments


The Agency wishes to acknowledge AERMIC (the American Meteorological
Society/Environmental Protection Agency Regulatory Model Improvement Committee),
members of which have given a considerable amount of time, energy and dedication over the last
10 years to develop the AERMOD air dispersion modeling system:

W.D. Peters, U.S.  Environmental Protection Agency, OAQPS, EMAD,AQMG

A. Venkatram, College of Engineering, University of California at Riverside

J. C. Weil, Cooperative Institute for Research in Environmental Sciences,
 University of Colorado

R. B. Wilson,  U.S. Environmental Protection Agency, Region X

R. J. Paine,  ENSR Corporation

S.G. Perry1, Atmospheric Sciences Modeling Division,
Air Resources Laboratory, EPA/ NOAA

R. F. Lee, Consultant, Meteorologist

A. J. Cimorelli, U.S. Environmental Protection Agency, Region III
In addition, Mr.  Roger Erode of MACTEC Federal Programs, Inc. (formerly known as Pacific
Environmental Services Inc) has provided considerable talent and support to the AERMOD
project and has conducted many analyses over the last several years to test and develop this new
air dispersion model.
^n assignment to the Atmospheric Research and Exposure Assessment Laboratory, U.S.
Environmental Protection Agency.

-------
                        TABLE OF CONTENTS

1. INTRODUCTION	4

2. METEOROLOGICAL DATA BASES 	10

3. MODEL OPTIONS AND SOURCE DEFINITIONS	14

4. MODELING RESULTS  	21

5. DISCUSSION OF RESULTS  	  29

6. GENERAL CONCLUSIONS 	36

7. COMPUTER RUN TIMES	37

APPENDIX A
     SIDE BY SIDE COMPARISON OF MODEL FEATURES:
     AERMOD VS ISCST3	38

APPENDIX B
     FIGURES FOR DESCRIBING COMPLEX TERRAIN	44

APPENDIX C
     LOCATION OF SOURCE AND RECEPTORS FOR THE
     COMPLEX TERRAIN ANALYSIS	51

APPENDIX D
      FLAT AND SIMPLE TERRAIN MODELING RESULTS	67

APPENDIX E
     COMPLEX TERRAIN MODELING RESULTS  	80

-------
1. INTRODUCTION

1.1 Background information.
       This report is a final version of an earlier consequence analysis2, which was released to
support the proposal of AERMOD (the American Meteorological Society/Environmental
Protection Agency Regulatory Model Improvement Committee's Dispersion Model, version
99351) in a Federal Register notice3 on April 21, 2000. At that time, the EPA also proposed an
additional model, ISC-PRIME ( Industrial Source Complex -Short Term Model[Version 3] -
Plume Rise Model Enhancements), designed to be used in cases where building downwash was
significant;  AERMOD was to be used for air pollution source scenarios where downwash was not
an issue.  To support the ISC-PRIME proposal, there was a separate but similar building-
downwash-consequence analysis completed which compared ISC-PRIME to ISCST34 (Industrial
Source Complex -Short Term Model—Version 3).  Responding to the overwhelming reaction  from
the commenters on the proposal, the Agency decided to incorporate PRIME algorithms into
AERMOD and thereby eliminate the use of the ISC-PRIME model. The final results in this report
consider both downwash and non-downwash source scenarios since AERMOD now provides the
state of the science for modeling both types of source scenarios. Thus, this  report is designed to
supercede the two earlier consequence analyses.
       This analysis is based on the lastest version of AERMOD, version 022225, which includes
the PRIME algorithms and the proposed version of AERMOD (99351). The ISC-PRIME results
are based on version 99020;  the ISCST3 results are based on version 96113 (for the downwash
analysis) and version 97363 (for the point, area and volume sources) which are the same versions
of the models used in the earlier consequence analyses.
       The introduction includes the following additional sections: a description and purpose of a
consequence analysis;  a description of the 3  components to this study;  and, a brief description of
the air dispersion models of interest - AERMOD (including a list of AERMOD changes since the
proposal), ISCST3, ISC-PRIME,  and, CTDMPLUS (the Complex Terrain Dispersion Model-
Plus)

 1.2 What is a consequence analysis?
       The purpose of this report, often called a consequence analysis, is to give the user
community a sense of how regulatory design concentrations from a new air dispersion model
compare to those from an established model via a series of  "representative" examples. After the
release of a new model for regulatory applications, the user community will want to know:  "What
does this mean to my modeling projects?". This analysis is designed to answer that question by
2Peters,W.D. et al, "Comparison of Regulatory Design Coi
draft document, April 1999, available on the EPA website:

'Federal Register notice, 65FR21506, April 21, 2000.
Concentrations: AERMOD versus ISCST3 and CTDMPlus",
  : www.epa.gov/scramOO 1.
4Paine, R.J. and Lew, F., "Consequence Analysis for Adoption of PRIME: an Advanced Building Downwash Model",
August 24, 1998, available on the EPA website: www.epa.gov/scramOO 1.

5Available on the EPA website: www.epa.gov/scramOO 1.

                                            4

-------
showing the effects of the new model as compared to the existing regulatory model which it
replaces.  For this study, the new model is AERMOD with the PRIME algorithms. The existing
regulatory models used in this report are ISCST3, ISC-PRIME, and CTDMPLUS. This
consequence analysis does not substitute for detailed comparative evaluations or sensitivity
analyses, but rather,  provides to the modeler some simple comparisons of regulatory design
concentration estimates from these air quality models for an extensive number of typical source
scenarios.

 1.3  The three components of this study.
       There are three parts to this study: the flat and simple terrain component; the building
downwash component; and, the complex terrain component. The building downwash component
has been added to the original report since AERMOD now contains the PRIME building
downwash feature and will  be used for sources near buildings. All of the study components use
source scenarios and meteorological data sets which remain unchanged from the earlier
consequence analyses.
       The flat and simple terrain consequence analysis is based on comparative runs made using
a composite of standard data sets. These data sets include a range of point sources with varying
stack parameters,  area and  volume sources, and two point sources in simple terrain6. All source
scenarios are evaluated with two meteorological data sets representing different climatic regimes
in the U.S.   For building downwash, a series of point sources with varying stack heights and
different building configurations are included in the data sets.  Only one of the meteorological
data  sets used in the previous description is used in this part of the analysis.  For the complex
terrain, the study  includes a number of stack heights, buoyancy regimes, distances from source to
hill,  and hill types along with its own meteorological data base (one site).
       After applying the model to all of the above source scenarios, the consequence analysis is
summarized by tabulating the important regulatory (design) concentrations for the new model
against those predicted by the existing regulatory models. Often, the concentrations of regulatory
interest are the high and the high-second-highest concentrations for 1-hour, 3-hour, 24-hour, and
annual averages,  and they are used in this study. The choice of averaging times is based on the
earlier consequence analyses, although this choice is not consistent across all three components of
this study.

1.4 A Brief Description of AERMOD7
       A committee, AERMIC (the American Meteorological Society/Environmental Protection
Agency Regulatory Model  Improvement Committee), was formed to introduce state-of-the-art
modeling concepts into the  EPA's local-scale air quality models.  AERMIC's focus was on  a new
platform for regulatory steady-state plume modeling; this platform would include air dispersion
6 Simple terrain includes receptors with elevations below the top of the stack and at elevations above or below the stack
base. Intermediate terrain includes receptors with elevations above stack top and below the plume centerline.
Complex terrain includes receptors with elevations above the top of the stack.

'User's Guide for the AMS/EPA Regulatory Model - AERMOD, US EPA, Office of Air Quality Planning and
Standards, Research Triangle Park, NC 27711, Report No EPA-454/B-03-001, July  2003. Available on the EPA
website: www.epa.gov/scramOO 1.

-------
fundamentally based on planetary boundary layer turbulence structure, scaling and concepts.
AERMOD is designed to treat both surface and elevated sources in simple and complex terrain.
       Special features of AERMOD include its ability to treat the vertical inhomogeneity of the
planetary boundary layer, special treatment of surface releases, irregularly-shaped area sources, a
three-plume model for the convective boundary layer, and limitation of vertical mixing in the
stable boundary layer. A treatment of dispersion in the presence of intermediate and complex
terrain is used that improves on that treatment currently in use in ISCST3 and other models, yet
without the complexity of a model such as CTDMPLUS.
       AERMOD incorporates, with a new simple approach, current concepts about flow and
dispersion in complex terrain. Where appropriate, the plume is modeled as either impacting and/or
following the terrain. This approach is designed to be physically realistic and simple to implement
while avoiding the need to distinguish among simple, intermediate and complex terrain, as is
required by present regulatory models. As a result, AERMOD removes the need for defining
complex terrain regimes; all terrain is handled in a consistent and continuous  manner that is
simple while still considering the dividing streamline concept in stably-stratified conditions.
       AERMOD is actually a modeling system with three separate components:  AERMOD
(AERMIC Dispersion Model), AERMAP (AERMOD Terrain Preprocessor), and AERMET
(AERMOD Meteorological Preprocessor).
       AERMET is the meteorological preprocessor for AERMOD.  Input data can come from
hourly cloud cover observations, surface  meteorological observations and twice-a-day upper air
soundings.  Output includes surface meteorological observations and parameters and vertical
profiles of several atmosheric parameters.
       AERMAP is a terrain preprocessor designed to simplify and standardize the input of terrain
data for AERMOD.  Input data include receptor terrain elevation data. The terrain data may be in
the form of digital terrain data that is available from the U.S. Geological Survey.  For each
receptor, the output includes a location and height scale, which is an elevation used for the
computation of air flow around hills.
       Additional information about AERMOD can be found in other documents.  The model
evaluation paper8 compares both AERMOD (proposed and current versions), CTDMPLUS,
ISCSTS's and ISC-PRIME's model predictions against measured ambient concentrations.  The
Model Formulation Document9 provides  a detailed explanation of the science behind the model.
8 USEPA, "AERMOD: Latest features and Evaluation Results." Office of Air Quality Planning and Standards,
Research Triangle Park, NC 27711, EPA Report No. EPA-454/R-03-003. July 2003. Available on the EPA website:
www.epa. gov/scramOO 1.


9 USEPA, " AERMOD: Description of Model Formulation (Version 02222) ", Office of Air Quality Planning and
Standards, Research Triangle Park, NC 27711, EPA Report No. EPA-454/R-03-004, October 2002. Available on the
EPA website: www.epa.gov/scram001.

-------
The AERMOD, AERMET and AERMAP  User's Guides7-10-11 inform the user community about
the various options and features of the model and its preprocessors.
1.5  Changes made to AERMOD since the proposal
       A summary of the changes made to the AERMOD in response to comments include the
following:

       * adding the PRIME algorithms to the model (response to public comments);
       * modifying the complex terrain algorithms to make AERMOD less sensitive to the
       selection of the domain of the study area (response to public comments);
       * modifying the urban dispersion for low-level emission sources, such as area sources, to
       produce a more realistic urban dispersion and, as a part of this change, changing the
       minimum layer depth used to calculate the effective dispersion parameters for all
       dispersion settings (scientific formulation correction which was requested by beta testers);
       and making an adjustment to the friction velocity and the Monin-Obukhov length for urban
       stable cases (improved scientific formulation);
       * upgrading AERMOD to include all the newest features that exist in the latest version of
       ISC such as FORTRAN 90 compliance and allocatable arrays, EVENTS processing and
       the TOXICS option (response to public comments).

       In doing the follow-up quality control checking of the model and the source code, the need
for additional changes were identified and the following changes have been made:
       * adding meander to: 1) the stable and unstable urban and 2) the rural unstable dispersion
       settings (only the rural, stable dispersion setting considered meander in the earlier version of
       AERMOD - this change provides a consistent treatment of air dispersion in all dispersion
       settings);
       * making some changes to the basic meander algorithms (improved scientific formulation);
       * making a correction to avoid elevated concentrations for terrain below stack base from the
       virtual image source (response to public comments about spurious results in complex
       terrain); and,
       * repairing miscellaneous coding errors.
A more detailed list of corrections are given in the model evaluation paper8.
1.6  Overview of ISCST312.
10USEPA, "User's Guide for the AERMOD Meteorological Preprocessor (AERMET)", US EPA, Office of Air
Quality Planning and Standards, Research Triangle Park, NC 27711, EPA Report No EPA-454/B-03-002, July 2003.
Available on the EPA website: www.epa.gov/scramOO 1.

"USEPA, "User's Guide for the AERMOD Terrain Preprocessor (AERMAP), US EPA, Office of Air Quality
Planning and Standards, Research Triangle Park, NC 27711, EPA Report No EPA-454/B-03-003, August 2002.
Available on the EPA website: www.epa.gov/scramOO 1.

12USEPA,"User's Guide for the Industrial Source Complex (ISC3) Dispersion Models", Office of Air Quality
Planning and Standards, Research Triangle Park, NC 27711, Report No. EPA-454/B-95-003a, September 1995.
Available on the EPA website: www.epa.gov/scramOO 1.

-------
       ISCST3 is especially designed to support the EPA's regulatory modeling programs. This
model is a steady-state Gaussian dispersion model with a number of options available to the user.
These options include the use of stack-tip downwash, buoyancy-induced dispersion, final plume
rise (except for sources with building downwash), a routine for processing averages when calm
winds occur, and default values for wind profile exponents and for the vertical potential
temperature gradients.  The Short Term model also incorporates COMPLEXl screening model
dispersion algorithms for receptors in complex terrain. The user may  select either rural or urban
dispersion parameters, depending on the characteristics of the source location.  A more detailed
side-by-side explanation and comparison of features between ISCST3  and AERMOD is given in
Appendix A.

1.7 Overview of PRIME.
       PRIME was developed by the Electric Power Research Institute to provide new and
improved plume rise and building downwash algorithms.  The PRIME set of algorithms was
incorporated into ISCST3 and the new model was called ISC-PRIME.  The improved algorithms
provided the following new features:
       . consideration of the location of the stack in relationship to the building;
       . consideration of the streamline deflection over the building;
       . inclusion of plume rise affected by the velocity deficit in the wake or vertical wind speed
shear;
       . a linkage between plume material captured by the near wake and far wake concentrations;
       . elimination of discontinuities at the interface between the  two downwash algorithms;
       . provision  of wind direction effects for squat buildings;
       . elimination of the large concentrations predicted by ISCST3 during light wind speed,
stable conditions that are not supported by observations.

A further, more detailed, description of the model13 and the evaluation results 14are available.

1.8 A Brief Description of CTDMPLUS15
       CTDMPLUS is a  refined Gaussian plume dispersion model designed to estimate hourly
concentrations of plume material from elevated point sources at receptors on or near isolated terrain
features. This model can assess stable and neutral atmospheric conditions as well as daytime,
unstable conditions. Its use of meteorological data and terrain information is different from other
regulatory models in that considerable detail for both types of input data is required and is  supplied
by preprocessors specifically designed for CTDMPLUS.
13  L.L. Schulman, D.G. Strimaitis, IS. Scire, "Development and Evaluation of the PRIME Plume Rise and Building
Downwash Model", Jounrnal of Air and Waste Management Association, 50: 378-390, March 2000.

14  R. J. Paine, F. Lew, "Results of of the Independent Evaluation of ISCST3 and ISC-PRIME", Electric Research
Institute, EPRI TR-2460026, November 1997.  Available at www.epa.gov/scramOO 1.

15 User's Guide to the Complex Terrain Dispersion Model Plus Algorithm for Unstable Situations, US EPA,
Atmospheric Research and Exposure Assessment Laboratory, Research Triangle Park, NC 27711, Report No.
EPA/600/8-89/041, March 1989. Available on the EPA website: www.epa.gov/scram001.

-------
       In modeling stable to neutral conditions, a central feature of CTDMPLUS is its use of a
critical dividing-streamline height to separate the flow in the vicinity of a hill into two separate
layers. Flow in the upper layer has sufficient kinetic energy to pass over the top of the hill, while
the streamlines in the lower layer are constrained to flow in a horizontal plane around the hill. In
modeling unstable or convective conditions, the model relies on a probability density function
(PDF) description of the vertical velocities to estimate the vertical distribution of pollutants.

       Hourly profiles of wind and temperature measurements are used by CTDMPLUS to
compute plume rise, plume penetration, convective scaling parameters.  In stable/neutral conditions,
the profiles of turbulence data are used to compute dispersion parameter values at plume height.

       The model calculates on an hourly basis how the plume trajectory is deformed by each hill.
The computed concentration at each receptor is then derived from the receptor position on the hill
and the resultant plume position and shape.

-------
2. METEOROLOGICAL DATA BASES

2.1 Flat and Simple Terrain
       One year of hourly data for two sites were retrieved and processed. The two sites selected
for this study are Pittsburgh, PA (WBAN [Weather Bureau-Air Force-Navy] station No. 94823),
representative of an urban eastern site; and Oklahoma City, OK (WBAN station No. 13967),
representative of a southwestern plains site. The 1964 data are used at the Pittsburgh site and 1984
data are used at the Oklahoma City site. ISCST3 meteorological data were preprocessed by
PCRAMMET and AERMOD meteorological data were preprocessed by AERMET.

       2.1.1  AERMET Overview.   AERMET '"provides a general purpose meteorological
preprocessor for organizing available meteorological data into a format suitable for use by
AERMOD.  National Weather Service (NWS) hourly surface observations and twice-daily upper
air soundings, plus site-specific data from a meteorological measurement program can be processed
in AERMET. There are three stages to processing the data.  The first stage extracts meteorological
data from archive data files and processes the data through various quality assessment checks. The
second stage merges all data available for 24-hour periods (NWS and site-specific data) and stores
these data together in a single file.  The third stage reads the merged meteorological data and
estimates the necessary parameters for use by AERMOD.  Two files are written for AERMOD: 1) a
file of hourly boundary layer parameter estimates;  and, 2) a file of multiple-level observations
(profiles) of wind speed and direction, temperature, and standard deviation of the fluctuating
horizontal and vertical components of the wind.

       Input data used in this part of the study include: 1) hourly specification of wind speed; 2)
hourly specification of wind direction;  3) hourly ambient temperature; 4) hourly solar radiation16;
5) hourly cloud cover values;  6) a quantification of surface characteristics (surface roughness,
albedo, Bowen ratio);  and 7) twice-daily upper air soundings17. Output includes hourly values for
mixing heights and Monin-Obukhov lengths, surface friction velocity, convective velocity scale,
and profiles of wind speed and direction, temperature and turbulence. Table 2-1 lists the albedo,
Bowen ratio, and  surface roughness that are assumed for this analysis.  Table 2-1 lists only the
rural settings for the meteorological data.  The urban analysis is accomplished by setting the urban
mode and urban source option in AERMOD and using the rural meteorological data for the model
inputs.
16Solar and Meteorological Surface Observation Network 1961-1990, Version 1.0, US Department of Commerce,
National Climatic Data Center, Asheville, NC / US Department of Energy, National Renewable Energy Laboratory,
Golden CO, September 1993.

"Radiosonde Data of North America 1946-1992, Version 1.0, Forecast Systems Laboratory, Boulder, CO and
National Climatic Data Center, Asheville, NC, August 1993.

                                             10

-------
Table 2-1. Albedo, Bowen Ratio, and Surface Roughness length assumed for AERMET
preprocessor.
Site
Pittsburgh
Oklahoma City
Option
Rural
Rural
Albedo
0.25
0.25
Bowen
Ratio
0.75
0.75
Surface roughness
(meters)
0.15
0.15
       2.1.2  PCRAMMET Overview.  The PCRAMMET 18model requires the twice-daily
mixing heights and NWS surface observations.  Prior to being made available, the data were
checked for blank fields (missing data) and filled by accepted procedures. A modification was
made to the data sets by setting the minimum mixing heights to 10 meters.  This change was made
to avoid spuriously high or low concentrations for the short stacks.   Only the meteorological data
used for the ISCST3 analysis was affected.

        For ISCST3, the minimum input data requirements to the PCRAMMET are the twice-daily
mixing heights and hourly surface observations of wind speed, wind direction, dry bulb
temperature, opaque cloud cover and ceiling height.  The operations performed by the
PCRAMMET include:  1) calculation of  hourly values for atmospheric stability from
meteorological surface  observations;  and, 2) interpolation of twice-daily-mixing heights to hourly
values. A brief description of the meteorological data for the two sites is given in Table 2-2.
Table 2-2.     Missing soundings and calm wind conditions by site and year.
Site Year
Pittsburgh 1964
Oklahoma City 1984
Anemometer
height (feet)
20
20
Hours/
year
8784
8784
Missing Soundings
0000 GMT19 1200 GMT
0 0
0 0
Calm wind
conditions
858
181
18 PCRAMMET User's Guide, US EPA, Office of Air Quality Planning and Standards, RTF, NC 27711, EPA-454/B-
96-001, October 1996. Available fromEPA's world-wide-web site atwww.epa.gov/scram001.


19 GMT = Greenwich Mean Time
                                             11

-------
2.2  Building Downwash.
       Only the meteorological data from Pittsburgh (1964), as described in the preceding section,
is used in the building downwash scenarios. No modifications were made to the data because there
are no short stacks, (i.e. less than 20 meters) in this part of the analysis.

2.3  Complex Terrain.
      The meteorological data base used in the complex terrain portion of this study is taken from
a project where site-specific data were collected20.  A 100-m tower, instrumented at 10, 50, and 100
meters and sodar equipment were used to gather the meteorological data. The sodar data was
collected at 50-meter intervals, and the 150 - 400 meter sodar data were used with the tower data to
construct the  meteorological  profiles. The use of sodar turbulence data is limited to vertical
turbulence values only. All of the tower  and sodar levels are used in AERMOD and CTDMPLUS
runs.  Only the 100-m tower data (wind speed and wind direction) are used in ISCST3 runs (see
Figure 2-1 for the 100 meter wind rose).  The atmospheric turbulence and dispersion for ISCST3
are addressed by applying atmospheric stability classifications which are estimated by the solar
radiation/delta-T (SRDT) stability scheme21.
       To confine the differences between CTDMPLUS and AERMOD to differences in the
dispersion algorithms, the METPRO22 output used for CTDMPLUS (including the boundary layer
parameters) is reformatted in a mode compatible with AERMOD meteorological data requirements.
However, the predicted concentrations are not sensitive to these boundary layer values because
profiled meteorological data are available at several levels straddling the stack release heights.
20  The data came from an unnamed source.

21 "An Evaluation of A Solar Radiation/Delta-T Method for EstimatingPasquill-Gifford (P-G) Stability Categories",
EPA-454/R-93-055, October 1993.

22 "User's Guide to the CTDM Meteorological Preprocessor Program", EPA-600/8-88-004, 1988.  Available on the
EPA website: www.epa.gov/scram001.

                                             12

-------
                     Wind Rose for 100-m Tower Winds
W —
                                11.0
      1.8    3.3     5.4   8.f
 WIND  SPEED  CLASS  BOUNDARIES
           (METERS/SECOND)

 NOTES:
 DIAGRAM OF THE FREQUENCY  OF
 OCCURRENCE FOR  EACH WIND DIRECTION.
 WIND  DIRECTION  IS THE DIRECTION
 FROM WHICH THE WIND  IS BLOWING.
 EXAMPLE - WIND IS  BLOWING FROM THE
 NORTH  3.7  PERCENT OF THE TIME.
                                           WINDROSE
2-1.
                               13

-------
3. MODEL OPTIONS AND SOURCE DEFINITIONS

3.1 Modeling Options for Flat, Simple and Complex Terrain.
       The regulatory dispersion model used in this study for the flat and simple terrain is the
ISCST3 model. The model was run in the "regulatory mode" which uses the option settings as
described in Table 3-1. Table 3-1 also shows the parallel settings or options used for AERMOD
setup.
Table 3-1. Model Options Used in Consequence Analysis.
ISCST3
* Use stack tip downwash
* Use buoyancy -induced dispersion
* Do not use gradual plume rise (gradual plume
rise is used in complex terrain)
* Use the calms-processing routines
* Use default wind profile exponents
* Use default vertical potential temperature
gradients
AERMOD
* Use stack tip downwash
* Use buoyancy-induced dispersion (not
an option)
* Use gradual plume rise (not an option)
* Use the calms-processing routines(not an
option)
* Calculate wind profiles (not an option)
* Calculate vertical potential temperature
gradients (not an option)
       The results reported in these 2 components of the study are the high and the highest second-
high concentrations averaged over 1-hr, 3-hr and 24-hr short term averages and the high annual
average.

3.2 Source Characteristics for Flat, Simple Terrain
       Ten source types are processed  for the flat terrain part of this study: seven point sources,
one area source and two volume sources.   Source characteristics for each source type are presented
in Table 3-2.  The very buoyant 35 meter stack source and the 200 meter stack source in Table 3-2
are used in the simple terrain part of this study.  All these sources are evaluated using:  1) both the
rural and urban settings;  and  2) both sets of meteorological data.  Thus, there are 48 scenarios
[(10 flat terrain sources + 2 simple terrain sources) x 2 land use settings (rural, urban) x 2
meteorological sites] and 7 different maximum concentration values for a total of 336 cases .
                                            14

-------
                 Table 3-2. Source characteristics for flat and simple terrain.
Point sources
Stack height
(m)
5
10
20
35
(moderately
buoyant)
35 (very
buoyant)23*
100
200*
X,Y location & base
elevation (m)
0,0,0
0,0,0
0,0,0
0,0,0
0,0,0
0,0,0
0,0,0
Emission rate
(g^1)
100
100
100
100
100
100
100
Exit velocity
(ms-1)
0
0
0
11.7
11.7
18.8
26.5
Stack diameter
(m)
2.4
2.4
2.4
2.4
2.4
4.6
5.6
Temperature
(K)
Ambient
Ambient
Ambient
293
432
416
425
Area source
Area (m2)
1,000,000
Length of side (m)
1000
Emission rate Height of emission
(gs^rn2) release (m)
0.0001 0.0
Volume sources
Emission rate
(gs'1)
100
100
Height of emision
release (m)
10
35
Length of side divided
by 4.3 (m)
14.
14.
Vertical dimension
divided by 4. 3 (m)
16.
16.

23* These sources are also used for the simple terrain part of the consequence analysis.




                                                15

-------
3.3 Source Characteristics for Complex Terrain.
       The complex terrain analysis examines a combination of four hills, two stack heights, two
buoyancies, and two source-hill distances.  The four hills are: 1) Piedmont, a hill near Keyser, WV;
 2) Montour Ridge - Crosswind, near Sunbury, PA;  3) Montour Ridge - Alongwind; and  4)
Cinder Cone Butte, located near Boise, ID.  Except for "Montour Crosswind", the sources are
located to the west of the hill centers, at distances of about 1 kilometer for the "close-in" case and
about 10 kilometers for the "far-out" case (See Appendix B for the figures describing the hills).
For "Montour Crosswind", the sources are located to the north of the east-west oriented ridge.  The
meteorological data base used in this study features a high percentage of winds from the northwest
quadrant (see Figure 2-1).  Therefore, the modeling results reflect a large number of cases of plume
transport from the hypothetical sources to these hills.  The source parameters for the complex
terrain analysis  are provided in Table 3-3.  Although there are 32 possible combinations of
hill/source/source-hill distances ( 4 hills x 2 stack heights x 2 buoyancies x 2 source hill distances),
the plume never significantly impacts the Cinder Cone Butte hill in 4 of the cases and are not
included in the analysis. Thus, the results are reported for a total of 28 complex terrain cases.

3.4 Source Characteristics for Building Downwash.
       A series of hypothetical scenarios involving single point sources and rectangularly shaped
buildings were chosen in an earlier work and these configurations are retained for this study.
ISCST3, ISC-PRIME and AERMOD are applied to each scenario.  The test cases include the
following situations:

    * a stack adjacent to a building structure, and also four building heights away from the northeast
    corner of the building;

    * stack height to building height ratios of 1.0 and 2.0;

    * squat, supersquat, and tall building shapes; and,

    * urban and rural settings.

A no-building set of cases is also used for "control" runs.  Not counting the no-building cases,
there are 20 source/building scenarios and three averaging times to provide a total of 60 cases in
this component of the study. The selection of this set of source configurations and averaging times
matches that of the earlier consequence analysis.

The stack parameters are listed in Table 3.4.

       One year of meteorological data (Pittsburgh, 1964) is employed in this analysis.  The results
for the highest second-highest 3-hour and 24-hour concentrations, as well as the highest annual
concentration are tabulated for each run.   The analysis also includes the model predictions for the
highest 1 hour cavity concentration.
                                             16

-------
Table 3-3. Complex Terrain Source Configurations.
Stack height -
Buoyancy

Low/Low
Low/High
High/Low
High/High
Emission rate
(g/s)

1.0
1.0
1.0
1.0
Stack Height
(m)

30.
30.
150.
150.
Stack Gas
Temperature
(K)

400.
500.
400.
500.
Exit
Velocity
(m/s)

10.
30.
10.
30.

Stack
Diameter
(m)

2.0
6.
2.0
6.
Table 3-4. Source characteristics for building downwash analysis - point sources.
Stack height (m)
35
100
Emission rate
(gs'1)
100
100
Exit velocity
(ms-1)
11.7
18.8
Stack diameter
(m)
2.4
4.6
Temperature
(K)
432
416
3.5  Receptor Configuration for Flat and Simple Terrain

       A gridded polar array of receptors is used in the flat terrain portion of the analysis. For the
point sources, there are  36 radials (beginning at 10 degrees from north and spaced every 10
degrees). The distance of the concentric rings are: 125m, 250m, 400m, 800m, 2000m, 4000m,
8000m, and 16000m. The volume and the area source polar grid is also set up for 10 degree radials
but uses concentric ring distances of 125m, 250m, 400m, 800m, and 2000m.

       A gridded polar array of receptors is used for the point sources in simple terrain settings.
There are 36 radials (beginning at 10 degrees from north and spaced every 10 degrees).  The
distance of the concentric rings were:  800m, 2000m, 4000m, 7000m, and 15000m.  The
elevations for the receptors are plotted (with isopleths) in Figures 3-1 (35 meter stack)  and 3-2
(200 meter stack).
                                           17

-------
 3.6  Receptor Configuration for Building Downwash.

       A cartesian receptor grid extending out to 10 kilometers is used in the building downwash
 analysis.  The receptor density varies, with 50-m spacing for the first 500 meters, 100 m spacing out
 to 1000 m, 200 m spacing out to 2000 m and 1000 m spacing out to 10000 m.  This spacing
 matches that used in the original ISC-PRIME consequence analysis.

 3.7  The Complex Terrain Receptor Locations.

       The Figures in Appendix B show the contours of the hills used in the analysis.  AERMOD,
ISCST3, and CTDMPLUS are run with the full year of data described above for 28 combinations of
sources, and source-hill distances (1 and 10 kilometers).  The CCB andMontour
longwind/crosswind setting includes a total of 140 receptors; the Piedmont Hill setting uses  a total
of 144 receptors; and, the Cinder Cone Butte setting uses 140 receptors. Appendix C contains the
input files used to run AERMOD and provides the location and elevations of all the receptor
locations for all runs. In all cases, each model estimates concentrations on single hills downwind
from the source.
                                            18

-------
                            Figure 3-1
         Elevation in Feet Around the 35 Meter Stack
    15000
    10000
     5000
 CD
    -5000
   -10000
   -15000
                                      §
      -15000  -10000   -5000      0      5000    10000    15000


                             Meters (m)


         Note: 25 ft. contours
Stack base = 797 feet, stack top = 911 feet
                             19

-------
                           Figure 3-2
        Elevation in Feet Around the 200 Meter Stack
    15000
 CO
    10000
     5000
       0
    -5000
   -10000
   -15000
                           r//
                                                   8
      -15000  -10000   -5000     0      5000   10000

                             Meters (m)

         Note: 100 ft. contours
Stack base = 797 feet,  stack top = 1453 feet.
15000
                              20

-------
4. MODELING RESULTS

       The results from the three components of the study are given in this section. The three
components are for the flat and simple terrain, the building downwash and the complex terrain
scenarios. The results compare the new (AERMOD version 02222) versus the old model's (ISCST3
or ISC-PRIME) predicted maximum concentrations.  The relationship between the current version
of AERMOD versus the old model is the focal point of this study.  Generally, the modeling
community is not concerned about the magnitude of the concentration predictions, but is interested
in those situations where the new model predicts higher or lower concentrations than the old model.
Thus, the parameter of choice to present the consequences of the new model is the concentration
ratio. The concentration ratio can be calculated by dividing the current version of AERMOD's
maximum predicted concentration by the old model's corresponding maximum concentration. The
concentration ratio parameter is convenient because a ratio greater than 1 occurs when the new
model predicts maximum concentrations higher than the old model and, conversely, concentration
ratios less than 1 occur when the new model predicts lower maximum concentrations.

       As additional information for those who are interested in the model changes since the
proposal, the results include concentration ratios which are based on the earlier consequence
analyses2.   That is,  the concentration ratios between the current version of the AERMOD versus the
proposed version of AERMOD (version 99351) are provided.  Often, this second set of
concentration ratios is redundant, but they directly  help readers who are familiar with the original
consequence analysis and who want to study the changes to the consequences subsequent to the
Federal Register proposal3.

4.1 The Flat and Simple Terrain Results.

       The results for the flat and simple terrain part of this study are found in Appendix D.  The
high and the highest second high (H2H) concentrations for the 1, 3, 24 hour and annual averaging
times are provided for the ISCST3 model (column  2) and, in parallel, the concentration ratios are
provided for the proposed and the current version of AERMOD (columns 3 and 4). The third
column reproduces information presented in the earlier consequence analysis, that is, the ratio of
maximum concentrations comparing the proposed version of AERMOD to ISCST3.  The last
column of numbers presents the ratio of the air quality concentrations as predicted by the current
versus the proposed  version of AERMOD.   Although redundant, this last column helps the reader
to quickly determine the  changes in the consequence analysis since the April 1999 report.  Each
modeling scenario is defined by a code and the code key is provided at the bottom of each page for
convenience.

       Because of the amount of data and complexity of the tables in Appendix D, a series of tables
are presented  to summarize the statistics of the relationship between the predicted concentrations
from AERMOD and ISCST3.  Tables 4-1, 4-2 and 4-3 have identical  structures.  There are 4
columns of numbers providing a distribution of concentration ratios. The second column provides a
distribution, an average and the maximum and minimum value of the concentration ratios, for the
proposed version of AERMOD versus ISCST3, as  reported in Appendix D. Again, this is a
                                            21

-------
reference point back to the earlier consequence analysis.  The third column shows the new
concentration ratios based on the current version of AERMOD (02222) and ISCST3.  The last
column is redundant but directly supplies information about the magnitude of the changes between
the earlier and the current version of the new air dispersion model since those ratios compare the
new version of AERMOD to the proposed version of AERMOD.  Table 4-1 provides the results for
all the modeling scenarios, while Tables 4-2 and Tables 4-3 break out the results by the rural and
urban settings.

       For example, to further explain the summary tables, refer to the third column in Table 4-1.
According to the number in the second row, there are 5 cases where the current version of
AERMOD predicts a maximum concentration that is a factor of 3 greater than ISCST's prediction.
The third row, column 3 indicates that there are 25 cases where AERMOD's predictions are a factor
of 2 greater than ISCST3's.  The fifth row indicates the total number of cases (336) in this
component of the study.  The fourth row and the sixth row values are the most significant.  These
entries  provide the number of cases where the AERMOD (version 02222) maximum concentrations
are higher than ISCST3's (116) or lower than ISCST3's maximum concentrations (220).  In the
seventh row, there are 46 cases where AERMOD concentrations are less than /^ of the ISCST3 's
maximum concentrations. The average concentration ratio in row 10, the highest ratio in row 11
and the lowest ratio are based on all 336 cases.
Table 4-1. Summary statistics based on the ratio of AERMOD predicted
concentration to ISCST3 predicted concentrations for flat and simple terrain
(see Appendix D) over all averaging times and all source types and both rural
and urban settings.



no. of ratios > 4
no. of ratios > 3
no. of ratios > 2
no. of ratios > 1

Total No.

no. of ratios < 1 .0
no. of ratios < 0.5
no. of ratios <0.33
no. of ratios < 0.25

average
high
low
99351AER/
ISC

1
12
43
157

336

179
44
11
2

1.14
4.25
0.22
02222AER/
ISC

0
5
25
116

336

220
46
14
3

0.96
3.82
0.20
02222AER/
99351AER

0
0
0
82

336

254
24
7
1

0.90
1.73
0.22
                                            22

-------
Table 4-2. Summary statistics based on the ratio of AERMOD  predicted concentration to ISCST3
predicted concentrations for flat and simple terrain (see Appendix D) overall averaging times and all
source types - FOR THE RURAL SETTING ONLY.
RURAL RESULTS
99351AER/ 02222AER/
ISC3 ISC3
number of ratios >4
>3
>2
>1
total
<1.0
<0.5
<0.33
<0.25
max
min
average
0
8
27
91
168
77
16
0
0
3.89
0.35
1.25
0
5
25
85
168
83
6
0
0
3.83
0.41
1.21
02222AER/
99351AER
0
0
0
48
168
120
0
0
0
1.73
0.73
1.00
Table 4-3. Summary statistics based on the ratio of AERMOD  predicted concentration to ISCST3
predicted concentrations for flat and simple terrain (see Appendix D) overall averaging times and all
source types - FOR THE URBAN SETTING ONLY.
URBAN RESULTS
99351AER/ 02222AER/
ISC3 ISC3
number of ratios >4
>3
>2
>1
total
<1.0
<0.5
<0.33
<0.25
max
min
average
1
4
16
66
168
102
28
11
2
4.25
0.22
1.02
0
0
0
31
168
137
40
14
3
1.49
0.20
0.71
02222AER/
99351AER
0
0
0
34
168
134
24
7
1
1.61
0.22
0.80
                                           23

-------
4.2  The Building Downwash Results.

       As mentioned in the introduction, this section was not in the April 1999 AERMOD
consequence analysis since AERMOD was not proposed as the model of choice for building
downwash analyses. Because PRIME has now been incorporated into AERMOD, the AERMOD
consequence analysis now contains comparisons of the building downwash models.  The results,
which  are patterned after the earlier ISC-PRIME consequence analysis4, are given in Table 4-4.

       Table 4-4 has four sets of columns.  The first set of three columns are the scenario
descriptions.  The second set of 4 columns include: the maximum annual concentrations from
ISCST3;   the ISC-PRIME to ISCST3 annual concentration ratios (which were reported in the earlier
ISC-PRIME consequence analysis); the AERMOD (with PRIME) to ISCST3 annual concentration
ratios;  and the AERMOD to ISC-PRIME annual concentration ratios (which are redundant).  The
third set of four columns are for the high-second-high 24 hour concentration rations, using the
column structure as for the annual results.   The fourth set of four columns are for the high-second-
high 3 hour concentration ratios. Table 4-5 presents the summary statistics for the building
downwash analysis, that is, the maximum, minimum and average concentration ratios for each  of the
three averaging times and over all the averaging times.

       Table 4-4 makes note of those source/building scenarios where building downwash is
significant. This criterion is based on cavity concentrations. Many source scenarios do not produce
an estimated cavity concentration (e.g. the 100 meter stack separated from the tall building in a rural
setting);  but, those that do are marked as significant downwash sources (e.g. the 35 meter stack next
to the squat building in a rural setting).  There are two examples where both models generate a
cavity concentration output, but the estimated cavity concentrations are very small (the 100 meter
stack next to the squat building in the urban and rural settings).  Table 4-6 presents the summary
statistics  of the maximum concentration ratios only for those cases where there is significant
building downwash.

       In addition to the downwind concentrations, cavity concentrations are calculated and Table
4-7 presents the results for each source/building scenario. The maximum cavity concentrations from
ISC-PRIME and AERMOD (with PRIME) are given respectively in columns 4 and 5 and the
concentration ratios are seen in column 6.  ISCST3 does not contain an algorithm to estimate the
cavity concentration and could not be included in this table. The summary statistics for the
maximum 1 hour cavity concentration ratios are given in Table 4-8.
                                            24

-------
                                                     Table 4-4.  Building downwash results.

Case Dispersion Stack ISC
(M)
No building (reference) Urban 35 27.1
No building (reference) Rural 35 5.8
No building (reference) Urban 100 3.3
No building (reference) Rural 100 0.4
Squat Building -Stack adjacent to NE of building]
Hb=34; 60x120m * Urban 35 232.6
Hb=34;60x120m* Rural 35 236.6
Hb=50; 60x120m Urban 100 4.1
Hb=50; 60x120m Rural 100 1.5
Squat Building -Stack at distance 4*Hb to NE of building]
Hb=34; 60x120m Urban 35 180.4
Hb=34; 60x120m Rural 35 180.7
Hb=50; 60x120m Urban 100 3.9
Hb=50; 60x120m Rural 100 1.4
Tall Building -Stack adjacent to NE of building]
Hb=34; 30x30m * Urban 35 243.9
Hb=34; 30x30m * Rural 35 242.5
Hb=50; 30x30m Urban 100 3.4
Hb=50; 30x30m Rural 100 0.4
Tall Building -Stack at distance 4*Hb to NE of building]
Hb=34; 30x30m Urban 35 150.8
Hb=34;30x30m Rural 35 181.5
Hb=50; 30x30m Urban 100 3.4
Hb=50; 30x30m Rural 100 0.4
Super Squat Building -Stack adjacent to NE of building]
Hb=34; 180x180m* Urban 35 244.2
Hb=34; 180x180m* Rural 35 243.2

ISCP/
ISC3
1.00
1.00
1.00
1.00

1.23
0.87
2.15
1.61

0.22
0.05
1.09
0.38

1.72
1.32
1.63
1.21

0.26
0.04
1.32
1.06

0.75
0.57
ANNUAL RATIOS
AERMOD/
ISC3
0.65
3.36
0.46
3.76

1.23
1.11
0.77
2.79

0.09
0.12
0.44
1.22

1.35
1.32
0.45
3.66

0.12
0.11
0.45
3.67

0.74
0.69
AERMOD/
ISCP
0.65
3.36
0.46
3.76

ISC

198.8
56.1
22.4
5.1

1.00
1.29
0.36
1.74
1439.8
1574.7
30.0
21.9

0.42
2.46
0.40
3.19
1097.4
1556.5
26.0
17.3

0.79
0.99
0.28
3.03
1508.3
1751.9
23.4
7.5

0.45
2.59
0.34
3.46
1023.5
1556.1
22.7
7.3

0.99
1.23
1508.3
1761.8

ISCP/
ISC3
1.00
1.00
1.00
1.00

1.48
1.03
2.07
1.80

0.25
0.11
1.17
0.32

2.24
1.36
1.60
0.76

0.28
0.05
1.42
0.69

0.90
0.66
24 H2H RATIOS
AERMOD/
ISC3
0.63
2.78
0.53
2.34

1.38
1.16
0.87
1.91

0.11
0.12
0.47
0.71

1.87
1.55
0.51
1.57

0.12
0.10
0.52
1.62

0.86
0.69
AERMOD/
ISCP
0.63
2.78
0.53
2.34

0.93
1.13
0.42
1.06

0.47
1.16
0.40
2.19

0.84
1.14
0.32
2.07

0.43
2.18
0.37
2.34

0.96
1.06
Super Squat Building -Stack at distance 4*Hb to NE of building]
Hb=34; 180x180m Urban 35 226.2
Hb=34; 180x180m Rural 35 240.6
0.17
0.05
0.08
0.10
0.45
2.08
1154.4
1556.5
0.20
0.14
0.14
0.15
0.70
1.04


ISC ISCP/ AERMOD/
ISC3 ISC3
362.4 1.00 0.85
174.7 1.00 2.06
58.6 1.00 0.60
22.6 1.00 1.56


3442.4 0.85 0.79
5662.6 0.39 0.45
62.4 1.93 1.29
59.2 1.62 1.41


3007.7 0.18 0.12
4292.8 0.11 0.10
58.6 1.38 0.65
59.2 0.41 0.64


3442.4 1.29 1.20
5662.6 0.64 0.71
58.6 1.58 0.60
31.3 0.99 1.12


2780.6 0.23 0.11
4206.6 0.07 0.09
58.6 1.45 0.60
31.3 0.76 1.12


3017.8 0.74 0.76
5614.5 0.34 0.38


3007.7 0.16 0.14
4292.8 0.12 0.12
3 H2H RATIOS
AERMOD/
ISCP
0.85
2.06
0.60
1.56

0.93
1.17
0.66
0.87

0.66
0.87
0.47
1.57

0.93
1.11
0.38
1.13

0.47
1.25
0.41
1.48

1.02
1.12

0.85
0.98
  Significant downwash
source
                                                                       25

-------
Table 4-5. Summary of the building downwash analysis.
ALL AVERAGING TIMES
ANNUAL RATIOS


ave
max
min
No cases
ISCP/
ISC3
0.88
2.15
0.04
20
AERMOD/
ISC3
1.03
3.67
0.08

AERMOD/
ISCP
1.38
3.46
0.28

ISCP/
ISC3
0.93
2.24
0.05
20
24 H2H RATIOS
AERMOD/
ISC3
0.82
1.91
0.10

AERMOD/
ISCP
1.06
2.34
0.32

3 H2H RATIOS
ISCP/ AERMOD/
ISC3 ISC3
0.76 0.62
1.93 1.41
0.07 0.09
20
AERMOD/
ISCP
0.92
1.57
0.38

ISCP/
ISC3
0.86
2.24
0.04
60
AERMOD/
ISC3
0.82
3.67
0.08

AERMOD/
ISCP
1.12
3.46
0.28

Table 4-6.  Summary of results for those sources with significant downwash.
ANNUAL RATIOS


ave
max
min
No cases
ISCP/
ISC3
1.08
1.72
0.57
6
AERMOD/
ISC3
1.08
1.35
0.69

AERMOD/
ISCP
1.05
1.29
0.79

ISCP/
ISC3
1.28
2.24
0.66
6
24 H2H RATIOS
AERMOD/
ISC3
1.25
1.87
0.69

AERMOD/
ISCP
1.01
1.14
0.84

3 H2H RATIOS
ISCP/ AERMOD/
ISC3 ISC3
0.71 0.71
1.29 1.20
0.34 0.38
6
AERMOD/
ISCP
1.05
1.17
0.93

ISCP/
ISC3
1.02
2.24
0.34
18
ALL AVE RAGING TIMES
AERMOD/
ISC3
1.01
1.87
0.38

AERMOD/
ISCP
1.03
1.29
0.79

                                                       26

-------
       Table 4-7.  Results of the PRIME cavity max hourly concentrations (ug/m3) in ISC-PRIME
and AERMOD (version 02222).   Summary statistics of the AERMOD to ISC-PRIME ratios are
included.
    Case
    No building
    No building
    No building
    No building
Dispersion  Stack
          (M)
          MAX 1 HR CAVITY CONG
                                  AERMOD/
                 ISCP    AERMOD     ISCP
Urban
Rural
Urban
Rural
35
35
100
100
    Squat Building -Stack adjacent to NE of building
    Hb=34; 60x120m  Urban     35               3202        3180      0.99
    Hb=34; 60x120m  Rural     35               2341        3180      1.36
    Hb=50; 60x120m  Urban     100                10*        0.08*      N/A
    Hb=50; 60x120m  Rural     100                 0*        0.08*      N/A
    Squat Building -Stack at distance 4*Hb to NE of building
    Hb=34; 60x120m  Urban     35
    Hb=34; 60x120m  Rural     35
    Hb=50; 60x120m  Urban     100
    Hb=50; 60x120m  Rural     100
    Tall Building -Stack adjacent to NE of building
    Hb=34;30x30m   Urban     35               5034        4388      0.87
    Hb=34; 30x30m   Rural     35               3963        4388      1.11
    Hb=50;  30x30m   Urban     100
    Hb=50;  30x30m   Rural     100
    Tall Building -Stack at distance 4*Hb to NE of building
    Hb=34; 30x30m   Urban     35
    Hb=34; 30x30m   Rural     35
    Hb=50;  30x30m   Urban     100
    Hb=50;  30x30m   Rural     100
    Super Squat Building -Stack adjacent to NE of building
    Hb=34; 180x180m Urban     35               2524        2498      0.99
    Hb=34; 180x180m Rural     35               2321        3276      1.41
    Super Squat Building -Stack at distance 4*Hb to NE of building
    Hb=34; 180x180m Urban     35
    Hb=34; 180x180m Rural     35	
    ' considered
     insignificant
                                             27

-------
Table 4-8.  Summary of the building downwash 1 hour cavity concentration ratios.


ave
max
min
No cases
AERMOD/
ISCP
1.12
1.41
0.87
6
4.3  The Complex Terrain Results.

       The complex terrain results which compare AERMOD to ISCST324 and to
CTDMPLUS, are presented in Appendix E.  The Appendix E table includes the highest
second high ratios along with the highest annual concentrations.  A summary of the complex
terrain results are provided below in Table 4-9. Table 4-9 provides a distribution of
AERMOD to other model concentration ratios.  There were no cases where AERMOD
predicted higher concentrations than either ISC or CTDMPLUS.   As seen before, the last
column in Table 4-9 provides data, presented in the earlier consequence analysis, as a
convenient reference for the reader.
Table 4-9. Summary statistics based on the ratio of AERMOD predicted regulatory design
concentration to ISCST3(COMPLEX1) and CTDMPIus predicted-regulatory-design concentrations
for complex terrain (see Appendix E).


NO OF RATIOS>4
NO OF RATIOS>3
NO OF RATIOS>2
NO OF RATIOS>1
TOTAL NO
NO OF RATIOS<1 .0
NO OF RATIOS<0.5
NO OF RATIOS<0.33
NO OF RATIOS<0.25
AVERAGE
MAX
MIN
02222AER/
ISC3
0
0
0
0
196
196
189
163
119
0.24
0.79
0.07
02222AER/
CTDM+
0
0
2
40
196
152
58
28
15
0.75
2.13
0.14
02222AER/
99351AER
0
0
1
74
196
77
0
0
0
1.01
2.67
0.61
24In complex terrain, the COMPLEXl portion of the ISCST3 model is playing a significant role. For receptors with
elevations above the plume height, the COMPLEXl concentration estimates are used. In intermediate terrain, the
highest estimates from the "ISCST3" simple terrain model and the COMPLEXl model are used.
                                          28

-------
5. DISCUSSION OF RESULTS

5.1 Discussion of Flat and Simple Terrain Results

5.1.1.  The current version of AERMOD (02222).  With all the flat and simple terrain
results viewed as a whole, the current version of AERMOD produces maximum
concentrations that are similar to ISCSTS's predicted concentrations.   The reported average of
0.96 ( Table 4-1,  row 10 column 3) indicates that AERMOD predicts concentrations only
about 4% lower than ISCST3. This average is taken over all source types, stack heights,
settings and concentration averaging times.  However, as expected when a new model is
developed, there are differences between the old and the new air dispersion model predictions.
Although about 80% of the AERMOD concentrations are within a factor of 2 (high or low)
from the ISCST3 concentrations Table 4-1, column 3 also indicates that, for certain situations,
the proposed AERMOD predictions are higher or lower than ISCST3 predictions by  a factor
of 3 or more.  Upon studying column 4 in the detailed Appendix D listing, one can find that
the most significant differences between the 2 models are found in the following scenarios:

1. in the rural, low level stacks (AERMOD is lower);

2. in the long-term concentrations for the rural, taller stacks (AERMOD is higher);

3. in the short-term, urban short stacks and urban area sources (AERMOD is lower);   and,

4. in all of the regulatory concentrations for urban very tall stacks in simple terrain
(AERMOD is lower).

These results are  consistent across the 2 meteorological data bases.  Because the 2 air
dispersion models are significantly different from one another (see Appendix A for side-by-
side comparison), such variation in the model differences is expected.

5.1.2.  Impacts of the changes to AERMOD.    Changes made to AERMOD produced changes
to the consequence analysis.  Column 2 in Table 4-1 shows that, over all source types, settings, and
averaging times, the proposed  version of AERMOD predicts an average concentration ratio of 1.14
or about 14% higher than ISCST3 (compared to 0.96  or about 4% lower based on the current
version of AERMOD).  Thus, there is about a 16% 25reduction in the overall average of
concentration ratios.  Table 4-1 column 3 also indicates that the differences in the AERMOD
predictions are not as extreme as seen in the earlier consequence analysis.  The number of cases
where the AERMOD to ISCST3 ratios are greater than a factor of 2 drops from 43 out of a total of
336 (proposed version) to 25 (current version). The number of cases where the concentration ratios
are greater than 3 drops from 12 to 5. On the other side of the distribution where AERMOD
predictions are less than ISCST3's predictions, the distribution of cases where the AERMOD to
  The 16% = (1.14-0.96) /1.14 x 100.

                                            29

-------
ISCST3 concentration ratios are less than one-half, less than one-third and less than one-fourth
remains almost the same as those from the April 1999 analysis.

       Reviewing the urban and rural breakout tables (Tables 4-2, 4-3), one sees that the largest
changes occur in the urban setting while the rural distribution and averages do not change
significantly.  In the urban setting results in Appendix D, the most significant changes are found in
the following scenarios:

1.  in the short-term concentrations for area sources (a decrease in concentrations due to the proper
urbanization of the dispersions parameters);

2.  in the short-term concentrations for low stacks (a decrease in concentrations due to the addition
of meander for both stable and unstable settings); and,

3.  in the overall urban category when comparing the averages in Table 4-3, row 11, columns 3 and
4 (a decrease due to the addition of stable and unstable urban meander).

In the rural setting, less significant changes are seen and those changes  are found:

1.  in the short-term concentration for low stacks (an increase due to the changing of the minimum
layer depth used to calculate the effective dispersion parameters- a  secondary effect from fixing the
urban dispersion problem); and,

2.  in the overall results (a slight overall reduction in concentration  predictions due to the addition
of meander in the rural, stable setting).

5.1.3.  Model Evaluation Study support.   Differences between models leads to the next topic for
discussion - how do these models perform when compared to measured data?  Do the differences
represent an improvement in model predictions? The model evaluation study (MES)8 provides
AERMOD's (including the proposed version and the current version) and ISCST3's predictions and
compares them to ambient air quality data.   Of the 5  available flat terrain data bases in the MES,
there is one MES site with a low-level release in a flat rural setting. In this scenario, the current
version of AERMOD's (version 02222) short-term concentration predictions (the Robust Highest
Concentrations26 [RHCs]) are about /^ of the ISCST3 estimates, with the AERMOD estimates more
closely matching the observed values.  The MES also includes 3 rural, tall stack scenarios in flat or
simple terrain locations which shows AERMOD predicting long-term concentrations (RHCs that
are almost twice as high as the ISCST3 predictions.  In all 3 cases, AERMOD predictions are
closer to the measured values.  The MES does not have any data bases which are representative of
the shorter stacks in urban settings.  The one urban data set (tall stack) in the MES is based on a
limited monitoring study.  The urban one-hour RHC for AERMOD is about 20% lower than the
ISCST3 with AERMOD predictions closer to the measured concentrations.  The MES supports the
26 W.M. Cox, J.A. Tikvart, "A Statistical Procedure for Determining the Best Performing Air Quality Simulation
Model", Atmospheric Environment, Vol.24A, No 9, pp 2387-2395, 1990

                                             30

-------
AERMOD concentration predictions over those provided by the ISCST3 model, i.e., AERMOD's
performance is better than ISCSTS's performance when compared to monitored concentrations.

       Also, the MES indicates a slight performance improvement when comparing the current
version of AERMOD (version 02222) to the proposed version of AERMOD.  The most notable
differences are in the short term concentrations with a tall stack in the urban area and with a tall
stack in moderate hilly terrain in a rural area.  In both settings, the current version of AERMOD
predicts lower RHCs than the proposed version (which is consistent with this study) and predicts
concentrations that are closer to the measured values.

5.2 Discussion of Building Downwash Results.

      The discussion in this section does not include the proposed version of AERMOD.  The
proposed version of AERMOD (99351) does not include the PRIME algorithms and does not play a
role in the building downwash analysis.  For this component of the consequence analysis,  the
current version of AERMOD (02222) is the new model,  ISC-PRIME is the proposed model, and
ISCST3 is the currently approved model.

5.2.1 AERMOD versus ISCST3.   Table 4-4 presents the results of the downwash analysis and
Table 4-5 presents a summary of the results.  The summary table indicates that AERMOD (with
PRIME) produces somewhat lower maximum concentration estimates, on average,  than ISCST3.
The AERMOD to ISCST3 ratios in column 12 displays an average concentration ratio of 0.82;
that is, on average over all 60 cases (10 source types, 2 settings and 3 averaging times),  AERMOD
predicts concentrations that are about 18% lower than ISCST3's predictions.  However, the range
of the concentration ratios is more significant; AERMOD predicted concentrations that are up to a
factor of 4 higher and up to a factor of 10 lower than ISCST3.

       The situation where AERMOD's predicts maximum concentrations much lower than the
ISCST concentrations is found in several cases.  For example,  the 35 meter stack separated from
the building for all averaging times in urban and rural settings for all 3 building types (total of 18
cases in Table 4-4) indicate lower AERMOD predictions.  Because the downwash algorithms in
ISCST3 ignore the separation between stacks and buildings and is designed to be environmentally
conservative27, relatively smaller AERMOD concentrations are expected for the shorter stacks.
This scenario with stack/building separation was chosen originally to highlight the differences in
the ISCST3 and PRIME models, thus, significant differences are expected.   The vast majority of
the remaining cases have AERMOD concentrations that are within a factor of 2 of the ISCST3
concentrations.

       There are situations where AERMOD is higher than ISCST3. For example, in the annual
concentration estimates for the rural setting, tall and squat building with a  100 meter stack (both
adjacent to the building and separated from the building), the AERMOD maximum concentrations
are larger than the ISCST3 estimates.  In these cases, building downwash is not important to the
27ISCST3 assumes that the stack is located in the center of the building which maximizes the impacts of the building
wake on the plume dispersion. Generally, this assumption will predict the highest concentrations.
                                           31

-------
calculation of maximum annual concentration estimates, so the difference between AERMOD and
ISCST3 concentration estimates are attributable only to the differences in the dispersion algorithms
within the two models.  The insignificance of building downwash for these cases is seen:  1) by the
lack of a calculated cavity concentration (Table 4-7); and, by comparing the ISCST3 maximum
concentrations to the matching non-downwash case for 100 meter stack in the rural setting (Table
4-4).  The maximum annual concentrations and the concentration ratios for the downwash cases
remains basically unchanged from the corresponding "No building" case.

       There is a second set of statistics prepared for those sources where cavity concentrations are
calculated, that is, building downwash is a known significant factor in the dispersion of the plume
(Table 4-6).  These sources are marked in Table 4-4.  The summary results in column 12, Table 4-
6, show that the AERMOD maximum concentration predictions are, on average, about the same as
those produced by ISCST3 (average concentration ratio of 1.01).  The concentration ratios range
from a maximum of 1.87 to a minimum of 0.38.  In all these significant downwash cases in Table 4-
4, the stack is close to the building and the discrepancies between the 2 models are due only to the
differences in the downwash algorithms.

5.2.2  AERMOD versus ISC-PRIME.   Table 4-5 provides the summary information about the
comparison between AERMOD and ISC-PRIME.  Because PRIME is in  both models, the
expectation is that the 2 models should be in reasonable agreement.  On average over the 60 cases,
the ratio of AERMOD to ISC-PRIME maximum concentration predictions was 1.12 (AERMOD
predictions are about  12% higher than ISC-PRIME predictions), which is  rather good agreement.
However,  the maximum concentration ratio (3.46) and the minimum concentration ratio ( 0.28) are
of initial concern.  When studying the concentration values in Table 4-4  for those cases with the
highest differences (rural case with 100 meter stack separated from a tall building) and lowest
differences (urban case with 100 meter stack near a tall building), one can see that there  are similar
differences between the two models in the no-building scenario (rows 3 and 4).  Thus, these
extreme cases, which are not significant downwash cases, are mostly explained by the differences
in the dispersion algorithms.

       For those cases where building downwash was more important, the two models should be in
closer agreement, because the PRIME algorithms are in both models and should be dominating the
dispersion calculations.   Although there are differences in the way that PRIME interacts with the
two dispersion models (the numerical plume rise, the plume capture criteria and blending of the
disturbed plume with the surrounding undisturbed atmosphere), there are a number of tests to check
AERMOD with the PRIME insertion.

       The first test involves the cavity concentrations which should be similar between the two
models. Tables 4-7 and 4-8 indicate that  this is so.  In Table 4-8, for the 6  cases with significant
maximum 1 hour cavity concentrations estimates, the average AERMOD/ISC-PRIME
concentration ratio is  1.12 with the concentration ratio ranging from a maximum of 1.41 to a
minimum  of 0.87.  These cavity concentration differences are attributable to the differences in the
plume rise equations.  There are two cases where a very small cavity concentration is calculated by
                                           32

-------
both models, but these results not included in the previously mentioned statistics.  These two cases
include the 100 meter stack with a squat building in the rural and urban setting.

      The second related test reviews the AERMOD versus ISC-PRIME maximum concentration
results for all averaging times for those  cases with significant downwash (Table 4-6).  The average
AERMOD/ISC-PRIME concentration ratio is 1.03 with a maximum value of 1.29 and a minimum
value of 0.79 (column  13).  As expected in these cases with a definite cavity, the two models agree
very closely (plus or minus 30%).

      The third test of model consistently compares those cases where building downwash is not
significant enough to produce concentrations different from the reference ("No-building") case (the
first 4 rows in Table 4-4). In this test, both models should predict maximum concentrations that are
essentially the  same as those predicted by the model in the corresponding "No building" scenario.
Table 4-4  indicates there are three source  scenarios (all averaging times) where the building has
little or no effects on the ISC-PRIME maximum concentrations: rural, urban 100 meter stacks
separated from a squat building;  and rural, urban 100 meter stacks near and separated from a tall
building.  In all these cases, AERMOD  produces a matching  result in that the building had almost
no effect on the estimated maximum concentrations.

      Conversely, there are cases where AERMOD shows no change from the no-building to the
with building scenario and ISC-PRIME modeling results predict some impact from the building.
Examples  of this result are seen in the 3 hour averaging columns: the urban and rural 100 meter
stack near and  separated from a tall building; the urban and rural 35  meter stack separated from the
tall building; and, the urban 100 meter  stack separated from the tall building.  This difference in
the two models is expected because of the critical angle of plume rise used for calculating the
amount of pollution that is caught in the wake of the building.  AERMOD development work
suggested  a change in this area and a different critical angle was implemented in the AERMOD.

5.2.3. Model  Evaluation Study support.   The original model evaluation reports13>14, in general,
support  the addition of PRIME to the ISCST3 model.  The reports conclude that the PRIME had a
statistically better performance result for each data base in the independent evaluation.  Although
these results support the implementation of PRIME into the regulatory models,  none of the
evaluation databases have examples of stacks significantly separated from the building.   So the
model differences cannot be confirmed  for this scenario.

      Although the AERMOD and ISC-PRIME maximum concentration estimates are in general
agreement, there are differences.  Some variations are expected because of the way that PRIME is
integrated  into AERMOD.  When reviewing the model evaluation results8, one finds that the
AERMOD performance is slightly better than ISC-PRIME.   There are four cases of slight
degradation of performance, four cases  of similar performance and, five cases with improved
performance.   Of the five cases with improved performance, there is one case with a rather
dramatic performance improvement.
                                           33

-------
5.3 Discussion of Complex Terrain Results.

5.3.1  The current version of AERMOD (02222).   The current version of AERMOD (02222)
consistently produced lower or significantly lower regulatory design concentrations estimates than
those generated from the ISCST3 (Table 4-9).  This result was expected since that portion of the
ISCST3 model that deals with complex terrain,COMPLEXl, is used for screening purposes and
has been designed to be conservative28.  The average AERMOD/ISCST3 concentration ratio,
overall scenarios and averaging times, is 0.24 (column 2, row 10). Thus, AERMOD produced
maximum concentrations that were, over all all cases, a factor of 4 lower than the ISCST3
predictions.  The concentration ratios ranged from 0.07 to 0.79.  In well over half of the cases (119
of the 196 cases),  AERMOD produced maximum concentrations that were a factor of 4 lower than
the ISCST3 estimates.  There were no cases where AERMOD predicted maximum concentrations
higher than ISCST3.  When examining Appendix E, the differences between the 2 models tended
to increase with averaging time, i.e. the largest differences were seen in the maximum annual
averages.

       Also, Table 4-9 and Appendix E contain the results for the AERMOD/CTDMPLUS
comparison.  As expected, the models agreed more closely, since CTDMPLUS, which is not a
screening model like COMPLEX1, is a more refined, site-specific complex terrain model.  Table 4-
9 indicates that the average AERMOD/CTDMPLUS ratio over all cases was 0.75 with a range of
0.14 (lowest value) to 2.13 (highest value).   Table 4-9 indicates that AERMOD predicted
maximum concentrations that were larger than CTDMPLUS in about 20% of cases (40 out of 196).

5.3.2.  Impacts of the changes to AERMOD.  In the complex terrain scenarios, changes to the
complex terrain algorithms do not produce much of a change to the consequence analysis.  Some
minor changes were made to the algorithm to make the model's concentration predictions less
sensitive to the domain selection. The critical value of hill height scale was not changed in the
input files. Thus, significant changes from the earlier consequence analysis are not expected.
Table 4-9  column 4 summarizes the numerical changes to the concentration ratios of the current
version of AERMOD to the proposed version of AERMOD.   The average concentration ratio over
all 196 cases is 1.01 which implies that the changes did not bias the model towards higher or lower
concentration estimates. The range of concentration ratios of new versus the old version of
AERMOD was 0.61 to 2.67.  The one high value of 2.67 was due to changes in the meander
algorithm  and was not a repercussion of the complex terrain changes.  Although not shown in the
summary table, about 94 % of the concentration ratios were within a factor of plus or minus 30%.

5.3.3.  Model evaluation study support. As mentioned above, the differences between
AERMOD and ISCST3 tend to increase with averaging time, i.e. the largest differences are seen in
the maximum annual averages. These results are confirmed by the complex terrain model
evaluation results. AERMOD consistently predicts lower maximum concentrations than ISCST3
28 That is, the model is designed to overestimate or produce concentration estimates which are larger than
concentrations that one would measure.

                                           34

-------
with the largest variations occurring in the annual averages. AERMOD provides much better
performance in the complex terrain data bases.

       AERMOD does not consistently predict lower maximum concentrations than
CTDMPLUS's estimates, as is seen in the comparisons.  There are cases where AERMOD is higher
and lower than the site-specific complex terrain model.   These results are consistent with the model
evaluation results as AERMOD produces RHCs higher and lower than  CTDMPLUS.  However, in
all but one of the 10 evaluation cases, AERMOD outperforms  CTDMPLUS.

       The model evaluation results indicate that the current version of the model performs slightly
better than the proposed version.  In all of the 10 complex terrain database cases, the change to the
performance is no greater than 9%.
                                           35

-------
6. GENERAL CONCLUSIONS

The following general conclusions are made.

1) For non-downwash settings in flat and simple terrain, the current version of AERMOD (version
02222): a) on average, tends to predict maximum concentrations that are similar to ISCST3; b) ,
on average, tends to predict concentrations closer to ISCST3 than the proposed version of
AERMOD; and, c) predicts maximum concentrations which are not as extreme in their differences
from ISCST3 as those seen when  applying the proposed version of AERMOD; and, on average,
tends to predict urban maximum concentrations that are lower than the proposed version of
AERMOD.

2) Where building downwash is a significant factor in the air dispersion analysis, the current
version of AERMOD predicts maximum concentrations and maximum cavity concentrations that
are very similar to ISC-PRIME.

3) In general, the consequences from using the current version of AERMOD instead of ISCST3 in
complex terrain are significant, the current version of AERMOD produces much lower maximum
concentrations than the screening technique in ISCST3.  Also, the  current version of AERMOD
produced results that are essentially unchanged from the results reported using the proposed version
of AERMOD.  When compared to CTDMPLUS, AERMOD tends to predict somewhat lower
maximum concentrations with examples of AERMOD predictions being higher and lower than the
CTDMPLUS predictions.

4) Where data are available, the model evaluation results support the differences identified in this
report when comparing the proposed version of AERMOD to ISCST3 and when comparing the
current version of AERMOD to the proposed version of AERMOD.  . The model evaluation report
indicates that the current version of AERMOD outperforms all the  other four models ( ISCST3,
ISC-PRIME, CDTMPLUS and the proposed version of AERMOD).

5) Because of the stability  of AERMOD model throughout the consequence analysis and because the
model evaluation study supports AERMOD (02222) when significant differences occur between the
current version of AERMOD and ISCST3 or the earlier version of AERMOD, it is appropriate for
the Agency to adopt the current version of AERMOD (02222) as a regulatory model and is a
suitable replacement for ISCST3  for many regulatory applications.
                                           36

-------
7. COMPUTER RUN TIMES

       As an additional feature to the consequence analysis, ISCST3 and AERMOD models were
compiled and run on a typical personal computer.  The purpose of this exercise was to provide the
user community with a sense of the potential changes in the amount of time to run typical source
configurations on their computer systems. The results are tabulated below in Table 7-1.   The
computer used to complete this table was a Pentium 2.4  Gigahertz computer with 256 megabytes  of
random access memory. Each source run evaluated 180 receptors in a polar grid (36 radials with  5
ring distances).  One full year (8784 hours) of meteorological data was used for each run.
Table 7-1. Computer Run Times for Typical Source Configurations.
Source Type
Point source
Point source w/ downwash
Volume source
Area source
AERMOD
20.1 seconds
179
12.8
2190
ISCST3
2.34 seconds
42
2.03
515
AERTSC
RATIO
8.5
4.1
6.3
4.3
                                           37

-------
       APPENDIX A





SIDE BY  SIDE COMPARISON





 AERMOD VERSUS ISCST3
            38

-------
Feature
ISCST3
AERMOD (version 02222)
Comments
Types of sources
modeled
Point, area, and
volume sources
Same as ISCST3
Models are comparable
Plume Rise
Uses Briggs equations
with stack-top wind
speed and vertical
temperature gradient
In stable conditions, uses Briggs
equations with winds and
temperature gradient at stack top
and half-way to final plume rise;
in convective conditions, plume
rise is superposed on the
displacements by random convective
velocities
AERMOD is better because in
stable conditions it factors in
wind and temperature changes
above stack top, and in
unstable conditions it accounts
for convective updrafts and
downdrafts
Meteorological Data
Input
One level of data
accepted
An arbitrarily large number of data
levels can be accommodated
AERMOD can adapt multiple
levels of data to various stack
and plume height	
Profiling
Meteorological Data
Only wind speed is
profiled
AERMOD creates profiles of wind,
temperature, and turbulence, using
all available measurement levels
AERMOD is much improved over
ISCST3 in this area
Use of
Meteorological Data
in Plume Dispersion
Stack-top variables
for all downwind
distances
Variables measured throughout the
plume depth (averaged from plume
centerline to 2.15 sigma-z below
centerline; changes with downwind
distance)
AERMOD treatment is far more
advanced than that of ISCST3;
accounts for meteorological
data throughout the plume depth
Plume Dispersion:
General Treatment
Gaussian treatment in
horizontal and
vertical
Gaussian treatment in horizontal
and in vertical for stable
conditions; non-Gaussian
probability density function in
vertical for unstable conditions
AERMOD's unstable treatment of
vertical dispersion is a more
accurate portrayal of actual
conditions	
                                                         39

-------
Feature
ISCST3
AERMOD (version 02222)
Comments
Urban Treatment
Urban option either
on or off; no other
specification
available; all
sources must be
modeled either rural
or urban
Population is specified, so
treatment can consider a variety of
urban conditions; sources can
individually be modeled rural or
urban
AERMOD provides variable urban
treatment as a function of city
population, and can selectively
model sources as rural or urban
Characterization of
Modeling Domain
Surface
Characteristics
Choice of rural or
urban
Selection by direction and month of
roughness length, albedo, and Bowen
ratio, providing user flexibility
to vary surface characteristics	
AERMOD provides the user with
considerably more options in
the selection of the surface
characteristics
Boundary Layer
Parameters
Wind speed, mixing
height, and stability
class
Friction velocity, Monin-Obukhov
length, convective velocity scale,
mechanical and convective mixing
height, sensible heat flux
AERMOD provides parameters
required for use with up-to-
date  planetary boundary layer
(PEL) parameterizations; ISCST3
does not
Mixed Layer Height
Holzworth scheme;
uses interpolation
based upon maximum
afternoon mixing
height
Has convective and mechanical mixed
layer height; convective height
based upon hourly accumulation of
sensible heat flux
AERMOD's formulation is
significantly more advanced
than that of ISCST3, includes a
mechanical component, and in
using hourly input data,
provides a more realistic
sequence of the diurnal mixing
height changes	
                                                         40

-------
Feature
ISCST3
AERMOD (version 02222)
Comments
Terrain Depiction
Elevation at each
receptor point
                                            Controlling hill  elevation  and
                                            point  elevation at  each  receptor,
                                            obtained  from  special  terrain pre-
                                            processor (AERMAP)  that  uses
                                            digital elevation model  (DEM) data
                                      AERMOD's terrain pre-processor
                                      provides information for
                                      advanced critical dividing
                                      streamline height algorithms
                                      and uses digital data to obtain
                                      receptor elevations
Plume Dispersion:
Plume Growth Rates
Based upon 6 discrete
stability classes
only; dispersion
curves (Pasquill-
Gifford)  are based
upon surface release
experiments (e.g..
Prairie Grass)
Uses profiles of vertical and
horizontal turbulence (from
measurements and/or PEL theory);
variable with height; uses
continuous growth functions rather
than a discrete (stability-based)
formulation
Use of turbulence-based plume
growth with height dependence
rather than that based upon
stability class provides AERMOD
with a substantial advancement
over the ISCST3 treatment
Plume Interaction
with Mixing Lid:
convective
conditions
If plume centerline
is above lid, a zero
ground-level
concentration is
assumed
Three plume components are
considered:  a "direct" plume that
is advected to the ground in a
downdraft, an "indirect" plume
caught in an updraft that reaches
the lid and eventually is brought
to the ground, and a plume that
penetrates the mixing lid and
disperses more slowly in the stable
layer aloft (and which can re-enter
The AERMOD treatment avoids
potential underpredictions
suffered by ISCST3 due to its
"all or nothing" treatment of
the plume; AERMOD's use of
convective updrafts and
downdrafts in a probability
density function approach is a
significant advancement over
ISCST3
                                                         41

-------
Feature

Plume Interaction
with Mixing Lid:
stable conditions
Building Downwash
ISCST3

The mixing lid is
ignored (assumed to
be infinitely high)
Combination of Huber-
Snyder and Scire-
Schulman algorithms;
many discontinuities
AERMOD (version 02222)
the mixed layer and disperse to the
ground)
A mechanically mixed layer near the
ground is considered. Plume
reflection from an elevated lid is
considered.
New PRIME downwash algorithm
installed
Comments

AERMOD' s use of a mechanically
mixed layer is an advancement
over the very simplistic ISCST3
approach
AERMOD benefits from the
technological advances offered
by the PRIME model
This page left intentionally blank.
                 42

-------
           APPENDIX B
FIGURES FOR DESCRIBING COMPLEX
            TERRAIN

-------
lished Qlipses.   Prgss  ENTER IPIEPHNTB - unedited |
      t  Kn.
                              Figure B-l. Plot of Piedmont Hill Contours.
                                       44

-------
                              TERRAIN INFORMATION
HILL NAMED PIEDMONT                 HILL TOP: 2240.0 feet
Figure B-l displays both actual and transformed contours for the Piedmont hill.  The
transformed contours were used in the modeling analysis. There are 2 sets of 13 contours;
the dotted lines are the actual contours and the solid lines are the transformed or modeled
contours.  The elevations  for the 13 contours are listed below:

  1000.0 feet
  1100.0
  1200.0
  1300.0
  1400.0
  1500.0
  1600.0
  1700.0
  1800.0
  1900.0
  2000.0
  2100.0

The actual receptor locations are listed in the AERMOD control files which are listed in
Appendix C.
                                          45

-------
Finished  QlipsQs.
                    CNTEP  IHONTDUR - Unedited)
1 KM
                        Figure B-2. Plot of Montour Ridge Contours.
                                  46

-------
                              TERRAIN INFORMATION
 HILL NAMED MONTOUR          HILLTOP:  1429.0  meters

Figure B-2 displays both actual and transformed contours for the Montour Ridge. The
transformed contours were used in the modeling analysis. There are 2 sets of 8 contours;  the
dotted lines are the actual contours and the solid lines are the transformed or modeled
contours.  The elevations for the 8 contours are listed below:
  700.0 meters
  800.0
  900.0
  1000.0
  1100.0
  1200.0
  1300.0
  1400.0
The actual receptor locations are listed in the AERMOD control files which are listed in
Appendix C.
                                         47

-------
u
Ol
W
D.


ol
                                                                                                            o
                                                                                                           U
                                                                                                           33
                                                                                                            o
                                                                                                           U
                                                                                                           •B
                                                                                                            o
                                                                                                           PH

                                                                                                           fi

                                                                                                           (C

-------
TERRAIN INFORMATION
 HILL NAMED CINDER CONE BUTTE (CCB)           HILL TOP:  100.0 meters
Figure B-3 displays both actual and transformed contours for  CCB.  The transformed
contours were used in the modeling analysis. There are 2 sets of 11 contours; the dotted lines
are the actual contours and the solid lines are the transformed or modeled contours.  The
elevations for the 11 contours are listed below:
   0.0 meters
   5.0
   10.0
   20.0
   30.0
   40.0
   50.0
   60.0
   70.0
   80.0
   90.0
The actual receptor locations are listed in the AERMOD control files which are listed in
Appendix C.

-------
   APPENDIX C
LOCATION OF SOURCE AND RECEPTORS
FOR THE COMPLEX TERRAIN ANALYSIS
                  50

-------
                               ** Cinder Cone Butte
                           ** Low Stack Height: 30 meters
                ** High Buoyancy Case: 6-m diameter, 30 m/s Exit Vel, 500 K
                             ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar

** Cinder Cone  Butte
** Low Stack Height:  30 meters
** High Buoyancy Case:  6-m diameter,  30 m/s Exit Vel, 500 K
** Close to Hill: 1 km  away
** Meteorology  from 100-m Tower  and  Sodar
**
CO STARTING
   TITLEONE  Cinder Cone Butte:  Source 1 km away
   TITLETWO  30-m Stack Height;  High Buoyancy Case
   MODELOPT  CONG             MSGPRO
   AVERTIME  1  3 24 Period
   POLLUTID  SO2
   RUNORNOT  RUN
   ERRORFIL  ERRORS.OUT
   TERRHGTS  ELEV
CO FINISHED
SO STARTING
   ELEVUNIT  FEET
   LOCATION  STACK1   POINT
** Point Source       QS
** Parameters:        	
   SRCPARAM  STACK1     1.
0.0  1500.0   3100.0
HS    TS    VS   DS

30.  500.  30.0  6.0
   SRCGROUP  ALL
SO FINISHED
RE
RE
**
RE
RE
RE
RE
RE
RE
RE
STARTING
ELEVUNIT
FEET
X (meters) y (meters)
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
2
-88
-164
-222
-258
-288
-287
.00
.98
.25
.04
.22
.26
.25
-302.
-269.
-208.
-130.
-40.
51.
147.
00
90
74
75
01
90
62
z (feet)
3198
3198
3198
3198
3198
3198
3198
.42
.42
.42
.42
.42
.42
.42
3395
3395
3395
3395
3395
3395
3395
.27
.27
.27
.27
.27
.27
.27
                                        51

-------
               ** Cinder Cone Butte
           ** Low Stack Height: 30 meters
High Buoyancy Case: 6-m diameter, 30 m/s Exit Vel, 500 K
             ** Close to Hill: 1 km away
     ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
-241.
-146.
-98.
-11.
78.
154.
211.
256.
196.
202.
253.
190.
103.
2.
-79.
-146.
-198.
-232.
-261.
-267.
-230.
-146.
-91.
-26.
55.
126.
179.
202.
162.
151.
201.
176.
97.
2.
-71.
-131.
-178.
83
77
96
03
05
00
91
99
67
01
73
44
16
00
33
86
80
68
87
17
33
35
21
62
30
01
34
47
42
72
61
30
36
00
09
43
05
223.
217.
288.
302.
263.
201.
123.
38.
-33.
-107.
-185.
-256.
-295.
-278.
-250.
-196.
-127.
-47.
34.
121.
194.
179.
229.
275.
248.
197.
129.
52.
-23.
-91.
-163.
-232.
-267.
-251.
-230.
-182.
-120.
90
43
04
58
18
72
27
11
00
64
45
50
08
00
67
29
07
04
64
16
77
65
83
52
31
64
32
18
57
57
01
57
36
00
09
39
85
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3198.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3231.
3264.
3264.
3264.
3264.
42
42
42
42
42
42
42
42
42
42
42
42
42
23
23
23
23
23
23
23
23
23
23
23
23
23
23
23
23
23
23
23
23
04
04
04
04
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
                        52

-------
               ** Cinder Cone Butte
           ** Low Stack Height: 30 meters
High Buoyancy Case: 6-m diameter, 30 m/s Exit Vel, 500 K
             ** Close to Hill: 1 km away
     ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
-201.
-230.
-237.
-214.
-140.
-81.
-31.
41.
105.
155.
163.
135.
113.
149.
151.
84.
2.
-61.
-114.
-153.
-180.
-206.
-215.
-195.
-129.
-74.
-32.
31.
86.
130.
136.
114.
85.
106.
125.
73.
3.
61
06
00
19
21
40
28
60
42
35
34
40
60
22
97
02
00
01
04
87
97
06
00
31
76
20
19
93
72
34
25
98
47
04
03
77
00
-47.
24.
100.
169.
150.
187.
243.
231.
187.
129.
61.
-9.
-71.
-140.
-204.
-240.
-231.
-208.
-168.
-115.
-54.
7.
73.
130.
120.
152.
204.
209.
172.
123.
61.
-2.
-49.
-113.
-175.
-215.
-206.
29
39
64
15
71
62
37
33
96
97
63
50
80
06
88
37
00
49
94
63
34
87
91
00
63
79
87
23
28
48
17
29
87
29
86
85
00
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3264
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3296
3329
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.04
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.85
.66
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
                        53

-------
               ** Cinder Cone Butte
           ** Low Stack Height: 30 meters
High Buoyancy Case: 6-m diameter, 30 m/s Exit Vel, 500 K
             ** Close to Hill: 1 km away
     ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
-53.
-99.
-134.
-159.
-179.
-187.
-156.
-99.
-55.
-17.
36.
83.
114.
104.
72.
50.
69.
88.
65.
2.
-51.
-94.
-127.
-143.
-161.
-125.
-72.
-32.
7.
59.
87.
83.
41.
-12.
-54.
-17.
36.
07
24
94
94
33
72
80
24
52
95
11
25
87
12
48
78
73
06
18
00
36
97
05
69
36
43
68
04
68
05
00
26
54
50
22
83
25
-188.
-153.
-107.
-54.
1.
58.
93.
95.
134.
179.
181.
147.
99.
42.
-4.
-31.
-83.
-138.
-187.
-181.
-161.
-124.
-78.
-23.
30.
62.
83.
122.
162.
141.
94.
38.
8.
25.
2.
-30.
-48.
80
39
29
13
32
87
46
55
37
52
51
75
51
87
28
31
64
92
19
00
23
59
21
72
37
70
20
40
36
43
58
04
84
13
65
31
00
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3329
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
3362
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.66
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
.46
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
                        54

-------
                                ** Cinder Cone Butte
                            ** Low Stack Height: 30 meters
                  High Buoyancy Case: 6-m diameter, 30 m/s Exit Vel, 500 K
                              ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
52
57
4
-22
-46
-67
-82
-88
-66
-40
-14
10
23
23
7
15
3
9
34
55
65
57
.65
.52
.00
.33
.82
.86
.70
.59
.55
.90
.47
.97
.83
.86
.13
.65
.80
.99
.87
.86
.52
.11
-100.
-155.
-146.
-143.
-133.
-117.
-95.
-70.
-64.
-69.
-72.
-79.
-101.
-127.
-144.
73.
96.
119.
124.
108.
85.
61.
60
07
00
15
11
36
60
00
01
61
14
29
15
28
96
30
84
90
11
49
47
70
3362
3362
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
.46
.46
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
3395
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
.27
RE FINISHED

ME STARTING
   SURFFILE   SURFACE  FREE
   PROFFILE   PROFILE  FREE
   SURFDATA   99999  1994  Tower/Sodar
   UAIRDATA   99999  1994  Tower/Sodar
   SITEDATA   00000  1994  Tower/Sodar
ME FINISHED

OU STARTING
   RECTABLE   ALLAVE  FIRST-SECOND
OU FINISHED
                                         55

-------
                             ** Montour Alongwind Hill
                           ** High Stack Height: 150 meters
                ** Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
                             ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar

** Montour Alongwind  Hill
** High Stack Height:  150  meters
** Low Buoyancy Case:  2-m  diameter,  10  m/s Exit Vel, 400 K
** Close to Hill: 1 km away
** Meteorology from 100-m  Tower  and Sodar
**
CO STARTING
   TITLEONE  Montour  Alongwind Hill:  Source 1 km away
   TITLETWO  150-m Stack Height;  Low Buoyancy Case
   MODELOPT  CONG              MSGPRO
             1 3 24 Period
             S02
             RUN
             ERRORS.OUT
               ELEV
   AVERTIME
   POLLUTID
   RUNORNOT
   ERRORFIL
**   TERRHGTS
CO FINISHED
SO STARTING
   ELEVUNIT  FEET
   LOCATION  STACK1
** Point Source
** Parameters:
   SRCPARAM  STACK1
                     POINT
                      QS
-3500.0
HS    TS
0.0
  VS
700.0
 DS
                        1.   150.   400.   10.0  2.0
   SRCGROUP  ALL
SO FINISHED
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
STARTING
ELEVUNIT
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART

FEET
1655
1075
481
-125
-747
-1367
-1332
-761


.00
.54
.15
.99
.97
.68
.89
.22


146
-98
-308
-475
-576
-518
-67
195


.00
.85
.48
.98
.12
.96
.94
.56


1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.


00
00
00
00
00
00
00
00


1400
1400
1400
1400
1400
1400
1400
1400
                                        56

-------
             ** Montour Alongwind Hill
          ** High Stack Height: 150 meters
Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
             ** Close to Hill: 1 km away
    ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
-149
465
1093
1712
2281
2905
3514
4142
4040
3413
2851
2224
141
-108
-358
-607
-857
-1083
-1074
-875
-637
-395
-152
92
342
592
843
1081
1102
886
649
420
-373
-287
-188
-89
9
.69
.45
.57
.61
.57
.52
.58
.81
.15
.23
.31
.11
.00
.32
.91
.97
.60
.05
.53
.96
.92
.29
.75
.72
.50
.42
.23
.96
.55
.71
.92
.08
.00
.24
.10
.24
.09
347
485
506
443
691
783
925
947
665
607
427
416
-299
-341
-376
-421
-453
-366
-137
10
95
168
240
301
342
378
366
287
95
-29
-118
-223
-31
27
58
90
119
.21
.10
.82
.97
.14
.48
.91
.59
.33
.18
.99
.49
.00
.14
.51
.69
.09
.88
.83
.86
.96
.10
.59
.83
.65
.00
.02
.07
.11
.07
.37
.10
.00
.50
.73
.31
.15
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1400.
1400.
1400.
1400.
1400.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
                        57

-------
                              ** Montour Alongwind Hill
                           ** High Stack Height: 150 meters
                  Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
                              ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
113
214
311
410
405
337
244
152
53
-48
-151
-253
-354
-450
-460
.04
.03
.84
.80
.63
.71
.77
.08
.15
.48
.63
.71
.07
.53
.63
125
144
177
157
69
-8
-55
-102
-132
-155
-169
-189
-217
-194
-106
.22
.82
.82
.11
.59
.29
.12
.54
.88
.28
.15
.72
.34
.01
.35
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
RE FINISHED

ME STARTING
   SURFFILE   SURFACE  FREE
   PROFFILE   PROFILE  FREE
   SURFDATA   99999  1994  Tower/Sodar
   UAIRDATA   99999  1994  Tower/Sodar
   SITEDATA   00000  1994  Tower/Sodar
   PROFBASE  0
ME FINISHED
OU STARTING
   RECTABLE
OU FINISHED
ALLAVE  FIRST-SECOND
                                         58

-------
                             ** Montour Crosswind Hill
                           ** High Stack Height: 150 meters
                ** Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
                             ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar

** Montour Crosswind  Hill
** High Stack Height:  150  meters
** Low Buoyancy Case:  2-m  diameter,  10  m/s  Exit Vel, 400 K
** Close to Hill: 1 km away
** Meteorology from 100-m  Tower  and Sodar
**
CO STARTING
   TITLEONE  Montour  Crosswind Hill:  Source 1 km away
   TITLETWO  150-m Stack Height;  Low Buoyancy Case
   MODELOPT  CONG              MSGPRO
             1 3 24 Period
             S02
             RUN
             ERRORS.OUT
               ELEV
   AVERTIME
   POLLUTID
   RUNORNOT
   ERRORFIL
**   TERRHGTS
CO FINISHED
SO STARTING
   ELEVUNIT  FEET
   LOCATION  STACK1
** Point Source
** Parameters:
   SRCPARAM  STACK1
                     POINT   0.0    2000.0   700.0
                      QS     HS     TS    VS   DS

                        1.   150.   400.   10.0  2.0
   SRCGROUP  ALL
SO FINISHED
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
STARTING
ELEVUNIT
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART

FEET
1655
1075
481
-125
-747
-1367
-1332
-761


.00
.54
.15
.99
.97
.68
.89
.22


146
-98
-308
-475
-576
-518
-67
195


.00
.85
.48
.98
.12
.96
.94
.56


1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.


00
00
00
00
00
00
00
00


1400
1400
1400
1400
1400
1400
1400
1400
                                        59

-------
             ** Montour Crosswind Hill
          ** High Stack Height: 150 meters
Low Buoyancy Case: 2-m diameter,  10 m/s Exit Vel, 400 K
             ** Close to Hill: 1 km away
    ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
-149
465
1093
1712
2281
2905
3514
4142
4040
3413
2851
2224
141
-108
-358
-607
-857
-1083
-1074
-875
-637
-395
-152
92
342
592
843
1081
1102
886
649
420
-373
-287
-188
-89
9
.69
.45
.57
.61
.57
.52
.58
.81
.15
.23
.31
.11
.00
.32
.91
.97
.60
.05
.53
.96
.92
.29
.75
.72
.50
.42
.23
.96
.55
.71
.92
.08
.00
.24
.10
.24
.09
347
485
506
443
691
783
925
947
665
607
427
416
-299
-341
-376
-421
-453
-366
-137
10
95
168
240
301
342
378
366
287
95
-29
-118
-223
-31
27
58
90
119
.21
.10
.82
.97
.14
.48
.91
.59
.33
.18
.99
.49
.00
.14
.51
.69
.09
.88
.83
.86
.96
.10
.59
.83
.65
.00
.02
.07
.11
.07
.37
.10
.00
.50
.73
.31
.15
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1200.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1300.
1400.
1400.
1400.
1400.
1400.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
                        60

-------
                              ** Montour Crosswind Hill
                           ** High Stack Height: 150 meters
                  Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
                              ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
113
214
311
410
405
337
244
152
53
-48
-151
-253
-354
-450
-460
.04
.03
.84
.80
.63
.71
.77
.08
.15
.48
.63
.71
.07
.53
.63
125
144
177
157
69
-8
-55
-102
-132
-155
-169
-189
-217
-194
-106
.22
.82
.82
.11
.59
.29
.12
.54
.88
.28
.15
.72
.34
.01
.35
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
1400.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
1400
RE FINISHED

ME STARTING
   SURFFILE   SURFACE  FREE
   PROFFILE   PROFILE  FREE
   SURFDATA   99999  1994  Tower/Sodar
   UAIRDATA   99999  1994  Tower/Sodar
   SITEDATA   00000  1994  Tower/Sodar
   PROFBASE  0
ME FINISHED
OU STARTING
   RECTABLE
OU FINISHED
ALLAVE  FIRST-SECOND
                                         61

-------
                                 ** Piedmont Hill
                           ** High Stack Height: 150 meters
                  Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
                             ** Close to Hill: 1  km away
                      ** Meteorology  from 100-m Tower and Sodar
** Piedmont Hill
** High Stack Height:
** Low Buoyancy Case:
                      150 meters
                      2-m diameter,  10 m/s Exit Vel, 400 K
** Close to Hill: 1 km away
** Meteorology from 100-m Tower and  Sodar
CO STARTING
   TITLEONE
   TITLETWO
   MODELOPT
   AVERTIME
   POLLUTID
   RUNORNOT
   ERRORFIL
**   TERRHGTS
CO FINISHED
             Piedmont Hill:  Source 1 km away
             150-m Stack Height;  Low Buoyancy Case
             CONG              MSGPRO
             1 3 24 Period
             S02
             RUN
             ERRORS.OUT
               ELEV
SO STARTING
   ELEVUNIT  FEET
   LOCATION  STACK1
** Point Source
** Parameters:
   SRCPARAM  STACK1
                     POINT   0.0    1000.0   1000.0
                      QS     HS     TS    VS   DS

                        1.   150.   400.   10.0  2.0
   SRCGROUP  ALL
SO FINISHED
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
STARTING
ELEVUNIT
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART








FEET
-612
-632
-495
-433
-106
245
413
553
.40
.69
.55
.60
.68
.80
.78
.12
-1792
-1395
-1025
-668
-522
-585
-230
138
.00
.44
.07
.18
.37
.23
.98
.16
1500
1500
1500
1500
1500
1500
1500
1500
.00
.00
.00
.00
.00
.00
.00
.00
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
00
00
00
00
00
00
00
00
                                        62

-------
                 ** Piedmont Hill
          ** High Stack Height: 150 meters
Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
             ** Close to Hill: 1 km away
    ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
882.
1211.
1277.
1032.
-582.
-561.
-479.
-341.
-133.
242.
448.
592.
918.
1228.
1088.
885.
-505.
-480.
-465.
-310.
-138.
231.
469.
618.
944.
1164.
962.
783.
-357.
-437.
-369.
-285.
-166.
192.
486.
638.
945.
78
30
72
05
60
02
65
53
80
75
16
43
60
23
95
35
50
41
41
06
32
42
14
92
57
28
35
63
70
57
41
45
77
72
38
29
65
337.
132.
-198.
-511.
-2015.
-1616.
-1223.
-863.
-612.
-647.
-311.
59.
238.
4.
-350.
-698.
-2155.
-1760.
-1358.
-987.
-654.
-710.
-408.
-32.
168.
-117.
-467.
-831.
-2305.
-1945.
-1566.
-1172.
-788.
-801.
-562.
-189.
31.
97
16
82
55
00
37
82
47
84
07
09
94
86
15
46
68
00
47
34
10
32
58
24
23
79
76
47
20
00
98
25
66
53
51
30
10
54
1500
1500
1500
1500
1600
1600
1600
1600
1600
1600
1600
1600
1600
1600
1600
1600
1700
1700
1700
1700
1700
1700
1700
1700
1700
1700
1700
1700
1800
1800
1800
1800
1800
1800
1800
1800
1800
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
                        63

-------
                 ** Piedmont Hill
          ** High Stack Height: 150 meters
Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
             ** Close to Hill: 1 km away
    ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
1040.
845.
700.
-312.
-352.
-284.
-219.
-103.
231.
526.
675.
957.
856.
728.
646.
-211.
-208.
-103.
-133.
-159.
-71.
67.
296.
524.
573.
607.
677.
99.
52.
34.
93.
-9.
2.
127.
276.
443.
421.
60
76
78
10
85
39
41
65
54
84
48
97
22
44
54
00
98
15
27
02
30
68
79
56
38
64
99
85
53
82
40
02
67
68
83
05
80
-246.
-601.
-978.
-2275.
-1926.
-1572.
-1204.
-845.
-848.
-634.
-292.
-77.
-432.
-787.
-1155.
-2248.
-2019.
-1812.
-1588.
-1357.
-1142.
-957.
-912.
-859.
-1049.
-1272.
-1493.
-2515.
-2279.
-2043.
-1816.
-1599.
-1394.
-1190.
-1010.
-1102.
-1337.
59
19
55
00
16
42
49
19
17
35
95
44
11
98
69
00
81
17
07
53
46
10
52
69
63
62
23
00
43
20
61
99
18
55
06
06
79
1800
1800
1800
1900
1900
1900
1900
1900
1900
1900
1900
1900
1900
1900
1900
2000
2000
2000
2000
2000
2000
2000
2000
2000
2000
2000
2000
2100
2100
2100
2100
2100
2100
2100
2100
2100
2100
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
00
                        64

-------
                                  ** Piedmont Hill
                           ** High Stack Height: 150 meters
                  Low Buoyancy Case: 2-m diameter, 10 m/s Exit Vel, 400 K
                              ** Close to Hill: 1 km away
                      ** Meteorology from 100-m Tower and Sodar
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
RE
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
DISCCART
613
643
348
323
375
308
340
409
387
296
440
552
606
666
.92
.13
.70
.75
.69
.07
.42
.93
.04
.38
.59
.40
.70
.27
-1451.
-1689.
-2665.
-2501.
-2350.
-2194.
-2038.
-1883.
-1722.
-1584.
-1521.
-1613.
-1773.
-1932.
98
91
00
71
23
93
55
47
33
62
57
29
86
25
2100
2100
2200
2200
2200
2200
2200
2200
2200
2200
2200
2200
2200
2200
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
.00
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
2200.
00
00
00
00
00
00
00
00
00
00
00
00
00
00
RE FINISHED
ME STARTING
   SURFFILE   SURFACE  FREE
   PROFFILE   PROFILE  FREE
   SURFDATA   99999   1994  Tower/Sodar
   UAIRDATA   99999   1994  Tower/Sodar
   SITEDATA   00000   1994  Tower/Sodar
   PROFBASE  0
ME FINISHED
OU STARTING
   RECTABLE
OU FINISHED
ALLAVE  FIRST-SECOND
                                         65

-------
APPENDIX D
FLAT AND SIMPLE TERRAIN MODELING RESULTS -





LISTING OF REGULATORY DESIGN CONCENTRATIONS
                      66

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
R05NDNBO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R05NDNBP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R10NDNBO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R10NDNBP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

1.477E+06
1.477E+06
9.403E+05
6.626E+05
1.678E+05
1.355E+05
1.537E+04
ISC3

2.215E+06
1.477E+06
8.294E+05
6.729E+05
3.350E+05
1.503E+05
9.266E+03
ISC3

8.679E+05
8.679E+05
5.524E+05
3.918E+05
1.018E+05
8.792E+04
1.147E+04
ISC3

1.283E+06
8.679E+05
4.871 E+05
3.966E+05
2.044E+05
9.117E+04
7.044E+03
99351AER/
ISC3
0.434
0.432
0.454
0.461
0.539
0.591
0.388
99351AER/
ISC3
0.357
0.433
0.608
0.427
0.347
0.472
0.500
99351AER/
ISC3
0.442
0.440
0.484
0.574
0.694
0.645
0.515
99351AER/
ISC3
0.355
0.441
0.706
0.542
0.399
0.550
0.615
02222AER/
ISC3
0.673
0.668
0.688
0.574
0.596
0.730
0.434
02222AER/
ISC3
0.619
0.671
0.783
0.629
0.510
0.687
0.593
02222AER/
ISC3
0.504
0.502
0.541
0.591
0.592
0.605
0.529
02222AER/
ISC3
0.410
0.503
0.764
0.546
0.438
0.613
0.661
02222AER/
99351AER
1.550
1.546
1.516
1.244
1.106
1.235
1.119
02222AER/
99351AER
1.734
1.548
1.288
1.473
1.470
1.455
1.185
02222AER/
99351AER
1.139
1.140
1.118
1.030
0.852
0.937
1.027
02222AER/
99351AER
1.155
1.140
1.082
1.009
1.098
1.115
1.075
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               67

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
R20NDNBO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R20NDNBP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R35MFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R35MFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

7.518E+04
5.041 E+04
3.027E+04
1.975E+04
8.005E+03
5.218E+03
8.692E+02
ISC3

7.595E+04
4.889E+04
3.242E+04
2.192E+04
7.737E+03
6.170E+03
6.952E+02
ISC3

1.394E+03
1.257E+03
1.091E+03
9.733E+02
3.291 E+02
3.022E+02
4.046E+01
ISC3

1.394E+03
1.267E+03
9.518E+02
7.229E+02
3.393E+02
2.322E+02
2.781 E+01
99351AER/
ISC3
0.510
0.551
0.735
0.875
1.076
1.496
2.277
99351AER/
ISC3
0.524
0.709
0.745
0.835
1.207
1.254
1.926
99351AER/
ISC3
1.093
1.103
0.971
0.932
1.624
1.723
2.005
99351AER/
ISC3
1.215
1.193
1.207
1.181
1.206
1.295
1.835
02222AER/
ISC3
0.480
0.602
0.621
0.850
1.052
1.417
2.269
02222AER/
ISC3
0.508
0.618
0.682
0.767
1.166
1.242
1.931
02222AER/
ISC3
1.652
1.192
0.989
0.932
1.581
1.689
2.002
02222AER/
ISC3
1.823
1.326
1.154
1.087
1.262
1.260
1.815
02222AER/
99351AER
0.941
1.092
0.846
0.971
0.978
0.947
0.996
02222AER/
99351AER
0.970
0.873
0.916
0.918
0.966
0.990
1.002
02222AER/
99351AER
1.511
1.081
1.019
1.000
0.974
0.980
0.998
02222AER/
99351AER
1.500
1.111
0.956
0.920
1.046
0.97
0.989
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               68

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
R35BFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R35BFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R100FO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R100FP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

2.651 E+02
2.551 E+02
1.832E+02
1.807E+02
8.005E+01
6.871 E+01
8.103E+00
ISC3

2.626E+02
2.608E+02
1.951 E+02
1.394E+02
6.414E+01
5.184E+01
5.464E+00
ISC3

6.150E+01
5.391 E+01
3.820E+01
2.791 E+01
6.271 E+00
5.537E+00
6.320E-01
ISC3

5.296E+01
5.003E+01
2.751 E+01
1.826E+01
7.965E+00
4.170E+00
4.000E-01
99351AER/
ISC3
1.589
1.623
2.142
2.081
2.534
2.721
3.121
99351AER/
ISC3
1.575
1.537
2.006
2.411
2.650
2.995
3.388
99351AER/
ISC3
0.917
0.901
1.153
1.413
2.107
2.347
2.696
99351AER/
ISC3
0.931
0.925
1.474
1.907
1.759
2.628
3.253
02222AER/
ISC3
1.589
1.622
2.142
2.081
2.463
2.714
3.054
02222AER/
ISC3
1.575
1.537
2.006
2.411
2.618
2.936
3.322
02222AER/
ISC3
0.816
0.805
1.020
1.290
2.128
2.313
2.663
02222AER/
ISC3
0.815
0.816
1.309
1.724
1.622
2.395
3.203
02222AER/
99351AER
1.000
1.000
1.000
1.000
0.972
0.998
0.979
02222AER/
99351AER
1.000
1.000
1.000
1.000
0.988
0.980
0.981
02222AER/
99351AER
0.889
0.894
0.884
0.913
1.010
0.986
0.988
02222AER/
99351AER
0.875
0.882
0.888
0.904
0.922
0.911
0.985
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               69

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
R200FO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R200FP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R35TO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R35TP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

3.441 E+01
2.074E+01
1.335E+01
1.043E+01
2.594E+00
1.907E+00
1.480E-01
ISC3

3.479E+01
3.1 31 E+01
1.160E+01
1.098E+01
3.174E+00
1.714E+00
8.200E-02
ISC3

2.863E+02
2.522E+02
1.606E+02
1.471E+02
5.868E+01
5.074E+01
6.466E+00
ISC3

2.984E+02
2.835E+02
2.087E+02
1.668E+02
6.858E+01
5.750E+01
6.193E+00
99351AER/
ISC3
1.036
0.989
1.212
1.044
1.414
1.781
3.162
99351AER/
ISC3
0.836
0.583
1.015
0.886
1.223
1.824
3.890
99351AER/
ISC3
1.111
1.255
1.903
2.039
3.151
3.411
3.394
99351AER/
ISC3
1.040
1.085
1.366
1.694
2.141
2.227
2.289
02222AER/
ISC3
0.941
0.922
1.091
0.929
1.328
1.716
3.088
02222AER/
ISC3
0.721
0.495
0.938
0.818
1.154
1.630
3.829
02222AER/
ISC3
1.045
1.149
1.608
1.550
2.414
2.494
2.526
02222AER/
ISC3
1.025
1.050
1.356
1.544
2.132
2.182
2.265
02222AER/
99351AER
0.908
0.932
0.900
0.889
0.939
0.963
0.976
02222AER/
99351AER
0.863
0.850
0.924
0.923
0.944
0.894
0.984
02222AER/
99351AER
0.940
0.915
0.845
0.760
0.766
0.731
0.744
02222AER/
99351AER
0.986
0.967
0.993
0.912
0.996
0.979
0.990
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               70

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.     The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
R200TO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R200TP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R10VOLFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R10VOLFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

3.441 E+01
2.226E+01
1.723E+01
1.446E+01
5.957E+00
3.802E+00
5.140E-01
ISC3

3.479E+01
3.130E+01
1.512E+01
1.299E+01
4.775E+00
4.005E+00
3.440E-01
ISC3

6.407E+04
4.611E+04
3.675E+04
3.094E+04
1.342E+04
1.102E+04
2.496E+03
ISC3

8.044E+04
8.004E+04
4.150E+04
3.422E+04
1.719E+04
1.191E+04
1.656E+03
99351AER/
ISC3
1.004
0.910
0.923
0.746
0.626
0.931
0.930
99351AER/
ISC3
0.830
0.577
0.755
0.774
0.852
0.812
1.073
99351AER/
ISC3
0.909
1.249
1.154
1.296
1.142
1.150
1.059
99351AER/
ISC3
0.900
0.768
1.276
1.168
1.070
1.025
1.026
02222AER/
ISC3
0.918
0.888
0.831
0.677
0.549
0.828
0.860
02222AER/
ISC3
0.720
0.495
0.705
0.723
0.811
0.728
1.058
02222AER/
ISC3
0.909
1.249
1.154
1.206
1.054
1.109
1.054
02222AER/
ISC3
0.900
0.768
1.276
1.168
1.063
1.010
1.039
02222AER/
99351AER
0.915
0.976
0.900
0.907
0.878
0.890
0.925
02222AER/
99351AER
0.867
0.859
0.934
0.933
0.952
0.897
0.986
02222AER/
99351AER
1.000
1.000
1.000
0.931
0.923
0.964
0.995
02222AER/
99351AER
1.000
1.000
1.000
1.000
0.994
0.986
1.012
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               71

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
R35VOLFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
R35VOLFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
RAREFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
RAREFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

1.347E+04
8.979E+03
6.060E+03
4.795E+03
2.022E+03
1.728E+03
4.123E+02
ISC3

1.457E+04
1.024E+04
6.079E+03
5.010E+03
2.178E+03
1.807E+03
2.988E+02
ISC3

1.374E+04
1.374E+04
1.317E+04
1.218E+04
6.654E+03
5.904E+03
2.573E+03
ISC3

2.047E+04
1.983E+04
1.532E+04
1.257E+04
7.766E+03
7.353E+03
3.044E+03
99351AER/
ISC3
0.683
0.962
0.906
0.932
0.890
0.889
0.925
99351AER/
ISC3
0.905
0.999
1.026
0.973
0.810
0.967
0.915
99351AER/
ISC3
2.040
1.382
1.379
1.201
1.169
0.995
0.646
99351AER/
ISC3
1.509
1.515
1.280
1.324
1.057
1.050
0.715
02222AER/
ISC3
0.619
0.875
0.807
0.880
0.897
0.840
0.920
02222AER/
ISC3
0.818
0.847
0.956
0.868
0.807
0.930
0.913
02222AER/
ISC3
1.954
1.312
1.305
1.132
1.019
0.832
0.614
02222AER/
ISC3
1.441
1.441
1.252
1.292
1.010
0.989
0.682
02222AER/
99351AER
0.906
0.910
0.891
0.945
1.008
0.945
0.994
02222AER/
99351AER
0.904
0.848
0.932
0.892
0.996
0.961
0.998
02222AER/
99351AER
0.958
0.950
0.946
0.943
0.872
0.836
0.950
02222AER/
99351AER
0.955
0.951
0.978
0.975
0.955
0.942
0.954
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               72

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
U05NDNBO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U05NDNBP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U10NDNBO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U10NDNBP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

1.668E+05
1.668E+05
1.284E+05
9.909E+04
3.038E+04
2.785E+04
5.129E+03
ISC3

2.498E+05
1.668E+05
1.382E+05
1.106E+05
5.385E+04
2.891 E+04
3.033E+03
ISC3

1.588E+05
1.588E+05
1.223E+05
9.443E+04
2.899E+04
2.661 E+04
4.942E+03
ISC3

2.375E+05
1.588E+05
1.315E+05
1.053E+05
5.135E+04
2.774E+04
2.923E+03
99351AER/
ISC3
2.096
2.078
1.869
2.130
2.103
1.818
0.992
99351AER/
ISC3
2.743
3.152
2.318
1.841
1.537
1.502
1.260
99351AER/
ISC3
1.127
1.125
1.061
1.322
1.497
1.163
0.941
99351AER/
ISC3
1.255
1.696
1.319
1.136
0.949
1.344
1.153
02222AER/
ISC3
0.887
0.852
0.804
0.813
1.091
1.041
1.261
02222AER/
ISC3
0.592
0.873
0.857
0.849
0.801
1.013
1.412
02222AER/
ISC3
0.533
0.516
0.523
0.634
0.864
0.860
1.187
02222AER/
ISC3
0.356
0.525
0.575
0.598
0.570
0.809
1.342
02222AER/
99351AER
0.423
0.410
0.430
0.382
0.519
0.573
1.271
02222AER/
99351AER
0.216
0.277
0.370
0.461
0.521
0.674
1.120
02222AER/
99351AER
0.473
0.459
0.493
0.479
0.577
0.739
1.262
02222AER/
99351AER
0.284
0.310
0.436
0.526
0.601
0.602
1.164
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               73

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
U20NDNBO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U20NDNBP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U35MFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U35MFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

5.789E+04
5.078E+04
3.907E+04
3.145E+04
1.145E+04
9.885E+03
2.119E+03
ISC3

7.568E+04
5.516E+04
4.189E+04
3.361 E+04
1.730E+04
1.089E+04
1.398E+03
ISC3

2.222E+03
1.859E+03
1.462E+03
1.403E+03
8.936E+02
8.090E+02
1.415E+02
ISC3

2.358E+03
1.859E+03
1.644E+03
1.401E+03
9.331 E+02
7.561 E+02
1.037E+02
99351AER/
ISC3
0.522
0.422
0.485
0.432
0.658
0.713
1.013
99351AER/
ISC3
0.365
0.482
0.443
0.462
0.468
0.682
0.997
99351AER/
ISC3
0.691
0.690
0.692
0.683
0.595
0.644
0.795
99351AER/
ISC3
0.694
0.820
0.644
0.604
0.466
0.477
0.611
02222AER/
ISC3
0.582
0.588
0.384
0.431
0.511
0.544
0.670
02222AER/
ISC3
0.438
0.535
0.528
0.500
0.346
0.469
0.799
02222AER/
ISC3
0.567
0.638
0.567
0.558
0.456
0.441
0.401
02222AER/
ISC3
0.833
0.676
0.574
0.516
0.272
0.303
0.399
02222AER/
99351AER
1.116
1.392
0.792
0.996
0.777
0.763
0.662
02222AER/
99351AER
1.200
1.110
1.192
1.081
0.739
0.687
0.801
02222AER/
99351AER
0.820
0.925
0.820
0.818
0.766
0.684
0.504
02222AER/
99351AER
1.201
0.824
0.891
0.854
0.583
0.635
0.653
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               74

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
U35BFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U35BFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U100FO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U100FP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

6.987E+02
5.931 E+02
4.513E+02
3.648E+02
2.849E+02
2.716E+02
3.748E+01
ISC3

5.718E+02
4.429E+02
3.615E+02
3.473E+02
2.1 71 E+02
1.927E+02
2.568E+01
ISC3

8.344E+01
7.354E+01
5.392E+01
4.886E+01
3.015E+01
2.984E+01
4.082E+00
ISC3

8.015E+01
7.214E+01
5.665E+01
4.963E+01
2.329E+01
2.022E+01
2.602E+00
99351AER/
ISC3
0.686
0.754
0.925
1.111
1.027
0.997
0.975
99351AER/
ISC3
0.787
0.990
1.128
1.094
1.046
1.066
1.027
99351AER/
ISC3
0.644
0.709
0.912
0.953
0.601
0.519
0.549
99351AER/
ISC3
0.664
0.703
0.812
0.820
0.736
0.670
0.648
02222AER/
ISC3
0.539
0.621
0.724
0.841
0.548
0.477
0.580
02222AER/
ISC3
0.615
0.773
0.866
0.894
0.593
0.614
0.649
02222AER/
ISC3
0.601
0.591
0.722
0.737
0.383
0.367
0.392
02222AER/
ISC3
0.571
0.607
0.636
0.634
0.555
0.494
0.479
02222AER/
99351AER
0.786
0.823
0.782
0.758
0.533
0.479
0.594
02222AER/
99351AER
0.782
0.781
0.768
0.817
0.567
0.576
0.631
02222AER/
99351AER
0.934
0.833
0.792
0.773
0.638
0.707
0.714
02222AER/
99351AER
0.859
0.864
0.783
0.773
0.754
0.738
0.739
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               75

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
U200FO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U200FP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U35TO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U35TP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

3.134E+01
2.972E+01
2.152E+01
1.978E+01
1.107E+01
9.360E+00
1.276E+00
ISC3

3.104E+01
2.726E+01
2.028E+01
1.941E+01
8.419E+00
6.496E+00
8.800E-01
ISC3

4.808E+02
4.364E+02
3.433E+02
2.541 E+02
1.271E+02
1.107E+02
1.826E+01
ISC3

6.314E+02
4.922E+02
3.071 E+02
2.784E+02
1.516E+02
1.333E+02
2.224E+01
99351AER/
ISC3
0.952
0.614
0.651
0.562
0.382
0.449
0.490
99351AER/
ISC3
0.804
0.570
0.550
0.559
0.524
0.537
0.527
99351AER/
ISC3
0.556
0.604
0.749
0.994
1.295
1.467
1.520
99351AER/
ISC3
0.420
0.522
0.805
0.884
1.053
0.961
0.762
02222AER/
ISC3
1.033
0.643
0.760
0.576
0.314
0.350
0.378
02222AER/
ISC3
0.935
0.920
0.644
0.550
0.688
0.527
0.435
02222AER/
ISC3
0.465
0.477
0.533
0.705
0.883
0.898
0.759
02222AER/
ISC3
0.367
0.461
0.698
0.731
0.664
0.749
0.554
02222AER/
99351AER
1.085
1.048
1.167
1.025
0.822
0.779
0.771
02222AER/
99351AER
1.163
1.613
1.171
0.984
1.312
0.981
0.825
02222AER/
99351AER
0.838
0.790
0.711
0.710
0.682
0.612
0.499
02222AER/
99351AER
0.875
0.884
0.866
0.826
0.631
0.780
0.728
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               76

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
U200TO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U200TP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U10VOLFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U10VOLFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

7.411E+01
6.340E+01
5.143E+01
4.184E+01
1.332E+01
1.160E+01
1.855E+00
ISC3

7.576E+01
7.495E+01
4.611E+01
4.378E+01
1.884E+01
1.631E+01
1.673E+00
ISC3

2.722E+04
2.722E+04
2.225E+04
2.007E+04
8.807E+03
7.552E+03
1.775E+03
ISC3

4.076E+04
3.631 E+04
2.595E+04
2.219E+04
1.102E+04
7.792E+03
1.106E+03
99351AER/
ISC3
0.385
0.320
0.260
0.267
0.328
0.369
0.343
99351AER/
ISC3
0.323
0.279
0.275
0.264
0.238
0.218
0.304
99351AER/
ISC3
1.140
1.132
1.104
1.189
1.261
1.159
1.319
99351AER/
ISC3
1.070
1.114
1.148
1.062
1.043
1.108
1.388
02222AER/
ISC3
0.426
0.312
0.326
0.282
0.262
0.296
0.277
02222AER/
ISC3
0.389
0.334
0.266
0.244
0.330
0.203
0.231
02222AER/
ISC3
1.040
1.027
1.007
1.070
1.146
1.151
1.348
02222AER/
ISC3
0.850
0.861
1.031
0.990
0.974
1.109
1.492
02222AER/
99351AER
1.106
0.974
1.253
1.055
0.798
0.803
0.807
02222AER/
99351AER
1.205
1.195
0.966
0.925
1.384
0.933
0.760
02222AER/
99351AER
0.913
0.908
0.912
0.900
0.909
0.993
1.022
02222AER/
99351AER
0.794
0.773
0.898
0.932
0.934
1.001
1.075
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               77

-------
Appendix D.    Results of  the  flat  and simple  terrain analysis.    The ISC3
concentrations are  in ug/m3 and the last  3 columns  are  concentration ratios.
U35VOLFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
U35VOLFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
UAREFO

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
UAREFP

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
ISC3

1.079E+04
8.430E+03
5.600E+03
5.198E+03
2.736E+03
2.149E+03
5.942E+02
ISC3

1.079E+04
9.609E+03
6.187E+03
5.280E+03
2.839E+03
2.302E+03
4.048E+02
ISC3

4.902E+03
4.902E+03
4.787E+03
4.411E+03
2.651 E+03
2.276E+03
1.095E+03
ISC3

7.322E+03
7.202E+03
5.631 E+03
4.691 E+03
3.005E+03
2.882E+03
1.287E+03
99351AER/
ISC3
0.768
0.919
0.890
0.783
0.594
0.681
0.653
99351AER/
ISC3
1.114
0.981
0.939
0.868
0.554
0.634
0.680
99351AER/
ISC3
4.254
2.856
2.820
2.405
2.168
1.823
1.151
99351AER/
ISC3
3.157
3.150
2.576
2.668
1.893
1.954
1.272
02222AER/
ISC3
0.773
0.932
0.873
0.790
0.593
0.639
0.593
02222AER/
ISC3
1.105
0.903
0.940
0.824
0.533
0.610
0.650
02222AER/
ISC3
1.383
1.247
1.224
1.246
1.281
1.373
1.434
02222AER/
ISC3
0.958
0.932
1.015
1.178
1.201
1.237
1.417
02222AER/
99351AER
1.007
1.014
0.980
1.009
0.998
0.938
0.909
02222AER/
99351AER
0.992
0.921
1.000
0.949
0.962
0.963
0.955
02222AER/
99351AER
0.325
0.437
0.434
0.518
0.591
0.753
1.246
02222AER/
99351AER
0.303
0.296
0.394
0.441
0.635
0.633
1.114
The first letter indicates urban (U) or rural (R) settings: the number indicates the stack height (meters): the next set of
letters - VOL, volume source; ARE, area source; ND, no downwash; NB, no building; M, moderately buoyant; B,
bouyant plume; F, flat terrain; T, simple terrain: the last letter indicates the met data site - O, Oklahoma City; P,
Pittsburgh
                                               78

-------
APPENDIX  E
COMPLEX TERRAIN MODELING RESULTS:





      LISTING OF REGULATORY DESIGN CONCENTRATIONS (UG/M3)
                           79

-------
Appendix E.  Complex terrain results.  The ISC3 concentrations are in ug/m3 and the last 3
columns are concentration ratios.
PHLLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MCLLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MALLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
CCLLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
88.349
87.987
78.591
73.306
24.103
18.047
2.703

ISC3
51.612
50.919
48.584
41.564
12.846
10.534
1.348

ISC3
43.626
43.351
37.516
34.43
9.473
8.015
0.564

ISC3
48.416
48.416
42.68
32.291
9.707
7.256
0.48
02222AER/
ISC3
0.39
0.32
0.20
0.19
0.18
0.20
0.19
02222AER/
ISC3
0.51
0.41
0.21
0.24
0.16
0.17
0.17
02222AER/
ISC3
0.26
0.25
0.11
0.10
0.11
0.12
0.16
02222AER/
ISC3
0.30
0.24
0.20
0.19
0.17
0.16
0.20
02222AER/
CTDM+
0.67
0.55
0.54
0.48
0.38
0.44
0.41
02222AER/
CTDM+
1.04
0.97
0.63
1.08
0.53
0.63
0.52
02222AER/
CTDM+
0.41
0.56
0.47
0.44
0.79
0.74
0.62
02222AER/
CTDM+
1.23
1.22
1.77
1.78
2.06
1.57
1.18
02222AER/
99351 AER
0.87
0.74
0.61
0.61
0.75
0.89
0.89
02222AER/
99351 AER
0.85
0.96
0.87
0.93
0.91
1.01
1.03
02222AER/
99351 AER
0.91
0.91
0.92
0.91
0.99
0.99
1.09
02222AER/
99351 AER
1.00
1.00
1.01
1.05
0.99
1.02
1.18
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                       8 0

-------
Appendix E.  Complex terrain results.  The ISC3 concentrations are in ug/m3 and the last 3
columns are concentration ratios.
PHLLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MCLLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MALLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24 Hour H2H
Annual
CCLLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
8.293
7.971
4.908
3.841
1.122
0.739
0.067

ISC3
8.024
7.649
5.108
4.486
1.133
0.979
0.109

ISC3
7.78
7.662
4.03
3.327
0.842
0.57
0.038

ISC3
8.073
7.663
4.604
3.016
0.805
0.762
0.037
02222AER/
ISC3
0.40
0.33
0.25
0.27
0.18
0.26
0.27
02222AER/
ISC3
0.47
0.37
0.26
0.24
0.32
0.32
0.33
02222AER/
ISC3
0.19
0.14
0.17
0.17
0.20
0.20
0.21
02222AER/
ISC3
0.22
0.21
0.22
0.22
0.17
0.12
0.24
02222AER/
CTDM+
1.17
1.13
0.87
0.84
0.80
0.96
1.20
02222AER/
CTDM+
1.23
1.19
0.73
0.84
1.00
0.88
0.88
02222AER/
CTDM+
0.74
0.59
1.02
0.92
1.61
1.33
0.89
02222AER/
CTDM+
1.22
1.40
1.50
1.68
1.30
1.07
1.13
02222AER/
99351 AER
0.94
0.81
0.76
0.90
0.80
0.91
1.00
02222AER/
99351 AER
0.90
0.92
0.90
0.95
1.01
0.98
1.01
02222AER/
99351 AER
0.99
1.02
1.00
1.00
1.00
1.01
1.19
02222AER/
99351 AER
1.00
1.10
1.00
1.01
1.01
1.12
1.36
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                       81

-------
Appendix E.  Complex terrain results.  The ISC3 concentrations are in ug/m3 and the last 3
columns are concentration ratios.
PHLHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MCLHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MALHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
CCLHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
18.709
18.672
17.274
16.654
5.399
4.246
0.581

ISC3
9.41
9.385
9.208
9.164
3.195
2.775
0.249

ISC3
5.684
5.684
5.534
4.72
1.284
0.836
0.052

ISC3
0.878
0.869
0.788
0.779
0.304
0.143
0.01
02222AER/
ISC3
0.31
0.27
0.20
0.17
0.13
0.15
0.11
02222AER/
ISC3
0.38
0.28
0.23
0.20
0.12
0.13
0.10
02222AER/
ISC3
0.22
0.19
0.10
0.11
0.08
0.09
0.13
02222AER/
ISC3
0.38
0.32
0.15
0.13
0.09
0.18
0.20
02222AER/
CTDM+
0.66
0.66
0.77
0.75
0.70
0.80
0.97
02222AER/
CTDM+
0.55
0.58
0.97
1.10
0.89
0.95
0.75
02222AER/
CTDM+
0.34
0.30
0.23
0.39
0.31
0.28
0.39
02222AER/
CTDM+
0.22
0.22
0.14
0.18
0.18
0.34
0.67
02222AER/
99351 AER
1.01
1.01
1.00
1.00
1.00
1.01
1.03
02222AER/
99351 AER
1.20
1.02
1.01
1.01
1.02
1.02
1.05
02222AER/
99351 AER
0.97
1.00
1.00
1.01
0.98
0.98
1.00
02222AER/
99351 AER
0.71
0.85
0.76
0.95
0.87
1.29
1.60
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                       8 2

-------
Appendix E.  Complex terrain results.  The ISC3 concentrations are in ug/m3 and the last 3
columns are concentration ratios.
PHLHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24 Hour H2H
Annual
MCLHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MALHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
CCLHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
2.309
2.309
2.037
1.867
0.445
0.36
0.026

ISC3
1.547
1.547
1.477
1.164
0.408
0.284
0.028

ISC3
1.381
1.381
1.321
0.908
0.307
0.207
0.012

ISC3
0.27
0.262
0.219
0.217
0.088
0.044
0.003
02222AER/
ISC3
0.35
0.27
0.16
0.16
0.16
0.14
0.12
02222AER/
ISC3
0.54
0.29
0.20
0.19
0.12
0.12
0.14
02222AER/
ISC3
0.22
0.17
0.08
0.11
0.07
0.10
0.25
02222AER/
ISC3
0.79
0.73
0.41
0.37
0.16
0.27
0.33
02222AER/
CTDM+
0.45
0.39
0.51
0.57
0.52
0.63
0.75
02222AER/
CTDM+
0.60
0.52
0.66
0.59
0.78
0.62
0.80
02222AER/
CTDM+
0.31
0.34
0.31
0.34
0.42
0.40
0.75
02222AER/
CTDM+
0.35
0.47
0.39
0.41
0.45
0.48
1.00
02222AER/
99351 AER
1.23
1.09
1.01
1.39
1.05
1.22
1.17
02222AER/
99351 AER
1.80
1.18
1.12
1.23
1.07
1.01
1.12
02222AER/
99351 AER
1.00
0.91
0.84
1.01
0.84
0.90
1.20
02222AER/
99351 AER
0.86
1.00
0.97
0.95
0.85
0.97
1.46
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                       8 3

-------
Appendix E.  Complex terrain results.  The ISC3 concentrations are in ug/m3 and the last 3
columns are concentration ratios.
PHHLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MCHLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MAHLC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
PHHLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
49.86
49.715
28.962
26.97
9.578
8.275
1.022

ISC3
28.886
28.812
19.196
16.396
4.287
3.845
0.579

ISC3
15.289
15.289
8.119
6.722
2.004
1.815
0.128

ISC3
7.064
6.271
2.727
2.338
0.64
0.44
0.037
02222AER/
ISC3
0.35
0.26
0.31
0.28
0.30
0.27
0.25
02222AER/
ISC3
0.41
0.36
0.32
0.33
0.33
0.31
0.17
02222AER/
ISC3
0.43
0.28
0.36
0.28
0.25
0.16
0.19
02222AER/
ISC3
0.33
0.33
0.32
0.30
0.20
0.27
0.19
02222AER/
CTDM+
0.90
0.79
0.96
1.03
0.94
1.08
0.82
02222AER/
CTDM+
0.86
2.13
1.33
1.76
1.63
1.49
0.92
02222AER/
CTDM+
0.67
0.78
0.74
0.88
1.02
0.71
0.71
02222AER/
CTDM+
0.80
0.82
0.85
0.78
0.78
0.90
0.88
02222AER/
99351 AER
1.00
1.00
1.00
1.00
1.00
1.00
1.01
02222AER/
99351 AER
1.00
1.00
1.00
1.00
1.01
1.00
1.01
02222AER/
99351 AER
1.00
1.00
1.00
1.00
1.00
1.00
1.01
02222AER/
99351 AER
1.00
1.00
1.00
1.00
1.01
1.00
1.18
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                       8 4

-------
Appendix E.  Complex terrain results.  The ISC3 concentrations are in ug/m3 and the last 3
columns are concentration ratios.
MCHLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MAHLF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
PHHHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MCHHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
6.74
6.254
4.232
2.611
0.61
0.566
0.058

ISC3
6.509
6.371
2.17
2.16
0.43
0.423
0.022

ISC3
8.606
8.578
6.329
6.16
2.065
1.635
0.205

ISC3
1.029
1.026
0.945
0.935
0.389
0.371
0.022
02222AER/
ISC3
0.26
0.28
0.18
0.22
0.19
0.16
0.12
02222AER/
ISC3
0.22
0.14
0.22
0.19
0.17
0.14
0.27
02222AER/
ISC3
0.23
0.14
0.11
0.08
0.18
0.09
0.10
02222AER/
ISC3
0.66
0.28
0.24
0.15
0.17
0.11
0.18
02222AER/
CTDM+
0.85
1.58
1.09
1.39
1.30
1.39
1.00
02222AER/
CTDM+
0.81
0.82
0.67
0.92
0.83
0.72
0.86
02222AER/
CTDM+
0.32
0.29
0.31
0.27
0.91
0.43
1.25
02222AER/
CTDM+
0.40
0.25
0.31
0.25
0.60
0.44
0.50
02222AER/
99351 AER
1.00
1.00
1.08
1.00
1.00
1.23
1.00
02222AER/
99351 AER
1.00
1.00
1.00
1.00
1.00
0.99
1.10
02222AER/
99351 AER
1.25
1.00
1.07
0.99
1.00
1.00
0.99
02222AER/
99351 AER
0.95
0.91
0.95
0.98
1.00
0.95
1.02
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                       8 5

-------
Appendix E.    Complex  terrain results.   The ISC3 concentrations  are  in ug/m3
and the  last  3 columns are  concentration ratios.
MAHHC

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
PHHHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24 Hour H2H
Annual
MCHHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual
MAHHF

1 HourHIH
1 HourH2H
3 HourHIH
3HourH2H
24 HourHIH
24HourH2H
Annual

ISC3
0.659
0.643
0.614
0.525
0.124
0.077
0.006

ISC3
1.625
1.625
1.083
0.964
0.204
0.16
0.011

ISC3
0.274
0.272
0.235
0.233
0.095
0.053
0.007

ISC3
0.276
0.266
0.217
0.193
0.053
0.034
0.003
02222AER/
ISC3
0.38
0.32
0.19
0.16
0.27
0.38
0.33
02222AER/
ISC3
0.27
0.11
0.14
0.08
0.11
0.12
0.09
02222AER/
ISC3
0.59
0.51
0.33
0.30
0.21
0.25
0.29
02222AER/
ISC3
0.42
0.42
0.29
0.24
0.23
0.35
0.33
02222AER/
CTDM+
0.20
0.23
0.27
0.24
0.28
0.35
0.40
02222AER/
CTDM+
0.33
0.16
0.23
0.21
0.24
0.23
0.50
02222AER/
CTDM+
0.56
0.71
0.63
0.85
1.11
0.72
1.00
02222AER/
CTDM+
0.69
0.86
0.78
0.73
0.48
0.67
0.50
02222AER/
99351 AER
1.00
0.93
0.94
0.93
0.97
0.96
0.97
02222AER/
99351 AER
2.66
1.31
1.71
0.98
1.50
1.56
0.99
02222AER/
99351 AER
0.93
0.82
0.94
0.96
1.01
1.02
1.22
02222AER/
99351 AER
0.75
0.88
0.74
0.85
0.85
1.00
0.66
The code for each scenario evaluated is as follows: The first 2 letters refer to the name of the hill: PH = Piedmont,
MC = Montour Crosswind, MA = Montour Alongwind, CC = Cinder Cone Butte. The 3rd letter refers to the stack
height: L = 30 meters, H = 150m. The fourth letter refers to the buoyancy: L = low, H = high. The fifth letter stands
for the distance between the source and the hill: C = close or 1 kilometer, F = far or 10 km.                     8 6

-------
                                         TECHNICAL REPORT DATA

                                    (Please read Instructions on reverse before completing)
  1. REPORT NO.
   EPA-454/R-03-002
                                                                            3. RECIPIENT'S ACCESSION NO.
  4. TITLE AND SUBTITLE
  Comparison of Regulatory Design Concentrations: AERMOD
  vs. ISCST3, CTDMPLUS, ISC-PRIME
                                                                            5. REORT DATE
                    June 2003
                                                                            6. PERFORMING ORGANIZATION CODE
  7. AUTHOR(S)
                                                                            8. PERFORMING ORGANIZATION REPORT NO.
  9. PERFORMING ORGANIZATION NAME AND ADDRESS
                                                                            10. PROGRAM ELEMENT NO.
  See below.
                                                                            11. CONTRACT/GRANT NO.
  12. SPONSORING AGENCY NAME AND ADDRESS

  U.S. Environmental Protection Agency

  Office of Air Quality Planning and Standards

  Emissions Monitoring and Analysis Division

  Research Triangle Park, NC 27711	
                                                                            13. TYPE OF REPORT AND PERIOD COVERED
                    Final technical report.
                    14. SPONSORING AGENCY CODE
                    EPA/200/04
  15. SUPPLEMENTARY NOTES
  16. ABSTRACT     This report is a consequence analysis for the implementation of a new air dispersion model.  A consequence
  analysis is designed to give the user community a sense of how regulatory design concentrations from the new model compare to
  those from an established model via a series of "representative" examples.  For this study, the new model is AERMOD with the
  PRIME algorithms.  The existing regulatory models used in this report are ISCST3, ISC-PRIME, and CTDMPLUS. This analysis
  shows, for an extensive number of typical source scenarios,  the effects of the new model as compared to the existing regulatory
  model which it replaces.  Overall, except for complex terrain, the models produce rather similar results, although there are
  significant differences for individual source types and settings. As expected in complex terrain, AERMOD typically produces
  concentrations that are lower than ISCST3.  The model evaluation results, where available,  support the differences seen between
  the models.
  17.
                                           KEY WORDS AND DOCUMENT ANALYSIS
                      DESCRIPTORS
                                                        b. IDENTIFIERS/OPEN ENDED TERMS
                                                                                                 c. COSATI Field/Group
  consequence analysis, regulatory design concentrations,
  AERMOD, ISCST3, CTDMPLUS, PRIM E,ISC-PRIM
Air Pollution models
  18. DISTRIBUTION STATEMENT
   Release Unlimited
                                                        19. SECURITY CLASS (Report)   Unclassified
                                         21. No of pages  91
20. SECURITY CLASS (Page)   Unclassified
EPA Form 2220-1 (Rev. 4-77)   PREVIOUS EDITION IS OBSOLETE

-------
United States                       Office of Air Quality Planning and Standards             Publication No. EPA-454/R-03-002
Environmental Protection            Emissions Monitoring and Analysis Division                                    June 2003
 Agency                             Research Triangle Park, NC

-------