Analysis of CAMx and UAM-V Modeling in an Indirect Technique to Assess Contribution to Ozone Exceedances in the Eastern United States

Page created by Laurie Castillo
 
CONTINUE READING
Analysis of CAMx and UAM-V Modeling in an Indirect
Technique to Assess Contribution to Ozone Exceedances
in the Eastern United States
Paper presented at the 93rd Annual Meeting of the Air and Waste Management
Association, Salt Lake City, UT, June 2000.

Saravanan Arunachalam
Environmental Programs
MCNC–North Carolina Supercomputing Center
P.O. Box 12889, 3021 Cornwallis Road
Research Triangle Park, NC 27709-2889

ABSTRACT
Regional-scale photochemical air quality simulation models (PAQSMs) are increasingly
being used to address the complex nature of the tropospheric ozone problem. The U.S.
Environmental Protection Agency recently promulgated a NOx State Implementation Plan
(SIP) Call to address the problem of transport of ozone and its precursors across the
eastern United States. The NOx SIP Call includes emissions budgets that are designed to
eliminate specific amounts of NOx, one of the precursors to ozone. In developing the NOx
SIP Call, EPA performed numerous air quality simulations using the Variable-grid Urban
Airshed Model(UAM-V) and the Comprehensive Air Quality Model with
Extensions(CAMx) to identify the amounts of NOx emissions that contribute significantly
to nonattainment of the National Ambient Air Quality Standard (NAAQS) for ozone in
downwind areas. The analysis described here, which focuses on the “zero-out” technique,
and the “source-apportionment” technique, specifically compares the quantification of
state-by-state contributions to ozone exceedances using two modeling systems. The zero-
out method estimates the downwind impacts on ozone nonattainment by comparing the
model predictions from a base case simulation to the predictions from a simulation in
which all anthropogenic emissions are removed from a specific state. EPA had applied
this “zero-out” technique for three different states using both UAM-V and CAMx for
four different episodes in which high levels of ozone were observed in various parts of
the eastern United States. To assess the results from each of these simulations, we
identified various metrics designed to provide information on the fundamental factors for
evaluating whether emissions in an upwind state can contribute to downwind
nonattainment. To further corroborate these results, the outputs from CAMx simulations
with source apportionment were also analyzed and compared. The results from this
analysis should assist the air quality modeling community in the scientific and regulatory
world in comparing the performance of two different modeling systems, and identify
advantages, or even potential limitations, in the choice of a given PAQSM for addressing
regional-scale air quality issues.

                                                                                             1
INTRODUCTION
The problem of regional transport of ozone and its precursors is being studied for more
than two decades1 . The scales of such transport could be inter-city, inter-state, or even
inter-regional. The movement of elevated regional ozone concentrations may reflect
impact of large-scale weather systems. Exceedances of the ozone standard often occur
when the ozone precursor emissions from major metropolitan areas interact with the
elevated regional concentrations resulting from transport. Emission areas that are
potentially responsible for the violations of the air quality standard at downwind Areas of
Violations (AOV) are termed Areas of Influence (AOI). Various tools that could help in
establishing the domains of AOIs and AOVs are: Analysis of Monitor Data, Air Parcel
Trajectory Modeling and Photochemical Air Quality Modeling. Though each have their
own advantages and limitations, grid-based Photochemical Air Quality Simulation
Modeling (PAQSM) has the potential to be a most comprehensive tool for source
attribution processes. Of course, this requires that multi-scale grid is designed in such a
way that it is able to “catch” the transport phenomena of importance; unfortunately, this
is not the case yet with most existing models. Major metropolitan areas like New York,
Chicago, and Atlanta are often dominated by overwhelming transport of ozone and its
precursors from upwind sources. While local emissions within the domain are expected
to be more important with respect to 1-h concentrations of ozone, larger areas and longer
distance scales can affect the system when longer averaging times are considered for
reducing the regional pollutant levels.

In developing the NOx State Implementation Plan (SIP) Call2 , the U.S. EPA performed
various air quality analyses3 to support the multifactor approach used to identify amounts
of NOx emissions that contribute significantly to nonattainment in downwind areas.
These analyses include subregional and state-by-state modeling intended to

1. quantify the air quality contributions from emissions in upwind states to both 1-h and
   8-h nonattainment in downwind areas, and
2. determine whether these contributions are significant.

The modeling for assessing contributions consisted of state-by-state zero-out modeling
using the Variable-grid Urban Airshed Model (UAM-V)4 and Comprehensive Air Quality
Model with Extensions (CAMx)5 , and state-by-state source apportionment modeling
using the CAMx Anthropogenic Precursor Culpability Assessment (APCA)6 technique.
For assessing regional strategies, EPA performed a series of additional UAM-V model
runs. This study focuses on the state-by-state zero-out modeling performed with UAM-V
and CAMx, and source apportionment modeling performed with CAMx. The motivation
for this study is to see how two very different modeling systems behave in assessing air
quality contributions when the models are subject to similar emissions perturbations.

APPROACH
To assess the effect of transport of ozone and its precursors, various geographical
sensitivity simulations have been designed and performed: source-apportionment
technique, rollout simulations, inverse-rollout7 simulations, zero-out technique, direct

                                                                                              2
emissions reductions approach, etc. While all of them have their limitations, it is
important to understand the differences between them. The zero-out modeling is quite
different from the source apportionment modeling for evaluating upwind contributions to
downwind receptors. This is because the zero-out modeling is an indirect technique to
quantify the contributions by the differences of two model-runs (i.e. base case versus
zero-out), while the source apportionment technique directly estimates the contributions
of upwind sources to receptor areas from a single simulation. This might lead to
differences in both magnitude and frequency of contributions for individual upwind-to-
downwind linkages. From the various modeling runs, EPA found that the source
apportionment modeling tended to show greater contribution magnitudes than the zero-
out modeling. While both modeling techniques have their usefulness, they also have their
limitations and there is no strong technical evidence to show that one is clearly superior
to the other.

The two models used by EPA, UAM-V (Version 1.24) and CAMx (Version 1.13), are
three-dimensional photochemical air quality simulation models (PAQSMs) designed for
integrated assessment of photochemical air pollution over regional and urban scales. Both
the models simulate the emission, dispersion, and removal of inert and chemically
reactive pollutants in the lower troposphere by solving a set of species mass continuity
equations that represent the relevant processes for each chemical species on a system of
nested three-dimensional grids. While there are many similarities between the two
modeling systems, there are major differences between them: the treatment of plume-in-
grid, the numerical methods for solving chemistry and vertical diffusion, and the
calculations of dry deposition. A key feature of CAMx is its ability to perform source
apportionment to assess the contributions of multiple source regions and source
categories to specific receptor areas from a single simulation. While giving a clear picture
of the likely distribution of ozone and its precursors by source category, this also
indicates whether the ozone at the selected time and location would be more likely
respond to upwind NOx or VOC controls. Both models are suitable for evaluating the air
quality effects of emission control scenarios because they account for spatial and
temporal variations as well as differences in the reactivity of emissions. Traditionally, the
modeling applications involve first replicating historical ozone episodes to establish a
base year simulation. Once the model results have been evaluated and the model
determined to perform within prescribed levels, the same base year meteorological inputs
for each episode are combined with projected emissions inventories to simulate a
projected future base case and alternative emission control strategies as well.

In support of the NOx SIP Call, for purposes of evaluating air quality impacts, EPA
modeled all four episodes that were used by the Ozone Transport and Assessment Group
(OTAG)8 . The OTAG effort is one of the most comprehensive three-dimensional gridded
regional ozone modeling projects to date, and the results from OTAG are likely to be
used as boundary conditions for future urban-scale ozone modeling studies. These
episodes correspond to periods of elevated ozone concentrations in the eastern United
States. The first period was July 1-15, 1988, which had elevated ozone in the Midwest,
Northeast, and Southeastern United States. The second period was July 13-21, 1991,
which had high ozone in the Midwest and Northeastern areas. The third period was July

                                                                                                3
20-30, 1993, which had elevated ozone levels throughout the Southeast and in portions of
the Southwest, Midwest, and Northeast. The fourth period was July 7-18, 1995, which
had elevated ozone in the Midwest, Northeast, and Southeast.

Each of the above episodes included a 2- or 3-day “ramp-up” period to initialize the
model (the first three simulation days for the 1988, 1991, and 1995 simulations and the
first two days for the 1993 simulations). Predictions from these the ramp-up periods were
excluded from the evaluation and analyses.

After projecting the emissions from these episodic years to the future year 2007 and
establishing a 2007 base case for each episode, EPA performed the state-by-state zero-out
modeling with UAM-V for the following states: Alabama, Georgia, Illinois, Indiana,
Kentucky, Massachusetts, Michigan, Missouri, North Carolina, Ohio, South Carolina,
Tennessee, Virginia, West Virginia, and Wisconsin, for the combined states of Illinois
and Wisconsin, and the combined Southeastern states of Alabama, Georgia, North
Carolina, South Carolina, and Tennessee. CAMx, on the other hand, was used to perform
zero-out modeling only for the states of Illinois, Michigan, and North Carolina.
Therefore, this study focuses on analyses of the zero-out CAMx and UAM-V simulations
for Illinois, Michigan, and North Carolina only; these are referred to as zeroIL, zeroMI,
and zeroNC.

CAMx simulations with “source-apportionment” technique were performed by EPA for
the 2007 base year for all the four episodes. The analyses of model comparisons from
UAM-V and CAMx zero-out simulations are then corroborated with outputs from CAMx
simulations performed with source apportionment.

MODEL SET-UP
The modeling domain used in all these simulations was the same as that used for OTAG.
It includes portions or all of 37 states in the eastern United States, the District of
Columbia, and part of Canada. EPA used the same OTAG configuration here for initial
and boundary conditions, meteorological inputs, and grid configuration in both the
horizontal and vertical resolutions. Horizontally, the modeling domain consisted of a
coarse grid with 64 x 63 cells at a 36 km grid resolution and a fine grid with 137 x 110
cells at a 12 km grid resolution. The meteorological inputs for the 1991 and 1995
episodes were derived from applications of the Regional Atmospheric Modeling System–
Version 3a (RAMS)9 . For the 1988 and 1993 episodes, the meteorological inputs were
derived from the Systems Application International Mesoscale Meteorological model
(SAIMM)10 .

The base year emissions inventory was a combination of 1995 and 1996 data, intended to
be as representative of the 1994-1996 time period as possible. To produce the 2007 base
case inventory, the 1995/1996 base year inventory was projected to 2007 and certain
controls were applied. The 2007 base case emissions reflect Clean Air Act mandated
controls as well as certain Federal measures that EPA has promulgated. Across the entire
modeling domain, the 2007 base case NOx emissions and VOC emissions are 1% and
19% less than the 1995/1996 base year emissions, respectively.

                                                                                            4
MODEL PERFORMANCE EVALUATION
As part of the OTAG study, an objective evaluation of UAM-V model performance11 was
conducted for each of the four episodes to determine how well the modeling system
represented regional ozone concentration levels. Various statistical metrics comparing
predicted ozone to observed ground-level ozone concentrations were computed. While
the results indicated generally good agreement between observed and predicted values, a
few relatively minor concerns were found. One of them was that the model tended to
underestimate ozone aloft, especially overnight into the early morning hours. Thus, the
contribution of upwind source regions to ozone levels in downwind areas may actually be
greater than that estimated by the model.

Several modeling comparative studies have evaluated the model performance of CAMx
and compared it with that of UAM-V. Two such studies independently focused on the
1991 and 1995 episodes.12,13 In general, these comparative studies found that the
performance of CAMx for ozone was similar to the performance of UAM-V. CAMx,
with a lower mean normalized error, performed slightly better on average for 1995 than
did UAM-V. Comparisons with aircraft data14,15 revealed that both models may
somewhat underestimate the average ozone concentrations aloft. In support of the NOx
SIP Call, EPA also compared the two models in terms of the predictions of ozone and its
precursors for the 1995 episode alone. In addition, the comparisons by EPA found that
both models tend to underpredict the mean ozone concentration above 80 and 120 ppb,
and that CAMx tends to predict higher ozone domain-wide than UAM-V does (especially
in the high-ozone hours) by 2-5 ppb.

The model performance evaluations and comparisons discussed above show that the two
models are technically equivalent based on various statistical measures.

ANALYSES
In using the two models to assess the relative contribution from zeroing out individual
states, we focused on the 1-h metrics alone for nine nonattainment regions. The major
urban areas included in these regions are Atlanta, Baltimore, New York, Philadelphia,
Chicago, Milwaukee, Western Massachusetts (Springfield), Greater Connecticut
(Hartford), and Washington, D.C.

We analyzed the model-predicted ozone concentrations using the following three metrics,
calculated by EPA for the 2007 base case and then the three zero-out scenarios from both
models.

•   Metric 1: Maximum 1-h Concentration, which is the maximum concentration
    predicted by the model for the entire episode
•   Metric 2: Exceedances, which is the number of 12-km grid cells with predictions of
    1-h ozone =125 ppb
•   Metric 3: Ozone reduced (ppb), which is the magnitude and frequency of the “ppb”
    contributions in ozone concentrations =125 ppb for the 1-h NAAQS. These ozone

                                                                                           5
impacts were quantified by calculating the difference in predicted ozone between
    each scenario and the 2007 base case. Contributions are found in terms of the
    frequency of impacts from the zero-out simulation on ozone concentrations in each
    downwind area for six concentration ranges: >2 to 5 ppb, >5 to 10 ppb, >10 to
    15 ppb, >15 to 20 ppb, >20 to 25 ppb, and >25 ppb.

These three metrics are each calculated twice, first for all the grid cells within a state and
then for just the designated 1-h-exceedance grid cells within each state. The latter is
intended to assess the impacts seen only in the nonattainment counties, as opposed to the
entire state.

To assess the outputs from the CAMx simulations with source apportionment, the
following metrics alone calculated by EPA are used in the analyses:

•   Metric 4: Total ppb 1-h ozone, contributed by a single source state to various
    downwind receptor states
•    Metric 5: Total ppb 1-h ozone, in a single downwind receptor state contributed by
    various upwind source states
•   Metric 6: Total ppb 8-h ozone, contributed by a single source state to various
    downwind receptor states
•   Metric 7: Total ppb 8-h ozone, in a single downwind receptor state contributed by
    various upwind source states

These four metrics are averaged across all the four episodes mentioned earlier. It should
be noted that while the calculation of Metrics 4 and 5 was done only in the nonattainment
counties of each state, Metrics 6 and 7 were calculated in all grid-cells in each state.

RESULTS
The metrics as defined above are presented for both CAMx and UAM-V simulations for
the 2007 base case, and then for the three zero-out simulations for all four episodes.

Metric 1:
Figure 1 shows in the form of bar charts, the episodic maximum 1-h ozone concentration
predicted within the nonattainment counties in each state for each of the four episodes
from CAMx and UAM-V simulations. The top portion of Figure 3 shows the maximum
1-h ozone concentration from all four episodes, calculated within nonattainment counties.

The only significant contribution to this metric seems to be from the zeroIL simulation,
and the effects are seen in the maxima predicted within the states of both IL and WI. The
decreases range from 10% with respect to the base case in the 1998 episode to 37% in the
case of the 1995 episode. The other two cases, zeroMI and zeroNC, seem to have almost
negligible impacts on all the states considered. Both CAMx and UAM-V show similar
behavior for all four episodes.

                                                                                                 6
Metric 2:
Figure 2 shows metric 2 predicted within nonattainment counties in each state for each of
the four episodes from CAMx and UAM-V simulations. The bottom portion of Figure 3
shows the same metric calculated cumulatively for all four episodes.
The zeroIL case seems to show almost 100 % reduction in the number of exceedances in
IL itself in 1988 and 1995 from both CAMx and UAM-V. While most of the other
reductions are comparable between the two models, the zeroMI case seems to show a
100% reduction in NY (CAMx 1993) and the zeroNC case seems to show a 100%
reduction in PA (UAM-V 1991) cases. These reductions are not shown by the other
models in the corresponding cases.

Metric 3:
This metric is presented in Figures 4-7 as stacked bar charts representing the various bins
of the “ppb” contributions to ozone concentrations =125 ppb at different threshold ranges
for each state.

Figures 4-6 show the contributions in the nonattainment counties of each state from
zeroIL, zeroMI, and zeroNC, from both CAMx and UAM-V modeling for each of the
four episodes. Figure 4 shows that the zeroIL case for the 1988 episode, while exhibiting
impacts in just the 2-5 ppb range in both CT and NY, has impacts in the 20-25 ppb range
in IL from both CAMx and UAM-V simulations. In the case of the zeroMI case (Figure
5), both CAMx and UAM-V seem to show reductions in the states of MD and NY in the
2-5 ppb range. The zeroNC case in Figure 6 shows impacts in MD in ranges up to 10-
15 ppb as predicted by CAMx for the 1988 episode. Marginal impacts are predicted by
both models for the other three episodes

Figure 7 shows the same metric plotted for the composite of the four episodes for zeroIL,
zeroMI, and zeroNC, calculated in the nonattainment counties of each state. While these
plots do not show very large bias from either model, it should be noted that the
percentage reductions predicted by CAMx are somewhat higher than those predicted by
UAM-V in most cases.

Metric 4 and 6:
The percent contribution of ozone from each of the three source states – IL, MI, and NC
– to downwind states are ranked, and the top 10 states affected are shown as pie charts in
Figure 8 (1-h ozone on the left and 8-h ozone on the right).

Metric 5 and 7:
The percent contribution of ozone in six downwind receptor states – CT, DE, MD, NJ,
NY, and PA – from various source states are ranked, and the top 10 contributing states
are shown as pie charts in Figure 9 (for 1-h ozone) and Figure 10 (for 8-h ozone).

There are some very interesting results seen from these pie charts. While sources in IL
contribute 10-12% of 1-h ozone in both receptors NJ and NY, of all the 1-h ozone in NJ

                                                                                              7
and NY, sources in IL constitute less than 5 %. Also, while sources in MI contribute
around 17% to 1-h ozone in NJ and NY, sources in MI constitute only less than 2% in the
NJ and NY receptor charts. However, the zeroIL and zeroMI CAMx simulations in
Figure 2 for the 1993 and 1995 episodes seem to show almost 50 to 65% reductions in
the number of 1-h exceedances in NJ and NY.

Similarly, sources in IL, MI and NC contribute anywhere between 14.6 to 32.5% to 1-h
ozone in MD. However, of all the 1-h ozone in MD, these three combined contribute by
less than 1%. From Figure 2, there was only less than 10% reduction in the number of 1-
h exceedances from any of the three zero-out cases from both models.

                                                                                          8
Figure 1. Episodic maximum 1-h ozone concentration (ppb) predicted within non-
attainment counties from CAMx (left) and UAM-V (right) for 1988, 1991, 1993, and
1995 episodes (top to bottom).

                                                                                   9
Figure 2. Number of exceedances of 1-h ozone concentration (ppb) predicted within non-
attainment counties from CAMx (left) and UAM-V (right) for 1988, 1991, 1993, and
1995 episodes (top to bottom).

                                                                                         10
Figure 3. Maximum 1-h ozone concentration (top) and total number of exceedances of 1-
h ozone concentration (ppb) (bottom) predicted within non-attainment counties from all
four episodes from CAMx (left) and UAM-V (right)

                                                                                         11
Figure 4. Magnitude and frequency of “ppb” contributions to ozone concentrations
=125 ppb in nonattainment counties due to Illinois zero-out from CAMx (left) and
UAM-V (right) for 1988, 1991, 1993, and 1995 episodes (top to bottom).

                                                                                   12
Figure 5. Magnitude and frequency of “ppb” contributions to ozone concentrations
=125 ppb in nonattainment counties due to Michigan zero-out from CAMx (left) and
UAM-V (right) for 1988, 1991, 1993, and 1995 episodes (top to bottom).

                                                                                   13
Figure 6. Magnitude and frequency of “ppb” contributions to ozone concentrations
=125 ppb in nonattainment counties due to North Carolina zero-out from CAMx (left)
and UAM-V (right) for 1988, 1991, 1993, and 1995 episodes (top to bottom).

                                                                                     14
Figure 7. Magnitude and frequency of “ppb” contributions to ozone concentrations
=125 ppb in nonattainment counties for composite of the four episodes from CAMx (left)
and UAM-V (right) for Illinois, Michigan, and North Carolina zero-out simulations (top
to bottom).

                                                                                         15
Figure 8. Percent of the ozone contribution from each of the source states (IL, MI, and
NC) to downwind states for 1-h ozone concentrations (left) and 8-h ozone concentrations
(right) from CAMx Source Apportionment modeling (composite of all the four episodes).

                                                                                          16
Figure 9. Percent of the 1-h ozone contribution in receptor states (CT, DE, MD, NJ, NY,
and PA) from various source states from CAMx Source Apportionment modeling
(composite of all the four episodes).

                                                                                          17
Figure 10. Percent of the 8-h ozone contribution in receptor states (CT, DE, MD, NJ,
NY, and PA) from various source states from CAMx Source Apportionment modeling
(composite of all the four episodes).

                                                                                       18
CONCLUSIONS
We used the zero-out technique to analyze state-by-state contributions to ozone
exceedances, based on outputs from two different modeling systems, CAMx and
UAM-V. The results from CAMx simulations with source apportionment were also used
to compare with findings from the zero-out simulations. Various metrics calculated from
zero-out simulations for the states of Illinois, Michigan, and North Carolina for four
different episodes were used to assess downwind impacts on ozone nonattainment. The
model performance evaluations that were performed as part of earlier studies for both
CAMx and UAM-V indicated that both models are technically equivalent, and that both
models may underestimate the average ozone concentrations aloft. The results of the
zero-out simulations analyzed for all the grid cells seem to show similar patterns to those
analyzed for the nonattainment counties alone within each state, for both CAMx and
UAM-V. Both models predicted very small changes in the peak 1-h ozone concentration
within each downwind area for all cases, except in IL and WI due to zeroIL. In most
cases, the other two metrics behaved fairly similarly in assessing upwind contribution to
downwind ozone nonattainment. Comparisons with assessment from source-
apportionment seemed to show somewhat different contributions, especially to some
states in the Northeastern U.S. Given the nature of the two techniques used here (direct
versus indirect approach to assessing ozone nonattainment), some of these differences
can be placed in perspective. However, the results from this study still present an
opportunity to compare the behavior of two different modeling systems with similar
inputs in addressing regional air quality issues. The results presented here are limited to
the extent that the zero-out simulations were available only for three states for both
modeling systems. Additional zero-out simulations performed for other states might help
to validate the findings discussed here.

ACKNOWLEDGMENTS
The author gratefully acknowledges the U.S. EPA’s Office of Air Quality Planning and
Standards, which provided the results from all the model simulations that were used for
these analyses, and the immense contributions from the various modeling centers that
participated in the OTAG study. The author would also like to acknowledge the editorial
assistance provided by Jeanne Eichinger in reviewing this manuscript and providing
useful comments.

DISCLAIMER
The opinions presented in this paper are those of the author alone and do not reflect those
of any other agency or organization. Mention of model names does not constitute
endorsement or recommendation for use in any manner.

REFERENCES
1. Wolff, G.T.; Lioy, P.J. “Development of an ozone river associated with synoptic scale
   episodes in the Eastern United States”, Environ. Sci. and Technol. 1980, 14(10),
   1257-1260.

                                                                                              19
2. U.S. Environmental Protection Agency. Air Quality Modeling Technical Support
   Document for the NOx SIP Call, Office of Air and Radiation: Research Triangle Park,
   NC, September 1998.
3. U.S. Environmental Protection Agency Regional Modeling Center Home Page.
   http://www.epa.gov/ttn/scram/regmodcenter/t28.htm (accessed December 1999).
4. Systems Applications International. User’s Guide to the Variable-grid Urban Airshed
   Model (UAM-V), San Rafael, CA, April 1995; SYSAPP-95/027.
5. Environ International Co. User’s Guide to the Comprehensive Air Quality Model with
   Extensions (CAMx) Version 1.1, Novato, CA, April 1997.
6. Environ International Co. User’s Guide to the Ozone Tool: Ozone Source
   Apportionment Technology. Novato, CA, March 1997.
7. Arunachalam, S.; Georgopoulos, P.G.; Sistla, G. A study of transport of ozone and
   precursor emissions for estimating Areas of Influence (AOI) in the Eastern United
   States, In Proceedings of the A&WMA Emissions Inventory Conference: Planning for
   the Future, Research Triangle Park, NC, October 1997.
8. Ozone Transport and Assessment Group. Ozone Transport and Assessment Group
   Technical Support Document, June 1997.
9. Pielke, R.A.; Cotton, W.R.; Walko, R.L.; Tremback, C.J.; Lyons, W.A.; Grasso, L.D.;
   Nicholls, M.E.; Moran, M.D.; Wesley, D.A.; Lee, T.J.; Copeland, J.H. Meteorol.
   Atmos. Phys. 1992, 49, 69-91.
10. Systems Applications International. User’s Guide to the SAI Mesoscale
    Meteorological Model (SAIMM), San Rafael, CA, 1995.
11. Lurman, F.W.; Kumar, N.; Londergan, R.; Moore, G. Evaluation of the UAM-V
    Model Performance in the Northeast Region for OTAG Episodes, Sonoma
    Technology, Inc.: Santa Rosa, CA, March 1997; STI-996133-1716-WD2.1.
12. Sonoma Technology, Inc. Comparison of CAMx and UAM-V Model Performance for
    Two Ozone Episodes in the Eastern United States, Santa Rosa, CA, 1997; STI-
    996203-1733-FR.
13. Alpine Geophysics. Modeled Effects of Indiana Point-Source NOx Emissions
    Reductions on Local and Regional 1-h and 8-h Ground Level Ozone Concentrations
    in 1995 and 2007 Using Two OTAG Episodes, Covington, KY, January 1998.
14. Earth Tech, Inc. Assessment of UAM-V Model Performance in Northeast Region for
    OTAG Episodes, Concord, MA, March 1997.
15. Tesche, T.W. Evaluation of the MM5, SAQM, UAM-IV and UAM-V Models over
    the Northeast US for four ozone episodes using routine and intensive NARSTO-NE
    and LMOS Aerometric Datasets. In Proceedings of the First NARSTO-Northeast
    Data Analysis Symposium and Workshop, Norfolk, VA, December 1996.

                                                                                         20
You can also read