CALIBRATING AND VERIFYING THE SSARR MODEL -MISSOURI RIVER WATERSHEDS STUDY1

The Streamflow Synthesis and Reservoir Regulation (SSARR) model was calibrated and verified on the Madison and Gallatin watersheds in the upper Missouri River drainage. The study was performed to determine if the SSARR model could simulate snowmelt-runoff volumes to effect better operation of six multipurpose reservoirs on the Missouri River. Physical watershed characteristics and parameter sensitivity are incorporated into a procedure which expedites model calibration. Criteria are established to facilitate parameter development and to objectively evaluate calibration and verification results. A ratio of simulated to observed snowmelt-runoff volumes of the Madison River averaged 1.00 and 1.02 for calibration (N = 8 years) and verification (N = 6 years) with corresponding standard deviations of 0.08 and 0.13. Gallatin volume ratios averaged 0.99 and 0.95 for calibration (N = 7 years) and verification (N = 5 years) with respective standard deviations of 0.08 and 0.28.     

THE APPLICATION OF HYDROLOGIC MODELS TO SMALL WATERSHEDS HAVING MILD TOPOGRAPHY1

ABSTRACT: The application of hydrologic models to small watersheds of mild topography is not well documented. This study evaluates the applicability of hydrologic models described by Huggins and the Soil Conservation Service to small watersheds by comparing the simulated and actual hydrograph for both gaged and ungaged situations. The annual maximum rainfall events plus storms exceeding 2.5 inches from 25 years of rainfall and runoff data for two small watersheds were selected for the model evaluations. These storms had a variety of patterns and occurred on many different watershed conditions. Simulated and actual hydrographs were compared using a parameter which contained volume, peak, and shape factors. One-half of the selected storms were used to calibrate the models. For both models, there were no significant differences between the simulated and actual runoff volumes and peak runoff rates. Parameters obtained during the calibration process and relationships developed to estimate antecedent moisture and to modify tabulated runoff curve numbers were used to simulate the runoff hydrograph from the remaining storms. These remaining storms or test storms were simulated only once in order to imitate an ungaged situation. In general, both the Huggins and SCS model performed similarly on the test storms, but the level of model performance was lower than that for the calibration storms. For both models, the two-day antecedent rainfall was more important than the five-day in determining antecedent moisture and modifying tabulated curve numbers. The time of concentration which resulted in good hydrograph simulations was about three times larger than that estimated using published empirical relationships.     

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 STORM PATTERNS ON RUNOFF HYDROGRAPHS1

ABSTRACT: As an alternative to the conventional single-peak design storms commonly used in hydrologic practice, a large number of Southeastern Pennsylvania storm events were selected from hourly U.S. National Oceanographic and Atmospheric Administration (NOAA) records, and their temporal distributions were analyzed. From these recorded events, design storms of a typical distribution were developed for storm durations between 6 and 18 hours. All of these generated design storms have two or more peaks. The conventional single peak as well as the “typical” multi-peak storms were then applied to a simulated watershed. It was found that the multi-peak storms consistently produced more dispersed hydrographs with lower runoff peaks than the conventional single peak storms.     

HYDROLOGIC CALIBRATION OF THE SWAT MODEL IN A WATERSHED CONTAINING FRAGIPAN SOILS1

ABSTRACT: The Soil and Water Assessment Tool (SWAT) model, designed for use on rural ungaged basins and incorporating a GRASS GIS interface, was used to model the hydrologic response of the Ariel Creek watershed of northeastern Pennsylvania. Model evaluation of daily flow prior to calibration revealed a deviation of runoff volumes (D_v) of 68.3 percent and a Nash-Sutcliffe coefficient of-0.03. Model performance was affected by unusually large observed snowmelt events and the inability of the model to accurately simulate baseflow, which was influenced by the presence of fragipans. Seventy-five percent of the soils in the watershed contain fragipans. Model calibration yielded a D_v of 39.9 percent and a Nash-Sutcliffe coefficient of 0.04, when compared on a daily basis. Monthly comparisons yielded a Nash-Sutcliffe coefficient of 0.14. Snowmelt events in the springs of 1993 and 1994, which were unusually severe, were not adequately simulated. Neglecting these severe events, which produced the largest and third largest measured flows for the period of record, a D_v of 4.1 percent and Nash-Sutcliffe coefficient of 0.20 were calculated on a daily comparison, while on a monthly basis the Nash-Sutciffe coefficient was 0.55. These results suggest that the SWAT model is better suited to longer period simulations of hydrologic yields. Baseflow volumes were accurately simulated after calibration (D_v= -0.2 percent). Refinements made to the algorithms controlling subsurface hydrology and snowmelt, to better represent the presence of fragipans and snowmelt events, would likely improve model performance.     

ASPEN CLEARCUTTING INCREASES SNOWMELT AND STORM FLOW PEAKS IN NORTH CENTRAL MINNESOTA1

ABSTRACT: Clearcutting aspen from the upland portion of an upland peatland watershed in north central Minnesota caused snowmelt peak discharge to increase 11 to 143 percent. Rainfall peak discharge size increased as much as 250 percent during the first two years after clearcutting, then decreased toward precutting levels in subsequent years. Storm flow volumes from rain during the first two years increased as much as 170 percent but declined to preharvest volumes in the third year. Snowmelt volumes did not significantly change. Snowmelt peak discharge occurred about four to five days earlier after clearcutting, but the timing of storm flow from rainfall was not changed. Snowmelt peaks remained above precut size for nine years after clearcutting on an area undergoing natural regeneration to aspen saplings. Partial cutting - up to approximately one-half of the watershed - reduced peak snowmelt discharge because melt was desynchronized in cleared and forested parts. Clearing more than 2/3 of the watershed caused snowmelt flood peak size to double during years with snow packs in excess of seven inches of water that remained until a day when maximum air temperatures exceeded 60d?F.     

FLOOD PEAKS AS MODIFIED BY DAM SIZE AND LOCATION1

ABSTRACT. The ability to replicate storms on a rainfall simulator over an iconic watershed model permits the study of the effectiveness of dams of different sizes and locations in reducing flood peaks. This study empirically confirms that, when storms uniformly and completely cover the watershed, large dams are nearly twice as effective as a network of small and medium dams in reducing the flood.     

Calibration and Verification of SWMM for Low Impact Development

The Storm Water Management Model was used to simulate runoff and nutrient export from a low impact development (LID) watershed and a watershed using traditional runoff controls. Predictions were compared to observed values. Uncalibrated simulations underpredicted weekly runoff volume and average peak flow rates from the multiple subcatchment LID watershed by over 80%; the single subcatchment traditional watershed had better predictions. Saturated hydraulic conductivity, Manning's n for swales, and initial soil moisture deficit were sensitive parameters. After calibration, prediction of total weekly runoff volume for the LID and traditional watersheds improved to within 12 and 5% of observed values, respectively. For the validation period, predicted total weekly runoff volumes for the LID and traditional watersheds were within 6 and 2% of observed values, respectively. Water quality simulation was less successful, Nash–Sutcliffe coefficients >0.5 for both calibration and validation periods were only achieved for prediction of total nitrogen export from the LID watershed. Simulation of a 100‐year, 24‐h storm resulted in a runoff coefficient of 0.46 for the LID watershed and 0.59 for the traditional watershed. Results suggest either calibration is needed to improve predictions for LID watersheds or expanded look‐up tables for Green–Ampt infiltration parameter values that account for compaction of urban soil and antecedent conditions are needed.     

THE EFFECTS OF A MOVING RAINSTORM ON DIRECT RUNOFF PROPERTIES1

ABSTRACT: The effects of a moving rainstorm on flood runoff characteristics were investigated. A flood hydrograph simulation model called “FH-Model” and a natural watershed were used. A hypothetical rainstorm of 50 years recurrence interval, 75 mm depth, and 4 hours duration was used to show the effects of velocity and direction of the moving rainstorm on the runoff characteristics. Compared with an equivalent stationary rainstorm (ESRS), the peak flow caused by a rainstorm moving in a downstream direction with a speed equal to channel velocity, V, was 27.5 percent higher and the peak flow caused by the same rainstorm moving in an upstream direction was 21.7 percent smaller. These percentages reduced to 10.5 percent and 8.6 percent for storms moving downstream and upstream, respectively, at three times the channel velocity, 3V. There were negligible differences in the time of peak, T_p between runoff caused by storms moving downstream and runoff produced by ESRS. However, T_p for a storm moving upstream at V velocity was 82 percent higher than that produced by ESRS, but was reduced to 27 percent higher when the storm velocity was 3V.     

THE USE OF FLOOD POTENTIAL INDICES FOR FLOOD PEAK ESTIMATION ON UNGAGED WATERSHEDS1

: The construction of a flood peak index map was attempted for use by hydrologists in the simple format of rainfall maps. Since flood peaks are highly dependent on watershed area, the effect of area was removed. By regression analysis flood peaks of 2.33 and 100-year return periods were found to be proportional to watershed area to the 0.8 and 0.7 powers, respectively. Therefore, indices C_{2 33}= Q_{2 33}/_A0.7 were completed at each gage and plotted on a Pennsylvania map. It was attempted to further remove some of the scatter by regression of C with several other watershed parameters like slope, percent forest cover, and watershed shape, but no significant correlation could be found. The index maps, drawn without attenuation of the scatter, can be used by hydrologists to compute flood peaks as Q = CAⁿ (with n = 0.8 and 0.7 for the 2.33 and 100-year flood peaks, respectively). Flood peak safety factors can be based on visual observation of the index variation in the vicinity of the location for which the flood peak estimate is needed.     

RECHARGE ESTIMATES USING A GEOMORPHIC/DISTRIBUTED-PARAMETER SIMULATION APPROACH,AMARGOSA RWER BASIN1

ABSTRACT: Average-annual volumes of runoff, evapotranspiration, channel loss, upland (interchannel) recharge, and total recharge were estimated for watersheds of 53 channel sites in the Amargosa River basin above Shoshone, California. Estimates were based on a water-balance approach combining field techniques for determining streamflow with distributed-parameter simulation models to calculate transmission losses of ephemeral streamflow and upland recharge resulting from high-magnitude, low-frequency precipitation events. Application of the water-balance models to the Amargosa River basin, Nevada and California, including part of the Nevada Test Site, suggests that about 20.5 million cubic meters of water recharges the ground-water reservoir above Shoshone annually. About 1.6 percent of precipitation becomes recharge basinwide. About 90 percent of the recharge is by transmission loss in channels, and the remainder occurs when infrequent storms yield sufficient precipitation that soil water percolates beyond the rooting zone and reaches the zone of saturation from interchannel areas. Highest rates of recharge are in headwaters of the Amargosa River and Fortymile Wash; the least recharge occurs in areas of relatively low precipitation in the lowermost Amargosa River watershed.     

HYDROLOGIC MODELING OF TWO GLACIATED WATERSHEDS IN NORTHEAST PENNSYLVANIA1

ABSTRACT: A hydrologic modeling study, using the Hydrologic Simulation Program - FORTRAN (HSPF), was conducted in two glaciated watersheds, Purdy Creek and Ariel Creek in northeastern Pennsylvania. Both watersheds have wetlands and poorly drained soils due to low hydraulic conductivity and presence of fragipans. The HSPF model was calibrated in the Purdy Creek watershed and verified in the Ariel Creek watershed for June 1992 to December 1993 period. In Purdy Creek, the total volume of observed stream-flow during the entire simulation period was 13.36 × 10⁶ m³ and the simulated streamflow volume was 13.82 × 10⁶ m³ (5 percent difference). For the verification simulation in Ariel Creek, the difference between the total observed and simulated flow volumes was 17 percent. Simulated peak flow discharges were within two hours of the observed for 30 of 46 peak flow events (discharge greater than 0.1 m³/sec) in Purdy Creek and 27 of 53 events in Ariel Creek. For 22 of the 46 events in Purdy Creek and 24 of 53 in Ariel Creek, the differences between the observed and simulated peak discharge rates were less than 30 percent. These 22 events accounted for 63 percent of total volume of streamflow observed during the selected 46 peak flow events in Purdy Creek. In Ariel Creek, these 24 peak flow events accounted for 62 percent of the total flow observed during all peak flow events. Differences in observed and simulated peak flow rates and volumes (on a percent basis) were greater during the snowmelt runoff events and summer periods than for other times.     

CALIBRATION OF STORM LOADS IN THE SOUTH PRONG WATERSHED,FLORIDA, USING BASINS/HSPF1

ABSTRACT: The South Prong watershed is a major tributary system of the Sebastian River and adjacent Indian River Lagoon. Continued urbanization of the Sebastian River drainage basin and other watersheds of the Indian River Lagoon is expected to increase runoff and nonpoint source pollutant loads. The St. Johns River Water Management District developed watershed simulation models to estimate potential impacts on the ecological systems of receiving waters and to assist planners in devising strategies to prevent further degradation of water resources. In the South Prong system, a storm water sampling program was carried out to calibrate the water quality components of the watershed model for total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN). During the period of May to November 1999, water quality and flow data were collected at three locations within the watershed. Two of the sampling stations were located at the downstream end of major watercourses. The third station was located at the watershed outlet. Five storm events were sampled and measured at each station. Sampling was conducted at appropriate intervals to represent the rising limb, peak, and recession limb of each storm event. The simulations were handled by HSPF (Hydrologic Simulation Program‐Fortran). Results include calibration of the hydrology and calibration of the individual storm loads. The hydrologic calibration was continuous over the period 1994 through 1999. Simulated storm runoff, storm loads, and event mean concentrations were compared with their corresponding observed values. The hydrologic calibration showed good results. The outcome of the individual storm calibrations was mixed. Overall, however, the simulated storm loads agreed reasonably well with measured loads for a majority of the storms.     

FIRST-YEAR EFFECTS OF CLEARCUTTING AN OAK-HICKORY WATERSHED ON WATER YIELD1

ABSTRACT: Two intermittent streams on oak-hickory watersheds in southern Illinois were gaged with a V-notch weir and sampled with an automatic water sampler. Baseline data was collected for a period of three years. Flow volume showed large variations between years and watersheds. Water samples were analyzed for Na, K, Ca, Mg, ortho-P, and NO₃-N. Water quality was consistently high, but there were significant differences between the watersheds during the calibration period. One watershed was clearcut in November 1979. One year of postharvest data has been analyzed. Flow volume increased 95 percent, but there was no evidence of increased sedimentation. There were significant increases in the stream water concentrations of K, Mg, and NO₃-N of 18 percent, 8 percent, and 274 percent, respectively. Nutrient budgets for the site were not adversely affected by the harvest. The clearcutting operation appears to have had a small impact on the watershed due to minimal disturbance during the logging and below normal precipitation the first year following the harvest.     

NONPARAMETRIC STATISTICAL METHODS IN URBAN HYDROLOGIC RESEARCH1

ABSTRACT. In urban hydrologic studies, it is often necessary to determine the effect of changes in urban land use patterns on such runoff characteristics as flood peaks and flow volumes. Nonparametric statistical methods have certain properties that make them a valuable tool for detecting hydrologic change caused by a treatment, such as urbanization, that changes watershed over a period of time. As many hydrologists do not have a working familiarity with nonparametric methods, a number of them are used for illustrative purposes to analyze the effect of urbanization on 24 years of annual flood peaks for a Louisville, Kentucky, watershed. In the example, urbanization was found to increase the central tendency, but not the dispersion of the peaks. Peak flows modeled by holding watershed parameters constant were also found to be increasing because of an upward trend in precipitation. By following the numerical examples in the paper and looking up test statistics in referenced sources, the reader can easily apply these methods to other situations.     

DEVELOPMENT AND EVALUATION OF A DIMENSIONLESS UNIT HYDROGRAPH1

ABSTRACT: A generalized unit hydrograph method is developed and evaluated for ungaged watersheds. A key component in this method is the value of a dimensionless storage coefficient. Procedures to estimate this coefficient are given using calibrated values from 142 rainfall-runoff events gaged in watershed located mainly in the Eastern US. Only limited success was obtained in predicting this storage coefficient. Thirty-seven, independent rainfall-runoff events were used to test the proposed technique. The generalized unit hydrograph predicted the observed runoff hydrographs fairly well with considerable improvement in accuracy over the SCS dimensionless unit hydrograph. Approximately one-half of test storms had percent errors in predicted peak flow rates that were less than 34 percent compared to percent error of 88 percent with the SCS method.     

BMP IMPACTS ON WATERSHED RUNOFF,SEDIMENT, AND NUTRIENT YIELDS1

ABSTRACT: To quantify the effectiveness of best management practice (BMP) implementation on runoff, sediment, and nutrient yields from a watershed, the Nomini Creek watershed and water quality monitoring project was initiated in 1985, in Westmoreland County, Virginia. The changes in nonpoint source (NPS) loadings resulting from BMPs were evaluated by comparing selected parameters from data series obtained before, during, and after periods of BMP implementation. The results indicated that the watershed-averaged curve number, sediment, and nutrient (N and P) concentrations were reduced by approximately 5, 20, and 40 percent, respectively, due to BMP implementation. The nutrient yield model developed by Frere et al. (1980) was applied to the water quality parameters from 175 storms, but it failed to adequately describe the observed phenomena. Seasonal changes in nutrient availability factors were not consistent with field conditions, nor were they significantly different in the pm- and post-BMP periods. An extended period of monitoring, with intensive BMP implementation over a larger portion of the watershed, is required to identify BMP effectiveness.     

APPLICATION OF THE GREEN-AMPT INFILTRATION EQUATION TO WATERSHED MODELING1

ABSTRACT: A computer model was developed, based on the Green-Ampt infiltration equation, to computed rainfall excess for a single precipitation event. The model requires an estimate of parameters related to hydraulic conductivity, wetting front section, and fillable porosity of the soil layers. Values of parameters were estimated from soil textural averages or regression equations based on percent sand, percent clay, and porosity. Average values of effective porosity and wetting front suction were largely acceptable due to the relatively low variability and low model sensitivity to the parameters. Hydraulic conductivity was the most erratic constituent of the loss rate computation due to the high variability and the high sensitivity of the computed infiltration to the parameter. The performance of the Green-Ampt infiltration model was tested through a comparison with the SCS curve number procedure. Seven watersheds and 23 storms with precipitation of one inch or greater were used in the comparison. For storms with less than one inch of rainfall excess, the SCS curve number procedure generally gave the best results; however, for six of the seven storms with precipitation excess greater than one inch, the Green-Ampt procedure delivered better results. In this comparison, both procedures used the same initial abstractions. The separation of rainfall losses into infiltration, interception, and surface retention is, in theory, an accurate method of estimating precipitation excess. In the second phase of the study using nine watersheds and 39 storms, interception and surface retention losses were computed by the Horton equations. Green-Ampt and interception parameters were estimated from value sin the literature, while the surface retention parameter was calibrated so that the computed runoff volumes matched observed volumes. A relationship was found between the surface retention storage capacity and the 15-day antecedent precipitation index, month of year, and precipitation amount.     

NORMALIZED UNIT HYDROGRAPH AGGREGATION1

ABSTRACT: The unit hydrograph is a common tool in hydraulic design. Used correctly, it allows a design engineer to estimate a runoff hydrograph from a drainage basin given a rainfall event. The typical method for estimating a unit hydrograph for a gaged watershed is by deconvolution. However, distinct storms produce different unit hydrographs for a single watershed. Consequently, a design engineer usually develops a composite, or average, unit hydrograph based on several recorded storm events. Common methods for estimating this composite unit hydrograph include curve fitting, simple aggregation, and multistorm optimization techniques. This paper introduces a new method to perform aggregation of unit hydrographs. The method is an extension to the simple averaging technique, in which prior to averaging, the individual unit hydrograph time ordinates are normalized with respect to the average time to peak. The normalization method is compared to a simple averaging technique and two multistorm aggregation techniques at six rural watersheds in Alabama. The results indicate that on average the normalization method predicts runoff nearly as accurately as the multistorm techniques, and displays improvement for 60 percent of the storms tested when compared with the simple averaging technique.     

PATTERNS OF SUSPENDED SEDIMENT TRANSPORT IN A COASTAL ALASKA STREAM1

Suspended sediment data from a 154 ha watershed on northeast Chichagof Island, Alaska, were collected over three fall storm seasons from 1980 to 1982. Sediment rating curves for nine pooled storms explained less than 34 percent of the variation in total suspended solids (TSS). Significantly higher concentrations of suspended sediment occurred during the rising limb of storm hydrographs than for similar flows on the falling limb, accounting for hysteresis loops in TSS versus streamflow plots for individual storms. These hysteresis loops were wider during early season storms, indicating that easily transportable fine sediment may have been flushed from the upper portion of channel banks and from behind large organic debris during early season peak flows. Regression relationships (TSS versus Q) developed for the highest stormflows (> 1 m³/s) had steeper slopes than the lower stormflows (< 1 m³/s). Turbidity correlated well (r=0.94) with TSS for all storm-flow data combined. Organic matter constituted an average of 35 percent (by weight) of TSS for all water quality samples.