CHANNEL STABILITY DOWNSTREAM FROM A DAM ASSESSED USING AERIAL PHOTOGRAPHS AND STREAM‐GAGE INFORMATION1

ABSTRACT: The stability of the Neosho River channel downstream from John Redmond Dam, in southeast Kansas, was investigated using multiple‐date aerial photographs and stream‐gage information. Bankfull channel width was used as the primary indicator variable to assess pre‐ and post‐dam channel change. Five sin‐mile river reaches and four stream gages were used in the analysis. Results indicated that, aside from some localized channel widening, the overall channel change has been minor with little post‐dam change in bankfull channel width. The lack of a pronounced post‐dam channel change may be attributed to a substantial reduction in the magnitude of the post‐dam annual peak discharges in combination with the resistance to erosion of the bed and bank materials. Also, the channel may have been overwidened by a series of large floods that predated construction of the dam, including one with an estimated 500‐year recurrence interval.     

ASSESSING STREAM CHANNEL STABILITY THRESHOLDS USING FLOW COMPETENCE ESTIMATES AT BANKFULL STAGE1

ABSTRACT: Stream channel stability is affected by peak flows rather than average annual water yield. Timber harvesting and other land management activities that contribute to soil compaction, vegetation removal, or increased drainage density can increase peak discharges and decrease the recurrence interval of bankfull discharges. Increased peak discharges can cause more frequent movement of large streambed materials, leading to more rapid stream channel change or instability. This study proposes a relationship between increased discharge and channel stability, and presents a methodology that can be used to evaluate stream channel stability thresholds on a stream reach basis. Detailed surveys of the channel cross section, water surface slope, streambed particle size distribution, and field identification of bankfull stage are used to estimate existing bankfull flow conditions. These site specific stream channel characteristics are used in bed load movement formulae to predict critical flow conditions for entrainment of coarse bed material (D₈₄ size fraction). The “relative bed stability” index, defined as the ratio of critical flow condition to the existing condition at bankfull discharge, can predict whether increased peak discharges will exceed stream channel thresholds.     

STREAM-CHANNEL INCISION FOLLOWING DRAINAGE-BASIN URBANIZATION1

ABSTRACT: Urbanization of a drainage basin results in pervasive hydrologic changes that in turn initiate long-term changes in stream channels. Increases in peak discharges and in durations of high flows result in either quasi-equilibrium channel expansion, where cross-section area increases in near-proportion to the discharge increase, or catastrophic channel incision, where changes occur far out of proportion to the discharge increases that initiated them. Field data and hydrologic modeling of rapidly urbanizing basins in King County, Washington, define conditions of flow, topography, geology, and channel roughness that identify streams susceptible to incision. Channel slope and geologic material are particularly critical; thus simple map overlays, nearly irrespective of contributing drainage area, provide a valuable planning tool for identification of susceptible terrain. Where such conditions exist, basal shear stress provides a quantifiable parameter for predicting likely problems, although knickpoints are typical in such settings and confound simple calculation of sediment-transport rates. Where urbanization proceeds in such areas, effective mitigation of the incision hazards requires a degree of stormwater control far in excess of standards typically applied to present development activity.     

EFFECTS OF DAM IMPOUNDMENT ON THE FLOOD REGIME OF NATURAL FLOODPLAIN COMMUNITIES IN THE UPPER CONNECTICUT RIVER1

ABSTRACT: Understanding the effects of dams on the inundation regime of natural floodplain communities is critical for effective decision making on dam management or dam removal. To test the implications of hydrologic alteration by dams for floodplain natural communities, we conducted a combined field and modeling study along two reaches in the Connecticut River Rapids Macrosite (CRRM), one of the last remaining flowing water sections of the Upper Connecticut River. We surveyed multiple channel cross sections at both locations and concurrently identified and surveyed the elevations of important natural communities, native species of concern, and nonnative invasive species. Using a hydrologic model, HEC‐RAS, we routed estimated pre‐and post‐impoundment discharges of different design recurrence intervals (two year through 100 year floods) through each reach to establish corresponding reductions in elevation and effective wetted perimeter following post‐dam discharge reductions. By comparing (1) the frequency and duration of flooding of these surfaces before and after impoundment and (2) the total area flooded at different recurrence intervals, our goal was to derive a spatially explicit assessment of hydrologic alteration, directly relevant to natural floodplain communities. Post‐impoundment hydrologic alteration profoundly affected the subsequent inundation regime, and this impact was particularly true of higher floodplain terraces. These riparian communities, which were flooded, on average, every 20 to 100 years pre‐impoundment, were predicted to flood at 100 ? 100 year intervals, essentially isolating them completely from riverine influence. At the pre‐dam five to ten year floodplain elevations, we observed smaller differences in predicted flood frequency but substantial differences in the total area flooded and in the average flood duration. For floodplain forests in the Upper Connecticut River, this alteration by impoundment suggests that even if other stresses facing these communities (human development, invasive exotics) were alleviated, this may not be sufficient to restore intact natural communities. More generally, our approach provides a way to combine site specific variables with long term gage records in assessing the restorative potential of dam removal.     

THE CONCEPT OF HYDROLOGIC LANDSCAPES1

ABSTRACT: Hydrologic landscapes are multiples or variations of fundamental hydrologic landscape units. A fundamental hydrologic landscape unit is defined on the basis of land‐surface form, geology, and climate. The basic land‐surface form of a fundamental hydrologic landscape unit is an upland separated from a lowland by an intervening steeper slope. Fundamental hydrologic landscape units have a complete hydrologic system consisting of surface runoff, ground‐water flow, and interaction with atmospheric water. By describing actual landscapes in terms of land‐surface slope, hydraulic properties of soils and geologic framework, and the difference between precipitation and evapotranspiration, the hydrologic system of actual landscapes can be conceptualized in a uniform way. This conceptual framework can then be the foundation for design of studies and data networks, syntheses of information on local to national scales, and comparison of process research across small study units in a variety of settings. The Crow Wing River watershed in central Minnesota is used as an example of evaluating stream discharge in the context of hydrologic landscapes. Lake‐research watersheds in Wisconsin, Minnesota, North Dakota, and Nebraska are used as an example of using the hydrologic‐land‐scapes concept to evaluate the effect of ground water on the degree of mineralization and major‐ion chemistry of lakes that lie within ground‐water flow systems.     

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.     

Modeling Hydrologic Benefits of Low Impact Development: A Distributed Hydrologic Model of The Woodlands,Texas

Low Impact Development (LID) is alternative design approach to land development that conserves and utilizes natural resources to minimize the potential negative environmental impacts of development, such as flooding. The Woodlands near Houston, Texas is one of the premier master‐planned communities in the United States. Unlike in a typical urban development where riparian corridors are often replaced with concrete channels, pervious surfaces, vegetation, and natural drainage pathways were preserved as much as possible during development. In addition, a number of detention ponds were strategically located to manage runoff on site. This article uses a unique distributed hydrologic model, Vflo?, combined with historical (1974) and recent (2008 and 2009) rainfall events to evaluate the long‐term effectiveness of The Woodlands natural drainage design as a stormwater management technique. This study analyzed the influence of LID in The Woodlands by comparing the hydrologic response of the watershed under undeveloped, developed, and highly urbanized conditions. The results show that The Woodlands drainage design successfully reflects predeveloped hydrologic conditions and produces peak flows two to three times lower than highly urbanized development. Furthermore, results indicate that the LID practices employed in The Woodlands successfully attenuate the peak flow from a 100‐year design event, resulting in flows comparable to undeveloped hydrologic conditions.     

Measuring the Impact of Urbanization on Channel Widths Using Historic Aerial Photographs and Modern Surveys1

Abstract: Land use in a watershed is commonly held to exert a strong influence on trunk channel form and process. Land use changes act over human time‐scales, which are short enough to measure effects on channels directly using historic aerial photographs. We show that high‐resolution topographic surveys for the channels of paired watersheds in the Lehigh Valley, Pennsylvania, are comparable, but have channel widths that have changed dramatically in the past five decades. The two watersheds, Little Lehigh Creek and Sacony Creek, are similar in most aspects except in their respective amount of urban land use. Aerial photographs of the urbanized Little Lehigh Creek show that a majority of the measured widths (67 of 85) were statistically wider in 1999 than in 1947. In contrast, the measured widths from the agricultural Sacony Creek are more evenly distributed among those that widened (18), narrowed (28), and those that were statistically unchanged (6) from 1946 to 1999. From 1946 to 1999 the only section of Sacony Creek that widened was that reach downstream of the only sizable urban area in the watershed. The current land use in Sacony Creek watershed resembles that of 1946, while the Little Lehigh Creek watershed has more than tripled its urban area. These data, in concert with other recent hydrologic data from the watersheds suggest that the increase in urban area‐generated peak discharges is the mechanism behind the widening that occurred in the Little Lehigh Creek. These wider channels can affect water quality, aquatic habitat, suspended sediment loads, and river esthetics.     

Impact of Changing Climate and Land Cover on Flood Magnitudes in the Delaware River Basin,USA

Changing climate and land cover are expected to impact flood hydrology in the Delaware River Basin over the 21st Century. HEC‐HMS models (U.S. Army Corps of Engineers Hydrologic Engineering Center‐Hydrologic Modeling System) were developed for five case study watersheds selected to represent a range of scale, soil types, climate, and land cover. Model results indicate that climate change alone could affect peak flood discharges by ?6% to +58% a wide range that reflects regional variation in projected rainfall and snowmelt and local watershed conditions. Land cover changes could increase peak flood discharges up to 10% in four of the five watersheds. In those watersheds, the combination of climate and land cover change increase modeled peak flood discharges by up to 66% and runoff volumes by up to 44%. Precipitation projections are a key source of uncertainty, but there is a high likelihood of greater precipitation falling on a more urbanized landscape that produces larger floods. The influence of climate and land cover changes on flood hydrology for the modeled watersheds varies according to future time period, climate scenario, watershed land cover and soil conditions, and flood frequency. The impacts of climate change alone are typically greater than land cover change but there is substantial geographic variation, with urbanization the greater influence on some small, developing watersheds.     

SPATIALLY EXPLICIT HYDROLOGIC MODELING OF LAND USE CHANGE1

ABSTRACT: Using a Geographic Information System (GIS), a method is presented to develop a spatially explicit time series of land use in an urbanizing watershed. The method is prefaced on the existence of independent observations of land use at different times and data that describes the spatial‐temporal land use transition characteristics of the watershed between these two points in time. A method is then presented to generalize the TR‐55 graphical method, a common lumped hydrologic model for estimating peak discharge, for use in a spatially explicit scheme. This scheme predicts peak discharge throughout a watershed, rather than at a single selected watershed outlet. Coupling these two methods allows the engineer to model both the temporal and spatial evolution of peak discharge for the watershed. An illustrative watershed in a suburban area of Washington, DC is selected to demonstrate the methods. The model results from these analyses are presented graphically to highlight the complex features in peak discharge behavior that exist both spatially, as a function of position within the watershed drainage network, and temporally, as the watershed undergoes urbanization. These features are not commonly noted in most hydrologic analyses but are captured in these analyses because of the high spatial and temporal resolution of the methods presented. The physical implications of the modeled results are discussed in the context of the information content of a stream gauge located at the overall outlet of the illustrative watershed. This work shows that the common practice of transposition of gauge information to locations internal to the watershed would neglect internal variability in peak discharge behavior, and could potentially lead to the determination of inappropriate design discharges.     

A Framework to Develop Nationwide Flooding Extents Using Climate Models and Assess Forecast Potential for Flood Resilience

The methods used to simulate flood inundation extents can be significantly improved by high‐resolution spatial data captured over a large area. This paper presents a hydraulic analysis methodology and framework to estimate national‐level floodplain changes likely to be generated by climate change. The hydraulic analysis was performed using existing published Federal Emergency Management Agency 100‐year floodplains and estimated 100‐ and 10‐year return period peak flow discharges. The discharges were estimated using climate variables from global climate models for two future growth scenarios: Representative Concentration Pathways 2.6 and 8.5. River channel dimensions were developed based on existing regional United States Geological Survey publications relating bankfull discharges with channel characteristics. Mathematic relationships for channel bankfull topwidth, depth, and side slope to contributing drainage area measured at model cross sections were developed. The proposed framework can be utilized at a national level to identify critical areas for flood risk assessment. Existing hydraulic models at these “hot spots” could be repurposed for near–real‐time flood forecasting operations. Revitalizing these models for use in simulating flood scenarios in near–real time through the use of meteorological forecasts could provide useful information for first responders of flood emergencies.     

LONG‐TERM CHANGES IN REGIONAL HYDROLOGIC REGIME FOLLOWING IMPOUNDMENT IN A HUMID‐CLIMATE WATERSHED1

ABSTRACT: We analyzed the type of hydrologic adjustments resulting from flow regulation across a range of dam types, distributed throughout the Connecticut River watershed, using two approaches: (1) the Index of Hydrologic Alteration (IHA) and (2) log‐Pearson Type III flood frequency analysis. We applied these analyses to seven rivers that have extensive pre‐and post‐disturbance flow records and to six rivers that have only long post‐regulation flow records. Lastly, we analyzed six unregulated streams to establish the regional natural flow regime and to test whether it has changed significantly over time in the context of an increase in forest cover from less than 20 percent historically to greater than 80 percent at present. We found significant hydrologic adjustments associated with both impoundments and land use change. On average, maximum peak flows decrease by 32 percent in impounded rivers, but the effect decreases with increasing flow duration. One‐day minimum low flows increase following regulation, except for the hydro‐electric facility on the mainstem. Hydrograph reversals occur more commonly now on the mainstem, but the tributary flood control structures experience diminished reversals. Major shifts in flood frequency occur with the largest effect occurring downstream of tributary flood control impoundments and less so downstream of the mainstem's hydroelectric facility. These overall results indicate that the hydrologic impacts of dams in humid environments can be as significant as those for large, multiple‐purpose reservoirs in more arid environments.     

Bankfull Regional Curves for North and Northwest Florida Streams1

Abstract: Regional curves, which relate bankfull channel dimensions and discharge to watershed drainage area, are developed to aid in identifying the bankfull stage in ungaged watersheds, and estimating the bankfull discharge and dimensions for river studies and natural channel design applications. This study assessed 26 stable stream reaches in two hydro‐physiographic regions of the Florida Coastal Plain: the Northwest Florida Coastal Plain (NWFCP) and the North Florida Coastal Plain (NFCP). Data from stream reaches in Georgia and Alabama were also used to develop the Florida regional curves, since they are located in the same hydro‐physiographic region. Reaches were selected based on the presence of U.S. Geological Survey gage stations and indicators of limited watershed development (e.g., <10% impervious surface). Analyses were conducted to determine bankfull channel dimensions, bankfull discharge, average channel slope, and Rosgen stream classification. Based on these data, significant relationships were found between bankfull cross‐sectional area, width, mean depth, and discharge as a function of drainage area for both regions. Data from this study suggested that bankfull discharges and channel dimensions were larger from NWFCP streams than from Coastal Plain streams in North Carolina and Maryland. Bankfull discharges were similar between NFCP and Georgia coastal plain streams; therefore, the data were combined into one regional curve. In addition, the data were stratified by Rosgen stream type. This stratification strengthened the relationships of bankfull width and mean depth as a function of drainage area.     

EFFECTS OF SCALE ON LAND USE AND WATER QUALITY RELATIONSHIPS: A LONGITUDINAL BASIN‐WIDE PERSPECTIVE1

ABSTRACT: Human land use is a major source of change in catchments in developing areas. To better anticipate the long‐term effects of growth, land use planning requires estimates of how changes in land use will affect ecosystem processes and patterns across multiple scales of space and time. The complexity of biogeochemical and hydrologic interactions within a basin makes it difficult to scale up from process‐based studies of individual reaches to watershed scales over multiple decades. Empirical models relating land use/land cover (LULC) to water quality can be useful in long‐term planning, but require an understanding of the effects of scale on apparent land use‐water quality relationships. We empirically determined how apparent relationships between water quality and LULC data change at different scales, using LIJLC data from the Willapa Bay watershed (Washington) and water quality data collected along the Willapa and North Rivers. Spatial scales examined ranged from the local riparian scale to total upstream catchment. The strength of the correlations between LTJLC data and longitudinal water quality trends varied with scale. Different water quality parameters also varied in their response to changes in scale. Intermediate scales of land use data generally were better predictors than local riparian or total catchment scales. Additional data from the stream network did not increase the strength of relationships significantly. Because of the likelihood of scale‐induced artifacts, studies quantifying land use‐water quality relationships performed at single scales should be viewed with great caution.     

Habitat Improvements and Fish Community Response Associated with an Agricultural Two‐Stage Ditch in Mower County,Minnesota

Water quality and stream habitat in agricultural watersheds are under greater scrutiny as hydrologic pathways are altered to increase crop production. Ditches have been traditionally constructed to remove water from agricultural lands. Little attention has been placed on alternative ditch designs that are more stable and provide greater habitat diversity for wildlife and aquatic species. In 2009, 1.89 km of a conventional drainage ditch in Mower County, Minnesota, was converted to a two‐stage ditch (TSD) with small, adjacent floodplains to mimic a natural system. Cross section surveys, conducted pre‐ and post‐construction, generally indicate a stable channel with minor adjustments over time. Vegetation surveys showed differences in species composition and biomass between the slopes and the benches, with changes ongoing. Longitudinal surveys demonstrated a 12‐fold increase in depth variability. Fish habitat quality improved with well‐sorted gravel riffles and deeper pool habitat. The biological response to improved habitat quality was investigated using a Fish Index of Biological Integrity (FIBI). Our results show higher FIBI scores post‐construction with scores more similar to natural streams. In summary, the TSD demonstrated improvements in riparian and instream habitat quality and fish communities, which showed greater fish species richness, higher percentages of gravel spawning fish, and better FIBI scores. This type of management tool could benefit ditches in other regions where gradient and geology allow.     

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.     

Analyses of Urban Drainage Network Structure and its Impact on Hydrologic Response1

Meierdiercks, Katherine L., James A. Smith, Mary Lynn Baeck, and Andrew J. Miller, 2010. Analyses of Urban Drainage Network Structure and Its Impact on Hydrologic Response. Journal of the American Water Resources Association (JAWRA) 1-12. DOI: 10.1111/j.1752-1688.2010.00465.x Abstract: Urban flood studies have linked the severity of flooding to the percent imperviousness or land use classifications of a watershed, but relatively little attention has been given to the impact of urban drainage networks on hydrologic response. The drainage network, which can include storm pipes, surface channels, street gutters, and stormwater management ponds, is examined in the Dead Run watershed (14.3 km²). Comprehensive digital representations of the urban drainage network in Dead Run were developed and provide a key observational resource for analyses of urban drainage networks and their impact on hydrologic response. Analyses in this study focus on three headwater subbasins with drainage areas ranging from 1.3 to 1.9 km² and that exhibit striking contrasts in their patterns and history of development. It is shown that the drainage networks of the three subbasins, like natural river networks, exhibit characteristic structures and that these features play critical roles in determining urban hydrologic response. Hydrologic modeling analyses utilize the Environmental Protection Agency’s Stormwater Management Model (SWMM), which provides a flexible platform for examining the impacts of drainage network structure on hydrologic response. Results of SWMM modeling analyses suggest that drainage density and presence of stormwater ponds impact peak discharge more significantly in the Dead Run subbasins than the percent impervious or land use type of the subbasins.     

Hydrologic Impact Assessment of Land Cover Change and Stormwater Management Using the Hydrologic Footprint Residence

Urbanization impacts the stormwater regime through increased runoff volumes and velocities. Detention ponds and low impact development (LID) strategies may be implemented to control stormwater runoff. Typically, mitigation strategies are designed to maintain postdevelopment peak flows at predevelopment levels for a set of design storms. Peak flow does not capture the extent of changes to the hydrologic flow regime, and the hydrologic footprint residence (HFR) was developed to calculate the area and duration of inundated land during a storm. This study couples a cellular automata land cover change model with a hydrologic and hydraulic framework to generate spatial projections of future development on the fringe of a rapidly urbanizing metropolitan area. The hydrologic flow regime is characterized for existing and projected land cover patterns under detention pond and LID‐based control, using the HFR and peak flow values. Results demonstrate that for less intense and frequent rainfall events, LID solutions are better with respect to HFR; for larger storms, detention pond strategies perform better with respect to HFR and peak flow.     

SIMULATION OF FRESHWATER DISCHARGES FROM UNGAGED AREAS TO THE SEBASTIAN RIVER,FLORIDA1

ABSTRACT: Due to alterations in the natural drainage system over the past several decades and intensified agricultural practices, freshwater discharges to the Sebastian River, Florida, have increased substantially. As a result, salinity patterns in the Sebastian River and adjacent Indian River lagoon have been disrupted and the influx of nutrients has increased. Recently, the St. Johns River Water Management District has developed a 3‐D hydrodynamic and salinity model for the Sebastian River and adjacent Indian River to study the effects of freshwater inflows, and to set guidelines for management of future freshwater discharges. Freshwater inflows to the Sebastian River are part of the input data of the hydrodynamic model. Except for the downstream drainage areas, inflows are gaged, and the data were used for calibration of the hydrologic simulations. Collectively, the downstream ungaged areas constitute about 16 percent of the total drainage area. Because of the significant contribution to the total drainage area, reliable estimates of freshwater discharges from the ungaged areas to the Sebastian River are needed. This case study illustrates the development of a set of model parameters, reflecting the hydrologic and physiographic characteristics of the entire region. In this context region applies to the watersheds located in the coastal area along the Indian River from Titusville in the north to Vero Beach in the south. The parameter set was first tested on a number of gaged drainage basins in the region, and was then applied to the ungaged areas.     

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.