<br />
<br /> FL Ogdell el al./ JowrMl of Hydrology 228 (2000) 82-]00 87
<br /> 120
<br /> 0
<br /> 100
<br />E
<br />S 80 ,
<br />.. "
<br />'E 0
<br />'.
<br />a: 80
<br />~ 0
<br />-
<br />, 0
<br />0
<br />J:
<br />" 40
<br />~
<br />.
<br />a:
<br /> 20
<br /> + +
<br /> + +
<br /> 0 100 120
<br /> 0 20 40 80 80
<br /> Gage Hourly Rainfall (mm)
<br />
<br />86
<br />
<br />FL Ogdell el a/./ JOWrMl of Hydrology 228 (2000) 82-100
<br />
<br />Bras (1995); BJoschl et al (1995), and Seyfried and
<br />Wilcox (1995).
<br />Horsetooth Reservoir, shown in Fig. 2. provides
<br />irrigation storage and serves as the water supply for
<br />the City of Fort Collins. II is a permanent water
<br />storage facility that receives most of its supply from
<br />trans-basins transfers from the upper Colorado river
<br />via the Colorado-Big Thompson (CSn diversion
<br />project. Horsetooth lake is formed by four dams on
<br />a ridge in the foothills of the Rocky Mountains. The
<br />sUIface area of the lake on 28 July 1997 was 7.4lcm2,
<br />with a contributing area of 37.4 km~. This mountai-
<br />nous watershed contains six different soil types, five
<br />of which are common to the Spring Creek watershed.
<br />As shown in Fig. 2. the Horsetooth watershed is adja-
<br />cent to the Spring Cre~k and was affected by the same
<br />stonns as Spring Creek on 27-28 July 1997, Interest-
<br />ingly, when the Horsetooth Reservoir was constructed
<br />in the early 19505, the Spring Creek watershed was
<br />hydraulically cut in half.
<br />Soil texture maps for both Spring Creek and Horse-
<br />tooth watersheds were obtained and digitized. Flood
<br />plain maps at I :100 scale. developed by others as pan
<br />of the Spring Creek master drainage plan, were used
<br />to develop channel cross-sections and 10 obtain details
<br />of the major detention basins. Land-use data reveal
<br />that 23% of the Spring Creek watershed area is
<br />covered by impervious surfaces. The location of the
<br />only stage-recording station on the main channel of
<br />the Spring Creek was considered the watershed outlet.
<br />UnfoI1unately, this slage recorder was damaged by
<br />floodwaters, and no data were collected after 22:30
<br />MDT on 28 July 1997.
<br />
<br />7, Radar-rainfaU
<br />
<br />The spatial and temporal resolution of rainfall dala
<br />has a great impacl on the accuracy of physics-based
<br />distributed hydrologic models as verified by many
<br />researchers. Woolhiser (1996) adeptly stated that
<br />expecting a physically based hydrologic model to
<br />correctly predict storm runoff using erroneous rainfall
<br />data is akin to expecting a beam flexure equation to
<br />predict the correct beam displacement under an erro-
<br />neous loading. The magnitude of runoff prediction
<br />errors resulting from rainfall spatial and temporal
<br />resolution errors depends on various factors due to
<br />
<br />the complexity of the flood producing system, i.e.
<br />the interaction of the watershed topography and
<br />geomorphology, the runoff production mechanism,
<br />the ratio of the rainfall rate to the saturated hydraulic
<br />conductivity of the soil (Saghafian et al., 1995) and
<br />the spatial_and temporal variability of the rainfall
<br />field. This wide range of factors might be the reason
<br />for varied conclusions among researchers who sfudjed
<br />the impacts of spatial and temporal variability on
<br />runoff. Although some researchers have downplayed
<br />lhe importance of averaging in certain situations,
<br />many researchers have shown the opposite in other
<br />studies. Ogden and Julien (1994) identified twO
<br />diSlinct causes of error due to precipitation spatial
<br />aggregation, "storm smearing" and "watershed smear-
<br />ing". Stann smearing occurs because spatial aggrega-
<br />tion tends to reduce areal-average rainfall rates.
<br />Watershed smearing causes uncertainty in the location
<br />of precipitation over the catchment due to aggrega-
<br />tion. Winchel et al. (1998) showed that spatial and
<br />temporal averaging of rainfall affects modeling results
<br />more for Hortonian runoff production than for satura-
<br />tion excess runoff production. Radars and rain gages
<br />have very different sampling characteristics in both
<br />space and time. As an example, it is possible that
<br />instantaneous radar measurements taken every
<br />30 min can produce more erroneous results than lhe
<br />use of hourly integrated averages of rainfall rates
<br />(Sharif and Ogden, 1998).
<br />Only four recording rain gages cover lhe Spring
<br />Creek watershed area when using Thiessen poly-
<br />gons to determine the area of influence of each
<br />rain gage, Two of lhese four gages cover 90%
<br />of the watershed. Furthermore, extremely high
<br />rainfall intensity gradients were observed on 28
<br />July 1997. Given the likely inadequacy of the
<br />rain gage network, investigation of the accuracy of
<br />radar-rainfall estimates was necessary for this study
<br />(Pessoa et al.. 1993). Both available single- and dual-
<br />polarization rainfall rate estimates were verified using
<br />rain gage measurements.
<br />
<br />
<br />
<br />Fi;. 3. Radar-raiofall estimates vs. rain galle observations for uncalibJaled polarimetric and WSR-88D radu.rainfall algorithms.
<br />
<br />radar located near Cheyenne, WY, which is 85 km
<br />north-northeast of Fort Collins. The WSR-88D is a
<br />single-polarization radar which observes the reflec-
<br />tivity factor, Z, which is theoretically given by:
<br />
<br />Z= V-I iD~
<br />1'=1
<br />
<br />where V is the volume of the sample element (cubic
<br />meters), 11 is the number of raindrops in the sample
<br />element, and Dj is the diameter (in millimeters) of the
<br />jth drop. The reflectivity factor is converted to rainfall
<br />rate R through an empirical R-Z power function of the
<br />form:
<br />
<br />R=aZb
<br />
<br />The R-Z parameters used are a = 0.017 and b =
<br />0.714 for R in mmlh and Z in mmt-/m3. These para-
<br />meters correspond to the relation Z = 300R1.4 (Batlan,
<br />1973) which was in use at the Cheyenne. Wyoming,
<br />WSR-88D site on 28 July 1997. WSR-88D openlfors
<br />have the option of using a tropical rainfall l-R rela-
<br />tion (Fulton et a1. 1998), however, this relation was
<br />not in use during !he Fort Collins stonn.
<br />The operational WSR-88D rainfall products (Smith
<br />et al., 1996) underestimated the total rainfall for this
<br />storm by a factor of 2.0 compared to rain gage
<br />accumulations, Hourly WSR-88D radar-rainfall esti-
<br />mates compared with hourly measurements at the
<br />
<br />7.1. WSR-88D raillfa" estimates
<br />
<br />The precipitation field was estimated from the
<br />three-dimensional volume scan reflectivity fields
<br />observed by the S-band US National Weather Service
<br />WSR-88D (Weather Surveillance Radar, 1988-Doppler)
<br />
<br />(I)
<br />
<br />seven rain gages nearest to Spring Creek are shown
<br />on Fig. 3. WSR-g8D underestimation of rainfall can
<br />be attributed to the distinctive doud microphysics of
<br />extreme storms, e.g. anomalous drop size distribution,
<br />the vertical structure of the storm (Smith et aL. 1996),
<br />abnonnal propagation of the radar beam (ground
<br />clutter) and radar calibration error. Petersen et al.
<br />(1999) present a detailed discussion of radar-rainfall
<br />estimation errOrs for the Fort Collins storm.
<br />
<br />7.2. Polarimetric radar-rainfall estimates
<br />
<br />(2)
<br />
<br />The polarimetric weather radar observations used
<br />in this study were recorded by the CSU-CHILL
<br />research radar facility located near Greeley, Colorado,
<br />which is approximately 35 Ian east-southeast of Fort
<br />CoHins. The CSU-CHILL S-band radar is dual, line-
<br />arly polarized with one-degree beam width. The radar
<br />measures the reflectivity in twO orthogonal (horizontal
<br />and vertical) polarizations, ZIl and lv. respectively.
<br />When the two reflectivities are measured in an
<br />approximately simultaneous fashion, the differential
<br />reflectivity (in decibels) can be derived by lOR =
<br />10 10glO(ZH/Zv), Differential reflectivity ZOR provides
<br />an indicator of raindrop oblateness and thus raindrop
<br />size (Seliga and Bringi, 1976). For an assumed rain-
<br />drop size distribution (RSD). lOR provides an estimate
<br />of one of RSD parameters, and when combined with
<br />ZH results in a more accurate estimate of rainfall rate
<br />
|