t *j<"y 11 %i
e* 1983
EPA COMPLEX TERRAIN MODEL DEVELOPMENT
Description of a Computer Data Base
from Small Hill Impaction Study No. 1
Cinder Cone Butte, Idaho
ENVIRONMENTAL SCIENCES RESEARCH LABORATORY
OFFICE OF RESEARCH AND DEVELOPMENT
U.S. ENVIRONMENTAL PROTECTION AGENCY
RESEARCH TRIANGLE PARK, NORTH CAROLINA 27711
-------
EPA COMPLEX TERRAIN MODEL DEVELOPMENT
Description of a Computer Data Base
from Small Hill Impaction Study No. 1
Cinder Cone Butte, Idaho
by
Lawrence E. Truppi
and
George C. Holzworth
Meteorology and Assessment Division
Environmental Sciences Research Laboratory
U.S. Environmental Protection Agency
Research Triangle Park, North Carolina 27711
ENVIRONMENTAL SCIENCES RESEARCH LABORATORY
OFFICE OF RESEARCH AND DEVELOPMENT
U.S. ENVIRONMENTAL PROTECTION AGENCY
RESEARCH TRIANGLE PARK, NORTH CAROLINA 27711
-------
DISCLAIMER
This report has been reviewed by the Environmental Sciences Research
Laboratory of the U.S. Environmental Protection Agency and approved for
publication. Mention of trade names or commercial products does not
constitute endorsement or recommendation for use.
-------
ABSTRACT
As part of the U.S. Environmental Protection Agency's effort to develop
and demonstrate a reliable model of atmospheric dispersion for pollutant
emmissions in irregular mountainous terrain, the Complex Terrain Model
Development Program was initiated. The first phase, a comprehensive tracer
field study, was carried out on Cinder Cone Butte, Idaho, during the autumn
of 1980. Eighteen quantitative tracer experiments were conducted, each
lasting 8 hr at night or early morning. The main tracer gas was sulfur
hexafluoride; a second tracer, Freon 13B1 was used in ten of the eighteen
experiments. Averaged meteorological data were recorded from six towers
near and on the slopes on the hill. Data consisted of direct and derived
measures of temperature, wind, turbulence, solar and net radiation, and
nephelometer coefficent of scattering. Hourly wind profiles were obtained
from pilot balloon observations; tethersonde observations recorded profiles
of wind and temperature.
Tracer gas concentrations were detected by a network of approximately
100 samplers located on the slopes of the hill. The system used to collect
the data, the operational procedures used to run the system, and its
performance record are described. Tables of tracer gas release data have
been included to assist in any modeling effort. All meteorological and
i n
-------
tracer concentration data have been edited and recorded on magnetic tape and
are now available, upon request, at the National Computer Center, Research
Triangle Park, North Carolina, either as copies or by interactive computer
access.
-------
CONTENTS
Abstract i i i
Figures vii
Tables viii
List of Symbols and Abbreviations x
Acknowledgements xi i i
1. Introduction 1
1.1 EPA program 1
1.2 Objective 2
2. Field Study at Cinder Cone Butte 5
2.1 Geographic and meteorological settings 5
2.2 Experimental design 6
3. Tower Meteorological Data 11
3.1 Fixed meteorological network 11
3.1.1 Data acquisition system 21
3.1.2^ Quality assurance plan 23
3.1.3 Data editing 26
3.1.4 Periods of data collection 31
3.2. Meteorological Data Tape Files 32
3.2.1 Tape file index 33
3.2.2 Tape file records 37
3.2.3 Data record types 38
-------
CONTENTS (Continued)
4. Tracer Gas Data 43
4.1 Tracer gas release system 43
4.2 Tracer gas sampling system 44
4.3 Tracer analysis system 53
4.4 Tracer gas data tape files 59
4.4.1 Tape file index 60
4.4.2 Tape fi 1 e records 60
4.5 Gas chromatograph calibration data tape files 63
4.5.1 Tape file index 64
4.5.2 Tape file records 65
5. Pilot Balloon Wind Data 67
5.1 Pilot balloon wind system 67
5.2 Pilot balloon wind data tape files 69
5.2.1 Tape file index 70
5.2.2 Tape file records 71
6. Tethersonde and Minisonde Data 74
6.1 Tethersonde and minisonde data systems 74
6.2 Tethersonde and minisonde data tape files 75
6.2.1 Tape file index 75
6.2.2 Tape file records 76
7. EPA Complex Terrain Model Development
SHIS #1 Modelers' Data Archive - 1982 81
7.1 Modelers' data archive 81
7.2 Tracer concentration data 82
7.3 Tracer release information 83
7.4 Meteorological data 84
7.5 Archive structure 91
7.6 Modelers' data tape files 92
7.6.1 Tape file i ndex 93
7.6.2 Tape file records - meteorological data .............. 95
7.6.3 Tape file records - tracer concentration data .... 95
7.7 Conclusion 99
References 101
vi
-------
FIGURES
Number Page
1 Topography of Cinder Cone Butte 6
2 Field experiment layout 12
3 Data collection system configuration 25
4 Tracer gas sampler locations 49
5 Reflection mast design 50
6 Tracer gas chromatograph 55
7 Bag sampling and analysis procedures 57
8 Procedures to obtain tracer gas concentrations 58
-------
TABLES
Number Page
1A Tower Instrumentation and Measures 14
IB Climatronics Instruments Used for Tower
Instrumentation 17
2 Identification of Measurement Codes 18
3 Measurement Range and Resolution Due to Integer
Communication 24
4 Periods of Meteorological Data Collection 32
5 Meteorological Data Set No. 1:
Edited Data with Uncorrected Wind Data Tape
File Numbers - Tape File Numbers 34
6 Meteorological Data Set No. 2:
Edited Data with Corrected Wind Data:
Tape File Numbers 36
7 Data Record Format 37
8 Meteorological Record Types 38
9 Meteorological Data Set No. 1: Sample Printout -
Edited Data with Uncorrected Wind Data 41
10 Meteorological Data Set No. 2: Sample Printout -
Edited Data with Corrected Wind Data 42
11 Tracer Rel ease Data 45
12 Tracer Gas Sampler Network 51
13 Tracer Gas Tape Fi 1 es 60
\/m
-------
TABLES (Continued)
14 Tracer Data Format 61
15 Tracer Concentration Data - Sample Printout 62
16 Gas Chromatograph Calibration Data Tape Files 64
17 Gas Chromatograph Calioration Data Format 65
18 Gas Chromatograph Calibration Data - Sample Printout ... 66
19 Pilot Balloon Wind Tape Files 70
20 Pi 1 ot Bal 1 oon Wind Data Format 71
21 Pilot Balloon Wind Data -
Sample Printout 73
22A Tethersonde Data Tape Fi 1 es 76
223 Minisonde Data Tape Files 76
23 Tethersonde Data Format 77
24 Tethersonde Meteorological Data -
Sample Printout 78
25 Mi ni sonde Data Format 79
26 Minisonde Meteorological Data -
Sample Pri ntout 80
27 Modelers' Data Tape Files 93
28 Modelers' Meteorological Record Data Format - Tower A .. 94
29 Modelers' Meteorological and Tracer Release Data -
Sample Printout 96
30 Modelers' Tracer Concentration Record Data Format 97
31 Modelers' Tracer Concentration Data -
Sample Printout 98
ix
-------
LIST OF SYMBOLS AND ABBREVIATIONS
SYMBOL
b . Scattering coefficient - Nephelometer
C Concentration
Fr Froude number
g Acceleration caused by gravity
h Hill height
H Height of the plume center!ine above the ground over
flat terrain
H Critical dividing streamline height
IX, IY, IZ Turbulence intensities alongwind, crosswind, and
vertical
N Brunt-Vaisala frequency
Q Tracer emission rate
r, 9, z CCB polar coordinate system coordinates
a. Standard deviation of horizontal wind direction
B
a Standard deviation of vertical velocity fluctuations
w
a Standard deviation of alongwind velocity fluctuations
a Standard deviation of crosswind velocity fluctuations
T Average temperature
u Wind speed at source
U, u Uniform wind speed of flow approaching hill
-------
LIST OF SYMBOLS AND ABBREVIATIONS (Continued)
z. Mixed layer height
z Plume release height
ABBREVIATION
AFB Air Force Base
AID Analytical Instrument Development, Inc.
CCB Cinder Cone Butte
CTMD Complex Terrain Model Development Program
ECL Executive Control Language
EPA U.S. Environmental Protection Agency
ERT Environmental Research & Technology, Inc.
FAA Federal Aviation Administration
FMF Fluid Modeling Facility
GC Gas chromatograph
MSL Mean Sea Level
MST Mountain Standard Time
NOAA National Oceanic and Atmospheric Administration
NRTS National Reactor Testing Station
ppb Parts per billion by volume
xi
-------
LIST OF SYMBOLS AND ABBREVIATIONS (Continued)
ppt Parts per trillion by volume
RTD Resistance Thermometric Device
SHIS Small Hill Impaction Study
TRC TRC Environmental Consultants, Inc.
WPL Wave Propagation Laboratory
-------
ACKNOWLEDGEMENTS
This report is partly composed of excerpts from publications and
documents produced by Environmental Research and Technology, Inc. (ERT) the
prime contractor for the Complex Terrain Model Development project, who
compiled the computer data base of magnetic tapes. As referenced in the
text, the First Milestone Report - 1981 by Lavery et al. (1982), was another
important source, as was the Quality Assurance Project Report for Small Hill
Impaction Study No. 1 by Greene and Heisler (1982). Section 7 of this
report is a reproduction of a document by Strimaitis and DiCristofaro of ERT
describing a special Modelers' Data Base they developed from data acquired
at Cinder Cone Butte. The magnetic tape files containing observed and
computed values derived to assist any modeling effort were included in the
computer data base and made accessible along with all other tape files.
All credit for creation of the computer data base and documentation of
the effort must go to the scientists and investigators at ERT. The purpose
of this report was to condense available documentation into one volume that
would serve as a convenient handbook for any investigators who might acquire
and use this valuable computer data base.
XT 11
-------
SECTION 1
INTRODUCTION
1.1 EPA PROGRAM
The extensive development of energy resources, especially in the
mountainous terrain of the western United States, has generated concern
about the resulting impact on air quality (as well as on water and land).
Even in relatively simple situtations, it has been difficult to produce
reliable calculations of atmospheric transport and diffusion. In complex
terrain, mathematical modeling is confounded because the physical processes
are more complicated and meteorological measurements are less
"representative" than in level terrain settings. Responding to this
fundamental problem, the U.S. Environmental Protection Agency (EPA) has
embarked upon the Complex Terrain Model Development (CTMD) Program, a major
program to develop and demonstrate reliable models of atmospheric dispersion
for emissions in mountainous terrain.
An early step in the development of this program was the convening of a
workshop on problems in modeling atmospheric dispersion over complex
terrain. In concert with recommendations of the workshop report (Hovind et
al., 1979), EPA's CTMD Program involves a coordinated effort in mathematical
-------
model development, field experimentation, and scaled physical modeling. The
Program's basic objective is the production of practical models with
demonstrated reliability. Initially the CTMD Program has focused on the
problem of stable plume impaction/interaction with elevated terrain. This
phenomenon was singled out because of the likelihood of relatively high
concentrations and because models that are in use have been challenged
extensively. The approach has been to study stable plume interactions first
in relatively simple terrain settings and subsequently in more complex
situations.
EPA's prime contractor for carrying out the CTMD program is
Environmental Research and Technology, Inc. (ERT). Significant
contributions are also being provided by EPA's Fluid Modeling Facility (FMF)
and by the National Oceanographic and Atmospheric Administration's Wave
Propogation Laboratory (WPL) through their sophisticated measurement
capabilities. A comprehensive tracer field study was carried out on Cinder
Cone Butte (CCB), near Boise, Idaho, during the autumn of 1980'(Small Hill
Impaction Study No. 1, SHIS #1). Based on those data, several models of
plume impaction have been tested and some relatively new modeling concepts
have been introduced (Lavery et al., 1982; Strimaitis et al., 1983).
1.2 OBJECTIVE
The purpose of this report is to describe the data that were collected
in the tracer field study on CCB and to publicize their availability. These
data offer a wealth of information for model development/testing, which is
-------
continuing in the EPA Program. This report is based on three other reports
which provide further details and trace the history and refinement of the
data. The first two are program milestone reports (Lavery et a!., 1982;
Strimaitis et al. , 1983), while the third is a very thorough report on
quality assurance (Greene and Heisler, 1982). In spite of the publication
dates, these documents were written in the order mentioned above. They
should be consulted for details beyond those provided here.
This report describes the setting of CCB, the experimental approach,
and the following data archived on magnetic tape in seven sets of data
files:
Tower (six) wind and temperature measurements (unaltered but
flagged), solar and net radiation at one location, and
nephelometer data;
Tower wind data refined by applied quality assurance procedures;
Tracer gas concentrations and release data;
Winds based on pilot balloon data;
Winds, temperatures, and moisture measured from tethersondes;
Winds, temperatures, and moisture measured from balloon-borne
minisondes; and
The modelers' archive of derived wind and temperature values at
tracer release locations and measured tracer concentrations
(tracer values in this file differ from those in data file 3; here
averages are taken of colocated samplers, reanalyzed samples, and
10-min samples during a given hour).
-------
In the first set of data files, the nephelometer measurements were taken at
three locations near the top of CCB. These data (5-min averages of
backscatter) are listed with Tower B data. A preliminary evaluation of
these measurements indicates that they are qualitatively useful for
determining when and where plume impact occurs. Although lidar measurements
(by WPL) and extensive photography were made of the oil fog plumes, those
data are not available for publication at this time. Pertinent scaled
physical modeling studies by EPA's FMF are being published as they become
available.
-------
SECTION 2
FIELD STUDY AT CINDER CONE BUTTE
2.1 GEOGRAPHIC AND METEOROLOGICAL SETTINGS
Cinder Cone Butte is an isolated small hill in the Snake River Valley,
located about 50 km SSE of Boise and 25 km NW of Mountain Home Air Force
Base, Idaho. Near CCB the main axis of the valley is SE-NW. Mountains
begin to rise sharply at roughly 20 km NE and SW of CCB. Immediately east
of CCB some summer farming is done for potatoes, sugar beets, and grain.
The butte is on the eastern boundary of a National Guard training range.
The soil is mostly sandy and rocky with some grass as tall as 0.5 m and
sparse scrub brush that rarely reaches 2 m.
Figure 1 illustrates the topography of CCB. Contours are shown in 10-m
intervals. The base (zero) contour is at an elevation of 945 m (3100 ft)
above sea level. Both rectangular and polar coordinates, with origins at
the center of CCB, have been used. Cinder Cone Butte is a two-peaked,
roughly axisymmetrical hill, about 95 m high. The nearly circular base is
about 1 km in diameter. The upper part of the butte has slopes that are
generally around 25 degrees. There are several roads near the base of CCB
and one to the top to service a permanent FAA tower on the northern peak.
Several photographs of CCB are presented in the report by Lavery et al.
(1982).
-------
CINDER CONE BUTTE, ID
TOWER "A" 3S2DEG 2.3KM
W I
: TOWER'
... , I -/F A A TOWER / I I TOWER
i'1
\ \
\ TOWER "E" TOWER "B"
CONTOUR INTERVAL = 10M
OATUM=3100FT CONTOUR (945MIMSL
METERS
150
300 450
600
Figure 1. Topography of Cinder Cone Butte
-------
General meteorological information indicated a high frequency
throughout the year of nights with marked temperature inversions. The
autumn season was targeted for the field study in order to operate during
longer nights but before the normal onset of winter storms. Summaries of
hourly wind observations from nearby Mountain Home AFB indicated predominant
surface air flow from either the southeast or northwest, along the valley's
axis.
2.2 EXPERIMENTAL DESIGN
The field study was designed so that tracer releases from a platform
lifted by a large crane would impact on or interact with CCB. Sixty
sequential bag samplers were deployed over the butte to collect up to twelve
1-h samples each during each experiment. Twenty samplers were programmed to
collect 10-min samples. Seventy of these samplers were in fixed locations
and ten were moved according to meteorological conditions. In addition, a
small number of samplers were used for background determination, for
colocation, and for vertical sampling at two or three levels on four 6-m
"reflection" masts (operated during experiments 203 to 218). Except for 3-
or 6-m levels on masts, the sampler intakes were nominally at 1 m above
local ground.
During Phase I, September 16 to 27, 1980, ten "shakedown" and strictly
visualization experiments were conducted. During Phase II, October 16 to
November 12, 1980, eighteen quantitative tracer experiments (numbers 201 to
-------
218) were carried out. Each lasted 8 h and was conducted at night or early
morning. The main tracer was sulfur-hexafluoride (SFfi). An oil fog, for
visualization of the plume, was generated within 1 m (horizontally) of the
SFg release point. A second tracer, Freon 13B1 (CF-Br) was released
simultaneously at a different elevation (usually lower) in nine experiments
(numbers 208, 210, 211, and 213 to 218). The tracer gas concentrations in
the bags were measured by gas chromatography. The lowest quantifiable
limits of detection were 20 parts per trillion (ppt) of SFg and 220 ppt of
Freon. More than 14,000 bag samples (1-h and 10-min) were analyzed. In
cases where it was clearly demonstrated that the tracer gas(es) did not
reach the bag samplers, most of those bags were not analyzed. In one tracer
experiment (number 212) the wind was so variable that it was not possible to
align the release system upwind of the sampler array. The visible oil fog
plume was never observed to hit the hill. No concentrations of SFC were
b
above 5 ppt, indicating a lack of background contamination. No
meteorological data are included in the files for experiment 212. The
tracer release and analysis procedures and quality assurance program are
discussed extensively in the three documents mentioned previously and in
other appropriate sections of this report.
A lidar (operated by WPL) took backscatter measurements through
vertical sections of the oil fog plume but those data are not included in
the archive described here. Nephelometer measurements were collected at
three locations near the top of CCB. These data (5-min averages) are
qualitatively useful for determining plume impact. An attempt was made to
collect nephelometer data from instruments hung from a crane, but the data
were not usable and are not included in the archive.
-------
The arrangement of the meteorological monitoring equipment is shown in
Figure 1. There were six towers. Tower A, 150 ni tall, was instrumented at
eight levels with a resistance thermometric device, at five of those levels
with triaxial propeller anemometers, and at three of the latter levels with
cup-and-vane wind instruments. Fast-response thermistors were operated at
two levels for determining standard deviations of temperature. Tower A was
located 2 km north of CCB, where scaled modeling studies at EPA's FMF had
shown it to be outside the region of flow disturbed by the butte in the
frequent stable southeasterly or northwesterly flows. A pyranometer and a
net radiometer were located 40 m WSW of the tower base.
Tower B, 30 m tall, was located on top of the south knoll of CCB. It
was instrumented at three levels with resistance thermometric devices and
triaxial propellers and also at the two upper levels with cup-and-vane
instruments.
Towers C through F, 10 m tall, were located on CCB near the 70-m
contour, except F was about 10 m lower. Each was instrumented at two levels
with resistance thermometric devices and triaxial propellers and also at the
upper level with cup-and-vane instruments. The sampling frequency for all
tower-mounted meteorological instruments was 4 Hz, which was used in real
time to calculate 5-min and 1-h average values for automatic storage in the
computer and display at the command post. Differences in the operating
characteristics of the two wind systems and refinement of the measurements
are discussed in detail in the quality assurance document mentioned earlier
(Greene and Heisler, 1982). The base of the FAA tower was used to install a
nephelometer instrument.
-------
A tethersonde was operated at one of two sites (depending on wind
direction) located about 1.3 km NW and SE from the butte center; the sondes
were released from the more upwind of the two sites. Ascent and descent
measurements were made of temperature, pressure, wind direction and speed,
and moisture, normally at 1-h intervals during tracer releases. When winds
were too fast to operate the tethersonde or the tethersonde system was
inoperable, data were obtained from minisondes tracked by theodolite. In
addition, double-theodolite pibals were taken at hourly intervals between
tethersonde soundings.
Temperature and wind instruments 1 m above ground were deployed at five
locations along the east-facing draw on CCB in order to document low-level
drainage winds. However, the records and some instruments were stolen from
the site and have not been recovered.
In addition to the lidar, the WPL also operated a frequency-modulated,
continuous-wave radar to determine winds aloft and two monostatic acoustic
sounders to determine stability structure aloft. However, these data are
not yet available.
Further details on the CCB data archive are given in other appropriate
sections of this report and in the reports by Lavery et al. (1982),
Strimaitis et al. (1983), and Greene and Heisler (1982).
10
-------
SECTION 3
TOWER METEOROLOGICAL DATA
3.1 FIXED METEOROLOGICAL NETWORK
Six meteorological towers, designated A through F, were deployed at the
CCB site as shown in Figure 2 (Lavery et al., 1982). Tower A, the 150-m
"profile" tower, was located about 2 km north of the butte, where tow-tank
experiments at EPA's FMF had shown it would be outside the region of flow
disturbed by the butte in the predominantly NW and ESE nocturnal winds of
the Snake River Basin. Tower B, 30 m high, was located on the top of the
south peak of the butte. Towers C through F were 10 m h'igh, based within
the 62- to 78-m contours (see Figure 1), and deployed approximately at 90
degree intervals from the center of the butte.
Meteorological data are identified in the computer data base by a code
of four to six characters, the first of which identifies the tower, A
through F, the second and third the level on the tower, and the remainder
the type of data. Because the heights of the nine instrument levels on
Tower A extend to 100 m and beyond, the levels are coded as simple
identifiers, LO to L8. For the other towers, the level codes are the level
heights in meters above the tower base. Table 1 (Lavery et al. , 1982)
11
-------
S- .
O 0)
i
o
ur>
0)
-C V
P J-
(0
Ul
3 O
O
&U
r- 0)
P O
C 4->
E E
i- i
S- O
0) PO
Q.
X
-------
illustrates the arrangement of instrumentation on each tower, and Table 2
(Lavery et al., 1982) identifies the codes used for each type of
meteorological data measured.
Tower A was a Unarco-Rohn Model 80 with a Tower Systems elevator that
lifted the instrument booms to 10, 40, 80 and 150 m. It was instrumented
with temperature sensors at eight levels and wind sensors at five levels. A
net ^adiometer and pyranometer were located about 1 m above the ground and
roughly 40 m WSW of the tower base. The tower was surrounded by fairly
sparse prairie grass about 50 cm high except near its base, where the grass
had been trodden down by people and vehicles during installation and
maintenance.
There were two temperature sensors on Tower A at both the 10- and 150-m
levels. One was a standard platinum resistance thermometric device (RTD),
giving temperature differences from the 2-m level, and one fast response
thermistor bead, installed for the purpose of calculating the standard
deviation cf temperature. Both sensors at either level were mounted in the
same aspirated shield. All aspirators were installed with the intake end
facing north.
Tower B, a Unarco-Rohn Model 55 on top of the butte, was instrumented
with wind and temperature sensors at the 2-, 10-, and 30-m levels. The area
around its base was fairly smooth, soft rock and soil without vegetation.
Very sparse grass, 10 to 20 cm long, and small stones covered the slopes of
the nil! away from this tower.
The other four towers, C through F, were identically instrumented
Unarco-Rohn Model 45's; wind and temperature sensors were mounted at 2- and
13
-------
TABLE 1A. TOWER INSTRUMENTATION AND MEASURES
Site
Instrument*
Direct Measures**
Derived Measures*
Tower A
Level 0 (1 m)
Level I (2 m)
Level 2 (10 m)
Level 3 (20 m)
Level 4 (40 m)
Level 5 (60 m)
Level 6 (80 m)
Level 7 (100 m)
Level 8 (150 m)
Pyranometer
Net radiometer
Triaxial props
Cup-and-vane
RTD
Insolation
Net radiation
U, V, W, IX, IY, IZ
UX, VX
T
Triaxial props U, V, W, IX, IY, IZ
Cup-and-vane UX, VX, MS, MD, SD
RTD TD
Fast bead thermistor T, ST
RTD
Triaxial props
RTD
RTD
Triaxial props
RTD
RTD
Triaxial props
Cup-and-vane
RTD
U, V, W, IX, IY, IZ
TD
U, V, W, IX, IY, IZ
TD
U, V, W, IX, IY, IZ
UX, VX
TD
Fast bead thermistor T, ST
WS, WD
SP, DR
WS, WD
SP, DR
TC
WS, WD
TC
WS, WD
TC
WS, WD
SP, DR
TC
(continued)
14
-------
TABLE 1A. (Continued)
Site
Instrument*
Direct Measures**
Derived Measures*
Tower B
Level 2 (2 m) Triaxial props
RTD
Level 10 (10 m) Triaxial props
Cup-and-vane
RTD
Level 30 (30 m) Triaxial props
Cup-and-vane
RTD
(1.5 m) Nephelometer
U, V, W, IX, IY, IZ WS, WD
T
U, V, W, IX, IY, IZ WS, WD
UX, VX SP, DR
TD TC
U, V, W, IX, IY, IZ WS, WD
UX, VX SP, DR
TD TC
ANEPH
Towers C, D, E, F
Level 2 (2 m) Triaxial props
RTD
Level 10 (10 m) Triaxial props
Cup-and-vane
RTD
Tower D
(1.5 m) Nephelometer
FAA Tower (north peak)
(1.5 m) Nephelometer
U, V, W, IX, IY, IZ WS, WD
T
U, V, W, IX, IY, IZ WS, WD
UX, VX SP, DR
TD TC
CNEPH
BNEPH
(continued)
15
-------
TABLE 1A. (Continued)
*A11 temperature sensors were mounted in aspirated radiation shields; an RTD
is a Resistance Thermometric Device.
**"Direct" measures are those calculated by the data station microprocessors
from the outputs of the instrument translators. The turbulence data (IX,
IY, IZ, SD, ST) were calculated for both 5-min and 1-h data periods. All
direct measures were updated once per min by the data stations. UX and
VX are the westerly and southerly components of the wind calculated from
the cup-and-vane outputs at the 4 Hz sampling frequency.
+"Derived" measures are those calculated by the data collector computer from
the 5-min averages provided by the data stations. These derived measures
comprise 5-min average values of the measures indicated as well as 1-h
averages of all direct and derived measures except the turbulence data,
1-h averages of which were calculated by the data stations.
Tower elevations - datum = 944.9 m MSL
Tower A = -03.5 m
Tower B = 98.0 m
Tower C = 69.0 m
Tower D = 69.4 m
Tower E = 78.8 m
Tower F = 61.8 m
FAA Tower = 94.9 m
16
-------
TABLE IB. CLIMATRONICS INSTRUMENTS USED FOR TOWER INSTRUMENTATION
Type
Model no.
UVW component anemometer - Triaxial props
(23 cm styrofoam propellers)
Cup-and-vane windset
Platinum resistance thermometric device - RTD
Fast response bead thermistor probe
Motor aspirated temperature probe
Solar radiation sensor - pyranometer
Net radiation sensor - net radiometer
WC-13
F460
100826
100093-4
TS-10
1000484
100881
17
-------
TABLE 2. IDENTIFICATION OF MEASUREMENT CODES
Code Units
SR ly/nri n
NR ly/min
U, V, W m/s
WS m/s
WD degrees
IX, IY, IZ percent
Instrument
pyranometer
net radiometer
propellers
propellers
propellers
propel lers
Measurement
solar radiation
net radiation
westerly, southerly
vertical wind
components
vector resultant wi
speed
vector resultant wi
direction
>
nd
nd
alongwind, crosswind,
UX, VX
SP
DR
T
TD
TC
TF
m/s
m/s
degrees
°C
°C
°C
°C
cup-
cup-
cup-
RTD
RTD
T &
bea<
bead thermistor
and vertical
intensity of
turbulence
westerly and southerly
wind components
vector resultant wind
speed
vector resultant wind
direction
dry bulb temperature
temperature difference
from 2-m level
calculated temperature
dry bulb temperature
(continued)
18
-------
TABLE 2. (Continued)
Code
ST
MS
MD
SD
ANEPH,
BNEPH,
CNEPH
Units
°C
m/s
degrees
degrees
bSCAT'
nT1
Instrument
bead thermistor
cups
vane
vane
nephelometer
Measurement
standard deviation of
dry bulb temperature
scalar mean wind
speed
scalar mean wind
direction
standard deviation of
wind direction
scattering coefficient
10 m. Tower C, to the NE of the butte's center, was based on a steep slope
covered with sparse, short grass, small stones and soil. Tower D was to the
SE of the center of the butte and located approximately 20 m north of the
crest of the ridge coming down from Tower B. Here the grass was somewhat
longer (20 to 30 cm) and thicker than it was around the other towers on the
butte, particularly on the steep slope NE of this tower. The ground around
Tower E, on the other hand, was nearly free of vegetation except for some
sagebrush and large rocks up to 0.5 m in diameter and half that in height.
Tower F sat on a small, fairly level spot on the ridge to the S of the NW
draw. The ground around it was quite variable, mostly short sagebrush below
it and to the N on the slope of the draw, with short grass and small rocks
in other directions.
19
-------
Because the sites of the four 10-m towers were sloped, and the
different wind instruments at the 10-m level were separated from one another
by crossarms, the heights of the sensors above the ground must be regarded
as "nominal".
Sites were established for three nephelometer instruments to detect the
impact of the oil-fog plume on CCB. Nephelometer-A was located at the base
of Tower 8 on the south peak, nephelometer-B on the north peak near the FAA
Tower, and nephelometer-C at the base of Tower D. All instrument sample
intakes were approximately 1.5 m above the ground. Data were recorded as a
scattering coefficient (b t) and placed in Tower B data tape files.
Certain of the measures in the CCB data base require explanation. With
the exception of measurements made at 10 m on Tower A, all wind speeds and
directions are horizontal vector resultant values whether derived from
triaxial propeller anemometers or from cup-and-vane sets. These measures
are derived from the propellers simply by averaging the u (westerly) and v
(southerly) component values from the respective propellers. The
measurements from the cup-and-vane sets were resolved into westerly and
southerly components as the data were taken at the 4-Hz sampling frequency.
Only at the 10-m level on Tower A were the cup mean speed (MS) and vane mean
direction (MD) recorded in the traditional, single-sensor fashion. Standard
deviation of vane direction (SD) was also calculated only for this site.
The other measures requiring explanation are the intensities of
turbulence, IX, IY, and II, respectively, the alongwind, crosswind, and
vertical intensities. These are approximately au/0, av/U, and ow/U, where U
is the mean horizontal wind speed. The CCB measures, however, are
20
-------
calculated from the signals of the propeller anemometers by means of the
following algorithms:
TY - l r Iy + 2ZUZVZUV -, .
IX - ( [ 2 J} - U
(lu) + (Iv)
IY = {1 [Zu2 + Zv2 - (^)2 lu2 + iZv)2 Zv2 + 2ZUZVZUV
+ Zv^
IZ =
2 2
where U = ^ N Vis the vector resultant mean wind speed, N is the
number of samples in the calculations, and u, v, and w are the instantaneous
wind component speeds from the triaxial propeller anemometers.
3.1.1 Data Acquisition System
Each of the 89 instrument transmitters in the fixed meteorological
network was sampled four times a second. This sampling was done by nine
microprocessors (ERT Data Station Model No. DS-00) installed with the
instrument translators and power supplies in the two shelters provided, one
to the west of the 150-m tower and one between the two peaks of the butte,
which served the five shorter towers. The microprocessors converted the 0-
to 5- volt instrument outputs to digital values, calculated the intensities
of turbulence and standard deviations of temperature and wind direction,
resolved the cup-and-vane data into westerly and southerly components, and
-------
made 5-min averages of 142 measurements. All calculations were updated once
per minute by the microprocessors.
The 5-min data were collected by a minicomputer located in ERT's base
trailer approximately 3 km E of the butte. This data collector was a Data
General Nova 3 equipped with two disk drives, a tape drive, and a Digital
Equipment Corporation Decwriter as interactive console. The data collector
polled each microprocessor channel every 5 min for the 5-min data and every
hour for the turbulence quantities. It did parity and range checks on the
data received (flagging faulty data), resolved the vector component wind
data into speeds and directions, calculated temperatures at elevated tower
locations by adding AT's to 2-m T's, wrote a selected subset of 5-min data
on the console for experiment control, accumulated hour averages of the 142
measurements provided by the microprocessors, and wrote the data to disk.
On the hour it also requested the 1-h turbulence quantities from the
microprocessors and made the same calculations from its accumulated hour
averages that it did from each 5-min scan of data.
The scan of the microprocessors' 142 channels took longer than 1 min so
that all the data associated with a particular time are not necessarily
representative of the same 5-min period. Consequently a calculated
temperature (TC) may sometimes not equal the sum of the 2-m temperature (T)
and the temperature difference (TD) because the T value was requested again
by the data collector just before the calculation of the higher level TC.
At least once per hour (every 20 minutes during the later experiments)
the Data General M-600 at ERT's office in Concord, MA, requested the 5-min
data (and 1-h data if appropriate) from the Nova and wrote it to disk for
archive.
22
-------
The microprocessors in the instrument shelters and the data collector
in the trailer communicated via cables. The microprocessors transmitted
their data as integers between 0 and 1,023. The resolution of the
transmitted data was therefore slightly better than 0.1% of full range. The
range and resolution of each of the measurements is listed in Table 3
(Lavery et a!., 1982). The data collector and the computer in Concord
communicated via telephone line. Figure 3 (Lavery et al. , 1982) is a
schematic of the data acquisition system.
3.1.2 Quality Assurance Plan
ERT's quality assurance plan (Greene and Heisler, 1982) for the fixed
meteorological measurements was based on careful, documented calibration of
the instruments before installation; a calibration check shortly before
Phase 2, the complete tracer experiments that started in the middle of
October; a calibration check at take-down in mid-November; and the
performance audit by TRC. Additional automatic quality control (QC) checks
for parity and range were done by the data collector computer in real time
(see Section 3.1.1 above).
The instruments could not be completely checked out and calibrated, as
planned in the ERT laboratory in Ft. Collins prior to installation because
the period between receipt of the instruments from the manufacturer and
startup was too brief. The field calibration done at installation, however,
incorporated most of the laboratory procedures. The output of the cup and
propeller anemometers was checked when spun with a synchronous motor and
when held stationary, the bearings of the vanes and anemometers were checked
23
-------
c
o
^
p
Z 3
O i
r-< O
h- co
< CU
1 1
z
o
O
0*
UJ
CD
UJ
H-
Z
i i
0
r-
UJ
13 CU
o en
c
z re
0 CC
1
0
co
UJ
Q£
Q
Z
^
UJ
CD
Z
^
f^
CD
rH
O
O
C
r-
E
>^
!_
O
CM
O
P
O
o
c
o
F
p
re
*-
73
re
5-
S-
re
^
o
10
Q£
CO
c
^
E
*^
r~~
^"
CM
o
o
c
f
E
^(^
fmmm
O
*
CM
O
P
m
.
o
1
c
o
"^
-p
re
*(_
T3
re
P
CU
c
Q£
Z
U)
^^
E
G\
^"
O
CO
^*x.
E
lO
CM
^-
O
-P
m
CM
1
/^^
10
Q.
O
i-
Q.
v^
(/)
+J
C
CU
c
o
a.
E
O
u
TJ
C
r
3.
r
re
c
o
N
r»
^
O
.c
>
M
p
c
CU
u
CO S_
v. cu
E a.
O 00
CM CD
O O
P
C
CU
u
CO 5-
V. CU
E a.
0 0
rH 0
4- rH
O O
0 0
1
(A
CU
r^. 4J
CL -r-
O 10
s- c
a. cu
c
P -r-
c
CU CU
c u
o c
a. cu
E
O 3
U J3
i.
TJ 3
C +J
r-
5 r~-
re
r- -P
re c
U 0
r- N
P T-
S- i.
CU O
> J=
5.
*""
M
X
S >~l
C
Q
U
5-
CU
a.
en
co
o
4-*
c
CU
u
5-
CU
a.
o
^3"
o
.p
o
^^
^->
r-
10
C
01
p
c
^
CU
u
c
CU
^3
^
i.
3
-P
r
re
u
r~
P
S-
0)
>
M
CO
\
E 0
o
00 CO
cr* r^*
o o
CO
"s^
E 0
o
o o
m m
+ -r
O O
P P
O O
m co
i i
/*^
^>
od
CJ
CO
^J
c
CU
c
o
a. cu
E S-
O 3
U P
re
c cu
r- a.
3 E
cu
re
P ja
C r-
O 3
M J3
^ ^^
O S-
JT -a
x:
A
x:
3 h-
f ^
O
0
CM
O
CJ
o
m
rH
O
P
m
0)
M-
.^
a
cu
3
+J
re
cu
Q.
cu
Q
r-
O
O
CO
o
CJ
o
o
m
o
-p
o
co
i
a>
^
p
re
cu
a.
cu
P
cu
CO
c
o
o.
W)
cu
f-
p
CO
re
14-
u_
t~
CJ
o
rH
O
CJ
O
CO
rH
O
O
cu
i-
3
-P
re
i.
cu
o.
E
CU
p
<4_
o
c
o
+J
re
r-
>
CU
a
T!3
fi-»
re
a
c
re
to
r-
CO
^
C71
O
(A
>X,
E
O
in
o
o
TJ
cu
cu
a.
IA
o
c
^>
3
c
re
cu
E
i.
re
re
u
ifl
co
CO (A
cu cu
cu cu
i. i.
O) O5
0) CU
o -a
00 CO
CM rH
m o
(/> (A
CU CU
CU CU
i. i.
O) O)
cu cu
a -a
O CO
* o
in rH
o o
p p
o o
c
o
p
CJ
cu
i.
I
XJ
c cu
o c
r- re
p >
u
CU
-------An error occurred while trying to OCR this image.
-------
for proper performance, and vane output was checked when pointed towards and
away from known landmarks. Temperature sensors were immersed in cold and
warm stirred water baths whose temperatures were measured with NBS-traceable
precision thermometers. All instrument translator cards were tuned to give
correct output voltages at various calibration points.
3.1.3 Data Editing
ERT has performed two types of editing on the data from the fixed
meteorological network. The first was an examination of each 5-min value to
flag data that were identified as incorrectly transmitted to the data
collector or that were taken from a malfunctioning sensor. The second
editing process-was correction of the data for consistent and significant
calibration errors, for misorientation of wind sets to true north as
determined by TRC's audit, for instrument response characteristics, and for
effects of tower wakes.
The first data editing, or validation, was largely necessitated by two
types of equipment failures. Shielded cable had been specified by ERT for
the data links from the microprocessors in the instrument shelters to the
data collector, but it could not be obtained from vendors in time for the
startup of the experiments, and unshielded cable was used. This resulted in
frequent communication errors in spite of parity checks done on both the
data requests from the data collector and the data transmissions from the
microprocessors. The second type of equipment failure was a loss of
response of a UVW component propeller. Sometimes a propeller would stop
turning almost completely, in which case the fault was easy to identify.
26
-------
More frequently, however, che failure was more subtle - a slight
"stiffening" of the instrument identified by an unusually low ratio of
crosswind intensity of turbulence, IY, to alongwind intensity of turbulence,
IX, or by changes in the relative values of U or V components with respect
to an adjacent cup-and-vane set in similar conditions of wind speed and
direction but separated by an hour or more.
The communication problems resulting from the unshielded cable caused
two major types cf errors. The first type was a miscommunication from a
data station to the data collector that was not always picked up by the
parity checking on the transmission. Such an error resulted in a value that
looked peculiar in the time series of values for the measure affected. From
the redundant wind measurements (both cup-and-vane and UVW propellers),
errors of this type could be fairly easily verified for wind speed and
direction except at the 40- and 80-m levels of Tower A, where propellers
were alone and vertically separated by 30 m or more from the nearest source
of data for comparison. Because of the strong thermal layering during many
experiments, it was often unfeasible to verify a communications "hit" by
comparison to these levels, and the determination that a value was suspect
or in error depended entirely on whether it was unreasonable or out of place
in the time series. Calculated temperatures at 10 m and 150 m could be
validated by comparison with the values from the fast-response thermistors
at these sites. Temperatures and temperature differences at other heights
on Tower A were validated by comparison with temperatures above and below
the height being validated.
27
-------
Fortunately, few errors of communication from the data stations to the
data collector resulted in values that were in the range of possible values.
Most were recognized as faulty by the data collector and identified by an
"M" flag.
A more difficult communication problem occurred in the data requests
from the data collector to the data stations. A request for a wind
component might, for example, be received as a request for a temperature,
and the temperature would therefore be returned to the data collector and
put into the data base as the wind value. All measured values were
transmitted to the data collector as integers (called "counts") between 0
and 1,023 inclusively. The data collector converted them to proper
engineering units by interpolation in the range of the measure. A
temperature transmitted in error as a wind component would therefore not
appear in the data base as the value of temperature that was sent but rather
as the value of the wind component appropriate to the number of counts
corresponding to the temperature. Consequently, one could look through all
the data for the 5-min scan in which the suspect value occurred for another
measure value that had the same associated counts. If such a measure was
found, the suspect value was regarded as bad.
To expedite the time-consuming error check through all measurements
taken during the experiments, the 5-min data were retrieved from the data
base in time-series files for each experiment. In general, each of these
files included all the 5-min measures for one level on a tower; there are 19
files per experiment in this "edit" format. A set of data flags for
identifying the quality of the data was established as follows:
28
-------
" " (blank): Both the editor and the data system concur that the
value is valid.
"M" (missing): Both the editor and the data system concur the value
is invalid.
"U" (unavailable): The value is unavailable because of data collector
or data station failure.
"B" (bad): The editor believes the value is invalid but the data
system did not catch the error; this flag is therefore associated
either with instrument malfunction or communication problems.
"R" (restored): The editor believes the value is valid although the
data system had flagged it "M".
"C" (calculated): The editor calculated a derived measure (WD, WS,
SP, DR, TC), usually from "R" flagged values; the only exceptions
are nonzero solar radiation values in a string of zeros at night,
for which a 0.000 was inserted and flagged "C".
"S" (suspect): The editor believes the data are somewhat in error but
cannot confirm either an instrument malfunction or communication
failure.
"L" (at limit): The measure is at the upper limit of its range and
the "true" value exceeds that shown. The instrument ranges were
not themselves exceeded during the experiment, and this flag is
necessary only for the turbulence data (IX, IY, IZ, SD) in very
light and variable winds.
29
-------
No data have been estimated and inserted into the data base.
In this validation editing, ERT tried to maintain a balance between the
premise that all data are potentially valid and the premise that no data are
above suspicion. Consequently, if no instrument failure or communications
error could be verified, a value was regarded as valid unless it appeared to
be unreasonable with respect to comparable values adjoining it in time and
space. This is generally not a difficult judgement to make, but in some
situations a value may look peculiar but not completely unreasonable and
might indicate a significant phenomenon. Such data are left unflagged if
they will not be misleading and are flagged "S" if they are substantially
removed from the general trend.
The different characteristics of propeller wind sets and cup-and-vane
systems are well demonstrated in the CCB data. In general, the vector
resultant wind speed (WS) from the propellers was less than the vector
resultant wind speed (SP) from a cup-and-vane set at the same location. The
ratio of WS to SP decreased from 0.8 to 0.9 in high-speed smooth flows down
to 0.5 or less in light and variable winds. In near-calm conditions, the
props were observed to be more responsive to gentle puffs than the vanes, so
that a 5-min wind direction (WD) and wind speed (WS) resolved from the props
might be 175° at 0.2 m/sec, whereas the corresponding cup-and-vane direction
and speed might be 245° at 0.5 m/sec. Both these pairs of wind measurements
may appear in the validated data .vithout any error flag because there was no
indication of instrument malfunction or communication error. The
differences between the measurements are attributable to the differences in
the instruments.
30
-------
Similarly, the response of propeller sansors is direction-dependant.
Often the difference between WD and DR at the site changed markedly when WD
passed through a cardinal direction such as 0°, 45°, 90°, or 135°. Again,
the measures were both retained as valid in the data base.
The differences between the speeds and directions from the two kinds of
instruments show general consistency with the differences anticipated as a
result of the departure of the UVW systems from the cosine response curve,
as discussed by Horst (1973). Furthermore, the horizontal intensities of
turbulence IX and IY tend to become more nearly equal when the average angle
of attack of the wind is approximately equal on both propellers (i.e.,
directions near 45°, 135°, 225°, 315°), whereas IX tends to exceed IY when
the average angles of attack are substantially different (i.e., directions
near 0°, 90°, 180°, 270°). This consistency suggested that the quality of
the UVW data might significantly improve if corrections were applied similar
to those described by Horst (1973), which were derived from comparisons of
R. M. Young propeller data and sonic anemometer data. Although ERT was
unable to find any similar comparative analysis of data for the Climatronics
is! "*
system, corrections for .moncosine response were applied to all data from
wind component propeller sensors, and a separate file of wind data was
produced.
3.1.4 Periods Of Data Collection
Table 4 illustrates the dates and times of the experiments and the
concurrent periods of meteorological tower data collection. No
meteorological tower data collected during experiment 212, November 11,
31
-------
1980, 1700 to 2400 MST have been included in the data base because the
tracer gas never hit the sampler array on the butte.
TABLE 4. PERIODS OF METEOROLOGICAL DATA COLLECTION
Experiment
no.
201
202
203
204
205
206
207
208
209
210
211
213
214
215
216
217
218
1980
Date
10/16
10/17
10/20
10/21
10/23
10/24
10/25
10/27
10/28
10/30
10/31
11/04
11/05
11/06
11/09
11/10
11/12
Experiment
hours (MST)
1700 to 2300
1700 to 2300
0000 to 0800
0000 to 0800
0000 to 0800
0000 to 0800
0000 to 0800
1700 to 0100
1700 to 0100
0000 to 0740
0000 to 0800
0000 to 0800
0200 to 1000
0000 to 0600
0000 to 0700
0200 to 1000
0200 to 1000
Stability
E
E
E to F
E to F
E
E
E to F
E to F
F
E to F
E to F
E to F
E to F
E to F
E to F
E to F
E
NOTE: NO DATA COLLECTED FOR EXPERIMENT 212.
3.2 TOWER METEOROLOGICAL DATA TAPE FILES
Data are stored at the National Computer Center, Environmental Research
Center, Research Triangle Park, North Carolina on Sperry UNIVAC 1100/83
systems magnetic tape, nine track, odd parity, ASCII-quarter word mode,
density 6250 BPI, tape number 004700. Record length is 132 characters, and
the block size is 1320 words or 40 records per block. Each file has three
blocks. UNIVAC users may assign the tape, @ASG,T CCB.U9S//////Q,004700 using
32
-------
UNIVAC Executive Control Language (ECL). Upon request, copies can be
furnished and translated into formats acceptable to any computer using nine
track tape drives.
3.2.1 Meteorological Data Tape File Index
Two sets of meteorological data files are recorded on tape number
004700. The first set of files, numbers 1 to 323, are edited but
uncorrected data; data editing procedures and flags were performed as
described. Table 5 illustrates hew individual tape files are related to
tower sites, record types, and experiment number in the first set of files.
The second set of files, numbers 324 to 612, are derived from the same
wind speed and direction observations as set number one, except that the
data have been corrected to account for audited misalignment of wind sets,
for consistent errors in instrument calibration, for noncosine response of
the wind component propeller sensors, and for the effect of tower wakes on
wind speeds.
Wind speed and directions from the Climatronics F460 cup-and-vane
anemometers were corrected for erroneous calibration, misalignment to true
north, and mean nonlinearity in vane response. Wind speeds and directions
derived from the UVW propeller anemometers were corrected for noncosine
response, misalignment to true north, and consistent calibration errors that
were greater than the resolution of the measurement provided by the data
acquisition system. In addition, corrections were applied to wind speeds
derived from both types of wind instruments to account for tower wakes.
These corrections result in substantially improved correspondence between
speed and direction data from the two types of wind sets.
33
-------
LU
m
z
3
Z
LU
1
HH
u_
LU
0_
^
1
|
^3
^
0
u
C£
.
o
c
J_)
c
0)
E
r
S-
cu
Q.
X
LU
LU
S- S
3 0
O rH
r-
a
^ E
CD
3 o
0 (H
r-
0
S- E
OJ
3 0
O rH
r-
ca
i- E
CD
3 o
O CO
r
rH rH
CM CO
rH
rH CSJ
rH
«d- rH
rH
CO O
rH
CM cn
rH 00
p*^ r*^"
to to
in m
«3- «*
co co
CM CM
rH rH
rH
O
CSJ
oor-»toin*j-cocMrHocnoor tom*}-co
comr*»cnrHrOLnr«..cnocM<3"tocoocM
rHrHrHrHrHCMCMCSJCSlCSJCOCO
f-»tr>in«3-cocsi-Hocnoor-'ioin*3-cocM
coinr»«cnrHcoinr^ooocM'^-toooocM
r-lrHrHrHrHCMCMCSJCMCMCOCO
tOinJCOCO
m^-cocsjr-iocnQOr^ioin^rcocsjrHo
COinr>.CnrHCO«tlOOOO'.'M<3-tOCOOCM
rHrHrHr-trHCMCMCMCSJCSJCOCO
<4-cocsjrHocnoor^-toin^tOLn
co«*tocoocM^t-tocoocMcoinr>.cnrH
rHrHr-lrHrHCMCMCMCMCMCMCO
cnoor*-tom«!i-cocMrHocj'icor»-toin'3-
CM«i-tOCOOCSJ^lOOOOrHCOinr^CnrH
rHrHr-lrHr-tCMCSJCMCMCMCMCO
cOf-»iOin^'cocMrHocncoPN-tom*j-co
CM«^-tOCOOCM^tDOOCnrHCOinr~«-Cnr i
r-lr Ir-tr-trHrHCsJCsJCSJCSJCSJCO
p*.LDin«tcocMrHocnoors«-toin«sfcocsj
CM^J-tOOOOCM^tDI^-CnrHCOini ^CnrH
rHrHrHrHrHrHCMCMCSJCsJCSJCO
toin^-cocsjr-tocnoor^toin'srcocMr-t
csj«3-toooocM«3-inr-«cnrHcomr>.cnrH
rHrHrHrHrHrHCMCMCMCMCMCO
in«a-COCSlr-IOCnoOr>.U3in<3-cOCMrHO
csj«4-tocoocsicoinr~>cn Hcoir)r^cnrH
rHrHrHrHrHrHCMCMCMCMCMCO
Nf-COCSJrHOCnCOr^tDin^COCMrHOCn
csj<*toooorHcoinr*»cnrHcoinr^cno
rHrHrHrHr-tr-tCMCMCMCMCsjoo
COCMrHOCnoOr^tDin^'COCMrHOCnCO
CM«ttooocnrHcomr>-cnrHcoLnr»-ooo
rHr-tr-trHr-t'NJCMCMs'MCMCO
CMr-tocnoor>-toin'^-cocsjrHocnoor^
cM^-tor^cnrHcoLnr^cnrHcomtocoo
rHrHrHrHrHCSJCMCMCMCMCO
rHocnoor^toin<*cocMrHocnoor«vto
csj^ini^cnr-icoinp^cnrHco^-toooo
rHrHrHrHrHCMCMCSJCSjCMCO
ocncor~>tom<3-cocMr-tocnoor^toin
csjcoinr^cnr-icoinr^cnrHCM
-------
No corrections were made to the temperature data because of
inconsistencies between results of two independent audits. The differences
between audit or calibration values and output measurements were quite small
for most of the temperature sensors. All but two indicated errors were less
than 0.2 °C in magnitude, and it appeared that the magnitude of the errors
was close to the resolution of the auditing procedures.
No corrections were made to the turbulence intensity data. Certainly
these data are in error -because of the response characteristics of the
propeller sensors, but no satisfactory justifiable corrections were found to
apply. In most cases, the intensity of turbulence data is probably too low.
Other errors that remain in the corrected data are the effects of the
wake of one instrument on another and the effects of tower wakes on
turbulence measurements and wind direction. Installation of propeller
sensors on the north side of the towers resulted in a region between 90° and
125° in which the propeller speed was about 40% of cup anemometer speed.
Data users should be advised to give precedence where possible to wind
measurements from instruments that are clearly out of the wakes. A
comprehensive discussion of this problem is presented in Greene (1982).
The second set of meteorological data files, numbers 324 to 612,
contains a corrected version of wind data that are in the first set of tape
files. Since no corrections were applied to temperature data, record type 1
from Tower A and record type 4 from Tower B were omitted from the second set
of files. Table 6 illustrates the relationship of files to experiments.
35
-------
cc
LU
CO
z
""^
z
LU
_l
r (
u_
LU
Q.
-COOCM«^LOP~-CnOCM«TLOP^Cnr-i
COCOCOCO^t«*-«3-^J-<^-^J-LOLOLOLOLOLOLO
P^SfrHOOLOCMCnLOCOOP^'COO
COCOCOCO<3-«T«^-'!l-<5j-«}'LOLOLOLOLOLOLD
C\JCnLDCOOP~-^-rHCOLOCMCnLDCOOP«->OOO
COCOCOCO<3-«3-«3-«d-*J-«*'LOLf>LOLOLOlOLD
rHCOLOCMCnLDCOOP^^l-rHCOLOCMCnLDCO
fJ^j-LOOOCTlrHrOLOLDOOOrHrOtOLOOOO
COCOCOCOCO«j-«^-«}-«J-«*'LOLOLOLOLOLOLO
OP^^rrHOOLOCMCnLOCOOl ^j-rHCOLOCM
CO^LOOOCnrHCO^tLDOOOrHCOLOLDCOO
COCOCOCOCO«*'.^-rHCOLOCMCnLDCOO
cj«^-LOP^cniHco«3-LDOOcni ICO^J-LOOOO
COCOCOCOCO«J-^t-^^-^-«^-LOLOLOLOtOLO
P^'^iHCOLOCMCnLOCOOP^'i-r-ICOLOeMCn
CM«l-LDP«-CnrHCM^}-LDCOCnrHCO«^-LOCOcn
cococococo^-LDCOOP~
cM«^-LOP*-cniHCM^fLOP^cniicM>xj-LOCOcn
COCOCOCOCO«J-^J-'^-«^-^l-«J'LOLOLOLOLOLO
«d-i ICOLOCMCnLDCOOP^'^-r-IOOLOCMCnLD
CM^Lor^cnccM«!j-LOP^cnrHCM'*'LDr~-cn
cococococo^rj-^-^^-^j-LOLOLOLOLOm
rHCMCO^l-LOLOP>.COCnOrHCO^-LOLDP~-CO
OOOOOOOOOrHrHrHiHiHrHrHr-l
CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCNJCM
36
-------
3.2.2 Tape Data Records
The first record of the first block of each file contains alphabetic
ASCII characters of column headings for the data fields in the records that
follow. This first record may be considered to have a FORTRAN statement,
format (132A1). Column headings are coded in four or five characters, such
as AL2TD, where the first three characters are tower and level identifiers,
and the last one or two identify the meterological data. Here the "A"
identifies Tower A, "L2" identifies level 2 (10 m), and the "TD" specifies
the meteorological data, the value of temperature difference between the
10-m level and the 2-m level reference temperature. Table 7 lists the codes
for towers, levels, and meteorological measurements.
All data records following the first record have data fields arranged
as illustrated in Table 7.
TABLE 7. DATA RECORDS FORMAT
Position
1 to 4
6 to 8
10 to 11
13 to 16
17 to 24
25
26 to 132
Contents FORTRAN format
Year 1980
Day (Julian) 290 to 315
Hour (MST) 00 to 23
Seconds 0, 300, 600, 1500,
1800, 2100, 2400,
2700, 3000, 3300
Meteorological data
Data quality flag
Meteorological data
& data flag fields
14
13
12
14
F8.3
Al
F8.3, Al
37
-------
3.2.3 Data Record Types
Meteorological data acquired during each experiment was classified into
19 separate files, of 13 different record types, according to tower, level
and data type.
Table 8 classifies the record types of the first set of files, numbers 1 to
323, edited but uncorrected data.
TABLE 8. METEOROLOGICAL RECORD TYPES
No. of Meteorological
Type data fields data Data codes
Tower A
1 12 Temperature profile, 2 to 150 m T, TD, TC
2 9 Solar radiation, fast response NR, SR, TF,
temperature, cup-and-vane wind ST, MS, MD,
statistics SD
3 12 Wind data, 2-m level U, V, W, WS,
WD, IX, IY, IZ,
UX, VX, SP, DR
4 12 Wind data, 10-m level U, V, W, WS,
WD, IX, IY, IZ,
UX, VX, SP, DR
5 12 Wind data, 150-m level U, V, W, WS,
WD, IX, IY, IZ,
UX, VX, SP, DR
6 8 Wind -.'ata, 40-m level U, V, W, WS,
WD, IX, IY, IZ
7 8 Wind data, 80-m level U, V, W, WS,
WD, IX, IY, IZ
(continued)
38
-------
TABLE 8. (Continued)
No. of Meteorological
Type data fields data Data codes
Tower B
1 8 Wind data, 2-rn level U, V, W, WS,
WD, IX, IY, IZ
2 12 Wind data, 10-m level U, V, W, WS,
WO, IX, IY, IZ,
UX, VX, SP, DR
3 12 Wind data, 30-m level U, V, W, WS,
WD, IX, IY, IZ,
UX, VX, SP, DR
4 8 Temperature profile, T, TO, TC, ANEPH,
2 to 30 m BNEPH, CNEPH
Tower C, 0, E, F
1 11 Temperature profile, 2 to 10 m & U, V, W, WS,
Wind data, 2-m level WD, IX, IY, IZ,
T, TD, TC
2 12 Wind data, 10-m level U, V, W, WS,
WD, IX, IY, IZ,
UX, VX, SP, DR
For the second set of files, numbers 324 to 612, corrected wind data,
only 17 separate files were required for each experiment since no
temperature data were corrected, thus record type 1 from Tower A and record
type 4 for Tower B were omitted. Otherwise, the classification of files by
record type is the same as in the first set of files as shown in Table 8.
39
-------
Table 9 is a printout of the first five records from the first seven
filss-from the first set of meteorological data files, numbers 1 to 323. It
illustrates how the first alphabetic record of each file identifies with
headings of the meteorological data that follows in the data records in the
remainder of the file. The seven files shown represent the seven different
record types used to record data from Tower A.
Table 10 is a printout of the first five records from the first seven
files from the second set of meteorological data files, numbers 324 to 612.
Since temperature profile data were omitted from this set of files, only six
record types were used for Tower A, and the seventh file shown is wind data
from Tower B. If wind data are required from Cinder Cone Butte, only the
second set of files should be used if ERT's corrections are considered
satisfactory. The first set of files, edited but uncorrected data, are
available for application of other corrections and for comparative purposes.
Only the temperature profile data are considered reliable.
Further refinements of propeller anemometer data may/be possible when
results of more wind tunnel studies are analyzed and when comparisons are
made between propeller and sonic anemometer data from colocated systems
operated at Small Hill Impaction Study No. 2 at Hogback Ridge, New Mexico
during 1982.
40
-------
a
o
a
UJ
o
ce.
o
u
a
LU
a:
a.
LU
_J
a.
00
00
O
O
Q£
O
f-
-T
U 0
O NO
03 rH
_J 1
O
r».
in
<
in
co
L) O
eg
IT*
O O
F-* *
O 1
_J
<
fH
t?
in 0-
<
jpj
^
U i-l
^ o
-J iH
r-t
O
Q sO
H-
* 1
^ j
<
o
fH 10
*
II tO 0
-I rH
O <
UJ eg
w o m
U r~
eg o
<
tH
i-.
0 rH
eg i
_j
<
to
*
^m ,
rH O
J rH
<
0
U
UJ
10
3Z rH
>- 0
Q C"
~) CJ
O
Gj
K ->
3
o -o <£
o -o co
o o o
00-
o cj eg
o in NT
in rH rH
1 1 1
3
0 O O
o in in
i
o o o
O Is* (7*
O Is* 1^
3
o -o co
O I-l rH
O CO f-
m i i
1
3
O rH fH
O -O NO
0 l> e>
1
3
O fH fH
O rH iH
* 0 0
3
o m r-
o eg NT
o -I to
in i i
3
o eg eg
o in in
o to ro
o o o
ro fH rH
3
o to
r-t fH rH
0 O O
rj cj eg
o o o
cl . j CO
u 0 ON
o
Q 0
in
o
p o>
c
CJ rH
-J rH
< to
rH
V) 0
CJ NO
<
iH
to
(- 0
to
03
rH
O
U. rH
eo o
.j
<
O
eg IH
t- 0
n eg
. _j
O <
LU CO
M u. in
eg o
<
in
o
re eg
in
o
_j
<
to
rH
QT O
^ .
O
<
O
(J
UJ
in
E rH
>- O
O O'
-> cj
C
o
a: c
3
0 CO rH
O CJ in
3
O CO rH
O CO CO
0 rH rO
NO O
rH iH
ro to
3
O O iH
o o eg
o in NO
o in
3
O O NO
O O rH
0 O O
3
o eg eg
o eg eg
o o o
o o o
to
i
3
O O NO
O O >H
0 O O
3
o eg eg
o in in
o to to
o o o
i
3
o NO in
0 rH rH
3
0 -C IO
O rH tO
in o o
1 1 1
o o o
O 0 O
ro NO ON
rH rH rH
0 O O
c- c- o
CJ CJ CJ
o o o
>- O 'O
C ." O
ro
co
a: co
3*
-I 0
o
in
a -o
rH -T
_l
CO
X cf
rH eg
_i i
CO
o
X in
3to"
_j
<
CO
fH
M -T
r-t -T
_J
<
CO
tH
rH CO
<
O
in
x in
_J rH
**
CO
to
§m
H in
_J O
< to
^
ro r-
en in
n IH ro
O <
UJ CO
-j ro
W O
rH
<
fN.
fs.
O
rH eg
_i i
<
03
o
ON
^ .
rH CJ
_J
<
O
0
UJ
in
S. rH
V 0
Q 0
1 CJ
O
ei
Q: *>
3 U
o NO o
o in >o
o -o in
rH O
to ro
3
0 rH -T
o r-. r-
o f^ to
f
co r-
rH rH
3
o t>- eg
O v3- o
rH rH
3 u r
o < o
o m o
O CO O
eg
rH
to
3 r
0 rH 0
O O O
O t-- O
' ro
3
o co r-
o o eg
O rH rH
0
1
3
o r*. ro
o rn ro
o in 03
in eg rH
eg i i
I
3 s:
o to o
O rH O
o r- o
in eg in
CJ CJ
1 1
o o o
O 0 0
tO *O C7*
rH r4 rH
0 Cl O
eg cj cj
o o o
co cl co
& a- - Pi
eg s3-
J rH
in
to
x in
CJ rH
-J rH
to
O
D eg
cj r-
I O
< to
CO
* ro
in -a
~r ,
II CJ vT
O <
Ul rH
J CO
M to
eg
<
o
03
eg eg
_j i
<
o
ON
NO
3
eg to
_j
<
o
(J
LU
in
I rH
>- O
a o-
-) C-J
o
rj
a: c--
3 u
0 CJ CJ
O rH »0
in o
to ro
3 U
O O- 0
o in o
o to eo
NO m
3
O in rH
o ro
rH fH
3
o co eg
o eg r*.
o in rH
eg in
rH O
to to
3
o 03 in
o o in
O CO N3-
' -* -3-'
3
o ro rH
o eg to
o ro to
o
i
3
o in o
o eg in
in to eg
CJ 1 1
i
3
0 Nj- rH
O in o NO eg
CO O Is- CO
<
3
f o in eg
x in o co in
03 -o in in
_i
3 U
in o eg to
NO O NO NO
P rH O rH eg
CU rH rH r*»
-J rH rH O
< to to to
3 U
O O NT CO
W 0 O -O 1-4
in eg o rH N}
II CO NO NO NO
_J
o <
Z 3
UJ NO O O CO
_j 03 o N^ o
M rH O O rH
03 O
J rH
< 1
3
rH 0 CO -O
CO O CO O3
O 0 t^ - 0 O 0 O
O O C7- CN O
-> cj cj ej eg
o o o o
co a co co
rN- (J. (J- C.-N O
CO
in
N in
«3- 10
<
^
CO
>- ro
M
N± O
<
«
in
X rH
_l
O
0
Q eg
NT eg
_J rH
< to
fN»
NO NT
in o/
n -T in
o<
**
UJ CO
J O
IH fH
3-
<
in
to
o-
«T 10
i i
<
NO
CJ
to
^ .
NT NT
_l
<
o
u
UJ
in
X -H
> O
~> eg
o
CO
or o
3
o m o
0 f- 0
ro NT
3
o eo f»
0 |H fN.
0 ON NO
r» o
3
0 rH >O
o to co
o to eg
(N, e>
3 3
O O «T
0 0 NO
O O CJ
o
rH
10
3 3
O 0 rH
0 O >O
0 O O
m
3
O 03 «T
O NO NT
o o eg
o
rH
1
3 n
0 0 (O
O 0 CO
O O I*N.
in in ro
CJ CJ 1
1 1
3
o in eg
CD Kl fN»
O l> -T
in to NT
eg
1
o o o
000
ro NO (TN
rH rH rH
O O O
^- C*" CT*
cj cj eg
O 0 0
CO CO CO
a- v> f
3
o o o in
N O rH r-
M fN. O IT1 O
NO NT to ro
<
3
CO o o ro
NT O f^ CO
>- IT- 0 NO O
M
NO o in in
< ""*
3
ro o eo o
NT 0 in CM
X 03 0 r4 CO
_1
3 U
r- o NT -o
NT O tO CJ
P O» O Is- tJN
NO in NO eg
_J O O O
< 10 to to
3 O
& O rH ^*
N o o in o
(rt Kl O CJ -O
3 *
U *0 D >O
. _J
0<
3
UJ «O O N Ut
J CO O ^ O
M rH 0 r-t CJ
< i
3
l> o l> cj
to o to o
t> o r- in
43 ro in to ro
_j i eg i i
< i
3
in O rH NT
rH O O in
3 ....
NO m m in in
-i eg
< 1
o o o o
U 0 O 0
uj ro NO l>
in
X rH rH rH rH
y- o o o o
~> eg eg eg eg
o o o o
Q eo co co
c: o- » ON CTN
41
-------
<-c
<.
o
UJ
c_>
UJ
O£
O
O
o
Q
a:
Q-
O-
- o o o o
Q O O* O
> CJ CJ CM CJ
ry o o o
- c*> o o to
r£ CT C- ON O
3 U
ro o NO o
co o in NO
a: co o NO in
a ....
rH rH CM CO rH
-J O in O O
< ro ro to ro
3
to NT in in
a. -T cj in rH
t/i * . *
rH NT 1 NT «T
-1
3
ro CM o NO
NO rH r-- O
X ro cj co rH
rH CJ CJ CM
-II II
3
co o co in
o- ro co r^
x r* o in in
3 ....
IH ro i ro ro
_i
3
eo o ro CM
rH O O NO
NJ NT 0 FN, CJ
H NT 10 NT
<
3
CO 0 rH NO
NO O fN" CO
>- r» o NO o
rH CO CO f*-
-J rH rH rH
<
3
O 0 C» CM
in o NO NO
X in o NT o
rH O NT O
-J rH rH rH
3 CJ C
rH r-» ro r*
- ro o o o
CM -T ro ro
J rH iH f^
<
3
in o in o
ro o in o
X in o ro o
CM rH O fO
1 rH rH rH
in r>. ro CM
r*» CM P-» NO
a to r-. to r<-
cj co r-. CM NO
_l O KI rH 0
< ro CJ ro ro
3
Ift 0 CJ .T
NO cj o ro co
CJ LO O O rH NO
ii cj m NO in
i
o <
3
LU rH O IO rH
1 CO 0 CM CO
rH tO O IO IO
CM O
-J rH
< 1
3
co o ro CM
fN, O Kl O
NO 0 rH NT
CJ IO NT tO
II It
<
3
in o o NT
NT O Kt in
NO o in in
3 ....
CM NT NT NT
_J
<
O 0 O O
CJ 000
LU tO NO ON
cn
>- O O O O
: cj cj CM cj
O 0 0 0
o no cn co
;y o o c- o-
3 CJ
s3 0 CO NT
CO O O rH
o: NT o o NO
a
CO IO CO rH O
_J rH LT» rH rH
< to to to to
3 U
co NV to .r
0 rH 0 in
O. O CM CO rH
to ....
co co i r» co
_J
3 U
rH NT 0 - CO O
rH 0 CO O
X CO 0 CO rH
3 ....
to in i m NO
_i
3
CO 0 IO NO
f*- O rH CO
M rH 0 in NO
co NO in NO
<
3
0 0 in rH
CO 0 CM rH
>- in o NO cj
CO O rN. CO
1
<
3
& o in CM
NT o NO r~
x in o co in
co NO in in
^
Ro
NO IO
CM CJ CJ in
O ro r*. ro ro
co in rH in CM
J rH NT rH rH
< to CM ro ro
3 U
NO O O C^
fN, rH O r^ O
'CM CO O- 0 CO rH
II CO Is- Is- CO
O <
2 3
uj NO o o co
i co a NT o
M rH O 0 rH
CO O
-J rH
< 1
3 CJ
O O NO tO
CM O O CM
NO o m in
ND o ro in
in o in o
3 ....
co in m NO
.j
<
o o o o
CJ O O O
uj ro NO o
CO
>- 0 0 O O
"1 CJ CJ CM CJ
O 0 O O
CC '0 CO CO
CK C* o* o cr
3
co o in o
in o r^ o
N in o NO co
NT ro ro - ro o o NO
NT O f^ 17-
_i
<
3
CO O rH NO
in o ro co
X rH O tO CM
NT NO P* O
^
3 3
O CM CM rO
O rH fN> r** N0
NT CM f-- ps o
J rH IO K) rH
< ro CM cj to
3 3
O 0 0 NO
co in o o co
CJ CO NT O o NT
ii NT r» r-
> _j
o <
3
UJ CO 0 CO -T
1 O O NO NT
W i-l O O CM
NT O
-1 IH
< 1
3 3
O 0 0 NO
o o o r-
0 0 O CO
NT in O O O O
O O C" CT" O*
) CJ cj cj cj
0 O O O
CO T CO CO
py Q«\ Q. ^ CT-
3
o o o in
r- o rH r~
M r- o o NO
NO -T ro ro
i
<
3
CO o o ro
-T o r* co
>- o o NO o
NO o in in
-J rH
<
3
ro o co o
NT O Ul CM
X CO 0 rH CO
NO Nfl NO f»
_J
3 U
CM r* 0 rH
in cj in in
a to r^ o* o
NO CM CM CM O
1 rH NT rH rH
^ ro cj to ro
3 CJ
rH 0 -T i-<
O IO 0 CO rH
CJ CO rH O O NT
II NO CO f- CO
* _J
O <
z
3
uj NO o f*. in
J CO O NT O
M rH O rH CJ
NO O
-t rH
< 1
3 CJ
CO O O CM
tN. 0 -T rH
NT O NT -T
NO m in in
ii ii
<.
3 CJ
O* cj -T CO
o o NT ro
o o co NT
3 ....
NO NO m NO
<
O O 0 0
U O 0 0
LU ro NO O
CO
>- o o o o
a c- o a- t>
~> CJ CM CM CJ
O O G O
CO CO Cj co
CC O- (7- Cf- N>
3
NT 0 CM NT
CM a NT CJ
hsl NT 0 in NT
CJ CJ CM O
O
CO
3
o^ o r*. o
NT O CM O
>- in o in r-
CJ NO f** CO
0
ca
3
m o CM o
to o to ro
X in o NO r»
CJ rH rH rH
0 rH rH rH
CJ 3 U
O CO O- -T
o cj ro ro
U 3 0
in o co o-
o in o NO o
10 in NO o r~ NT
u CM co r» co
* o
O CQ
Z
3
UJ PN, Q fNt fNj
; NT o CM NT
M rH O rH rH
CJ O
0 rH
at i
030
Kl 0 NT CM
NT 0 0 NT
o a o in
^ a ...
CM NO NO NO
O 1 II
ca
O 3 O
co co o in
NO O CM CM
rH O O NT
0 ....
CM m NT m
0
CO
o o o o
O o o o
LU tO NO O
tn
>- o o o o
Q cr c?1 ^ o*
> Cd CJ CJ CM
O 0 0 O
CT CO CO co
EC O- O O 0-
42
-------
SECTION 4
TRACER GAS DATA
4.1 TRACER GAS RELEASE SYSTEM
Two tracer gases, SFg and Freon, were released at different heights
from the boom of a mobile crane. The mobility of the release system
resulted in a higher number of successful hours per test (normally six or
seven hours out of eight) in which significant tracer concentrations were
recorded on the hill. In only one experiment (212) were the wind patterns
so variable that it was not possible to align the release system upwind of
the hill.
The SFg and Freon tracer gases were stored in individual compressed gas
cylinders kept at ground level; flexible Tygon tubing, approximately 100 m
long, led from the gas cylinders to different release heights on the crane
boom. For the first nine experiments (201 to 209), the tracer release tube
was attached to the smoke generator platform at the smoke release height but
from 0.5 m to 1 m away, horizontally. For the last nine experiments (210 to
218), the tracer release tube was on a separate pulley system independent of
the smoke generator platform and about 1 m away, horizontally, from the
smoke release. The gas flow was monitored by separate rotameters on the SF/-
43
-------
and Freon cylinders, and each cylinder's weight loss was monitored by a
separate electronic digital scale.
Because of the difficulty in calibrating rotameters with 100 m of
tubing attached, the rotameters were used simply to monitor a constant
tracer flow rate; the weight loss of the cylinders (as recorded by the
digital scale) was used to determine the emission rate of each tracer. The
scales could be read accurately only to the nearest 0.05 Ib, and because the
SFfi flow rate was initially as low as 0.06 g/sec (0.5 Ib/h), the possible
uncertainty in the hourly emission rate determination could be up to 10%.
This problem was alleviated in the later experiments by increasing the SF..
b
flow rate to about 0.18 g/sec (1.5 Ib/h), thus reducing the emission rate
uncertainty to about 3%. Table 11 presents the average tracer release rates
in each experiment; release rates ranged from 0.06 g/sec to 0.20 g/sec for
SFg and from 0.86 g/sec to 0.98 g/sec for Freon.* Table 11 (Greene and
Heisler, 1982) also identifies locations of tracer release points. Release
point range (r) and azimuth (0) are determined with respect to the center of
CCB. , Elevation at release point (z) is determined with respect to CCB datum
(3100 ft. or S44.9-m MSL), and the heipht of release (Ht) is determined with
respect to local ground elevation.
4.2 TRACER GAS SAMPLING SYSTEM
Tracer sampling was accomplished by means of approximately 100
individual battery-operated samplers capable of either 10-min or 1-h
sequential operation. Each sampler contained 12 individual pumps, each of
*In analysis by electron-capture gas chromatography, Freon is about 20 times
less sensitive than SF_.
44
-------
= 3
£ =
OOO OOOOO OO^OO OOO OOOOO ^0 >£ i/
* -Si OO OOOO OOOOO
<
O
I
I OOO
10 O9>OO OOO OOOOO
OOOO OOOOO
LU
Qi
LU
o
OOO OOOOO
«n«r»«n OOOOO
OO
o
Ol
3
c
c
o
(J
OOOOO
CO
<
§0 c* 00*000
OO Ot^OOO OOO-
~ £ £ 2 2 ^ £ £ 3ooo
III I I I I I I I I I
-»O~ OO O «-« Or*
iflO^t r* BO ff1 *N O -T "^ -^>
*N «N C1 «N OOOO
I O O 9> O
000 OOOOO OOOO
i-t f- O O «"« O
«/< *n OOOO
o o S M ?3 o
- o o
oooo
000
OOOO OO
-------
v e
aa *j
*O -J
ss
1 1
o &
O <-i
o S
Or-COCOO O O OO » ^
O ^0 "^ W
in u^ t"> -a
oooo
lit!
^ o -o w*
-T O
09
7 -
1 1 1
cr» ^.
*M (*»
3S£
3 0 r^O
-T
^
O
1
CO
("1
0
r-
300
333
o o *-*
9 ft
3 a a
o o o
«^ M^ irt
a
r-* r- co cc .3- .» «o op eo oo eo o> «p
r-» <>* n O O O <
MOO O O i
O t/*O -J»I
ocooo
1 1 1 1 1
OOJ'O-?
oo
>^> NO »* IN
ooooo
1 1 1 1
<*-ioo\OOfOff*
oo o<->-ro
(~*>jt/^ioiOP^C'
ooooooo
-------
O O O O O
sO -^ O O O
sj ;N O (S C
*^ r* 03 CO O
0000
v e
XI -^
a =
Ch
o
o
f^
o
o
0 O
n O
M -J
o o
o o
0
s
oo «;
o
o
o
0
5 «- 0
CO
M
sO
o
l/>
r>
O
O
o
^ ^
0 0
00 93
0 0
o o
CO
U-l
(J
O
O O
0
-J
a
o
o
.*
C4
O
irt
O NO
O f^
CO CO
o o
w< o
O)
3
C
C
O
u
o o o o o o
o o o o
LLJ
C_5
e* a* <
a o
i o o
o o
OQ
<
u-i O O O
-T O O O
r* -J- i/* cc
C 000
S*0 O ct
*^ O O
>o ^ ^
o o o o
47
-------
which intermittently** filled a Tedlar bag over the time period of interest.
Thus, each sampler could take sequential 1-h samples over a 12~h period or
sequential 10-min samples over a 2-h period. Normally, 1-1 bags were used
for both hourly and 10-min samples. Except for samples taken from
reflection masts (described below), all samples were taken at 1 m above
ground level.
Figure 4 (Lavery et al., 1982) shows the locations of the 70 fixed
samplers and also the 10 movable samplers that were placed on either the NW
or SE side of the hill, depending on the prevailing wind direction. Of
these 80 samplers, typically 60 were used for 1-h average samples and 20
were used for 10-min average samples. Another 20 samplers were used for
reflection masts, for background ambient air samples, and for colocated
samplers. Table 12 presents the locations of samplers with regard to range
and azimuth and Cartesian coordinates using the center of CCB as origin.
Elevations of samplers are also presented with regard to a datum of 944.9-s-i
MSL.
The design of a reflection mast is shown in Figure 5 (Lavery et al.,
1982). Air samples were drawn in from 3 m and 6 m (in addition to the
normal 1-m height) and also at an uphill site equal in elevation to the 6-m
height. The purpose of this sampling strategy was to determine if tracer
concentrations would "reflect" off the surface as predicted by some
disperson models. Four of these reflection masts were used during Cases 203
to 218. Normally, the 3-m height was sampled on only one of the reflection
masts; the other masts were sampled at 6 m and 1 m, in addition to the
uphill site.
**For a 1-h average sample, a pump sampled intermittently for about 1 sec
every 15 sec.
48
-------
o o
g V)
II
,^.- ;- -^_ .y,.- , ,-- i
^ ,' \x i»vg^5^L '^^ ,. '^^-~/-' _ " ' / _ _
1.
t-
0)
s_
3
at
49
-------
c
en
«r-
in
cu
o
in
03
C
o
o
O)
0)
-------
ce
0.
to
CO CO O
^^ L^^ 00 0^ '^ CO G^ *-.P ^^ ^«O VJD uO tO ^t" i-O ^5* 00 ^P C3 f^^ 00 ^*^ t"^ (^^ ^^"
f^> ^^ 00 ^^ rH f^^ ^J* r^* OO O^ O"^ *-O C^ 00 ^P (-O CO r**"t CO uO O^ CO CO ^^
'S- CO CM CM rH iI rH rH rH i-HCMCMCMrHrH
I I I I I I I I I I I I I I I I I I I I I I I I
rHr-t romCMCMrH OJ
I I I I I I I I I I I I I I I I I I I
ooooooooooooooooooooooooo
pH CO ^^ r^ ^O ^^ f^ f^^ ^H ^^ QQ f^. p^ f^^ 00 ^^ f^^ cO 00 G^ C3 r~^ C^^ *^* ^O
^^ LU
E U
ro O
s~\
M E
I I
un-^-COOOCMCMCMrHrH rH LD
X E
COCMrHCM <*'rHrH'HrH
I
OOOOOOOOOOOOOOOOOOOOOOOOO
E
fOCMr-HCM rHcdforo<^tv5oodaDo6o6odo!Dooa3adodcc)o6i^
a
0)
3
C
C
o
U
51
-------
o
01
c
r- j
-P
C
o
u
UJ
_l
Q.
<
CO
CO
a:
OJ
o
<
cc
I-
csj
LU
CO
- S
a: e
O a)
"O
S-
01
"a. >
e u
(0 O
CO r
s-
0)
(DOTtHrHOOOOOOOO DOOOOOOOO
CM C\J r-i i-l rH i I rH
rH Csl i I -P
i- C
P O
to o
Ol -P
in o
Ol O
i r-i
O) «
C PO
TJ _1
C 1
TO 2
W> I
O! <7>
O5
C «*
03
-------
Background air was sampled during each experiment by at least one
sampler upwind of the tracer release point. At two locations during each
test, an extra sampler was placed next to the normal sampler and set to
sample at the same time. These colocated samplers were used to assess the
variability in the sampling technique.
The mechanical reliability of the samplers was relatively poor, with a
typical pump breakdown rate of about 20% during each test. During the
earlier experiments, the mechanical breakdown problems, when combined with
sampler crew mistakes in setting the sampler times, resulted in fairly
low data capture for some of the experiments. However, as the sampler crew
gained experience, the data capture during the later experiments was limited
mainly by mechanical problems.
The sampling system design is proved to be a good compromise between
total system flexibility and personnel endurance. For example, it was not
possible to operate many more 10-min samplers because bags had to be
manually changed by the sampler crew every two hours for each 10-min
sampler. The utility of the reflection mast system cannot be properly
assessed at present because a more detailed study of the results is
necessary.
4.3 TRACER ANALYSIS SYSTEM
The analysis of the bag samples was done by means of chromatographic
separation and electron capture detection at the NAWC gas analysis
laboratory in Boise. Seven gas chromatographs (GCs) were originally used in
53
-------
3
the laboratorythree Baseline Industries units, three S Inc. units, and
one AID (Analytical Instrument Development, Inc.) unit. The output of all
but the AID GC was evaluated by electronic integraters (with strip chart
backup) to give the area under tracer gas peaks. The AID's output was
recorded on a strip chart and evaluated by measuring peak height. The AID
differed from the other GCs in that it operated with several attenuation
factors ranging up to 64 x 10 . It was the only instrument whose molecular
sieve column could separate SFg, Freon, and oxygen successfully;
consequently all Freon analyses were done on the AID GC.
The GCs were numbered 1 through 8 for simplicity of identification,
although no instrument was designated number 4, due to the requirements of
the electronic integrators, each of which processed the output from several
instruments. The Baseline units were removed from service after the first
five experiments; two more S instruments replaced them after the twelfth
experiment as GCs numbers 1 and 2 because of drift problems. The detection
limit of the GCs was about 5 parts per trill on (ppt) for SFg and about 100
ppt for Freon.
A chromatograph showing a good separation of the tracer gases using a
5A molecular sieve column is illustrated in Figure 6 (Lavery et al. , 1982).
The SFfi and Freon separate before the large oxygen peak, with a total
analysis time of about 4 min per sample. The SF& areas were calculated by
an electronic digital integrator (the area under the peaks is directly
proportional to concentration). With six chromatographs and an average of 4
min per sample, a total of 90 samples per hour could be analyzed.
54
-------
T
CN
O
in
co
U
O
* T
. (U I
5-5
-------
For quality assurance, about 3% of the bag samples from each experiment
were analyzed again, usually on a chromatograph different from that used in
the first analysis; all Freon analysis, of course, had to be done on GC #8
(AID). These "recount" data have been used for estimates of the precision
of the analytical procedures. Most analyzed bags were then flushed two
times with nitrogen and returned to the field. The exceptions were bags
that contained high tracer concentrations (greater than 1 part per billion
SFg or greater than 10 part per billion Freon); these bags were discarded to
prevent possible contamination caused by tracer desorption from the bag
walls. Figure 7 (Lavery et al., 1982) illustrates the flow of procedures
followed in the bag sampling and analysis.
Calibrations were performed on each GC at the start and finish of each
analysis day. Nine calibration gases, ranging from about 10 ppt to 40 ppb
SF- and from about 200 ppt to 800 ppb Freon, were used to calibrate each GC
in the early experiments. The calibration points were reduced to seven (10
ppt to 10 ppb SFg) in later experiments because no SFg tracer concentration
greater than 10 ppb was ever detected in the field studies. A check with
one calibration gas (usually.100 ppt SFg) was performed every four hours on
every GC; information (date, time, and integrator area) was then written on
the same data sheet. Because of the large number (more than 14,000) of
tracer analyses performed w'th these data, the actual calculation of
concentrations, Figure 8 (Lavery et al. , 1982), from GC responses was done
by computer at ERT's office in Concord, MA. Experiment, bag, and sampler
numDers, sampler location code, and sampling time were entered into the
computer system from the sampler log sheets by means of a remote terminal in
NAWC's laboratory in Boise. GC calibration data (GC number, time and date
-------
Field
Samplers
1
GC Analysis
Recounts
Bag Flushing
High
Concentrations
Figure 7. Bag sampling and analysis procedures.
57
-------
Field
Sampler
Data Sheets
Recounts
Computer
i
Concentrations
GC Calibrations
Calibration
Data Entry
Figure 8. Procedures to obtain tracer gas concentrations.
58
-------
of calibration, calibration gas concentrations, and GC responses) were
similarly entered from calibration log sheets. After sample analysis, GC
number, bag number, experiment number, analysis time and date, and GC
responses were entered and merged with sampler data. The concentration in
the bag was then calculated from the GC response (peak height or peak area)
by means of curves fit to the calibration data.
In view of the huge number of tracer samples and the operation of the
gas chromatographs for 16 h per day, the tracer analysis system worked quite
well. All samples were analyzed within 48 h of sample collection. The main
deficiency was that only one chromatograph could analyze both SFg and Freon.
The only major instrument problem occurred during the early experiments when
it was difficult to obtain reproducible results from three of the
chromatographs. These chromatographs were subsequently replaced, and the
analysis proceeded smoothly. The preliminary recount statistics show good
reproducibility of the tracer analysis system.
4.4 TRACER GAS DATA TAPE FILES
Data are stored at the National Computer Center, Environmental Research
Center, Research Triangle Park, North Carolina on Sperry UNIVAC 1100/83
systems magnetic tape, nine track, odd parity, ASCII-quarter word mode,
density 6250 BPI, tape number 004700. Record length is 132 characters, and
the block size is 1320 words, or 40 records per block.
UNIVAC users may assign the tape, @ASG,T CCB.U9S//////Q,004700 using
UNIVAC ECL. Upon request, copies can be furnished and translated into
formats acceptable to any computer using nine-track tape drives.
59
-------
4.4.1 Tape File Index
There are 18 tape files, one for each experiment, numbered 613 to 630
following the corrected meteorological tower data on tape number 004700.
Table 13 shows how tape files are related to experiments and dates of
operation.
TABLE 13. TRACER GAS TAPE FILES
File no.
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
Exp. no.
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
Date 1980
10/16
10/17
10/20
10/21
10/23
10/24
10/25
10/27
10/28
10/30
10/31
11/02
11/04
11/05
11/06
11/09
11/10
11/12
Exp. hours
(MST)
1700 to 2300
17CO to 2300
0000 to 0800
0000 to 0800
0000 to 0800
0000 to 0800
0000 to 0800
1700 to 0100
1700 to 0100
1700 to 0100
0000 to 0800
1700 to 2400
0000 to 0800
0200 to 1000
0000 to 0700
0000 to 0700
0200 to 1000
0200 to 1000
Stability class
E
E
E to F
E
E
E
E to F
E to F
F
E to F
E to F
E to F
E to F
E to F
E to F
E to F
E to F
E
4.4.2 Tape File Records
Table 14 shows the tracer data formats on each data record within the
files. The first two records of the first block of each file contain
alphabetic characters for the column headings for the data records that
follow. Table 15 is a printout of the first block (40 records) of the first
tracer data file (Experiment 201). Tracer concentrations of SFg are
60
-------
TABLE 14. TRACER DATA FORMAT
Position
1 to
5 to
14 to
20 to
24 to
32 to
38 to
44
46 to
56 to
63
66 to
77 to
82 to
93 to
102 to
113
118 to
129
3
12
16
22
29
35
41
53
59
75
80
91
96
110
126
Contents
Experiment number
Sample collection date
Collection location
Sampler ID
Bag Number
Sampling start time
Sampling end time
Sample flag
Analysis date
Analysis time
Gas chromatograph
GC response
to SF6
GC attenuation
for SF6 analysis
GC response
to Freon
GC attenuation
for Freon analysis
SF6 concentrations
Data reduction flag
Freon concentration
Data reduction flag
Heading
EXP
SAMPLING DATE
SITE
SID
BAG
SAMPLING START
SAMPLING END
QF
ANALYSIS DATE
ANALYSIS TIME
GC #
SF6 RESPONSE
GC AT
FREON RESPONSE
GC AT
SF6 PPT
QF
FREON PPT
QF
Comments
mm/dd/yy
characters
characters
6 digits
hhmm (MST)
hhmm (MST)
G = valid sample
R = recount sample
Q = questionable
sample
B = bad sample
S = suspect sample
mm/dd/yy
hhmm (MST)
F10.2 FORTRAN
format
FORTRAN integer
zero if GC is
not AID
F10.2 FORTRAN
format
FORTRAN integer
F9.1 FORTRAN format
Always "E"
F9.1 FORTRAN format
Blank if no
Freon analysis;
"E" otherwise
Notes:
1. Headings occupy first two records of the file.
2. Records are sorted by sampling location, sample collection date, sampler ID
and sample flag, in that order.
61
-------
or u.
o
Z H O O O O O
O Q.
UJ Q.
CK
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOi
UJ UJ UJ UJ UJ Itl In til nt yj in QJ nj ii I i» i
n i in 11 1 ii i ii i
yj .it ii i ii t it I ii I til yj n t n I
u. a
in a. 1-1 r^ r-
co -o
o o ro (
in in 1-1
i-icoinc\iroinroeo
r*» rn in o f-i f-H
o r*. in HI
i-l -O <^^if^'«fi^*fOOOOOOOO
ujooooooooooooooooooor^oooooa-ooo-ooooooaoo
inooooooooooooooooooocjaoinroorHor-to-ooooooooo
o o a) i-t m ooinm cjo^>Hin>o>o i-iin -H -3- O* tO «3* 9^ CJ f4 pH F^> CO O &*
iH r^rttOrt coi-iioinr^Hi HCOHI-J'
tKuifi r-ii-t in o vi- cj in rt co
i-l rt
Ut*CJMiHlOrO«\li-IIOi-ICMr-IH1lOCMfHi-l«OCOcOCO<0HtOK1l-ICMr-ll
>
i oooooooooooooooooooooooooooooooooooooo
< CQCOcocOcOcOcOcOcocOcocOcocococOcOcOcQcOcocOcOcOcOcOcocOcOcOcocococococOcoco
oooooooooooooooooooooooooooooooooooooo
Qoooooooooooooooooooomoinutuiinininininoooooooo
_J rt IH 1-1 rl
a, »-
z:o:ooooooooooooooooooooinoinininintninininoooooooo
oooooooooooooooooooooooooooooooooooooo
Ul ^
O oooooooooooooooooooooooooooooooooooooo
MIUXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
coxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx>.xxxx
Ok OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
XOi^OOOCOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO
UJ CJ (J CJ CJ CJ CJ CJ CJ CJ CJ CJ CJ CJ CJ
-------
contained in record positions 102 to 110 (F9.1 FORTRAN format), and Freon
concentrations are located in positions 118 to 126 (F9.1 FORTRAN format). A
quality flag accompanies each concentration that must be "E" to validate the
value. Since no Freon tracer was released until Experiment 208, the "0"
values listed in the records in Table 15 do not have an "E" flag indicating
no measured value rather than a zero Freon concentration.
An overall evaluation of the sample quality is indicated in record
position 44 where an appearance of "G", "Q", "B" or "S" indicates quality.
An "R" in this flag position indicates the sample is one that has been
analyzed twice or "recounted".
Four reflection mast sampler systems were used during Experiments 203
to 218. Air samples were drawn in from 3 m and 6 m (in addition to the
normal 1-m height) and also from an uphill site equal in elevation to the
6-m height. Normally, the 3-m height was sampled on only one of the
reflection masts; the other masts had two samplers, 1-m and 6-m, in addition
to the uphill sampler. All reflection mast systems are recorded under one
collection location identification in tape record position 14 to 16; each
sampler in the mast system is recorded with a separate number in record
position 20 to 22. The uphill sampler is denoted "900", 1-m height "901",
3-m height "903", and 6-m height "906".
4.5 GAS CHROMATOGRAPH CALIBRATION DATA TAPE FILES
Calibration data observed on GC's during all experiments are stored on
eight files, 631 to 638, immediately following the tracer gas data tape
63
-------
files on the same tape reel and with the same tape block and record
specifications.
4.5.1 Tape File Index
There are eight tape files, one for each GC employed in tracer gas
analysis. Nine calibration gases were used to determine responses in each
GC in early experiments, but this number was reduced to seven or less in
later experiments when no SFg tracer concentration greater than 10 ppb was
ever detected. GC number 8, the AID instrument, has calibration data on two
files, 637 and 638; the first has calibrated responses for SFg and the
second for Freon. No GC was assigned number 4. Table 16 shows how tape
files 631 to 638 related to calibration data.
TABLE 16. GAS CHROMATOGRAPH CALIBRATION DATA TAPE FILES
File no.
631
632
633
634
635
636
637
638
GC
1
2
3
5
6
7
8
8
no.
(SF6)
(Freon)
64
-------
4.5.2 Tape File Records
Table 17 shows the formats for each data record. The first record of
each calibration procedure contains alphabetic characters to identify and
supplement the calibration data that follows. Table 18 is a printout of the
first block, 40 records, of the first file.
TABLE 17. GC CALIBRATION DATA FORMAT
Position
Contents
Comments
1 to 5
7 to 16
18 to 27
Attenuation
GC response
Calibration gas
concentration
15 FORTRAN format, always
"1" except for GC #8
where measured value is
presented
F10.2 FORTRAN format, area
under tracer gas peaks
on GC except GC #8 where
peak height is measured
F10.2 FORTRAN format
Notes:
1. The header record that precedes each calibration is in the format:
"GC" - GC number
"COL" - molecular sieve column number
"GAS" - 1=SF6, 2= Freon
"mm/dd/yy" - date of calibration
"hhmm" - hour and minute of calibration (MST)
2. The last record for each calibration contains a value of -1 in each field.
65
-------
to
i
o
z
o
r-t
CM
O
CO
o o o o
I O O O O
o o o
o o o
O O '
o o <
CM -J- iH O
r-< O O rl
o o o
in o o
tfj vT1 Kl
i-t K> O
r-l
N
O O ^H O CM
O O 1 *-4 rH
oooeoooor-.ooo.>ooooeocooooooeoo
OOOOOOOOi-tOOCTOOOOOOOi-IOOOOOOOO
'pHOOOOfHOCJ^i-tOOOOOOiHOCNJ^i-IOOOOO
lOPtlDOO I rHrHOOKlLIOOOO I r-lr-IOOK)LnoOO
Wooooooooootnooooooootnoooooooooo cioooooooo
o f~ O O
(J
1
(J
0
K^Ot^^TTLn^-CM *J K1OCJ-^C
fH CU \T \T\ CJ O C^J O CM \£\ I
r-t i-l U
1
M vO CJ
i-t
66
-------
SECTION 5
PILOT BALLOON WIND DATA
5.1 PILOT BALLOON WIND SYSTEM
North American Weather Consultants (NAWC) operated pilot balloon
systems from the more upwind of two locations about 1.3 km NW and SE of the
center of CCB (see Figure 2). Wind profiles are derived from double
theodolite measurements of trajectories from pilot balloons (pibals) or
minisonde balloons released approximately once an hour.
Two theodolites were positioned a known distance apart, and both
theodolites were aligned to true north. The positions of the theodolites
are determined by the azimuth and elevation angles of theodolite station 2
as observed by theodolite station 1. After the balloon is released, both
theodolite stations take simultaneous measurements of the balloon's
position, which is recorded as an azimuth and an elevation angle as observed
from each station. Thus, at each data point, two rays are defined. The
first ray is the ray Rl from Theodolite 1 to the balloon described by the
angles AZ1 and ELI. The second ray R2 is described by the angles AZ2 and
EL2 and is analogous to the ray Rl. Theoretically, these two rays will
intersect at the exact location of the balloon but experimental errors
generally cause the two rays to be skew and not intersect at all. Based on
67
-------
two such nom'ntersecting rays, the position of the balloon is analyzed as
follows.
It is necessary to find the line segment AB connecting a point A on Rl
with a point B on R2 such that the segment AB is the shortest possible line
segment with endpoints on the two rays. This constrains the segment AB to
be perpendicular to both Rl and R2. Now let Rl no longer be a ray of
infinite length, but be a vector originating at Theodolite 1 and ending at
the Point A. Similarly, R2 becomes a vector originating at Theodolite 2 and
ending at the Point B. Because the experimental error in determining the
balloon's position is an error in the measurement of an angle, the resulting
linear error in balloon position is directly proportional to the distance
from the point of observation to the balloon. Therefore, the point chosen
as the most probable location for the balloon is the Point C lying on the
line segment AB such that the ratio of the distances AC/BC is equal to the
ratio R1/R2. Thus, if the origin of the coordinate system is taken as
Theodolite 1 and the vector AB is taken as the vector originating at the
Point A (Point A is now synonymous with Rl) and ending at the Point B, the
balloon position C can be expressed by the following equation:
C = Rl + (--) A3
R1+R2
The quantity reported in the column labeled "error" in the heading records
of each tape file is actually the length of the line segment AB.
The values listed as "direction correction" in the heading records have
been applied to both the observed angles and the computed wind directions.
68
-------
The wind directions and speeds are midpoint averages. For example, the wind
at min 3.5 is taken from the midpoint of the balloon position at min 3 and
min 4. If min 4 was missing and it was necessary to use min 3 and min 4.5,
then the resultant wind would really be representative of min 3.75 instead
of min 3.5.
In situations where only one theodolite was operative or where one
theodolite lost track of the balloon, an assumed balloon ascent rate was
employed to determine vertical distance. If one theodolite lost track near
the end of the run, then calculations were based on a continued constant
vertical velocity. In all such cases, the assumed ascent rate is indicated
in the heading records and the requirement of a continued vertical velocity
noted by a special record at the end of wind profile data records.
b.2 PILOT BALLOON WIND DATA TAPE FILES
Data are stored at the National Computer Center, Environmental Research
Center, Research Triangle Park, North Carolina on Sperry UNIVAC 1100/83
systems magnetic tape, nine track, odd parity, ASCII-quarter word mode,
density 6250 BPI, tape number 004700. Record length is 132 characters, and
the block size is 1320 words, or 40 records per block.
UNIVAC users may assign the tape, @ASG,T CCB,U9S//////Q,004700 using
UNIVAC ECL. Upon request, copies can be furnished and translated into
formats acceptable to any computer using nine track tape drives.
69
-------
5.2.1 Tape File Index
*
There are 21 tape files, 9 containing wind profiles for 9 days
preceeding the days with tracer gas release, and 18 with data from
experiments 201 to 218. Files are numbered 639 to 665 following the gas
chromatograph calibration data on tape number 004700. Table 19 shows how
tape files are related to experiments and dates of operation.
TABLE 19. PILOT BALLOON WIND TAPE FILES
File no.
Exp. no.
Date 1980
Release location
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
102
103
104
105
106
107
108
109
110
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
9/17
9/19
9/20
9/22
9/23
9/24
9/25
9/26
9/27
10/16
10/17
10/20
10/21
10/23
10/24
10/25
10/27
10/28
10/30
10/31
11/02
11/04
11/05
11/06
11/09
11/10
11/12
NW
SE
SE
SE
SE
SE
SE
SE
SE
NW
NW
SE
SE
SE
SE
SE
NW
NW
SE
SE
NW
SE
SE
SE
NW
SE
SE
70
-------
5.2.2 Tape File Records
The ffrst seven records of each pilot balloon ascent are alphabetic
characters that identify and describe the observation. The next two records
are alphabetic characters of column headings for the data record that
follow. A blank record separates each balloon ascent during an experiment.
Table 20 illustrates the data formats on each data record within the files.
TABLE 20. PILOT BALLOON WIND DATA FORMAT
Position Contents
Comments
Heading (FORTRAN Format)
1 to 5 Time from release
10 to 15 Azimuth angle,
theodolite 1
19 to 24 Elevation angle,
theodolite 1
29 to 34 Azimuth angle,
theodolite 2
38 to 43 Elevation angle,
theodolite 2
46 to 51 Linear error of balloon
position
56 to 60 E-W coordinate of balloon
position
65 to 69 N-S coordinate of balloon
position
74 to 78 Ht. of balloon above
datum, 944.9 m
82 to 86 E-W component of wind
speed
90 to 94 N-S component of wind
speed
98 to 102 Vertical speed
106 to 110 Wind speed
115 to 117 Wind direction
TIME (MIN)
AZ.l (DEC)
EL.l (DEC)
A2.2 (DEC)
EL.2 (DEC)
ERROR (M)
X (M)
Y (M)
Z (M)
U (MRS)
V (MRS)
W (MPS)
SPEED (MPS)
DIRN (DEC)
F5.2, min.sec
F6.2, degrees.tenths
F6.2, degrees.tenths
F6.2, degrees.tenths
F6.2, degrees.tenths
F6.1, meters.tenths
15, meters
15, meters
15, meters
F5.1, meters per sec
F5.1, meters per sec
F5.1, meters per sec
F5.1, meters per sec
13, degrees
71
-------
Table 21 is a printout of the first block, 40 records, of the first
file. It illustrates a situation where only one theodolite was operative
and an assumed ascent rate was used.
72
-------An error occurred while trying to OCR this image.
-------
SECTION 6
TETHERSONDE AND MINISONDE DATA
6.1 TETHERSONDE AND MINISONDE DATA SYSTEMS
North American Weather Consultants (NAWC) operated tethersonde and
minisonde systems from the same two locations NW and SE of CCB from which
the pilot balloons were released. The tethersonde was operated in an
ascent-descent sequence yielding profiles of temperature, pressure, wind
speed and direction at intervals of once an hour (or more frequently) to
heights of at least 200 m above the local terrain. Release times were
scheduled 30 min after the pilot balloon soundings to obtain wind profiles
every half hour.
The principal quality control check performed by NAWC, other than
routine operational checks before each flight, was a comparison of the
output from the sonde package with data taken concurrently from the 150 m
tower. Three profiles taken when the sonde was flown from the NW sounding
site, about 1.5 kn SSW of Tower A, have been compared with tower data taken
during experiment 201. Averaged data for the period of 15 min at a constant
altitude near 80 m are also compared with tower data at that level. The
temperature data from the sonde appear to be about 1° C high, although the
profiles are still useful in defining gradients. The wind speed data
compared favorably with the now corrected tower data.
74
-------
Minisonde flights were conducted when wind speeds were too high to
allow tethersonde operation, or when the tethersonde system was not working.
The minisonde was operated as a free release balloon sounding to heights as
high as 500 mb or 5.5 km. Sometimes the balloon was tracked to serve as a
pilot balloon for wind profiles. Minisonde data consisted of profiles of
temperature, wet bulb temperature and pressure.
6.2 TETHERSONDE AND MINISONDE DATA TAPE FILES
Data are stored at the National Computer Center, Environmental Research
Center, Research Triangle Park, North Carolina on Sperry UNIVAC 1100/83
systems magnetic tape, nine track, odd parity, ASCII-quarter word mode,
density 6250 BPI, tape number 004700. Record length is 132 characters, and
the block size is 1320 words, or 40 records per block.
UNIVAC users may assign the tape, @ASG,T CCB.U9S//////Q,004700 using
UNIVAC ECL. Upon request, copies can be furnished and translated into
formats acceptable to any computer using nine track tape drives.
6.2.1 Tape File Index
There are 17 tape files, 10 containing tethersonde data, and 7 with
minisonde data. Files are numbered 666 to 675 for tethersonde data, 676 to
682 for minisonde data. Tables 22A and 22B show how tape files are related
to experiments and dates of operation.
75
-------
TABLE 22A. TETHERSONDE DATA TAPE FILES
File no.
666
667
668
669
670
671
672
673
674
675
Exp. no.
201
202
203
204
205
206
207
208
209
212
Date 1980
10/16
10/17
10/20
10/21
10/23
10/24
10/25
10/27
10/28
11/02
Release location
NW
NW
SE
SE
SE
SE
SE
NW
NW
NW
TABLE 22B.
MINISONDE DATA TAPE
FILES
File no.
676
677
678
679
630
681
682
Exp. no.
211
213
214
215
216
217
218
Date 1980
10/31
11/04
11/05
11/06
11/09
11/10
11/12
Release location
SE
SE
SE
SE
SE
SE
SE
Note: No valid tethersonde or minisonde data from Experiment 210, 10/30.
6-2.2 Tape File Records
The first six records of each tethersonde ascent-descent are alphabetic
characters that identify and describe the observation. The next two records
are alphabetic characters of column headings for the data records that
follow. A blank record separates each tethersonde ascent-descent during an
experiment. Table 23 contains the data formats on each data record within
the files.
76
-------
TABLE 23. TETHERSONDE DATA FORMAT
Position
1
11
19
26
33
40
47
54
62
70
to
to
to
to
to
to
to
to
to
to
7
15
23
30
37
43
51
58
66
73
Contents
Time of observation
Barometric pressure
Height of obs. AGL
Temperature
Rel. humidity
Mixing ratio
Wind direction
Wind speed
Potential temperature
Voltage
Heading
TIME
PRES
HT.
TEMP
(MIN)
. (MBS)
(M)
(O
RH (%)
M.R.
DIRN
SPD.
P.T.
. (DEC)
(MPS)
(K)
VOLTS
Comments
(FORTRAN Format)
F7.
F5.
F5.
F5.
F5.
F4.
F5.
F5.
F5.
F4.
4,
1,
1,
1,
1,
1,
1,
1,
1,
1,
HH.MMSS
millibars. tenths
meters. tenths
degrees Celsius. tenths
percent. tenths
ratio. tenths
degrees. tenths
meters per second.
tenths
degrees kelvin. tenths
volts. tenths
Table 24 is a printout of the first block, 40 records of the first
file. It illustrates a tethersonde ascent from the surface to 309 m, with
the first block of records containing the identification and heading records
and data records up to 151.2 m.
The first four records of each minisonde ascent are alphabetic
characters that identify and describe the observation. The next two records
are alphabetic characters of column headings for the data records that
follow. A blank record separates each minisonde ascent during an
experiment. This instrument was used when high winds prevented use of the
tethersonde, and in all experiments after 212, when the tetersonde system
did not operate. Data recorded were not subject to quality assurance
procedures as were the pilot balloon data or the tethersonde data. Table 25
contains the data formats on each data record within the files.
77
-------
Q.
1-lr-trHi-ii-lr-lr-lrHrH
U O> I
-» U. O O I
UJ (O IO i-l I
(Q 2 ^" O* O* ff* ^ ^* O* Iff* 9* O*
UJ 13 O UJ UJ
Z UJ 2: o n
O CO UJ Z Z __
bj uj o (ft
o: ao n.Q,K>tnr^ x:
§P f- OJ «^
M HI (Q
M H M U- £
~ <<-K1Z(3
_...."
h- - S X
uj M M .. in [!: oO>OOOOOOO
>-ZC3C3 cV5h~ r^i^f^Mi-ifHM
§>- 5 g a 10
c: 5 5 = uj
t o o o uj o:
HE
< O O
a. M co
CHS
-JON. tf> ir>
_/ M -.
-J O UJ
< n i- z: r^ ' ' '
UJCJ ^t^J.HH^^^p^f-I^^Mi
78
-------
TABLE 25. MINISONDE DATA FORMAT
Position
Contents
Headi ng
Comments (FORTRAN Format)
1 to 6 Time from release
Time (SEC)
9 to 13 Temperature TEMP. (C)
18 to 22 Wet bulb temperature WET BULB (C)
F6.1, MM.SEC(tenths);
8.9=8:54; 9.0=9:00;
9.2=9:12
F5.1, degrees Celsius.tenths
F5.1, all data unusable
29 to 33 Barometeric pressure PRESSURE (MBS) F5.1, mil 1ibars.tenths
Note: Altitude of observation must be developed from pilot balloon profile
that tracked minisonde ascent. Balloon ascent rate and time from
release will yield altitude.
Table 26 is a printout of the first block, 40 records, of the first
file. It illustrates the beginning of a minisonde ascent with observations
of temperature and pressure at 6- or 12- second intervals. The presence of
erroneous wet bulb data is also evident.
79
-------
OS
a.
UJ
_l
a.
o
CD
O
O
LU
UJ
K
5
UJ
O
o
to
uj
s
K
UJ
9
UD
CM
CQ
-» co uj
r:»- \ 3
U rO
_i O V
! _j o a. i
M r-( ^ U
I I I I I I I I I 1 I I I I I I I I I I I I F I I I I I I t I I I I I
80
-------
SECTION 7
EPA COMPLEX TERRAIN MODEL DEVELOPMENT
SHIS #1 MODELER'S DATA ARCHIVE - 1982
by
David G. Strimaitis*
Sonald C. DiCristofaro*
7.1 MODELERS' DATA ARCHIVE
The present modelers' data archive for the first Small Hill Impaction
Study (SHIS #1) of the CTMD program contains observed 1-h average tracer
concentration data, tracer release information, and meteorological variables
and derived parameters estimated at release height for each of the hours
during SHIS #1 in which either SFg or Freon tracer gas was released. The
method of estimating meteorological data appropriate to the release heights
of the tracer gases relies on a few central assumptions, and is objectively
applied to all the data with few exceptions. Those assumptions and a
description of the procedures are presented in the following sections.
Because the meteorological data contained in the archive are spatial
estimates of the meteorological conditions affecting the transport and
dispersion of the tracer gases, they should be viewed as approximate. At
*Environmental Research & Technology, 696 Virginia Road, Concord, MA 01742.
Section 7.1 to 7.5 and Section 7.7.
81
-------
best, the data could adequately represent conditions at the release point.
At worst, the data could be misleading. This is particularly evident in the
estimated wind direction. For example, the wind set at the 40-m level on
Tower A was partially operational during all of Experiments 205, 206, 210,
and 211. Some data from this instrument are recoverable, but winds are
estimated at this level only by making an assumption about the variability
in wind directions between 40 m and 80 m. At times, the assumption may not
represent the real situation. With the scale of CCB, a resulting wind
direction error of 10° could at times cause the modeled plume to miss the
hill entirely during a period when the actual plume produced significant
concentrations on the hill. Potential users of this archive should be aware
that subjective estimates of the most appropriate wind directions as derived
from evidence of actual plume transport directions are in preparation, and
will be included in a second version of this archive when available.
7.2 TRACER CONCENTRATION DATA
Methods of collection and analysis of SFg and Freon tracer gas data are
described by Lavery et al. (1982), and revised calibration procedures are
described by Strimaitis et al. (1983). Procedures and results of quality
assurance analyses of the tracer data are described by Greene and Heisler
(1982).
Tracer gas concentration data from the SHIS #1 data base available
through EPA include data from 1-h samplers, colocated 1-h samplers, 10-min
samplers, and reflection mast samplers (a subset of the 10-min samplers),
82
-------
and "recount" concentration data resulting from re-analyzing a subset of all
bag samples as part of the quality assurance program. The data base also
provides additional information on the sampler location and the quality of
each individual sample.
Tracer gas data contained in the modeler's data archive have been
assembled from these data. All values are reported as 1-h averages. These
averages include concentrations from all 1-h sampling bags labeled as good,
the average of good concentrations obtained from the two samples collected
at co-located sampler sites, the average of the two concentrations obtained
from samples included in the recount analyses, the average of good 10-min
concentrations at standard 10-min sampler sites, and the average of good
10-min concentrations obtained at the foot of the sampling mast locations.
In the case of the 10-min samples, hourly concentration averages are
included only if no more than one 10-min period was missing from the hour.
The position of each sampler is included in the archive along with the
1-h average tracer gas concentration. The coordinate system is a Cartesian
system with origin at the center of Cinder Cone Butte, x-axis oriented
toward the E, and y-axis oriented toward the N. The sampler position
identification code for each concentration is also included. A map of CCB
identifying each sampler position is presented in Figure 4.
7.3 TRACER RELEASE INFORMATION
The modeler's data archive contains the average emission rate of the
tracer gas source, the polar coordinates of the source position is the
33
-------
elevation of the base of the source crane, the height of the source above
the ground, and the times at which the tracer gas was turned on and turned
off.
As discussed in Greene and Heisler (1982), the emission rate is an
average mass release rate (g/s) from the time at which the release valve was
opened to the time at which it was shut off. In some cases, this period of
time was less than 1 h, but in most cases it was several hours. The start
and stop times for the release are referenced to the beginning and ending
time of each experiment hour, respectively. A start time of -10 (min)
indicates that the tracer was released 10 min before the start of the
sampling hour, and a time of -5 (min) indicates a release ended 5 min before
the end of the hour.
Coordinates of the source position are expressed in the hill coordinate
system, a polar grid centered on CCB. The zero height contour in this
system corresponds to the 3100-ft elevation MSL (944.9 rn). Release
elevations are presented in meters above the ground, and the elevation of
the ground at the release position is given as the difference in meters from
944.9 m MSL. A topographical map of CCB is shown in Figure I.
7.4 METEOROLOGICAL DATA
Meteorological data contained in the modelers' data archive differ from
those contained in the SHIS #1 data base in that all quantities apply to the
release height of the tracer gas rather than to the height of fixed
instrument levels, derived parameters computed from the meteorological data
84
-------
base are included, and 1-h averages are constructed. Background information
on the design of the SHIS #1 meteorological data system can be found in
Lavery et al. (1982), and information on the adjustments applied to the data
in preparing the refined data base can be found in Strimaitis et al. (1983)
and Greene and Heisler (1982).
A "spline under tension" method is used to interpolate meteorological
variables between instrument levels on Tower A (Cline, 1974). This method
produces a linear interpolation when a tension factor of 50 is specified,
and a cubic spline curve through the data when a tension factor of zero is
specified. The suggested nominal tension factor of 1.0 produces a smooth
curve through all data points in a profile without the cusps and regions of
high curvature between data points common with the cubic spline. After
inspecting a number of profiles produced with tension factors between 0.5
and 3.0, the factor 1.2 was selected. This factor produces slightly greater
curvature near instrument levels than does a factor of 1.0, but it also
reduces the magnitude of local maxima/minima between data points in regions
where the vertical gradient of the profile quantity must change sign.
Meteorology representative of release height is assumed to be
equivalent to data taken from the Tower A vertical profile at the same
height above the surface as the height of release even though the surface
elevation at the release point generally differs from the surface elevation
at Tower A. This approximation is consistent with the spatial resolution of
the meteorological instrumentation because the release locations lie between
1 and 2.5 km from Tower A, because differences in surface elevations between
the base of Tower A and release locations vary between -6.1 m and 1.5 m with
85
-------
a mean difference of -3.7 m, and because the vertical resolution of wind
measurements on Tower A is 30 m or greater above 10 m (wind sets are located
at 2, 10, 40, 80, and 150 m).
The 5-min sequence of data in the modeler's data archive is constructed
as follows:
Tower A wind speeds, wind directions, and temperatures contained
in the refined data base are scanned for missing data. If missing
5-min values are found, they are replaced with values estimated
using linear interpolation in time. Only UVW prop wind data are
used to develop wind information because the F460 cup-and-vane
instruments were placed at the 2-, 10-, and 150-m levels while
most release heights are between 20 m and 60 m.
The temperature, the vertical component of the wind speed, and the
horizontal wind speed and direction are estimated at release
height by "spline under tension" interpolation with a tension
factor of 1.2. Horizontal speeds and directions are first broken
into wind components, and the components are interpolated to
obtain the wind direction at release height. The speeds are
interpolated directly. The -"rO-m level wind data are incomplete in
four of the experiments (205, 206, 210, 211) leaving the
undesirable prospect of having to interpolate the profile between
10 m and 80 m. However, the u-component (E) of the wind data from
40 m is available, and the mean winds were approximately SE during
86
-------
each of these experiments rather than nearly N or S. Assuming
that the directional wind shear was small (between 40 m and 80 m)
at Tower A during these four experiments, estimates of the wind
speed at 40 m are made with the 40-m u-component and the 80-m wind
direction. The resulting wind speed profiles indicate that the
speeds estimated at 40 m are generally reasonable, and so these
speed and direction estimates of 40-m level winds are used in the
interpolation procedure for nearly all 5-min periods when the 40 m
v-component was unavailable. Because the wind directions near 40
m appear to differ substantially from those at 80 m during the
first hour and twenty minutes and the last hour of Experiment 211,
the spline interpolation is used for these periods of time.
Turbulence data are estimated at release height by employing a
linear interpolation rather than the spline interpolation. The
turbulence velocity scales a , a , and a are obtained from the
turbulence intensity values contained in the refined data base by
multiplying by the wind speed. Unlike wind and temperature data,
missing turbulence data are not filled in by interpolating in
time. However, estimates of a are prepared for those experiments
W
in which one of the horizontal wind speed props malfunctioned at
the 40-m level. In these instances (205, 206, 210, 211) reported
values of I in the refined data base are flagged as bad because
the reported wind speed is incorrect, but because the w-prop was
working, the a values are not necessarily deficient. a is
w vv
S7
-------
recovered by multiplying the original value of I by the original
value of the wind speed. A prop response correction derived from
the work of Horst (1973) is applied. More information on its
application to CCB data can be found in Strimaitis et al. (1983).
No prop response corrections are applied to a or a .
The Brunt-Vaisala frequency, N, is estimated at source height by
interpolating the temperature profile in the immediate vicinity of
the release height to obtain the local temperature gradient.
The critical streamline height, H , is obtained from the splined
C*
profiles of temperature and wind speed by means of the integral
formula presented in Lavery et al. (1982). A bulk Hill Froude
number is calculated for the layer between H and the top of the
tower, 150 m, and also for the layer between 2 m and 150 m. The
hill height in both calculations is the difference between 95 m
and the height of the bottom of the layer.
The Turner dispersion stability class is calculated from net
radiation and wind speed data by means of the method of Williamson
and Krenmayer (1980). Wind speeds measured by the cups at the
10-m level on Tower A (reported as scalar averages) and the net
radiation data are interpolated in time whenever missing values
are encountered. The stability class is calculated as a number
between 1 and 6, where 1 denotes stability class A. Both the
88
-------
stability class and the 10-m wind speed are included in the data
archive.
Most 1-h average data in the modelers' data archive are obtained from
this sequence of 5-min average data interpolated to release height. Only
stability class data are not obtained in this way. The 1-h stability class
is found from the 1-h average net radiation and 10-m wind speed. In a
second method for computing the hourly stability class, the 10-m 1-h average
wind speed is combined with cloud information according to the Turner
objective method. Because cloud observations were not recorded as part of
the SHIS #1 data base, we are providing the 1-h average 10-m wind speeds as
part of the modelers' archive so that other users of the archive may obtain
the cloud cover data from Mountain Home AFB or Boise, Idaho and determine
the stability class.
The remaining 1-h average data in the modeler's data archive are
constructed as follows:
The wind speed and direction are calculated as both vector and
scalar averages. Two versions of the scalar wind direction are
calculated. One is a scalar average of the 5 min vector resultant
wind directions. The second is a vector average of unit vectors
along each 5-min vector resultant wind direction so that all
directions have equal weighting, as in the scalar average, but the
averaging is performed with vector arithmetic.
89
-------
Temperatures and parameters calculated from the splined profiles
of the 5-min temperature and wind data (N, H , Fr) are simply
averaged to provide 1-h average values.
Horizontal turbulence data, a and a , are computed for 1-h
periods by adding the contribution due to the 5-min turbulence
values to the contribution due to the variability of the 5-min
average winds. Let (a),.,, ,. denote the total 1-h value of the
ou, u
standard deviation, let (a)5 Q denote the total 5-min standard
deviation, and let (d)60 5 denote the standard deviation of the
5-min average winds over a 1-h period. Then:
Note that although prop response corrections are not applied to
the (a)g Q values, they are implicity contained in Men c because
the 5-min wind data in the refined data base include corrections
for prop response and wake effects. Also, because no time inter-
polation is performed on the 5-min turbulence data, less than a
full set of 12 values may be available during some hours. In
2
these cases the average of the (o")r n values will be incomplete.
o, u
The number of (a),- 0 values contained in each hour are denoted by
N(su), N(sv), and N(sw) in the archive.
The vertical turbulence a is computed for 1-h averaging periods
W
in the same way as a and a , but the prop response correction
90
-------
suggested by Horst (1973) is also applied. Neither the (a ),- n or
W 3 U
the (a )cn r. data used in the formula for (a )cri n have the cor-
w bU,o w bU,u
rection already applied. Rather, the correction is applied
directly to (a ),-0 Q. Also, the construction of (
-------
metaorology file. Concentrations are given in ppt, followed by a data
quality flag (they should all be G) and a column indicating whether the
concentration is from a 1-h sample (0), or made up of five or six individual
10-min samples. The Cartesian coordinates of the sampler position follow
(*, y, z), expressed in meters.
7.6 MODELERS' DATA TAPE FILES
Data are stored at the National Computer Center, Environmental Research
Center, Research Triangle Park, North Carolina on Sperry UNIVAC 1100/83
systems magentic tape, nine track, odd parity, ASCII-quarter word mode,
density 6250 BPI, tape number 002689. Record length is 132 characters, and
the block size is 1320 words, or 40 records per block.
UNIVAC users may assign the tape, @ASG,T CCBTR.U9S//////Q,002689
using UNIVAC ECL. Upon request, copies can be furnished and translated into
formats acceptable to any computer using nine track tape drives.
7.6.1 Tape File Index
There are four tape files; for each tracer gas, there is one file of
combined meteorological data and tracer release data and one file of tracer
gas concentrations. Table 27 shows how the tape files are an anged on the
tape.
92
-------
TABLE 27. MODELERS' DATA TAPE FILES
File no.
1
2
3
4
Experiment no.
201 to 218
201 to 218
208 to 218
208 to 218
Comments
SFg meteorological and tracer
release data
SFg tracer concentration
Freon meteorological and
release data
Freon tracer concentrati
data
tracer
on data
7.6.2 Tape File Records - Meteorological and Tracer Release Data
The first record of the file is composed of alphabetic characters that
identify the experiment, hour and tracer. The second record has alphabetic
characters that reveal the tracer release data with regard to emission rate,
position of the mobile crane, height of emission release, and start and stop
time of the release. The third record is a blank, and the fourth has
alphabetic characters that are column headings for 12 records of
meteorological data that follow. A blank record follows the meteorological
records, and eight alphabetic records follow, containing hourly averages of
meteorological data. Table 28 illustrates the data formats of the
meteorological records.
Table 29 is a printout of the first block, 40 records, of file no. 1,
the first file of combined meteorological and tracer release data related to
SFfi gas. The first block contains data from the first hour of the first
experiment, 201, as well as data from part of the second hour. File no. 3
93
-------
TABLE 28. MODELERS' METEOROLOGICAL RECORD DATA FORMAT - TOWER A
Position
1 to
5 to
8 to
20 to
27 to
37 to
44 to
53 to
63 to
76 to
89 to
3
6
11
15
23
30
33
40
48
50
57
59
70
72
83
85
96
Contents Heading
Day (Julian)
Hour
Seconds
Stability class
10-m wind speed, C&V
Average temp, (splined)
Number of data points
in Tower A profile
Average wind speed
(splined)
Average wind direction
(splined)
Number of data points
in Tower A profile
Average vertical wind
(splined)
Number of data points
in vertical wind
Tower A profile
Turbulence velocity
u-component (along-
wind)
Number of data points
interpolated in
Tower A profile
Turbulence velocity
v-component (cross-
wind)
Number of data points
interpolated in
Tower A profile
Turbulence velocity
w-component
(vertical)
DAY
HR
SEC
SC
MS(10)
TEMP
#
WS
WD
#
W
#
Sigma-U
#
Sigma-V
#
Sigma-W
Comments (FORTRAN Format)
13
12
14
11
F4.1,
F4.1,
11
F4.1,
F5.1,
11
F5.3,
11
F8.3,
11
F8.3,
11
F8.3,
meters per second
degrees eel si us
meters per second
degrees
meters per second
meters per second
meters per second
meters per second
(continued)
94
-------
TABLE 28. (Continued)
Position
Contents
Heading
Comments (FORTRAN Format)
98
102 to 105
110 to 114
Number of data points
interpolated in
Tower A profile
Height of critical
streamline
Hill Froude number
118 to 122 Froude number
126 to 130 Brunt-Viasala
frequency
II
HC F4.1, meters
FR(HC) F5.1, calculated for the
layer HC to 150 m
FR F5.1, calculated for the
layer 2 m to 150 m
N F5.4
is similar to file no.l except the tracer release data are associated with
another gas, Freon, and meteorological parameters are estimated for
different release heights.
7.6.3 Tape File Records - Tracer Concentration Data
The first record of the file is composed of alphabetic characters that
identify the experiment, hour, and tracer concentrations. Data records
follow the first record and continue until another alphabetic identification
record is encountered to indicate another hour begins. Table 30 shows the
data formats for the tracer concentration records.
95
-------An error occurred while trying to OCR this image.
-------
TABLE 30. MODELERS' TRACER CONCENTRATION RECORD DATA FORMAT
Position
1 to 3
4 to 6
7 to 10
12 to 14
16 to 17
20 to 23
27 to 34
37
43
46 to 53
56 to 63
67 to 73
Contents Comments (FORTRAN
Sampler ID
Tracer ID
Year 1980
Day (Julian)
Hour
Second-ending
Concentration, PPT
Quality flag (G)
Number of samples in 1-h aver.
X-coordinate (E-W)
of sampler relative to
center of CCB
Y-coordinate (N-S)
of sampler relative to
center of CCB
Z-Height of sampler
above datum, 944.9 m
A3
A3
A4
13
12
14
F8.3
Al
11
F8.3,
F8.3,
F7.3,
Format)
meters
meters
meters
Table 31 is a printout of the first block, 40 records, of file no. 2,
the first file of tracer concentration data. The first block contains
tracer concentrations, SFg, as 1-h averages collected as 1-h samples (0), or
averaged from 5 or 6 individual 10-min samples. File no. 4 is similar to
file no. 2 except the tracer gas concentrations are of Freon and the tracer
release was started at Experiment 208.
97
-------
o
I
II
0=
0.
o
o
z
LU
O
o
o
a:
LU
o
K
a
o
ro
ca
i-tcooinoooo-oooLntoinooo-Too i-toainoor-iiooirioior^ooo^o
1-4 IO in to co v> tn t") tD o to to to to to to vt in vi tn v? en in to irt co or
-------
7.7 CONCLUSION
The present contents of the modelers' data archive for SHIS #1 cover
all experiment hours in which either SFC or Freon tracer concentrations are
o
quantifiable at locations on Cinder Cone Butte, and in which these gases are
released for a significant part of an hour. Meteorological data are
estimated at tracer gas release height by means of data measured at Tower A,
the 150-m tower erected approximately 2 km north of the hill. The method of
estimation is applied objectively, and relies on few assumptions. Tracer
gas concentrations for all available samplers are reported as 1-h averages.
Users of this archive should pay particular attention to the starting
and stopping times included with the tracer gas release information.
Observed tracer gas concentrations are not equivalent to modeled 1-h tracer
gas concentrations if the release began well after the start of the hour, or
if the release terminated well before the end of the hour. However, in some
cases adjustments to the modeled 1-h tracer gas concentrations could be
designed to take account of the actual release period, and the travel time
from the source to the sampler array. In any case, it is the user's
responsibility to screen out those periods which are inappropriate for
driving his model.
The user should also consult the modeling work presented in the first
and recond CTMD milestone reports. These will give the user some idea of
the representativeness of the interpolated wind directions contained in the
archive for 45 of the experiment hours. Note, however, that the other
meteorological data presented in those reports were not derived explicitly
from the data contained in this archive. Present and future CTMD modeling
will make use of the archive data.
99
-------
Finally, two hours of the original 45 hours contained in the CTMD first
milestone report have since been considered inappropriate for Gaussian model
development. Experiment-hour 205-5 displays an SF, concentration pattern
which appears to be inconsistent with the release height when compared to
the pattern during the previous hour. Experiment-hour 209-7 contains a
slowly propagating abrupt wind shift which clearly passes Tower A well
before it passes the tracer gas release location, so Tower A data are not
representative of the meteorological conditions at the release location.
100
-------
REFERENCES
1. Cline, A.K. Scalar- and planar-valued curve fitting using splines under
tension. Comm. of Association for Computing Machinery. 17: 4, 218-220
1974.
2. Greene, 8.R. and Heisler, S. EPA Complex terrain model development:
Quality assurance project report for small hill impaction study no. 1.
ERT Document no. P-B348-350, Environmental Research and Technology,
Inc., Concord, MA, 1982. 72 pp.
3. Horst, T.W. Corrections for response errors in ,- th1"-*-1?-component
propeller anemometer. J. Appl. Meteorol. 12: 716-725, 1973.
4. Hovind, E.I. , Edelstein, M.W., and Sutherland, V.C. Workshop on
atmospheric dispersion models in complex terrain. EPA-600/9-79-041.
U.S. Environmental Protection Agency, Research Triangle Park, NC. 1979.
5. Lavery, T.F., Bass, A., Strimaitis, D.G., Venkatram, A., Greene, B.R.
Drivas, P. J. , and Egan, B.A. EPA Complex terrain model development:
first milestone report-1981. EPA-600/3-82-036. U.S. Environmental
Protection Agency, Research Triangle Park, NC, 1982. 304 pp.
101
-------
REFERENCES (Continued)
6. Strimaitis, D.G,, Venkatram, A., Greene, B.R., Hanna, S. , Heisler,
Lavery, T.F., Bass, A. and Egan, B.A. EPA Complex terrain model
development: Second milestone report-1982. EPA-600/3-83-015. U.S.
Environmental Protection Agency, Research Triangle Park, NC, 1982. 375
PP-
7. Williamson, H.J. and Krenmayer, R.R. Analysis of the relationship
between turner's stability classifications and wind speed and direct
measurements of net radiation. Paper presented at the 2nd Joint
Conference on Applications of Air Pollution Meteorology, Am. Meteorol.
Soc., Boston, MA, March 24-27, 1980.
102
-------
TECHNICAL REPORT DATA
(Please read Instructions on the reverse before completing)
1. REPORT NO.
2.
3. RECIPIENT'S ACCESSION NO.
4. TITUE AND SUBTITLE
I I I UK f+IVLJ ^UD 1 I I U.&
EPA COMPLEX TERRAIN MODEL DEVELOPMENT. Description
of a Computer Data Base from Small Hill Impaction
Study No. 1 Cinder Cone Butte, Idaho
5. REPORT DATE
6. PERFORMING ORGANIZATION CODE
7. AUTHOR(S)
Lawrence E. Truppi and George C. Holzworth
8. PERFORMING ORGANIZATION REPORT NO.
9. PERFORMING ORGANIZATION NAME AND ADDRESS
10. PROGRAM ELEMENT NO.
CDTA1D/09-3062 (FY-84)
Same as 12.
11. CONTRACT/GRANT NO.
12. SPONSORING AGENCY NAME AND ADDRESS
Environmental Sciences Research Laboratory - RTP, NC
Office of Research and Development
U.S. Environmental Protection Agency
Research Triangle Park, North Carolina 27711
13. TYPE OF REPORT AND PERIOD COVERED
Final
14. SPONSORING AGENCY CODE
EPA/600/09
15. SUPPLEMENTARY NOTES
ABSTRACT
As part of the U.S. Environmental Protection Agency's effort to develop and
demonstrate a reliable model of atmospheric dispersion for pollutant emissions in
irregular mountainous terrain, the Complex Terrain Model Development Program was
initiated. The first phase, a comprehensive tracer field study, was carried out on
Cinder Cone Butte, Idaho, during the autumn of 1980. Eighteen quantitative tracer
experiments were conducted, each lasting 8 hr at night or early morning.' The main
tracer gas was sulfur hexafluoride; a second tracer, Freon 13B1 was used in ten of
the eighteen experiments. Averaged meteorological data were recorded from six towers
near and on the slopes of the hill. Data consisted of direct and derived measures of
temperature, wind, turbulence, solar and net radiation, and nephelometer coefficient
of scattering. Hourly wind profiles were obtained from pilot balloon observations;
tethersonde observations recorded profiles of wind and temperature.
Tracer gas concentrations were detected by a network of approximately 100 sam-
plers located on the slopes of the hill. The system used to collect the data, the
operation procedures used to run the system, and its performance record are described.
Tables of tracer gas release data have been included to assist in any modeling effort.
All meteorological and tracer concentration data have been edited and recorded on
magnetic tape and are now available, upon request, at the National Computer Center,
R.T.P., NC, either as copies or by interactive computer access.
17.
KEY WORDS AND DOCUMENT ANALYSIS
DESCRIPTORS
b.lDENTIFIERS/OPEN ENDED TERMS C. COSATI Field/Group
18. DISTRIBUTION STATEMENT
RELEASE TO PUBLIC
19. SECURITY CLASS (ThisReport/
UNCLASSIFIED
21. NO. OF PAGES
20. SECURITY CLASS (This page)
UNCLASSIFIED
22. PRICE
EPA Form 2220-1 (9-73)
------- |