-------
Local terrain elevations near the release point are estimated to
be 17.4 m above the base elevation of the hill coordinate system, so
the net release height corresponds to the 37.4 m height level on the
ridge.
Meteorological Information
Figure 44 contains plots of the 5-minute Hc and Fr values for
this hour. Hc decreases from 23 m in the beginning of the hour to
5 TO half-way through the hour, then rises to 34 m during the latter
part of the hour before decreasing to 10 m at the end of the hour.
The average of the 5-minute values over the hour is 18 m. Fr
gradually rises from 1.3 to 1.8 during the hour. The one hour average
value of Fr above Hc is 1.5, which is an indication that
stratification has only a weak effect on the flow over the top of the
hill.
The hourly averaged vector wind direction and wind speed at the
tracer release height are estimated to be 59.3° and 1.2 m/s,
respectively. Figure 45 shows the trend in wind speeds and directions
between 5 m and 40 m during the hour. There is little directional
wind shear for approximately two-thirds of the hour, but there is up
to 100° of directional wind shear during the remaining third of the
hour. Large directional wind shear between 5m and 40 m occurs in
combination with increasing Hc above 5 m. The flow direction is
more variable well below Hc than above and this is apparently due to
the blocking effect of the nearly two-dimensional ridge. The wind
speed time series varies inversely with the Hc time series. As the
wind speed increases during the first half of the hour, Hc
decreases; conversely, as the wind speed decreases during the second
half of the hour, Hc increases.
The vertical turbulent intensity values estimated at the tracer
release height vary from 8 to 27% and the av values estimated at
the tracer release height vary from 0.18 to 0.29 m/s during the hour.
The one hour value for aw estimated at the release height is
0.22 m/s (iz = 18%).
The vertical profiles of the hourly averaged wind direction, wind
speed, and temperature are shown in Figure 46. For the layer 10 m
above and below the release height, the wind speed shear is 0.5 m/s
and the directional shear is 32°. The wind data from the propeller
anemometers at 5 m and 40 m compare favorably with the sonic
anemometer wind data at the same levels. The wind direction profile
between 5 m and 40 m is approximately linear, so the MDA value for the
wind direction should be reasonable. The hourly average of the
5-minute propeller wind direction and wind speed values measured at
the tracer gas release height are 50° and 1.7 m/s which corroborate
the one-hour average MDA values.
In summary, the plume was released above Hc during the first
half of the hour and below Hc during the second half of the hour.
The hourly average value for Hc is 18 m which is close to the
release height. The flow is characterized by substantial wind
129
-------
u
Release Height
Time (Hour)
2.0T-
1.5-
o
5- i .0-
(Experiment 14, 10/26/82, 0300-0400 MDT).
130
-------
m IM
(s/ui) peodg PUIAA
3
o
I
o
0)
tn
*£ 3
r-f
lj (d
jf
O
fl •
<*-s
e H
0 Q
£ s
O
(D O
4J *^
n) o
•O 1
o
o o
••-I CO
G O
o
w -
CM
Oi 00
-*J ^.
D \o
C CM
•H \
e o
1 H
in
e
o
Ll
14J
•a
H ^ o e
60
IM IM
-------
t- :
O.
o
H
O
U
8 8
8 ? 8
8 8 S
S
JR R
S 8
CO
•o
o
•a
to
•o
a
o
00
o
— o
M (1)
\ ^
1 I
o
o
O
o
CO
O
•_ .
3 CM
O 00
(4-1 CM
O >»
O
14-1
o
- §
« -rM
O U
• r-t 0)
*•» ^
vo
-5T
Q
t-i
00
132
-------
meander, wind speeds varying from 0.7 to 2.2 m/s, and a hill Froude
number of 1.5.
Cjf.3§r Concentrations
The distribution of the observed hourly averaged CP^Br
concentrations over the surface of the ridge is shown in Figure 47.
The concentration pattern is consistent with the variable wind
directions observed during the hour. The largest concentrations are
found near the bottom-half of the sampler array at an elevation that
is less than Hc. The maximum observed concentration (117 ys/m^)
is found near the estimated hourly-averaged plume centerline. There
is a sharp decrease of observed tracer concentrations above the Hc
surface (1635 m).
Model Performance
Results from the HER (Flat) model are displayed in Figure 48.
Overall, the mean of the concentration estimates is in close agreement
with the mean of the observed concentrations (Co/Cp - 1.04,) but
the spatial correlation is poor (r^ = 0.31). The maximum observed
concentration (117 ps/nP) is more than a factor of two greater
than the maximum modeled concentration (52 ps/m3).
Estimates from the HBR (Terrain) model, shown in Figure 49, are
somewhat larger than the flat-terrain estimates. The maximum modeled
concentration (79 v>s/m3) is less than the maximum observed
concentration (117 ys/rn-*), although on average, the model is over
estimating (cfo/cl = 0.74). Because the distribution of the peak
observed concentrations in space is so unlike that estimated by the
model using one-hour average data, the temporal variations in the
meteorology (particularly iz) may be responsible for the poor model
performance in estimating the peak concentrations.
»
4.4.3 Experiment 6, Experiment-Hour 9 (0700-0800)
Release Description
The CF3Br tracer gas was released from position 203 for the
entire hour at a release height 20 m above the ground. The release
was continuous from the previous three hours. The CF3Br release
rate is calculated to be 1.30 g/s.
Local terrain elevations near the release point are estimated to
be 17.4 m above the base elevation of the hill coordinate system, so
the net release height corresponds to the 37.4 m height level on the
ridge.
4
Meteorological Information
The plots of the 5-minute Hc and Fr values for this hour are
displayed in Figure 50. The 5-minute values of Hc are all greater
than the release height, except for the second 5-minute period when
Hc drops to 14 m. This sudden drop in Hc corresponds to a 1.0 m/s
133
-------
u>
I
I E
co
CO
:D
10 T
-. r- in
X
•
O uj M C3 CN
u o
-4-1 i-H
0) -
O -3-
C rH
o
o -P
ti a)
CD e
CO-H
ii U
u a)
CJ *~*
mm
oo -
cd
H
O
(U
C
O
t>0
O
C O
o <•
.rH O
CO I
en o
134
-------
135
-------
1
I o
I <-<
" <•> in c; f-
«) w « sr •
(VI (V (U — -•
a a. a. — o
o u o a. x
*» *» o i a
o o o
- D JC
7 oo <>• 19 ar
ULU N O t/1
CO
e -
(N
tn ^
c o
O iH
CO
c e
O -H
o u
(U
*o o*
3 J=
O *1
0) O O
C [j co
O "4-1 O
ao
bu
136
-------
o
Z
Release Height
Time (Hour)
2.8-
1.5-1
.0-
Time (Hour)
E'igure 50. Time series of 5-minute calculated dividing-streamline heights
(Hc) and bulk hill Froude numbers above Hc (Fr(Hc))
(Experiment 6, 10/31/82, 0700-0800 MDT).
137
-------
increase in Tower-A propeller wind speeds. This increase in speed is
not corroborated by the sonic anemometer wind data. The average of
the 5 minute Hc values over the hour is 35 m. Because Hc is less
than the release height for only one 5-minute period and this value is
suspect, this hour is considered representative of flow below Hc.
Fr remains steady during the hour varying from 1.0 to 1.3, with the
one hour value of Fr above Hc being 1.2.
The hourlr averaged wind speed and direction estimated at the
tracer release height are 1.0 m/s and 67.1°, respectively. Figure 51
shows the trend in wind speeds and directions between 5 m and 40 m
during the hour. The linearly interpolated values estimated at the
tracer release height are also shown. The 5-m wind directions vary
from 1 to 51° during the hour and the 40-m wind directions vary from
72 to 164°. The wind speed measured at 5 m oscillates from a high of
1.5 m/s during the beginning of the hour to a low of 0.6 m/s in the
middle of the hour, then increases to 1.2 m/s towards the end of the
hour. The wind speed measured at 40 m steadily decreases from 1.6 to
0.6 m/s during the hour, with the exception of a 1.7 m/s peak measured
during the seventh 5-minute period.
The trend in iz and aw values during the hour are shown in
Figure 51. The values of iz and aw estimated at the tracer
release height vary from 7 to 10% and 0.06 to 0.13 m/s, respectively.
The one hour value for
-------
(M
(s/iu) peeds
O
(9
(O
a>
U)
—r—
s>
m
CB
(M
IZI
o
o
H
(B
I
CD
W
3
r-l
CQ
o
H
e
o
u
y
•r-l
c
o
w
o
o
oo
o
o
o
CO
r-l
o
o
i m
co
At
u
(Baa)
111
(s/iu) M/>
139
-------
o
I
I
0)
o
H
O
8 5 S 8
co
•o
8 5 S
I S 8
3 3
a
O
CJ i
g o
,-t O
u
O CM
.c oo
in n
O r-<
-V
w o
P. C
-1
CO -H
o
.H
o
s-
CM
LO
ao
140
-------
141
-------
release height. This is the largest observed concentration found in
the Ct^Br modeling data set. Also, this sampler is collocated with
sampler 809 where the observed concentration is 425 ys/m-*. A
sharp decrease of observed tracer concentrations is found above the
mean Hc surface (1652 m).
Model Performance
The hourly averaged scaled concentrations estimated from the HER
(Flat) model are shown in Figure 54. The model grossly underestimates
the mean of the observed concentrations by more than a factor of 4
with~Co/Cp - 4.71 and r^ = .004. The maximum observed
concentration (445 ps/m^) is found towards the bottom-half of the
ridge at an elevation of 26 m; whereas, the maximum estimated
concentration (24 ys/m-*) is found towards the top of the ridge at
75 m.
Estimates from the HBR (Terrain) model, shown in Figure 55, are
considerably larger than the flat-terrain estimates with Co/cT =
1.02. However, the maximum estimated concentration (105 ys/m^)
found towards the center of the sampler array at 56 m is still much
less than the observed maximum concentration, and occurs at a location
too far up on the hill. The model performance apparently suffers
during this hour because the plume seems to have "dropped" in height
over a short distance from the release and HBR (Terrain) does not
account for this.
4.4.4 Experiment 8, Experiment-Hour 7 (0500 - 0600 MDT)
Release Description
The CF3Br tracer gas was released from position 215 at 25 m
above the ground for the entire hour. The CF3Br release rate is
computed to be 0.97 g/s. >
Local terrain elevations near the release point are estimated to
be 12.9 m above the base elevation of the hill coordinate system, so
the net release height corresponds to the 37.9 m height level on the
ridge.
Meteorological Information
Figure 56 contains plots of the 5-minute Hc and Fr values for
this hour. The tracer gas release height is less than Hc for the
first eight 5-minute periods of the hour. The average of the 5-minute
values over the hour is 35 m. Fr is steady during the first eight
5-minute periods of the hour ranging from 0.8 to 1.0. For the last
third of the hour, Fr rises to 1.5. The one hour average value of Fr
above Hc is 1.0.
Time series plots of the sonic anemometer data from Tower A for
this experiment-hour are presented in Figure 57. Large directional
wind shear is found between 5 m and 40 m until the tenth 5-minute
period when the 5-m wind direction shifts from 56 to 152°. There is
142
-------
en
^v
r/l
3.
CM
00
IJ
£ -u
O C
O
•O
4) ,-H
4-> 0)
O 73
.* O
*o c
0)
a -P
a)
00 b
0)
>
ca
u
o
£
o
01
Jj
3
00
Dd Q
m Z
o
-> o
I
e o
o o
143
-------
W CM
a oo
in r-t
C x
O O
•H rH
-P
n) -
c *>
a) c
o a)
c e
o -^
o u
01
*O &•
a) 5<
i-t W
B) ^
O
01 ^1
•a -a
<]) O
+> e
o
•a c
o 4-1 o
e o
o o
&o
•H
fc.
144
-------
1.5-
U
I
Release Height
Time (Hour)
~1
Time (Hour)
Figure 56. Time series of 5-minute calculated dividing—streamline heights
(Hc) and bulk hill Froude numbers above Hc (Fr(H))
(Experiment 8, 10/15/82, 0500-0600 MDT).
145
-------
(s/iu) peedg PUJM
3
O
(%) Zl
3
O
\
u
0)
o
H
o
Ll
14-1
o
cn
cu
-P
C
• H
e
H
Q
£
O
O
O
I
o
o
0)
-P •
CO 6
*—1
O /-N
& :
-------
little wind speed shear for the first eight 5- minute periods of the
hour. During the final third of the hour, the wind speeds between 5 m
and 40 m differ by as much as 2.0 m/s. The hourly averaged vector
wind direction and wind speed estimated at the tracer release height
are 83.7° and 1.2 m/s, respectively.
The iz and aw values are fairly steady and low for the
first half of the hour with iz values less than 10% and crw
values less than 0.1 m/s. During the second half of the hour, there
is a sharp increase in measured turbulence. The estimated value of
aw at the release height peaks at 0.4 m/s and iz peaks at 30%.
The one hour value of ow estimated at the tracer release height is
0.21 m/s (iz - 17.5%).
The vertical profiles of the hourly averaged wind direction, wind
speed, and temperature are shown in Figure 58. There is approximately
60° of directional shear in the vicinity (+10 m) of the tracer release
height. The wind speed decreases for the 10-m layer below the release
height. The wind data from the propeller anemometers at 5 m and 40 m
are in fair agreement with the sonic anemometer data at the same
levels. However, the vertical profiles of the propeller wind data
between 5 m and 40 m are not linear. The hourly average of the
5-minute propeller wind direction and wind speed values interpolated
to the tracer gas release height are 130° and 1.3 m/s, respectively.
The wind speed value is in close agreement with the hourly averaged
MDA value. However, the wind directions differ by 46° which indicates
that the linear interpolation is questionable. In this case, the 40 m
sonic wind direction (128°) is better than the interpolated value.
In summary, the plume was released on average below Hc into a
flow characterized by significant directional wind shear, 0.9 to
2.2 m/s wind speeds, and a hill Froude number above Hc equal to
unity.
»
Concentrations
The distribution of the observed hourly averaged CF3Br
concentrations over the surface of the ridge is shown in Figure 59.
The largest concentrations are found towards the middle of the sampler
array at an elevation that is near or less than the hourly averaged
value for Hc. The maximum observed concentration (72 ps/m^) is
found away from the estimated hourly-averaged plume centerline at
sampler 205. This large concentration may be associated with the last
third of the hour when Hc was near the tracer release height and the
wind direction was towards this sampler. The second highest observed
concentration (69 vis/m3) is located near the estimated hourly
averaged plume centerline at an elevation that is less than Hc. It
is most likely that the estimated hourly-averaged plume centerline is
incorrect. The concentration pattern is more representative of flow
from 130°, as indicated by the propeller anemometer data.
147
-------
IB
tn
CD
(M
S>
CT
(LU)
o
¥
8 <
CD
Q. P
I s
* £
e
o
IB
ID
(M
(B
a
IT
CO
•a
CO
O
f
TJ
Q.
(A
TJ
00 H
O Q
-< s
O
u o
o o
0) vO
-P O
Q) I
e§
X m
r-t O
O CSJ
^r oo
tn
(B
(9
c
O
O
£
a
•o
c
e-
w
00
in
!iO
148
-------
149
-------
Model Performance
Scaled concentrations from the HBR (Flat) model, shown in
Figure 60, are much less than the observed concentrations with
CQ/CP = 2.17 and r2 = 0.05. The largest estimated
concentrations (20 ys/m^) are found near the crest of the ridge at
a total of nine sampler locations.
Scaled concentrations estimated from the HBR (Terrain) model are
shown in Figure 61. On average, the model does well with Tfg/Cp =
0.88 and MCO/MCP = 1.29, although the correlation is low (r" =
.08). The observed concentration pattern exhibits a bi-modal tendency
which mirrors the wind direction pattern. This appears to contribute
to the low correlation with the concentration pattern estimated with a
Gaussian distribution about the mean wind direction.
4.4.5 Summary
The flat terrain model and the empirical HBR model have been
tested against CF3Br observations for four case study hours. These
case study hours represent three classes of meteorology (zr < Hc,
zr ~ Hc, zr > Hc) and include the highest concentrations
observed within the*three classes. Also, the MDA wind directions are,
in general, consistent with the observed concentration patterns during
these hours.
Comparisons were made between the sonic and propeller anemometer
data to evaluate the method used to construct the preliminary MDA for
the release height. In general, data from the sonic anemometers
compare favorably with data from the propeller anemometer for these
case study hours. However, there are many hours in the preliminary
CE^Br MDA that contain uncertain meteorological data, probably
because the linear interpolation between 5 m and 40 m is
inappropriate. The greater vertical resolution provided by the Tower
A propeller anemometer data will most likely provide more
representative meteorological data which should improve the modeling
results.
For each case study hour, the ratio of the maximum observed to
the maximum predicted concentration (unpaired in space) is closer to
unity for the empirical HBR model than the flat-terrain model.
Significant discrepancies between observed and modeled concentrations
(from the empirical HBR model) appear to result in part from
variations in the meteorology during the hour that are not adequately
represented in the "hourly" Gaussian plume formulation.
150
-------
O 0. 0. — I)
o u u a
~ ~ i u N
e b o i o
<-> a u
- o r
151
-------
-• OD r* ai
• • • "
U) (D O) ~'
If II M H It
o a. a. — u
u u o a. z
— ~ I u> x
» k o i o
u a u
» —
t
-« 0
• • • * OB CD r**
u> m 10 • «»
(U Al I -
n ii i i it H •
o a. a. ~ a. ai z
u u o a. u w
i u N w
o i o or
u a o
u
f)
)K
m
CO
CO
Of)
:co
x
a: •>:
CD oj rv
ro (s -. • o
U. (3 CD IT) (VI
o co —• ro •
•• •• u u ii
(K UJ
co :DCO
CJ LU M C5 r/i
'f, CM
3. 00
v-^ V.
LO
01 r-(
C ^
O O
•H H
-P
Cd -
U 03
C
0)
O
CD
O -H
O Ui
CD
•a a,
4) X
•—I [T1
cd ^
CJ
m ^^
0)
•o -a
6
o
•O "c
CD -H
'M Cd
O. !J
CD CD
OB 4-" •
'<-• H
CD « Q
> 33 2
Ct) X
3 .C
-------
SUCTION 5
THE PRELIMINARY TRACY EXPERIMENT
5.1 Geographic and Meteorological Setting
The Tracy Power Plant (TPP) near Reno, Nevada has been selected
as the site for the Full Scale Plume Study (FSPS), the third field
experiment in the CTMD project. The Tracy plant was tentatively
selected in May 1983 after (1) a review of operating power plants in
the western United States and of available topographic, meteorological
and air quality data, and (2) site visits to three power plants. The
final selection was made after a preliminary experiment that was
conducted at Tracy in November 1983. This section summarizes the
results from the preliminary experiment and discusses the plans for
the FSPS.
The Tracy station is operated by Sierra Pacific Power Company.
The Sierra Pacific personnel had agreed to participate in the program
and were very cooperative in the design of the FSPS. The power plant
is located about 27 km (17 mi) east of Reno, Nevada in the Truckee
River Valley. It has three units--53 MW, 80 MW and 120 MW, although
only the latter has been used recently. The 120-MW unit is serviced
by a 91.4-m (300-ft) stack. This stack was used to release the
oil-fog and SF$ during the preliminary experiment.
Figure 62 shows the location of the Tracy station on a 1:250,000
scale topographic map. The plant is located east, of the Reno-Sparks
metropolitan area and about 40 km (25 mi) east of* the Sierra Nevada
Mountains, the source of the Truckee River. The river runs eastward
near the plant and eventually drains into Pyramid Lake.
The 91.4-m stack is located near the southwest corner of the
plant (Figure 63) south of the Truckee River. The plant site is in a
relatively narrow valley with mountains surrounding the plant on all
sides. Figure 64 shows the location of the plant on a 1:62,500 scale
map. Mountain peaks rise to elevations of 900 m above the stack base
elevation within 6.5 km of the plant. The area is characterized by a
sparse vegetative cover of shrubs and grasses. Much of the
uncultivated area in the immediate plant environs is covered by small
boulders.
The Truckee River enters the valley through a narrow opening near
Patrick. It flows eastward just north of the plant and then takes an
abrupt turn to the north about 4 km east of the plant. The river
flows between two mountains at its northward bend. These two
mountains are the primary "target" areas for the dispersion
experiments.
153
-------
'
n.
Ui
o
&>
o
(0
0)
5
•o
I
tJ
CO
c
o
•H
00
0)
CM
\D
O
u
00
154
-------
155
-------
c
1-1
a.
o
Du
U
-------
In addition to the Tracy station, there is a diatomaceous earth
plant run by Eagle Picher Industries Inc. in the valley near Clark.
There are also a few ranches in the valley. Interstate 80 runs north
of the plant along the river.
No historical meteorological record is available from anywhere
in the valley. However, a previous field program (Kapsha et al.
1976), which included aircraft and mobile van measurements of SO2 as
well as pilot balloon data, suggested plume transport from the Tracy
stack to locations that would produce ground- level concentrations on
mountains 5229 (called beacon hill) and 5764 (called target mountain)
east of the plant and around the mountain 5610 complex northwest of
the plant. This field program was conducted in December 1975 and
experienced easterly wind flows* associated with transistory
anticyclones, as well as the more usual stable drainage winds. In
July 1983, as part of their initial feasibility testing, ARLF'RD
released oil-fog from the Tracy stack, and the smoke plume was
transported to and interacted with target mountain and beacon hill.
Although there is currently no corroborative data base, it is expected
that during the late summer, nighttime winds in the valley will be
dominated by terrain effects--probably producing a prevailing westerly
wind down the Truckee River with superimposed local katabatic effects.
5.2 Experimental Design
The preliminary flow visualization and tracer experiment that was
conducted during the period November 7-20, 1983 was co-sponsored by
EPA and the Electric Power Research Institute (EPRI). The EPR1
participation was in anticipation of a complex terrain field
experiment as part of their Plume Model Validation and Development
project. The results of the preliminary Tracy experiment will be used
to guide the design of the next PMV&D field experiment and to provide
the PMV&D modelers preliminary information on the relationship among
emissions, meteorological conditions, and observed concentrations in a
complex terrain setting.
The EPA objectives of the November experiment were (1) to assess
the feasibility of the Tracy site for the FSPS and (2) to obtain
sufficient information to design and plan the full scale experiment.
These objectives were satisfied and the site was selected for the
FSPS.**
The experimental methods were similar to those used and tested at
CCB and HBR and at the two previous EPRI field sites. The experiment
included:
• Releases of SV& and oil-fog from the 91.4-m stack;f
*Easterly winds were also experienced during the 1983 Tracy
experiment. See Section 5.5
**A detailed work plan for the FSPS was prepared. See also
Section 5.6.
t In the two small hill experiments, releases were made from mobile
cranes.
157
-------
• Ground- level SF^ concentration measurements at up to 53
sites;
• Fixed meteorological measurements:
a 150-m tower instrumented at four levels (5, 10, 100
and 150 m),
two 10-m towers instrumented at one level,
two monostatic acoustic sounders,
a doppler acoustic sounder, and
two optical crosswind anemometers;
• Two tethersondes;
• T-sonde releases with double theodolite tracking at two
locations;
• Two solar-powered electronic weather stations;
• Airborne lidar; and
• Photographs and videotapes.
The participants and their principal responsibilities were:
• ERT (EPA prime contractor)
(1) directed experiment operations in consultation with
other participants;
(2) operated the command post (Sierra Pacific Power
provided office space at the Tracy station);
(3) provided smoke candles, two carbon arc lamps, cameras
and personnel for flow visualization experiments and
scientific observations;
(4) provided, installed and operated two electronic weather
stations to measure winds and temperature at two
locations;
(5) provided, installed, and operated a 150-m tower
instrumented at four levels to measure winds,
turbulence, and temperature;
(6) provided for other site logistics as needed; and
(7) produced the master data archive ayd disseminated it to
participants.
• NOAA ARLFRD (via an interagency agreement with EPA)
(1) provided, installed, and operated two 10-m towers with
wind and temperature instruments telemetering to
command post for display and storage;
(2) provided data logging equipment; archived, reduced and
disseminated the meteorological data;
(3) adapted fogger to Tracy flue, provided oil and operated
fogger; provided SFg and injected into flue (with oil
fog);
(4) provided one tethersonde and operator;
(5) provided radios and repeater station;
(6) provided for photography contract (two photographers,
about four time-exposures/hr and one video during
daylight).
• NOAA WPL (via an interagency agreement with EPA)
(1) provided one tethersonde and operator,
(2) provided Doppler acoustic sounder and operator,
158
-------
(3) provided two monostatic acoustic sounders,
(4) provided two optical crosswind anemometers.
• SKI International (KVRl contractor)
(1) provided and operated airborne lidar (approximately two
3 hour missions per experiment).
• Rockwell. International (KPR1 contractor)
(1) provided 1-hour sequential (9 hrs) syringe samplers
(including spares) to operate on consecutive nights at
53 locations, and provided necessary deployment crews
and vehicles (including helicopter);
(2) determined SFg concentrations (by GO in all syringes
and bags with a turnaround time of about 24 hours;
(3) provided equipment and took T-sondes to 2.5-3 km above
ground at 1-hr intervals at two locations during tracer
releases; and
(4) provided survey/identification of all SFg samplers,
meteorological instruments, and fixed photography
locations.
• Research Triangle Institute (EPRI "External Audits"
contractor)
(1) provided an independent review of QC plans of other
EPRI contractors; conducted onsite systems audits of
the field measurements and the data handling activities;
(2) provided independent checks of the precision and
accuracy of the field measurements and data handling
results. Onsite performance audits were performed at
the Tracy Power Plant for the following measurements
systems: tracer, T-sonde, tethersonde, and two 10-tn
tower? with wind speed, direction, turbulence, and
temperature measurement systems; and
(3) provided reports to EPRI, through TRC, on the results
of systems and performance audit results.
• TRC Environmental Consultants (EPRI technical management
contractor)
(1) represented EPRI in the field and coordinated Kf'Rl
contractors, and
(2) undertook scientific observations and analyses as
appropriate.
5.2.1 Oil-fog and Tracer Gas Release System
ARLFRD provided an oil fog generator and a SFg release system
to inject oil-fog and SFg directly into the 91.4 m stack flue at the
Tracy station. The injections were made through a "door" into the
ducting leading to the 91.4-m stack (Figure 65). The SFg tracer gas
was stored in two compressed gas cylinders at ground level. Piping
carried the gas through a linear mass flow meter (LKM) system to the
point of discharge into the stack. The LFM measured and displayed the
rate of gaseous tracer discharge via real- time digital
159
-------
T60
-------
display, the total amount of gas discharged via a digital counter, and
the analog output voltage directly proportional to the flow rate. The
voltage was logged and monitored on a strip chart recorder. Pre and
post-test release weights of gas tracer cylinders were measured by
certified scales. Similarly, the oil consumption rate was logged to
document the quantity of oil-fog injected through the 91.4-m stack.
A nominal SF$ release rate of 1.26 g/s (10 Ib/hour) was used
for the Tracy experiment. The SFg and oil-fog releases commenced
approximately 30 minutes prior to the start of the sampling to ensure
that the tracer gas had actually reached the sampling grid when the
samplers were turned on. During the course of the experiment, the
91.4-m stack was vented through the use of a fan located at the bottom
of the stack. Only occasionally during the experiment was the unit
used to generate electricity. During these times the generation rate
was typically 20 MW.
5.2.2 Tracer Sampling and Analysis
Fifty-three syringe samplers were deployed by Rockwell
International (Cher 1984) to sample hourly concentrations. Sampling
sites were selected from 63 sampler locations (Figure 66).
Forty-three locations were specified to sample concentrations during
westerly winds and were used during the first few experiments. During
the course of the experiment, easterly winds were experienced
frequently so it was decided to select additional sites west of the
plant.
The samplers used in the program were sequential syringe samplers
manufactured by D&S Instruments. Nine 30 cm-^ syringes were housed
in each sampler. The syringe samplers functioned over a 9-hour period
with syringes sampling for consecutive 1-hour periods. Because of the
rough terrain, deployment of samplers required the use of either a
helicopter or a 4-wheel drive vehicle. t
In order to identify each sample by date, location, and time, a
numbering system was devised whereby each sampler, stake, and syringe
was given a unique number code. The numbers were printed on special,
double sided labels using a computer generated bar code system.
Labels were attached to the samplers, stakes, and syringes. Upon
deployment of each sampler, stake labels and sampler labels were
affixed to s. data sheet, which also contained information on the
syringe codes associated with the sampler. During analysis, a section
of the double-sided label from the stake, sampler, and each syringe
was transferred to the strip chart record containing the corresponding
trace. In this way, the strip chart record and the deployment data
sheet contained complete, redundant records of the deployment history
of each sample. When the strip chart records were read, all labels
were scanned with a bar code reader, and therefore all identifying
information was automatically transferred to the computer with
essentially no transcription errors.
A typical experiment ran from midnight to 9 A.M. Samplers were
loaded with syringes two days before the scheduled test. Loading
consisted of labelling the new plastic syringes, attaching the
161
-------
01
C
o
01
o
o
-------
needles, installing the units in the sampler, replacing any missing
septa, and setting the timers for the desired timing sequence and
starting time.
Deployment of the samplers occurred the day before the test and
exposed samplers were retrieved at the end of the test. For this
program a sufficient number of samplers were available so that
deployment and retrieval of samplers from the previous test could be
done simultaneously. Deployment and retrieval times were entered in
the data sheets for the two tests, and the stake identifying label was
also transferred to the deployment sheet at this time. To avoid
mixing samplers from different tests, the exposed samplers were marked
upon retrieval with an identifying tag. Analysis of the air samples
was usually performed on the same day as the test.
In addition to the syringe samplers, ARLFRD provided five
sequential bag samplers. These were collocated with five syringe
samplers for quality assurance purposes.
The analyses of the air samples were performed by gas
chromatography using an electron capture detector. Three Varian 3700
gas chromatographs were used, of which two were fitted with dual
columns to speed up the analyses. The columns were 3.175 mm x 1.829 m
stainless steel packed with molecular sieve 5A. The oven temperature
was 50°C, the detector temperature was 150°C, and the carrier gas was
nitrogen flowing at the rate of 50 cm3/min. Under these conditions,
the elution time of SFg and 02 were approximately 25 and 40
seconds, respectively. The entire analysis time (with backflushing
starting as soon as the oxygen peak began to elute) was approximately
2.5 min.
Each column-detector combination was calibrated using standards
obtained from compressed gas cylinders (Scott Specialty Gases)
containing manufacturer-certified SFg concentrations of 17, 83, 505,
and 1145 ppt SF6. Calibrations were repeated every two hours. The
response of the gas chromatographs remained constant within 10%.
Because the response of one of the gas chromatographs was slightly
non-linear, the response curves for all calibrations were determined
by fitting calibration data to a quadratic equation of the form:
concentration = ah(l+bh)
where a and b are calibration constants, and h is the peak height.
For the highest concentration used (1145 ppt) and the most non- linear
response gas chromatograph, the ratio of the quadratic term to the
linear term bh was approximately 0.15-0.20 or 15-20%. All calibration
constants were stored on disk using a DEC Professional 350 Computer.
The SFg peak heights were digitized electronically using a Science
Accessories Corporation sonic digitizer Model GP6-40 connected
directly to the computer. Peak heights were converted to
concentrations by applying the appropriate calibration constants.
163
-------
5.2.3 Plume Photographs
Two dedicated photographers took 5-minute exposures of the
oil-fog plume approximately every 15-minutes during the course of the
experiment. Their locations are depicted in Figure 67. The plume was
illuminated by two carbon arc lamps- one located south and one north
of the river. In addition, the ERT scientific observers took plume
photographs of interest from a number of locations.
5.2.4 Airborne Lidar Sampling
SRI International provided the ALPHA-1 airborne lidar to document
the three dimensional interaction of the oil-fog plume with the
terrain. The aircraft flew a "creeping-ladder" pattern from the stack
to approximately five to ten kilometers downwind. The flight pattern
legs were approximately at right angles to the mean plume direction.
The aircraft typically flew two three-hour missions during each
experiment. Facsimile lidar cross sections of the plume were
available at the end of each experiment for analysis by the project
scientists.
A separate report on the Alpha 1 observations was prepared by
Uthe and Morley (1984). Some of their material is presented in
Appendix C.
5.2.5 Meteorological Measurements
Table 15 lists the meteorological instruments used during the
November experiment. Their approximate locations are depicted in
Figure 67. The minisondes were released and tracked hourly on the
hour. The tethersondes were used to obtain vertical profiles of winds
and temperature up to 600 m. The tethersonde ascents occurred during
the first half-hour of each hour and the descents during the second
half. The tethersonde data were processed continuously so that the
data were available for near real-time analysis and operational
planning.
All data (including SFg) were achived and identified by the end
time (PST) of the averaging period, e.g., an hourly average wind
direction for the period 0200-0300" is identified as a 1-hour average
ending at 0300. Shorter-term measurements were also identified by the
end time of the measurement, e.g., Doppler wind data at 0200, 0220,
0240.
5.3 Preliminary Field Study Results
5.3.1 Summary of Data Base
Ten experiments were conducted for 73 hours during the period
7-20 November 1983. The use of 43 sampling sites was planned for the
first four days and 53 for the last six. Because of high winds and
snow, only 37 and 30 samplers were deployed for two experiments.
During the first four days, unfavorable weather conditions resulted in
164
-------
.
! :vg> \ W 7p--W: m^r
'!
h-7
'WV/'*f}\/SS&jg$^ \
\J'.v^oH;;jo:\ \f li ] .,
3
O
C
0)
a,
x
la
o
(-.
a,
•£>
0)
•x)
165
-------
TABLE 15. METEOROLOGICAL INSTRUMENTS AT THE
PRELIMINARY TRACY EXPERIMENT
Instrument
T-sonde (two)
Tethersonde (two)
Doppler acoustic sounder
To provide wind and temperature data
throughout the boundary layer (up to
3000 m) near the west and east
boundaries of the experiment-region.
To provide wind and temperature data up
to 600 m in (1) the flow upwind of the
stack and at (2) the flow upwind of the
bend to the north of the Truckee River.
To provide real-time information on the
vertical profile of winds near the
source location; data are archived
every 10 minutes.
Fixed meteorological sensors To provide
two 10-m towers
(winds and temperature)
two electronic weather
stations (winds and
temperature)
one 150-m tower east of
the stack (winds, tur-
bulence and temperature
at four levels -5m,
10 m, 100 m, 150 m)
Monostatic acoustic
sounders (two)
Optical anemometers (two)
(1) real-time information on
meteorological conditions in the
high terrain,
(2) data on the drainage flow, and
(3) data (historical only) rep-
resentative of plume conditions.
To provide information on the structure
of the boundary layer (1) upwind of the
source, and (2) upwind of the bend in
the river.
To provide information on the drainage
flow on the south side of target
mountain.
166
-------
few ground- level SFg concentrations. Table 16 Summarizes the
available concentration data from the entire Tracy experiment. A
total of 3167 SF$ concentrations are in the data base. All the
remaining measurements have also been delivered to ERT and are now
part of the Tracy Preliminary data base, which is summarized in
Table 17.
The complete data base was delivered to EPA accompanied by a
descriptive data report. The data base will be available from the EPA
Project Officer.
Maps of the SFg concentrations and the various meteorological
data were examined to assess which hours are sufficient for modeling.
From the 68-hour tracer data base, it was judged that approximately 34
"good" hours and 13 "marginal" hours* are available for modeling
purposes. Table 18 summarizes the modeling data base.
5.3.2 Overview of Results
The November experiment did achieve the program objectives:
(1) the Tracy site is feasible for the FSPS; and (2) the tracer gas,
meteorological and photographic data base are sufficient to design the
FSPS (see Section 5.6). Although the data base from the November
experiment is small compared to the CCB and HBR data bases, the
experiment did capture a wide variety of dispersion conditions.
Tracer gas concentrations were observed in the primary target areas
during stable conditions and also during windy, neutral conditions.
The data can be used to evaluate the CTMD modeling approaches and to
extend the hill and ridge data bases.
5.4 Example Results from Specific Experiments
To illustrate some of the experimental results, a few hours from
three experiment days are described qualitatively in terms of the
relationship among emissions, observed meteorological conditions and
subsequent ground- level tracer gas concentrations. Predominantly
stable atmospheric conditions occurred during the first two days and
windy, neutral conditions on the third day.
Experiment 5 (November 12, 1983) 0000-0500 PST
concentrations were observed in the high terrain of the
primary target areas during the course of Experiment 5. Figure 68
shows the geographical distribution of hourly concentrations observed
during the first five hours. In this figure, the dotted lines
indicate the position of Interstate- 80 and the Eagle Picher haul
road. The river flows from the lower left to the upper right and it
parallels Interstate-80. The X in the middle of the diagram
represents the location of the stack. The concentration range of
is indicated by the numerical symbols, defined as follows:
*The criteria for good hours are qualitative and subjective.
167
-------
TABLE 16.
SUMMARY OF ACQUISITION OF SF, CONCENTRATION DATA
o
Experiment
1 (A)
2 (B)
3 (C)
4 (D)
5 (E)
6 (G)
7 (H)
8 (1)
9 (J)
10 (K)
Date
Nov 7
Nov 8
Nov 9
Nov 10
Nov 12
Nov 14
Nov 15
Nov 16
Nov 18
Nov 19
Number of
Sampling
Sites
43
37
43
43
53
30
53
53
53
53
Sampling
Hours
9
9
9
9
9
9
9
9
9
9
Samplers
Analyzed
43
37
43
10
53
30
53
53
53
53
Hours
Analyzed
6
4
2
2
9
9
9
9
9
9
Total Number
of Samples
Analyzed
258
148
86
20
477
270
477
477
477
477
Total
68
3167
168
-------
TABLE 17. DATA BASE, PRELIMINARY TRACY EXPERIMENT
Information
Concentration
emission rate
150-m tower met data:
WS, WD, o0, w, ow,
T, AT from 5, 10, 100,
150-m, levels
10-m tower met data:
WS, WD, T
Doppler sounder data:
WS, WD
electronic weather station
data:
WS, WD
monostatic sounder
optical anemometers
tethersonde data:
WS, WD, T
T-sonde data:
WS, WD, T
Lidar data
Photographs
Base Log
Observer notes
Averaging Time
1-hour
15-min
5-min & 1-hour
5-min
10-min
1-hour
continuous
5-min
instantaneous
Comments
Available for 68 hours
instantaneous
continuous during each
mission
5-min
Two towers
T available from one
of the two stations
Facsimile record from
two locations
Path average wind
speed for two paths
Two profiles per hour
to 600 m at two
sites. RH
available from
Clark site,
pressure from site
west of plant.
One profile per hour
to 3,000 m at two
sites
Facsimile records and
aircraft position
(See Uthe & Morley
1984)
Five minute exposures
every 15 minutes
Available from
Project Officer
Definitions
WS wind speed
WD wind direction
-------
CO
H
E
t—'
K
a,
X
o
X
OS
a*
s
ly winds, neutral conditions
^J
c
CD
4J
CO
o
60
c
0
Li
j->
M
d variable winds; considerable
c
cd
jj
f*
bO
•H
•— i
^
Li
O
0)
cn
•rl
Li
CD
£
(X
sterly winds; no impact on samplers
ca
CD
^J
.c
60
.i-l
i— I
^
Li
a)
>
y and southerly winds throughout
!—(
Li
0)
^j
cn
cd
o>
4J
.c
60
uj
4J
C
0)
s
.1-i
Li
0)
&
CD
on primary target area as well as new
-f-
u
fcj
o
*O
CD
Li
3
Cfl
cd
CD
£
-U
cn
i
0)
-tj
0
CQ
0)
P,
I
CQ
60
.r<
C
r— t
t— 4
cd
vD
w
•d
o
o
erywhere after 0200
^
CD
j_)
cn
o
g
cd
fcf
w
•o
.1-1
S
•o
c
cd
^
•o
3
O
O
o
0
CM
O
.1-1
jJ
C
cn
•o
.r-l
3
gj
3
^J
C
jj
CO
.H
cn
u
0)
ex.
\o
GO
3
00
w
CD
4->
CD
60
C
•H
X
Li
cd
>
w
•o
c
.1-1
3
o
o
0
1
o
o
CO
o
0)
Li
CD
5j
•>^
Li
CD
CD
O
—4
H
c-1
O
cn 60
L, c
3 -H
o -*
X ffl
•o
•c o
cd u
> o
£ £
co o
O O
O O
£ £
O O
o o
£
CM
O
£ £ £ £ £
H r-l CO CM H
O O O O O
00 vO *«T CO vO
00
T-J
5
H
cn
Li
I
\O CO
o o
o o
o o
CO CO
oo oo
oo
o o
i i
CM O
O O
CO
CO 00
oo -v
-v O
ON r-l
o
I
o
O
CO
00
CM
o o o o o
i I I I I
O O r-l CM O
O O CM CM O
CO CO
CO 00
CO CO \D 00 CO
00 00 r-l H 00
>» I
ID IT)
I
i-l r-l H r-t
o
..-4
4J
CO
Li
-J
CO G)
e -i
PI -a
a)
ri
Ul
9?!
CQ OQ H OXr-ir>Sx{
V—' *—' N_^ V_^ \^*• \^' V^*- S^/ x^x
0^ CO *^ ^^ ^D r*"» 00 ^1^ ^?
o
•- u
•O 60
o
O CO
60 -i-l
III 0
CO J
04 O O
* * +
170
-------
Data for Hour Ending 0100 on 12-November-1983
2 1 *
J .3--"-^4
/' 9
.-&*'
xVQ.,---«'
^_/«,_
o o
Data for Hour Ending 0200 on 12-November-1983
0
2
30
,' 0
o o
Figure 68. Hourly SF6 concentrations November 12, 1983 0000-0500.
171
-------
Data for Hour Ending 0300 on 12-November-1983
1 1
...-5V p 0 o
~
/ °
o o
o S o
y
f
y
o
0
Data for Hour Ending 0400 on 12-November-1983
0 0 5 2 -
2 V'V
3
4 i
" ~: ~ 5 4 0 0
f"1 3 . <• *•
0 1
1 0
0
0
Figure 68 (Continued). Hourly SF6 concentrations November 12, 1983
0000-0500.
172
-------
Data for Hour Ending 0500 on 12-November-1983
z
>""-s ..............
5
00
.'*-.
— ^ ™
5 0
0
0
0
Figure 68 (Continued). Hourly SFg concentrations November 12, 1983
0000-0500.
173
-------
Concentration Range (ppt)
Symbol
0
1
2
3
4
5
6
7
8
9
During the first three hours, the highest concentrations were measured
at elevations typically 60-80 m above the top of the Tracy stack.
During the period 0300-0500, the highest concentrations were observed
along the valley floor in the gorge where the Truckee River bends to
the north. The highest concentrations, their location and elevation
above the stack base are (Cher 1984):
Time
(Ending Hour)
0100
Min
0
20
40
60
80
100
125
150
175
>200
Max
19
39
59
79
99
124
149
174
199
Sampler
Site
8
34
31
Cone.
ppt
211
178
159
R
km
6.2
3.6
3.2
e
deg
68
82
62
AZ
m*
169
108
170
0200
31
28
34
6
176
154
149
146
3.2
4.5
3.6
5.3
62
71
82
77
170
113
108
159
0300
1
2
6
310
169
119
4.6
5.0
5.3
95
90
77
158
170
159
0400
25
26
37
27
241
216
207
206
5.1
5.0
6.2
4.8
73
73
63
71
20
-15
-13
-8
174
-------
0500 26 190 5.0 73 -15
37 188 6.2 63 -13
25 171 5.1 73 20
*AZ measured from the base of stack. Stack height is about 91 m.
Figure 69 gives a 5-minute exposure of the oil-fog plume as taken
from Prospect hill (photography position #1, see Figure 67) at 0000
PST. Notice the dispersion of a stable plume from the Tracy stack
toward the east at an elevation near the top of the 150-m tower.
Figures 70 and 71 show 5-minute exposures, again taken front prospect
hill, at 0015 and 0030, respectively. The plume is apparently
producing ground-level concentrations on beacon hill. A 5-minute
exposure taken from Old Lonesome at the west end of the valley
(photography position #3) is shown in Figure 72. It illustrates the
extensive crosswind growth of the plume and shows some plume material
reaching the northwest sections of target mountain.
Hourly average wind directions and speeds taken at the 150-m
tower during the five hours are given in Figures 73 and 74. Winds
near plume height during the first three hours have a westerly
component, consistent with the observed concentrations, photographs
and observer comments. During the last two hours the winds have an
easterly component. Figure 75 gives 5-minute values of wind direction
and speed and ov from the 150-m le-'el of the tower. Again, the
winds are primarily westerly until about 0300, and then primarily
easterly. Notice the increase in turbulence after 0300. Figure 76
gives doppler profiles of wind direction and speed taken at 0020, 0130
and 0400.
A five-minute exposure of the oil-fog plume taken at 0130 from
Prospect Hill is shown in Figure 77. The plume is still being
transported toward the east, but notice that the top of the smoke
plume is below the top of the 150-m tower. Figure 78 shows the
geographic distribution of winds at the approximate elevation of the
plume above the base elevation of the 150-m tower for each of the five
hours. This figure and the doppler wind data illustrate the
complexity of the wind fields in the valley.
In any event, it is clear that the SFg plume was transported
directly to the samplers during the first three hours. Values of Hc
were calculated based on the elevation of beacon hill and using the
150-m tower data. A time series of 5-minute values of Hc are listed
along with other meteorological data in Table 19. The three hourly
Hc averages are 218,222 and 219 m. The elevations of the higher
concentrations are always less than Hc during these three hours.
Plume material was not transported to the samplers at the highest
elevations.
175
-------
§
o
o
00
CM
lib
176
-------
o
o
+J
fl
CO
00
u
o
I
o
00
•H
fc,
177
-------
o
o
+3
ffl
co
00
CM
0)
.Q
O
SB
cd
u
a)
2
CO
o
0>
4J
•r-l
•T-l
fc
0)
u
00
•H
Et.
1/8
-------
o
o
CO
oo
CM
a)
.o
o
s
I
O
u
u
o
d)
(D
tM
O
Li
Ed
179
-------
CM
II
II
K
.•....
'
I !
• 2
i !
a
1
I J
i '
i
8
R
I
It
u
o
f
O
2
0)
a)
ca
C
o
•o
co
0)
u
J!
> u>
180
-------
OL
I
£C
X
7 C
«- CM co 3 in
O O ^5 ^5 ^?
D
o
« o
" I
•
i
10
rg o
u>
s a IB
ea IA
3)
O
O
in
o
-o
-------
«.* i
LO
O
CO
d
q
LO
Q "ro I
> ^ i
> -o •
L_
O
(O
CO
O
CO
CO
8
CO
g
Is:
o
r-
CM
q
cvi
o
in
O
cn
o
to
o
CO
o
o
o
lJ
o
O
I
o
is
o
in
en
•o
Q
SJ
aO
•H
182
-------
400
375
350
325
300
275
250
225
200
175
150
125
100
75
50
-9.9
l.B
1. I
I ..3
-9.9
1.8
1.4
1.7
1.9
2.1
2.0
2.1
1.4
0.6
1.0
j. a. / a. a
999
265
247
249
999
203
222
208
213
229
235
236
214
157
82
O t S>0
400
180
270
360
0.0
2.0
HEIGHT
CM)
400
375
350
325
300
275
250
225
200
175
150
125
100
75
50
SPEED
(M/S)
-9.9
-9.9
-9.9
1.6
0.0
1.6
1.4
1.4
1.7
1.3
0.9
0.6
0.5
1.0
1.2
DIRECTION
(DEC FM N)
999
999
999
169
189
168
163
150
147
153
171
178
220
240
221
400 -
400
375
350
325
300
275
250
225
200
175
150
125
100
75
50
. x :i. :;>
-9.9
-9.9
4.3
4.2
4.5
4.0
2.9
3.0
2.5
2.2
2.1
2.4
2.4
1.6
1.5
t ISO
4.0 6.0 8.0 10.0 12.0
MIND SPEED (M/S)
90 180 270 360
0.0 2.0
999
999
183
179
182
172
151
126
109
82
88
89
88
72
67
AS O
400 _"
~" iTo' 6To ' 5to
WIND SPEED (M/S)
90 180 270
10.0 ' !?.(
360
350 ••
300 ..
250 ..
200 . .
150 ..
100 ..
•4-
t
50 f
?To 4.0 ?To 8.0 10.(
WIND SPEED (M/S)
17 o
Figure 76. Doppler sounder wind profiles.
183
-------
-u
cd
oo
er>
o
a)
U
0)
U
U
W
O
Q)
4)
I
•>
•rl
h.
2
2,
184
-------
, ;: • . • / y,7^>, • " \ „ • ~ '( Mf&. 'p V ^>J S i") Vv^'/VV-''! f' I *r "77^' v~f "~ ' ; ^ 'J\
•/] ^sn-> ^A4i'>w%^ ^F' ',"T/
.' -< ' vi V" ,.VsJ- V* i '(^'-^vV--^K-Mr r^" " • ^'>H^ /S2 - i; )s
1.' ' •: s' v; 1 A -1 '^ '' ^ - ^":7J//ii"^"^Tf V'7y "-vv - i RI;$V'
^^v^>l^>:;^ft|K^5W^
^:.. ,v. v. fit ^.^!/»rid J
•O^.-'j'i C'M /i:(((0 <^M:f^f U^sll^
w«y»^\f ^' - Ki" ^M^^'
^M^;Jiwm&^{ i' M IA1 ;< A ^ ^l% • r A'
W^wte^H^TT^rrrr^
J>A 0! \.'JZ\ Vil\^^v$v\ ok/vU^^HX, (^- / •' ( ] \
l^.p^^i^^
li|l«S!rIv* -i-VMi
-)? JJ-///('u^j ^(/^^S^^cv^ v ( p/^"^"""' _' l"J ! - ^"^-,< <
^>^:^4^41^|A:--f^^^^~ "^v" ^'' "~~ 1 * ^ ''^'" "^
[?^>i^4>^?»^
fctoibmw;
V'^'fi^'r-^rx''),)}^.;!'/ - , r"~i/^7 F '
^^T^ffi-f/,1'--;"
wov'^v/ (iv •••:
•Sr-^Si '-^^^ '"; '
{(l^y^^i^.S .
J'ii V/^flV'-V^" ^
^00
hi
,-vc
'{ '''
=^
Vf-^Lr. v-' ^-:- / ^'!' "
-^' ] --.j,^' •< f/c;:^^^1, «
)^.'K.- ;/'V^\ -' -•
J /V •— :/, i ^ /> 1 .,.
.-- .. «,;-.-. .— =---i , -,- -. r , -„,.
iSS'
^ «*
-'**\
".•''k-
.'••&'<
^ JJ ^ ^/ i
LT
0)
4J
cd
g
• iH
X
o
L<
o.
a
U)
• H
T3
y
•*-<
x;
a •
cd i -*
tj ,^
60 60
o -^
0) 0)
o jc
-------
T-=.--vfm-M--- "i*'-m,-\-i 4^A"iY^--?--"'-v- - -
\:O3tf'/-'-»;, ?'j: fU^MK ^i !•
- /r^-^/-/ ! , \U (|VA>jk-V(r- v' J V,
*-^T- •'•' ' m^^-i ""i r^
f|i;W-|telft^M^>
n.^y^v./g Lfefen\sr) T- r-^ H
;^: Jf^'DRli)PK^ ?' "= •-h,L': ;j
^»»fe>K^^^ai
' /,
X3i
i ,J*f - -
\' i
*y\
••\r
•a
c
U-l
o
•r* D.
TJ
(1)
O -1-1
•f-4 Cd
•s e
Q> •<->
a) X
U O
60 t-l
O O,
0) (X
O f?
186
-------
tn
•O
c
c x:
o oo
in F-I
.- O.
•O
0)
O -t->
.-<
-------
4J
of
en
•o
c
C A
O 60
•H *tH
4J 0)
a .c
.o
OS i-l
-rH Q.
•o
<]>
y *>
•rt a)
O. -r<
O) X
u o
60 U
O Q.
0) Oi
O 0)
•O
0)
-u
c
o
00
u
60
188
-------
4J
in
•O
c
o •
•4J
C J=
O 60
•*-4 -f-l
•U 4>
.»< 0)
^ E
•»J 3
01 r-l
•1 O.
•o
0)
O -P
•^ «
je 6
Q, .^i
-------
Hr
Sec
TABLK 19. FIVE-MINUTE METEOROLOGICAL DATA*
NOVKMBER 12, 1983, 0000 0300 PST
ws
WD
Ft-
Bulk H
0
0
0
0
0
0
0
0
0
0
0
1
1
1
1
1
1
1
1
1
1
1
1
2
2
2
2
2
2
2
2
2
2
2
2
3
300
600
900
1200
1500
1800
2100
2400
2700
3000
3300
0
300
600
900
1200
1500
1800
2100
2400
2700
3000
3300
0
300
600
900
1200
1500
1800
2100
2400
2700
3000
3300
0
1.1
.9
.5
.3
.3
.9
1.1
1.3
1.1
1.1
.5
.7
1.6
1.7
1.9
.7
.7
.1
.9
1.1
.7
.7
2.1
1.5
1.3
1.5
3.2
2.5
1.7
2.3
1.7
1.6
1.1
.7
1.7
1. 7
237.0
249.0
253.0
253.0
10.0
70.0
112.0
70.0
87.0
54.0
7.0
310.0
291.0
296.0
288.0
347.0
237.0
305.0
270.0
249.0
277.0
281.0
319.0
317.0
314.0
282.0
249.0
253.0
218.0
256.0
323.0
298.0
29.0
128.0
77.00
38.0
-.1
.1
.1
.1
.1
.1
-.1
-.1
.1
— 1
-.1
-.1
-.1
.1
.1
.1
.1
„ f 2
.1
.1
-.1
.1
-.1
.1
.1
.1
.1
-.1
.1
-.1
.1
— 1
.1
-.1
.3
- .1
.253
.291
.093
.046
.202
.330
.215
1.021
.167
.129
.286
.247
.360
.258
.345
.445
.097
.081
.245
.394
.349
.127
.451
.217
.145
.749
.743
.489
.383
.711
1.397
.360
1.092
.651
.730
1.311
.200
.120
.130
.140
.130
.130
.200
.200
.130
.200
.130
.200
.130
.130
.120
.120
.120
.130
.120
.120
.130
.280
.200
.130
.200
.200
.200
.280
.120
.280
.280
.360
.510
.200
.200
.200
178.0
198.0
143.0
166.0
187.0
246.0
254.0
240.0
242.0
259.0
253.0
267.0
247.0
238.0
229.0
211.0
230.0
253.0
220.0
188.0
230.0
159.0
226.0
233.0
227.0
238.0
206.0
227.0
216.0
219.0
224.0
182.0
222.0
179.0
260.0
213.0
1.1
1.2
.5
.6
.8
1.2
1.3
1.3
1.4
1.8
1.4
1.2
.9
1.0
1.3
.3
.9
1.0
.4
.1
. 7
.1
.8
1.2
1.1
1.2
1.7
1.5
1.1
1.3
1.0
.5
1.0
.3
1.2
.6
.2
.2
.1
.1
.1
.1
.1
.2
.2
.2
.1
.1
.1
.1
.1
. 1
. 1
.0
. I
. L
. L
.0
.1
.1
.]
.]
. 2'
.1
.1
.1
.2
.2
.2
.1
.1
.2
.0328
.0297
.0318
.0363
.0260
.0256
.0253
.0345
.0333
.0325
.0303
.0317
.0325
.0327
.0288
.0320
.0339
.0315
.0350
.0349
.0358
.0402
.0311
.0312
.0308
.0297
.0220
.0250
.0319
.0303
.0353
.0333
.0297
.0374
.0323
.0339
.0376
.0381
.0381
.0333
.0360
.0371
,0349
.0381
.0386
.0409
.0423
.0385
.0365
.0370
.0360
.0376
.0360
.0400
,0365
.0309
.0355
.0309
.0349
.0355
.0349
.0338
0327
.0390
.0395
.0400
.0360
.0321
.0338
.0326
.0284
.0327
*Collected at the upper level of the 150-m tower.
190
-------
During the last two hours, 0300-0500, the highest SFg
concentrations were measured in the gorge at the valley floor level,
yet the meteorological data suggest the plume was transported toward
the west. Figure 79, a 5-minute exposure taken from the west end of
the valley at 0315, shows the plume near the stack being transported
toward the west with considerable plume material remaining in the
target area. The Alpha-1 lidar observations* (Figure 80) also show
plume material at valley level in the gorge.
How did elevated plume material enter the gorge and produce high
concentrations on the valley floor? ERT observer comments suggest the
turbulent transport by drainage winds off beacon hill and target
mountain. The monostatic sounder at Clark also showed an increase in
turbulence after 0300. See also the 150-m tower ov data in
Figure 75. Figure 81 depicts a photograph of the monostatic sounder
facsimile record for 12 November. During 0000-0300 the sounder data
show decoupled layers with waves in the lower atmospheric boundary
layer. At 0300 the record shows the interaction of layers—implying
transport of elevated plume material from aloft—and complete mixing
by about 0440.
In summary, the first three hours of Experiment 5 illustrate
stable plume impingement conditions. The SF$ plume evidently
produced ground-level concentrations at samplers whose elevations were
below the calculated Hc. The hours 0300-0500 illustrate the
occurrence of high concentrations on the valley floor. The
concentrations evidently resulted from the turbulent transport by
drainage winds of elevated plume material.
Experiment 7 (November 15. 1983) 0300-0800
SF$ was observed by samplers located east and west of the TPP
stack during the morning of the 15th. Figure 82 shows maps of hourly
concentrations over the five hours from 0300-0800. Plan view maps of
the winds representative of 100 m and 150 m are given in Figures 83
and 84, Evidently, the winds in the valley switched back and forth
from westerly to easterly and then back to westerly during the five
hours. Figure 85 gives time series of wind direction and speed and
ow measured at the 150-m level of the tower. Notice the very low
values of ov. The highest concentrations, their locations and
elevations are (Cher 1984):
*See Appendix C for an explanation of the Alpha-1 observations.
191
-------
<*>
00
CM
H
U
I
co
CD
u
I
(ft
O
g
2
-------
193
-------
November 12
Time(PST)
600
• 400
200-
0200
0100
0000
•Time(PST)
600
0600
0500
0400
0300
600
01
0900
- Time (PST)
0800
0700
0600
Figure 81. Eagle Picher acoustic sounder records.
194
-------
Data for Hour Ending 0400 on 15-November-1983
0
0
1 ,' 0
o o / r, o
2 o cAi ° o
; 6
-..-'V 8| 5
0''
0 0
Data for Hour Ending 0500 on 15-November-1983
o
o o . . ' o
°
o
o o
. .....
Figure 82. Hourly SF6 concentrations November 15, 1983 0300 0800.
195
-------
Data for Hour Ending 0600 on 15-November-1983
° ° ° 0 -
0 o V o"
0 Ov' < 0
..
0 --" °
X
p
/
/ 1
/
o o
o
t
0
0
Data for Hour Ending 0700 on 15-November-1983
o o >--*
0 0 / „ 0
2 14£ °
3
1 6
0
0
Figure 82 (Continued). Hourly SF6 concentrations November 15, 1983
0300-0800.
196
-------
Data for Hour Ending 0800 on 15-November-1983
o o o ° '
0 ft 1 ,-1'" """*
0 o: i 0
° 2 I2?£l !
-'-- .............
2 3
0
Figure 82 (Continued). Hourly SFg concentrations November 15, 1983
0300-0800.
197
-------
M
•o
c
e
o
o •
r-t O
O
0) O
U 8
oo
-------
0)
T>
6
o
o
c
o
XI •
.1-1 <*)
U 00
*> 0>
01 r-l
•*4
•o -
in
O rH
a. a>
a) u
u e
tO 4)
o >
0) O
u z
•o
01
.1-1
,J
c
o
o
00
Ll
bO
199
-------
m
•G
C
o
o
.o •
•r< CO
u oo
at r-t
•«-l
T> -
J= U
O. (U
o >
-------
u
•O
• «-4
*
e
o
o
c
o
u oo
4J O>
01 r-1
•i-(
T) -
in
4) O
0 2
•a
0)
g
••-I
4->
c
o
ro
CO
a)
»-l
3
201
-------
tn
C
-»-*
»
e
o
o
I-l
c
o
••-I
4J
3
J3 •
•rH CO
[J 00
jj O>
Ul r-l
.«-*
•O -
in
O r-l
.1-*
J= fcl
a
-------
n
•O
E
i
o
u->
rH O
O
-
in
Q. 0>
td JD
t. e
60 0)
o >
0) O
O Z
00
u
10
.1-1
203
-------
: t ,%w >^f«Wf^ W-^
;-^ :-;',- ^>i/t 4:v^4^^^Kl/a^ ^ft
'-.'•^ '- ?"!«/ ^>%^»!^^^
f/a f ^^H-'^Wjnii/^
& ^^^Sft^
SIS4/-- --l-W^yT*-^-- -K+ • : / ,f—iv\-J \\,"_~, I
J)
W
TJ
g
o
m
r-i O
O
u co
JJ 0»
o> H
•*-i
•O -
m
O r-l
•VH
JS it
o. a)
u g
40
0) O
•O
0)
3
O
o
co
cu
60
204
-------
r—*
10
*r4
•o
0) O
O 2
•o
a
3
O
O
00
s
oo
205
-------
01
•o
e
i
o
in
r-l O
O
(w 00
o o
I
c o
o o
«* n
*> o
OO
•o -
in
a) A
U E
60 0)
o >
0) O
O 2
•o
a)
3
O
o
03
-------
in
TJ
C
f
o
1/1 •
r-l O
O
<4-i CO
O O
I
G O
O O
.1-1 CO
•^ o
.0 -
••-< in
(j 00
•o -
in
o tH
»iH
A U
CU V
0) A
t4 6
60 0)
o >
0) O
O 2
•o
a)
3
.•-I
*)
c
o
o
00
il
I,
-------
CN
d
o
c>
00
q
ci
q
d
SCM
O
O ci
I I
J I
s
o
o
e
I
O
ffl
-P
(0
•a
a)
o
oo
o
o
'*
o
.(J
o
oo
0)
u
•H
fa
Q "3 I O
> ® i n
5 T3 ! CO
O
o
CO
-------
Time
(Ending Hour)
Sampler
Site
Cone.
ppt
R
km
e
deg
AZ
ni*
0400
10
3
2
648
216
180
5.3
5.4
5.0
89
84
90
241
154
170
0500
2
3
10
417
313
269
5.0
5.4
5.3
90
84
89
170
165
241
0600
28
1
31
25
92
83
67
65
4.5
4.6
3.2
5.1
71
95
62
73
113
158
170
20
0700
56
28
1
22
53
126
96
77
76
76
5.3
4.5
4.6
4.1
5.8
231
71
95
103
246
254
113
158
56
224
0800
58
56
60
462
211
139
4.8
5.3
8.7
265
231
225
208
254
370
*AZ naasured from the base of stack. Stack height is about 91 m.
During the hour 0300-0400 the plume was transported directly (but
slowly) from the stack to the southwest corner of target mountain.
Figures 86 and 87 show photographs of the oil-fog plume taken from
Prospect hill at 0315 and 0330. The first five-minute exposure shows
a very stable plume being transported from the stack toward the east
at an altitude near the top of the 150-m tower. The second exposure
shows plume material in the valley south of beacon hill being advected
into target mountain. Figure 88 gives a photograph of the plume
209
-------
PI
O
CO
00
o
55
I
o
Q)
U
Dl
O
•H
0)
•H
fe
oo
u
ao
210
-------
o
PI
f)
O
-p
R)
oo
en
u
0)
1
00
o
u
•H
fa
211
-------
o
CO
CO
o
CO
oo
I
CO
0)
0)
Jj
w
O
0)
•rl
CO
00
212
-------
taken from position #3 at the vest end of the valley at 0300. The
exposure shows a plume with extensive crosswind diffusion being
transported toward the east.
Sampler 10, at an elevation of 241 m above the stack base,
measured a concentration of 648 ppt, the highest concentration
observed during the entire 10-day experiment. This peak
concentration, normalized by the SFg emission rate, is
3.1 ysec/m3. The calculated value of Hc is 245 m.
SFg tracer gas remained on target mountain during the hour
ending at 0500 despite the fact that winds gradually shifted to
generally easterly. A photograph (Figure 89) taken at 0430 from
Prospect Hill shows plume material overhead. A photograph (Figure 90)
taken at 0615 from the west end of the valley shows plume material
approaching from the east, eventually producing concentrations above
100 ppt at samplers west of the plant. By 0900 concentrations above
200 ppt were again observed in the target area.
Experiment 9 (November 19. 1983) 0000-0100
Figure 91 shows hourly tracer gas concentrations observed during
a windy neutral period. The photograph in Figure 92 illustrates a
persistent, coherent plume from the Tracy stack to the haul road south
of target mountain. The winds in the valley (Figure 93) were all
west-northwesterly at about 8 m/sec.
5,5 Summary of the Preliminary Experiment
The analysis of the three case study experiments and a
preliminary analysis of the other data show the occurrence of stable
plume impingement conditions. SF$ concentrations were observed
during stable conditions on target mountain and beacon hill. Although
the sampler coverage was relatively sparse, concentrations were also
observed during stable conditions in the hill 5610 complex northwest
of the plant.
The highest SFg concentrations were observed on the southwest
corner (samplers 10, 3, 2) of target mountain during Experiment 7.
The elevations of the samplers that captured plume material were a few
meters below the calculated hourly values of Hc. During other hours
of stable plume impingement conditions, plume material was observed to
stay below Hc. In short, it appears that the concept of a dividing-
streamline height will be useful to distinguish flow regimes and to
help simulate observed tracer gas concentration patterns in the Tracy
area.
Drainage winds and katabatic effects were seen to produce
ground-level concentrations on the valley floor in the gorge where the
Truckee River bends to the north. Observer comments, photographs and
acoustic sounder records all suggest the turbulent transport of "old"
plume material from aloft to the valley floor. The fumigation of
oil-fog by drainage winds was also observed on the south side of
target mountain. These katabatic effects were not observed at CCB and
HER and must be accounted for in the final design of the FSPS.
213
-------
o
rfi
•*
O
CO
CO
1
O
S3
cd
0)
u
s
o
•r-l
e
0)
00
0)
u
aO
214
-------
in
\o
o
a)
CO
00
O
s
en
CO
u
3
in
o
&
-------
Data for Hour Ending 0100 on 18-November-1983
0 °0°°
o o .-«••'"•« ....... ........
— kff' 1 0 0
/ °
0 0
o
0
0
Figure 91. Hourly SF& concentrations November 18, 1983.
216
-------
o
o
-p
(0
PI
00
00
u
4)
o
S5
CO
o
u
3
01
o
0
.1-1
Cx
2
217
-------
u
•o
-t-1
f
O
i/l •
H O
O
W-J t-t
O O
C O
O O
.rt o
4J O
.0 -
.t-l fO
b oo
JJ CT>
W t-t
•iH
•o -
oo
O r*
•f-4
£ Ll
Oi 0)
0) O
O 2
00
218
-------
The meteorological measurements depicted very complicated wind
flows during the November experiment. Horizontal and vertical wind
shears were common. These were probably caused by the combined
effects of the complex terrain and migratory anticyclones and cyclones
moving over the area in November. In summer and early fall we expect
the flows in the valley to be more dominated by drainage flows,
principally drainage down the Truckee River. In any event, during the
FSPS there must be a sufficient number of meteorological measurement
systems to provide information on the three-dimensional structure of
winds in the valley.
Finally, the November experiment produced a data base with a wide
variety of dispersion conditions--from windy, neutral cases to stable
plume impingement. About 34 hours were judged good for modeling. The
CTDM concepts will be tested using this data base.
5.6 Plans for the Full Scale Plume Study
The FSPS* will be conducted at the Tracy plant during the period
August 6, 1984 through approximately August 27, 1984. The
participants will include ERT, NOAA WPL, NOAA ARLFRD, and SRI
International. Table 20 provides a schedule of the daily experiments.
The FSPS will commence with two 4-hour shakedown experiments on
the 6th and 7th. These will start around 0300 PST and end around 0700
PST and will be conducted to test the equipment and the experiment
protocol. The shakedown experiments will be followed by 12 10-hour
tracer and simultaneous flow visualization experiments. These will
run primarily during the nighttime hours to capture stable conditions
and are scheduled to take photographic advantage of the full moon on
the llth.
The experimental methods of the FSPS are based on the preliminary
experiment conducted in November and to a large extent are similar to
the methods used and tested at CCB and HBR. The FSPS at the Tracy
station will include:
• Releases of SFg and oil-fog from the 300-ft stack and
CF$Br from various heights on the 150-m tower;
• Operation of 110 tracer gas samplers at 107 locations (four
will be on the 150-m tower);
*The detailed plans are described in the Work Plan for the Full Scale
Plume Study. ERT document P-B876-625, May, 1984.
219
-------
TABLE 20. FSPS SCHEDULE*
M 6 Aug 0300 - 0700**
T 7 Aug 0300 - 0700
T 9 Aug 2000 - 0600
F 10 Aug 2000 - 0600
S 11 Aug 2000 - 0600
W 15 Aug 2200 - 0800
T 16 Aug 2200 - 0800
F 17 Aug 2200 - 0800
M 20 Aug 2200 - 0800
T 21 Aug 2200 - 0800
W 22 Aug 2200 - 0800
S 25 Aug 0000 - 1000***
S 26 Aug 0000 - 1000
M 27 Aug ,0000 - 1000
*Schedule subject to change depending on weather conditions.
**0nly the fan for the 300-ft stack will be operating during the
first 11 experiments.
***The 120-MW unit will be operated during the last three experiments.
220
-------
• Fixed meteorological measurements:
a 150-m tower instrumented at six levels (instruments
will include sonic, propeller and cup-and-vane
anemometers and temperature and radiation sensors),
three 10-m towers instrumented at one level,
one 10-m tower instrumented at two levels,
three monostatic acoustic sounders,
two doppler acoustic sounders, and
two solar-powered electronic weather stations;
• Two tethersondes:
one operated at plume elevation to document
meteorological conditions representative of the
effective source height, and
one operated to measure vertical profiles of
meteorological parameters upwind of the source;
• Two radar balloon tracking systems:
one located near the west end of the valley to measure
the approach flow, and
one located north of the plant to document the winds
near potential impact areas;
• Ground-level tracer gas concentrations;
• Airborne and ground-based lidar measurements;
• Photographs (from five locations) and movies; and
• A command post near the 150-m tower that includes:
real-time display of data from the 150-m and 10-m
towers,
radio base station, and
- facsimile output of weather maps.
The meteorological data will be archived and displayed in real-time by
a system of onsite minicomputers. Real-time information on ambient
meteorological conditions will aid in understanding the dispersion
phenomena and will help the field managers maintain real-time
experimental control. The real-time operations management will be
supplemented by near real-time lidar data, a 48-hour turnaround on the
tracer gas concentrations, and a 48-hour turnaround on the
photographs. Provisions will be made for data analysis in the field
to help guide the experiment. Figure 94 illustrates the layout of the
FSPS.
Oil-fog and SFg will be released from the 91.4-m stack.
CF3Br will be released from one of the three levels (100, 120 or
140 m) on the 150-m tower. The tracer gases will be sampled at about
107 locations in the valley and on the mountains.
ARLFRD will operate 110 samplers during each approximately
10-hour experiment. All samplers will be used to get one-hour
averages. Four samplers will be operated at the 150-m tower. Each
10-hour experiment could produce 1,100 bags and 2,200 concentrations.
The samples will be analyzed and concentration maps produced within
about 48 hours after collection of the bags.
221
-------
illiillll
o
01
On
H
w
222
-------
The sampler network (see Figure 94) was designed principally from
the results of the preliminary experiment that was held in November
1983. It will measure ground-level tracer concentrations during
several plume dispersion conditions:
1) transport along known (observed) plume paths,
2) stable plume impaction,
3) channeling by major terrain features,
4) lee side phenomena,
5) katabatic fumigation,
6) recirculation in the valley, and
7) flat terrain versus hill effects.
Samplers on target mountain and beacon hill and the string of
samplers west of the plant are located to measure concentrations
during events similar to those observed in November and during
previous field experiments. The samplers on target mountain, beacon
hill and on the hill 5085 - hill 5610 complex are located to measure
concentrations during stable plume impaction conditions.
The meteorological towers, the acoustic sounder systems, the
radar wind systems and the tethersondes will provide information to
characterize the three dimensional structure of the wind and
temperature fields in the valley. The 150-m tower, which will be
instrumented with sonic, propeller and cup-and-vane anemometers, will
provide data on the winds and turbulence near the elevation of the
Tracy plume and data to define the meteorological conditions near the
height of the Freon emissions.
The lidar systems and photographs will document the path and
growth of the oil-fog (and coincidental SFg) plume upwind of the
major terrain elements. They will provide information to calculate
-------
SECTION 6
SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS FOR FURTHER STUDY
This Fourth Milestone Report documents the further evolution of
CTDM and presents a detailed mathematical description of the model
components, including those that explicitly account for the effects
that terrain has on plume dispersion. The latest version of CTDM has
been evaluated using an 80-hour subset of the SHIS #1 data base. The
report also describes an empirical modeling approach to the
Hogback Ridge data base. It provides an overview of the preliminary
dispersion experiment that was conducted at the Tracy Power Plant in
November, 1983 and presents plans for the Full Scale Plume Study.
6.1 Principal Accomplishments and Conclusions
The Complex Terrain Dispersion Model (CTDM)
Substantial progress has been made in the development of CTDM as
a method for simulating tracer gas concentrations observed at CCB and
as a practical regulatory model. The central feature of the model
still is its use of the concept of a dividing streamline to separate
the flow into two discrete layers. The upper weakly stratified layer
is handled by the LIFT component and the lower stable layer by the
WRAP component. Both model components include explicit mathematical
expressions that account for the important phenomena that control
dispersion in mountainous terrain:
• streamline contraction in the vertical,
• streamline distortion in the horizontal,
• acceleration of the flow, and
• changes in lateral and vertical diffusivities.
CTDM also includes a method to simulate the transition between the
upper and lower flows and subsequent ground-level concentrations and
to simulate the temporal variability of the Hc interface. The model
also includes new formulations for c*z and
-------
EPA Project Officer. The MDA was evaluated by comparing interpolated
winds and turbulent intensities to those estimated from lidar data and
plume photos. The analysis indicates that the MDA values are
generally appropriate for modeling, although photo and lidar estimates
were substituted in some cases when warranted. A subset of 80 hours
of SE'g and coincidental meteorological data was selected to evaluate
CTDM.
Investigations of Plume Growth
An analysis of SHIS #1 Tower A turbulence intensity measurements
with values of oz estimated from plume photographs and lidar data
suggests the following model of vertical dispersion of elevated
releases in the stable boundary layer:
a t
w _
(1 -H t/2T_ )
LJ
The dispersion time scale TL is given by
TL - a/aw
where the mixing length 9, is
_ _ , _
l I I
n s
and the neutral length scale 8,n and the stable length scale S.s
are given by
J,n = Tzr ; !LS = y2 ow/N
The variance of the vertical velocity fluctuations <*w and the
Brunt-Vaisala frequency N are evaluated at the release height zr.
The constants y and f were derived from surface layer
flux-gradient relationships and were found to have values of 0.27 and
0.32.
Evaluation of CTDM
The current version of CTDM was evaluated by comparing model
calculations to (1) observed SFg concentrations, (2) concentration
estimates based on a flat terrain model, and (3) concentration
estimates based on the COMPLEX I/II plume path assumptions. The
80-hour subset of the MDA was divided into four classes: neutral,
weakly stratified, impingement, and very stable. Performance
statistics were generated for each class as well as the entire
subset. The results show that CTDM simulates the observations better
than the other two approaches. CTDM has a weak overall bias toward
225
-------
overestimating the larger concentrations (m~ = 0.84), the COMPLEX
modeling approach generally overestimates by more than a factor of two
(mg = 0.42), and the flat terrain model generally underestimates by
more than a factor of two (m~ = 2.2). The noise in the CTDM
calculations is lower for the time- and space-paired statistics, which
indicates that CTDM is better able to simulate the observed
distribution of concentrations. These comparisons illustrate the
importance of including Hc and the stagnation streamline in a
modeling framework such as CTDM.
The performance of CTDM in each of the four classes indicates
that the model tends to overestimate the larger observed
concentrations for Hc greater than 1.25 times the release height
(the two "more stable" classes). These classes include the hours in
which the largest concentrations (scaled by the emission rate) were
observed at CCB. CTDM tends to underestimate in the other two "less
stable" classes, producing estimates that are 70% to 80% of the larger
observed concentrations. These results indicate that more work is
needed to further improve CTDM performance within each of the classes.
Model Performance for Various Quantities of Onsite MeteoroloRical
Data
The ability of CTDM to simulate the CCB observations with various
quantities of onsite meteorological data was evaluated by constructing
three alternate model input data sets: (1) 5-minute meteorological
data contained in the MDA; (2) an hourly data set constructed from
one-hour wind and temperature data measured at 10-m and 150-m; and (3)
an hourly data set constructed from wind and temperature data measured
at 10-m only. The results suggest that the 5-minute simulations
improve some individual hours but do not substantially change the
overall performance statistics. The results from the simulations
using the simplified meteorological input are inferior to those based
on the full MDA even though MDA wind directions were used for each.
The conclusion of this analysis is that onsite measurements of
turbulence intensity near release height and detailed vertical
profiles of wind and temperature are essential for accurately
simulating concentrations at a complex terrain site such as CCB.
Modeling the Hogback Ridge Data Base
To help understand the phenomena that control disperssion at the
ridge site a simple empirical modeling approach was taken. A model
was constructed by modifying the effective plume height as a function
of Hc. Model simulations were performed using a subset of the
SHIS #2 CF3Br data base. These empirical model calculations were
compared to observed CF3Br concentrations and to calculations made
with a flat terrain model and a model based on the half-height plume
path assumption. The empirical model performed better than the other
two models. Furthermore, for both CCB and HBR the flat terrain model
underestimated concentrations while the COMPLEX (or half-height) model
overestimated concentrations by roughly a factor of two.
226
-------
Fluid Modeling Facility Simulations
A series of tows was conducted in the EPA FMF salt-water-
stratified towing tank using, as the basic hill shape, the
fourth-order polynomial hill as used by Hunt and Snyder (1980).
Linear density gradients were established in the tank and the hill was
towed at a speed such that the Froude number Fr was 0,5. Since the
density gradient was linear and Fr = 0.5, the dividing-streamline
height was 0.5h. After that series was completed, a second series of
tows was conducted wherein the entire model was raised out of the
water to the point where the water surface was precisely at the
dividing-streamline height, i.e., the water surface was at half the
hill height. The model was towed at the same speed as in the
full-immersion tows, so that the Froude number with this now
half-height hill was unity, and all streamlines passed over the hill
top. The flat water surface thus forced a flat dividing-streamline
surface. The resulting surface concentration patterns were then
compared with the corresponding full-immersion tows. These
simulations, which are described in Appendix A, were performed to
answer the question: how good is the assumption of a flat
dividing-streamline surface?
The results suggest that this assumption is a reasonable
approximation to make, at least with regard to predicting the
locations and values of maximum concentrations and areas of coverage
on the windward side of the hill. When the stack heights are
relatively close to the dividing-streamline height, the lee-side
concentrations are also predicted reasonably well. These results were
used directly in the formulation of CTDM.
The FMF staff also conducted wind-tunnel studies (Appendix B) to
investigate the influence of an idealized three-dimensional hill on
ground-level concentrations from upwind sources in a neutral
atmospheric boundary layer and to locate the source locations where
this influence is the greatest.
The presence of the hill was found to influence the dispersion of
the plume to increase the maximum concentration in three ways. For
low sources, at moderate distances from the hill, the reduction in
mean wind speed allows the plume to reach the ground surface closer to
the source, producing higher concentrations than in the absence of the
hill. Plumes from higher sources can be thought of as being
intercepted by the hill. That is, the hill penetrates the plume where
the concentrations are greater than those that would occur at ground
level farther downstream over flat terrain. For yet higher sources,
the streamline convergence over the hilltop and the corresponding
downward flow in the lee of the hill brings the plume to the ground
more rapidly than over flat terrain. The maximum concentration for
these three regimes occurs upwind of the hilltop, near or on the
hilltop, or downwind of the hilltop, respectively. Terrain
amplification factors ranged from near 1.0 to 3.63. The region of
source locations that produced an amplification factor of 1.4 or more
extended to an upwind distance of 14 hill heights.
227
-------
6.2 Recommendations for Further Study
6.2.1 The SHIS #2 Data Base
Further refinements are recommended in screening and adjusting
the meteorological data obtained at the three primary towers at HER.
Once this is completed, the Quality Assurance Report for SHIS #2 can
be completed, and the modelers will have better guidance in
interpreting the data.
In addition to the sonic anemometer data used in the preliminary
modeling reported in this Milestone Report and the acoustic profile
data, the meteorological measurements from the 150-m and 30-m towers
(towers A and B) are essential to investigating the flow properties
below Hc. The temperature patterns and changes in the turbulence
properties will be investigated. A comparison of profiles of these
data within 30 m of the surface at these two towers should essentially
provide a description of the flow field within the "blocked" region
below Hc for nearly all releases below Hc.
For releases above Hc, comparison and integration of
tethersonde data, acoustic profile data, and measurements from
towers A and C are particularly important. Also, because of the
oil-fog plume rise so characteristic of many of the experiment-hours,
analysis of the photographs and the lidar scans of the visible plume
will be pursued to document and model the plume rise, and to model the
initial growth of the visible plume. Analysis of lidar scans and
photographs will also document patterns of streamline deflections in
the flow over the crest of HBR, and possibly the distribution of plume
material near the surface.
6.2.2 CTDM
Application of CTDM to the SHIS #1 data base has reached a stage
where the modeling framework is largely complete. Some effort will be
devoted to investigating the LIFT/WRAP transition zone, and fluid
modeling studies at the EPA FMF will help refine the formulation of
the flow in this region. In addition, it appears that the performance
of the model will be improved by including directional wind shear,
which so far has been ignored. Also ignored has been the variation of
the turbulence with height (away from the plume centerline).
Aside from these enhancements to CTDM, most of the work remaining
in applying CTDM to the SHIS #1 data base will focus on the
terrain-effects factors T^, T^, Tu, and T^. The present formulation
will be compared with computations performed using thin aerofoil
theory as well as results from fluid modeling simulations,. Also,
particular model parameters will be selected for optimization in order
to infer the terrain effect description which best matches the
observed concentrations and meteorological data. Presumably,
alternate choices of the terrain-effect parameters will be needed to
improve the model performance. The specification of T^ (either
T^z or T^y) is a good example. Our present calculations are made
with Tj_ = 1. Rapid distortion theory will be used to infer non-zero
values for T^. Optimized values of T^ will be compared with those
computed from the rapid distortion theory.
228
-------
Application of CTDM to other sites will require some
generalization of the code. A good test-case for such a
generalization is HER. The model will be altered so that the present
CTDM "assumptions" can be tested at HER. The experiments at the Tracy
Power Plant (FSPS preliminary and FSPS) also need to be considered in
structuring the general CTDM. New algorithms will surely be needed in
some circumstances, but we also need to know when the detail in the
present CTDM is most needed. A sensitivity analysis will help define
the conditions in which the terrain effects contained in CTDM have the
greatest and least impact on the magnitude of the concentration
estimates compared to estimates made by means of a simple
"flat-terrain" model. Key input parameters could then be designed to
signal a "CTDM" calculation or a "FLAT" calculation when the model is
applied in a regulatory permitting mode.
The development of a theoretical rather than an empirical model
for HBR will receive much attention in the next year. We first need
to identify the circumstances in which CTDM does as well as the
empirical model, and circumstances in which CTDM does much worse. The
cases in which CTDM fails must be analyzed in detail to develop a
rationale (theory) for modifying CTDM. In this way CTDM will increase
its range of applicability, and the need for empiricism will be
reduced. This process will begin by applying CTDM to the subset of 35
hours modeled in this Milestone Report. Although the meteorological
data are incomplete, we expect to learn something from this
preliminary modeling. We expect to develop a framework to simulate
the effects of an unsteady, "blocked" flow upwind of a ridge.
6.2.3 The FSPS Data Base
The FSPS at Tracy Power Plant has been designed to provide a data
base that will augment our understanding of how a plume in large-scale
stable flow interacts with topography and local slope flows to reach
the surface. The phenomena observed at this site will be described
and documented, and the relationship between these phenomena and
observed ground-level concentrations will be investigated. This
investigation will include comparisons with both the SHIS #1 and #2
results. We expect that the similarities among the FSPS and the SHIS
observations will identify those aspects of the FSPS data base that
can be modeled with CTDM directly. The dissimilarities will focus
attention on aspects that will require modification to CTDM.
229
-------
REFERENCES
Briggs, G.A. 1973. Diffusion Estimation for Small Emissions, ATDL
Contribution File No. 79, Atmospheric Turbulence and Diffusion
Laboratory.
Briggs, G.A. 1975. Plume Rise Predictions. Lectures on Air Pollution
and Environmental Impact Analyses. AMS, Boston.
Brighton, P.W.M. 1978. Strongly Stratified Flow Past Three
Dimensional Obstacles. Quarterly Journal of the Royal
Meteorological Society. 104; 289-307.
Businger, J.A. 1973, Turbulent Transfer in the Atmospheric Surface
Layer. In Workshop on Micrometeorology, AMS, Boston, 67-100.
Csanady, G.T. 1973. Turbulent Diffusion in the Environment.
D. Reidel Publishing Co., Dordrecht, Holland.
Deardorff, J.W. and G.E. Willis 1975. A Parameterization of
Diffusion into the Mixed Layer. J. Appl. Met.. 14: 1451-1458.
Drazin, P.G. 1961. On The Steady Flow of a Fluid of Variable Density
Past an Obstacle. Tellus. 13.: 239-251.
Gifford, F.A. 1980. Smoke as a Quantitative Atmospheric Diffusion
Tracer. Atmospheric Environment. 14: 1119-1121.
Holzworth, G.C. 1980. The EPA Program for Dispersion Model
Development for Sources in Complex Terrain. Second Joint
Conference on Applications of Air Pollution Meteorology, New
Orleans, LA. AMS, Boston.
Hunt, J.C.R., and R.J. Mulhearn 1973. Turbulent Dispersion from
Sources Near Two-Dimensional Obstacles. J. Fluid Mech., 61:
245-274.
Hunt, J.C.R. and W.H. Snyder 1980. Experiments on Stably and
Neutrally Stratified Flow Over a Model Three Dimensional Hill.
J. Fluid Mech.. 96: 671-704.
Hunt, J.C.R. 1981. Diffusion in the Stable Boundary Layer.
Atmospheric Turbulence and the Pollution Modelling. F.T.M.
Nieustadt and H. Van Dop (Ed.). D. Reidel, Holland.
230
-------
REFERENCES (Continued)
Hunt, J.C.R. 1982. Diffusion in the Stable Boundary Layer. In
Atmospheric Turbulence and Air Pollution Modeling. D. Reidel
Publishing Company, Dordrecht, Holland, 231-274.
Hunt, J.C.R., J.C. Kaimal, J.E. Gaynor and A. KorrelL 1983.
Observations of Turbulence Structure in Stable Layers at the
Boulder Atmospheric Laboratory. In Studies of Nocturnal Stable
Layers at BAO. Report Number Four, Jan. 1983. Available from
NOAA/ERL, Boulder, CO 80303, U.S.A., 1-52.
Kapsha, T.P., et al. 1976. Behavior of S02 Plumes from Tracy and
Fort Churchill Generating Stations Under Stable Atmospheric
Conditions in the Vicinity of Complex Terrain. Westinghouse
Electric Corporation, Pittsburgh.
Lavery, T.F., A. Bass, D.G. Strimaitis, A. Venkatram, B.R. Greene,
P.J. Drivas, and B.A. Egan 1982. EPA Complex Terrain Model
Development: First Milestone Report - 1981. EPA-600/3-82-036,
Research Triangle Park, NC. 304 p.
Lauwerier, H.A. 1954. Diffusion from a Source in a Skew Velocity
Field. Appl. Sci. Res.. 4, p!53.
Overcamp, T.J. 1983. A Surface-Corrected Gaussian Model for Elevated
Sources. J _._ of Climate and Applied Met. . 22: 1111-1115.
Pasquill, F., and F.B. Smith 1983. Atmospheric Diffusion, 3rd
Edition. Ellis Horwood Ltd, England.
Pearson, H.J., J.S. Puttock and J.C.R. Hunt 1983. A Statistical
Model of Fluid Element Motions and Vertical Diffusion in a
Homogeneous Stratified Turbulent Flow. JA Fluid Mech., 129:
219-249.
Riley, J.J., Liu, H.T. and Geller, E.W. 1976. A Numerical and
Experimental Study of Stably Stratified Flow Around Complex
Terrain. EPA Report No. EPA-600/4-76-021, Res. Tri. Pk., NC, 41p.
Sheppard, P.A. 1956. Airflow Over Mountains. Quart. J. R. _Met«jor.
Soc., 82: 528-529.
Stiyder, W.H. , R.E. Britter and J.C.R. Hunt 1980. A Fluid Mo_deling
Study of the Flow Structure and Plume Impingement on a
Three-Dimensional Hill in Stably Stratified Flow. Proc. Fifth
Int. Conf. on Wind Engr. (J.E. Cermak, ed.), 1: 319-329, Pergamon
Press, NY, NY.
231
-------
REFERENCES (Continued)
Snyder, W.H. and J.C.R. Hunt 1984. Turbulent Diffusion from a Point
Source in Stratified and Neutral Flows Around a Three Dimensional
Hill; Part II - Laboratory Measurement of Surface Concentrations.
Submitted to Atmospheric Environment.
Uthe, E.E. and B.M. Morley 1984. Alpha-1 Observations of Plume
Behavior for PMV&D Tracy Site. Electric Power Research Institute,
Palo Alto, California.
Venkatram, A. and R. Paine 1984. Development of a Model to Estimate
Dispersion of Elevated Releases in the Shear-Dominated Boundary
Layer. ERT Report No. P-B846-500. Prepared for the Maryland
Power Plant Siting Program.
Venkatram, A., D. Strimaitis, D. DiCristofaro 1984. A Semiempirieal
Model to Estimate Vertical Dispersion of Elevated Releases in the
Stable Boundary Layer. Atmospheric Environment. 18: 923-928.
232
-------
APPENDIX A
STABLE PLUME DISPERSION OVER AN ISOLATED HILL
RELEASES ABOVE THE DIVIDING-STREAMLINE HEIGHT
233
-------
STABLE PLUME DISPERSION OVER AN ISOLATED HILL
Releases above the Dividing-Streamline Height
by
William H. Snyder *
and
Robert E. Lawson, Jr. *
Meteorology and Assessment Division
Atmospheric Sciences Research Laboratory
U.S. Environmental Protection Agency
Research Triangle Park, NC 27711
December 1983
*On Assignment from the National Oceanic and Atmospheric Administration,
U.S. Department of Commerce.
23A
-------
ABSTRACT
A series of tows was conducted in a stably stratified salt-water
towing tank wherein the density gradient was linear and the
dividing-streamline height was half the hill height. Effluent was
released at three elevations above the dividing-streamline height.
Pairs of tows were made such that, in one tow, the hill (upside-down)
was fully-immersed in the water and the towing speed was adjusted to
provide a "natural" dividing-streamline surface. In the second tow of
the pair, the hill was raised out of the water to the point where only
the top half of the hill was immersed, thus, forcing a flat dividing-
streamline surface, while all other conditions were maintained
identical. Concentration distributions were measured on the hill
surface and in the absence of the hill. Concentration distributions
from each pair of tows were compared to ascertain effects of an
assumed flat dividing-streamline surface as is used in some
mathematical models. The results suggest that the assumption of a
flat dividing-streamline surface is a reasonable approximation to
make, at least with regard to predicting the locations and values of
maximum concentrations and areas of coverage on the windward side of
the hill. When the stack heights are relatively close to the
dividing-streamline height, the lee-side concentrations are also
predicted reasonably well. The apparent cause of the relatively poor
agreement between lee-side concentration patterns in the higher stack
cases is the presence of a hydraulic jump at the downwind base of the
hill in the full-immersion case which was absent in the half-immersion
case.
235
-------
1. INTRODUCTION
The structure of strongly stratified flows over a three
dimensional hill has been envisioned as composed of two layers: a
lower layer of essentially horizontal flow wherein plumes from upwind
sources impinge directly on the hill surface, and an upper layer
wherein plumes from upwind sources may pass over the hill top. This
basic concept was suggested by theoretical arguments of Drazin (1961)
and Sheppard (1956) and was demonstrated through laboratory
experiments by Riley et al. (1976), Brighton (1978), Hunt and Snyder
(1980), Snyder et al. (1980) and Snyder and Hunt (1983). Complex
terrain diffusion models utilizing this concept have been developed by
Hunt, Puttock and Snyder (1979) and, more extensively, by Lavery
et al. (1982), Strimaitis et al. (1982), and Lavery et al. (1983).
The basic parameter characterizing the flow structure is the Froude
number, F (see Snyder and Hunt, 1984).
We are concerned in this report with the upper-layer flow. The
basic assumption in this approach is that, in strongly stratified
flows (0
-------
hypothesis (balance of kinetic and potential energy), terminates on
the surface at the top of the hill, again suggesting that the
dividing-streamline surface cannot be flat. On the other hand,
concentration measurements on hill surfaces (Snyder and Hunt, 1984)
suggested that a flat.surface approximation may yield reasonable
estimates. From a practical viewpoint, the mathematical models are
vastly simplified if such an assumption yields reasonable estimates of
surface concentration. Hence, we attempt to answer the question not
from a detailed analysis of the shape of such a dividing streamline
surface, but from the more practical comparison of surface
concentration patterns.
A series of tows was conducted in the salt-water-stratified
towing tank using, as the basic hill shape, the fourth order
polynomial hill as used by Hunt and Snyder (1980). Linear density
gradients were established in the tank and the hill was towed at a
speed U such that the Froude number F ( = U/Nh, where N is the
Brunt Vaisala frequency) was 0.5. Effluent was released at heights of
Hs = 0.6h, 0.7h and 0.8h, and the resulting hill-surface
concentration patterns were measured. Since the density gradient was
linear and F-0.5, the dividing-streamline height was also 0.5h (Hc/h
= 1 - B'); since the effluent was released above Hc, the plumes did
not impinge directly on the hill, but instead pollutants reached the
hiil surface by the combination of streamline displacement and
diffusion. After that series was completed, a second series of tows
was conducted wherein the entire model (hill, baseplate and stack, as
a unit) was raised out of the water* to the point where the water
surface was precisely at the dividing-streamline height, i.e., the
water surface was at half the hill height. The model was towed at the
same speed as in the full-immersion tows, so that the Froude number
with this now half-height hill was unity, and all streamlines passed
over the hill top. The flat water surface thus forced a flat
dividing-streamline surface. The resulting surface concentration
patterns were then compared with the corresponding full-immersion tows
to ascertain the effects of a flat dividing-streamline surface.
*The model is routinely mounted upside down such that the baseplate is
submerged a few millimeters below the water surface. In discussion
of flow structure and plume behavior, however, we discuss the results
as if the model were right-side up.
237
-------
2. EXPERIMENTAL APPARATUS AND TECHNIQUES
Most of the details of the experimental apparatus and techniques
were given by Hunt and Snyder (1980) and in a laboratory report by
Hunt, Snyder and Lawson (1978). The basic method of making
concentration measurements was described by Snyder and Hunt (1984).
Only a brief overview is given here, but changes in the techniques and
apparatus as well as special features of these experiments are
described in detail.
A fourth-order polynomial hill (z(r) = h/(l+(r/L)4)) of height
24.3 cm was used in a stratified towing tank. The tank, 1.2 m in
depth, 2.4 m in width and 25 m in length, was stably stratified with
layered mixtures of salt water. This dye mixture was emitted at four
times the isokinetic rate from a bent-over "stack" of 0.635 cm o.d.
The stack exit was located 84.8 cm (3.5 h) upstream of the hill
center. The non-isokinetic effluent release rate was used to obtain a
realistic plume size and shape, i.e., in the isokinetic releases used
previously by Snyder and Hunt (1984), because of the nonexistent
approach-flow turbulence, the plumes tended to be exceptionally thin
and narrow. With the four-times-isokinetic rate, a weak but turbulent
jet was formed at the stack exit. This jet grew in size with
downstream distance to provide a plume with dimensions not
insignificant in comparison with, say, the hill height and hence, a
more realistic simulation of a typical atmospheric situation. The
maximum jet velocity at the hill center was estimated to be about
2 cm/s relative to the hill (Townsend, 1956), i.e., small compared
with the towing speed. Characteristics of the plume in the absence of
the hill were measured with horizontal and vertical rakes of sampling
tubes, and are also presented herein.
One hundred sampling ports were fixed on the hill surface,
distributed as shown in Figure A-l. In some preliminary tows, the
sampling ports (2.4 mm o.d.) protruded 2.5 mm above the smooth hill
surface, but as the plumes were spread broadly to cover most of the
hill surface (but very thinly in the direction normal to the surface),
narrow, clear wakes were observed down-stream of the protuberant
sampling tubes. This disturbance caused strong reductions in
concentrations measured at ports directly downstream from others,
e.g., ports along the 0° line (Figure A-l). The sampling ports were
then cut to be flushed with the hill surface, but some amount of
interference was still observed in the concentration distributions,
apparently due to the withdrawal of sample fluid through the ports.
After considerable experimentation, a final configuration was found
that displayed no interference: the ports were raised to the original
2.5 mm above the smooth hill surface and the hill was covered with
238
-------
312. 90°
-34-0,
180° |
.113.U211HimQ91.Q81lIl7iQei05.10ii.lQ3.
FIJOW '
DIRECTION \
2C.?U'ii01.
12. 101. 0 501. 502. 5Q3.5Qu.5055.065Q75a£S.0951Q5U512.513.5'.i-l.
801. 601. I
701> / I
802. 602.
612.
-45°
E
Pfl
POLYNOMIflL
HILL
FLUID
MODELING
SECTION
Figure A-l. Port locations and numbering system for lift study.
(Precise location is at decimal point. Dashed circle
marks half the hill height.)
239
-------
sharp-edged gravel of grain size (longest dimension) 2 to 4 mm. Also,
the sample withdrawal rate was reduced by a factor of about 4, to
13 cm-Vmin. This corresponds to a sampling "stream tube" diameter
of 1.3 mm'at the typical tow speed. The rough surface provided the
additional benefit of eliminating a viscous sublayer on the hill
surface, hence, avoided problems of molecular diffusion through this
viscous sublayer.
The concentration of dye in the collected samples were analyzed
on a Brinkman Model PC-600 probe colorimeter. The fiber-optics probe
was immersed sequentially into the sample test tubes. The wavelength
used was 570 nm. The output of the colorimeter, a voltage related to
the opacity of the solution being tested, was fed to a POP 11/44
minicomputer, where it was converted to a concentration in percent
dye. The conversion utilized a calibration curve formed by recording
the output voltage versus concentration for at least 12 "standards"
which consisted of accurately known dilutions of the same dye used for
the effluent source. A "Beer's Law" type of curve was best-fit to the
standards for use in converting the voltage from the unknown sample
into a dye concentration. Although the instrument required care in
use (by frequently checking the "zero transmittance" and "zero
absorbance" controls), frequent checks of the calibration showed
excellent repeatability. A typical calibration curve is shown in
Figure A-2. The concentration of dye in the effluent was occasionally
adjusted such that the majority of samples would be within the most
reliable range of the calibration curve, i.e., the vast majority of
samples analyzed were in the range of 0.005 to 0.25%.
The sample lines leading from the ports on the hill surface (or
rake) to the sample test tubes had to be filled prior to the beginning
of a tow in order for the vacuum-sampling system to work properly.
Because of the reduced sample flow rate, the volume of dye-free water
stored in the lines (and which was thereafter drawn into the test
tubes along with the dye samples) was approximately 25% of the total
sample collected. A first-order correction to the measured
concentrations was made to account for this dilution of the samples.
However, the precise volume of sample collected in each test tube
varied slightly from one tube to the next, so that the accuracy of the
concentration-measurement system is estimated to be in the range of
to 15%.
240
-------
A fOLCRL. it C 2. 1) ^.I^r 3:UD1 . POLrHFLL.31.
" 3CERFT. '4 t ?., 5) QTCON7 MGV-25-77
27-DEC-83
IP fOCD DT£.(^-J S3. 3231 . BRINKMHM . IPCGOOi
»&r!/SEC,OH..H2CCEP.=:0.5aCM.?'r-V:, «••
-------
3. PRESENTATION AND DISCUSSION OP RESULTS
3.1 Plumes in the Absence of the Hill
Vertical and horizontal concentration profiles measured at the
downstream position of the center of the hill (but in its absence) are
presented in Figures A-3 and A-4, respectively. Profiles were
measured for each of three stack heights, 0.1 h, 0.2 h and 0..6 h. For
the two lowest stack heights, the baseplate was raised out of the
water; for the highest stack height, it was submerged as it was for
the full-immersion tows with the hill (approximately 6 mm).
Concentrations have been nondimensionalized as x = CUooh^/Q,
where C is the dimensional dye concentration, Um is the towing
speed, h is the hill height, and Q is the volume flow rate of dye in
the effluent. The calculated plume parameters are listed in
Table A-l. The two higher-level plumes are essentially identical and
are very nearly Gaussian in shape. (Best-fit Gaussian curves; are
shown for comparison in the figures. The parameter Xmxg shown in
Table A-l is the maximum concentration in the best-fit Gaussian
distribution with the same standard deviation as the measured data.)
The lower-level plume, however, has clearly diffused to the water
surface (Figure A-3); a reflected-Gaussian profile assuming ain
effective stack height of 0.09 h appears to fit the data quite well.
The standard deviation of the vertical distribution az (calculated
assuming perfect reflection at the water surface), however, is 40%
larger than those of the higher level plumes. The standard deviation
of the lateral distribution cry of the lower level plume is also
somewhat larger (10%) than that of the upper level plumes. The
probable cause of this increased diffusion is the reduced stability
near the water surface, as shown in Figure A-5. Ideally, the density
gradient would have been perfectly constant all the way to the water
surface, so that the plume released at the lowest level (0.1 h) in the
half-immersion case would have been essentially identical in shape to
that released at the upper level (0.6 h) in the full-immersion case.
However, as will be shown later, these differences in approaching
plume shapes appeared to have little effect on the concentration
distributions measured on the hill surface.
3.2 Distributions on the Hill Surface
Figure A-6 shows side and top views of the plume surmounting the
hill, in this case with the hill fully immersed and the effluent
emitted at 0.6 h. These photographs show the plume to be deformed to
cover essentially the entire surface of the hill above half the hill
242
-------
LIFT .25 t 3. 1) LIFT STUDT, NO HILL2. VERT. PROF.. STK. HT.-2...
09-DEC-83 DTE=U5.6ML/MIN.DIL=865.8ML/MIN,HCRIT=. 5M, ci i: X P. . L'LK, T-0, 2= J 4. 58CMPBLCP
SCflL 1,0000, 0.000, 2*. ViO, 0.000
Q LIFT .21 (3, 1) LIFT STLIDT. NO HILL2, VERT. PROF. . 5TK. HT. =u. 8K K. i /?= 1 S. 77CK/SEC
07-DCC-83 DrC^U7. 7ML/MIN, DIL-905. 8ML/MIN, IICRIT-G. C-!!. n-vT: X-OH. GCt-'. T-0, Z-lPCLCfiL.
-------
A LIFT .26 (1. 3) LIFT STUDT.NO HI LL2. VERT. PROF. . STK. HT . -2 A3CK , T /:,- 1 :,. i Cf./L^C
scfiL ' 211.3000. o.'ooo, ' i'.'ooo', '"o'.ooo"""
LIFT STUDT.NO HILL2. LflT. PROF.,STK.HT.=U.96CM.T/S=15.77CK/SEC
DTE-i47.7ML/MIN,OIl.-905.8ML/MIN,IICRIT-0.5i:.riKt:i X-OU.CCM.T-D.2-17.
SCfiL 2U.3000. 0.000, 1.000, 0.000
O LIFT . ?li [ 1. 3) LIFT STUDT.NO HILL.VERT.PROF.,STK.HT.=1U.58CK,T/S=15.83CM/SfC
08-DEC 30
D LIFT .22 I 1. 3)
07-DCC-83
X
u
10 -
5
-.03
Y/H
Figure A-4. Horizontal concentration profiles measured at the
downstream position of the center of the hill, but in
its absence.
244
-------
File
No.
21
22
23
24
25
26
Horizontal/
Vertical
V
H
V
V
V
H
Hs/ll
0.2
0.2
0.6
0.6
0.1
0.1
Xmx
25.8
23.6
23.4
21.6
21.1
22.4
a
(cm)
1.07
2.36
1.07
2.28
1.50'
2.55
inixg
28.7
27.2
25.3
25.1
22.2
24.0
(a) Calculated assuming perfect reflection at the water surface and an
effective stack height of 2.2 cm (0.09 h).
245
-------
PL5RS .30 [3.
09-DEC-83
POLTHILL LIFT STUDT , WHTCR IIT . = 1 00. 8CK, MCTTLCH (PLE'OO)
TIME: 0610
SCRL 1.0000, 0.000. 21.500, 0.000
.8
.7
.6
.5
Z
\
N
«-»
1.06 1.07
1.08 1.09
SPECIFIC GRAVITY
1.1
1.11
Figure A-5. Density profile measured prior to Tow #25.
(Concentration distributions shown as lower plume in
Figures A-3 and A-4.)
246
-------
Figure A-6. side and top-views of plume released at 0.6 H.
F = 0.5, fully submerged.
247
-------
height and yet to be quite thin in the direction normal to the hill
surface (the dotted line on the photographs marks half the hill
height). Notice that, on the lee side, the plume sweeps to much lower
elevations, so that we may expect to see substantial surface
concentrations on the lee side below half the hill height.
Figure A-7 presents a scatter diagram showing the repeatability
of surface concentration measurements from one tow to the next. Tn
this case, both tows (numbers 10 and 11) were made with the hill
fully-immersed and the stack height fixed at 0.6 h. From the; diagram,
it may be seen that, for any port, repeatability is well-within a
factor of two. At large values of concentration, repeatability is
generally within 10%. The repeatability is less good at the lower
concentrations. A few points with the largest percentage differences
between the tows are marked with the port number (see Figure A-l);
these show that the largest differences occurred below half the hill
height, i.e., on the lower edges of the plume.
Figures A-8 and A-9 show the concentration distributions measured
on the hill surface in the full-immersed and half-immersed cases,
respectively. The most obvious difference between the two ca,5es is
the absence of lee-side concentrations below half the hill height in
the half-immersed case. Of course, in the half-immersed case,
concentrations at positions below half the hill height are zero,
because that portion of the hill is outside the water. In the
fully-immersed case, the plume diffused to some extent below half the
hill height around the upwind side, but also, as mentioned eai-lier,
this plume "hugged" the hill surface as it was swept down the lee side
to much lower elevation than the release height. Just beyond the
downwind base of the hill (not shown on the photographs), the flow
appeared to separate from the surface and rise abruptly in a hydraulic
jump (see Hunt and Snyder, 1980).
Side and top views of this plume in the half-immersed case are
shown in Figure A-10. These are to be compared with the
full-immersion case of Figure A-6. The comparison shows similar
behavior as discussed in the above paragraph.
Figure A-ll presents a scatter diagram comparing, on a
port-to-port basis, the surface concentrations measured in the half-
and fully-immersed cases. Measurements at points below half the hill
height are not included here because, of course, in the half-immersion
case, these ports were out of the water. Within the region of large
concentrations, the two cases compare quite favorably, the
half- immersed case yielding concentrations approximately 10 to 20%
larger than the fully-immersed case. In the region of low
concentration, quite large differences occur (worst case, a factor of
10). However, a close examination shows that in all cases where
concentrations differed by more than a factor of 2, the port locations
were very close to half the hill height i.e., either at 0.505 h or
0.59 h (port numbers are marked on Figure A-ll for those cases where
concentrations differed by more than a factor of 2).
248
-------
L1FOIF. 7 ( 3. 2) POLTHIL LIFT STUDY. COL 1. PGfU NO.: COL.2 TO U. CGMI. PRO1; i ITT.010
22-NOV-83 TO .012
100
I
s
o
cc
u.
X
O
1 10
CHI FROM TOW 11
100
Figure A-7. Repeatability of surface concentration measurements.
F « 0.5. HS/H =0.6, fully submerged.
249
-------
Figure A-8. Concentration distributions measured on the hill
surface. F = 0.5. HS/H = 0.6, fully submerged.
250
-------
Figure A-9. Concentration distributions measured on the hill
surface. F = 0.5, HS/H =0.6, half submerged.
251
-------
Figure A-10. Side and top views of plume released at 6.6 H.
F = 0.5, half submerged.
252
-------
LIFDIF. 7 ( li. 3) PGLTHIL LIFT STUDT. COL 1, PCRT NO.: COL.2 TO 4. CGNC. FfiCiM LIFT.OiC
22-NCV-83 TO .012
100
1=1
tu
o
-------
Note that the maximum hill-surface concentration is in the range
of 25 to 30, whereas the maximum in the plume in the absence of the
hill was in the range of 20 to 25. This suggests that even though the
plume shape was highly contorted through the streamline deformation
process, the basic diffusion process was limited by the presence of
the surface; this "reflection" at the hill surface resulted in a
maximum surface concentration that was 20 to 25% larger than that
observed at the center of the plume in the absence of the hill.
However, as mentioned earlier, the accuracy of the concentration-
measurement system was estimated to be in the range of +10 to 15%, so
that the above estimates of increased concentration should be taken
with caution.
Figures A-12 and A-13 show the concentration distributions
resulting from sources elevated at 0.7 h. As was the case with a
source height of 0.6 h, the contours on the windward side of the hill
are roughly circular, but, in this case, somewhat more elongated in
the streamwise direction. Again, of course, no concentrations were
observed below half the hill height in the half-immersion case; in the
full-immersion case, the plume "hugged" the lee side of the hill, but
was spread less broadly in the crosswind direction and lifted off the
surface into the hydraulic jump somewhat sooner.
Also, as was the case at a source height of 0.6 h, the maximum
concentration is about 10% larger in the half- than in the
full-immersion case. The location of the maximum, however, has moved
from the windward side to the lee side of the hill.
A scatter diagram comparing concentrations on a port-by-port
basis for the half- and fully-sumberged hills with a source height of
0.7 h is shown in Figure A-14. The correspondence between the full
and half-depth immersions is not as good as was the case with a source
height of 0.6 h. For the large concentrations, the agreement is
excellent, but for small concentrations, the scatter is quite large.
Port numbers where the concentrations differed by more than a factor
of two are marked on the figure. These show that the comparisons were
poor only in the elevation range of 0.5 h to 0.67 h, i.e., n&ar the
water surface in the half-immersion case. In the half-depth case,
concentrations were larger on the windward line of ports, apparently
related to the wider vertical diffusion of the plume to the water
surface, whereas concentrations elsewhere round the sides of the hill
were considerably smaller than in the full-depth case.
Figures A-15 and A-16 show the concentration distributions
resulting from sources elevated at 0.8 h. In this case, the
correspondence between the two distributions is not nearly as good as
was the case at the lower stack heights. whereas the locations of the
maximum concentrations are approximately the same, the values in the
full-immersion case are 55% larger and the area of plume contact on
the hill surface is nearly quadrupled over the half-immersion case.
Note that the location of the maximum concentration is well to the lee
side of the hill and that the area of plume contact is much smaller
than for the lower stack heights.
254
-------
Figure A-12. Concentration distribution measured on the hill
surface. F = 0.5, HS/H =0.7, fully submerged.
255
-------
\
\
\
\
\
Figure A-13. Concentration distribution measured on the hill
surface. F = 0.5, HS/H =0.7, half submerged.
256
-------
LITDir. B ( 2, 3) POLTHIL LIFT STUDY. COL 1, PORT NO. ; COL 2 TO 7, CONIC. T'l'ul', LITT.OiJ.
Ol-DEC-83 .017. .014. .015. .016. RND .018. RESPECTIVELY.
100
o
UJ
O
OC
Ul
5
0)
D
09
3
O
X
o
1 10
CHI, TOW 13, HALF SUBMERGED
100
Figure A-14.
Comparison of concentrations. Half submerged versus
fully submerged. F = 0.5, HS/H =0.7.
257
-------
\
\
Figure A-15. Concentration distributions measurea on the hill
surface. F = 0.5, HS/H = 0.8, fully submerged.
258
-------
\
\
Figure A-16. Concentration distributions measured on the hill
surface. F = 0.5, HS/H = 0.8, half submerged.
259
-------
A scatter diagram comparing these two cases is presented in
Figure A-17. The correspondence between the full- and half-depth
immersion cases is clearly poor, with the half-depth immersion tow
exhibiting surface concentrations that are generally less than half
those from the full-depth tow. Both of these tows were repeated and
the repeatability was found to be excellent, i.e., similar to that
shown in Figure A-7 for the lower stack height of 0.6 h.
The poor correspondence in this case may be partially understood
by examining the side- and top-view photographs of the plume as shown
in Figures A-18 and A-19. The top-views show that, in the half-depth
tow, the plume was much narrower on the lee side of the hill. The
side views show that, in the full-depth tow, the plume "hugged" the
hill surface on the lee side to elevations considerably lower than
half the hill height, then lifted off the surface and rose somewhat in
elevation with the hydraulic jump at the downwind base of the hill
(not seen in the photograph). This hydraulic jump was not present in
the half-depth tow, and the plume remained essentially at the water
surface. The presence of the hydraulic jump evidently affects the
flow structure on the lee side of the hill and, therefore, the
concentration distributions on the lee side.
260
-------
LJFDIF. C I U. G) POLTIIIL LIFT STUDT. COu 1, PORT NO.;CO^ 2 1C 7, CC'.MC. rf,3!
01-DEC-S3 .017. .014. .015, .015, fiN'D .018, RESPECT! V£LT.
100
Q
UJ
O
oc
111
2
CD
3
M
UL
<0
T-
I
»
I
U
1 10
CHI, TOW 14. HALF SUBMERGED
100
Figure A-17. Comparison of concentrations. Half submerged versus
fully submerged. F = 0.5, HS/H =0.8.
261
-------
Figure A-18.
Side views of plume released at 0.8 H. F = 0.5; upper,
fully submerged; lower, half submerged.
262
-------
Figure A-19.
Top views of plume released at 0.8 H. F=0.5; upper,
fully submerged; lower, half submerged.
263
-------
4. CONCLUSIONS
A series of tows was conducted in a stably stratified, salt-water
towing tank wherein the density gradient was linear and the
dividing-streamline height was half the hill height. Effluent was
released at 0.1 h, 0.2 h and 0.3 h above the dividing-streamline
height, i.e., at elevations above the base of the hill of 0.6 h, 0.7 h
and 0.8 h. Pairs of tows were made such that, in one tow, the hill
war fully immersed ir the water and the towing speed was adjusted to
provide a "natural" dividing-streamline surface (i.e., not flat) of
half the hill height (far upstream). In the second tow of the pair,
the hill was raised out of the water to the point where only the top
half of the hill was immersed; this, in effect, forced a flat
dividing-streamline surface (the water surface itself). In both tows
of each pair, all other conditions were maintained identical, e.g.,
towing speed, effluent release height and flow rate, and density
gradient. Concentration distributions were measured on the hill
surface and in the absence of the hill. The distributions from each
pair of tows were compared to ascertain the effects of an assumed flat
dividing-streamline surface.
Conclusions that may be drawn from this study are as follows:
1. Repeatability of concentration patterns from one tow to the
next, under ostensibly identical conditions, is regarded as
excellent, generally within +10%. The few points outside
this range (but generally well within a factor of two) were
found to be on the edges of the plume where the
concentrations were quite small and the turbulent
fluctuations naturally gave rise to this kind of variability.
2. At a stack height of 0.6 h (0.1 h above Hc), the maximum
concentration occurred on the windward side of the hill; the
plume was spread broadly to cover essentially the entire top
half of the hill, but its thickness in the direction normal
to the hill surface was small. Isoconcentration contours on
the windward side were roughly circular. Comparison of the
concentration patterns from the half- and full-immersion
tows showed remarkable similarities over the top half of the
hill. The locations of the maxima were essentially the
same, although the value of the maxima was 10 to 20% larger
for the half-immersion case. Regions of poor agreement
included a thin band just above the half-hill height (where
concentrations differed, in the worst case, by a factor of
about four) and, of course, everywhere below the half-hill
height.
264
-------
3. At a stack height of 0.7 h (0.2 h above Hc), the maximum
surface concentration occurred just downwind of the hill
center, with a value nearly the. same as with the stack.
height of 0.6 h. The plume was again spread broadly to
cover most of the top half of the hill. Isoconcentration
contours were again roughly circular on the windward side,
but somewhat elongated in the flow direction. Comparison of
the concentration patterns from the full- and half-depth
tows again showed remarkable similarities and quite good
agreement in the areas of higher concentration. Regions of
poor agreement covered a broader band just above the
half-hill height and, in that region, concentration values
differed by as much as a factor of 8.
4 At a stack height of 0.8 h (0.3 h above Hc), the maximum
concentration occurred well to the lee side of the hill and
its value was substantially smaller than the maxima observed
in the lower stack cases. Also in comparison with the lower
stack cases, this plume was quite narrow, and the area of
coverage on the hill surface was quite small. Comparison of
the concentration patterns from the full- and half-depth
tows showed that, in areas of large concentration, values of
concentration differed by roughly factors of two, In areas
of low concentration, values of concentration differed by up
to a factor of 10.
5. Observations of the flow structure revealed that a hydraulic
jump occurred near the downwind base of the fully immersed
hill, but was absent in the case of the half-immersed hill.
This hydraulic jump appeared to have a rather minor effect
on the streamline pattern near the surface on the lee side
of the hill, but, as pointed out by Snyder and Hunt (1984),
surface concentrations are extremely sensitive to the ratio
of the plume width to the normal distance of the plume
centerline from the hill surface when this ratio is near
unity. Hence, even a seemingly minor change in the
streamline patterns can have a major effect on the
concentration patterns.
6. Maximum concentrations on the surface under the "worst"
conditions were equal to or somewhat greater than those
observed at the plume centerline in the absence of the hill.
7. These results suggest that the assumption of a flat
dividing-streamline surface in a mathematical model is a
reasonable approximation to make, at least with regard to
predicting the locations and values of maximum
concentrations (and areas of coverage) on the windward side
of the hill. When the stack heights are relatively close to
the dividing-streamline height, the lee-side concentrations
would also be predicted reasonably well.
265
-------
Future work should include an improved concentration-
measurement technique to resolve the question of whether or
not hill-surface concentrations can exceed those at the
plume centerline in the absence of the hill (and by how
much).
266
-------
REFERENCES
Bass, A., Strimaitis, D.6. and Egan, B.A., 1981: Potential Flow Model
for Gaussian Plume Interaction with Simple Terrain Features, Rpt.
to Envir. Prot. Agcy. under Contract No. 68-02-2759, Res. Tri.
Pk., NC, p. 201.
Brighton, P.W.M., 1978: Strongly Stratified Flow Past Three-Dimensional
Obstacles, Quart. J. Roy. Meteorol. Soc., v. 104, p. 289-307.
Drazin, P.G., 1961: On the Steady Flow of a Fluid of Variable Density
Past an Obstacle, Tellus, v. 13, no. 2, p. 239-51.
Hunt, J.C.R., Puttock, J.S. and Snyder, W.H., 1979: Turbulent
Diffusion from a Point Source in Stratified and Neutral Flows
around a Three-Dimensional Hill: Part I: Diffusion Equation
Analysis, Atmos. Envir., v. 13, p. 1227-39.
Hunt, J.C.R., Richards, K.J. and Brighton, P.W.M., 1984: Stratified
Shear Flow over Low Hills: II. Stratification Effects in the
Outer Flow Region, To be submitted to Quart. J. Roy. Meteorol.
Soc.
Hunt, J.C.R. and Snyder, W.H., 1980: Experiments on Stably and
Neutrally Stratified Flow over a Model Three-Dimensional Hill, J.
Fluid Mech., v. 96, pt. 4, p. 671-704.
Hunt, J.C.R., Snyder, W.H., and Lawson, R.E. Jr., 1978: Flow Structure
and Turbulent Diffusion around a Three-Dimensional Hill: Fluid
Modeling Study on Effects of Stratification; Part I: Flow
Structure, Rpt. No. EPA-600/4-78-041, Envir. Prot. Agcy., Res.
Tri. Pk., NC.
Lavery, T.F., Bass, A., Strimaitis, D.6., Venkatram, A., Greene, B.R.,
Drivas, P.J. and Egan, B.A., 1982: EPA Complex Terrain Modeling
Program: First Milestone Report - 1981, Rpt. No. EPA-600/3-82-036,
Envir. Prot. Agcy., Res. Tri. Pk., NC, p. 304.
Lavery, T.F., Strimaitis, D.G., Venkatram, A., Greene, B.R.,
DiCristofaro, D.C. and Egan, B.A., 1983: EPA Complex Terrain
Model Development: Third Milestone Report-1983, Rpt. No.
EPA-600/3-83-101, Envir. Prot. Agcy. Rpt., Res. Tri. Pk., NC, p.
271.
267
-------
Riley, J.J., Liu, H.T. and Geller, E.W., 1976: A Numerical and
Experimental Study of Stably Stratified Flow Around Complex
Terrain, Rpt. No. EPA-600/4-76-021, Envir. Prot. Agcy., Res. Tri.
Pk., NC, p. 41.
Sheppard, P.A., 1956: Airflow over Mountains, Quart. J. Roy. Meteorol.
Soc., v. 82, p. 528-9.
Snyder, W.H., Britter, R.E. and Hunt, J.C.R., 1980: A Fluid Modeling
Study of the Flow Structure and Plume Impingement on a
Three-Dimensional Hill in Stably Stratified Flow, Proc. Fifth
Int. Conf. on Wind Engr. (J.E. Cermak, ed.), v. 1, p, 319-29,
Pergamon Press, NY, NY.
Snyder, W.H. and Hunt, J.C.R., 1984: Turbulent Diffusion from a Point
Source in Stratified and Neutral Flows around a Three-Dimensional
Hill; Part II: Laboratory Measurements of Surface Concentration,
Atmos. Envir. (to appear).
Strimaitis, D.G., Venkatram, A., Greene, B.R., Hanna, S., Heisler, S.,
Lavery, T.F., Bass, A., and Egan, B.A., 1982: EPA Complex Terrain
Model Development: Second Milestone Report - 1982, Rpt. No.
EPA-600/3-83-015, U.S. Envir. Prot. Agcy. Rpt., Res. Tri. Pk., NC.
Townsend, A.A., 1956: The Structure of Turbulent Shear Flow, Cambridge
Univ. Press, Cambridge, England, p. 315.
268
-------
APPENDIX B
DISPERSION FROM A SOURCE UPWIND OF A
THREE-DIMENSIONAL HILL OF MODERATE SLOPE
269
-------
DISPERSION FROM A SOURCE UPWIND OF A
THREE-DIMENSIONAL HILL OF MODERATE SLOPE
Roger S. Thompson
and
William H. Snyder*
Meteorology and Assessment Division
Atmospheric Sciences Research Laboratory
Environmental Protection Agency
Research Triangle Park, NC 27711
May 1984
*0n assignment from the National Oceanic and Atmospheric
Administration, U.S. Department of Commerce.
270
-------
1. INTRODUCTION
The series of experiments performed in 1980 at Cinder Cone Butte
(CCB), Idaho, by Environmental Research and Technology, Inc.
(Strimaitis et al., 1983 and Lavery et al., 1983) provide a data base
for development and evaluation of dispersion models for complex
terrain situations. As a follow-up and companion effort, fluid
modeling studies have been performed at the EPA Fluid Modeling
Facility using both a stratified towing tank and a meteorological wind
tunnel. The work described in this report was inspired by an earlier
case study of a one-hour period of the field study for which the
atmosphere was neutral. That wind-tunnel study (see Appendix A of
Lavery et al., 1983) provided good agreement between laboratory and
field measurements of concentrations. However, the maximum observed
concentration on the hill surface was found to be only about 10%
greater than that observed from the same source over flat terrain;
that is, the influence of CCB in increasing the maximum observed
ground-level concentration (glc) was not dramatic. Thus, a set of
experiments to determine the influence of the hill on the maximum glc
and to locate the source positions where this influence is the
greatest was a natural continuation of the fluid modeling effort. The
data base for mathematical modeling evaluation is extended by these
measurements. All measurements were made under simulated atmospheric
conditions.
The nearly-axisymmetric CCB shape was replaced with a truly
axisymmetric hill having a simple mathematical formula. The average
silhouette of CCB was best fit by the formula to maintain a close
relationship with the previous efforts: the maximum slope of CCB was
matched in the model.
Since, for upwind sources, the maximum ground-level concentration
can be expected to occur from sources on the centerline, this study
was restricted to only those source positions. The maximum
concentrations from these sources are expected to occur along the
centerline through the source and hill center; thus, measurements were
limited to centerline sampling port locations.
The measure of the influence of the hill on the maximum glc was
computed as a "terrain amplification factor" A. This factor is
defined as the ratio of the maximum glc observed with the hill present
to the maximum observed in the absence of the hill. The location of
the maxima are not considered in this evaluation.
271
-------
2. EXPERIMENTAL DETAILS
The model selected for this study is an idealization of Cinder
Cone Butte (CCB), an isolated three-dimensional hill about 100 meters
high. As described in the First Milestone Report (Lavery et al.,
1982), CCB has a double peak but is roughly axisymmetric. Typical
maximum slopes are about 25 degrees.
Profiles of the shape of CCB along radial lines for azimuths of
0, 10, 20, .... and 350 degrees were graphed. Several functional
forms were considered as approximations to the average silhouette of
CCB. The function that was found to best fit the hill shape is
h(r) = {(H+«)/(l+(r/L)4) - & 0 < r < R
0 r > R
where: r is the radius,
h is the local height,
H is the height at the center,
6 is a vertical offset for a finite model,
L is the radius for which h = H/2 if 6 =0,
and R is the radius of the modeled area.
The vertical offset was included in the formula to avoid a
discontinuity at the edge of the model. For CCB, H = 100 m,
L = 250 m, R = 500 m, and 4 = 6.25 m were chosen. This gives a
maximum slope angle of 24.4 degrees. Comparison of this function with
CCB silhouettes at intervals of 90 degrees in the azimuth are shown in
Figure B-l.
Another appealing reason to use this hill shape is that Hunt et
al. (1978) previously used the same functional form (however, with
H/L = 1 which gives a maximum slope angle of 45 degrees) in
experiments that established a large data base, which is available for
comparison. Also, the Fourier transform of this function should be
obtainable for use in mathematical models.
A scale model of the idealized CCB was vacuum molded of acrylic
plastic over a wood form. The height of the model was 15.5 cm which
corresponds to a scale of 1:640 based on the 100 m height assumed for
CCB. The radius of the model was 0.78 m. The vertical offset in the
shape function prevented a discontinuity in the surface at the edge of
the model. If the model had followed the formula exactly there would
have been a discontinuity in the slope at the model edge of 2.7
degrees. But the edge was faired during construction.
272
-------
IOC
?no 300
R, METERS
U'lO
Figure B-l. Silhouettes of Cinder Cone Butte and the profile selected
for this axisymmetric hill study.
273
-------
The model was covered with gravel of the same size used in
roughening the test section floor for boundary layer generation. The
same boundary layer generating devices used in the study described in
the Third Milestone Report (Lavery et al., 1983) were again used for
this study. The model was located with its center at a distance of
13.8 m downstream of the entrance to the test section.
Mean velocity and turbulence intensity measurements were made
with hot-film anemometers (TSI model 1241 probes and model 1053B
anemometers). A bent-over brass tube (o.d. =0.32 cm) was used to
emit a non-buoyant tracer gas (methane, 99% pure). The emission rate
was 2950 cm^/min which gave an exit velocity of about two times the
local speed. Samples were drawn through ports located on the hill
surface or through sampling rakes, positioned with the instrument
carriage, for locations upwind or downwind of the hill or in the
absence of the hill. Flame ionization detectors (Beckman model 400
Hydrocarbon Analyzers) were used to measure the sampled
concentration. Two-minute averages were obtained of values sampled at
rates of one sample per second for the flame ionization detectors and
500 samples per second for the anemometers with the laboratory's
minicomputer.
When necessary to consider actual concentrations rather than
terrain amplification factors or other ratios of concentrations, a
non-dimensional concentration x was calculated according to
X = CU(H)H2/Q
Where C = the measured concentration,
U(H) = the approach wind speed at z = H,
H = the hill height, and
Q = the emission rate of tracer.
For visualization, plumes were generated with cotton tipped glass
tubes of liquid titanium tetrachloride (marketed as "smoke sticks" by
E. Vernon Hill, Inc.). Photographs of these plumes were taken with a
Graphlex 4X5 camera fitted with a Polaroid back and type 55 P/N film.
Two 500 w photolamps were directed at the hill from elevated positions
downwind of the hill on either side of the test section. A
"time-averaged" photograph was obtained by making eight one-second
exposures with an interval of a few seconds between the exposures.
274
-------
3. MEASUREMENTS AND RESULTS
Vertical and lateral profiles of mean velocity and turbulence
intensity were obtained in the wind tunnel test section before
installing the model. The vertical profiles were measured at
distances of 10 and 14 m from the beginning of the floor roughness
(11.6 hill heights upstream and 14.2 hill heights downstream of the
position selected for the hill center). These profiles define the
approach wind field and demonstrate the uniformity of the wind tunnel
flow over the portion of the test section used for this study (Figures
B-2 and B-3). Lateral traverses were made at the 10 m distance at
heights of one-half and one times the hill height. Along these
traverses, the mean velocity and turbulence intensities varied by less
than 4%.
Fitting a power-law profile to the approach flow mean velocity
gives an exponent of 0.166, that is, u(z)« z°-166. A logarithmic
fit gives z0 = 0.2 mm or 0.0013 times the hill height and
u*/k =0.51 m/s.
As required for computation of A, ground-level concentrations
were measured for sources in the absence of the hill. Stack heights
of 0, 40, 80, 120, 160, and 200 mm were selected to be integer
multiples of 1/4 of the anticipated hill height of 160 mm; however,
the true height of the model came out to be 155 mm. These same stack
heights were used with the model in place and are referred to as stack
heights of 1/4, 1/2, 3/4, 1, and 1 1/4 hill heights. The distances
from the hill center to the source were also scaled by the hill
height.
The reference (no hill) data base consists of longitudinal
profiles of glc for the various stack heights (Figure B-4) normalized
(Berlyand, 1975) using the maximum glc (Xmax^ an<* tlie distance to
that maximum (xmax') as similarity parameters (Figure B-5). The
dependence of these parameters on the source height Hs was
determined to be
X aH1'15.
max s
The maximum concentrations from these profiles were used to calculate
the terrain amplification factors.
275
-------
7 --
6 --
A 10 M FROM STflRT OF ROUGHNESS
D l"i M
5 --
T"
S
6
3 -•
2 --
.25
.5 75
U IZ) /U (H)
A D
AD
A
1.2E
Figure B-2. Vertical profiles of the mean air speed in the test
section of the wind tunnel. H = Hill height.
276
-------
7 --
5 --
BE
so
O
A UVU, 10 M
A i.,1 Vu, i 4 y
D V'/V, iO K
• V'/v, 114 f
O uVU, 10 M
* W'/U, 114 M
A
Figure B-3. Vertical profiles of local turbulence intensities in the
test section of the wind tunnel.
277
-------
c
t\i
LJ
A
—T£i i r
2i
A
O
i
j.
A
o
Figure B-4. Longitudinal profiles of ground-level concentration from
point sources in the absence of the model.
278
-------
o
A
O
O
O
L1
O
a
A
A
A
A ,
A
A
c
A
Figure B-5. Longitudinal profiles of ground-level concentration from
point sources in the absence of the model plotted in
similarity coordinates.
279
-------
With the hill in place, glc's were measured on the hill via
surface ports on the hill and both upwind and downwind of the hill
with sampling rakes. Sources were located at distances of 4 to 16
hill heights upwind of the center of the model and at heights of 1/4
to 1 1/4 times the hill height. Complete longitudinal profiles of the
glc were obtained for selected source positions; but for most
positions, only sufficient measurements were made to find the maximum
glc.
The presence of the hill may be considered to influence the
transport and dispersion of the plume to increase the maximum glc in
three ways. For low sources, at moderate distances from the hill, the
reduction in mean wind speed allows the plume to reach the ground
surface closer to the source, producing higher concentrations than in
the absence of the hill. Plumes from higher sources can be thought of
as being intercepted by the hill. That is, the hill penetrates the
plume where the concentrations are greater than those that would occur
at ground level farther downstream over flat terrain. For yet higher
sources, the streamline convergence over the hilltop and the
corresponding downward flow in the lee of the hill brings the plume to
the ground more rapidly than over flat terrain. The maximum
concentration for these three regimes occurs upwind of the hilltop,
near or on the hilltop, or downwind of the hilltop, respectively. The
distance from the source to the maximum glc was always observed to be
less in the presence of the hill. No amplification factors less than
1.0 were observed in this study; that is, the hill's presence never
influenced the dispersion from these upwind sources to produce a
maximum glc that was less than would occur over flat terrain.
Hunt et al. (1979) presented a theoretical argument for the
presence of a "window" of upwind source locations that produce a
terrain amplification factor of 2 or greater for three-dimensional
terrain objects. They show that for sources within this window, the
maximum concentrations occur on or near the hill.
A series of longitudinal profiles of glc for all source heights
at a. distance of six hill heights upwind of the hill center exhibit
these types of influence (Figure B-6). The Hs = 0.25 H source has
its maximum glc occur near the upwind base of the hill; a low plume
being slowed by the presence of the hill. The maximum glc for the
Hs = 0.75 H stack occurs on the hill top; the hill intersecting the
plume. The Hs = 0.5 H stack produces a large region of nearly
constant concentration over the upwind face and crest of the hill; a
combination of the first two types of influence. Two nearly equal
maxima occur for the Hs = 1.0 H stack, one on the hill top and one
downwind of the hill. The maximum glc for the Hs = 1.25 H stack
clearly occurs downwind of the hill. Figure B-7 gives photographs of
the Hs = 0.5 H and Hs = 1.25 H plumes to illustrate their behavior.
All of the above concentration profiles exhibit a rapid reduction
in concentration at approximately 1 1/2 to 2 hill heights downwind of
the hill center. The flow separated from the hill surface at this
point and concentrations just downstream were determined by the
concentrations in the recirculating flow.
280
-------
A-
A
A.
A
A
X/L'
r
^ _ C"
Figure B-6. Longitudinal profiles of glc for sources located six hill
heights upwind of the hill center. Hill profile shown
below; note exaggerated scale.
281
-------
a) H /H = 0.50
b) HS/H = 1.25
Figure B-7. Visualization of plumes from sources at x/H = -5.
282
-------
Terrain amplification factors for all source positions are
presented in Figure B-8. Lines of constant A (= 1.4 and 2.0) are
shown as dashed lines. The solid lines divide the area into three
regions, where the source produces its maximum glc upwind of the
hilltop, on the hilltop (between the peak and the separation point),
and downwind of the hill.
Terrain amplification factors were computed for a hill of similar
functional form but much greater slope (45°) using data from Snyder
and Britter (1979) who report on a continuation of the work by Hunt
et al. (1978). That study used a boundary layer of approximately the
same depth but it was generated with a rougher surface (larger size of
gravel coating) and a smooth hill that was one and a half times the
height of the hill used in the present effort. The rougher surface
produced turbulence intensities a factor of two larger than in the
present study. Surface concentration profiles were obtained for two
source positions: release heights of one-half and one times the hill
height at a distance of 3.62 hill heights upwind of the hill center.
The terrain amplification factors for these source locations were
computed to be A = 4.00 and 1.27 for the stacks one-half and one times
the hill height, respectively. Values for the present study were
calculated from releases at a slightly greater upwind distance (4 hill
heights). They were 3.63 and 1.27 for the same stack heights,
respectively.
The earlier work with the scale model of Cinder Cone Butte (see
Appendix A, Lavery et al., 1983) also provides a point for
comparison. The source was located 10 hill heights upwind, if we
consider the hill height to be 100 m. However, it was not directly
upwind of the hill center but on a line 3° from the mean wind
direction. The source height was 0.42 H. The maximum concentration
occurred at a position toward the side of the hill about halfway up.
The terrain amplification factor was computed to be 1.08. This falls,
between the values for the present study of 1.0 and 1.55 for stacks at
the same upwind distance but heights of 0.25 H and 0.50 H,
respectively.
283
-------
\ •
l
284
-------
4. SUMMARY AND CONCLUSIONS
A wind-tunnel study was performed to investigate the influence of
an idealized three-dimensional hill on ground-level concentrations
from upwind sources in a neutral atmospheric boundary layer and to
locate the source locations where this influence is the greatest.
The presence of the hill was found to influence the transport and
dispersion of the plume to increase the maximum glc in three ways.
For low sources, at moderate distances from the hill, the reduction in
mean wind speed allows the plume to reach the ground surface closer to
the source, producing higher concentrations than in the absence of the
hill. Plumes from higher sources can be thought of as being
intercepted by the hill. That is, the hill penetrates the plume where
the concentrations are greater than those that would occur at ground
level farther downstream over flat terrain. For yet higher sources,
the streamline convergence over the hilltop and the corresponding
downward flow in the lee of the hill brings the plume to the ground
more rapidly than over flat terrain. The maximum concentration for
these three regimes occurs upwind of the hilltop, near or on the
hilltop, or downwind of the hilltop, respectively. Terrain
amplification factors ranged from near 1.0 to 3.63. The region of
source locations that produced an amplification factor of 1.4 or more
extended to an upwind distance of 14 hill heights.
285
-------
REFERENCES
Berlyand, M.G., 1975: Contemporary Problems of Atmospheric Diffusion
and Air Pollution, Hydromet. Press, Leningrad, USSR, 448p.
Hunt, J.C.R., Puttock, J.S. and Snyder, W.H., 1979: Turbulent Diffusion
from a Point Source in Stratified and Neutral Flows Around a
Three Dimensional Hill: Part I: Diffusion Equation Analysis,
Atmos. Envir., v. 13, p. 1227-39.
Hunt, J.C.R., Snyder, W.H, and Lawson, R.E., 1978: Flow Structure and
Turbulent Diffusion around a Three-Dimensional Hill: Fluid
Modeling Study on Effects of Stratification: Part I: Flow
Structure, Rpt. No. EPA 600/4-78-041, Envir. Prot. Agcy., Res.
Tri. Pk., NC, 96p.
Lavery, T.F., Bass, A., Strimaitis, D.G., Venkatram, A., Greene, B.R.,
Drivas, P.J. and Egan, B.A., 1982: EPA Complex Terrain Model
Development: First Milestone Report - 1981, Rpt. No. EPA
600/3-82-036, Envir. Prot. Agcy., Res. Tri. Pk., NC, 304p.
Lavery, T.F., Strimaitis, D.G., Venkatram, A., Greene, B.R.,
DiCristofaro, D.C. and Egan, B.A., 1983: EPA Complex Terrain
Model Development: Third Milestone Report - 1983, Rpt. No. EPA
600/3-83-101, Envir. Prot. Agcy., Res. Tri. Pk., NC, 271p.
Snyder, W.H. and Britter, R.E., 1979: Aspect Ratio Study: Unpublished
In-house Data Report, Envir. Prot. Agcy., Res, Tri. Pk., NC.
Strimaitis, D.G., Venkatram, A., Greene, B.R., Hanna, S.R.,
Heisler, S., Lavery, T.F., Bass, A. and Egan, B.A., 1983: EPA
Complex Terrain Model Development: Second Milestone Report -
1982, EPA 600/3-85-015. Envir. Prot. Agcy., Res. Tri. Pk., NC
375p.
286
-------
APPENDIX C
ALPHA-1 OBSERVATIONS OF THE
TRACY OIL-FOG PLUME*
*This appendix is a condensation of a report by Uthe and Morley (1984)
287
-------
C.I Introduction
During the November experiment the ALPHA-1 airborne lidar mapped
aerosol plume behavior in the vicinity of the Tracy Generating
Station. The ALPHA-1 was based at the Reno airport and was operated
by SRI International under the sponsorship of the Electric Power
Research Institute.
An oil-fog tracer aerosol was generated and emitted froit. the
Tracy Station with a plume sufficiently dense for photographic
purposes. The ALPHA-1 mapped the plume behavior by making cross-plume
traverses above the plume at various downwind distances from the Tracy
Station. Typically, the plume traverses were made at downwind
increments of 0.5 nmi from the Tracy Station to 6 nmi downwind. The
aircraft typically operated at an altitude of 10,000 ft MSL, although
lower altitudes were flown during cloudy conditions. The mountainous
terrain and nighttime conditions prevented data collection during
times when cloud bases were below 7500 ft. The experimental data
collection typically began at about midnight and extended to
0900 PST. The ALPHA-1 normally made two flights, each between three
and four hours in duration.
The ALPHA-1 operation consisted of transmitting 1.06 ym
wavelength laser pulses at a repetition rate of 5/second (horizontal
resolution of about 12 m depending on aircraft ground speed). See
Uthe et al. (1980a and 1980b) for a description of the ALPHA-1
system. Received energy was detected, logarithmically amplified, and
digitized at 20 ns (3 meter) intervals. Lidar signatures were stored
on nine-track 1600 bpi magnetic tape and were also used to generate
real-time facsimile pictorial displays of plume, atmosphere and
terrain structure. Aircraft location was determined both by a LORAN-C
and RNAV system. The RNAV was programmed to provide distance and
bearing to the Tracy Station. Aircraft location and time information
were recorded on magnetic tape for each laser firing. The RNAV data
also were plotted on the facsimile recorder to provide the information
needed to direct the lidar data collection. The ALPHA-1 logarithmic
amplifier saturates for signals greater than about four orders of
magnitude (40 dB) . The ground returns are normally saturated,,
Because of the dense aerosol plume, plume returns near the source were
near saturation of the receiver amplifiers. Therefore, near-surface
plume returns were difficult to separate from ground returns. To
solve this problem, the transmitted energy was decreased by a factor
of 2 (3 dB), providing more dynamic range between plume and ground
returns. However, the reduced transmitted energy also decreased
sensitivity to low-density atmospheric aerosol features. A log of
data collection times is presented in Table C-l.
The ALPHA-1 data records are stored on nine-track 1600 bpi
magnetic tape in binary forms. The data consists of 1000 8-bit lidar
backscatter values for each laser firing. Also recorded in each block
of data are the time and location. The magnetic tape records are
maintained at SRI for future analysis.
288
-------
TABLE C-l
LOG OF ALPHA-1 DATA COLLECTION
TRACY GENERATING STATION, NOVEMBER 1983
Date
(Nov 83)
5
6
7
8
9
12
12
14
14
15
15
15
16
17/18
18
19
19
Time
(PST)
— —
1125-1140
0155-0420
0015-0325
0240-0530
0035-0340
0500-0720
0040-0335
0435-0800
0025-0340
0435-0735
2135-2330
0235-0240
2225-0130
0320-0345
0025-0130
Flight Time
(hr)
1.2
1.0
3.5
3.7
3.4
3.6
2.8
3.4
3.8
3.7
3.4
2.5
0.5
3.5
0.8
1.5
1.8
Remarks
ferry
test flight
2nd flight
2nd flight
2nd flight
clouds and turbulence
2nd flight/clouds
rain and snow
ferry
Total
44.1
289
-------
The ALPHA-1 magnetic tapes were used to generate facsimile
displays that depict plume, atmosphere, and terrain structure. The
displays also provide time and aircraft location. Tables presenting
information needed to interpret the lidar pictures are given in
Section C.2 and the lidar pictures are reproduced in Section C.3. The
reproduced pictorial records are reduced in size and the originals
provide a better indication of plume behavior. These original
facsimile records are available on a loan basis from SRI.
C.2 Date, Time, Altitude, and Location Data
This section provides a listing (Table C-2) of date and time
information needed to interpret the ALPHA-1 facsimile data plots
presented in Section C.3.
The listing provides the page number of Appendix B of the
original report by Uthe and Morley (1984) that presents facsimile
plots for each of the 87 data tapes collected during the Tracy
program. Tape number is identified on the lower left hand side of the
facsimile plots. Picture number refers to the data presented within
vertical solid lines and each picture presents an individual ALPHA-1
cross section. The ALPHA-1 was turned off during times between the
individual cross sections. Time (PST) is given for the first time
mark of each picture. Time marks are identified in the key to reading
the lidar cross sections presented in Figure C-l. In Section C.3 only
the facsimile plots representative of the three case studies described
in Section 6 of this Milestone Report are given.
Aircraft location is determined by an area navigation (RNAV)
system that provides information on bearing and distance to a selected
waypoint position. The waypoint was established at the Tracy Station,
and distance and bearing to the station are plotted on the facsimile
data displays. Figure C-2 presents a map of the Tracy area with
ALPHA-1 coordinates of distance (nmi) and bearing to the Tracy Station
superimposed. LORAN-C data on longitude and latitude were also
recorded but not plotted on the displays. Figure C-l is an example of
ALPHA-1 backscatter signatures processed as a facsimile record for
Tape 45 (page 8-23). Terrain, plume, and atmospheric features are
identified. Location of time marks, altitude scale, and aircraft
location plots are given as a key for interpreting the cross sections
presented in Section C.3. The altitude scale is 600 m/division for
all data.
290
-------
TABLE C-2
DATE AND TIME DATA FOR ALPHA-1 DATA COLLECTION
AT THE TRACY GENERATING STATION, NOVEMBER 1983
Gray Scale Depiction
Gray Scale Depiction
Page * Tape
B-2
1
B-2 2
B-3 3
B-3 4
B-4 5
Picture
1
1
2
3
4
5
6
-7
8
1
2
3
4
5
6
7
8
1
2
3
4
5
6
7
1
2
3
4
5
6
7
L/clUC/
Time (PST)
6 Nov 1983(1)
1125
7 Nov 1983(2)
0156
0159
0202
0204
0207
0212
0216
0219
0223
0225
0227
0230
0233
0236
0242
0245
0253
0256
0259
0303
0307
0310
0314
0317
0321
0325
0329
0332
0338
0342
Page
B-4
B-5
B-5
B-6
B-6
Tape Picture
6 1
2
3
4
5
6
7
8
7 1
2
3
4
5
8 1
2
3
4
5
6
7
8
9
9 1
2
3
4
5
6
7
10 1
2
3
Time (PST)
0346
0349
0354
0356
0359
0401
0404
0406
0409
0411
0414
0416
0419
8 Nov 1983
0014
0017
0020
0023
0026
0029
0031
0034
0037
0042
00^5
0048
0054
0059
0102
0105
0109
0113
0116
(1) Location from Reno Vortac (test flight).
(2) All following locations from Tracy Power Plant.
* Page number of the original report by Uthe and Morley (1984)
291
-------
TABLE C-2 (Continued)
Gray Scale Depiction
Picture
B-6 10
Time (PST)
B-7 11
B-7 12
B-8 13
B-8 14
B-9 15
B-9 16
4
5
6
7
1
2
3
4
5
1
2
3
4
5
6
1
2
3
4
5
6
7
1
2
3
4
5
6
7
8 Nov 1983
0119
0123
0128
0131
0138
0141
0144
0149
0154
0200
0203
0206
0216
0219
0222
0227
0331
0235
0243
0246
0251
0253
0257
0300
0303
0306
0310
0312
0316
1
2
3
4
5
6
0326
9 Nov 1983
0244
0246
0250
0252
0255
0258
Page Tape
B-9 16
B-9 17
B-10 18
B-10 19
B-ll 20
B-ll 21
Picture
7
8
9
10
1
2
3
4
5
6
7
8
1
2
3
4
5
6
1
8
1
2
3
4
5
6
7
8
1
2
3
4
5
6
1
2
3
4
uaue/
Time (PST)
9 Nov 1983
0302
0304
0307
0309
0312
0314
0317
0319
0322
0325
0328
0330
0336
0342
0345
0348
0352
0354
0358
0401
0406
0409
0412
0414
0417
0420
0422
0425
0434
0438
0442
0446
0449
0453
0459
0506
0509
0513
292
-------
TABLE C-2 (Continued)
Gray Scale Depiction
Page Tape
B-ll 21
B-12 22
B-12 23
B-13 24
8-13 25
B-14 26
B-14 27
B-15 28
Picture
5
1
2
1
2
3
4
C
_>
1
2
3
4
5
6
1
2
3
4
5
6
1
2
3
4
5
6
1
2
3
4
5
1
2
3
4
uate/
Time (PST)
9 Nov 1983
0516
0522
0528
12 Nov 1983
0035
0038
0042
0046
0049
0057
0102
0106
0110
0113
0116
0120
0124
0128
0131
0135
0138
0143
0146
0150
0153
0157
0203
0207
0213
0216
0220
0223
0228
0231
0235
0242
: :
Page Tape Picture
Brl5 28 5
B-15 29 1
2
3
4
5
B-16 30 1
2
3
4
5
B-16 31 1
2
3
4
B-17 32 1
2
3
4
5
B-17 33 1
2
3
4
5
6
B-18 34 1
2
3
4
5
B-18 35 1
2
3
4
5
Date/
Time (PST)
12 Nov 1983
0247
0253
0258
0302
0307
0311
0317
0322
0329
0334
0338
0503
0506
0512
0519
0525
0532
0536
0539
0543
0549
0554
0558
0602
0605
0609
0613
0617
0621
0625
0630
0636
0641
0644
0647
0650
293
-------
TABLE C-2 (Continued)
Gray Scale Depiction
Page Tape
B-18 35
B-19 36
B-19 37
B-20 38
B-20 39
B-21 40
B-21 41
B-22 42
Picture
6
1
2
3
4
1
2
3
4
5
6
7
1
2
3
4
5
1
2
3
4
5
1
2
3
4
5
1
2
3
4
5
1
2
3
Date/
Time (PST)
12 Nov 1983
0655
Gray Scale Depiction
Picture
1
2
3
4
1
2
3
4
5
6
7
1
2
3
4
5
1
2
3
4
5
1
2
3
4
5
1
2
3
4
5
1
2
3
0702
0709
0712
0718
14 Nov 1983
0043
0046
0049
0052
0055
0100
0105
0110
0114
0119
0125
0130
0134
0139
0142
0147
0151
0200
0204
0208
0213
0216
0222
0226
0231
0234
0240
0246
0250
0252
B-22 42
B-22 43
B-23 44
B-23 45
B-24 46
B-24 47
4
5
6
Date/
Time (PST)
Nov 1983
0256
0259
0303
B-25 48
I
2
3
4
5
6
7
8
1
2
3
4
5
6
1
2
3
4
5
6
1
2
3
4
5
1
2
3
4
5
6
1
2
3
0313
0317
0320
0323
0326
0330
0334
0336
0439
0441
0444
0467
0451
0454
0459
0502
0506
0509
0513
0516
0521
0524
0528
0539
0547
0551
0555
0558
0602
0605
0608
0624
0626
0629
294
-------
TABLE C-2 (Continued)
Gray Scale Depiction
Gray Scale Depiction
Page Tape
B-25 48
B-25 49
B-26 50
B-26 51
B-27 52
B-27 53
Picture
4
5
6
7
8
1
2
3
4
5
6
7
8
1
2
3
4
5
6
7
8
9
1
2
3
4
5
6
7
1
2
1
2
3
4
5
Date/
Time (PST)
14 Nov 1983
0632
0634
0637
0640
0642
0647
0650
0652
0655
0658
0701
0704
0707
0712
0715
0717
0720
0722
0724
0726
0728
0730
0734
0738
0741
0743
0746
0750
0753
0757
0800
15 Nov 1983
0026
0028
0031
0033
0036
Page Tape Picture
B-27 53 6
7
8
9
B-28 54 1
2
3
4
5
6
7
B-28 55 1
2
3
4
5
6
7
B-29 56 1
2
3
4
5
6
7
B-29 57 1
2
3
4
5
6
B-30 58 1
2
3
4
5
6
Date/
Time (PST)
15 Nov 1983
0039
0042
0044
0047
0051
0055
0057
0100
0105
0108
0111
0115
0119
0122
0124
0127
0130
0136
0139
0143
0146
0149
0152
0155
0158
0202
0205
0209
0212
0215
0219
0226
0229
0233
0237
0241
0245
295
-------
TABLE C-2 (Continued)
Gray
Page
B-30
B-31
B-31
B-32
B-32
B-33
B-33
Scale Depiction
Tape Picture
59 1
2
3
4
5
6
60 1
2
3
4
5
6
7
8
61 1
62 1
2
3
4
5
6
63 1
2
3
4
64 1
2
3
4
5
65 1
2
3
4
5
6
Date/
Time (PST)
15 Nov 1983
0249
0253
0256
0301
0305
0309
0313
0317
0320
0323
0326
0328
0332
0334
0339
0437
0440
0444
0448
0452
0455
0459
0505
0510
0515
0521
0525
0530
0534
0538
0542
0546
0552
0557
0600
0603
Gray Scale Depiction
Page Tape Picture
B-34 66 1
2
3
4
5
6
7
B-34 67 1
2
3
4
5
6
7
B-35 68 1
2
3
! 4
5
6
7
B-35 69 1
:f 2
3
4
5
B-36 70 1
i 2
3
4
5
6
7
8
B-36 71 1
2
3
Date/
Time (PST)
15 Nov 1983
0606
0609
0612
0615
0618
0621
0624
0633
0636
0638
0641
0644
0649
0651
0656
0701
0704
0707
0710
0713
0716
0719
0722
0725
0729
0732
2137
2139
2142
2146
2151
2154
2157
2201
2206
2209
2213
296
-------
TABLE C-2 (Continued)
Gray
Page
B-36
B-37
B-37
B-38
B-38
B-38
B-39
Scale Depiction
Tape Picture
71 4
5
6
72 1
2
3
4
5
6
7
73 1
2
3
4
5
6
7
74 1
2
75 1
2
76 1
2
3
4
5
6
7
77 1
2
3
4
5
6
7
Date/
Time (PST)
15 Nov 1983
2219
2222
2229
2233
2238
2240
2243
2245
2249
2252
2259
2302
2306
2309
2313
2317
2320
2325
2328
16 Nov 1983
0235
0237
17 Nov 1983
2228
2232
2236
2239
2243
2247
2250
2254
2257
2300
2304
2307
2311
2315
Gray Scale Depiction
Page Tape Picture
Date/
Time (PST)
B-39
78
B-40
79
B-41 81
B-41 82
1
2
3
4
5
1
2
3
4
17
Nov 1983
2319
2324
2331
2335
2340
2345
2349
2353
2357
B-40 .79
B-40 80
5
6
1
2
3
4
5
6
7
8
9
18 Nov 1983
0001
0004
0008
0011
0013
0016
0018
0024
0026
0029
0031
1
2
3
4
5
6
7
8
9
1
2
3
4
5
6
7
8
0034
0038
0040
0043
0046
0049
0054
0057
0059
0103
0105
0108
0110
0112
0115
0117
0123
297
-------
TABLE C-2 (Continued)
Gray Scale Depiction n t / Gray Scale Depiction n t- /
Page Tape Picture Time (PST) Page Tape Picture Time (PST)
18 Nov 1983
B-41 82 9 0126*
B-42 83 1 0131*
B-42 84 1 0322*
2 0327*
3 0330*
4 0336*
5 0339*
B-42 85
B-43 86
B-43 87
1
2
3
4
5
6
7
8
1
2
3
4
5
6
7
8
1
2
3
4
5
19 Nov 1983
0023
0025
0028
0030
0033
0036
0038
0042
0054
0059
0101
0105
0107
0110
0112
0114
0118
0122
0125*
—
0130
*
Cross section along axis of the plume.
298
-------
O
o
O
(H
O
(X
o
•p
X
I
u
00
299
-------
§
4->
I/I
X
O
•P
nj
•H
O
O
O
D.
rt
rt
-------
C.3 ALPHA-1 Facsimile Records for the Case Study Experiments
Discussed in Section 5
301
-------
302
-------
303
-------
304
-------
305
-------
306
-------
307
-------
308
-------
309
-------
310
-------
311
-------
Ill
312
-------
313
-------
314
-------
315
-------
316
-------
317
-------
318
-------
REFERENCES FOR APPENDIX C
Uthe, E.E., N.B. Nielsen, and W.L. Jimison 1980a: "Airborne Lidar
Plume and Haze Analyzer (ALPHA 1)," Bu 11. Amer. Meteoro1. Soc.,
61, 1035-1043.
Uthe, E.E., W.L. Jimison, and N.B. Nielsen 1980b: "Development of an
Airborne Lidar for Characterizing Particle Distribution in the
Atmosphere," Final Report, EPRI No. EA-1538, SRI International,
Menlo Park, CA.
Uthe, E.E. and B.M. Morley 1984. Alpha-1 Observations of Plume
Behavior for PMV&D Tracy Site, prepared for Electric Power
Research Institute, Palo Alto, California.
U S GOVERNMENT PRINTING OFFICE-559-013/20008
319
-------