EPA-450/3-74-085
EMISSIONS INVENTORY
OF AGRICULTURAL TILLING,
UNPAVED ROADS
AND AIRSTRIPS,
AND CONSTRUCTION SITES
by
Chattcn C. Cowherd, Jr. ,
Christine M. Guenther, and Dennis D. Wallace
Midwest Research Institute
425 Volker Boulevard
Kansas City, Missouri 64110
Contract No. 68-02-1437
EPA Project Officer: Charles O. Mann
Prepared for
ENVIRONMENTAL PROTECTION AGENCY
Office of Air and Waste Management
Office of Air Quality Planning and Standards
Research Triangle Park, N. C. 27711
November 1974
-------
This report is issued by the Environmental Protection Agency to
report technical data of interest to a limited number of readers.
Copies are available free of charge to Federal employees, current
contractors and grantees, and nonprofit organizations - as supplies
permit - from the Air Pollution Technical Information Center, Environmental
Protection Agency, Research Triangle Park, North Carolina 27711;
or, for a. fee, from the National Technical Information Service, 5285
Port Royal Road, Springfield, Virginia 22161.
This report was furnished to the Environmental Protection Agency
by Midwest Research Institute, Kansas City, Missouri 64110, in fulfillment
of Contract No. 68-02-1437. The contents of this report are reproduced
herein as received from Midwest Research Institute. The opinions,
findings, and conclusions expressed are those of the author and not
necessarily those of the Environmental Protection Agency, Mention
of company or product names is not to be considered as an endorsement
by the Environmental Protection Agency.
Publication No. EPA-450/3-74-085
11
-------
ACKNOWLEDGEMENTS
This report was prepared for the Environmental Protection Agency's
Office of Air Quality Planning and Standards under EPA Contract
No. 68-02-1437. Mr. Charles 0. Mann served as EPA Project Officer.
The program was conducted in MRl's Physical Sciences Division
under the supervision of Dr. Larry J. Shannon, Head, Environmental
Systems Section. Dr. Chatten Cowherd, Jr., was the Principal Investi-
gator for MRI. Dr. Cowherd was assisted by Ms. Christine Guenther and
Mr. Dennis Wallace.
Approved for:
i
TIJTE
H. M. Hubbard, Director
Physical Sciences Division
8 January 1975
ill
-------
CONTENTS
Page
List of Figures vii
List of Tables viii
Introduction. ........... ... 1
Unpaved Road Emissions 5
Source Extent 5
Correction Factors ... .... 7
Silt Content 7
Dry Days Per Year 9
Emission Factor. 9
Heavy Construction Emissions. .. 11
Source Extent. ...... 11
Emission Factor. .. ......... 11
Agricultural Tilling Emissions. .. 17
Source Extent. 17
Correction Factors 17
Silt Content 17
Precipitation-Evaporation Index ........... 22
Emission Factor. 24
Unpaved Airstrip Emissions 27
Source Extent 27
Correction Factors 29
Emission Factor. ............... 29
-------
CONTENTS (Concluded)
Page
Data Tabulations and Accuracies ..... ..... 31
References. .............. 39
VI
-------
LIST OF FIGURES
No. Page
1 Annual vehicle miles on unpaved roads 6
2 Unpaved roads silt content (percent) 8
3 Emission factors for unpaved roads (Ib/vehicle mile). . . 10
4 Annual acres of construction 12
5 Emission factors for construction (tons/acre) 14
6 Annual acres of land tilled 18
7 Agricultural silt content (percent) 20
8 Precipitation-Evaporation Index 23
9 Emission factors for agricultural tilling (Ib/acre) ... 25
10 Annual LTO cycles on dirt airstrips 28
11 Emission factors for dirt airstrips (Ib/LTO cycle). ... 30
-------
LIST OF TABLES
No. Page
1 Fugitive dust source parameters 2
2 Construction dollars-to-acres conversion factors 13
3 Estimated annual tillings by crop 19
4 Coded NEDS area source data for Alabama 32
5 Correction factors and corrected emission factors for
Alabama 34
6 Smallest geographical areas assigned single values. ... 36
7 Estimated error ranges for tabulated data 37
Vlll
-------
INTRODUCTION
Area dust emission sources contribute substantially to the
atmospheric particulate burden in many parts of the country. The
Environmental Protection Agency has identified several fugitive source
categories for inclusion in the National Emissions Data System (NEDS)
area source file. Among them are the following categories of fugitive
dust sources:
1. Unpaved Roads (dirt and gravel),
2. Heavy Construction Sites (such as road and building
construction),
3. Agricultural Land Tilling,
4. Unpaved Airstrips.
To determine the impact of these sources, it is necessary to
develop a national emissions inventory of these sources on a county-by-
county basis. Calculation of county emission totals for each source
category requires, in addition to the basic emission factor, (1) an
appropriate measure of the extent of the source type within the county
and (2) correction factors which adjust the emissions estimates for
local climatic conditions and properties of the emitting surface.
The basic emission factors with associated correction terms were
developed by MRI^/ under EPA Contract No. 68-02-0619. Table 1 lists
the measures of source extent and the correction parameters which are
required for the calculation of corrected emission estimates.
The objective of the program reported herein was (1) the calcula-
tion of source extent and emission factor correction terms on a county-
by-county basis for the source categories designated above, and
(2) the documentation of the methodology used in these calculations,
including procedures used to estimate missing data.
-------
I/I
Pi
w
H
H
S
T
PM
<4
PM
w
u
Drf
O
>
H
P
W
£>
M
H
M
O
;=>
pn
1 1
Ol
r-l
,0
cfl
H
CO
S-i
Ol
4~)
0)
B
td
j_j
CO
O-i
c
0
H
4-1
O
0)
J_)
}_l
o
u
/--*.
4J
d
01
4J
C
o
o
4-1
r-l
l-l
CO
\_s
0)
^
£j
4-1
X
01
4-J
Ol
U
Cfl
SJ
3
CO
4J
d
o>
4-)
X
0)
M-l
O
0)
S-i
3
co
cfl
01
^
£»-
S-i
o
M
01
4-1
cfl
U
01
o
J-l
3
o
CO
TD
0)
r-l
Ol
^
cfl
S-i
4-1
CO
0)
T 1
l-l
6
1
0)
T 1
0
l-l
fj
Ol
CO
"O
cfl
o
S-I
T3
01
[>
cfl
C
tD
^-^
S-i
cd
01
j>^
j_i
o>
ft
CO
^i
ct)
T3
K^
J_|
13
v^^
01
S-i
3
4J
CO
H
0
e
0)
O
ct)
M-l
J^l
3
C/3
j>^
4J
.,-1
^>
H
4J
O
cfl
c:
o
H
4J
O
3
^-i
1 t
CO
d
o
o
M-l
O
d
o
r-l
4-J
Cfl
^
3
Q
O
H
4-1
U
3
4J
CO
d
o
o
cu
^>
H
4-1
O
Cfl
M-l
0
CO
CU
^4
0
d
0
r-l
4->
O
3
4J
CO
d
o
o
r**t
^
CO
01
ffi
y^
4J
d
O)
4J
d
o
0
4-J
r 1
H
CO
\^s
0)
^1
3
4J
X
4-1
01
O
cfl
M-4
J^
3
C/D
T)
01
r-l
r-l
r-l
4-1
CO
ft
O
}_|
O
M-l
O
CO
01
S-i
o
TTJ
d
CO
i t
i i
cd
S-i
3 00
M d
r-l -H
3 r-l
O r-l
l-l T-1
S-I 4J
M
<
1
ct)
S-i
O
ft
cfl
£>
w
1
d
o
H
4J
cfl
4-1
r-l
ft
r-l
CJ
O)
S-i
pu|
V '
Ol
^
3
4-1 ^
co X
r-l 01
O 13
B d
M
O>
o d
cfl O
M-l -H
S-I 4J
3
CO
^*^
4J
d
CD
4-1
d
o
u
4-J
r i
H
CO
CU
3
4-1
X
01
4-1
0)
O
Cfl
M-l
j_l
3
C/l
/J-^S
O
H
V_X
M-4
M-4
O
0)
cfl
4-)
^- CO
c>0 Ol
d T-^
iH 0
13 >^
d o
cfl
Cfl
ft
H
j_j
4-1
CO
^j
r-l
cd
4-1
S-I
i-t
Q
^^
S-i
cfl
01
K^
J_j
0>
ft
CO
^»
cfl
T3
j>.
S-<
13
\^f
Ol
S-i
3
4J
CO
r-t
O
B
0)
a
cfl
M-4
S-i
3
-------
The following sections of this report present, for each source
category, a comprehensive description of the calculation methodology.
The final data tabulations including the coded NEDS area source data
forms, have been submitted as a separately bound attachment to this
report.
-------
UNPAVED ROAD EMISSIONS
SOURCE EXTENT
The basic equation for the calculation of annual traffic flow on
unpaved roads in a specified county is given by:
v = 365 (ADT) m (1)
where v is the vehicle miles traveled per year; ADT is average daily
traffic on unpaved roads in the county; and m is mileage of unpaved
roads in the county. The procedure used to determine ADT and m for
each county is depicted in Figure 1.
Regression analysis of statewide traffic counts for unpaved roads
in Kansas yielded the following equation:
ADT = 15 + 2.8 I Z ) (2)
where p is the county rural population?.' and a is the county
area (sq mile).^-/ Kansas was the only state which was found to have
actual ADT data for unpaved roads. For this reason and because Kansas
is thought to be fairly representative of areas of the country with
substantial mileage of unpaved roads, the above ADT correlation was
applied to all of the other states.
Tabulations of the mileage of unpaved roads (surface types A
through E) in each state are prepared annually by the Department of
o /
Transportation.' However, the county statistics must be obtained
individually from each state.
-------
1
1 1
1 1
1 »
11 .
- -s
2 °
o
«*
1 » "S
1 -C >
I HI D
1 > =
1 rvi ~~>
t>J . /
1 fc c
1 - °
1
(.
-
.
1
<.
_
"D
4) X
0) > "$Z
O c
0 CL 2.
.- C O
-s^u
> c >>
_ o -0
5 « «
2 m -o
Ji §
< -<; oi
*
X
E
^
0
J
X
0
-D
(LI
>
o
o
c
^ ^
c
0
*/»
"D
5 8
< c§
i
<
o!
00
CN
+
IO
1
Q
X
- j
3
X
D
I
|
i
o o 5
Q) " O I
< "1 \
"? R i
o o o
_l - 0.
1/1
1 HI
D Q- 1 "o
| 4-
«
1
-5
&
D
Ji
^5
5
[
X
"c
o
U
X
_Q
O
Cl
CD
D
HI
2
0)
cn
o
01
2
(LI
-»-
P
t/^
X
-*
c
^
o
u
X
^Q
\f>
-r-l
V
o
o
a:
D
(LI
|
D
)
~o
o
o
c£
-o
(LI
^
O
Q.
C
3
,
O
tj
0
LJ_
O)
c
'c
_o
u-
8.
9-
-------
Unpaved road mileage data by county for five states (Kansas,'
Nebraska,' Arkansas,' North Dakota,U and South Dakota^/) were
analyzed by plotting mileage density (m/a) versus rural population
density (p/a) for each state. It was observed that for p/a > 8 persons
per square mile, the mileage density becomes essentially independent of
rural population density. For the counties with p/a < 8, further
analysis leads to the conclusion that the dependence of mileage density
on population density was too small to justify development of complex
correlations to be applied to the relatively few sparsely populated
counties. Therefore, for all states (other than the five mentioned
above) the state unpaved road mileage totals were apportioned among
the respective counties on the basis of county area.
CORRECTION FACTORS
The emission factor for dust emissions from unpaved roads-i' con-
tains three correction factors: (1) average vehicle speed; (2) sur-
face silt content; and (3) rainfree days per year.
Based on previous field testing experience,' the average vehicle
speed on unpaved roads was taken to be 40 mph.
Silt Content
The average surface silt content of unpaved roads for each state
was calculated using the following equation (see Figure 2):
s = 3 (DE)(15) + (NS)s
r 3(DE) + NS
where sr is the weighted surface silt content; DE are miles of roads
with surface types D (soil-surfaced) and E (slag, gravel, or stone-
surfaced); NS are the miles of nonsurfaced* roads in the state; and
s~ is the average soil silt content for the state. The value for s~ ,
which represents the silt content for vehicles traveling on dirt roads,
was determined by averaging the county soil silt contents determined
for calculation of agricultural tilling emissions (see Agricultural
Tilling Correction Factors). The constant 15 represents the percent
silt for D- and E-surfaced roads.' The factor 3 is an estimate of
the ratio of vehicle miles traveled on D- and E-surfaced roads compared
to that on nonsurfaced roads.
* Nonsurfaced roads include primitive (type A), unimproved (type B),
and graded and drained (type C) roads.
-------
o c o
06 S. o5
~O C Q.
O *
co co
I
C .^-^
JU CO
"c Z
o >
U +
CO '
11
^^
co
z
Q
CO
- >-"- (11
D S, 30)
« g3 2 o
1 S i ^
o - U-" ^
co uj
Z Q
1 t
,^
0)
_i_
"c
£ "S
c o
o o
u i
CO uj
^ Q -0
S-, o o
ii: az
V)
> -o
°1I
1 1- D
Q D «
^ irt 1 i/>
X uj 0 0
CO Q Z S
\ '
b
CO
_D
"c
_!»_
c
o
U
_
M
CO
'o
co
-o
D
o
O
o -o »
1 «J 0)
^ o
CN O (0
CV C "
»^ >^^ -^
Y*
-,
D
)
T3 ^
.- O)
U- C
^^ -1
^ 0)
^ 1
I
t- -o
^^
l^J Q
"^^ Q^
O 0
*" D
A
"a c
_E D
jj_
UJ O
C
tfl
a.
a
01
l-l
3
-------
Dry Days per Year
The starting point for the determination of the number of dry days
per year for each county was a national map of average annual frequency
of measurable precipitation.&/ If the number of dry days in a given
state varied by less than 20, an average was taken for the state;
otherwise, the state was divided into regions for which the difference
was less than 20, and an average for each region was estimated. Finally,
the dry days for each county were taken to be the average value for the
region or the state in which that county is located.
EMISSION FACTOR
The emission factor for dust emissions from unpaved roadsi' is
given by (see Figure 3):
EFr = 0.49(sr) I < (4)
The factor gives the pounds of dust particles smaller than 30 urn in
diameter* (based on a particle density of 2 to 2.5 g/cnr) emitted by a
vehicle traveling at a speed of 40 mph over a distance of 1 mile.
* The approximate effective cutoff diameter of a standard high-volume
particulate sampler.!/
-------
X
_Q
1 £
c S
o> £
t £
o **
U «
"D
J± O
CO £^
P §
E o §
LU
Q.
4)
0) II
^"g
JJ O
o E
OO
Mo"
74-
x
|
o
U
X
-Q
CD
_Q
E «
3 ^
Zo
Q
1 ^
§Q
c
a
W
CO
0)
S-i
3
60
10
-------
HEAVY CONSTRUCTION EMISSIONS
SOURCE EXTENT
The most comprehensive available data on the extent of heavy
construction* in the United States are state construction receipts
as derived from the 1972 Census of Construction.2/ Values are broken
down by SIC subcategory. Figure 4 diagrams the procedure which was
used to go from state construction receipts to acres of construction
by county.
After consultation with construction statistics experts, it
was decided that the most reasonable technique for estimation of
county acres would consist of:
(1) conversion of state construction receipts to equivalent
acres of construction,
(2) apportionment of state construction acreage to counties
on the basis of the fraction of the state construction employment
assigned to each county.2J
The conversion factors for step 1, as presented in Table 2, were
developed by MRI for each SIC subcategory. These were applied
separately to the state construction dollars in each SIC code to
determine the acreage of active construction during the base year
(1972). It should be noted that construction dollars for certain
SIC categories were missing for a few states, and therefore were not
included in the computation of total construction acreage.
EMISSION FACTOR
To determine a state-wide emission factor for dust emissions from
construction activities, it was necessary to multiply the previously
Although heavy construction normally is defined as nonbuilding
construction, both building (SIC Code 15) and nonbuilding (SIC
Code 16) construction were included in this section.
11
-------
o .o x
t3 "c
3D
O £ O
3 */> C)
c c
c o -^
< u Jr
X
_Q
0)
O)
£. '" x
P < c
00 D
^" O
u- C ( j
O 0) v-'
- ^^
C >- O
o _o o
O Q.OJ
-------
Table 2. CONSTRUCTION DOLLARS-TO-ACRES CONVERSION FACTORS
Factor
SIC code (acres/$106)
1521 5
1522 5
1531 5
1541 5
1542 5
1611 25
1622 25
1623 5
1629 150
determined emission rate factor (1 ton/acre/month)'by an average dura-
tion of construction projects within each state. As indicated in
Figure 5, the average duration was determined from the relative propor-
tions of acreage differentiated by project type and the average duration
of construction for each project type. MRI economists estimated the
average duration of construction to be:
6 months for residential buildings,
11 months for nonresidential buildings,
18 months for nonbuilding construction.
Therefore, the emission factor for heavy construction can be
written as follows:
EFC = D ton/acre (5)
where D is the weighted average duration of construction within a
given state. Note that this factor describes emissions of particles
13
-------
i_ C
o .2
"o "*""
LL- E
+-
0)
"o
>w
vt v
£- -Q
«
o i
.2L ^
'o '55
>
0)4 >
P
LO
-°
.2oo
$ u <^
*i S x
UJ Lji J3
4
._
4
(J
*o .2
* ^-^
III
2 c o
Q^^
*-l \^ ^^^
1
X
_Q
r- <"
c ~D
.2 o
<*. V, ^
^ o
0 Su
to *- -^
CD « co
i_ C
"Ox
< U _Q
, t ,
0)
"o
LO
_Q"
0)
M
O
-~-,
en
C
o
4-J
N»X
c
o
r-(
4J
f)
3
S-i
4J
CO
C!
o
o
o
i [ |
jj
o
J-J
o
cfl
M-l
d
o
CO
en
i-i
g
W
u~i
(-1
1 1
w>
i*. So
0 So
vi ._ E m -
^ LLJ ts«*» r
14
-------
smaller than 30 urn in diameter* (based on a particle density of 2 to
2.5 g/cm3).
The approximate effective cutoff diameter for a standard high-
volume particulate sampler.'
15
-------
AGRICULTURAL TILLING EMISSIONS
SOURCE EXTENT
The annual acres of land tilled is the designated measure of
source extent for fugitive emissions from agricultural tilling. Data
used for this determination were (see Figure 6):
(1) harvested cropland by county, in acres; and
(2) an estimation of the number of annual tillings, by crop.
The acres of harvested cropland for all farms on a county basis
was obtained from the 1969 Census of Agriculture.A2/ The number of
annual tillings for major crops was estimated by knowledgeable MRI
personnel (see Table 3). An overall value of three tillings per year
was determined to be representative for all cropland. Therefore,
the acres of land tilled in each county was calculated to be three
times the annual harvested cropland.
CORRECTION FACTORS
Two correction factors were calculated for agricultural tilling:
(1) agricultural soil silt content; and (2) Precipitation-Evaporation
Index (a measure of average surface soil moisture content).
Silt Content
Three soils maps were used in determining agricultural silt
content by county. Figure 7 shows the procedure for this calculation.
A map of the soils of the North Central United States-^!/ was
used as the main source of data for the agricultural belt. This map
classifies soils according to their soil series (the most specific
soil classification unit). Map numbers indicate a predominant soil
series and one or two secondary soil series for each respective
geographical area.
17
-------
c
D
C
< X
_Q
C
< P
-D °
.
_£!
TJ
_O
Q.
O
i~
U
~O
0)
ft
(D
^
c
o ^_
:n «
-------
Table 3. ESTIMATED ANNUAL TILLINGS BY CROP
Crop
Barley
Corn
Cotton
Oats
Sorghum
Soybeans
Wheat
Number of tillings
per year
3
3
4, 3 (East, West)
3
2, 3 (East, West)
3
3, 2 (East, West)
19
-------
f ,
I 1
~o *"
O *
LU O
l_ CO
o
<*- (1)
"c **~ *-
c (j
E t
£ o o
00 U Z
t
-5 m
O 4-
LU O
i« CO
^o
*o D
1 s 1
c *£ (j
O (1)
U <" _c
«^ i
CO CO Z
t
I
1/1 r-
V C
£ o
3. '= £
oo 0 .-
_ .- <1>
w. o it ">
£ | 8
Q- « r? oo
0 U U
*> -s~
"^ -^~ n
oil i-
00 Z_ LL. O
<5 .E ^
c
t/1
5 x
i
.. OO
M- (D
C "
J-
^_
^ c
c ^5
Q_
3 0)
o^ U 6
13
-c u
u
o
I
_
/1\
II 1 l_ H*
1_
J «-
i£ ~E ^
*- -2 ««
J-
j_ -
^J
i -h-
-
-------
The following steps were taken to convert soil series information
to silt content values:
i 2 /
(1) The soil series were redefined into soil families ' (the
second most specific classification of soils, indicating the texture
of the soil);
(2) The soil texture triangle' was used to estimate silt
content for each family classification;
(3) A representative value of percent silt content was
determined for each map number by weighting the silt content of each of
the one to three soil series silt content values.
The silt content for each county in the 12-state North Central
region was then determined by the following procedure:
(1) The one or two soil map numbers covering the largest
areas within the county were identified and the fraction of the county
that each covered was estimated.
(2) If the map numbers identified in step 1 covered more
than two-thirds of the area, step 3 was skipped.
(3) The soil map number covering the third largest area
within the county was identified and the fraction of the county that
it covered was estimated.
(4) From the fractions of county area determined above,
and the silt content for each soil map number, a weighted value for
soil silt content was calculated.
The second map encompassed the Great Plains region of the United
States-"' and contained several states not included in the first map.
In this case, map numbers specified the soil family classification.
This information was used along with the soil texture triangle I^/
to assign an estimated silt content to each map number. The silt
content for each county was then determined in the manner described
above for the North Central states.
The third map was a soils map of the entire United States. i'
This map indicated only broad soil classifications: orders, suborders,
and great groups. No general procedure was available to determine
family (or texture) classifications from either suborders or great
groups; for this reason, it was not possible to use the soil texture
21
-------
triangle to estimate silt content. The methods detailed in the para-
graphs below were used to estimate silt content for each great group
(where possible) or suborder indicated on the national soils map.
First, for each soils area outside the states covered by the
first two maps, an attempt was made to locate the designated great
group (or suborder) in one of the previously defined soils areas. If
the desired great group (or suborder) was located, the soils area in
question was assigned the corresponding silt content.
If it was not possible to locate a similar soil for which the
silt content had been previously estimated, a set of scattered soil
texture profiles^-/ was searched. Each of these measured profiles
characterizes a soil suborder, i.e., the dominant soil suborder at the
location where that soil sample was obtained. If a profile for the
corresponding suborder was located, the soil in question was assigned
the silt content from the profile.
Finally, if neither of these procedures was possible, an
estimate of the silt content was made for the suborder based upon the
silt content of the surrounding area and the silt content of other
suborders within the same order.
After each of the suborders was assigned a silt content, the
silt content for each county was determined using the method described
previously for the North Central states.
It should be noted that the method for estimation of silt content
from the national soils map is less accurate than the methods which
utilized the first two maps. Thus, the confidence level of estimates
of silt content for areas not covered by the North Central states
map and Great Plains map is lower. However, agricultural tilling in
areas outside those maps is also less significant.
Precipitation-Evaporation Index
Thornthwaite's Precipitation-Evaporation index' is used to
correct emissions for geographical differences in soil mixture. A
map of PE-index by state climatic division was generated from an
earlier MRI study.!/ A value of the PE-index for each county was
determined by assigning all counties in each state climatic division
the value assigned to that state climatic division. Weighted values
of the PE-index were determined for those counties which were part of
more than one state climatic division. State maps -Is/ indicating both
the state climatic division and the counties were used for assigning
values and weighting functions. This procedure is outlined in Figure 8,
22
-------
c
£
i-
E
8.
o
llj |v
1 D
C O
.2 *kJ
"o x
o o>
0) -D
x
s §5
-------
EMISSION FACTOR
The emission factor for dust emissions from agricultural tilling
is given by the following equation:
EFt = iS (6)
t (PE/50)2
where the symbols are defined in Figure 9 and a value of 5.5 mph has
been substituted for average implement speed.
Equation (6) estimates the total emissions of dust particles
smaller than 30 urn in diameter"" (based on a particle density of 2 to
2.5 g/cm3).
* The approximate effective cutoff value for a standard high-volume
particulate sampler.
24
-------
~O 3
.2
?
*E X
O)
10
IT)
II
u_
X
4
3
O
U
X
_Q
1 ,
CD
"c
o
U
^-x
"c
U
0)
*-*
i!.
3
~3
U
-1- U.
~ 0)
oo <.
U)
-o
0)
(U
Q.
LO
,
C
E
0)
Q.
_E
OO
-------
UNPAVED AIRSTRIP EMISSIONS
SOURCE EXTENT
The measure of source extent for fugitive emissions from un-
paved* airstrips is annual LTO (landing/take-off) cycles. Figure 10
illustrates the procedure used to determine LTO cycles on unpaved
airstrips by county.
Contacts were made with the Federal Aviation Administration to
determine availability and accessibility of data on:
(1) the number of LTO cycles at small airport facilities,
in relation to the number of based aircraft: and
(2) the number of aircraft based at unpaved airstrips
in the United States, by county.
A computer tape with data on each airporti^-' was obtained from
the Washington offices of FAA. Data on this tape included the follow-
ing information for each airport in the United States: site number,
city, state, airport name, county code, latitude, longitude, airport
type, number of total based aircraft, number of multi-engine based
aircraft, runway pavement, runway length, population served, ownership
type and usage type.
The computer tape was processed and punched cards were generated
containing the pertinent information for those airports listed as
pavement type 5. (This pavement type includes dirt, water, sand, and
gravel pavement.) The punched cards were then sorted to eliminate air-
ports with no based aircraft, airports no longer in operation, and
heliport or seaplane bases. Data for over 1,000 airports remained, and
these cards were sorted by county and state codes. It was necessary to
convert the state and county codes from the FAA system to the SAROAD
coding system. This was manually accomplished by code comparison using
an IBM manual2.0-/ and the SAROAD Station Coding Manual./
* Excluding grass (turf) airstrips.
27
-------
O
§ Q - "-
Z < <
8.
<
X
_Q
O
U
L.
<
"D
V
v>
,2
x t- c~5
U " fo
.b CN
O < ii
jq o
g ^ E
QJ LLI
CO
(-1
rH
fi
o
o
>,
o
O
H
3
60
H
28
-------
Local FAA officials estimated the number of operations per based
aircraft at small airport facilities to be in the range of 400 to 800
operations per year with a typical value being 500, i.e., 250 LTO
cycles. Thus, the total number of LTO cycles on unpaved airstrips in
each county was calculated by multiplying the total number of aircraft
based at unpaved airstrips by 250.
CORRECTION FACTORS
The emission factor for dust emissions from unpaved airstrips
contains five correction parameters (see Figure 11):
(1) the surface texture, measured as percent silt content;
(2) the average LTO speed, estimated to be 40 mph;
(3) the surface soil moisture as determined by annual
number of dry days;
(4) the length of runway used for one complete LTO cycle,
estimated to 1 mile; and
(5) a wind-erosion multiplier, estimated to be equal to 2.
The silt content on a state basis and the annual number of dry
days on a county basis were assumed to be the same as those developed
as correction factors for unpaved roads. The estimates for average
LTO speed and length of runway used" in one LTO cycle were derived
from conversations with local FAA officials. The wind-erosion multi-
plier is an estimated value which accounts for the emissions generated
by the propeller wash.
EMISSION FACTOR
The emission factor for unpaved airstrips (Ib/LTO cycle) was
derived by analogy to the equation for unpaved roads. The equation for
unpaved airstrips is given by:
EFa = 2
0.49 si.
(7)
Predominate use of unpaved airstrips is limited to single-engine
aircraft.
29
-------
> 52
o
LL- < X
C *. §
°£ o
s Q u
ko
LL.
I
I1
< v-»
"S
0)
CL
OO
0
J
0)
0)
E
0)
OO
X
f-
c
3
O
u
X
_Q
p
O
sr
n
-o
0)
"o
E
4-
V)
LU
L_
0)
_Q
Z
O
c
c
°
X
0
Q
X
Q
M-
o
-£ < "^
c ~o ~
cu X p
-i U
0 0 II
^ -" "S
°^ .E o
o « LU
CN
.1 »
irt s- ~O
O OJ 0)
^1.1
.E "5 -
^s^
0)
r-(
O
O
O
H
in
cx
CO
S-i
H
CO
H
13
S-i
O
M-l
CO
^
O
4-1
O
d
o
H
CO
CO
H
I
3
60
30
-------
DATA TABULATIONS AND ACCURACIES
Two types of data tabulations were prepared as part of this
study:
(1) coded NEDS area source data forms, listing source
extent for each designated source category, on a county-by-county basis,
(2) tabulations of the correction factors and the corrected
emission factors for each source category, on a county-by-county basis.
Example completed forms (for Alabama) are shown in Tables 4 and 5.
A listing of data specificity of the submitted tabulation is pre-
sented in Table 6. As indicated, for some correction factors, single
values were assigned to multi-county regions or to states as a whole
rather than to individual counties.
The annual number of dry days and silt content for unpaved roads
and airstrips (mostly D- and E-surfaced) do not vary sharply enough
to justify calculation of separate values for each state. In the
case of duration of construction activity, county construction data
were not available.
Table 7 presents estimates for possible error in the calculated
data. These values correspond to a 90% confidence level. They were
determined by a progressive analysis of the errors associated with each
calculation step. Separate values are presented for the source
extent and the corrected emission factor for each source category.
31
-------
^
ll
3
3
K
o
<
3
fit
o
OS
o
- **
.9- $
£ o
± u
^0
i |
5 -d-
^^
CN
=
O
»
ao
rv
0
^
»
r»)
~
-
>
O
*^
ro
s9
Q
0
IV
>
0-
rl
0
0
0
*
ft
o
**,
0
0
r-
>
0
HI
m
0
0
O
«o
ro
ro
I/)
0
0
0
xS
4
o
o
4-
-
o
0
0
W
-X
*-
t-
0
Q
0
J
"1
o
r-
o
-
OQ
CN
O
o
0
te
J
~
ir
0
0
o
0
*
0
Oo
^0
^
vS
o
0
r-
vj
0
0
0
CO
o
^v^
m
0
O
O
Oo
ir\
o
rf
r>
^
sS
Q
^-
0
0
o
*v
'V
-N^
K)
o
0
o
o
-9
0
o-
^>
O
v5
w
0
0
o
o
^
^^
-1-
o
o
0
-S
4
0
rl
^
lo
\s
o
o
0
o
>
-X.
-
0
0
0
t»
sS
0
rl
-
^x
O
O
o
^
Ifl
>s
D
0
0
d
N
o
P~
-
to
5-
O
0
o
l»1
0
V.
IY^
0
0
Q
>
t~-
O
0
On
:S
o
0
s
N
o
-r
r
-
H
r-
^
0
0
0
r<
O
0
o
vB
^
"X.
r
0
o
o
0
OQ
0
^~
en
vS
N
o
0
0
sS
0
»»
0
0
o
H
0°
0
-9
O
10
o
0
0
4
9-
n
O
o
o
*-
oo
0
J
Oo
-
0
o
o|
C
0
0
0
")
**.
J
0
o
o
00
00
O
rS
o
(>-
Ul
O
0
o
N
Q
"S^
<\
O
O
0
o
ff-
0
«0
Vo
V-
r-
m
O
0
o
o
H
4
r^
o
o
0
-t-
v
0
0
o
o
>0
>-
J)
o
o
o
s.
H
^
o
0
o
xfi
-^
sO
V
lo
O
O
Q
00
(O
00
^v
^
>
O
0
Q
^
«h
^>
i\\
O
C
o
o
H
O
^3
«)
0
Q
0
-
i-
10
O
0
o
"0
r-
^~
o
0
o
o
4)
r
r
to
r-
o
0
0
>0
0
^
~r
0
o
o
O
r-
7~
00
-~
f j
o
X.
o
0
vB
OQ
vo
-J>
0
O
o
<»
t~
^
*
-
«0
^.
rl
>
O
0
0
O
U
0
i-H
r-l
3
PQ
rx
1-1
CD
r-i
4-1
3
PQ
CO
d
3
O
J3
i-H
CO
o
0
CO
rl
d
3
0)
i-H
O
0
01
cu
14H
>4H
O
a
K
4J
0)
l-H
O
00
S
U
I
o
o
o
tfl
en
o
0
U
o
fVJ
vN
d
0
ec
d
rl
o
o
Jj
J
U
0)
o
CN
d
to
i i
3
CO
4)
i-H
CO
a
*
a
i-i
O
IT)
i 1
sS
01
Q
"O
a)
|H
l-H
N
FVI
vN
CO
H
.0
8
CO
U
ca
w
CO
CN
43
CO
§
4-1
O
4)
4-1
CU
CO
fa
O
d
rl
t-H
d
CO
J-i
CO
0)
0)
O
CN
01
d
0)
0)
rl
CO
0)
rH
CO
S3
*
CO
rl
d
CU
m
CO
d
o
4-1
en
3
O
33
32
-------
c
^ 0)
i_
o
»- <
1°
1 *^s
.2 "^
S* o
t <
4/> ^-* '
O
11
«/> >***
°- v
1
<^ £
.t u
^0
5 ^
^-*»
0)
,/>
"%*
&
o -5
2>
^o
Z) O
Qi
u
O
<
tE
D
O
u
O)
*-
to
rs
("1
-O
S
s
<">
a
t-i
Si
R
00
IN
s
Z
00
^
*Q
£
*
2
£
-
0
o*
CD
r-.
0
'
o
~
-
*-
^
H
Ml
t-
r(
0
0
o
<0
^
^
h
0
O
O
N
n
ry
In
3-
o
O
0
do
v
7-
O
O
0
rt
O
H
m
-Jk
^
in
(^
+
o
o
0
*
*
d
r-
c
o
0
00
0
1
>s
o(
Of
1/1
r^
r|
0
o
o
X
0
o/
N
o
o
0
0
~
H
9^
U\
t-
lo
rt
c
O
o
M
*-
«J
O
O
0
^
-
H
>
-
"}
^
(T
r|
O
O
o
^
T
r^
0
o
o
^9
-
r(
(N
0-^
vS
«i
>*
o
0
0
ff-
0
j
o
o
0
H
r|
a
*v
X
X
0
0
o
Q
*
r^
0
0
o
H
rfi
ft
f)
1-
-
Oo
IV)
"1
O
Q
o
-
Q
.V
r-
0
o
o
-*
i»>
p<
t-i
in
-
o
*
M
N.
O
o
o
U)
-
>
10
0
o
0
o
*
If
4)
to
-v
r
m
X
O
O
0
o
*.
X.
-
0
o
0
r/
J-
t
0
^
-
0
a
Q
'^
Q
0
O
:S
~
ij
c^
O
O
0
H
00
r^
(0
to
O
0
Oo
r*
o~-
o
0
o
H
^
t
U)
rO
-
_
(T
O
O
0
«0
^
-
r
0
O
0
-3
t-
rt
V)
I"
~-
*^
~-
«
o
o
o
^
Q
^
rxl
O
O
O
0
00
0
v
^
0
o
O
s
«^
^
-
in
i^
0-
t
0
O
o
>
")
i
T
O
O
Q
,9
0
m
CO
O
-v
^
t^
O
0
0
W)
r/
^
T
0
o
o
00
O
«>
-a.
w^
^
H
M
O
O
O
*
^
^
-
0
o
0
,8
~*-
K~>
^
00
^
*0
X
X
.',
O
O
"0
o
X
r~
0
o
O
o
l/>
^
^
0
HI
o>
0)
r-l
U
C!
O
U
g
§
J3
rt
T-l
£
00
X
c
o
(J
Jackson
CM
d
o
Jeffers
CO
I
*
0)
Lawrenc
^o
Monroe
0
^
M
4J
M
AJ
a
s
IX
Morgan
00
M
M
0)
O
Pickens
o
CN
01
-1
PH
CN
r;
Randolo
CM
CN
T-l
i 1
01
w
M
co
CN
M
H
rt
tH
U
4-J
^J
>
a
j:
c/:
LO
tx
)-
a
4.
E
c/:
a
ct
"ej
E-
IX
CN
rt
CO
Q
Tallapoi
00
CN
rt
CO
Q
Tuscalo
Ov
Walker
o
CO
d
o
Washing
Wilcox
CM
Winston
CO
CO
*»
c^
in
CO
33
-------
1
C6
0
t*J
&
O
^J
to
0
H
C/3
c/3
M
1-4
w
s
§>
tJ
1
CO
O
%
M-i
5=5
B
y
a
cu
r-l
^
H
CO
Oi
0
1
u_
o
CO
co
CO
Qf
0
tj
U.
z
U
t
U
LU
CxL
o
U
OJ
a >.
^ - U
- 1
< \
"^
o ^~*
5 a, S
.± c u
3 ^^0
C P .a
05 ZI-
-^
'g
Sl-s
§-<£<
~
O 4>
Q >
>- ^
Q a!
c c
o o
3 0 X
i_ LU
Q- 1
C
-2? ^
O 0 ^cj
*- 1
t/>
U 3
J 3#
= <
en I
00
t
Jn
to
^
^
0
J
ll
r
0
r
0
V3
O
^
>
in
^
o
£;
J
^
t
0-
.0
O
r/
n
rJ
N
r-
-
xS
**
N
^
\n
")
^
r-
^
(V)
rJ
r-
?
^
t
r-
-
r-
1
-
i^
M
r-
(T-
r
^
rl
r-
i/)
^
^
r-
r*
">
^
^
s
(T)
J
r-
ifl
N
*
*>
>
V)
O
r-
IA
(Y)
.1
r,
O
"0
H
r-
ir>
^
.5
r-
Ul
^
r-
o
r-
O
rl
IV
rV
r-
Q-
^
J
r-
in
rr
vO
r
a-
n
rJ
r-
o
r-
^
c-
>
£
ro
r-
r
-
>
ft
*
-
X
>
0
r-
vS
^
IX
1-
\fl
^
*
^
^
In
n
s
v
^
^
^
^
V
O)
o
Q_
e
§
JQ
(0
i-l
CO
"c
o
U
Autauga
CN
Baldwin
CO
3
O
CO
"t
f-l
pfl
CO
Blount
^o
Bullock
IX
Butler
oo
Calhoun
o-
Chambers
o
Cherokee
Chilton
CM
Choctaw
CO
Clarke
<<
i-H
U
IT)
Cleburne
o
Coffee
IX
Colbert
CO
Conecuh
CK
Coosa
0
c
Covingto
CN
Crenshaw
CM
CN
Cullman
CO
c\
o
CV
Elmore
rx
CN
Escambia
00
CN
Etowah
cs
CM
Fayette
o
CO
Franklin
CO
Geneva
CM
CO
Greene
CO
CO
0)
rH
ta
-
-------
13
01
§
r-l
u
o
o
,0
tfl
H
CO
C*.
o
1
u
z
O
CO
CO
LU
CO
O
h-
u
ii
z
0
I
u
u 1
o;
O
u
4>
Q. X
^_ - u
O ' *
5 o> S!
.£ c u
D ^?
u \
'= ^ £
~
u E
S "8 _c'
S-j S
;§ ">
S-I
c c
0 O
^5 O X
"i Hf
u °
ff J5
0- |
c
V ^
C "D - -
tj O o"
^ a:
£ '
e S
U 3
U .- 3s
»/> i
00
^
0
f^
(t%
fO
^
0
si
H
f
r
Q
n
rv-
a
r/
1"
s?
(T
tr
^
^
00
Vo
N
o
IVi
vS
X
^
o
1-
>
0
r
r
^
X
r-
o
r-
IT
^
r-
o
r-
00
*
r-
r-
vS
X
^
o
^
r-
J?
o
*~
r-
J
-S
^
I-
tc
0
rJ
r
J
j
r^
r^
"i
rJ
r>
fi
Ov
,
^
tv>
o
^
rJ
J
xS
^
rJ
rJ
rv
rv
,
,x
^
r
o-
rJ
r
IX
^
ff-
^
T
r
o»
j
rl
4
Cr-
X.
,9
r-
^
Q
X
\
/
\
7
r~.
"~^
^
0
0)
O)
o
O-
^o
OJ
"B
o
o
o
0)
0
s
-O
cd
rH
2
>^
"c
D
O
U
Jackson
CN
c
o
Jeffers
c~
Lamar
TT
o>
r-l
CO
a
01
3
>J
IT)
o
aJ
rx
0)
c
Limes to
00
Lowndes
o*
Macon
o
Madison
Marengo
CM
Marion
CO
1-1
Marshal
<*
Mobile
w>
Monroe
>o
^,
(^
^
T4
cu
f
CN
j±
Randolp!
c>
CN
Russell
CO
CN
(-1
H
eg
CJ
CO
^
CN
X
a
a
"O
CN
Sumter
S3
CN
a
e
a
Cl
IX
CN
CO
CO
n
Tallapoc
00
CN
CO
CO
Q
Tuscaloi
o>
CN
Walker
o
CO
c
o
4J
a1
fl
JS
CO
CO
co
Wilcox
CN)
CO
Winston
CO
CO
-
-------
CO
W
>
w
,-1
o
2
rH
CO
O
w
jx^
o
r-l
CO
CO
~
trt
t-^
4J
t-1
O
M-l
H
O
Q)
O.
CO
cd
4_j
cO
Q
0)
CO
4->
CO
G
CU
4-1
G
O
o
4-1
r-l
H
CO
G
0
00
CU
^
cu
4-)
CO
[i
CO
CO
cO
*0
r^
}_)
T>
j_j
m cti
O 0)
K^
M
0) ^
i g.
3
2
4J
G
0
CQ
cu
i t
H
S
1
0)
t 1
o
1-1
r^
0)
>,
0
60
0)
4J
CO
O
O
J-l
3
O
CO
CO
ro
CO
o
**
*^J
C}
^>
cy
ft
c
CD
4-1
O
tfi
d
o
1-1
4-1 4-1
O O
3
G H
O 4J
H CO
U G
cfl o
t-l O
p
4J
G
01
4-1
c
O
O
CO
TJ
O
!-!
0>
cu
^
H
4-1
O
ctf
M-t
O
CO
CU
J_l
O
<3
G
o
r-l
4-1
o
J3
}_)
4J
CO
C
O
o
tilled
CO
ft
o
M
O
U-l
O
CO
0)
r-l
O
-------
Table 7. ESTIMATED ERROR RANGES FOR TABULATED DATA
Estimated relative error
Source category Extent of source Corrected emission factor
±
Unpaved roads
Heavy construction ± 40% ± 50%
Agricultural tilling ± 15% ± 30%
Dirt airstrips ± 25% ± 30%
37
-------
REFERENCES
1. Cowherd, C., Jr., K. Axetell, Jr., C. M. Guenther, and G. A.
Jutze, Development of Emission Factors for Fugitive Dust
Sources, prepared for the U.S. Environmental Protection
Agency, Office of Air and Waste Management, Office of Air
Quality Planning and Standards, Contract No. 68-02-0619,
Publication No. EPA-450/3-74-037, June 1974.
2. County and City Data Book 1972, a Statistical Abstract Supplement,
U.S. Department of Commerce, Bureau of the Census, U.S. Govern-
ment Printing Office, Washington, D.C. (1973).
3. Highway Statistics 1972, U.S. Department of Transportation, Fed-
eral Highway Administration, Washington, D.C., Table M-3,
October 1973.
4. "Number of Miles and Annual Average Daily Vehicle Miles of Travel
on Local Rural Roads in Kansas by County for the Year Ending
December 31, 1972," State Highway Commission of Kansas Planning
and Development Department, August 10, 1973.
5. Nebraska Highway Statistics State and Local Construction Mileage
for 1972, Nebraska Department of Roads, Office of Engineering
Services, Planning Division, Highway Statistics Unit, September
1973.
6. Arkansas Road and Street Mileages, Arkansas State Highway Depart-
ment, Division of Planning and Research, in cooperation with
the Federal Highway Administration, U.S. Department of Trans-
portation, January 1, 1973.
7. "Miles of Unpaved Roads by County," estimated for North Dakota
and South Dakota, personal communication from Charles Mann,
U.S. Environmental Protection Agency, National Air Data Branch,
July 23, 1974.
39
-------
8. Climatic Atlas of the United States, U.S. Department of Commerce,
Environmental Science Services Administration, Environmental
Data Service, U.S. Government Printing Office, Washington, D.C.,
June 1968.
9. 1972 Census of Construction Industries, Preliminary Report, U.S.
Department of Commerce, Bureau of the Census.
10. 1969 Census of Agriculture, County Summary, Table 2, U.S. Depart-
ment of Commerce, U.S. Government Printing Office, Washington,
B.C.
11. "Major Soils of the North Central Region, U.S.A.," a map from
Soils of the North Central Region of the United States,
North Central Regional Publication No. 76, Bulletin 544,
published by the Agricultural Experimental Station, University
of Wisconsin, in cooperation with the U.S. Department of
Agriculture, June 1960.
12. Soil Series of the United States, Puerto Rico, and Virgin Islands:
Their Taxonomic Classification, Soil Conservation Service,
U.S. Department of Agriculture, pp. 1-1 through 1-228, April
1972.
13. "Guide for Textural Classification in Soil Families," supplement
to Soil Classification: A Comprehensive System, Seventh Approxi-
mation, Soil Survey Staff, Soil Conservation Service, U.S.
Department of Agriculture, p. 40, March 1967.
14. "Soils of the Great Plains," soil map, copyright A.R. Aandahl, 1972,
P. 0. Box 81242, Lincoln, Nebraska.
15. "Distribution of Principal Kinds of Soils: Orders, Suborders, and
Great Groups," a map from the National Atlas, Sheet No. 85/86,
U. S. Geological Survey, Soil Conservation Service, Washington,
D.C. (1967).
16. Soil Classification: A Comprehensive System, Seventh Approximation,
Soil Conservation Service, U. S. Department of Agriculture,
August 1960.
17. Thornthwaite, C. W., "Climates of North America According to a New
Classification," Geograph. Rev., 21:633-655 (1931).
40
-------
18. State maps, National Climatic Center, Asheville, North Carolina.
19. "Airport Services Tape," Federal Aviation Administration, Public
Information Center, AIS 230, Washington, D.C. 20591.
20. Numerical Code for States, Counties, and Cities of the United
States, IBM Manual C20-8073-0.
21. Fair, Don H., SAROAD Station Coding Manual for Aerometric Sampling
Networks, U.S. Environmental Protection Agency, Office of Air
Programs, Research Triangle Park, North Carolina, Publication
No. APTD-0907, February 1972.
41
-------
TECHNICAL REPORT DATA
(Please read Instructions on the reverse before completing)
1 REPORT NO.
EPA-450/3-74-085
3. RECIPIENT'S ACCESSION-NO.
4. TITLE AND SUBTITLE
Emissions Inventory of Agricultural Tilling,
Unpaved Roads and Airstrips, and Construction Sites
5. REPORT DATE
November 1974
6. PERFORMING ORGANIZATION CODE
7 AUTHOR(S)
Chatten C. Cowherd, Jr., Christine M. -luenther, and
Dennis D. Wallace
8. PERFORMING ORGANIZATION REPORT NO.
9. PERFORMING ORGANIZATION NAME AND ADDRESS
Midwest Research Institute
425 Volker Boulevard
Kansas City, Missouri 64110
10. PROGRAM ELEMENT NO.
11. CONTRACT/GRANT NO.
68-02-1437
12. SPONSORING AGENCY NAME AND ADDRESS
U.S. Environmental Protection Agency
Office of Air Quality Planning and Standards
Research Triangle Park, North Carolina 27711
13. TYPE OF REPORT AND PERIOD COVERED
Final Report
14. SPONSORING AGENCY CODE
15. SUPPLEMENTARY NOTES
16. ABSTRACT
To determine the air pollution impact of selected fugitive dust sources, a
national emission inventory on a county-by-county basis was developed. Source
categories included were unpaved roads, unpaved airstrips, heavy construction
sites (road and building construction), and agricultural land tilling. Emission
factors, which include correction factors to adjust for local climatic conditions
and properties of emitting surfaces were calculated for each county based m
the results from EPA Contract No. 68-02-0619. Measures of the extent of activity
for each category were derived from available data. The methodology including
procedures used to estimate missing data, is documented in the report. The surce
data were coded on National Emission Data System (NEDS) area source forms.
17.
KEY WORDS AND DOCUMENT ANALYSIS
DESCRIPTORS
b.lDENTIFIERS/OPEN ENDED TERMS C. COSATI Field/Group
Emission Factors
Correction Factors
NEDS
13. DISTRIBUTION STATEMEN1
Release Unlimited
21. NO. OF PAGES
41
20. SECURITY qLA.SS (Thispage)
Unclassified
22. PRICE
EPA Form 2220-1 (9-73)
42
------- |