EFFECTS OF INTERANNUAL VARIATION OF PRECIPITATION ON STREAM DISCHARGE FROM RHODE RIVER SUBWATERSHEDS1

ABSTRACT: The purpose of this study was to determine the relationships between precipitation at the seasonal and annual scale and water discharge per surface area for seven contiguous first - and second-order tributaries of the Rhode River, a small tidal tributary to Chesapeake Bay, Maryland, USA. The goal was to quantify the effects of a wide range of precipitation, representative of inter-annual variations in weather in this region. The discharges measured included both overland storm flows and groundwater, since the aquifers were perched on a clay aquiclude. Precipitation varied from 824 to 1684 mm/yr and area-weighted Rhode River watershed discharge varied from 130 to 669 mm/yr with an average of 332 mm/yr or 29.1 percent of average precipitation. Average annual dis. charges from three first-order watersheds were significantly lower per surface area and varied from 16.0 to 21.9 percent of precipitation. Winter season precipitation varied from 125 to 541 mm. Area-weighted Rhode River winter discharge varied from 26.3 to 230 mm with an average of 115 mm or 43.9 percent of average precipitation. Spring season precipitation varied from 124 to 510 mm and watershed discharge varied from 40.0 to 321 mm with an average of 138 mm or 46.9 percent of average precipitation. In the summer and fall seasons, watershed discharge averaged 40.6 and 40.9 mm or 13.5 and 14.3 percent of average precipitation, respectively. Except in winter, the proportion of precipitation discharged in the streams increased rapidly with increasing volume of precipitation. Stream order showed a higher correlation with volume of discharge than vegetative cover on the watershed.     

EVALUATION OF A STREAM AQUIFER ANALYSIS TEST USING ANALYTICAL SOLUTIONS AND FIELD DATA1

ABSTRACT: Considerable advancements have been made in the development of analytical solutions for predicting the effects of pumping wells on adjacent streams and rivers. However, these solutions have not been sufficiently evaluated against field data. The objective of this research is to evaluate the predictive performance of recently proposed analytical solutions for unsteady stream depletion using field data collected during a stream/aquifer analysis test at the Tamarack State Wildlife Area in eastern Colorado. Two primary stream/aquifer interactions exist at the Tamarack site: (1) between the South Platte River and the alluvial aquifer and (2) between a backwater stream and the alluvial aquifer. A pumping test is performed next to the backwater stream channel. Drawdown measured in observation wells is matched to predictions by recently proposed analytical solutions to derive estimates of aquifer and streambed parameters. These estimates are compared to documented aquifer properties and field measured streambed conductivity. The analytical solutions are capable of estimating reasonable values of both aquifer and streambed parameters with one solution capable of simultaneously estimating delayed aquifer yield and stream flow recharge. However, for long term water management, it is reasonable to use simplified analytical solutions not concerned with early‐time delayed yield effects. For this site, changes in the water level in the stream during the test and a varying water level profile at the beginning of the pumping test influence the application of the analytical solutions.     

SPACE-TIME MODELING OF VECTOR HYDROLOGIC SEQUENCES1

Stochastic modeling of vector hydrologic sequences is examined with a general class of space-time autoregressive integrated moving average (STARIMA) models. The models describe spatial and temporal autocorrelatjon, through dependent variables lagged both in space and time. The model structures incorporate a hierarchical ordering scheme to map the vector of observations into a network configuration. The neighboring structure used introduces a physical/geographical hierarchy to enable the model identification procedures to assist in determining appropriate correlative relationships. The three-stage iterative space-time model building procedure is illustrated using average monthly streamfiow data for a four-station network of the Southeastern Hydropower System.     

USE OF STAGE DATA TO CHARACTERIZE HYDROLOGIC CONDITIONS IN AN URBANIZING ENVIRONMENT1

ABSTRACT: This paper presents the results of a study on the use of continuous stage data to describe the relation between urban development and three aspects of hydrologic condition that are thought to influence stream ecosystems—overall stage variability, stream flashiness, and the duration of extreme‐stage conditions. This relation is examined using data from more than 70 watersheds in three contrasting environmental settings—the humid Northeast (the metropolitan Boston, Massachusetts, area); the very humid Southeast (the metropolitan Birmingham, Alabama, area); and the semiarid West (the metropolitan Salt Lake City, Utah, area). Results from the Birmingham and Boston studies provide evidence linking increased urbanization with stream flashiness. Fragmentation of developed land cover patches appears to ameliorate the effects of urbanization on overall variability and flashiness. There was less success in relating urbanization and streamflow conditions in the Salt Lake City study. A related investigation of six North Carolina sites with long term discharge and stage data indicated that hydrologic condition metrics developed using continuous stage data are comparable to flow based metrics, particularly for stream flashiness measures.     

CHANGES IN HYDROLOGIC RESPONSE DUE TO STREAM NETWORK EXTENSION VIA LAND DRAINAGE ACTIVITIES1

ABSTRACT: The objective of this work is to determine the effects of extension of a stream network through land drainage activities during the late 1800s on the hydrologic response of a watershed. The Mackinaw River Basin in Central Illinois was chosen as the focus and the pre‐land and post‐land drainage activity hydrologic responses were obtained through convolution of the hill slope and channel responses and compared. The hill slope response was computed using the kinematic wave model and the channel response was determined using the geomorphologic instantaneous unit hydrograph method. Our hypothesis was that the hydrologic response of the basin would exhibit the characteristic effects of settlement (i.e., increases in peak discharges and decreases in times to peak). This, indeed, is what occurred; however, the increase in peak discharges diminishes as scale increases, leaving only the decrease in times to peak. At larger scales, the dispersive effects of the longer hill slope lengths in the pre‐settlement scenario seem to balance the depressive effects of the longer path lengths in the post‐settlement scenario, thus the pre‐settlement and post‐settlement peak discharges are approximately equivalent. At small scales, the dispersion caused by the hill slope is larger in the pre‐settlement case; thus, the post‐settlement peak discharges are greater than the pre‐settlement.     

CONTINENTAL SCALE SIMULATION OF THE HYDROLOGIC BALANCE1

This paper describes the application of a continuous daily water balance model called SWAT (Soil and Water Assessment Tool) for the conterminous U.S. The local water balance is represented by four control volumes; (1) snow, (2) soil profile, (3) shallow aquifer, and (4) deep aquifer. The components of the water balance are simulated using “storage” models and readily available input parameters. All the required databases (soils, landuse, and topography) were assembled for the conterminous U.S. at 1:250,000 scale. A GIS interface was utilized to automate the assembly of the model input files from map layers and relational databases. The hydrologic balance for each soil association polygon (78,863 nationwide) was simulated without calibration for 20 years using dominant soil and land use properties. The model was validated by comparing simulated average annual runoff with long term average annual runoff from USGS stream gage records. Results indicate over 45 percent of the modeled U.S. are within 50 mm of measured, and 18 percent are within 10 mm without calibration. The model tended to under predict runoff in mountain areas due to lack of climate stations at high elevations. Given the limitations of the study, (i.e., spatial resolution of the data bases and model simplicity), the results show that the large scale hydrologic balance can be realistically simulated using a continuous water balance model.     

RUNOFF SIMULATION USING RADAR RAINFALL DATA1

ABSTRACT: Rainfall data products generated with the national network of WSR-88D radars are an important new data source provided by the National Weather Service. Radar-based data include rainfall depth on an hourly basis for grid cells that are nominally 4 km square. The availability of such data enables application of improved techniques for rainfall-runoff simulation. A simple quasi-distributed approach that applies a linear runoff transform to grid-ded rainfall excess has been developed. The approach is an adaptation of the Clark conceptual runoff model, which employs translation and linear storage. Data development for, and results of, an initial application to a 4160 km² watershed in the Midwestern U.S. are illustrated.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

ANALYSIS AND MODELING OF LONG-TERM STREAM TEMPERATURES ON THE STEAMBOAT CREEK BASIN,OREGON: IMPLICATIONS FOR LAND USE AND FISH HABITAT1

ABSTRACT: Steamboat Creek basin is an important source of timber and provides crucial spawning and rearing habitat for anadromous steelhead trout (Oncorhynchus mykiss). Because stream temperatures are near the upper limit of tolerance for the survival of juvenile steelhead, the possible long-term effect of clear-cut logging on stream temperatures was assessed. Twenty-year (1969–1989) records of summer stream temperature and flow from four tributaries and two reaches of Steamboat Creek and Boulder Creek (a nearby unlogged watershed) were analyzed. Logging records for the Steamboat Creek basin and air temperature records also were used in the analysis. A time-series model of the components of stream temperature (seasonal cycle of solar radiation, air temperature, streamflow, an autoregressive term of order 1, and a linear trend variable) was fitted to the water-temperature data. The linear trend variable was significant in all the fitted models except Bend Creek (a tributary fed by cool ground-water discharge) and Boulder Creek. Because no trends in either climate (i.e., air temperature) or streamflow were found in the data, the trend variable was associated with the pre-1969 loss and subsequent regrowth of riparian vegetation and shading canopies.     

LARGE AREA HYDROLOGIC MODELING AND ASSESSMENT PART I: MODEL DEVELOPMENT1

ABSTRACT: A conceptual, continuous time model called SWAT (Soil and Water Assessment Tool) was developed to assist water resource managers in assessing the impact of management on water supplies and nonpoint source pollution in watersheds and large river basins. The model is currently being utilized in several large area projects by EPA, NOAA, NRCS and others to estimate the off-site impacts of climate and management on water use, non-point source loadings, and pesticide contamination. Model development, operation, limitations, and assumptions are discussed and components of the model are described. In Part II, a GIS input/output interface is presented along with model validation on three basins within the Upper Trinity basin in Texas.     

AN OBJECTWE TEST FOR HYDROLOGIC SCALE1

ABSTRACT: Improving the reliability of parametric hydrologic models (sometimes called cenceptual rainfall-runoff models) in the continuous simulation of runoff from ungaged catchments has been frustrated by difficulties in estimating model parameters from catchment characteristics. An underlying problem is that these models use parameters to represent catchments as a whole, whereas data on catchment characteristics are collected at multiple field locations and are difficult to transform into one measure of collective impact. Subdividing the catchment and calibrating a stochastic parametric model to estimate distributions for the parameters that covered the range of observed streamflow values was found to improve the simulations. This paper presents an optimization of the amount of subdivision to use in simulation with a version of the Stanford Watershed Model using available climatological data. The calibration process assumes that catchment heterogeneity introduces errors that can be reduced by calibrating parameters as spatial distributions rather than single values. Calibrations for three diverse small gaged catchments located in California and in Virginia found the optimal number of subdivisions to range from 4 to 25 and the optimal scale to range from 0.3 to 2.1 mi².     

REGIONAL SCALE MODELING OF SURFACE AND GROUND WATER INTERACTION IN THE SNAKE RIVER BASIN1

ABSTRACT: Changes in irrigation and land use may impact discharge of the Snake River Plain aquifer, which is a major contributor to flow of the Snake River in southern Idaho. The Snake River Basin planning and management model (SRBM) has been expanded to include the spatial distribution and temporal attenuation that occurs as aquifer stresses propagate through the aquifer to the river. The SRBM is a network flow model in which aquifer characteristics have been introduced through a matrix of response functions. The response functions were determined by independently simulating the effect of a unit stress in each cell of a finite difference groundwater flow model on six reaches of the Snake River. Cells were aggregated into 20 aquifer zones and average response functions for each river reach were included in the SRBM. This approach links many of the capabilities of surface and ground water flow models. Evaluation of an artificial recharge scenario approximately reproduced estimates made by direct simulation in a ground water flow model. The example demonstrated that the method can produce reasonable results but interpretation of the results can be biased if the simulation period is not of adequate duration.     

MODEL ACCURACY IN SNOWMELT-RUNOFF FORECASTS EXTENDING FROM 1 TO 20 DAYS1

ABSTRACT: This paper examines the performance of snowmelt-runoff models in conditions approximating real-time forecast situations. These tests are one part of an intercomparison of models recently conducted by the World Meteorological Organization (WMO). Daily runoff from the Canadian snowmelt basin Illecille. waet (1155 km², 509–3150 m a.s.l.) was forecast for 1 to 20 days ahead. The performance of models was better than in a previous WMO project, which dealt with runoff simulations from historical data, for the following reasons: (1) conditions for models were more favorable than a real-time forecast situation because measured input data and not meteorological forecast inputs were distributed to the modelers; (2) the selected test basin was relatively easy to handle and familiar from the previous WMO project; and (3) all kinds of updating were allowed so that some models even improved their accuracy towards longer forecast times. Based on this experience, a more realistic follow-up project can be imagined which would include temperature forecasts and quantitative precipitation forecasts instead of measured data.     

EVALUATION OF METHODS FOR ESTIMATING STREAM WATER QUALITY PARAMETERS IN A TRANSIENT MODEL FROM STOCHASTIC DATA1

ABSTRACT: In this study the estimation of parameters in water quality models represented by linear first order partial differential equations is investigated. Two sets of simulated input-output data, one with input noise and the other with output measurement error, were used. The parameters were estimated by a gradient technique (Bard's method) and a pattern search technique. The results indicate that the output measurement error significantly affects the values of parameter estimates as compared to the noise added to the input. Bard's method consistently gave results with a smaller sum of square value.     

CHARACTERIZATION OF HYDROLOGIC CONDITIONS TO SUPPORT PLATTE RIVER SPECIES RECOVERY EFFORTS1

ABSTRACT: Efforts are under way to recover habitat for several threatened and endangered species in and along the Platte River in central Nebraska. A proposed recovery program for these species requires a means of characterizing “wet” versus “normal” versus “dry” hydrologic conditions in order to set corresponding Platte River instream flow targets. Methods of characterizing hydrologic conditions in real time were investigated for this purpose. Initially, 10 watershed variables were identified as potentially valuable indicators of hydrologic conditions. Ultimately, six multiple linear regression equations were developed for six periods of the year using a subset of these variables expressed as frequencies of nonexceedence. The adequacy of these equations for characterizing conditions was assessed by evaluating their historic correlation to subsequent flow in the central Platte River (1947–1994). These equations explained 54 to 82 percent of variability in the observed flow exceedences in the validation datasets, depending upon the period of year evaluated. These equations will provide initial criteria for setting applicable flow targets to determine, in real time, whether water regulation projects associated with the species recovery effort can divert or store flows without conflicting with recovery objectives.     

USING FIELD SCALE MODELS TO PREDICT PEAK FLOWS ON AGRICULTURAL WATERSHEDS1

The goal of this study was to develop a methodology for generating storm hydrographs at a watershed scale based on daily runoff estimates from a field scale model. The methodology was evaluated on a small agricultural watershed using the ADAPT field scale process model. A comparison of observed and predicted peak flows for 11 of the largest events that occurred in a three year period gave r² values of 0.84, 0.82, and 0.81 when the watershed was subdivided into 1, 5, and 10 sub watersheds. However, all other statistical measures improved when the watershed was subdivided into at least five sub watersheds. Guidelines need to be developed on the use of the procedure but it first needs to be evaluated on several watersheds that exhibit a range in sizes, land uses, slopes, and soil properties.     

SIMULATION OF FRESHWATER DISCHARGES FROM UNGAGED AREAS TO THE SEBASTIAN RIVER,FLORIDA1

ABSTRACT: Due to alterations in the natural drainage system over the past several decades and intensified agricultural practices, freshwater discharges to the Sebastian River, Florida, have increased substantially. As a result, salinity patterns in the Sebastian River and adjacent Indian River lagoon have been disrupted and the influx of nutrients has increased. Recently, the St. Johns River Water Management District has developed a 3‐D hydrodynamic and salinity model for the Sebastian River and adjacent Indian River to study the effects of freshwater inflows, and to set guidelines for management of future freshwater discharges. Freshwater inflows to the Sebastian River are part of the input data of the hydrodynamic model. Except for the downstream drainage areas, inflows are gaged, and the data were used for calibration of the hydrologic simulations. Collectively, the downstream ungaged areas constitute about 16 percent of the total drainage area. Because of the significant contribution to the total drainage area, reliable estimates of freshwater discharges from the ungaged areas to the Sebastian River are needed. This case study illustrates the development of a set of model parameters, reflecting the hydrologic and physiographic characteristics of the entire region. In this context region applies to the watersheds located in the coastal area along the Indian River from Titusville in the north to Vero Beach in the south. The parameter set was first tested on a number of gaged drainage basins in the region, and was then applied to the ungaged areas.     

REGULATED RIVER MODELING FOR CLIMATE CHANGE IMPACT ASSESSMENT: THE MISSOURI RIVER1

ABSTRACT: The Great Plains of the United States, drained primanly by the Missouri River, are very sensitive to shifts in climate. The six main stem dams on the Missouri River control more than one‐half of the nearly 1.5 million square kilometer basin and can store three times the annual inflow from upstream. The dams are operated by the U.S. Army Corps of Engineers using a Master Manual that describes system priorities and benefits. The complex operational rules were incorporated into the Soil and Water Assessment Tool computer model (SWAT). SWAT is a distributed parameter rainfall‐runoff model capable of simulating the transpiration suppression effects of CO₂ enrichment. The new reservoir algorithms were calibrated using a 25‐year long historic record of basin climate and discharge records. Results demonstrate that it is possible to incorporate the operation of a highly regulated river system into a complex rainfall‐runoff model. The algorithms were then tested using extreme climate scenarios indicative of a prolonged drought, a short drought, and a ten percent increase in basin‐wide precipitation. It is apparent that the rules for operating the reservoirs will likely require modification if, for example, upper‐basin precipitation were to increase only ten percent under changed climate conditions.