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.     

MACROSCALE HYDROLOGIC MODELING FOR REGIONAL CLIMATE ASSESSMENT STUDIES IN THE SOUTHEASTERN UMTED STATES1

ABSTRACT: A macroscale hydrologic model is developed for regional climate assessment studies under way in the southeastern United States. The hydrologic modeling strategy is developed to optimize spatial representation of basin characteristics while maximizing computational efficiency. The model employs the “grouped response unit” methodology, which follows the natural drainage pattern of the area. First order streams are delineated and their surface characteristics are tested so that areas with statistically similar characteristics can be combined into larger computational zones for modeling purposes. Hydrologic response units (HRU) are identified within the modeling units and a simple three‐layer water balance model, Soil and Water Assessment Tool (SWAT), is executed for each HRU. The runoff values are then convoluted using a triangular unit hydrograph and routed by Muskingum‐Cunge method. The methodology is shown to produce accurate results relative to other studies, when compared to observations. The model is used to evaluate the potential error in hydrologic assessments when using GCM predictions as climatic input in a rainfall‐runoff dominated environment. In such areas, the results from this study, although limited in temporal and spatial scope, appear to imply that use of GCM climate predictions in short term quantitative analyses studies in rainfall‐runoff dominated environments should proceed with caution.     

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.     

PARAMETER UNCERTAINTY IN RAINFALL-RUNOFF MODELING FOR POLDER SYSTEM DESIGN1

ABSTRACT: An approach is developed for incorporating the uncertainty of parameters for estimating runoff in the design of polder systems in ungaged watersheds. Monte Carlo Simulation is used to derive a set of realizations of streamflow hydrographs for a given design rainstorm using the U. S. Soil Conservation Service (SCS) unit hydrograph model. The inverse of the SCS curve number, which is a function of the antecedent runoff condition in the SCS model, is the random input in the Monte Carlo Simulation. Monte Carlo realizations of streamfiow hydrographs are used to simulate the performance of a polder flood protection system. From this simulation the probability of occurrence of flood levels for a particular hydraulic design may be used to evaluate its effectiveness. This approach is demonstrated for the Pluit Polder flood protection system for the City of Jakarta, Indonesia. While the results of the application indicate that uncertainty in the antecedent runoff condition is important, the effects of uncertainty in rainfall data, in additional runoff parameters, such as time to peak, in the hydraulic design, and in the rainfall-runoff model selected should also be considered. Although, the SCS model is limited to agricultural conditions, the approach presented herein may be applied to other flood control systems if appropriate storm runoff models are selected.     

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.     

SIMULATION OF SPATIAL AND TEMPORAL CHANGES IN WATER QUALITY WITHIN A HYDROLOGIC UNIT1

Water quality must be considered in the development and planning aspects of water resource management. To accomplish this, the decision-maker needs to have at his disposal a systematized procedure for simulating water quality changes in both time and space. The simulation model should be capable of representing changes in several parameters of water quality as they are influenced by natural and human factors impinging on the hydrologic system. The objective of this work is two-fold. The first goal is to demonstrate the feasibility of developing and utilizing a water quality simulation model in conjunction with a hydrologic simulation model. The model represents water quality changes in both time and space in response to changing atmospheric and hydrologic conditions and time-varying waste discharges at various points in the system. This model has been developed from and verified with actual field data from a prototype system selected for this purpose. The second aim is to set forth procedural guidelines to assist in the development of water quality simulation models as tools for use in the quality-quantity management of a hydrologic unit.     

Hydrograph Separation by Incorporating Climatological Factors: Application to Small Experimental Watersheds 1

Abstract: Evaluating the relative amounts of water moving through the different components of the hydrological cycle is required for precise management and planning of water resources. An important aspect of this evaluation is the partitioning of streamflow into surface (quick flow) and base‐flow components. A prior study evaluated 40 different approaches for hydrograph‐partitioning on a field scale watershed in the Coastal Plain of the Southeastern United States and concluded that the Boughton’s method produced the most consistent and accurate results. However, its accuracy depends upon the proper estimation of: (1) the end of surface runoff, and (2) the fraction factor (α) that is function of many physical and hydrologic characteristics of a watershed. Proper identification of the end of surface runoff was accomplished by using a second derivative approach. Applying this approach to 12 years of separately measured surface and subsurface flow data from a field scale watershed (study area) proved to be accurate for 87% of the time. Estimation of the α value was accomplished in this study using two steps: (1) alpha was fitted to individual hydrographs: and, (2) a regression equation that determines these alpha values based on climatological factors (e.g., rainfall, evapotranspiration) was developed. Using these strategies improved the streamflow partitioning method’s performance significantly.     

CALIBRATION OF SNYDER COEFFICIENTS FOR PENNSYLVANIA1

ABSTRACT Unit hydrograph theory is one of the most widely used techniques to predict surface runoff. The present study is concerned with the Snyder unit hydrograph and the calibration of the Snyder coefficients for Pennsylvania. Twenty-seven study basins were selected, located randomly across the state. With the rainfall and runoff recorded for several events for each basin (more than 500 events were analyzed) unit hydrographs were calculated and the Snyder coefficients determined. A map of the coefficients was drawn to illustrate the variability in the coefficients and two equations using multiple regression theory were developed. The unexplained variability of the coefficients suggests that upper and lower bounds on the peak flow might be placed on storm hydrographs developed for ungaged watersheds.     

APPLICATION OF INFILTRATION THEORY FOR THE DETERMINATION OF EXCESS RAINFALL HYETOGRAPH1

ABSTRACT: The purpose of this paper is to show through the use of numerical examples that modern infiltration theory can be used in everyday hydrologic practice. The actual use of four methods of calculation of infiltration rates and of excess rainfall rates is demonstrated for the case when simultaneous data of rainfall and stream flow are available for a watershed. The four methods are: (1) the well known Π-index method, (2) the traditional Horton's infiltration capacity formula, (3) the less traditional Green and Ampt infiltration capacity formula, and (4) a ponding time approach. It is recommended that hydrologists become at least familiar with the numerical procedures involved in the ponding time and postponding infiltration approach. This approach, though not flawless, should be preferred to the other three methods if use of the other three is at all considered.     

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.     

A STRIP MODEL APPROACH TO PARAMETERIZE A COUPLED GREEN‐AMPT KINEMATIC WAVE MODEL1

ABSTRACT: Infiltration processes at the plot scale are often described and modeled using a single effective hydraulic conductivity (Kg) value. This can lead to errors in runoff and erosion prediction. An integrated field measurement and modeling study was conducted to evaluate: (1) the relationship among rainfall intensity, spatially variable soil and vegetation characteristics, and infiltration processes; and (2) how this relationship could be modeled using Green and Ampt and a spatially distributed hydrologic model. Experiments were conducted using a newly developed variable intensity rainfall simulator on 2 m by 6 m plots in a rangeland watershed in southeastern Arizona. Rainfall application rates varied between 50 and 200 mm/hr. Results of the rainfall simulator experiments showed that the observed hydrologic response changed with changes in rainfall intensity and that the response varied with antecedent moisture condition. A distributed process based hydrologic simulation model was used to model the plots at different levels of hydrologic complexity. The measurement and simulation model results show that the rainfall runoff relationship cannot be accurately described or modeled using a single Kg value at the plot scale. Multi‐plane model configurations with infiltration parameters based on soil and plot characteristics resulted in a significant improvement over single‐plane configurations.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

Estimating Lag to Peak between Rainfall and Peak Streamflow with a Mixed‐Effects Model

We test the use of a mixed‐effects model for estimating lag to peak for small basins in Maine (drainage areas from 0.8 to 78 km²). Lag to peak is defined as the time between the center of volume of the excess rainfall during a storm event and the resulting peak streamflow. A mixed‐effects model allows for multiple observations at sites without violating model assumptions inherent in traditional ordinary least squares models, which assume each observation is independent. The mixed model includes basin drainage area and maximum 15‐min rainfall depth for individual storms as explanatory features. Based on a remove‐one‐site cross‐validation analysis, the prediction errors of this model ranged from ?42% to +73%. The mixed model substantially outperformed three published models for lag to peak and one published model for centroid lag for estimating lag to peak for small basins in Maine. Lag to peak estimates are a key input to rainfall–runoff models used to design hydraulic infrastructure. The improved accuracy and consistency with model assumptions indicates that mixed models may provide increased data utilization that could enhance models and estimates of lag to peak in other regions.     

Sensitivity of SCS Models to Curve Number Variation1

ABSTRACT: The Soil Conservation Service (SCS) models, including the TR-20 computer program and the simplified methods in TR-55, are widely used in hydrologic design. The runoff curve number (CN), which is an important input parameter to SCS models, is defined in terms of land use tretments, hydrologic, condition, antecedent soil moisture, and soil type. The objective of this study was to evaluate the sensitivity of the SCS models to errors in CN estimates. The results show that the effects of CN variation decrease as the design rainfall depth increases, such as for the larger storm events. The value and use of the sensitivity curves are demonstrated using a comparison of Landsat and conventionally derived curve numbers for three watersheds in Pennsylvania.     

BASE FLOW RECESSION RATES,LOW FLOWS,AND HYDROLOGIC FEATURES OF SMALL WATERSHEDS IN PENNSYLVANIA,USA1

This paper examines the relationships between measurable watershed hydrologic features, base flow recession rates, and the Q_7,10 low flow statistic (the annual minimum seven‐day average streamflow occurring once every 10 years on average). Base flow recession constants were determined by analyzing hydrograph recession data from 24 small (>130 km²), unregulated watersheds across five major physiographic provinces of Pennsylvania, providing a highly variable dataset. Geomorphic, hydrogeologic, and land use parameters were determined for each watershed. The base flow recession constant was found to be most strongly correlated to drainage density, geologic index, and ruggedness number (watershed slope); however, these three parameters are intercorrelated. Multiple regression models were developed for predicting the recession rate, and it was found that only two parameters, drainage density and hydrologic soil group, were required to obtain good estimates of the recession constant. Equations were also developed to relate the recession rates to Q_7,10 per unit area, and to the Q_7,10/Q₅₀ ratio. Using these equations, estimates of base flow recession rates, Q_7,10, and streamflow reduction under drought conditions can be made for small, ungaged basins across a wide range of physiography.     

GENERATING DESIGN HYDROGRAPHS BY DEM ASSISTED GEOMORPHIC RUNOFF SIMULATION: A CASE STUDY1

ABSTRACT: The geomorphic instantaneous unit hydrograph (GIUH) may be one of the most successful methodologies for predicting flow characteristics in ungauged watersheds. However, one difficulty in applying the GIUH model is determination of travel time, and the other difficulty is the large amount of geomorphologic information required in the study watershed. Recently, using the kinematic-wave theory Lee and Yen (1997) have analytically determined the travel times for overland and channel flows in watersheds. The limitation of using an empirical velocity equation to estimate the runoff travel time for a specified watershed is then relaxed. To simplify the time-consuming work involved in geomorphic parameter measurement on topographic maps, the GIUH model is linked with geographic information systems to obtain geomorphic parameters from digital elevation models. In this paper, a case study performed for peak flow analysis in an ungauged watershed is presented. The geomorphic characteristics of the study watershed were analyzed using a digital elevation model and were used to construct the runoff simulation model. The design storm was then applied to the geomorphic runoff simulation model to obtain the design hydrograph. The analytical procedures proposed in this study can provide a convenient way for hydrologists to estimate hydrograph characteristics based on limited hydrologic information.     

MANAGEMENT OF FLOOD CONTROL SUMPS AND POLLUTANT TRANSPORT1

ABSTRACT: Levee sump systems are used by many riverine communities for temporary storage of urban wet weather flows. The hydrologic performance and transport of stormwater pollutants in sump systems, however, have not been systematically studied. The objective of this paper is to present a case study to demonstrate development and application of a procedure for assessing the hydraulic performance of flood control sumps in an urban watershed. Two sumps of highly variable physical and hydraulic characteristics were selected for analysis. A hydrologic modeling package was used to estimate the flow hydrograph for each outfall as part of the flow balance for the sump. To validate these results, a water balance was used to estimate the total runoff using sump operational data. The hydrologic model calculations provide a satisfactory estimate of the total runoff and its time‐distribution to the sump. The model was then used to estimate pollutant loads to the sump and to the river. Although flow of stormwater through a sump system is regulated solely by flood‐control requirements, these sumps may function as sedimentation basins that provide purification of stormwater. A sample calculation of removals of several conventional pollutants in the target sumps using a mass balance approach is presented.     

A COMPUTERIZED DATA BASE FOR FLOOD PREDICTION MODELING1

ABSTRACT: A computerized geographic information system (GIS) was created in support of data requirements by a hydrologic model designed to predict the runoff hydrograph from ungaged basins. Some geomorphologic characteristics (i.e., channel lengths) were manually measured from topographic maps, while other parameters such as drainage area and number of channels of a specified order, land use, and soil type were digitized and manipulated through use of the GIS. The model required the generation of an integrated Soil Conservation Service (SCS) curve number for the entire basin. To this end, soil associations and land use (generated from analysis of Landsat satellite data) were merged in the GIS to acquire a map representing SCS runoff curve numbers. The volume of runoff obtained from the Watershed Hydrology Simulation (WAHS) Model using this map was compared to the volume computed by hydrograph separation and found to be accurate within 19 percent error. To quantify the effect of changing land use on basin hydrology, the GIS was used to vary percentages from the drainage area from forest to bare soil. By changing the basin runoff curve numbers, significant changes in peak discharge were noted; however, the time to peak discharge remained essentially independent of change in area of land use. The GIS capability eliminated many of the more traditional manual phases of data input arid manipulation, thereby allowing researchers to concentrate on the development and calibration of the model and the interpretation of presumably more accurate results.     

RUNOFF MODELING OF A MOUNTAINOUS CATCHMENT USING TOPMODEL: A CASE STUDY1

ABSTRACT: The rainfall‐runoff response of the Tygarts Creek Catchment in eastern Kentucky is studied using TOPMODEL, a hydrologic model that simulates runoff at the catchment outlet based on the concepts of saturation excess overland flow and subsurface flow. Unlike the traditional application of this model to continuous rainfall‐runoff data, the use of TOPMOEL in single event runoff modeling, specifically floods, is explored here. TOPMODEL utilizes a topographic index as an indicator of the likely spatial distribution of rainfall excess generation in the catchment. The topographic index values within the catchment are determined using the digital terrain analysis procedures in conjunction with digital elevation model (DEM) data. Select parameters in TOPMODEL are calibrated using an iterative procedure to obtain the best‐fit runoff hydrograph. The calibrated parameters are the surface transmissivity, T_O, the transmissivity decay parameter, m, and the initial moisture deficit in the root zone, S_r0. These parameters are calibrated using three storm events and verified using three additional storm events. Overall, the calibration results obtained in this study are in general agreement with the results documented from previous studies using TOPMODEL. However, the parameter values did not perform well during the verification phase of this study.