Example Area of Influence AOI Analysis for VISTAS Class I Areas

1
Example Area of Influence (AOI) Analysis for
VISTAS Class I Areas

• Gerry Mansell and Ralph Morris, ENVIRON
• Joe Adlhoch and Cassie Archuleta,
• Air Resource Specialists
• Greg Stella, Alpine Geophysics
• VISTAS Joint Work Group Meeting
• Raleigh, North Carolina
• September 6, 2006

2
Overview of AOI Approach

• Calculate Back Trajectories (Particle Paths) from
  VISTAS Class I Areas using NOAA HYSPLIT Model for
  all days from 2000-2004 Baseline
• Use Back Trajectories for Worst 20 Days from
  5-year Baseline
• Calculate Residence Time of Back Trajectories by
  binning them in lat/long grid cells
• Map Residence Time from 1? x 1? lat/long grid to
  36 km modeling LCC projection grid
• Weight them by total extinction (Bext) and
  extinction due to SO4/NO3
• Normalize Residence Time analysis to maximum
  value for plots Examples for 4 sites SIPS,
  CHAS, OKEF, ROMA
• Superimpose on emissions inventory and decay with
  distance

3
Summary of 2018 Base F4 Uniform Rate of Progress
Assessment
.
.
Hercules Glade, MO
CHAS, OKEF, ROMA and SIPS selected for example
AOI analysis
.
.
.
.
.
.
.
.
Current and Proposed New IMPROVE Equations and
2000-2004 Baseline
4
(No Transcript)
5
Back Trajectory Modeling

• The NOAA HYSPLIT model was used to generate
  multiple daily back trajectories for 4 sites
  SIPS, CHAS, OKEF, and ROMA
• The EDAS (Eta Data Assimilation System)
  meteorological data set was used
• The HYSPLIT parameters selected were
• Duration 72 hrs (3 days)
• End times 0600, 1200, 1800, 2400 EST
• End heights 100m, 500m
• Vertical motion option data

6
Example Back Trajectories, SIPS
7
Trajectory Processing (1)

• Trajectories were run for all days during the
  baseline period, 2000 2004
• Daily trajectories corresponding to the 20
  haziest days for each year were selected
• Individual hourly trajectory points were binned
  into 1 degree lat/long grid cells
• Individual trajectory points were weighted by a
  sites extinction data for the corresponding 20
  worst day in four ways
• By total extinction
• By ammonium sulfate extinction
• By ammonium nitrate extinction
• By organic mass extinction

8
Trajectory Processing (2)

• The concentration-weighted residence time field
  was reprojected to fit the 36km modeling domain
• Each grid cell was normalized to the maximum grid
  cell value and plotted
• Residence Time fields were overlayed on emissions
  with a distance decay factor
• 1/r distance decay to account for dispersion and
  deposition

9
(No Transcript)
10
Residence Time Analysis SIPS 100m AGL Bext
weighted
11
Residence Time Analysis SIPS 100m AGL Bext
weighted
12
Residence Time Analysis SIPS 500m AGL Bext
weighted
13
SIPS W20 Bext Weighted Residence Time Comparison
for 100m vs. 500m AGL Back Trajectory start height
100m
500m
14
Residence Time Analysis CHAS 100m AGL
15
CHAS W20 Bext Weighted Residence Time Comparison
for 100m vs. 500m AGL Back Trajectory start height
16
(No Transcript)
17
OKEF W20 Bext Weighted Residence Time Comparison
for 100m vs. 500m AGL Back Trajectory start height
18
(No Transcript)
19
ROMA W20 Bext Weighted Residence Time Comparison
for 100m vs. 500m AGL Back Trajectory start height
20
Next Steps

• Overlay Residence Time analysis on top of VISTAS
  2002 Typical Base G emissions
• SO2 and NOx Point NOx Low-Level Other
• Add decay with distance to account for dispersion
  and deposition
• 1/r, with cell containing monitor assumed to be ¼
  length from monitoring site
• VISTAS 2002 Typical Base G Emissions

21
Bext, SO2 Point Emissions and Distance Weighted
Residence Time Analysis for W20 Days at SIPS
Interim Non-VISTAS Emissions