SIMULATION OF RUNOFF FROM URBAN WATERSHEDS1

ABSTRACT. .A mathematical model for urban watersheds is being developed in stages at the Utah Water Research Laboratory, Utah State University at Logan. In verifying the watershed as a unit, watershed coefficients are determined on the computer, and related to the urbanization characteristics. The second stage of verification consists of dividing the watershed into subzones, and determining the urban parameters within each subzone. Each subzone is then individually modeled, and outflow hydrographs are routed through succeeding downstream subzones to the gaging point. The model thus makes it possible to: (a) develop runoff models for subzone hydrographs within the urban watershed, and (b) account for spatial variations of storm and watershed characteristics. An attempt was also made to analytically model the outflow hydrograph based on storm and watershed characteristics.     

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.     

ADAPTATION OF FLOOD PEAKS AND DESIGN HYDROGRAPHS FROM GAGED TO NEARBY UNGAGED WATERSHEDS1

ABSTRACT: Equations were developed to transform peak flows and to adapt design hydrographs and unit hydrographs from gaged watersheds to ungaged watersheds with similar hydrologic characteristics. Dimensional analysis was used to develop adjustment equations for peak flow and time base, and these equations were reinforced with results from regional flood frequency research. The authors believe that the use of these transformation equations should yield more reliable flood peak values and hydrogrphs than the common use of empirical flood estimating curves or equations.     

OBSERVATIONS ON KINEMATIC RESPONSE IN URBAN RUNOFF MODELS1

ABSTRACT: This paper first discusses the results of sensitivity analyses conducted on various parameters of the San Francisco Stormwater Model ta version of WREM) and the Penn State Runoff Model in terms of their impact on outflow hydrographs. The parameters considered within a idealized catchment include: basin shape, imperivous fraction, overland roughness and slope: deterntion depth; infiltration capacity; and hyetograph timing. Second, the results for the hypothetical catchment are extended to the lazzard laboratory surfaces (asphalt, grass, roofing material) as a mean of illustrating the need for changes in model structure, as opposed to continued parameter adjustment Finally the effect of altering the scale of hydraulic representation in the surface runoff and sewer transport calculations are demonstrated for two gaged watersheds in Hamburg, West Germany.     

Modeling Hydrology in a Small Rocky Mountain Watershed Serving Large Urban Populations1

Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km²) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.     

MODELING RAINFALL-RUNOFF PROCESS INCLUDING LOSSES1

ABSTRACT: Existing discrete, linear rainfall-runoff models generally require the effective rainfall of a given storm as the input for computing the runoff hydrograph. This paper proposes a method for estimating, simultaneously, the optimal values of model parameters and the rainfall losses frem the measured rainfall hyetograph and the runoff hydrograph. The method involves an ARMA model for the rainfall-runoff process and a nonlinear iterative technique. The number of model parameters to be estimated for the ARMA model is much less than the unit hydrograph model. Applications of the model to three different watersheds show that the computed runoff hydrographs agree well with the measurements.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.     

Environmental policy indicators: A systems model

This article describes and tests a systems theory-based policy indicators model. The framework is used to examine propositions about linkages between states' ecological-spatial characteristics and subsequent selected solid waste management (SWM) -related environmental policies. It was hypothesized that state characteristics of: (1) population density (used as a garbage-per-land area index), (2) population convergence within urban areas, and (3) percent population change in the interval 1980–1985, could jointly explain state variation in both the number and the vigor of SWM policy outputs. Greater levels of spatial pressure were proposed to be related directly to more numerous, more convincing policies. Proposals are grounded in the literature of organizational search theory, crisis stimulation, and technological pressure. Results revealed that the sociospatial model in fact could explain a reasonable proportion of policy variation across states. However, not all hypotheses are supported. Population change shows an indirect, rather than the anticipated direct, relationship with policy output levels. In addition, when used in the model as a pollution intensity index, population density failed to contribute significantly to an explanation of differences in state SWM policy levels. The analysis raises questions about changes occurring over time in the nature and direction of linkages between sociospatial measures and policy responses. This study suggests that strengthening policy indicator models may require questioning key assumptions and theoretical bases, conducting longitudinal studies, and factoring in political, economic, and other policy environment forces.     

DUAL URBAN AND RURAL HYDROGRAPH SIGNALS IN THREE SMALL WATERSHEDS1

ABSTRACT: Many studies can be found in the literature pertaining to the effects of urbanization on surface runoff in small watersheds and the hydrologic response of undeveloped watersheds. However, an extensive literature review yielded few published studies that illustrate differing hydrologic responses from multiple source areas within a watershed. The concepts discussed here are not new, but the methods used provide a unique, basic procedure for investigating stormwater hydrology in topographically diverse basins. Six storm hydrographs from three small central Pennsylvania watersheds were analyzed for this paper; five are presented. Two important conclusions are deduced from this investigation. First, in all cases we found two distinct peaks in stream discharge, each representing different contributing areas to direct discharge with greatly differing curve numbers and lags representative of urban and rural source regions. Second, the direct discharge represents only a small fraction of the total drainage area with the urban peak becoming increasingly important with respect to the rural peak with the amount of urbanization and as the magnitude of the rain event decreases.     

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.     

THE UNIT HYDROGRAPH: A SATISFACTORY MODEL OF WATERSHED RESPONSE?1

ABSTRACT: Some 96 flood events larger than the mean annual flood from 16 watersheds in the Commonwealth of Pennsylvania were used to derive unit hydrographs by the least-squares method. Analyses of the unit hydrographs were conducted to ascertain their response to watershed parameters, climatic and storm variables and locations within different hydrologic regions. Significant differences both within and among watersheds led to the formulation and testing of hypotheses stating that differences among watersheds are caused by physiographic differences while differences within watersheds result from climatic and storm differences. The analysis showed, that while many watersheds parameters strongly influence the shape of the unit hydrograph, only the storm variables duration and volume of precipitation excess produce significant differences. Seasonal differences were apparent but not proven statistically significant.     

EVALUATION OF NINE MODELS FOR UNGAUGED URBAN BASINS IN LOUISIANA1

ABSTRACT: Nine flood-estimation models used for ungauged urban watersheds in Louisiana were evaluated. Flood-quantile predictions from simple regression models calibrated by local data were found to be more reliable than those more complicated models or models with many parameters that may not be accurately estimated. Flood prediction from models developed by using regionalization techniques were found to be reasonably good. Finally, application of a model outside of its limitations or domain may lead to substantial prediction error.     

EFFECTS OF SIMULATED CLOUD SEEDING ON STREAMFLOW OF SELECTED WATERSHEDS IN PENNSYLVANIA1

A mathematical model was developed to simulate the hydrologic behavior of five small watersheds in central Pennsylvania. Continuous hydrographs for the 6-month period, April to September 1964, were simulated. Synthesized rainfall cycles consisting of increasing rainfall by 10, 20, and 30 percent to simulate the effects of cloud seeding were processed through the watershed model to determine the effects on low flow augmentation. Other rainfall cycles used consisted of increasing every third storm by 30 percent and of developing a rainfall cycle by processing daily radiosonde data through a mathematical cumulus cloud model to obtain a prediction of rainfall following seeding. A comparison of actual and predicted hydrographs indicated that simulated cloud seeding resulted in significant monthly and seasonal water yields. In general, the results of the study appear to indicate that on a theoretical basis cloud seeding would be a feasible method of augmenting low stream-flow during the summer months on watersheds in the northern Appalachian region.     

UNIT HYDROGRAPHS DERWED FROM THE NASH MODEL1

ABSTRACT: Unit hydrograph ordinates are often estimated by deconvoluting excess rainfall pulses and corresponding direct runoff. The resulting ordinates are given at discrete times spaced evenly at intervals equal to the duration of the rainfall pulse. If the new duration is not a multiple of the parent duration, hydrograph interpolation is required. Linear interpolation, piece-wise nonlinear interpolation and graphical smoothing have been used. These interpolation schemes are expedient but they lack theoretical basis and can lead to undesirable results. Interpolation can be avoided if the instantaneous unit hydrograph (IUH) for the watershed is known. Here two issues connected with the classic Nash IUH are examined: (1) how should the Nash parameters be estimated? and (2) under what conditions is the resulting hydrograph able to reasonably represent watershed response? In the first case, nonlinear constrained optimization provides better estimates of the IUH parameters than does the method of moments. In the second case, the Nash IUH gives good results on watersheds with mild shape unit hydrographs, but performs poorly on watersheds having sharply peaked unit hydrographs. Overall, in comparison to empirical interpolation alternatives, the Nash IUH offers a theoretically sound and practical approach to estimate unit hydrographs for a wide variety of watersheds.     

ESTIMATING PEAK RUNOFF FROM FIELD-SIZE WATERSHEDS1

ABSTRACT: A simple nonlinear runoff model was developed and tested for use on field-size agricultural watersheds. A Wooding idealization of the watershed topography was used. Kinematic wave equations were used with an assumed, instead of computed, overland flow, watersurface profile in order to simplify the numerical computations. The approach was used to synthesize runoff hydrographs for an agricultural watershed in Iowa. The accuracy of the synthesized hydro-graphs was judged by comparing the estimated and observed peak discharges and by comparing estimated and observed stages at the measuring weir. The mean errors were 0.01 in/hr and 0.05 ft, respectively. A qualitative comparison was also made with a detailed kinematic wave study. The largest variability occurred during the seedbed period for both models, which was attributed to changes in surface roughness. The roughness was more constant and the results more consistent for the canopy and ground residue periods.     

Effects of Changing Forest and Impervious Land Covers on Discharge Characteristics of Watersheds

Effects of changing patterns of forest and impervious land covers on hydrologic regimes of watersheds were evaluated for urban and rural areas of the lower Cedar River drainage near Seattle, Washington. Land cover characterizations were used in a spatially explicit hydrology model to assess effects of land covers on watershed hydrology during presettlement conditions (full forest cover), 1991 and 1998. For the presettlement to 1991 period, urban watersheds showed decreases in forest covers (range 63% to 83%) and increases in impervious surfaces (range 43% to 71%). Rural watersheds showed similar patterns but smaller changes, with forest covers decreasing (range 28% to 34%) and impervious surfaces increasing (range 8% to 15%). For the 1991 and 1998 period, changes in forest covers for urban and rural watershed were <24%, with losses in some watersheds and regeneration in others. Impervious surfaces continued to increase, but increases were larger in rural (range 38% to 60%) than in urban watersheds (range 4% to 27%). Flood-frequency curves indicated that discharge rates (m sec^–1) for all watersheds were higher in 1991 and 1998 than historical and suggested that chances for floods increase because of changing land covers. The largest increases in discharge rates were in urban watersheds, with rates for 2-year, 10-year, and 25-year recurrence intervals being more than two times greater than the rate during historical conditions. Changes in flow regimes were indicated by presettlement discharge levels of less frequent recurrence intervals (10-year and 25-year) occurring in posturbanization times (1991 and 1998) during more frequent intervals (2-year and 10-year). Normalized flows (m yr^–1) of watersheds for 2-year, 10-year, and 25-year recurrence intervals indicate how flow regimes in 1991 and 1998 can change as functions of different areas of land covers. During 1991 and 1998, abrupt increases in flows occurred when forest covers were low (range 17% to 37%) and impervious surfaces were >46%. In contrast, the lowest flows occurred when forest covers were most extensive (range 59% to 81%) and impervious surfaces were <23%. We conclude that our use of spatial characterizations of impervious surfaces and forested covers in a spatially explicit hydrology model provides a robust approach for revealing how variations in different types and spatial distributions of land covers can affect flood regimes and flows of different watersheds.     

APPLICATION OF AN INTEGRATED BASIN-SCALE HYDROLOGIC MODEL TO SIMULATE SURFACE-WATER AND GROUND-WATER INTERACTIONS1

ABSTRACT: Hydrologic models have become an indispensable tool for studying processes and water management in watersheds. A physically-based, distributed-parameter model, Basin-Scale Hydro-logic Model (BSIIM), has been developed to simulate the hydrologic response of large drainage basins. The model formulation is based on equations describing water movement both on the surface and in the subsurface. The model incorporates detailed information on climate, digital elevation, and soil moisture budget, as well as surface-water and ground-water systems. This model has been applied to the Big Darby Creek Watershed, Ohio in a 28-year simulation of rainfall-runoff processes. Unknown coefficients for controlling runoff, storativity, hydraulic conductivity, and streambed permeability are determined by a trial-and-error calibration. The performance of model calibration and predictive capability of the model was evaluated based on the correlation between simulated and observed daily stream discharges. Discrepancies between observed and simulated results exist because of limited precipitation data and simplifying assumptions related to soil, land use, and geology.     

A DISTRIBUTED DYNAMIC WATERSHED MODEL1

ABSTRACT: A distributed watershed model was developed to mathematically simulate overland and channel flow for a single-event storm. The modeled watersheds in the study were subdivided into rectangular grid elements. All hydrologically significant parameters, such as land slope, rainfall and precipitation excess, were assumed to be uniform within each element. The Green-Ampt method was adopted to generate precipitation excess for each element during the simulation period. A two-dimensional diffusion wave model was used for overland flow routing and an iterative Alternative Direction Implicit scheme was used to solve the simultaneous overland flow equations. Once the overland flow became inflow to the channel, a one-dimensional dynamic wave flood routing technique, based on a four-point, implicit, non-linear finite difference solution of the St. Venant equation of unsteady flow, was applied. A limited number of comparisons were made between simulated and observed hydrographs for areas of about one square mile. Given the appropriate parameters, the model was able to accurately simulate runoff for single-event storms. This paper describes a distributed watershed model developed to simulate overland and channel flow. Comparisons were made between simulated and observed hydrographs for three watersheds. The model was able to accurately simulate the runoff for single-event storms using 61-m by 61-m (200 ft by 200 ft) watershed grid elements.     

HYDROLOGIC MODELING USING A SIMULATION LANGUAGE BASED ON QUEUEING THEORY1

ABSTRACT: A novel approach has been developed by applying queueing theory to hydrologic modeling. A queueing situation is characterized by a flow of customers arriving randomly at one or more service facilities. In this case, the customers are represented by water and the service time is the time it takes to move through the soil and over the land. This approach has the potential to be simpler and more efficient than some previously developed models. Another important attribute of queueing theory is its ability to model almost unlimited detail. Many simulation languages based on queueing theory have been developed and tested and are available. Applications of the languages include manufacturing operations, transportation systems, computer systems, financial planning, and health care systems. A model was developed, using a simulation language, to predict runoff hydrographs from storms occurring on small homogeneous watersheds. With continued development, queueing theory could provide an effident, detailed approach to simulating many natural processes.     

Effects of watershed configuration and composition on downstream lake water quality

We evaluated the relationships between landscape characteristics and lake water quality in receiving waters by regressing four water quality responses on landscape variables that were measured for whole watersheds and three different buffer distances (30, 60, and 120 m). Classical percolation theory was used to conceptualize nutrient pathways and to explain nonlinear responses. The response variables were total nitrogen (TN), total phosphorus (TP), chlorophyll-a (Chl-a), and Secchi transparency (SD). Landscape data were obtained from satellite image-derived maps of 130 watersheds in Iowa using geographic information systems software. We developed regression models with a stepwise protocol selecting the optimal number of significant explanatory variables. Configuration variables such as contagion, the cohesion of cropland and urban land, and the aggregation index of forest were very important and more important than variables assessing landscape composition (e.g., percentage farmland). Whole watershed models predicted between 15 and 67% of the variability in TN, TP, Chl-a, and SD. Proximity-explicit data offered only slightly improved statistical power over land cover data derived from the entire watershed for variables TN, Chl-a. and SD, but not for TP.