SOIL MOISTURE ACCOUNTING COMPONENT OF THE USDAHL-74 MODEL OF WATERSHED HYDROLOGY1

ABSTRACT: Soil moisture in two layers of a soil near Chickasha, Oklahoma, was simulated, using USDAHL-74 Model of Watershed Hydrology. Weekly values computed for both layers compared well with those observed during the 15-month period. Certain key parameters required adjustments in the model which illustrate the need for accurate input information. The experiment demonstrates that the model, which has previously given good results in continuous streamflow prediction on watersheds up to 100 square miles, can also compute soil moisture continuously at a site. This capability suggests other model uses, for example, in monitoring the disposition of applied chemicals.     

HYDROLOGY OF A SMALL IMPOUNDMENT AND WATERSHED IN CENTRAL OKLAHOMA1

ABSTRACT: West Bitter Creek floodwater retarding structure site 3 in South Central Oklahoma was instrumented and records obtained and analyzed to obtain information concerning an impoundment water budget that is useful to landowners and designers of these impoundments. On-site loss of water from the impoundment was only 17 percent of the inflow during three years when the annual precipitation averaged 26 inches and the annual inflow averaged 1.4 inches. Runoff from an eroded area with no farm ponds was about 70 percent greater per unit area than from a portion of the watershed where 71 percent of the drainage area was controlled by farm ponds. A previous study indicated, however, that the ponds were reducing runoff only 13 percent. Loss of top soil increases runoff considerably. Only 24 percent of the total runoff into the impoundment was base flow. The flow rate into the impoundment was less than 0.05 cfs 70 percent of the time, and the inflow rate exceeded 10 cfs only 1 percent of the time. SCS runoff curve numbers varied between 57 and 96 for the impoundment watershed with an inverse relation between precipitation amount and curve number apprently caused by partial area runoff from impervious and semi-impervious areas. A comparison of measured event runoff versus event runoff computed by the SCS curve numbers gave an r² of only 0.44. However, the total computed surface runoff for eight years of record was less than 1 percent below the measured runoff which indicated the curve number method was a good tool for predicting long term runoff for the watershed.     

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 II: MODEL APPLICATION1

ABSTRACT: This paper describes the application of a river basin scale hydrologic model (described in Part I) to Richland and Chambers Creeks watershed (RC watershed) in upper Trinity River basin in Texas. The inputs to the model were accumulated from hydro-graphic and geographic databases and maps using a raster-based GIS. Available weather data from 12 weather stations in and around the watershed and stream flow data from two USGS stream gauge station for the period 1965 to 1984 were used in the flow calibration and validation. Sediment calibration was carried out for the period 1988 through 1994 using the 1994 sediment survey data from the Richland-Chambers lake. Sediment validation was conducted on a subwatershed (Mill Creek watershed) situated on Chambers Creek of the RC watershed. The model was evaluated by well established statistical and visual methods and was found to explain at least 84 percent and 65 percent of the variability in the observed stream flow data for the calibration and validation periods, respectively. In addition, the model predicted the accumulated sediment load within 2 percent and 9 percent from the observed data for the RC watershed and Mill Creek watershed, respectively.     

EFFECTS OF BASIN-SCALE TIMBER HARVEST ON WATER YIELD AND PEAK STREAMFLOW1

ABSTRACT: Streamflow changes resulting from clearcut harvest of lodgepole pine (Pinus contorta) on a 2145 hectare drainage basin are evaluated by the paired watershed technique. Thirty years of continuous daily streamflow records were used in the analysis, including 10 pre-harvest and 20 post-harvest years of data. Regression analysis was used to estimate the effects of timber harvest on annual water yield and annual peak discharge. Removal of 14 million board feet of lodgepole pine (Pinus contorta) from about 526 hectares (25 percent of the basin) produced an average of 14.7 cm additional water yield per year, or an increase of 52 percent. Mean annual daily maximum discharge also increased by 1.6 cubic meters per second or 66 percent. Increases occurred primarily during the period of May through August with little or no change in wintertime streamflows. Results suggest that clearcutting conifers in relatively large watersheds (> 2000 ha) may produce significant increases in water yield and flooding. Implications of altered streamflow regimes are important for assessing the future ecological integrity of stream ecosystems subject to large-scale timber harvest and other disturbances that remove a substantial proportion of the forest cover.     

MODELING THE FOREST HYDROLOGY OF WETLAND-UPLAND ECOSYSTEMS IN FLORIDA1

ABSTRACT: Few hydrological models are applicable to pine flat-woods which are a mosaic of pine plantations and cypress swamps. Unique features of this system include ephemeral sheet flow, shallow dynamic ground water table, high rainfall and evapotranspiration, and high infiltration rates. A FLATWOODS model has been developed specifically for the cypress wetland-pine upland landscape by integrating a 2-D ground water model, a Variable-Source-Area (VAS)-based surface flow model, an evapotranspiration (ET) model, and an unsaturated water flow model. The FLATWOODS model utilizes a distributed approach by dividing the entire simulation domain into regular cells. It has the capability to continuously simulate the daily values of ground water table depth, ET, and soil moisture content distributions in a watershed. The model has been calibrated and validated with a 15-year runoff and a four-year ground water table data set from two different pine flat woods research watersheds in northern Florida. This model may be used for predicting hydrologic impacts of different forest management practices in the coastal regions.     

FORMAL EDUCATION IN WILDLAND HYDROLOGY AND WATERSHED MANAGEMENT IN NORTH AMERICA1

ABSTRACT: A survey was made to determine the status of formal education in wildland hydrology by colleges and universities in the United States, Canada, and Mexico. As of December 1977 nine institutions offered the B.S. degree, 18 the M.S. degree, and 17 the degree of Ph.D. with a major or minor in watershed management, forest hydrology, or range hydrology. In addition, 8 other schools offer a minor in watershed management. The survey indicated 44 schools in the United States offer a total of 157 courses, five schools in Canada offer 24 courses, and 1 in Mexico offers one course in the related areas. The survey illustrated rapid growth in education programs and it is anticipated that growth will continue.     

SAGEBRUSH RANGELAND HYDROLOGY AND EVALUATION OF THE SPUR HYDROLOGY MODEL1

ABSTRACT: An excellent hydrologic record on sagebrush range-land has been developed at the Reynolds Creek Experimental Watershed in southwestern Idaho. The objectives of this paper were two-fold: (1) to analyze and describe the hydrologic record (8–18 years) from four sagebrush watersheds (1–83 ha); and (2) to evaluate the hydrology component of SPUR, a comprehensive rangeland model. The watersheds represent a gradient in elevation (1180–1658 m) and precipitation (240–350 mm/yr). Runoff was a small fraction (> 2 percent) of the total water budget for all of the watersheds. It occurred very infrequently at the three lower elevation watersheds: Summit, Flats, and Nancy Gulch. At Lower Sheep, the highest elevation watershed, runoff occurred most years for a period of 1 to 17 weeks in the winter. Frozen soil combined with rainfall or snowmelt was associated with most of the runoff from Flats and Nancy Gulch. At Summit summertime thunderstorms produced all of the runoff. The average annual sediment yield from all of the watersheds was low (17–950 kg/ha). It was highest from Summit, which had well developed alluvial channels and very steep slopes. SPUR was able to simulate runoff with reasonable accuracy only at Summit, where frozen soils were not a factor. There was poor correlation between predicted and actual annual 8ediment loss. The model tended to overpredict evapotranspiration early in the growing season and underpredict it in the late summer.     

APPLICATION OF THE RHESSys MODEL TO A CALIFORNIA SEMIARID SHRUBLAND WATERSHED1

ABSTRACT: Distributed hydrologic models which link seasonal streamflow and soil moisture patterns with spatial patterns of vegetation are important tools for understanding the sensitivity of Mediterranean type ecosystems to future climate and land use change. RHESSys (Regional Hydro‐Ecologic Simulation System) is a coupled spatially distributed hydroecological model that is designed to be able to represent these feedbacks between hydrologic and vegetation carbon and nutrient cycling processes. However, RHESSys has not previously been applied to semiarid shrubland watersheds. In this study, the hydrologic submodel of RHESSys is evaluated by comparing model predictions of monthly and annual streamflow to stream gage data and by comparing RHESSys behavior to that of another hydrologic model of similar complexity, MIKESHE, for a 34 km² watershed near Santa Barbara, California. In model intercomparison, the differences in predictions of temporal patterns in streamflow, sensitivity of model predictions to calibration parameters and landscape representation, and differences in model estimates of soil moisture patterns are explored. Results from this study show that both models adequately predict seasonal patterns of streamflow response relative to observed data, but differ significantly in terms of estimates of soil moisture patterns and sensitivity of those patterns to the scale of landscape tessellation used to derive spatially distributed elements. This sensitivity has implications for implementing RHESSys as a tool to investigate interactions between hydrology and ecosystem processes.     

GIS‐BASED HYIROLOGIC MODELING OF RIPARIAN AREAS: IMPLICATIONS FOR STREAM WATER QUALITY1

ABSTRACT: Riparian buffers have potential for reducing excess nutrient levels in surface water. Spatial variation in riparian buffer effectiveness is well recognized, yet researchers and managers still lack effective general tools for understanding the relevance of different hydrologic settings. We present several terrain‐based GIS models to predict spatial patterns of shallow, subsurface hydrologic flux and riparian hydrology. We then link predictions of riparian hydrology to patterns of nutrient export in order to demonstrate potential for augmenting the predictive power of land use/land cover (LU/LC) maps. Using predicted hydrology in addition to LUILC, we observed increases in the explained variation of nutrient exports from 290 sites across Lower Michigan. The results suggest that our hydrologic predictions relate more strongly to patterns of nutrient export than the presence or absence of wetland vegetation, and that in fact the influence of vegetative structure largely depends on its hydrologic context. Such GIS models are useful and complimentary tools for exploring the role of hydrologic routing in riparian ecosystem function and stream water quality. Modeling efforts that take a similar GIS approach to material transport might be used to further explore the causal implications of riparian buffers in heterogeneous watersheds.     

APPLICATION OF THE SNOWMELT RUNOFF MODEL USING MULTIPLE-PARAMETER LANDSCAPE ZONES ON THE TOWANDA CREEK BASIN,PENNSYLVANIA1

ABSTRACT: The Snowmelt Runoff Model (SRM) is designed to compute daily stream discharge using satellite snow cover data for a basin divided into elevation zones. For the Towanda Creek basin, a Pennsylvania watershed with relatively little relief, analysis of snow cover images revealed that both elevation and land use affected snow accumulation and melt on the landscape. The distribution of slope and aspect on the watershed was also considered; however, these landscape features were not well correlated with the available snow cover data. SRM streamflow predictions for 1990, 1993 and 1994 snowmelt seasons for the Towanda Creek basin using a combination of elevation and land use zones yielded more precise streamflow estimates than the use of standard elevation zones alone. The use of multiple-parameter zones worked best in non-rain-on-snow conditions such as in 1990 and 1994 seasons where melt was primarily driven by differences in solar radiation. For seasons with major rain-on-snow events such as 1993, only modest improvements were shown since melt was dominated by rainfall energy inputs, condensation and sensible heat convection. Availability of GIS coverages containing satellite snow cover data and other landscape attributes should permit similar reformulation of multiple-parameter watershed zones and improved SRM streamflow predictions on other basins.     

SPATIAL MODELLING WITH GEOGRAPHIC INFORMATION SYSTEMS FOR DETERMINATION OF WATER RESOURCES VULNERABILITY APPLICATION TO AN AREA IN MASSIF CENTRAL (FRANCE)1

ABSTRACT: A method for water resources protection based on spatial variability of vulnerability is proposed. The vulnerability of a water resource is defined as the risk that the resource will become contaminated if a pollutant is placed on the surface at one point as compared to another. A spatial modelling method is defined in this paper to estimate a travel time between any point of a catchment and a resource (river or well). This method is based on spatial analysis tools integrated in Geographical Information Systems (GIS). The method is illustrated by an application to an area of Massif Central (France) where three different types of flow appear: surface flow, shallow subsurface flow, and permanent ground water flow (baseflow). The proposed method gives results similar to classical methods of estimation of travel time. The contribution of GIS is to improve the mapping of vulnerability by taking the spatial variability of physical phenomena into account.     

TEMPORAL AND GEOMORPHIC VARIATIONS OF STREAM STABILITY AND MORPHOLOGY: MAHOGANY CREEK,NEVADA1

ABSTRACT: Detailed studies of long-term management impacts on rangeland streams are few because of the cost of obtaining detailed data replicated in time. This study uses government agency aquatic habitat, stream morphologic, and ocular stability data to assess land management impacts over four years on three stream reaches of an important rangeland watershed in northwestern Nevada. Aquatic habitat improved as riparian vegetation reestablished itself with decreased and better controlled livestock grazing. However, sediment from livestock disturbances and road crossings and very low stream flows limited the rate of change. Stream type limited the change of pool variables and width/depth ratio, which are linked to gradient and entrenchment. Coarse woody debris removal due to previous management limited pool recovery. Various critical-element ocular stability estimates represented changes with time and differences among reaches very well. Ocular stability variables tracked the quantitative habitat and morphologic variables well enough to recommend that ocular surveys be used to monitor changes with time between more intensive aquatic surveys.     

VARIABILITY IN MEASURES OF MACROINVERTEBRATE COMMUNITY STRUCTURE BY STREAM REACH AND STREAM CLASS1

ABSTRACT: Macroinvertebrate community data collected from streams in Wyoming were assessed at various scales: within one stream reach, between stream reaches within one stream, between streams, and between stream classes. Fourteen indices including number of individuals/m², biomass/m², number of taxa, Shannon's diversity index, and functional feeding group ratios were used to compare macroinvertebrates by stream reach and stream class. Statistical analysis indicated that for five of the 14 indices, significant variability occurred between macroinvertebrate communities within one reach. For two of the remaining nine indices there was significant variability between communities from several reaches within the same stream. For seven of the nine indices, there was significant variability among macroinvertebrate communities from streams of the same class. Variability among the macroinvertebrate communities from the three stream classes was significantly different for seven of the nine indices. ANOVA results suggest that macroinvertebrate communities from different samples within one reach and between reaches within one stream were more comparable than those from different streams and different stream classes.     

CALIBRATION OF A WATER-BALANCE MODEL FOR SMALL WATERSHEDS IN EASTERN OREGON1

ABSTRACT: The BURP water-balance model was calibrated for 13 small (0.46 to 7.00 mi²), forested watersheds in the Blue Mountains of eastern Oregon where snowmelt is the dominant source of runoff. BURP is the model name and is not an acronym. Six of the 16 parameters in BURP were calibrated. The subsurface recession coefficient and three subsurface water-storage parameters were most sensitive for simulating monthly flow. Calibrated subsurface recession coefficients ranged from 0.988 to 0.998. The subsurface-water storage parameters were calibrated at between 20 to 120 percent of their initial values obtained from a category III soil survey. That reconnaissance-level survey was apparently too broad to accurately reflect subsurface-water storage in small watersheds. Tests of model performance showed BURP is capable of producing accurate simulations of monthly flow for mountainous, snow-dominated watersheds with shallow (< 4 ft) soils when calibrated with 2 to 4 years of streamflow data. A regression of observed versus simulated monthly flows with data from all watersheds combined showed that BURP accounted for 85 percent of the variability in observed flows, which ranged from 0.01 to 20.8 inches, but underpredicted high flow months, with a slope of 1.15 that is significantly different from 1.0 (p = 0.05). Without prior calibration, subsurface-water storage parameters appeared to be the greatest source of potential error.     

PREDICTING TEMPORAL AND SPATIAL FLOOD DYNAMICS USING A PRECALIBRATED MODEL1

ABSTRACT: A large storm in December 1990 allowed the evaluation of flood predictions from a hydrologic model (TOPMODEL) that had been previously calibrated on the West Fork of Walker Branch Watershed, a gauged 37.5 ha catchment near Oak Ridge, Tennessee. The model predicts both hydrograph dynamics and the spatial distribution of overland flow using an index based on topography. Maximum extent of overland flow during the storm was determined from patterns of leaf litter removal from valley bottoms. Both the flood hydrograph and the extent of overland flow were accurately predicted using model parameters obtained from a three-month period of normal flow conditions during 1983.     

A CRITICAL APPROACH TO THE CALIBRATION OF A WATERSHED MODEL1

ABSTRACT: A complex watershed-scale water quality simulation model, the Hydrological Simulation Program-FORTRAN (HSPF) model, was calibrated for a 16 km² catchment. The simulation step size was 0.33 hours with predicted and recorded hydrologic flows compared on an annual and monthly basis during a total calibration period of four years. Unguided numerical optimization when applied alone did not yield a model parameter set with acceptable predictive capability; instead, it was necessary to apply a critical process that included sensitivity analysis, numerical optimization, and testing of derived model parameter sets to evaluate their performance for periods other than those for which they were determined. Using this critical calibration process, the model was proven to have significant predictive capability. Numerical optimization is an aid for model calibration, but it must not be used blindly.     

APPLICATION OF THE BASINS DATABASE AND NPSM MODEL ON A SMALL OHIO WATERSHED1

ABSTRACT: The purpose of this study was to evaluate the Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) watershed management system. BASINS data were used with the NPSM model to predict discharge and sediment concentrations at the outlet of a 103 km² Ohio watershed. It was concluded that the NPSM model should always be calibrated but only a few of the parameters provided with BASINS needed to be calibrated. For a three‐year study period, there was a 2 percent underestimation of discharge using area weighted precipitation values and a 25 percent overestimation using the single station data in BASINS. A comparison of observed and predicted monthly discharge resulted in an r² of 0.86 with area‐weighted precipitation and an r² of 0.74 with the single station data. Calibrating the model to substantially improve sediment predictions was unsuccessful and we concluded that a calibration period of one year was too short. For the three‐year study period, the r² for sediment was 0.36 with a slope of 0.37 and an intercept of 18.8 mg/l. The mean observed and predicted sediment concentrations were 27.1 mg/l and 22.6 mg/l, respectively.     

QUANTIFYING MODEL OUTPUT UNCERTAINTY DUE TO SPATIAL VARIABILITY OF RAINFALL1

Traditionally in the application of hydrologic/water quality (H/WQ) models, rainfall is assumed to be spatially homogeneous and is considered not to contribute to output uncertainty. The objective of this study was to assess the uncertainty induced in model outputs solely due to rainfall spatial variability. The study was conducted using the AGNPS model and the rainfall pattern captured by a network of 17 rain gauges. For each rainfall event, the model was run using the rainfall captured by each rain gauge, one at a time, under the assumption of rainfall spatial homogeneity. A large uncertainty in the modeled outputs resulted from the rainfall spatial variability. The uncertainty in the modeled outputs exceeded the input rainfall uncertainty. Results of this study indicate that spatial variability of rainfall should be captured and used in H/WQ models in order to accurately assess the release and transport of pollutants. A large uncertainty in the model outputs can be expected if this rainfall property is not taken into account.     

A CASE STUDY OF SHALLOW FLOW PATHS IN A STEEP ZERO-ORDER BASIN1

ABSTRACT: Soil water potentials, slope throughflow, runoff chemistry, and isotopic composition were monitored in a 97 m² zero-order basin within the Maimai 8 watershed on the South Island of New Zealand, for a natural rain storm and two artificial water applications. Contrary to results previously reported for other portions of the Maimai catchment, much of the runoff occurred as a shallow subsurface organic layer flow. For the 47 mm natural rain event, pre-storm soil matric potential ranged from ?60 to ?150 cm H₂O. No saturation was produced within the profile, and the majority of storm runoff emanated from flow within the organic horizon perched on the mineral soil surface. Hillslope applications corroborated this interpretation by showing >90 percent new water flushing with negligible mineral soil moisture response. Although the mechanisms cited in the text are not representative of the entire catchment, the study demonstrates: (1) the value of a combined physical-chemical-isotopic approach in quantifying slope processes, and (2) the heterogeneous nature and diversity of slope runoff pathways in a relatively homogeneous catchment.