DESCRIPTION OF THE MODEL 7
<br />random variables. The probability distribution assigned
<br />to a random variable represents the best-fit distribution
<br />(one of 17 possible probability distributions) as deter- i
<br />mined using the Chi-square goodness-of-fit criteria. The
<br />addition of future measurements may modify selection of ',
<br />the probability distribution and associated statistics. ',
<br />2. Measurement errors associated with random variables
<br />are second-order (less than 10 percent of the maximum)
<br />as compared to daily variability.
<br />Numerical simulation using random variables captures
<br />temporal variability thereby providing a realistic
<br />assessment of the effect reservoir operations have on
<br />streamflow quantity and quality in the Colorado River. ~,
<br />4. The stability and convergence of a prediction variable is
<br />accurate to within 1 percent when using the Monte Carlo
<br />approach. Using the bootstrap approach (Werckman and ',
<br />others, 2001), the minimum number of simulations
<br />required to achieve this level of stability and convergence ',
<br />was determined to be about 1,500.
<br />Measurements
<br />Measurements used in developing the stochastic mixing ',
<br />model were obtained from various sources including the
<br />USGS's National Water Information System (NWIS) (Kiner-
<br />ney, 2001); State of Colorado, Division of Water Resources
<br />HydroBase (Colorado Division of Water Resources, 2002), and ',
<br />Colorado State University, Orchard Mesa database (Harold
<br />Larsen; personal commun., 2002). For example, streamflow
<br />and dissolved-solids measurements were retrieved from the
<br />USGS NWIS database for the Colorado River streamflow '
<br />gaging stations near De Beque (09093700), Cameo (09095500), '
<br />and below the Grand Valley diversion near Palisade
<br />(09106000); Plateau Creek near Cameo (09105000); and Dry '
<br />Fork near De Beque, Colorado (09095400). Daily streamflow
<br />measurements were available fora 27-year period (1974-
<br />2001), whereas values of the dissolved solids and associated
<br />instantaneous streamflow were available at an approximately ',
<br />monthly frequency and a fewer number of years: Because i
<br />streamflow is regulated by upstream reservoirs, streamflow in
<br />the Upper Colorado River Basin system is seasonally periodic,
<br />and correlation (dependency) between annual streamflow
<br />records generally exceeds 0.70. A matrix of correlation coeffi-
<br />dents between total annual streamflow (integrated on a daily '~
<br />basis) time series for the Colorado station near Cameo I
<br />(09095500) is provided in table 1.
<br />Despite the persistent seasonal periodicities, correlations '',
<br />between annual streamflow records for different years are less ''
<br />than perfect (correlation coefficient less than 1.0) due to annual '~,
<br />variations in hydrologic conditions. Identifying the correspond-',
<br />ing hydrologic conditions is done by using a cumulative distri-.
<br />bution function (also called flow duration or exceedance curve) ',
<br />for annual flow (fig. 3). In the annual streamflow-exceedance
<br />curve, the 27-year period contains a range of hydrologic condi-'.
<br />lions that include 7 dry, 14 average, and 6 wet periods. The '
<br />average hydrologic period is defined as that period where flow
<br />exceedances are in the 25 to 75 percent range, whereas the
<br />wet and dry hydrologic periods are defined as those respective
<br />annual flows that are in the 0 to 25 percent and 75 to
<br />100 percent exceedance range.
<br />Discrete annual cumulative streamflow values, water year,
<br />rank, exceedance value, and corresponding hydrologic condi-
<br />tions are summarized by calendar year for the station near
<br />Cameo (09095500) in table 2. Based on information gathered at
<br />the time of this study, the annual streamflow of 657,210 ft3/s
<br />during 1977 was the minimum on record indicating a dry
<br />period, whereas streamflow of 2,859,690 ft3/s during 1984 was
<br />the maximum on record indicating a wet period. The 2002
<br />drought resulted in a slightly smaller (about 8-percent) annual
<br />minimum streamflow record of 604,026 ft3/s compared to 1977.
<br />Selection of this streamflow record provided a wide variety and
<br />balanced number of hydrologic conditions so as to minimize
<br />any potential model bias. At time scales on the order of days or
<br />weeks, streamflow is highly variable and can be characterized
<br />by using probability distributions.
<br />Daily Colorado River diversion records for water years
<br />1974 to 2001 were obtained for diverters downstream from
<br />De Beque and upstream from Palisade using the Colorado State
<br />HydroBase (Colorado Division of Water Resources, 2002). Pan
<br />evaporation measurements collected at the Orchard Mesa
<br />Research Station (about 15 mi west from the study area) were
<br />obtained from Harold Larsen (Colorado State University, writ-
<br />ten commun., 2002). Because evaporation measurements were
<br />limited to the growing season, a predictive equation was fit by
<br />using nonlinear regression so that the usefulness of these data
<br />could be extended to a full year. Because of the small number
<br />of salinity and evaporation measurements, these data also were
<br />fit to predictive equations by using nonlinear regression, and the
<br />probability distributions associated with residuals (difference
<br />between measurement and prediction) were used to generate
<br />random values for use in stochastic modeling. Nonlinear regres-
<br />sion also was used to fit functional relations for annual stream-
<br />flow exceedance probability and reservoir surface area as a
<br />function of reservoir volume.
<br />Instantaneous streamflow and salinity measurements were
<br />obtained from the USGS NWIS database for the Colorado River
<br />De Beque, Cameo, and Palisade and Plateau Creek stations, for
<br />Dry Fork and Upper Dry Fork stations near De Beque, and for
<br />a USGS National Water Quality Assessment (Spahr and others,
<br />1996) site on the Dry Fork at Upper Station near De Beque.
<br />Hydrograph separation was performed on the Dry Fork dis-
<br />charge measurements using published USGS software (Rut-
<br />ledge, 1998). By using the computed baseflow (low-flow)
<br />hydrograph, the streamflow hydrograph was adjusted to obtain
<br />a unit runoff hydrograph that was applied to the Sulphur Gulch
<br />drainage. Two water samples also were collected from a site
<br />(391607108153500) at the mouth of the Sulphur Gulch water-
<br />shed and analyzed for salinity and selenium during the spring
<br />2002 runoff. Instantaneous streamflow and basic parameters
<br />such as dissolved oxygen, alkalinity, pH, and temperature were
<br />measured concurrently with the sample collection.
<br />
|