WYOMINGUNITED STATES
DEPARTMENT OF THE INTERIOR GEOLOGICAL SURVEY
COLORADO
ANALYSIS OF WASTE-LOAD ASSIMILATIVE CAPACITY OF THEYAMPA RIVER, STEAMBOAT SPRINGS TO HAYDEN,
ROUTT COUNTY, COLORADO
WATER-RESOURCES INVESTIGATIONS 77-119
BIBLIOGRAPHIC DATA SHEET
1. Report No. 3. Recipient's Accession No.
4. Title and Subtitle
Analysis of Waste-Load Assimilative Capacity of the Yampa River, Steamboat Springs to Hayden, Routt County, Colorado
5. Report Date
March 19786.
7. Author(s)
____Daniel P. Bauer, Timothy D. Steele, and Richard D. Andersor8. Performing Organization Rept.
N°' USGS/WRI-77-1199. Performing Organization Name and Address
U.S. Geological Survey, Water Resources DivisionBox 25046, Mail Stop 415Denver Federal CenterLakewood. Colorado 80225 ________ __________
10. Project/Task/Work Unit No.
11. Contract/Grant No.
12. Sponsoring Organization Name and Address
U.S. Geological Survey, Water Resources Division Box 25046, Mail Stop 415 Denver Federal Center
____Lakewood. Colorado 8022S __ __ ______
13. Type of Report & Period Covered
_____Final______14.
15. Supplementary Notes
Prepared in cooperation with the Routt County Department of Environmental Health16. Abstracts
An analysis of the waste-load assimilative capacity of the Yampa River from Steamboat Springs to Hayden, Colo., was made to obtain information on the effects of projected waste loadings on this stream reach. Simulations of effects of waste loadings on streamflow quality were made using a steady-state water-quality model. The simulations were based on 7-day low-flow values with a 10-year recurrence interval and population projections for 2010. Model results for December and September stream- flow conditions indicated that a recommended water-quality standard of 0.02 milligram per liter for nonionized ammonia concentration would be exceeded. Model simulations also included the effect of a flow augmentation of 20 cubic feet per second (0.56 cubic meter per second) from a proposed upstream reservoir. The permissible ammonia loading in the study reach could be increased approximately 25 percent with this amount of flow augmentation. Simulations of concentrations of other water-quality variables indicated that water-quality goals proposed for 1978, 1983, or 1985 would not be exceeded.17. Key Words and Document Analysis. I7a. Descriptors
Colorado, Computer models, Data collection, Invertebrates, Low-flow augmentation, Model studies, Nutrients, Simulation analysis, Standards, Waste assimilative capacity, Water pollution control, Water Quality Act, Water storage, Water treatment
I7b. Identifiers/Open-Ended Terms
Benthic-invertebrate diversity, Cold-water fishery, Computational algorithms, Nonionized ammonia, Population projection, Seasonal conditions, Steady-state water-quality model, Steamboat Springs, Colo., Yampa River basin (Colorado-Wyoming)
I7c. COSATI Field 'Group
18. Availability Statement
No restriction on distribution.
19. Security Class (ThisReport)
__ UNCLASSIFIED20. Security Class (This
PageUNCLASSIFIED
21. No. of Pages
7622. Price
FORM NTis-as iREv. 1O-73) ENDORSED BY ANSI AND UNESCO. THIS FORM MAY BE REPRODUCED USCOMM-DC 8265-P74
ANALYSIS OF WASTE-LOAD ASSIMILATIVE CAPACITY OF THE YAMPA RIVER,
STEAMBOAT SPRINGS TO HAYDEN, ROUTT COUNTY, COLORADO
By Daniel P. Bauer and Timothy Doak Steele, U.S. Geological Survey,
and Richard D. Anderson, Colorado Department of Health
U.S. GEOLOGICAL SURVEY
Water-Resources Investigations 77-119
Prepared in cooperation with the
Routt County Department of Environmental Health
March 1978
UNITED STATES DEPARTMENT OF THE INTERIOR
CECIL D. ANDRUS, Secretary
GEOLOGICAL SURVEY
W. A. Radlinski, Acting Director
For additional information write to:
District Chief U.S. Geological Survey Box 25046, Mail Stop 415 Denver Federal Center Denver, Colorado 80225
CONTENTS
PageMetric conversion factors. ....................... VIAbstract ................................ 1Introduction .............................. 1
Purpose and scope ......................... 2Acknowledgments .......................... 2Previous investigations ...................... 4
Description of study reach ....................... 4Description of water-quality models used for analysis. ......... 11
U.S. Geological Survey model. ................... 11Pioneer-I model .......................... 11
Data used for model calibration. .................... 12Data-collection program ...................... 12Stream-reach and diel variations. ................. 14Data-reduction techniques ..................... 20
Calibration results. .......................... 20Factors incorporated into model simulations. .............. 35
Existing stream-reach classification. ............... 35Effluent standards of proposed regional wastewater-treatment plant. 38Flow augmentation of proposed Yamcolo Reservoir .......... 40Population projections for Steamboat Springs. ........... 40
Model simulations. ........................... 42Effluents from existing wastewater-treatment plants and 7-day,
10-year streamflow. ....................... 42Proposed regional wastewater-treatment plant with 1978 standards
for effluent. .......................... 44Proposed regional wastewater-treatment plant with 1985 standards
for effluent. .......................... 56Ammonia-nitrogen concentrations for varying conditions and 1978
standards for effluent. ..................... 56Stream biological conditions ...................... 61Summary. ................................ 65References ............................... 67
ILLUSTRATIONS
Figure 1. Map showing location and areal extent of the Yampa Riverbasin and location of study reach ............. 3
2-5. Graphs showing:2. Daily-mean discharges, Yampa River at Steamboat
Springs, 1975 water year. .............. 53. Daily-mean discharges, Yampa River downstream from
diversion near Hayden, 1975 water year. ....... 64. Annual-mean and mean-annual discharges, Yampa River at
Steamboat Springs (1905-06, 1910-77). ........ 75. Mean-monthly discharges, Yampa River at Steamboat
Springs (1905-06, 1910-77). ............. 8
III
Page Figure 6. Map showing location of sampling sites, study reach of Yampa
River, September 23-24, 1975 ................ 107-39. Graphs showing:
7. Mean temperature and mean concentrations of selected water-quality constituents in the Yampa River, September 23-24, 1975 ................ 15
8. Concentrations of total-coliform and fecal-coliformbacteria in the Yampa River, September 23-24, 1975. . 17
9. Diel variation in concentrations of total-coliform and fecal-coliform bacteria in the Yampa River downstream from campgrounds near Steamboat Springs, site YM-4, September 23-24, 1975 ................ 18
10. Diel variation in concentrations of dissolved oxygen at selected sites on the Yampa River, September 23-24, 1975 ..................... 19
11. Least-squares method of calculating decay rates and concentrations of biochemical oxygen demand, Yampa River downstream from Steamboat II near Steamboat Springs, site YM-5. ................. 21
12. Decay rate of carbonaceous biochemical oxygen demand for given subreaches downstream from reference location YM-4, Yampa River, Steamboat Springs area, September 23-24, 1975 ................ 22
13. Die-off rates of total-coliform and fecal-coliformbacteria for given subreaches downstream from refer ence location YM-4, Yampa River, Steamboat Springs area, September 23-24, 1975 ............. 23
14. Computed versus observed concentrations of ultimate- carbonaceous biochemical oxygen demand in the Yampa River, September 23-24, 1975. ............ 25
15. Computed versus observed concentrations of dissolvedoxygen in the Yampa River, September 23-24, 1975. . . 27
16. Computed versus observed concentrations of organicnitrogen in the Yampa River, September 23-24, 1975. . 29
17. Computed versus observed concentrations of ammonianitrogen in the Yampa River, September 23-24, 1975. . 30
18. Computed versus observed concentrations of nitrite plus nitrate as nitrogen in the Yampa River, September 23-24, 1975 ................ 31
19. Computed versus observed concentrations of totalnitrogen in the Yampa River, September 23-24, 1975. . 32
20. Computed versus observed concentrations of orthophos-phate in the Yampa River, September 23-24, 1975 ... 34
21. Computed versus observed concentrations of fecal- coliform bacteria in the Yampa River, September 23-24, 1975 ..................... 36
22. Computed versus observed concentrations of total- coliform bacteria in the Yampa River, September 23-24, 1975 ..................... 37
IV
Page Figures 23-39. Graphs Continued
23. Peak-day and permanent population projections for theSteamboat Springs area. ............... 41
24. Computed concentrations of ultimate-carbonaceous bio chemical oxygen demand for existing wastewater- treatment facilities and Q7,10-flow conditions in the Yampa River ................... 43
25. Computed concentrations of dissolved oxygen for existing wastewater-treatment facilities and Q7,10-flow conditions in the Yampa River. ...... 45
26. Computed concentrations of ammonia nitrogen for existing wastewater-treatment facilities and Q7,10-flow conditions in the Yampa River. ...... 46
27. Computed concentrations of fecal-coliform bacteria for existing wastewater-treatment facilities and Q7,10-flow conditions in the Yampa River. ...... 47
28. Computed concentrations of ultimate-carbonaceous biochemical oxygen demand using 1978 standards for effluent, Yampa River .............. 49
29. Computed concentrations of dissolved oxygen using1978 standards for effluent, Yampa River. ...... 50
30. Computed concentrations of ammonia nitrogen using1978 standards for effluent, Yampa River. ...... 51
31. Computed concentrations of total nitrogen using1978 standards for effluent, Yampa River. ...... 52
32. Computed concentrations of orthophosphate using1978 standards for effluent, Yampa River. ...... 53
33. Computed concentrations of nonionized ammonia nitrogen using 1978 standards for effluent, Yampa River, assuming 9 milligrams per liter (December) and 2.8 milligrams per liter (September) of ammonia nitrogen in effluent from wastewater-treatment plant ..... 54
34. Computed concentrations of ammonia nitrogen in the Yampa River, assuming 3.0 and 4.0 milligrams per liter of ammonia in effluent from wastewater- treatment plant ................... 55
35. Computed concentrations of fecal-coliform bacteriausing 1978 standards for effluent, Yampa River. ... 57
36. Computed concentrations of ammonia nitrogen using1985 standards for effluent, Yampa River. ...... 58
37. Computed concentrations of total nitrogen using 1985standards for effluent, Yampa River ......... 59
38. Mean diversity of benthic invertebrates and nutrientconcentrations, Yampa River ............. 63
39. Mean total number of organisms per square meter andnutrient concentrations, Yampa River. ........ 64
V
TABLES
Table 1. Existing point-source discharges from wastewater-treatmentplants ...........................
2. Sampling-site locations, September 23-24, 1975 ........3. Nonionized ammonia (NH3) for observed pH and temperature
values, Yampa River, September 23-24, 1975 ........4. Percentage of nonionized ammonia nitrogen in ammonia-water
solutions at various pH and water-temperature values . . . ,5. Suggested effluent standards for proposed regional wastewater-
treatment plant. .....................6. Maximum computed concentrations of ammonia nitrogen and non-
ionized ammonia nitrogen for permanent populations, and using 1978 standards for effluent, Yampa River, Steamboat Springs to Hayden, Colo. ..................
7. Maximum computed concentrations of ammonia nitrogen and non- ionized ammonia nitrogen for peak-day populations, and using 1978 standards for effluent, Yampa River, Steamboat Springs to Hayden, Colo. ..................
8. Periphyton data, Yampa River, Steamboat Springs to Hayden,Colo., August-September 1975 ................
Page
913
33
38
39
60
61
62
METRIC CONVERSION FACTORS
The English units used in the text of this report are followed by metric units; others may be converted to metric units by using the following conver sion factors:
To convert English units Multiply by
foot (ft) 0.3048
foot per mile (ft/mi) .1894
cubic foot per second (ft 3 /s) .02832
mile (mi) 1.609
square mile (mi2 ) 2.590
gallon per minute (gal/min) 6.309xlO~5
To obtain metric units
meter (m)
meter per kilometer (m/km)
cubic meter per second (m3 /s)
kilometers (km)
square kilometers (km2 )
cubic meter per second (m3/s)
VI
ANALYSIS OF WASTE-LOAD ASSIMILATIVE CAPACITYOF THE YAMPA RIVER, STEAMBOAT SPRINGS TO HAYDEN,
ROUTT COUNTY, COLORADO
By Daniel P. Bauer and Timothy Doak Steele, U.S. Geological Survey, and Richard D. Anderson, Colorado Department of Health
ABSTRACT
An analysis of the waste-load assimilative capacity of the Yampa River from Steamboat Springs to Hayden, Colo., a distance of 38 miles (61 kilo meters), was made during September 1975 to obtain information on the effects of projected waste loadings on this stream reach.
Simulations of effects of waste loadings on streamflow quality were made using a steady-state water-quality model. The simulations were based on 7-day low-flow values with a 10-year recurrence interval and population projections for 2010. Model results for December and September streamflow conditions indicated that the recommended 1978 Colorado and 1976 U.S. Environmental Protection Agency water-quality standard of 0.02 milligram per liter for nonionized ammonia concentration would be exceeded. Model simulations also included the effect of a flow augmentation of 20 cubic feet per second (0.56 cubic meter per second) from a proposed upstream reservoir. The permissible ammonia loading in the study reach could be increased approximately 25 percent with this amount of flow augmentation. Simulations of concentrations of dissolved oxygen, fecal-coliform bacteria, and nitrate nitrogen indicated that the State's water-quality goals proposed for 1978, 1983, or 1985 would not be exceeded.
A comparative analysis using the U.S. Geological Survey model and a modified Pioneer-I model, developed by Battelle Pacific Northwest Laboratories, was made with the observed data collected during September 1975. Some variation was shown between the two model results for the constituents considered.
INTRODUCTION
As a result of Public Law 92-500 (Federal Water Pollution Control Act Amendments of 1972), there are many studies now in progress for implementing areawide wastewater-management plans throughout the United States. Within the
State of Colorado, 12 regional councils of government have been established, and among their responsibilities are the coordination of studies dictated by section 208 of Public Law 92-500.
This report was prepared by the U.S. Geological Survey in cooperation with the Routt County Department of Environmental Health. The study reported herein was coordinated with the staff members of the North West Colorado Council of Governments. Both the Yampa River basin assessment by the U.S. Geological Survey (Steele and others, 1976a, 19762?) and the areawide wastewater-management study by the North West Colorado Council of Governments involve evaluating primary and secondary impacts of coal-resource development in their respective study areas (U.S. Department of the Interior, 1976). The study described herein was initiated because of impending population growth in the Steamboat Springs area and the concern for maintaining the regional envi ronmental esthetic quality in order to meet demands for various forms of recreation (including water-based uses) of the area. The specific area of interest for this study is a reach of the Yampa River, approximately from Steamboat Springs to Hayden, located in Routt County, Colo. (fig. 1).
Purpose and Scope
The study was conducted, in part, in fulfillment of guidelines for implementing section 208 of Public Law 92-500. The primary purpose of the study was to aid in evaluating the waste-load assimilative capacity of the study reach for a minimum mean 7-day low flow and 10-year recurrence interval (Q7,10). The results of this analysis will provide planners and managers in Routt County with information for determining the waste-load capacities of the study reach for possible design and operational alternatives of future wastewater-treatment plants.
The report presents an analysis of one set of data collected during a 24- hour period. These data were used to calibrate two water-quality models: A U.S. Geological Survey model and a modified Pioneer-I model developed by Battelle Pacific Northwest Laboratories. Calibration of the two models was done to evaluate the adequacy of the various model computational algorithms. The report also discusses simulation results using the U.S. Geological Survey model for alternative levels of waste loadings, Q7,10-flow conditions, and seasonal conditions during September and December. Discussion of possible relationships between measured concentrations of nutrients and observed populations of macroinvertebrates and periphyton also is included.
Acknowledgments
The authors are grateful to Phillip E. Stark and his staff, Routt County Department of Environmental Health, for assistance in the planning and execu tion of the field-data collection and in providing supplemental information; and to Ronald M. Eddy, Aquatic Biologist, Technical Investigations Branch, Surveillance and Analysis Division, U.S. Environmental Protection Agency,
u>
XLI
TTLE
ISN
AK
Ea.
/1
\
H{
r
Din
osa
ur
Natio
nal
Mon
umen
t
Bas
e fr
om U
.S.
Geo
logi
cal
Surv
ey
1:50
0 00
0 S
tate
bas
e m
ap,
1969
1 '
0
10
2
0
30 K
ILO
ME
TE
RS
40°<
Pro
pose
d Y
amco
lo
rese
rvo
ir s
ite
Fig
ure
1. L
ocat
ion
and
area
l ex
tent
of
the
Yam
pa R
iver
bas
in a
nd l
ocat
ion
of s
tudy
rea
ch.
Region VIII, for coordinating an extensive evaluation of macroinvertebrates in the study reach of the Yampa River. The authors also thank the many residents of Routt County for access to the Yampa River and tributary sampling sites during the study.
Previous Investigations
A study by McCall-Ellingson and Morrill, Inc. (1974), which was com pleted for the Colorado Department of Health, developed a water-quality management plan for that part of the Green River basin in Colorado. Their study was designed to fulfill the requirements of section 303(c) of Public Law 92-500. In their study, the reach of the Yampa River from Steamboat Springs to the confluence with the Elk River (fig. 1) was classified as an area with significant pollution problems, and an extensive computer-model analysis was made. Conclusions of the study by McCall-Ellingson and Morrill, Inc., were: Concentrations of ammonia nitrogen exceeded the recommended standards of the U.S. Environmental Protection Agency (National Academy of Sciences-National Academy of Engineering, 1973) for almost the entire study reach, concentrations of dissolved oxygen exceeded the minimum recommended standards of the State for the entire reach, and concentrations of fecal-coliform bacteria did not exceed recommended standards of the State. Their analysis was based upon data collected at relatively high flows during July 1972, and did not include data for low-flow conditions.
In 1974, Wright-McLaughlin Engineers conducted a study in the Mount Werner ski area, located northeast of Steamboat Springs. The purpose of the study was to evaluate the use of treated sewage effluent for snow-making on the ski slopes. Results of this study (Wright, 1976) were favorable toward the use of treated sewage in this manner. A separate summer land-treatment study was completed by Wright-McLaughlin Engineers (1974) for the Timbers Water and Sanitation District, also located near Steamboat Springs.
DESCRIPTION OF STUDY REACH
The 38-mi (61-km) study reach of the Yampa River begins approximately 3 mi (4.8 km) east of Steamboat Springs and ends 2 mi (3.2 km) west of Hayden (fig. 1). Daily-mean discharge of the Yampa River for the 1975 water year for the gaging station at Steamboat Springs is shown on figure 2 and for the gaging station near Hayden on figure 3. Mean-annual discharge for the Yampa River at Steamboat Springs for the period of record (1905-06 and 1910-77 water years) is 464 ft 3 /s (13.1 m3 /s) (fig. 4). Mean-annual discharge for the 1975 water year was 9 percent higher than the long-term average (fig. 4). The variability of mean-monthly discharges for the Yampa River at Steamboat Springs for the period of record is shown on figure 5. Low flows generally occur from September to February (fig. 5). The streamflow regime for the study reach of the Yampa River during low-flow periods is pool and riffle. The slope of the channel (11 to 26.4 ft/mi or 2 to 5 m/km) results in long
OcrP 1-1
era(T)C/3
2CD1-1
p
p 3 o* o p
OQ
DISCHARGE. IN CUBIC FEET PER SECOND - k» woo S
S So g o o
o o o
Q in -
8-
o
w
> ^
c _Z en m
C _» r~ en
> _,i"
I II I o o
I I I1 I
o o oO O
DISCHARGE, IN CUBIC METERS PER SECOND
9DISCHARGE. IN CUBIC FEET PER SECOND -j
- NJ 01 P-» NJ tn o o o o
g 888 888 8 8~i i n i i i i r
oc
3 <D P3
Q-
WOcr & -!
OQ
3
&
3>-*> -I
O3d-
O 33
>-i
sc
3
.. . . iiiiiiiiio O-*NJW o> rtfcJW 2 rtL, &, ?o b b b g P P P 8 0 o o o b
DISCHARGE. IN CUBIC METERS PER SECOND
Mea
n =4
64 c
ub
ic f
eet
per
seco
nd
(13
cub
ic m
eter
s pe
r se
cond
)
100 19
0519
1519
2519
3519
4519
5519
6519
75
Figu
re 4
. A
nnua
l-m
ean
and
mea
n an
nual
dis
char
ges,
Yam
pa R
iver
at
Stea
mbo
at S
prin
gs (
1905
06,
191
0 77
).
ME
AN
-MO
NT
HL
Y D
ISC
HA
RG
E.
IN C
UB
IC F
EET
PER
SE
CO
ND
oo
OQ c S n>
so o 3 Er £T CL
tf3 O
O
(W
ON
ft> Co r* C/3 3 or o
a r* 00
00
ME
AN
-MO
NT
HL
Y D
ISC
HA
RG
E.
IN C
UB
IC M
ETER
S PE
R S
EC
ON
D
riffles. Annual variability in stream temperatures in this reach of the Yampa River ranges from 0° to 20°C (Celsius) (T. D. Steele, D. P. Bauer, D. A. Wentz, and J. W. Warner, written commun., 1977). Minimum water temperatures generally occur in February and maximum temperatures in August. During the winter months, parts of the study reach may be completely covered with ice.
The locations of the sampling sites used for the modeling analysis of the study reach are shown on figure 6. Seventeen tributaries discharge into the Yampa River in this reach (fig. 6). Major tributaries include the Elk River, Trout Creek, and Walton Creek. In 1975, it was estimated that 8,000 people lived along the study reach, principally in the Steamboat Springs and Hayden areas (Gathers and Associates, Inc., 1976). Effluent discharges from municipal wastewater-treatment plants comprise the major source of organic pollution to the Yampa River in the study reach. A summary of the existing point-source discharges from wastewater-treatment plants in the study reach is presented in table 1. Until late 1975, a 180-MW (megawatt) powerplant upstream from the town of Hayden (fig. 6) discharged heated water into Sage Creek, a tributary to the Yampa River. During 1976, a second unit generating 250 MW was added to the powerplant. At about that time, a cooling-tower and evaporation-pond type of system was installed to dispose of the heated water resulting from the operation of the entire powerplant. No thermally heated waters are now being discharged into Sage Creek from the powerplant.
Table 1 .--Existing point-source discharges from wastewater-treatment plants
Reference code 1 Source System type Entities served
YE-1
YE-3-
YE-4-
YE-5-
YE-6-
YE-7-
Mount Werner Aerated lagoon
Town of Steamboat Springs
Aerated lagoon
KOA Package plant
Sleepy Bear Package plant
Steamboat I I Package plant
Mount Werner Water and Sanitation District, Tree Haus district.
Town of Steamboat Springs, Fish Creek Water and Sanitation Dis trict, West Steamboat Water and Sanitation District, Riverside Water and Sanitation District, and Mount Werner bypass.
KOA campground. Private system.
Sleepy Bear mobile-home park. Private system.
Steamboat II Water and Sanitation District.
Town of Hayden
Aerated lagoon Town of Hayden.
figure 6.
R.8
8W.
R.8
7W
.
107°
00'
R.8
6W.
R.8
5W.
R.8
4W.
T.7
N
40°3
0
T.6
N
YM
-6
Ste
ambo
at S
pri
ngs^
Bas
e fr
om U
.S.
Geo
logic
al S
urve
y C
raig
, 19
74
EX
PL
AN
AT
ION
V
SA
MP
LIN
G S
ITE
ON
YA
MPA
RIV
ER
V
SAM
PLIN
G S
ITE
ON
T
RIB
UT
AR
Y
WA
ST
EW
AT
ER
-TR
EA
TM
EN
T P
LA
NT
10 M
ILE
S
10
KIL
OM
ET
ER
S
Fig
ure
6.
Lo
cati
on
of
sam
plin
g si
tes,
stu
dy r
each
of
Yam
pa R
iver
, S
epte
mbe
r 23 2
4,
1975
.
DESCRIPTION OF WATER-QUALITY MODELS USED FOR ANALYSIS
Two models were used in the calibration phase of the analysis: A U.S. Geological Survey model and a modified Pioneer-I model developed by Battelle Pacific Northwest Laboratories. The purpose of this comparison was to show any differences in the computational algorithms used in the models for the same model-parameter values.
U.S. Geological Survey Model
The U.S. Geological Survey model is documented by Bauer and Jennings (1975). The model is based upon the oxygen-sag equation of Streeter and Phelps (1925). The model is intended for use under steady-state (constant flow rate of waste and stream discharges) conditions. When applying the model to a given reach of a river, the reach is divided into a given number of subreaches, generally defined by locations where a waste or tributary inflow mixes with the river water. The formulation assumes complete mixing in a stream cross section.
Some of the original programing concepts made for this form of model have been described by Shindala (1972). The first large-scale application of this model was in Arkansas (Jennings and Bryant, 1973).
The original model (Bauer and Jennings, 1975) was expanded to include processes for determining concentrations of fecal-coliform and total-colifora bacteria, the various nitrogen components, and orthophosphate. The equations for determining concentrations of coliform bacteria that were incorporated are those described by Mahloch (1973), equations for nitrification reactions are described by Thomann (1972), and the equation for determining concentrations of orthophosphate is described by Willis, Anderson, and Dracup (1975).
Pioneer-I Model
The Pioneer-I model (Waddel and others, 1973) is based upon the DOSAG-I model, originally developed by the Federal Water Pollution Control Administration and modified by the Texas Water Development Board (1970). The Pioneer-I model, which also utilizes the oxygen-sag equation of Streeter and Phelps (1925), was developed for stream studies relating to dissolved-oxygen and biochemical oxygen demand dynamics.
The mathematical basis of the Pioneer-I model is formulated around concepts of steady-state flow and uses Lagrangian coordinates to determine concentrations of various constituents for a given flow condition, waste- source concentration, water temperature, reaction coefficients, and water- diversion data. The water-quality dynamics usually are characterized by algorithms of first-order differential equations, although some algorithms of this model have the option for using second-order functions.
11
The Pioneer-I model was used in this study by the Colorado Department of Health to compare model results with those obtained with the U.S. Geological Survey model. The Pioneer-I model currently used by the Colorado Department of Health has been further modified with respect to the original Pioneer-I model developed by Battelle Pacific Northwest Laboratories (Waddel and others, 1973) for the U.S. Environmental Protection Agency. The additional modifications to the Pioneer-I model were made by the URS/Ken R. White Co. (1975) to facilitate the use of the model in situations where data are scarce. These modifications make it possible to study discrete segments of a river rather than an entire river basin, to add or remove discharges without redefining the entire model, and to allow streamflow in reaches or subreaches of the river to decrease to zero without physically removing these reaches or subreaches from the model input.
The Pioneer-I model recently has been modified again by the Colorado Department of Health to estimate the concentrations of the following constituents: Total nitrogen, dissolved solids, metal ions, orthophosphate, fecal-coliform bacteria, ammonia nitrogen, nitrite nitrogen, nitrate nitrogen, carbonaceous biochemical oxygen demand (CBOD), chlorophyll-a, and dissolved oxygen (DO).
DATA USED FOR MODEL CALIBRATION
Data-Collection Program
Field data used for calibrating the models were collected at 33 sites along the study reach during 24 hours on September 23-24, 1975. Samples were taken at approximately 3-hour intervals. The basic data collected for this study are tabulated in a report by Giles and Brogden (1977). A total of 16 main-stem, 11 tributary, and 6 wastewater-effluent sites were sampled (table 2 and fig. 6). Only 10 of the 17 tributaries in the study reach were sampled because the remainder had flows of less than 0.1 ft 3 /s (0.003 m3 /s). At all sites, field determinations were made for DO, water temperature, specific con ductance, fecal-coliform and total-coliform bacteria, and pH. In addition, water samples were collected for laboratory analysis of biochemical oxygen demand (BOD), organic nitrogen (nonfiltered form), ammonia nitrogen, nitrite- plus-nitrate nitrogen, Kjeldahl nitrogen, orthophosphate, total phosphorus, and total nitrogen. At least one discharge measurement was made at each sam pling site during the 24-hour sampling period.
As stated previously, data for this investigation were collected during a single 24-hour period. As noted by Rickert and Hines (1975), there is an inverse relationship between the complexity of a water-quality process and the ability to produce useful modeling assessments within a given time span. The Yampa River from Steamboat Springs to Hayden is not extremely complicated with respect to the flow regime or water-quality characteristics. Ideally, data for modeling would be collected during different seasons of a year; however, data collected during critical-flow conditions are sufficient to adequately calibrate a model (Hines and others, 1975). Both flow and water-quality char acteristics in the study reach normally vary considerably over a yearly cycle.
12
Table
2.--
Samp
ling
-sit
e lo
cati
ons,
September
23-24, 19
75
[YM=Yampa
Rive
r main-stem
sampling lo
cation;
YT=tributary sa
mpli
ng location;
YE=wastewater-effluent sampling lo
cati
on.
See
figu
re 6
for
rela
tive
location]
Reference
code
YM-0
Yj-2
YM-0
1
YE-1
YT-3
Yj-4
YM-1
YT-5
Yj-6
YM-2
YE-3
YM-3
YE-4
YE-5
YM-4
YE-6
YM-5
YM-6
YM-7
YT-8
YJ-9
YM-8
YT-1
0
YM-9
YM-1
0
YM-11
YT-15
YM-1
2
YM-1 3
YT-1
?a
YE-7
YT-1
?b
YM-1
4
Stat
ion
identifier
4025
4410
6493
600
4027
0010
6485
400
4027
3710
6493
700
4027
5310
6494
500
4027
5910
6493
100
4028
5710
6494
000
09239500
4029
4410
6495
900
4029
2010
6501
900
4029
3410
6505
400
403000106515700
402958106515200
403014106523300
4030
1710
6524
900
4030
1710
6525
800
4030
0410
6544
700
4030
0210
6545
500
402932106564900
402902106580000
4029
1310
6580
400
402903106584100
402840107004200
402816107003800
4028
5410
7020
500
402902107043600
09244410
402918107094400
4030
5110
7124
500
4030
0610
7154
800
402939107160100
402952107161500
4029
5210
7161
600
4029
3101
7174
200
_,.
...
Dist
ance
fr
omSt
at ion
na
me
. .
/ . ,
\
i mouth
(mi le
s) L
Yamp
a River
belo
w Oa
k Creek, near St
eamb
oat
Spri
ngs - - -
Walt
on Creek
near mouth, at U.S. Highway
40 br
idge
-- - -
Yamp
a River
below
Walt
on Creek -- - - --- -
Mt.
Wern
er se
wage
disposal po
nd--
-- ---
--- ________________ ___
____
_
Fish Creek
near
mouth,. at
St
eamb
oat
Springs-- -- -
Spring Cr
eek
near mouth, at St
eamb
oat
Springs-------------- ---
- ____
Yamp
a Ri
ver
at St
eamb
oat
Spri
ngs-
- ---
----
- -
Butcherknife Cr
eek
near
mouth, at St
eamb
oat
Spri
ngs-
- --------------
Soda
Creek
near
mouth, at St
eamb
oat
Springs----- - - -
Yampa
Rive
r above
wastewater-treatment plant, be
low
Stea
mboat
Spri
ngs-
Stea
mboa
t Springs
wastewater-treatment plant -- - - -
Yamp
a River
belo
w wastewater-treatment plant, below
Steamboat
Springs--
KOA
Campground pa
ckag
e plant - - -- -- -- -
Sleepy Bear Park pa
ckag
e pl
ant-
--- ---
-- - - -
Yamp
a Ri
ver
below
camp
grou
nds,
near St
eamb
oat
Spri
ngs -
Stea
mboa
t II
wastewater-treatment plant- ---
-- --- ---
----
----
--- -
Yamp
a Ri
ver
below
Steamboat
II,
near St
eamb
oat
Spri
ngs-
----
----
----
--
Yamp
a River
abov
e Elk
Rive
r, near Milner-- - -
Yamp
a Ri
ver
above
Elk
Rive
r co
nflu
ence
, near Milner-- - - --
Elk
River
near mouth, at U.S. Highway
40 br
idge
-- -
West
Fork Elk
Rive
r ne
ar mouth, at U.S. Hi
ghwa
y 40
bridge- - ---
Yamp
a Ri
ver
at Milner- - - - - - - - ---
Trout
Cree
k near Mi
lner
- - - - -
Yampa
Rive
r below
Trout
Creek, at
Mi
lner
- - - - - - -
Yamp
a Ri
ver
abov
e Tow Creek
oil
fiel
d---
- ---
----
----
- - ___
____
__
Yamp
a Ri
ver
belo
w diversion, near Hayden-------------------------------
Sage Creek
near Hayden - - ___
____
__
Yamp
a Ri
ver
below
Morg
an Creek, ne
ar Hayden - - - - -
Yamp
a Ri
ver
at Ha
yden
----
----
- ___
_ ____ ________ ___ ______
Dry
Cree
k above
wastewater-treatment plant, at Hayden- -
Hayd
en wastewater-treatment plant- - ---
-_ ___
_
Dry
Creek
below wastewater-treatment plant, at
Ha
yden
- - ---
Yamp
a Ri
ver
below
Hayd
en - - -- - - --
195.5
193.
1192.4
191.8
191-7
190.
619
0.5
190.3
189.
8189.3
188.3
188.3
187.
5187.4
186.7
185.5
185.3
183.4
181.
818
1.8
181.
1179.0
178.4
177-3
174.
6
169.2
169.
0163.8
159.
615
8.7
158.7
158.
7157.5
^or tribut
ary
and
wast
ewat
er ef
flue
nts,
"distance
from m
outh
" indicates
Yamp
a Ri
ver
distance for
the
trib
utar
y or
wastewater-effluent
junc
tion
.
The sampling was done during a 24-hour period to obtain the approximate diel variation for most of the constituents. Only two BOD and nutrient samples were collected at the 11 tributary locations because few diel varia tions were expected. The main-stem sites YM-9, YM-10, YM-11, YM-12, and YM-13 (fig. 6) are located downstream from the Elk River which contributes approx imately 40 percent of the total flow of the Tampa River at the confluence. These sites were sampled only for BOD and nutrients because water from the Elk River was expected to further dilute the wastewater effluent in the Yampa River below the confluence.
A second data-collection effort for a 24-hour period was made during De cember 1976. Because of the time constraints of this investigation, analysis of these data is not included in this report. The second effort was made to analyze the yearly variability of water-quality and flow characteristics and also to serve as a check of model results during winter (December) conditions.
Steady-state conditions can be assumed for most streams during low flows by dividing a reach into subreaches. Such was the case during the calibration sampling period for the study reach of the Yampa River. To help define depth, velocity, and traveltime data for each subreach of the Yampa River, 15 discharge measurements were made. These measurements were made so that average depth-velocity conditions of the stream in the general vicinity of each measurement could be determined. The discharge of each tributary and wastewater-treatment plant also was measured.
Stream-Reach and Diel Variations
Profiles of the mean values of DO concentrations, stream temperature, and 5-day CBOD (CBOD 5 ) are shown on figure 7. This summary is based on data collected during the 24-hour sampling period. The mean DO concentrations in the study reach were within ±5 percent of saturation, with the mean DO concen tration exceeding 8 mg/L (milligrams per liter) at each location. Mean DO concentrations were larger just upstream and downstream from the main wastewater-treatment plant at Steamboat Springs (sites YM-2 to YM-5, fig. 7 and table 2). These large concentrations may be due, in part, to the large number of algae and submerged vascular plants at the sites. Stream temperatures during the 24-hour sampling period averaged about 10° to 11°C from Steamboat Springs to Sage Creek. Downstream from Sage Creek (downstream from site YM-11), the stream temperature in the Yampa River increased to approximately 13°C. The discharge from Sage Creek (site YT-15), which was receiving thermal-heated effluent from the Hayden Power Plant, had an average temperature of 16°C and contributed to the increased temperature in the Yampa River downstream from Sage Creek (sites YM-12, YM-13, and YM-14, fig. 6). The slope of the Yampa River decreases between Sage Creek and Hayden. As a result, the mean river velocity decreases, creating larger pools in this part of the study reach. This results in a longer residence time per unit length of stream and allows the water in the stream to approach its equilibrium tem perature more quickly for a given length of stream reach.
14
14
FT
12 10
SE
LEC
TED
SIT
E N
UM
BER
S
I I
I^
CM
CO
*
T T7
T
S
2|U
2
>-
> >
- >
UJ
UJ
o. o
V)
MS
uj
<
UJ
oc o
e
Sat
ura
ted
jdjs
solv
ed o
xyge
n
Mea
n d
isso
lvo
fi o
xyge
n
UJ O O
O
4\ \
\
J\
195
190
185
180
175
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
MIL
ES
I__
__
__
__
__
__
__
I
170
165
160
310
300
290
280
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETER
S27
026
0
Fig
ure
1'. M
ean t
empe
ratu
re a
nd m
ean
conc
entr
atio
ns o
f se
lect
ed w
ater
qual
ity c
onst
itue
nts
in t
he Y
ampa
Riv
er,
Sep
tem
ber
23-2
4, 1
975.
The profile for average concentrations of CBODs also is shown on figure 7. Except for a subreach downstream from the main wastewater-treatment plant at Steamboat Springs (site YE-3) and another subreach downstream from the Elk River (site YT-9) near Milner, Colo. (table 2), the concentrations of CBOD5 are less than 3.0 mg/L. The larger values downstream from Steamboat Springs were expected because of the effects of waste discharges, but the relatively large values downstream from the Elk River are not readily explained. One possible cause of the larger concentrations of CBOD5 downstream from the Elk River is irrigation return flow in this area of the study reach. While the 24-hour survey was being completed, small volumes of return flow in this area were noted, but no samples for analysis of BOD were collected.
The concentrations of total-coliform and fecal-coliform bacteria along the study reach are shown on figure 8. As shown, the higher coliform-bacteria concentrations occur downstream from the Steamboat Springs area where the larger proportion of treated-wastewater effluents discharge into the Yampa River. The two profiles, in general, indicate a gradual die-off of coliform bacteria with increasing distance downstream from Steamboat Springs. The reason for the larger concentration of fecal-coliform bacteria in the vicinity of site YM-12 is not known.
Variations in concentrations of total-coliform and fecal-coliform bacte ria for the 24-hour sampling period at site YM-4 are shown on figure 9. Site YM-4 is located downstream from the effluent discharge of the main wastewater- treatment plant at Steamboat Springs (fig. 6). The large coliform-bacteria concentration at the site is the result of the upstream discharge of waste- water effluents. Concentrations of total-coliform bacteria varied consider ably during the 24-hour sampling period; for example, coliform bacteria at site YM-4 (fig. 9) ranged from 380 to 3,500 colonies/100 mL (milliliters) of water.
The diel variations of DO concentrations are shown on figure 10 for sites YM-2, YM-6, and YM-11 (fig. 6 and table 2). A general cyclical pattern is observed at all three sites, with highs during daylight hours typifying photosynthesis in the stream, and with lows at night indicating relative effects of respiration. The largest diel variation (7 to 10.5 mg/L) occurred at site YM-6, which is located approximately 5 mi (8 km) downstream from the discharge outlet of the main wastewater-treatment plant at Steamboat Springs. Site YM-2 is located upstream from this outlet. A diel variation of 8.2 to 10.2 mg/L in DO concentration occurred at site YM-2 during the 24-hour cycle. Data at this site include the effects of the discharge from the wastewater- plant at Mount Werner and any nonpoint sources from Steamboat Springs. Site YM-11 is located upstream from Sage Creek on the Yampa River (fig. 6). A diel variation of 7.3 to 9.7 mg/L in DO concentration occurred at this site during the 24-hour cycle. Along the entire length of the study reach, and especially in sections containing pools, varying concentrations of floating and attached plants were observed. In some subreaches, these plant concentrations were larger than in others, with the largest concentrations in the vicinity of Steamboat Springs.
16
T UJ uiS
iu
5
SE
LEC
TED
SIT
E N
UM
BE
RS
CM ^ S>
->-
>-
103
S IU O z
o
o
102
101
195
190
185
180
175
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
MIL
ES
I__
__
__
__
__
__
__
I
170
165
160
310
300
290
280
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETE
RS
270
260
Fig
ure
8.
Co
nce
ntr
atio
ns
of t
ota
l co
lifo
rm a
nd f
ecal
coli
form
bac
teri
a in
the
Yam
pa R
iver
, S
epte
mbe
r 23 2
4,
1975
.
8T
CONCENTRATION, IN COLONIES PER 100 MILLILITERS
_» o w
I I I I I I I I I I I I I
TO C 1-1 COVO
<ft pps 1-1 3 p"
3 | 3 °o »'
P o3 o >0 3
TO O1-1 CDO p
a O 3 Hft> ^pj W
" OC/3 ^
P r
3 51 a- V
Is £ I 2 1-13 3"B
^T ft) « OK-S1
^ 1
13
o 32 ^ 3 P a- oft> J*
12
IS
10
2
8
O i
UJ
4O
4
x
0 1200
1800
2400
TIM
E,
IN H
OU
RS
0600
1200
Fig
ure
10.-
-Die
l va
riat
ion
in c
on
cen
trat
ion
s of
dis
solv
ed o
xyge
n at
sel
ecte
d s
ites
on
the
Yam
pa R
iver
,S
epte
mbe
r 2
3-2
4,
1975
.
Data-Reduction Techniques
The CBOD in water samples collected during the 24-hour study on September 23-24, 1975, was analyzed using a procedure described by Hines, McKenzie, Rickert, and Renella (1977). Using this technique, DO concentrations in the sample bottles were measured initially and then remeasured 1, 2, 3 r 4, 5, 6, 8, 10, 12, 15, and 20 days after the initial measurement. These data were analyzed using a computer program developed by Jennings and Bauer (1976). Sample results of this form of analysis are shown on figure 11. By this technique, the ultimate CBOD (CBODU), CBOD 5 , and decay- rate coefficient (K^) are computed. In the example cited, CBODU=8.0 mg/L, CBOD5=2.9 mg/L, and #i=0.09 per day. An average-bottle decay-rate coefficient of about 0.1 per day for the entire study reach was computed using the average of the decay rates determined for the samples.
The decay-rate coefficients for CBOD computed by the decrease in DO in the sample bottles with time are somewhat different from those occurring in a stream environment. A technique described by Velz (1970) was applied to the CBOD data determined in the laboratory for the stream reach downstream from site YM-4, which is 1.6 mi (2.6 km) downstream from the main wastewater- treatment plant at Steamboat Springs. The stream reach included sites between YM-4 and YM-8. The log of the percentage of CBOD remaining relative to some initial concentration was plotted versus the estimated stream traveltime down stream from the starting point (fig. 12). By manually fitting a straight line through the data points, the arithmetic slope of the line will indicate the approximate decay rate, KI , of CBOD in the stream reach. A value of 0.4 per day was determined for KI by visual fit of the data points.
A similar analysis also was made for the coliform-bacteria data. Site YM-4 again was used as the reference site with sites between YM-4 and YM-8 considered. The results of this analysis are shown on figure 13. The die-off rates of coliform bacteria estimated by this visual fitting technique were 0.65 per day for fecal-coliform bacteria and 1.03 per day for total-coliform bacteria.
CALIBRATION RESULTS
The two water-quality models, U.S. Geological Survey and Pioneer-I, were calibrated using results of the data collected during September 23-24, 1975. In this section the calibration results are discussed and some variations in results of the two models are presented. To provide direct comparison, the same rate coefficients were used in both models.
Water-quality data used in the model calibration included the following variables: CBODU; DO; organic nitrogen, as nitrogen; ammonia nitrogen, as nitrogen; nitrite-plus-nitrate nitrogen, as nitrogen; total nitrogen, as ni trogen; orthophosphate; fecal-coliform bacteria; and total-coliform bacteria. For modeling analysis, the study reach was divided into a number of subreaches defined by the 10 tributaries, diversion sites, discharge points of wastewater-treatment plants, and geometry of the study reach (fig. 6).
20
I I I I I
Observed biochemical oxygen demand (BOD)
Decay-rate constant (Base E) = 0.091
Decay-rate constant (Base 10) =0.039
Ultimate BOD (mg/L)=S.O
5-day BOD (mg/L) = 2.9
I I1.9 3.8 5.7 7.6 9.5 11.4
TIME, IN DAYS
13.3 15.2 17.1 19.0
Figure 11. Least-squares method of calculating decay rates and concentrations of biochemical oxygen demand, Yampa River downstream from Steamboat II near Steamboat Springs, site YM 5.
21
100
90 80
e>70
z |
60
§
oc a
50 40
<
30
O
hO
(J O CO
20 10
YM
-4
YM
-5Y
M-7
YM
-6
YM
-8
0.0
0.1
0.2
0.3
0.4
TRA
VE
LTIM
E.
IN D
AY
S
0.5
0.6
0.7
0.8
Fig
ure
12
. D
ecay
rat
e of
car
bona
ceou
s bi
oche
mic
al o
xyge
n de
man
d fo
r gi
ven
sub
reac
hes
dow
nstr
eam
fro
m
refe
renc
e lo
cati
on Y
M 4
, Y
ampa
Riv
er,
Ste
ambo
at S
prin
gs a
rea,
Sep
tem
ber
23
24
, 19
75.
PERCENTAGE OF COLIFORM BACTERIA REMAINING
n ln -
o o3 ^
< CU
rff
-2- o -
c/a o
oOQ
o X3
f-> CLso O
Calibration results indicated nonmixed conditions in the stream cross sections at site YM-3. Visual inspection of this location in October 1975 also indicated nonmixed conditions. Site YM-3 is located approximately 500 ft (150 m) downstream from the discharge outlet of the main wastewater-treatment plant at Steamboat Springs. The stream at site YM-3 divides into two channels with approximately 75 percent of the flow in the left (southern) channel and 25 percent of the flow in the right (northern) channel. The right channel, where the samples were collected, receives the largest proportion of the treated-wastewater effluent. Because of this nonmixed condition, all observed data at site YM-3 were adjusted by the approximate flow ratio so that the data could be used in the model-calibration analysis.
Computed versus observed concentrations of CBODU for the study reach are shown on figure 14. The decay-rate coefficient for CBODU used was 0.1 per day for the entire study reach and corresponds to the average-bottle decay-rate coefficient previously computed. Even with using the 0.1-per-day decay-rate coefficient, the computed profile was low compared with the means of the observed data. One explanation could be that some additional BOD loading, such as irrigation return flows, could have affected concentrations of CBODU in the upstream subreaches. The values of CBODU computed by the two models were almost identical for the entire study reach.
The reaeration formula of Bennett and Rathbun (1972) was used in the DO computations. The reaeration formula has the following form:
K2 =8.76 F0 - 607 /!? 1 - 689 , (1) 20
whereK? =reaeration-rate coefficient at 20°C, in 1 per day;
20
F=mean velocity in the river, in feet per second; and
D-mean depth in the reach, in feet.
This reaeration formula expressed the reaeration-rate coefficient of the stream as a function of mean depth and velocity. The temperature-corrected reaeration-rate coefficients ranged from 3.8 to 13.1 per day with the larger coefficients occurring in the upstream subreaches between sites YM-01 and YM- 4. The temperature correction of the reaeration-rate coefficient is of the following form:
where
K2 =K2 (a) T-20, (2)1 20
#2 ,rrreaerat ion-rate coefficient at temperature, T, in 1 per day;
T=temperature , in degrees Celsius; and
(a)=empirical constant, set equal to 1.024 for the U.S. Geological Survey and Pioneer-I models.
24
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
N3
3
20
ffi o. OT I O
-J1
6 12
01 O
Z
4
O
310 I
- _ -
1
YM
-1
<
YM- (
YM
-2
»
<
-
YM
-01
T
-J 1
>
300 1
1
YM
-5
YM
-4
-3
(
>
)
<
YM
YM
-6
' 1
I
(
>
29
0 L 1
-7
YM
-8
l
<>
YM
<
1
280
270
260
1 1
11
YM
- <
-9 »
1 1
1-1
0
YM
-11
S
ITE
NU
MB
ER
-i-
MA
XIM
UM
OB
SE
RV
ED
C
ON
CE
NT
RA
TIO
N
< >
AV
ER
AG
E O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
MIN
IMU
M O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
+
ON
E S
AM
PLE
OR
NO
V
AR
IAT
ION
>
CO
MP
UT
ED
BY
PIO
NE
ER
H
MO
DE
L
CO
MP
UT
ED
BY
U.S
. G
EO
LOG
ICA
L S
UR
VE
Y M
OD
EL
YM
-11
YM~
13
YM
-12
^
YM.
J*~
'1 _
n- <
till
26
0 1 - - _ -
-14
195
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
14. C
om
pute
d v
ersu
s ob
serv
ed c
on
cen
trat
ion
s of
ult
imat
e ca
rbonac
eous
bioc
hem
ical
oxy
gen
dem
and
in t
he Y
ampa
Riv
er,
Sep
tem
ber
23 2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r I
mod
els)
.
Profiles of computed versus observed DO concentrations are shown on figure 15. The DO profiles computed by the two models are almost identical (fig. 15). To compare the DO concentrations computed by the two models required some consideration to the methods by which the data for mean depth and velocity of the stream are processed by each model. The U.S. Geological Survey model uses mean-depth and velocity data entered into the computer directly in the reaeration-rate computations; whereas, the Pioneer-I model computes the mean depth and velocity using the following equations:
(3)and ,
, (4)
where£?=mean flow in reach, in cubic feet per second; and
a, b, ct-i, Z?i=regression coefficients.
In the water-quality comparisons, the hydraulic coefficients a, b, a^, and b\ in equations 3 and 4 for the Pioneer-I model were adjusted until the computed values of mean depth and velocity approximately equaled those used by the U.S. Geological Survey model. In earlier computational comparisons of the two models, somewhat larger differences in the computed DO profiles were due to discrepancies in the traveltime computations in the Pioneer-I model. As indicated by the model computations (fig. 15), the DO sag in the study reach is barely discernible. The computed DO values remained within 95 percent of saturated DO concentrations. Relatively larger diel variations in DO concen trations were observed at several sites, especially in that part of the study reach downstream from the main wastewater-treatment plant at Steamboat Springs (sites YM-3 to YM-7, fig. 15) and downstream from Hayden (site YM-14, fig. 15).
In addition to total nitrogen, the following nitrogen components were modeled: Organic nitrogen (nonfiltered form), ammonia nitrogen, and nitrite- plus-nitrate nitrogen. The temperature-corrected rate coefficients using first-order exponential decay algorithms from the model calibrations were organic nitrogen, 0.2 per day; ammonia nitrogen, 0.3 per day; nitrite nitro gen, 1.0 per day; and nitrate nitrogen, 1.7 per day. Each nitrogen-rate coef ficient was corrected for water-temperature variation by an empirical coeffi cient of 1.047 using an expression similar to equation 2. The rate coefficient of 1.7 per day for nitrate nitrogen applies only to the U.S. Geological Survey model. The Pioneer-I model accumulates the nitrate nitrogen as the final step in its nitrogen-cycle computations.
The model-calibration results give estimates for nitrite and nitrate concentrations separately. These values then were added together and compared to the observed data. Total nitrogen, which was treated in both models as a conservative constituent, was modeled to provide some indication of where the nitrogen was being used or supplied in the study reach other than from point sources.
26
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
ro
310 I
12
oc
t oc Ul s'° oc
O 5 Z
8Z* o i o z 0 o z U
l > 4
o Ul J §
25
0
- -
YM
-
^£*
-
<
1
300
29
0
280
270
260
I 1
I I
I1
t 1
1 1
1
YM-
YM
-2
-01
-
YM
-1
(
>
^
>
\ ' {
26
0 1
_
YM
4
YM
14
~~
--Y
M-5
±
!4
YM
-6 Y
JS r
7
YM
-10
YM
-9
~ "
YM
-13
YIV
-8 --
YM
-11
Y
M-1
2
O
- -
-p-
) °
~~ i r
~~T
pi ___
-t ^
» ' * "
»
»
!! j.
,^^
j |
""""
<>
<>
<>
<>
_ _
-_- -
- -
YM
-11
S
ITE
NU
MB
ER
-
() ^
-r-
MA
XIM
UM
OB
SE
RV
ED
CO
NC
EN
TR
AT
ION
,
( >
AV
ER
AG
E O
BS
ER
VE
DC
ON
CE
NT
RA
TIO
N
_
MIN
IMU
M O
BS
ER
VE
DC
ON
CE
NT
RA
TIO
N
CO
MP
UT
ED
BY
PIO
NE
ER
-I_
MO
DE
L
CO
MP
UTE
D B
Y
-
U.S
. G
EO
LOG
ICA
L
195
| 190
SU
RV
EY
M
OD
EL
I I
I I
I I
185
180
175
170
165
160
15
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
Fig
ure
15
. C
ompu
ted
vers
us o
bser
ved
conc
entr
atio
ns o
f di
ssol
ved
oxyg
en i
n th
e Y
ampa
Riv
er,
Sept
embe
r 23
-24,
197
5 (U
.S.
Geo
logi
cal
Surv
ey a
nd P
ion
eer
I m
odel
s).
The profiles of computed organic-nitrogen concentrations (fig. 16) fit the average observed organic nitrogen within data and model accuracies. The computed ammonia-nitrogen concentrations (fig. 17) are larger than the observed concentrations except at site YM-3. The rapid disappearance of ammonia downstream from site YM-3 is a phenomenon observed on all Colorado streams surveyed by the Colorado Department of Health (R. D. Anderson, written commun., 1976). The computed profiles of concentrations of nitrite-plus- nitrate nitrogen (fig. 18) from sites YM-0 to YM-5 agree within 20 percent of the observed values, but downstream from site YM-5 the computed values were consistently larger than observed values. Computed versus observed total- nitrogen concentrations are shown on figure 19. Observed concentrations were 20 to 40 percent smaller than concentrations computed by both models downstream from site YM-8.
The comparison of the results of computed concentrations of organic nitrogen, ammonia nitrogen, nitrite nitrogen, and nitrate nitrogen agree closely for the U.S. Geological Survey and Pioneer-I models from approximately site YM-01 to site YM-8. For the reach downstream from site YM-8, some variations in the computed values are noted. The larger computed concentra tions of nitrite-plus-nitrate nitrogen for the Pioneer-I model can be explain ed by its accumulation of nitrogen in the nitrate form. The reason for the differences in computed organic-nitrogen and ammonia-nitrogen concentrations between the two models is not known. Some additional tests need to be made with each model to determine reasons for these differences. Computed total nitrogen, which was modeled as a conservative constituent, agreed closely be tween the two models along the entire study reach (fig. 19).
There are several possible explanations of the poor fit of computed versus observed ammonia nitrogen (fig. 17). Willingham (1976) reported that ammonia nitrogen in aqueous solutions exists in two states, ionized NH^ and nonionized NH3. In its nonionized (NH3) state, ammonia can escape as a gas from water. According to Willingham (1976), the partitioning of total ammonia between nonionized NH3 or NH^ forms is dependent primarily on pH and temperature conditions. The estimated ranges of ammonia as NH3 during the 24- hour sampling period are shown in table 3. The percentage of nonionized NH3 in most instances is less than 15 percent, which suggests that the loss of nitrogen as gaseous ammonia to the atmosphere is small. A second explanation is the use of the ammonia nitrogen by plants in the stream environment, for example, by the different algal forms. Kittrell (1969) reported that ammonia nitrogen can be assimilated by algae and changed to organic nitrogen by algae. It was noted that DO is not utilized in this process. Kittrell also pointed out that organic nitrogen, changed to ammonia nitrogen and oxidized to nitrate nitrogen, can be assimilated quickly by algae, which reconvert the nitrogen to protein as an organic-nitrogen form. Dissolved oxygen is used in this latter process. Because of the small oxygen sag in any part of the study reach (fig. 15), it is assumed for the Yampa River that nitrogen most probably is being lost as ammonia. The amount of DO used by the oxidation of organic and ammonia nitrogen to nitrate nitrogen is accounted for by both models, with a rate of 4.57 units of DO for each unit of nitrogen oxidized.
28
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
ro
1.5
BC
UJ £ _l ec
1.2
V) i 3 5 ?
0.9
» z o H t-
5 0 0
0.6
O Z UJ a o f£ z f
O
7
o ce O
0.0
310 I
300
I
2»0
2tO
27
0 260
2SO
T "i H
i ' i i h
YM
-6
YM
^-3
YM
-0 +
195
YM
-2
YM
-1
<
1 -I
4*
)
«
,
I19
0
) ~ Y
M-4
<
YW
__
,
> (
<
>
~ Y
M-1
1
SIT
E N
UM
BE
R
-r-
MA
XIM
UM
OB
SE
RV
ED
CO
NC
EN
TR
AT
ION
< >
AV
ER
AG
E O
BS
ER
VE
D
~C
ON
CE
NT
RA
TIO
N
MIN
IMU
M O
BS
ER
VE
DC
ON
CE
NT
RA
TIO
N
YM
-7
+
ON
E S
AM
PLE
OR
NO
» * ~
VA
RIA
TIO
N
-
CO
MP
UT
ED
BY
PIO
NE
ER
-M
OD
EL
CO
MP
UT
ED
BY
U.S
. G
EO
LO
GIC
AL
YM
-8
SU
RV
EY
MO
DE
L
( i ,
YM
-1 4
|_
u
-- ,Y
i.12
-T
°
YMT
9 Y
M-1
0
YM
-11
' - J
+
^.
--
I I II
I I
185
180
175
170
165
160
15
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
Fig
ure
16.-
-Com
pute
d ver
sus
obse
rved
co
nce
ntr
atio
ns
of o
rgan
ic n
itro
gen
in t
he Y
ampa
Riv
er,
Sep
tem
ber
23 2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r I
mod
els)
.
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETE
RS
1
A
TION.
IN
MILLIGRAMS
PER
LITER
P
~* r*oo b K> !
cc
CO
LU
0
^
0.6
0
0 LU
O fE
0.4
Z z 0
I
0,2
<
0.0 u
310
300
290
280
270
260
260
I I
I I
I I
11
1 1
1 1
1 1
YM
-11
S
ITE
NU
MB
ER
-r-
MA
XIM
UM
OB
SE
RV
ED
C
ON
CE
NT
RA
TIO
N
( »
AV
ER
AG
E O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
MIN
IMU
M O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
+
ON
E S
AM
PLE
OR
NO
V
AR
IAT
ION
-
YM-
3
CO
MP
UT
ED
BY
PIO
NE
ER
-I
MO
DE
L
_
<)
CO
MP
UT
ED
BY
_
U.S
. G
EO
LOG
ICA
L S
UR
VE
Y M
OD
EL
r ~~r
_
YM
-1
YM
-4
YM-
6 " ~~
"
TY
M-2
^"Y
M-5
V
M_,
YM
-9
k I
VJW
-14
YM
-01
JL
+
X
+
+
Y
M-7
Y
M-8
Y
M-1
0
YM
-11
Y
M-1
2
YM
-1 3
'
I "
I __
__
__
__
__
I +
I
+
+
I+ _
_______ L^f
c __
__
__
__
I +
U
+
15
180
185
180
175
170
165
160
1!
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
Fig
ure
17. C
om
pute
d v
ersu
s ob
serv
ed c
once
ntr
atio
ns
of a
mm
onia
nit
roge
n in
the
Yam
pa R
iver
, S
epte
mbe
r 2
3-2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r-I
mod
els)
.
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETER
S
0.10
310 I
30
029
0 I28
027
0 I26
0 I250
_I
YM
-7
n ii n
YM
-11
SIT
E N
UM
BE
R
-1-
MA
XIM
UM
OB
SE
RV
ED
C
ON
CE
NT
RA
TIO
N
AV
ER
AG
E O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
MIN
IMU
M O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
ON
E S
AM
PLE
OR
NO
V
AR
IAT
ION
CO
MP
UTE
D B
Y P
ION
EE
R-1
M
OD
EL
CO
MP
UTE
D B
Y U
.S.
GE
OL
OG
ICA
L
SU
RV
EY
MO
DE
L
YM
-10
+
0.00
YM
-11
+Y
M-1
2+
YM
-13
+
YM
-14
196
190
186
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
MIL
ES
165
160
155
Fig
ure
18.-
-Com
pute
d v
ersu
s ob
serv
ed c
on
cen
trat
ion
s of
nit
rite
-plu
s-nit
rate
as
nitr
ogen
in
the
Yam
pa R
iver
, S
epte
mbe
r 23-2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r-I
mod
els)
.
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
OJ
N3
310 I
2.0
DC
HI
H Hi
DC
HI
Q.
V) 5
1.6
DC o i Z Z
1.2
O H DC H o 8
0.8
Z HI
O
O Z <
0.4
O
0.0
30
0 I1
YM
-3
_ - _
+
(
YM
-2,
YM
-1
i.
(
(
)
1 1 1 j
)
)
LJK
£:
290
29
0
270
260
1 1
1 1
1 III
YM
-11
S
ITE
NU
MB
ER
-P
MA
XIM
UM
OB
SE
RV
ED
CO
NC
EN
TR
AT
ION
( »
AV
ER
AG
E O
BS
ER
VE
DC
ON
CE
NT
RA
TIO
N
MIN
IMU
M O
BS
ER
VE
DC
ON
CE
NT
RA
TIO
N
YM
YM
-4
<~Y
M
)
-5 )
-6
+
ON
E S
AM
PLE
OR
NO
VA
RIA
TIO
N
26
0 1
CO
MP
UT
ED
BY
PIO
NE
ER
-'M
OD
EL
CO
MP
UT
ED
BY
YM_
7 U
.S.
GE
OL
OG
ICA
L
-
SU
RV
EY
MO
DE
L
. -i_
Y
M-8
0
YM
-9
YM
-10
* '
*" *"
T r*
0
YM
-11
YM
~12
Y
M-1
3
J
T
"$"
(*
-
YM
-01
~*
~
1I
I I
I I
I
-
-14
>
195
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
19
. C
ompu
ted
vers
us o
bser
ved
con
cen
trat
ion
s of
tota
l ni
trog
en i
n th
e Y
ampa
Riv
er,
Sep
tem
ber
23 2
4,
1975
(U
.S,
Geo
logi
cal
Sur
vey
and
Pio
nee
r I
mod
els)
.
Table 3« Nonionized arrmoni-a (Nh^) for observed pH and temperature values* Yampa River, September 23-24, 1975
Location 1
YM-O YM-01 VM 1YM-l ---------- YM-2 YM-3
YM-lf VM C
VM £
YM-7 VM Q
YM-9 YM-10 YM-II VM-1 ) ---------
VM 1 '2 I PI I _j---------
VM 1 li
Average temperature, in degrees Cel sius
10 10 10 10 10
10 10 11 11 11
11 11 11 12 12 12
pH range
8 /1 Q -7. H-O. / 8.3-8.8 8.3-9.0 7.1-8.1 7.2-8.7
7.2-8.9 7.5-8.9 8.2-9.1 7.8-9.0 8.0-9.0
7.9-8.7 7.9-8.5 7.9-8.6 7.3-8.4 7.4-8.4 7.0-8.8
PH med Jan
8.6 8.5 8.6 7.8 7.8
8.1 8.28.7 8.2 8.2
8.18.3 8.17.5 7-7 8.2
MH3 range (percentage of total ammonia)
4.5- 8.5 3.6-10.5 3.6-15.7 .2- 2.3.3- 8.5
.3-12.9
.6-12.9 3.1-20.2 1.3-16.7 2.0-16.7
1.6- 9.1 1.6- 6.0 1.6- 7.4 .4- 5.2 .5- 5.2 .2-12.0
See figure 6
In the nitrogen-cycle fitting process by both models, it was the intent to fit the nitrite-plus-nitrate nitrogen and organic nitrogen more closely than the ammonia nitrogen. This was decided upon because there was no possible means to fit the observed ammonia-nitrogen profile with a first-order decay function with the models being used. As was noted earlier, the comparison between computed versus observed nitrite-plus-nitrate nitrogen and organic-nitrogen concentrations in the upstream subreaches of the Yampa River from sites YM-01 to YM-5 was thought to be adequate. Because the ammonia nitrogen is part of the nitrogen cycle, this procedure results in conservative estimates of ammonia-nitrogen concentrations.
A comparison of observed and computed orthophosphate concentrations is shown on figure 20. This constituent represents another critical nutrient for the plant community in the stream system. The orthophosphate-concentration profiles computed by the U.S. Geological Survey and Pioneer-I models differ more in the downstream part of the study reach. Both models appear to compute higher values than the observed average or single values. The orthophosphate- concentration profiles were computed assuming phosphorus to be a non- conservative substance and neglected the effects of nonconservative organic inputs of phosphorus. The orthophosphate computations of both models use two
33
0.!
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S31
0 I30
029
0 I28
0 I27
0 I26
0 I25
0_I
Z UJ Z
O
O
0 08
0.20
OJO
YM
-11
SIT
E N
UM
BER
-T-
MA
XIM
UM
OB
SE
RV
ED
C
ON
CE
NT
RA
TIO
N
C )
YM
-3
<)
YM
-01 +
\_
YM
-1
,-
YM
-6
YM
-7
~~\
0.0
0
...
YM
-S
<> L
___T
C.
I
O
YM
-9
YM
-10
I I
/Y
M-1
1
I
AV
ER
AG
E O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
MIN
IMU
M O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
ON
E S
AM
PLE
OR
NO
V
AR
IAT
ION
CO
MP
UTE
D B
Y P
ION
EE
R-I
M
OD
EL
CO
MP
UTE
D B
Y U
.S.
GE
OL
OG
ICA
L
SU
RV
EY
MO
DE
L
YM
-14
YM
-12
+
195
190
185
180
176
170
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
MIL
ES
166
180
156
Fig
ure
20. C
om
pute
d v
ersu
s ob
serv
ed c
once
ntr
atio
ns
of o
rtho
phos
phat
e in
the
Yam
pa R
iver
, S
epte
mbe
r 23-2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r-I
mod
els)
.
first-order decay coefficients, as described by Willis, Anderson, and Dracup (1975). The first coefficient relates phosphorus losses to benthos, and the second coefficient relates phosphorus use to chlorophyll-a. The temperature- corrected rate coefficients used in the models were 0.45 per day and 0.0 per day for the first and second orthophosphate-decay rates, respectively. The second coefficient relating to chloryphyll-a use was found to have little effect for the model calibrations and therefore was neglected.
The calibration results of the fecal-coliform and total-coliform bacteria are shown on figures 21 and 22, respectively. Considering the accuracy of the coliform-bacteria determinations and the range of observed diel variations (fig. 9), the modeling results give a general idea of profile variations. The results from the two models for the fecal-coliform bacteria were nearly identical (fig. 21). No comparison was made for the total-coliform bacteria because the Pioneer-I model did not include an algorithm for its computation.
One purpose of the model-calibration analysis was to evaluate the various computational algorithms in the models. A second purpose was to derive a unique set of model-parameter values. As shown in the preceding discussion of calibration results, some differences in the computed profiles were noted. The large scatter of the observed field data makes it difficult to derive the unique set of model parameters in the calibration phase for either model. As a result, calibration results are considered only fair. For the simulation phase of the analysis, described later, only the U.S. Geological Survey model was used. This choice was made because of previous experience with the model, the ease with which it can be used, and the form of required input data.
FACTORS INCORPORATED INTO MODEL SIMULATIONS
Existing Stream-Reach Classification
The study reach of the Yampa River-,- from Steamboat Springs to Hayden, has been classified as a cold-water-fishery secondary-body-contact type Bl (Colo rado Department of Health, 1974). Stream-quality requirements for this classification include the following: DO concentrations to equal or exceed 6 mg/L, concentrations of fecal-coliform bacteria to be less than 1,000 colonies/I00 mL, a pH range from 6.0 to 9.0, water temperature not to exceed 20°C, the maximum temperature increase from ambient water temperature not to exceed 1.1°C, and a recommended requirement that nonionized ammonia concentra tions be less than or equal to 0.02 mg/L (Colorado Department of Health, written commun., 1976; U.S. Environmental Protection Agency, 1976). Other variables considered in the standards of the State of Colorado and U.S. Environmental Protection Agency, such as taste, odor, color, and solids, are not included in this analysis.
The nonionized ammonia concentrations are not computed directly by the U.S. Geological Survey or Pioneer-1 models; instead, they are computed manually using reference tables. For the simulation phase of this study, tables in a report by Willingham (1976) were used. These tables express the
35
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETER
S
OJ
400
I
350
§
300
250
200
Z UJ O O O
150
100 50
310 I
300 I
29
0 128
0 I27
0 I26
0 I25
0_I
YM
~4
YM
-3
<>
YM
-S+
YM
-11
SIT
E N
UM
BE
R
MA
XIM
UM
OB
SE
RV
ED
C
ON
CE
NT
RA
TIO
N
<t>
AV
ER
AG
E O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
MIN
IMU
M O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
ON
E S
AM
PLE
OR
NO
V
AR
IAT
ION
-
CO
MP
UTE
D B
Y P
ION
EE
R-I
M
OD
EL
CO
MP
UTE
D B
Y
U.S
. G
EO
LO
GIC
AL
S
UR
VE
Y M
OD
EL
YM
-9
YM
-12
YM
-14
II
I19
519
018
518
0 17
5 17
0
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
21
, C
ompu
ted
ver
sus
obse
rved
co
nce
ntr
atio
ns
of f
ecal
co
lifo
rm b
acte
ria
in t
he Y
ampa
Riv
er,
'Sep
tem
ber
23 2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r I
mod
els)
.
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
U)
3500
£ LLJ U Hi i
3000
o
o cc LLJ
Q.
CO w
2500
z o _l o U z - 20
00o < £ II
I 0
15
00
z 0 o <
E LLJ o
1000
2 5 o i 50°
V 2 o H o
310
300
29
0
280
270
26
0I
______ J__ _________ L, I , I '
I Y
M-4
<
YM
-3 + IN>K
l\
I x
^J YM-2
^
YM
-1
++
YM
-01
-.._
T,
| I
I I
I
YM
-11
S
ITE
NU
MB
ER
MA
XIM
UM
OB
SE
RV
ED
CO
NC
EN
TR
AT
ION
( i
AV
ER
AG
E O
BS
ER
VE
DC
ON
CE
NT
RA
TIO
N
MIN
IMU
M O
BS
ER
VE
D
CO
NC
EN
TR
AT
ION
+
ON
E S
AM
PLE
OR
NO
VA
RIA
TIO
N
CO
MP
UT
ED
BY
U.S
. G
EO
LOG
ICA
L S
UR
VE
YY
M-5
S
UR
VE
Y M
OD
EL
+
)
[\
YM
-9!
x\ \s \
!Y
M-4
[
+
">..
YM
-8+
^^,^
-Y
M-7
^ ^
YM
-10
"~
"~
1^ _
_ _
V
M_
<<
.»
«
-T
'^' __
YM
-12
Y
M-1
3
»"|
TT
i.
i i
i i
i +
- i i u
.---
^
25
0 1 - - - -
*
195
190
1B5
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
22
. C
om
pu
ted
ver
sus
obse
rved
conce
ntr
atio
ns
of t
ota
l co
lifo
rm b
acte
ria
in t
he Y
ampa
Riv
er,
Sep
tem
ber
23-2
4,
1975
(U
.S.
Geo
logi
cal
Sur
vey
and
Pio
nee
r-I
mod
els)
.
nonionized ammonia nitrogen as a percentage of the ammonia-nitrogen concentration based on the water temperature and pH. A short summary of some of the nonionized ammonia-nitrogen concentrations from Willingham (1976) for a range of pH from 6 to 9, and a range of water temperature from 0° to 15°C is presented in table 4. It is noted that the percentage of nonionized ammonia nitrogen increases significantly with increases in pH with an approximate threefold change for each one-half unit change in pH for a constant water temperature. Studies by Skarheim (1973) also have shown the nonionized ammonia nitrogen as a function of total dissolved solids. Because water in the study reach does not contain more than 200 to 300 mg/L total dissolved solids (Wentz and Steele, 1976), the effect of total dissolved solids was not considered in the simulation phase of the analysis.
Table k.--Percentage of nonionized ammonia nitrogen in ammonia-water solutions at various pH and water-temperature values
[After Willingham, 1976]
pH(units)
9.08.58.07.5
7.06.56.0
Water temperature,
0
7.6k2.55.820.261
.083
.026
.008
3
9.603.251.05.335
.106
.034
.011
6
11.94.111.33Ml
.135
.043
.OH
in degrees
9
14.75.161.69.542
.172
.054
.017
Celsi us
12
17.96.442.13.68
.217
.069
.022
15
21.57.972.67.859
.273
.086
.027
Effluent Standards of Proposed Regional Wastewater-Treatment Plant
A regional wastewater-treatment plant is proposed in the Steamboat Springs area and is discussed in a wastewater-facilities plan (U.S. Environmental Protection Agency, 1977). The proposed regional plant would replace existing wastewater-treatment plants at Mount Werner, Steamboat Springs, the KOA campground, Sleepy Bear Trailer Park, and the Steamboat-II development (table 1).
As part of the wastewater-facilities plan for Steamboat Springs, four alternative treatment-plant proposals are discussed (U.S. Environmental Protection Agency, 1977). The alternatives are as follows:
1. Site in Steamboat Springs using extended aeration in lagoons and advanced biological treatment.
38
2. Site in Steamboat Springs using advanced wastewater treatment.
extended aeration in lagoons and
3. Midway regional site using extended aeration biological treatment.
in lagoons and advanced
4. Midway regional site using mixed advanced wastewater treatment.
The physical location of the facilities involved in alternatives 1 and 2 would be in the vicinity of the existing main wastewater-treatment plant at Steamboat Springs. The proposed location of the facilities involved in alternatives 3 and 4 is approximately 2 mi (3.2 km) downstream from the existing treatment plant at Steamboat Springs and 0.25 mi (0.40 km) upstream from the Steamboat-II development. Alternatives 1 and 3 would involve complete year-round plant treatment; whereas, alternatives 2 and 4 would involve land treatment during the summer months and biological treatment in the winter months. This report will consider only alternative 1. Land treatment, incorporated in alternatives 2 and 4, could produce nonpoint-source runoff, which this study framework has no means to evaluate. Alternative 3 would result in virtually the same analysis as alternative 1 displaced by the distance between the proposed locations. The suggested standards for certain constituents in the effluent to be discharged from the proposed wastewater- treatment plant (U.S. Environmental Protection Agency, 1977) would be upgraded with time, as indicated in table 5. These proposed limitations on the concen trations of constituents in the effluent were considered in the analysis and will be discussed later in this report. Because the suggested standards for 1978 and 1983 are the same, only one simulation for 1978 was done in the analysis. The results presented also apply to 1983.
Table 5. Suggested effluent standards for proposed regional wastewater-treatment plant
[U.S. Environmental Protection Agency, 1977; units are in milligrams per liter]
YearConsti tuent
1978 1983 1985
BOD (5-day)
Total N
Ammonia |sjl--------
Nitrate -----------
D0
10
10
2.8 (September) 9.0 (December)
9.0
6.0
10
10
2.8 (September) <9.0 (December) 9.0
6.0
5.0
5.0
1.01.0
<4.5
6.0
Allowable stream loading (1978 and 1983) of ammonia nitrogen varies over yearly cycle. Two values are shown for 2 months used in the modeling analysis.
39
Flow Augmentation of Proposed Yamcolo Reservoir
The Upper Yampa Water Conservancy District and the Colorado Water Conser vation Board have proposed the construction of the Yamcolo Reservoir on the Bear River, a tributary of the Yampa River, in Garfield County (fig. 1). An engineering feasibility and operational study of the proposed reservoir was conducted by Western Engineers, Inc. (1975). Of particular interest to the analysis of waste-load assimilative capacity is the possibility of streamflow augmentation from the Yamcolo Reservoir in the study reach of the Yampa River. The operating rules of the proposed reservoir would include possible augmenta tion of flow in the Yampa River by approximately 20 ft 3 /s (0.6 m 3 /s) during low-flow periods. On the basis of streamflow-data evaluation for the State of Colorado (Livingston, 1970), the computed Q7,10 flow of the Yampa River at Steamboat Springs was 28 ft 3 /s (0.78 m 3 /s) with a standard deviation of 27 ft 3 /s (0.76 m 3 /s) and a standard error of 58 percent. The regionalized Q7,10 flow estimated for tributaries flowing into the study reach should be consid ered only as fair to poor in accuracy (Livingston, 1970). The augmentation from the proposed Yamcolo Reservoir could increase the Q7,10 flow in the Yampa River at Steamboat Springs to 48 ft 3 /s (1.3 m 3 /s). The option of flow augmen tation of 20 ft 3 /s (0.6 m 3 /s) was used in evaluating the effects of treated wastewater effluent from the proposed regional treatment plant during Septem ber and December flow conditions on the Yampa River.
Population Projections for Steamboat Springs
A substantial seasonal variation in population in the Steamboat Springs area occurs each year because of summer and winter recreational activities. Because of this variation, two population indexes have been used in this study: peak-day and permanent population. The peak-day population is defined as the maximum daily population for a given year, and permanent population is defined as the average number of year-round residents in the Steamboat Springs area.
Population projections and extrapolations for 2010 based on studies by Gathers and Associates (1976) as revised by the Routt County Planning Commis sion and the U.S. Environmental Protection Agency (1977) (fig. 23) were used for the model-simulation phase of this study.
The projected population for 2010 was recommended by the Routt County Planning Commission (P. E. Stark, written commun., 1976) for use in the model- simulation phase of this report. The projected peak-day population for 2010 is approximately 26,000 and the projected permanent population is 11,500. These population estimates were determined by extending the graphical rela tionships in figure 23 to 2010. An assumed water use of 134 gal (0.509 m 3 ) per capita per day was used in conjunction with the projected peak-day popula tion of 26,000 and a maximum-design capacity of 3.5 Mgal (13.1xl0 3 m 3 ) per day for the proposed wastewater-treatment plant (F. J. Vogel, written commun., 1975). As noted in the wastewater-facilities plan (U.S. Environmental
40
i
IPiP
oo'
o5f-l
r*
cf
CT O
TO
Protection Agency, 1977), there has been a large variation in the population projections for the Steamboat Springs area, both in terms of population and of the per-capita use of water. For the study reported here, the peak-daily population was used for the December modeling condition and the permanent population projection was used for the September modeling condition.
The population projections considered in this analysis (Gathers and Asso ciates, 1976; U.S. Environmental Protection Agency, 1977) are considered only as indications of future growth. Because of the uncertainty of population projections, estimates for 2010 were used for the entire analysis. This pro moted a maximum loading viewpoint of the stream reach from the proposed re gional wastewater-treatment plant (U.S. Environmental Protection Agency, 1977).
MODEL SIMULATIONS
Effluents from Existing Wastewater-Treatment Plants and 7-Day, 10-Year Streamflow
Model simulations were conducted using the U.S. Geological Survey model for effluents from existing wastewater-treatment plants in the study reach, assuming low-flow conditions and population projections for 2010. The nonaugmented Q7,10 flow of 28 ft 3 /s (0.78 m3 /s) was considered. Also con sidered was a flow augmentation of 20 ft 3 /s (0.56 m 3 /s) from the proposed Yamcolo Reservoir. The analysis assumed the same effluent conditions measured during the September 23-24, 1975, sampling period. The results of this analysis using both flow conditions for selected water-quality variables are shown on figures 24 to 27. For the simulation phase, it was assumed that all model parameters were only temperature dependent. Therefore, for example, the simulation for December basically assumed the same biologic and pH conditions as in September. The amount of error induced because of this and other assumptions is not known.
The traveltime through the system, for the simulation phase of the study only, was estimated using a procedure described by Boning (1974). This procedure utilizes a technique of estimating traveltime as a function of discharge and mean streambed slope. After the traveltime data were computed, they were used in computations of subreach transit times and, also, of the stream-reaeration rates. A subsequent study of the traveltime and reaeration rates for the Yampa River has been completed (D. P. Bauer, R. E. Rathbun, and H. W. Lowham, written commun., 1977). Traveltime determined in the subsequent study is within 5 percent of the traveltime used in this study.
Simulated profiles for concentrations of CBODU are shown on figure 24. The largest concentrations, as expected, occur immediately downstream from the main wastewater-treatment plant at Steamboat Springs. No criteria concerning in-stream concentrations of this constituent have been proposed to date.
42
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
310 I
300 I
290 I
280 I
270 I
260 I
250 I
o I i o o o o
5.0 4.0
>i 0£
3.0
"I
to z
2.
0
O UJ O
1.0
0.0
NO
NA
UG
ME
NT
ED
FLO
W
AU
GM
EN
TE
D F
LOW
I19
519
018
518
0 17
5 17
0
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
24.-
Com
pute
d c
once
ntr
atio
ns
of u
ltim
ate-
carb
onac
eous
bioc
hem
ical
oxy
gen
dem
and
for
exis
ting
was
tew
ater
tre
atm
ent
faci
liti
es a
nd Q
7, 1
0 f
low
con
diti
ons
in t
he Y
ampa
Riv
er.
As shown on figure 25, the computed DO concentrations exceed by more than 2 mg/L the existing limit of 6.0 mg/L for a type Bl stream classification. Little increase in DO concentrations resulted from substituting augmented flow conditions in the model simulation (fig. 25).
Profiles for both computed ammonia-nitrogen and nonionized ammonia-nitro gen concentrations for the two flow conditions are shown on figure 26. The model computed the ammonia-nitrogen concentrations, from which the nonionized ammonia-nitrogen concentrations then were computed using the tables compiled by Willingham (1976). An average value of 13°C for water temperature was assumed throughout the stream reach, and the median pH values for each of the main-stem sampling sites observed during the 24-hour sampling period in September 1975 were used. The results (fig. 26) show a maximum nonionized ammonia-nitrogen concentration of 0.03 mg/L for augmented flow conditions and 0.05 mg/L for nonaugmented flow conditions. Both concentrations exceed the proposed 0.02-mg/L concentration for a type Bl stream classification for a short distance downstream from the main wastewater-treatment plant at Steamboat Springs. The exact length of stream in which the nonionized ammonia-nitrogen concentration would exceed the proposed 0.02-mg/L standard cannot be determined, because of the uncertainty of the decay-rate coefficient for ammonia nitrogen and, also, because of the possible large use of ammonia nitrogen by different algal forms, discussed earlier. The results presented above also can vary depending on the actual water temperature and pH values. For example, considering an acceptable pH range from 6 to 9 and a water temperature from 0° to 20°C, the maximum nonionized ammonia-nitrogen concentration could range from 0 to 0.11 mg/L for augmented flow conditions and from 0 to 0.16 mg/L for nonaugmented flow conditions (table 4).
The computed concentrations of fecal-coliform bacteria for existing conditions are shown on figure 27. The computed concentrations are signifi cantly less than the maximum of 1,000 colonies/100 mL allowed for type Bl stream classification (Colorado Department of Health, 1974). The greatest difference in computed coliform-bacteria concentrations based on augmented versus nonaugmented flow conditions occurs upstream from the Elk River.
Proposed Regional Wastewater-Treatment Plant with 1978 Standards for Effluent
The effects of discharge from the proposed regional wastewater-treatment plant on the Yampa River were simulated using the U.S. Geological Survey model for the following conditions: Permanent population for September 2010 and Q7,10 flow of 28 ft 3 /s (0.78 m3 /s), permanent population for September 2010 and Q7,10 flow with 20 ft 3 /s (0.56 m 3 /s) flow augmentation, peak-daily popula tion for December 2010 and Q7,10 flow, and peak-daily population for December 2010 and Q7,10 flow with 20 ft 3 /s (0.56 m3 /s) flow augmentation. Treatment- plant effluent was assumed to meet the proposed standards for 1978 (table 5). The Q7,10 flow was assumed for September and December because low-flow occurs in both months in the Yampa River at Steamboat Springs (figs. 5, 6). The median pH values obtained during the 24-hour data-collection effort were assumed for both the September and December conditions. Water temperatures of
44
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
u o
u HI o x
AU
GM
EN
TED
FLO
W
196
190
186
180
176
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
166
Fig
ure
25. C
om
pute
d c
once
ntr
atio
ns
of d
isso
lved
oxy
gen
for
exis
tin
g w
aste
wat
er t
reat
men
t fa
cili
ties
and
Q7,
1
0 f
low
con
diti
ons
in t
he Y
ampa
Riv
er.
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
280
270
1 _______ _
I26
0 I25
0_
I0.
06
NO
NA
UG
ME
NT
ED
FLO
W
-
AU
GM
EN
TE
D F
LOW
0.0
5
o
0.0
4 - z* o I z
Am
mo
nia
0.0
3
* A
mm
on
ia
Pro
po
sed
accep
tab
le c
on
cen
trat
ion
fo
r n
on
ion
ized
\
amm
on
ia n
itro
ge
n=
0.0
2 m
illi
gra
m p
er l
ite
r0.
02
0.0
1
O O
O O
O
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
MIL
ES
165
160
1550
.0
Fig
ure
26
. C
om
pu
ted
co
nce
ntr
atio
ns
of a
mm
onia
nit
roge
n fo
r ex
isti
ng
was
tew
ater
tre
atm
ent
faci
liti
es a
nd Q
7, 1
0 f
low
con
diti
ons
in t
he Y
ampa
Riv
er.
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETER
S
260
310 I
30
0 I29
0 128
0 I27
0 I26
0 I2
50 1
§
200
2
O O
O 2
15
0
ui
O 2
O
O10
0 50
l\l O
N A
UG
ME
NTE
D F
LOW
AU
GM
EN
TED
FLO
W
0 195
I19
018
518
0 17
5 17
0
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
27
. C
om
pu
ted
co
nce
ntr
atio
ns
of f
ecal
co
lifo
rm b
acte
ria
for
exis
ting w
aste
wat
er t
reat
men
t fa
cili
ties
and
Q7,
10-
flow
con
diti
ons
in t
he Y
ampa
Riv
er.
13°C for September and 2°C for December were assumed and were based on average yearly harmonic-mean water temperatures for Steamboat Springs (T. D. Steele, written commun., 1977). The regional treatment-plant discharge was based on projected population and 134 gal (0.509 m3 ) per capita per day for the plant location (see page 40) 188 river-miles (303 km) upstream from the mouth.
The computed profiles for CBODU are shown on figure 28 and for DO on figure 29. The profile patterns were distinct for the September and December conditions. The effects of flow augmentation were negligible in both instances. The simulated DO profiles indicated only a slight DO sag downstream from Steamboat Springs (fig. 29). The differences between the September and December DO profiles can be attributed to the different water- temperature conditions. The sharp increase of DO on the December profile at the start of the curve is a result of an assumed 8.5-mg/L initial concentration, much below the saturated DO concentration.
The nutrient species modeled included ammonia nitrogen (fig. 30), total nitrogen (fig. 31), and orthophosphate (fig. 32). The profiles for ammonia- nitrogen and orthophosphate concentrations indicate that large differences between September and December seasonal conditions and the effects of augmented flow conditions are discernible. The ammonia-nitrogen concentrations in the Yampa River for the two periods were computed assuming the treatment-plant effluent contained 9.0 mg/L ammonia nitrogen, the proposed standard for December 1978; and 2.8 mg/L, the proposed standard for September 1978 (table 5). For each of the four conditions modeled, the maximum concentrations of ammonia nitrogen, total nitrogen, and orthophosphate occurred downstream from the proposed regional wastewater-treatment plant to the confluence of the Elk River in December. The profiles of nonionized ammonia-nitrogen concentrations for the four assumed conditions are shown on figure 33. The concentrations of nonionized ammonia nitrogen shown of the two profiles for the augmented and nonaugmented flow conditions computed for December exceed the proposed 0.02-mg/L concentration for a type Bl stream classification with maximum concentrations of 0.042 mg/L for augmented flow and 0.065 mg/L for nonaugmented flow. These concentrations will vary depending on the actual water-temperature and pH values. For example, for the acceptable pH range of 6 to 9 and a water temperature from 0° to 20°C, the nonionized ammonia-nitrogen concentration could range from 0 to 0.27 mg/L for augmented flow and from 0 to 0.42 mg/L for nonaugmented flow in December.
The assumed ammonia-nitrogen concentration of the proposed treatment- plant effluent was varied in the model to evaluate its effect on the nonionized ammonia-nitrogen concentrations in the Yampa River. The results for December are shown on figure 34. An ammonia-nitrogen concentration for the plant effluent of 3 mg/L resulted in a maximum of 0.025 mg/L of nonionized ammonia nitrogen for nonaugmented flow. The results for a flow augmentation of 20 ft 3 /s (0.56 m 3 /s) in December also are shown on figure 34. An ammonia- nitrogen concentration of 4 mg/L in the effluent resulted in a maximum of about 0.02 mg/L of nonionized ammonia nitrogen in the study reach. The analysis indicates that flow augmentation may allow an approximate 25-percent increase in the permissible ammonia-nitrogen waste loading in the study reach, relative to the Q7,10-flow statistics.
48
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETE
RS
4.0
3.5
3.0
2.5
2.0
310 I
300 I
290 L
2tO I
270 I
260 I
260 1
UJ O Z
O o o z < QU
J ii 8 cc O IL
LI s f-
D05
OQ z
1.5
2 o1.
0
0.5
0.0
195
NO
NA
UG
ME
NT
ED
FLO
W
AU
GM
EN
TE
D F
LOW
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
28
. C
om
pu
ted
conce
ntr
atio
ns
of u
ltim
ate
carb
onac
eous
bioc
hem
ical
oxy
gen
dem
and
usin
g 19
78 s
tan
dar
ds
for
effl
uen
t, Y
ampa
Riv
er.
no 3
T3 Ci-* 03
o3t/3
O >->>
a.t/3 C/3
00033
C t/33'
00
oot/3r* 033 D- 03
03
3T303
2<'05
DISSOLVED-OXYGEN CONCENTRATION. IN MILLIGRAMS PER LITER
00 <O OXI i 7 ~-H-^_ r
w._>o
w §
g.8 3 o >
O m
-I
i
AMMONIA-NITROGEN CONCENTRATION. IN MILLIGRAMS PER LITERo o p o -»
_»o to b) <c fo
oc "i n>u> o'l
no 3
o00n>3cC/l5'
00
-4oo(/>(-> 03s a-03
03
3 o03
2<*
Ol
w
Proposed regional treatment plant
oN) M
W -IO >
Om
61m c5 i3 5 n H r~ mO O
O
s
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
Ul
ho
2.0
1.6
1.2
0.8
0.4
0.0
195
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
MIL
ES
310 I
I
- -
** 0 Q.
C B a t*
"a g 5)
£ g (0 1
1
300
29
0
280
270
260
250
1 1
1 1
1 1
ill
II
NO
NA
UG
ME
NT
ED
FLO
W
AU
GM
EN
TE
D F
LOW
-
1\ ^
\1
\_
Dec
embe
r J
1
.---I
f \
\ N \
Sep
tem
ber
j , _
^\
/
L
Sep
tem
ber
^< J~
-
Ill
II
165
160
155
Fig
ure
31. C
om
pute
d c
on
cen
trat
ion
s of
tot
al n
itro
gen
usin
g 19
78
stan
dard
s fo
r ef
fluen
t, Y
ampa
Riv
er.
Ul
U)
1.5
cc
1Q
J
I «
0. O S
Z O
0.9
O
|0.6
O o. V) O £
0.3 0.0 19
5
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
310 I
300 I
290 I
280 1
270 I
260 I
250
_I
«--
i
NO
NA
UG
ME
NT
ED
FLO
W
AU
GM
EN
TE
D F
LOW
Sep
tem
ber
Sep
tem
berip
tem
oer
= |
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
32
. C
ompu
ted
conc
entr
atio
ns o
f or
thop
hosp
hate
usi
ng 1
978
stan
dard
s fo
r ef
flue
nt,
Yam
pa R
iver
.
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
0.0
7
NO
NA
UG
ME
NT
ED
FLO
W
AU
GM
EN
TE
D F
LOW
Pro
po
sed
ac
ce
pta
ble
co
nce
ntr
atio
n f
or
no
nio
niz
ed
amm
on
ia n
itro
gen
= 0
.02 m
illig
ram
per
lit
er
Sep
tem
ber
S
epte
mb
er
190
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
155
Fig
ure
33
. C
om
pu
ted
co
nce
ntr
atio
ns
of n
onio
nize
d am
mon
ia n
itro
gen
usin
g 19
78 s
tan
dar
ds
for
effl
uen
t,
Yam
pa R
iver
, as
sum
ing
9 m
illi
gram
s pe
r li
ter
(Dec
embe
r) a
nd 2
.8 m
illi
gram
s pe
r li
ter
(Sep
tem
ber)
of
am
mon
ia n
itro
gen
in e
fflu
ent
from
was
tew
ater
tre
atm
ent
pla
nt.
DIS
TA
NC
E F
RO
M M
OU
TH
. IN
KIL
OM
ETE
RS
Ln
Ui
0.6
310 I
30
029
0 [28
0 I27
0 I26
0 I25
0 1
0.5
5 Z»
Z o HI
U Z o U O o oc I Z
O
0.4
I I
I I
Q7
,10
+ 2
0 cu
bic
fee
t pe
r da
y 4
mill
igra
ms
per
liter
Q7
,10
3 m
illig
ram
s pe
r lit
er
0.06
0.3
0.1
0.05
0.04
0,03
^
_A
mm
on
ia
nit
rog
en (
Dec
emb
er)
-\
Pro
pose
d ac
cep
tab
le c
on
cen
trat
ion
fo
r n
on
ion
ized
amm
onia
nit
rog
en=
0.0
2 m
illig
ram
per
lit
er
0.0
No
nio
niz
ed a
mm
onia
ni
trog
en (
Dec
embe
r)
0.02
0.01
O O Z o
o HI a
o
195
185
180
175
170
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
165
160
1550.
0
Fig
ure
34
. C
om
pu
ted
co
nce
ntr
atio
ns
of a
mm
onia
nit
roge
n in
the
Yam
pa R
iver
, as
sum
ing
3.0
and
4.0
mil
ligr
ams
per
lite
r of
am
mon
ia i
n ef
flue
nt f
rom
was
tew
ater
tre
atm
ent
pla
nt.
The simulated profile for concentrations of fecal-coliform bacteria is shown on figure 35. The maximum concentration of fecal-coliform bacteria estimated by the model was less than 50 colonies/100 mL. Substantial increases in concentration of fecal-coliform bacteria (fig. 35) are caused by discharge from the proposed regional wastewater-treatment plant at Steamboat Springs and by discharge from the existing wastewater-treatment plant at Hayden into Dry Creek, which is tributary to the Yampa River, 159 river-miles (256 km) upstream from the mouth. Even with these increases, the concentrations of fecal-coliform bacteria are significantly less than the acceptable 1,000 colonies/100 mL.
Proposed Regional Wastewater-Treatment Plant with 1985 Standards for Effluent
The effects of effluent discharge from the proposed regional treatment plant on the Yampa River were simulated by the U.S. Geological Survey model assuming compliance with the proposed 1985 standards for effluent (table 5). All other assumptions are the same as described in the previous section.
Some of the results of the model analysis are shown on figures 36 and 37. The figures show profiles for ammonia-nitrogen (fig. 36) and total-nitrogen (fig. 37) concentrations. Assuming the more stringent 1985 standards for effluent (table 5), the computed concentrations of all constituents were less than the standards for type Bl streams. The results for nonionized ammonia- nitrogen concentrations are not shown. The concentrations were computed but were less than 0.01 mg/L, significantly less than the recommended standard of 0.02 mg/L (Colorado Department of Health, written commun., 1976; U.S. Environmental Protection Agency, 1976) for a type Bl stream.
Ammonia-Nitrogen Concentrations for Varying Conditions and 1978 Standards for Effluent
The study reach was modeled assuming 1978 standards for effluent from the regional treatment plant (table 5), Q7,10 flow in the Yampa River at Steamboat Springs, and several projected populations and effluent-discharge rates. The results of the analyses are given in tables 6 and 7. The projection alterna tives included peak-day populations for December and permanent populations for September. The population estimates are from the U.S. Environmental Protec tion Agency (1977) for 1990, 1995, 2000, 2005, and 2010. The alternative effluent-discharge rates considered were 100, 125, and 150 gal (0.38, 0.47, and 0.57 m3 ) per capita per day.
For permanent populations (table 6), the maximum concentrations of ammonia nitrogen ranged from 0.15 to 0.24 mg/L and nonionized ammonia nitrogen from 0.015 to 0.023 mg/L. Thus, the maximum nonionized ammonia-nitrogen con centrations exceeded the recommended standard of 0.02 mg/L by 15 percent. The water temperature was assumed to be 13°C throughout the study reach and pH values were assumed to be similar to those measured during September 23-24, 1975.
56
FECAL-COLIFORM BACTERIA CONCENTRATION. IN COLONIES PER 100 MILLILITERS
o £ S & S
oc
no3
T3C
H< O 03 02
< 3
ct/55'
CTQ
D-
1-1 D-
AMMONIA-NITROGEN CONCENTRATION. IN MILLIGRAMS PER LITER
croc
nO3'Oc
oOQ
C C/35"
00
oo
I/)
3 O-
O-C/3
3c n>3
3 *o pPO
<
J«
I
0001
i 8
zom n
I
Z
p _A
Ok
_^
S
en
3 O
3 feI
--"^
Proposed regionaltreatment plant
^
-
_
^
1
c/>o oo
1
p-
ob00
1
o w
1
D
^,
\
1
OCD O O
1
I
c
N
J
Decem
8
> o p c» L» ro * 3 O> O J
1
~
J_
-""
-
1
1
1
1
> Zc o0 Zs >g ^H g
O 3 n H ~ r- mO 0
O
1
>3*
wo
w-§
o is> S5 -§|
0 m n
IN MOUTH.
280 1
z2Fi
-s"
-8O
M-s
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
Ul
vo
195
1.0
flC UJ -1 flC a. %
0.8
s (9 .j i 2 2
0.6
P U z o A
.,
UJ f <
0.2
"
0.0
310 I
I
-
- -
"5 g ?
, 1
30
0 II
._
______-
a
I
29
0
280
270
260
260
1 I
I I
1II
II
f
\ v
/ \
Dec
embe
r , ,
r-'
-
\ p
ece
mber
1 C
-.1
r -
K
r
1 /
' v
Sep
tem
ber
ft ? T
~
-- *
- r
j ! y
Sep
tem
ber
NO
NA
UG
ME
NT
ED
FLO
W
~
AU
GM
EN
TE
D F
LOW
II
119
018
518
0 17
6 17
0
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
186
160
166
Fig
ure
37. C
ompu
ted
conc
entr
atio
ns o
f to
tal
nitr
ogen
usi
ng 1
985
stan
dard
s fo
r ef
flue
nt,
Yam
pa R
iver
.
Table 6. Maximum computed concentrations of ammonia nitrogen and nonionizedammonia nitrogen for permanent populations3 and using 1978 standards
for effluent ., Yampa River> Steamboat Springs to Hoyden* Colo.
Maximum concentration,
Year
19901995200020052010
19901995200020052010
19901995200020052010
Projected permanent population
8,2009,0009,6009,75010,000
8,2009,0009,6009,75010,000
8,2009,0009,6009,75010,000
Effluent (gal Ions
per capita)
100100100100100
125125125125125
150150150150150
i it in i i i i y i a
Ammon i a ni trogen
0.15.16.18.18.19
0.19.19.21.21.21
0.21.22.23.23.24
1113 \J^\ 1 1 UGI
Nonionized ammonia ni trogen
0.015.016.017.017.018
0.018.018.020.020.020
0.020.021.022.022.023
Results of this analysis for the peak-day populations indicated that the maximum concentrations of ammonia nitrogen ranged from 0.81 to 1.55 mg/L and nonionized ammonia-nitrogen concentrations ranged from 0.036 to 0.070 mg/L (table 7). All nonionized ammonia-nitrogen concentrations for peak-day populations exceeded the recommended standard of 0.02 mg/L. The water temperature was assumed to be 2°C throughout the study reach and pH values were assumed to be similar to those measured during September 23-24, 1975.
60
Table 1.--Maximum computed concentrations of ammonia nitrogen and nonionized ammonia nitrogen for peak-day populations, and using 1978 standards
for effluent, Yampa River, Steamboat Springs to Hay den, Colo.
Projected Effluent
Maximum concentration, in milligrams per liter
Year
19901995200020052010
19901995200020052010
19901995200020052010
peak-day population
17,00020,50023,00024,50026,000
17,50020,50023,00024,50026,000
17,50020,50023,00024,50026,000
(gal Ions per capita)
100100100100100
125125125125125
150150150150150
Ammon i a ni trogen
0.81.94
1.041.091.11
0.991.141.231.301.33
1.161.331.441.531.55
Nonionized ammon i a nitrogen
0.036.042.047.049.050
0.044.051.055.058.059
0.052.059.064.068.070
STREAM BIOLOGICAL CONDITIONS
A concurrent assessment of periphyton and macroinvertebrates was con ducted cooperatively by personnel of the U.S. Environmental Protection Agency and U.S. Geological Survey for the 16 main-stem locations in the study reach. Sixteen of the sampling locations were at the main-stem locations used in the 24-hour study (table 2), and one additional site was located between YM-13 and YM-14. The macroinvertebrate data are discussed in a report by Eddy (1975). The diversity of the benthic communities identified from the samples was found to be related closely to the known point sources of effluent discharging to the Yampa River from wastewater-treatment plants. The mean diversity-index values (Slack and others, 1973) generally decreased in a downstream direction from the point-source discharges in the study reach.
The diversity-index values were compared to ammonia-nitrogen and ortho- phosphate concentrations estimated during the model-calibration phase. The
61
results of this comparison are shown on figure 38. A general inverse rela tionship between the various nutrient concentrations and mean diversity of the benthic invertebrates is indicated (correlation coefficient, r=0.7 for ammonia nitrogen and r=0.6 for orthophosphate). Other water-quality constituents that were compared with the mean diversity values included: CBODU, water temperature, total phosphorus, and Kjeldahl nitrogen. The correlations of these latter constituents to .the mean diversity data were relatively low (correlation coefficients, r=0.2 to 0.4).
A comparison of the mean total number of various benthic organisms and both orthophosphate and ammonia-nitrogen concentrations is shown on figure 39. There is no apparent relationship between the nutrient concentrations and the mean total number of organisms observed along the stream. Periphyton data that were collected at the 16 main-stem locations are tabulated in table 8. Periphyton populations were sampled using plastic strips suspended in the stream channel (Slack and others, 1973, p. 88-89). The strips for the 16 main-stem sites were placed in the stream during early September 1975 and were removed approximately 4 weeks later. These data also were plotted against the benthic diversity index and total numbers of organisms. No trend or bivariate relationship was indicated.
Table 8. Periphyton data, Yampa River, Steamboat Springs to Hayden, Colo., August-September 1975
Location 1
YM-n--------VM-fl1 ------.-TPIU 1
YM- 1 --------
YM-9--------
YM-3
YM-if VM_ CYM-i>--------VM_£ _
VM "7 YM-/ ---- VM Q
YM-Q--------
YM-10 YM-1 1 -------T PI 1 | VM 1 1
YM- 1 7-------Tn \ j YM- 1 li-------Trl |tf
Periphyton ash weight
(g/m2 )
k.6123.5 k.26.1
7.6 kk1.21.05.1
4-3.2
116 0
. O
7.0rOPO
Periphyton dry weight
(g/m2 )
6.51515 6.1o
13rOpo1.71 ft1 . O
7.6
5 Q y.5
27209.5
71
Periphyton chlorophyl 1-a
(mg/m2 )
5.0 (. ?bo k.\5.0L Q
8.9927 Q3-0
5-0 3 ft3.0
1 Q1 « J
.3109-2^.3
59
Periphyton chlorophyl \-b
(mg/m2 )
0 0.0
.6
.1.6.k
.5
.0
.3
.2
.3
.2
.0
.k
.5
.5 1.2
figure 6, and table 2.
62
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETE
RS
UJ
0.4
0
0.3
5
oc
0.3
0U
l
H
ec Ul
a. J*
0.2
5
2 0 d S
0.2
0Z
* z g 1 °
-15
Ul
O 8 0
.10
0.0
5
31
0
30
01
11
_ YM
-0
*?!
.v
YM
- 01 '. '
r f
L_ VI
YM
~1
i S
J :>.
^v I/I- ''.*. 2 Y
M-3
VIU
_ >
1
-i
.jjr^
t p
f I': r.:
..-\_
l
YIV/~ I-
29
0
280
270
I I
I260
250
I I
I I
I I
I
AM
MO
NIA
-NIT
RO
GE
N C
ON
CE
NT
RA
TIO
N
OR
TH
OP
HO
SP
HA
TE
CO
NC
EN
TR
AT
ION
_
fi\
ME
AN
DIV
ER
SIT
Y (
d)
~ A
dapte
dr~
^
from
Edd
y
s1
i
YIV
« YW
i"1^
-5
pT
>
"
1
1- *
7 YR
Y
M-8
\._ I
r: ;>I \_
.fl-.9
YW
_
_ ^
_. I
YW
-10 \i _
__ 1
1- .U«
11 _ _
_ _
_
I(1975)
YM
-13
YM
-12
\ i_«;V.
___ .
'i*iJ
~ YM
-14
._/
i^
- _
4.0 to b
AGRAMS;
o oe 3
2.0
£ CO oe U
l >1.
0 Q z i
o.o
0.0
0 _
.-
19
5
190
185
180
175
170
165
160
155
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
Fig
ure
38
. M
ean
div
ersi
ty o
f be
nthi
c in
ver
tebra
tes
and
nutr
ient
conce
ntr
atio
ns,
Yam
pa R
iver
.
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
KIL
OM
ETER
S
a o
310
0.4
0 '
0.3
5
0.3
0
0.2
5
0.2
0
0.1
0
0.0
5
0.0
0
30
0 II
290 L I
280
270
I I
1
26
0 I2
50 1
1 1
MIV
IIVIV
SIV
i/-\
~i^
1 1 nv^o
ciM
I,IH
VV
«C
:IV
i n
/A 1
1 \jr*
OR
TH
OP
HO
SP
HA
TE
CO
NC
EN
TR
AT
ION
_
px|
ME
AN
NU
MB
ER
OF
OR
GA
NIS
MS
-
- - _ ^
YM
-01
YM
-0 r19
5
/'
r^- YH
j
i A-
TP
-(
"L_
YM
-4~
~P
U
* |
\ __
I
1
YM
-3
YM
-2:
iPI
-!
F" v_
nEi
3 A
dapte
d f
rom
Edd
y (1
975)
s _
___ I
^ L_J
r
- -
YM
-14
YM
-6
L~\
YM
-5
190
^
185
\. i
Yl
YM
-7 1 ,
180
-\ i\~YW
8"
""
YM
-9
1
i- M«
175
ift Y v
-y{
YIV
~\_
_
__
_Y
M-1
2
YM
-11
[x
j
- >*
170
165
160
13/"
_J
*i*
__
40
Si
30 ,
* 52
OV
J/J
rf
g II
20 o
£**
U
J°s
£LU
oa
ccS
<
10 D
=Z
SZ
ee<
UJ
UJ
Q.
S
n15
5
DIS
TA
NC
E F
RO
M M
OU
TH
, IN
MIL
ES
Fig
ure
39
. M
ean
tota
l nu
mbe
r of
org
anis
ms
per
squa
re m
eter
and
nut
rien
t co
nce
ntr
atio
ns,
Yam
pa R
iver
.
A study of the aquatic insects in the Yampa and White Rivers also has been completed (Ames, 1977). The results of the study indicated an abundance of insects that live in relatively unpolluted water. Ephemeroptera (mayfly species) was found as the predominant insect with Trichoptera (caddisfly) and Diptera (two-winged insects) next in abundance.
SUMMARY
The analysis of the waste-assimilative capacity of the Yampa River from Steamboat Springs to Hayden indicated that nonionized ammonia-nitrogen concen trations may exceed proposed stream standards on peak-population days. Based upon this analysis, concentrations of dissolved oxygen, fecal-coliform bacte ria, total nitrogen, and nitrate nitrogen would not exceed current or proposed stream standards for the projected waste loading. The computed concentrations were based on population projections for 2010 and inflow to the proposed re gional wastewater-treatment plant at Steamboat Springs of 134 gal (0.509 m 3 ) per capita per day. Both December and September streamflow conditions and different population projections were considered. For December, a peak-day population of approximately 26,000 was used and, for September, a permanent population of approximately 11,500 was considered. Critical-period streamflow conditions in the study reach were assumed for all analyses to be related to the 7-day low flow with 10-year-recurrence interval (28 ft 3 /s or 0.78 m 3 /s). The regionalized Q7,10 flow estimated for tributaries flowing into the study reach should be considered only as fair to poor in accuracy (Livingston, 1970). A 20-ft 3 /s (0.56-m 3 /s) flow augmentation from the proposed Yamcolo Reservoir also was considered. A permissible increase in waste load of ap proximately 25 percent was indicated by the modeling analysis for this amount of flow augmentation.
Model calibration of this study reach was done using two models the U.S. Geological Survey and the Pioneer-I steady-state water-quality models. The primary purpose of applying two models was to show differences in the computa tional algorithms used by the models for the corresponding model-parameter values. Modeling comparisons included concentrations of: Total nitrogen, nitrite nitrogen, ammonia nitrogen, nitrate nitrogen, orthophosphate, carbon aceous biochemical oxygen demand, dissolved oxygen, and fecal-coliform bacte ria. The basic framework of both models was developed utilizing the Streeter and Phelps (1925) oxygen-sag equation.
The calibration results for this study are considered fair. The fitting procedure basically involved manually adjusting the nine model parameters in the U.S. Geological Survey and Pioneer-I models to fit the observed data with a minimum difference between the observed and modeled concentrations. Result ant profiles of concentrations for carbonaceous biochemical oxygen demand, dissolved oxygen, fecal-coliform bacteria, and total nitrogen agreed closely for the two models. The computed concentrations of organic nitrogen, nitrite nitrogen, nitrate nitrogen, and orthophosphate agreed closely for the U.S. Geological Survey and the Pioneer-I models for subreaches between sites YM-01 and YM-8. Downstream from site YM-8, there were some differences noted in the computed values. The differences in the computed concentrations of organic
65
nitrogen, ammonia nitrogen, and orthophosphate could not be explained. The higher computed concentrations of nitrite-plus-nitrate nitrogen by the Pioneer-I model were explained by its accumulation of nitrogen in the nitrate form.
For the simulation phase of the study, the U.S. Geological Survey model was used. This choice was made because of previous experience with the model, the ease with which it can be used, and the form of the required data.
Several major factors were considered in the model simulations for ex ample, existing stream-reach classifications, suggested standards for effluent from a proposed regional wastewater-treatment plant, and Steamboat Springs area population projections. The stream reach has been classified by the State of Colorado as a cold-water-fishery secondary-body-contact type Bl (Colorado Department of Health, 1974). With this classification, there are various water-quality requirements. For the simulation phase, it was assumed that all model parameters were only temperature dependent. Therefore, for example, the simulation for December basically assumed the same biologic and pH conditions as in September. The amount of error induced because of this and other assumptions is not known.
A regional wastewater-treatment plant is proposed in the Steamboat Springs area (U.S. Environmental Protection Agency, 1977). There are four alternative proposals for the location of the plant site and type of treatment to be used. Two alternative proposals involve extended aeration and advanced wastewater treatment. The other two alternative proposals involve mixed advance waste treatment and include land treatment during a part of each year. This report considered only the alternatives which included extended aeration and advanced waste treatment. The land treatment of the other alternatives can produce nonpoint surface runoff which this study framework has no means to evaluate. Suggested effluent standards for 1978, 1983, and 1985 were considered in the analysis.
The population projections considered in this analysis (Gathers and Asso ciates, 1976; U.S. Environmental Protection Agency, 1977) are considered only as indications of future growth. Because of the uncertainty of population projections, estimates for 2010 were used for the entire analysis. This promoted a maximum loading viewpoint of the stream reach from the proposed regional wastewater-treatment plant (U.S. Environmental Protection Agency, 1977).
A concurrent assessment of periphyton and macroinvertebrates was con ducted cooperatively for the 16 main-stem locations by the U.S. Environmental Protection Agency and U.S. Geological Survey. A possible relation between the benthic-invertebrate diversity index and several nutrient concentrations was compared with the mean diversity data and resulted in relatively high correla tion coefficients of r=0.7 for ammonia nitrogen and r=0.6 for orthophosphate. Periphyton data collected at the 16 main-stem locations indicated low correla tion with the nutrient concentrations observed in the stream reach.
66
REFERENCES
Ames, E. L., 1977, Aquatic insects in two western slope rivers, Colorado:Colorado State Univ., Dept. Zoology-Entomology, unpub. M.S. thesis, June1977, 95 p.
Bauer, D. P., and Jennings, M. E., 1975, Steady-state segmented dissolved-oxygen model: Bay St. Louis, Miss., U.S. Geol. Survey Computer Contr.,available only from U.S. Dept. Commerce, Natl. Tech. Inf. Service, PB-241779, 104 p.
Bennett, J. P., and Rathbun, R. E., 1972, Reaeration in open-channel flow:U.S. Geol. Survey Prof. Paper 737, 75 p.
Boning, C. W., 1974, Generalization of stream travel rates and dispersioncharacteristics from time-of-travel measurements: U.S. Geol. SurveyJour. Research, v. 2, no. 4, p. 495-499.
Colorado Department of Health, 1974, Water quality standards and streamclassification: Colorado Dept. Health, Water Quality Control Div., 52 p.
Eddy, R. M., 1975, The effect of point source discharges on the diversity ofbenthic invertebrates of the Yampa River, Steamboat Springs to Hayden,Colo., September 1975: Denver, Colo., U.S. Environmental ProtectionAgency, Region VIII, 35 p.
Gathers and Associates, Inc., 1976, Comprehensive plan for Routt County,Colorado: Steamboat Springs, Colo., Routt County Plan. Dept., 500 p.[In press.]
Giles, T. F., and Brogden, R. E., 1977, Selected hydrologic data, Yampa Riverbasin, and parts of the White River basin, northwestern Colorado andsouth-central Wyoming: U.S. Geol. Survey Open-File Rept. 78-23. [Inpress.]
Hines, W. G., McKenzie, S. W., Rickert, D. A., and Renella, F. A., 1977,Dissolved oxygen regime of the Willamette River, Oregon, under conditionsof basinwide secondary treatment: U.S. Geol. Survey Circ. 715-1. [Inpress.]
Hines, W. G., Rickert, D. A., McKenzie, S. W., and Bennett, J. P., 1975,Formulation and use of practical models for river-quality assessment:U.S. Geol. Survey Circ. 715-B, p. B1-B13.
Jennings, M. E., and Bauer, D. P., 1976, Determination of biochemical oxygendemand parameters: Bay St. Louis, Miss., U.S. Geol. Survey ComputerContr., available only from U.S. Dept. Commerce, Natl. Tech. Inf.Service, PB-253 739, 55 p.
Jennings, M. E., and Bryant, C. T., 1973, Water quality modeling for wasteload allocation studies in Arkansas, stream dissolved oxygen andconservative minerals: U.S. Geol. Survey open-file report, 16 p.
Kittrell, F. W., 1969, A practical guide to water-quality studies of streams:U.S. Dept. Interior, Federal Water Polution Control Adm., CWR-5, 135 p.
Livingston, R. K. , 1970, Evaluation of the streamflow data program inColorado: U.S. Geol. Survey open-file report, 72 p.
Mahloch, J. R., 1973, An investigation into modeling technique for coliformorganisms in fresh water streams: Mississippi State Univ., Dept. CivilEng., 36 p.
McCall-Ellingson and Morrill, Inc., 1974, Water quality management plan forthe Green River basin: Report for the Colorado Dept. Health, WaterQuality Control Div., Denver, Colo., April 1974, 142 p.
67
National Academy of Sciences-National Academy of Engineering, 1973, Water quality criteria, 1972: U.S. Environmental Protection Agency Rept. EPA- R3-73-033, 594 p.
Rickert, D. A., and Hines, W. G., 1975, A practical framework for river- quality assessment: U.S. Geol. Survey Circ. 715-A, 17 p.
Shindala, Adnan, 1972, Mathematical modeling for water quality management in streams and estuaries: Mississippi State Univ., Dept. Civil Eng., August 1972, 62 p.
Skarheim, H. P., 1973, Tables of the fraction of ammonia in the undissociated form for pH 6 to 9, temperature 0-30°C, IDS 0-3000 mg/L, and salinity 5- 35 g/kg: California Univ., Berkeley, College of Eng., School of Public Health, Sanitary Eng. Research Lab., SERL rept. 73-5, 33 p.
Slack, K. V., Averett, R. C., Greeson, P. E., and Lipscomb, R. G. , 1973, Methods for collection and analysis of aquatic biological and microbiological samples: U.S. Geol. Survey Techniques Water-Resources Inv., book 5, chap. A4, 165 p.
Steele, T. D., Bauer, D. P., Wentz, D. A., and Warner, J. W., 1976a, An environmental assessment of impacts of coal development on the water resources of the Yampa River basin, Colorado and Wyoming Phase-I work plan: U.S. Geol. Survey Open-File Rept. 76-367, 17 p.
Steele, T. D., James, I. C., II, Bauer, D. P., and others, 1976£, An environmental assessment of impacts of coal development on the water resources of the Yampa River basin, Colorado and Wyoming Phase-II work plan: U.S. Geol. Survey Open-File Rept. 76-368, 33 p.
Streeter, H. W., and Phelps, E. B., 1925, A study of the pollution and natural purification of the Ohio River: Washington, U.S. Public Health Service Bull. 146, 75 p.
Texas Water Development Board, 1970, Simulation of water quality in streams and canals, in Program documentation and users manual: Austin, Tex., DOSAG-I, 51 p.
Thomann, R. V., 1972, System analysis and water quality management: Envi ronmental Sci. Services, 286 p.
URS/Ken R. White Company, 1975, Arkansas River basin, water quality management plan: Report for the Colorado Department of Health, Denver, Colo.
U.S. Department of the Interior, 1976, Final environmental statement, north west Colorado coal: Regional analysis and 5 site-specific analyses, December 1976, 4 v.
U.S. Environmental Protection Agency, 1976, Quality criteria for water: U.S. Environmental Protection Agency Rept. EPA-440/9-76-023, 501 p.
____1977, Steamboat Springs Regional Service Authority, 201 Wastewater facil ities plan Final environmental impact statement: Region VIII, Denver, Colo., prepared by Weiner & Assoc., Rept. EPA-908/5-77-005, August 1977, 243 p.
Velz, C. J., 1970, Applied stream sanitation: New York, John Wiley & Sons, Inc., 619 p.
Waddel, W. W., Cole, C. R., and Baca, R. G., 1973, A water quality model of the South Platte River basin Documentation report: Richland, Wash., Battelle Pacific Northwest Labs., 44 p.
68
Wentz, D. A., and Steele, T. D., 1976, Surface-water quality in the Yampa River basin, Colorado and Wyoming An area of accelerated coal development: Paper presented at Conf. on Water for Energy Devel., Eng. Found., Asilomar Conf. Grounds, Pacific Grove, Calif., Dec. 5-10, 1976, 28 p.
Western Engineers, Inc., 1975, Yamcolo Reservoir project Feasibility report: Denver, Colo., 75 p.
Willingham, Tom, 1976, Ammonia toxity: Denver, U.S. Environmental Protection Agency, Region VIII.
Willis, Robert, Anderson, D. R., and Dracup, J. A., 1975, Steady-state water quality modeling in streams: Am. Soc. Civil Engineers, Jour. Environ mental Eng. Div., v. 101, no. EE2, p. 245-258.
Wright, K. R., 1976, Sewage effluent turned to snow, provides storage, removes pollutants: Am. Soc. Civil Engineers, Civil Eng., May 1976, p. 88-89.
Wright-McLaughlin Engineers, 1974, Sewage effluent renovation in mountain meadows: Steamboat Springs, Colo., Timbers Water and Sanitation Dist., 65 p.
GPO 844- 171