THE TWO‐DIMENSIONAL UPLAND EROSION MODEL CASC2D‐SED1

ABSTRACT: The two‐dimensional soil erosion model CASC2D‐SED simulates the dynamics of upland erosion during single rainstorms. The model is based on the raster‐based surface runoff calculations from CASC2D. Rainfall precipitation is distributed in time and space. Infiltration is calculated from the Green‐Ampt equations. Surface runoff is calculated from the diffusive wave approximation to the Saint‐Venant equations in two‐dimensions. Watershed data bases in raster Geographical Information System (GIS) provide information on the soil type, size fractions, soil erodibility, cropping management, and conservation practice factors for soil erosion calculations. Upland sediment transport is calculated for the size fractions (sand, silt, and clay), and the model displays the sediment flux, the amount of suspended sediment, and the net erosion and deposition using color graphics. The model has been tested on Goodwin Creek, Mississippi. The peak discharge and time to peak are within ± 20 percent and sediment transport rates within ?50 percent to 200 percent.     

RUNOFF SIMULATION USING RADAR RAINFALL DATA1

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

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.     

AN INTEGRATED ONE‐DIMENSIONAL AND TWO‐DIMENSIONAL URBAN STORMWATER FLOOD SIMULATION MODEL1

ABSTRACT: Flash flooding is the rapid flooding of low lying areas caused by the stormwater of intense rainfall associated with thunderstorms. Flash flooding occurs in many urban areas with relatively flat terrain and can result in severe property damage as well as the loss of lives. In this paper, an integrated one‐dimensional (1‐D) and two‐dimensional (2‐D) hydraulic simulation model has been established to simulate stormwater flooding processes in urban areas. With rainfall input, the model simulates 2‐D overland flow and 1‐D flow in underground stormwater pipes and drainage channels. Drainage channels are treated as special flow paths and arranged along one or more sides of a 2‐D computational grid. By using irregular computation grids, the model simulates unsteady flooding and drying processes over urban areas with complex drainage systems. The model results can provide spatial flood risk information (e.g., water depth, inundation time and flow velocity during flooding). The model was applied to the City of Beaumont, Texas, and validated with the recorded rainfall and runoff data from Tropical Storm Allison with good agreement.     

WATERSHED CONFIGURATION AND GEOGRAPHIC INFORMATION SYSTEM PARAMETERIZATION FOR SPUR MODEL HYDROLOGIC SIMULATIONS1

ABSTRACT: A grid cell geographic information system (GIS) is used to parameterize SPUR, a quasi-physically based surface runoff model in which a watershed is configured as a set of stream segments and contributing areas. GIS analysis techniques produce various watershed configurations by progressive simplification of a stream network delineated from digital elevation models (DEM). We used three watershed configurations: ≥ 2nd, ≥ 4th, and ≥ 13th Shreve order networks, where the watershed contains 28, 15, and 1 channel segments with 66, 37, and 3 contributing areas, respectively. Watershed configuration controls simulated daily and monthly sums of runoff volumes. For the climatic and topographic setting in southeastern Arizona the ≥ 4th order configuration of the stream network and contributing areas produces results that are typically as good as the ≥ 2nd order network. However both are consistently better than the ≥ 13th order configuration. Due to the degree of parameterization in SPUR, model simulations cannot be significantly improved by increasing watershed configuration beyond the ≥ 4th order network. However, a range of Soil Conservation Service curve numbers derived from rainfall/runoff data can affect model simulations. Higher curve numbers yield better results for the ≥ 2nd order network while lower curve numbers yield better results for the ≥ 4th order network.     

APPLICATION OF RORB MODEL TO A CATCHMENT IN SINGAPORE1

ABSTRACT: Runoff Routing model (RORB) is a general model applicable to both rural and urban catchments. The performance of the model is illustrated through its simulation of flood runoff hydrographs in an urban catchment in Singapore. The essential feature of the model is the routing of rainfall excesses on subareas through some arrangement of concentrated storage elements, which represent the distribution of temporary storage of flood runoff on the watershed. This nonlinear routing procedure of the storage elements has two common parameters, k_c and m. With the limited data available, these two parameter values were determined through calibration runs. The same set of values of k_c and m were then used in the model to determine the runoff hydrographs of five other storms selected from the rainfall events between 1979 and 1981. It was found that the simulated runoff hydrographs matched reasonably well with the recorded hydrographs.     

A GIS‐ASSISTED DISTRIBUTED WATERSHED MODEL FOR SIMULATING FLOODING AND INUNDATION1

ABSTRACT: A distributed watershed model combining kinematic wave routing, 1‐D dynamic channel‐flow routing, and 2‐D diffusive overland‐flow routing has been developed to simulate flooding and inundation levels of large watersheds. The study watershed was linked to a GIS database and was divided into an upstream mountainous area and a downstream alluvial plain. A kinematic wave routing was adopted at the mountainous area to compute the discharge flowing into the alluvial plain. A 1‐D dynamic channel routing solving the St. Venant equations by the Preissmann method was performed for the main channel of the alluvial plain, whereas a 2‐D overland‐flow routing solving the diffusion wave equation with the Alternating Direction Explicit scheme was used for floodplains. The above two routings were connected by weir‐link discharge formula. The parameters in the model were calibrated and independently verified by single‐event storms. An example application of flooding/inundation analysis was conducted for the Taichung station and the Woozi depot (Taiwan High Speed Rail). Suggested inundation‐proofing measures ‐ including raising ground surface elevation of the station and depot and building a waterproofing exterior wall and their combination ‐ were investigated. It was concluded that building the waterproofing exterior wall had a strong tendency to decrease peak inundation depth.     

STORM RUNOFF PREDICTION BASED ON A SPATIALLY DISTRIBUTED TRAVEL TIME METHOD UTILIZING REMOTE SENSING AND GIS1

ABSTRACT: In this study, remotely sensed data and geographic information system (GIS) tools were used to estimate storm runoff response for Simms Creek watershed in the Etonia basin in northeast Florida. Land cover information from digital orthophoto quarter quadrangles (DOQQ), and enhanced thematic mapper plus (ETM+) were analyzed for the years 1990, 1995, and 2000. The corresponding infiltration excess runoff response of the study area was estimated using the U.S. Department of Agriculture (USDA), Natural Resources Conservation Service Curve Number (NRCS‐CN) method. A digital elevation model (DEM)/GIS technique was developed to predict stream response to runoff events based on the travel time from each grid cell to the watershed outlet. A comparison of predicted to observed stream response shows that the model predicts the total runoff volume with an efficiency of 0.98, the peak flow rate at an efficiency of 0.85, and the full direct runoff hydrograph with an average efficiency of 0.65. The DEM/GIS travel time model can be used to predict the runoff response of ungaged watersheds and is useful for predicting runoff hydrographs resulting from proposed large scale changes in the land use.     

GIS‐Based Predictive Models of Hillslope Runoff Generation Processes1

Abstract: Successful nonpoint source pollution control using best management practice placement is a complex process that requires in‐depth knowledge of the locations of runoff source areas in a watershed. Currently, very few simulation tools are capable of identifying critical runoff source areas on hillslopes and those available are not directly applicable under all runoff conditions. In this paper, a comparison of two geographic information system (GIS)‐based approaches: a topographic index model and a likelihood indicator model is presented, in predicting likely locations of saturation excess and infiltration excess runoff source areas in a hillslope of the Savoy Experimental Watershed located in northwest Arkansas. Based on intensive data collected from a two‐year field study, the spatial distributions of hydrologic variables were processed using GIS software to develop the models. The likelihood indicator model was used to produce probability surfaces that indicated the likelihood of location of both saturation and infiltration excess runoff mechanisms on the hillslope. Overall accuracies of the likelihood indicator model predictions varied between 81 and 87% for the infiltration excess and saturation excess runoff locations respectively. On the basis of accuracy of prediction, the likelihood indicator models were found to be superior (accuracy 81‐87%) to the predications made by the topographic index model (accuracy 69.5%). By combining statistics with GIS, runoff source areas on a hillslope can be identified by incorporating easily determined hydrologic measurements (such as bulk density, porosity, slope, depth to bed rock, depth to water table) and could serve as a watershed management tool for identifying critical runoff source areas in locations where the topographic index or other similar methods do not provide reliable results.     

EFFECTS OF WATERSHED REPRESENTATION ON RUNOFF AND SEDIMENT YIELD MODELING1

ABSTRACT: This paper evaluates the effects of watershed geometric representation (i.e., plane and channel representation) on runoff and sediment yield simulations in a semiarid rangeland watershed. A process based, spatially distributed runoff erosion model (KINEROS2) was used to explore four spatial representations of a 4.4 ha experimental watershed. The most complex representation included all 96 channel elements identifiable in the field. The least complex representation contained only five channel elements. It was concluded that oversimplified watershed representations greatly influence runoff and sediment yield simulations by inducing excessive infiltration on hillslopes and distorting runoff patterns and sediment fluxes. Runoff and sediment yield decrease systematically with decreasing complexity in watershed representation. However, less complex representations had less impact on runoff and sediment‐yield simulations for small rainfall events. This study concludes that the selection of the appropriate level of watershed representation can have important theoretical and practical implications on runoff and sediment yield modeling in semiarid environments.     

MODEL PREDICTIONS OF WATERSHED HYDROLOGIC COMPONENTS: COMPARISON AND VERIFICATIION1

ABSTRACT: Model predictions of the relatively simple soil compartment model SESOIL are compared with those of the more data-intensive terrestrial ecosystem hydrology model AGTEHM. Comparisons were performed using data from a deciduous forest stand watershed, a grassland watershed, and two agricultural field plots. Good agreement was obtained between model predictions for annual values of infiltration, evapotranspiration, surface runoff, and groundwater runoff. SESOIL model predictions also compare well with empirical measurements at the forest stand and the grassland watersheds.     

A TRIANGULAR MODEL OF DIMENSIONLESS RUNOFF PRODUCING RAINFALL HYETOGRAPHS IN TEXAS1

ABSTRACT: A synthetic triangular hyetograph for a large data base of Texas rainfall and runoff is needed. A hyetograph represents the temporal distribution of rainfall intensity at a point or over a watershed during a storm. Synthetic hyetographs are estimates of the expected time distribution for a design storm and principally are used in small watershed hydraulic structure design. A data base of more than 1,600 observed cumulative hyetographs that produced runoff from 91 small watersheds (generally less than about 50 km²) was used to provide statistical parameters for a simple triangular shaped hyetograph model. The model provides an estimate of the average hyetograph in dimensionless form for storm durations of 0 to 24 hours and 24 to 72 hours. As a result of this study, the authors concluded that the expected dimensionless cumulative hyetographs of 0 to 12 hour and 12 to 24 hour durations were sufficiently similar to be combined with minimal information loss. The analysis also suggests that dimensionless cumulative hyetographs are independent of the frequency level or return period of total storm depth and thus are readily used for many design applications. The two triangular hyetographs presented are intended to enhance small watershed design practice in applicable parts of Texas.     

PARTIAL AREA RESPONSE ON SMALL SEMIARID WATERSHEDS1

ABSTRACT: Significant errors in estimating surface runoff and erosion rates are possible if a watershed is assumed to contribute runoff uniformly over the entire area, when actually only a portion of the entire area may be contributing. Generation of overland flow on portions of small semiarid watersheds was analyzed by three methods: an average loss rate procedure, a lumped-linear model, and a distributed-nonlinear model. These methods suggested that, on the average, 45, 60, and 50% of the drainage area was contributing runoff at the watershed outlet. Infiltrometer data support the partial area concept and indicate that the low infiltration zones are the runoff source areas as simulated with the distributed-nonlinear model.     

Hydrologic Modeling Uncertainty Resulting From Land Cover Misclassification1

Abstract: A stochastic, spatially explicit method for assessing the impact of land cover classification error on distributed hydrologic modeling is presented. One‐hundred land cover realizations were created by systematically altering the North American Landscape Characterization land cover data according to the dataset’s misclassification matrix. The matrix indicates the probability of errors of omission in land cover classes and is used to assess the uncertainty in hydrologic runoff simulation resulting from parameter estimation based on land cover. These land cover realizations were used in the GIS‐based Automated Geospatial Watershed Assessment tool in conjunction with topography and soils data to generate input to the physically‐based Kinematic Runoff and Erosion model. Uncertainties in modeled runoff volumes resulting from these land cover realizations were evaluated in the Upper San Pedro River basin for 40 watersheds ranging in size from 10 to 100 km² under two rainfall events of differing magnitudes and intensities. Simulation results show that model sensitivity to classification error varies directly with respect to watershed scale, inversely to rainfall magnitude and are mitigated or magnified by landscape variability depending on landscape composition.     

STORMFLOW SIMULATION USING A GEOGRAPHICAL INFORMATION SYSTEM WITH A DISTRIBUTED APPROACH1

ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.     

A QUASI-LINEAR SPATIALLY DISTRIBUTED CELL MODEL FOR THE SURFACE RUNOFF SYSTEM1

ABSTRACT: The paper presents a spatially distributed model consisting of cells that are interconnected in a pattern similar to the major drainage network of the watershed. Each cell receives as input the rainfall excess for the area represented by the cell as well as inflows from cells located upstream. Outflow from the cell is derived by routing the total input through the cell assuming it to be a linear reservoir during the storm. The time constant of the cells is however allowed to vary from storm to storm so that the model may be described as a quasi-linear model. The model was tested with rainfall excess and direct surface runoff data available for a medium size watershed with satisfactory results. The time constant was found to be related to the rainfall excess of the storms studied, its value decreasing with the increase in the total rainfall excess.     

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.     

ALAWAT: A spatially allocated watershed model for approximating stream,sediment, and pollutant flows in Hawaii,USA

The Ala Wai Canal Watershed Model (ALAWAT) is a planning-level watershed model for approximating direct runoff, streamflow, sediment loads, and loads for up to five pollutants. ALAWAT uses raster GIS data layers including land use, SCS soil hydrologic groups, annual rainfall, and subwatershed delineations as direct model parameter inputs and can use daily total rainfall from up to ten rain gauges and streamflow from up to ten stream gauges. ALAWAT uses a daily time step and can simulate flows for up to ten-year periods and for up to 50 subwatersheds. Pollutant loads are approximated using a user-defined combination of rating curve relationships, mean event concentrations, and loading/washoff parameters for specific subwatersheds, land uses, and times of year. Using ALAWAT, annual average streamflow and baseflow relationships and urban suspended sediment loads were approximated for the Ala Wai Canal watershed (about 10,400 acres) on the island of Oahu, Hawaii. Annual average urban suspended sediments were approximated using two methods: mean event concentrations and pollutant loading and washoff. Parameters for the pollutant loading and washoff method were then modified to simulate the effect of various street sweeping intervals on sediment loads.     

COMPOSITE VS. DISTRIBUTED CURVE NUMBERS: EFFECTS ON ESTIMATES OF STORM RUNOFF DEPTHS1

ABSTRACT: The U.S. Department of Agriculture Curve Number (CN) method is one of the most common and widely used techniques for estimating surface runoff and has been incorporated into a number of popular hydrologic models. The CN method has traditionally been applied using compositing techniques in which the area weighted average of all curve numbers is calculated for a watershed or a small number of sub-watersheds. CN compositing was originally developed as a time saving procedure, reducing the number of runoff calculations required. However, with the proliferation of high speed computers and geographic information systems, it is now feasible to use distributed CNs when applying the CN method. To determine the effect of using composited versus distributed CNs on runoff estimates, two simulations of idealized watersheds were developed to compare runoff depths using composite and distributed CNs. The results of these simulations were compared to the results of similar analyses performed on an urbanizing watershed located in central Indiana and show that runoff depth estimates using distributed CNs are as much as 100 percent higher than when composited CNs are used. Underestimation of runoff due to CN compositing is a result of the curvilinear relationship between CN and runoff depth and is most severe for wide CN ranges, low CN values, and low precipitation depths. For larger design storms, however, the difference in runoff computed using composite and distributed CNs is minimal.     

RUNOFF RESPONSE FROM TWO RECLAIMED WATERSHEDS1

ABSTRACT: Detailed studies of the surface hydrology of reclaimed surface-mined watersheds for both rainfall and snowmelt events are non-existent for central Alberta yet this information is crucial for design of runoff conveyance and storage structures. A study was initiated in 1992 with principal objectives of quantifying surface runoff for both summer rainfall and spring snowmelt events and identifying the dominant flow processes occurring in two reclaimed watersheds. Snowmelt accounted for 86 and 100% of annual watershed runoff in 1993 and 1994, respectively. The highest instantaneous peak flow was recorded during a summer rainfall event with a return period of greater than 50 years. Infiltration-excess overland flow was identified as the dominant flow process occurring within the Sandy Subsoil Watershed, whereas saturation overland flow was the principal runoff process occurring within the West Watershed.