SPATIALLY DISTRIBUTED MODELING OF STREAM FLOW DURING STORM EVENTS1

ABSTRACT: The purpose of this paper is to present a new approach for the spatially distributed modeling of water flow during storm events. Distributed modeling of flow during storm events is an important basis for any environmental modeling, including turbidity or sediment transport. During the initial phase of a rainstorm, surface runoff is the main contributor of flow. To provide the spatial components for distributed hydrological modeling a Geographic Information System (GIS) was used to map and visualize contributing areas around a stream channel. Stream segments were defined using the hydrologic response unit (HRU) concept. Lateral flows were derived from GIS output for each segment of the stream and at each time interval of the rain storm and were routed using the kinematic routing equation. This approach is new in hydrological modeling and can be used to enhance many existing simulations. The model is also unique in the fine time scale (i.e., intervals are on the order of minutes). Model results showed good correlation with measured discharge values; however, further studies of contributing area behavior, its relationship with soil types and slope categories, and the influence of watershed size are needed to improve model performance. This model will be used in the future as the basis to model turbidity in streams.     

Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed1

Macholl, Jacob A., Katherine A. Clancy, and Paul M. McGinley, 2011. Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed. Journal of the American Water Resources Association (JAWRA) 47(1):114‐125. DOI: 10.1111/j.1752‐1688.2010.00495.x Abstract: Glaciated watersheds are not easily delineated using geographic information systems’ elevation‐based algorithms, especially where stream networks are disconnected and there are large regions of internally drained areas. This paper presents the results of an analysis using the Potential Contributing Source Area (PCSA) model to identify potential contributing areas, defined as areas from which runoff is physically capable of reaching a drainage network. The investigation was conducted to define the potential contributing areas in a glaciated region of northwest Wisconsin. The curve number (CN) method was used to predict runoff volumes in the watershed. The streamflows of four tributaries were measured and the runoff portion of the hydrograph quantified to be compared with runoff estimates calculated using the potential contributing areas and the traditional catchment area. Runoff producing events occurred, but the use of area‐weighted CN values was unsuccessful in modeling runoff due to all precipitation depths during the study period falling below the initial abstraction. A distributed CN approach provided runoff estimates that were generally better using the potential contributing areas compared with using the traditional catchment area. The extent of the minimum contributing area, estimated for a range of precipitation events, was found to be substantially less than the potential contributing areas, suggesting that the PCSA model delimits the maximum boundary of potential contributing 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.     

HYDROLOGIC PARAMETERIZATION OF WATERSHEDS FOR RUNOFF PREDICTION USING SWAT1

ABSTRACT: The use of continuous time, distributed parameter hydrologic models like SWAT (Soil and Water Assessment Tool) has opened several opportunities to improve watershed modeling accuracy. However, it has also placed a heavy burden on users with respect to the amount of work involved in parameterizing the watershed in general and in adequately representing the spatial variability of the watershed in particular. Recent developments in Geographical Information Systems (GIS) have alleviated some of the difficulties associated with managing spatial data. However, the user must still choose among various parameterization approaches that are available within the model. This paper describes the important parameterization issues involved when modeling watershed hydrology for runoff prediction using SWAT with emphasis on how to improve model performance without resorting to tedious and arbitrary parameter by parameter calibration. Synthetic and actual watersheds in Indiana and Mississippi were used to illustrate the sensitivity of runoff prediction to spatial variability, watershed decomposition, and spatial and temporal adjustment of curve numbers and return flow contribution. SWAT was also used to predict stream runoff from actual watersheds in Indiana that have extensive subsurface drainage. The results of this study provide useful information for improving SWAT performance in terms of stream runoff prediction in a manner that is particularly useful for modeling ungaged watersheds wherein observed data for calibration is not available.     

Assessing Critical Source Areas in Watersheds for Conservation Buffer Planning and Riparian Restoration

A science-based geographic information system (GIS) approach is presented to target critical source areas in watersheds for conservation buffer placement. Critical source areas are the intersection of hydrologically sensitive areas and pollutant source areas in watersheds. Hydrologically sensitive areas are areas that actively generate runoff in the watershed and are derived using a modified topographic index approach based on variable source area hydrology. Pollutant source areas are the areas in watersheds that are actively and intensively used for such activities as agricultural production. The method is applied to the Neshanic River watershed in Hunterdon County, New Jersey. The capacity of the topographic index in predicting the spatial pattern of runoff generation and the runoff contribution to stream flow in the watershed is evaluated. A simple cost-effectiveness assessment is conducted to compare the conservation buffer placement scenario based on this GIS method to conventional riparian buffer scenarios for placing conservation buffers in agricultural lands in the watershed. The results show that the topographic index reasonably predicts the runoff generation in the watershed. The GIS-based conservation buffer scenario appears to be more cost-effective than the conventional riparian buffer scenarios.     

RASTER-BASED HYDROLOGIC MODELING OF SPATIALLY-VARIED SURFACE RUNOFF1

ABSTRACT: The proliferation of watershed databases in raster Geographic Information System (GIS) format and the availability of radar-estimated rainfall data foster rapid developments in raster-based surface runoff simulations. The two-dimensional physically-based rainfall-runoff model CASC2D simulates spatially-varied surface runoff while fully utilizing raster GIS and radar-rainfall data. The model uses the Green and Ampt infiltration method, and the diffusive wave formulation for overland and channel flow routing enables overbank flow storage and routing. CASC2D offers unique color capabilities to display the spatio-temporal variability of rainfall, cumulative infiltrated depth, and surface water depth as thunderstorms unfold. The model has been calibrated and independently verified to provide accurate simulations of catchment response to moving rainstorms on watersheds with spatially-varied infiltration. The model can accurately simulate surface runoff from flashfloods caused by intense thunderstorms moving across partial areas of a watershed.     

Assessing the Hydrologic and Water Quality Benefits of a Network of Stormwater Control Measures in a SE U.S. Piedmont Watershed

The hydrologic and water quality benefits of an existing engineered stormwater control measures (SCMs) network, along with the alternative stormwater control simulations, were assessed in the rapidly urbanizing Beaverdam Creek watershed located in SE U.S. Piedmont region through the use of distributed Model of Urban Stormwater Improvement Conceptualization stormwater model. When compared with predevelopment conditions, the postdevelopment watershed simulation without SCMs indicated a 2 times increase in total runoff volume, 3 times average increase in peak flow for 1.5‐3.2 cm 6‐h storm events, and 30 times, 12 times, and 3 times higher total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN) loadings, respectively. The existing SCMs network, in comparison with the postdeveloped watershed without SCMs, reduced the average peak flow rates for 1.5‐3.2 cm 6‐h storm events by 70%, lowered the annual runoff volume by 3%, and lowered TSS, TP, TN annual loads by 57, 51, and 10%, respectively. A backyard rain garden simulation resulted in minimal additional reduction in TSS (1.6%), TP (0.4%), and TN (4%). Model simulations indicate that mandatory 85% TSS and 70% TP annual load reductions in comparison with the predevelopment levels would require the diversion of runoff from at least 70% of the contributing drainage areas runoff into additional offline bioretention basins.     

WATERSHED SCALE MODELING OF CRITICAL SOURCE AREAS OF RUNOFF GENERATION AND PHOSPHORUS TRANSPORT1

ABSTRACT: A curve number based model, Soil and Water Assessment Tool (SWAT), and a physically based model, Soil Moisture Distribution and Routing (SMDR), were applied in a headwater watershed in Pennsylvania to identify runoff generation areas, as runoff areas have been shown to be critical for phosphorus management. SWAT performed better than SMDR in simulating daily streamflows over the four‐year simulation period (Nash‐Sutcliffe coefficient: SWAT, 0.62; SMDR, 0.33). Both models varied streamflow simulations seasonally as precipitation and watershed conditions varied. However, levels of agreement between simulated and observed flows were not consistent over seasons. SMDR, a variable source area based model, needs further improvement in model formulations to simulate large peak flows as observed. SWAT simulations matched the majority of observed peak flow events. SMDR overpredicted annual flow volumes, while SWAT underpredicted the same. Neither model routes runoff over the landscape to water bodies, which is critical to surface transport of phosphorus. SMDR representation of the watershed as grids may allow targeted management of phosphorus sources. SWAT representation of fields as hydrologic response units (HRUs) does not allow such targeted management.     

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.     

Assessing Watershed-Scale, Long-Term Hydrologic Impacts of Land-Use Change Using a GIS-NPS Model

Land-use change, dominated by an increase in urban/impervious areas, has a significant impact on water resources. This includes impacts on nonpoint source (NPS) pollution, which is the leading cause of degraded water quality in the United States. Traditional hydrologic models focus on estimating peak discharges and NPS pollution from high-magnitude, episodic storms and successfully address short-term, local-scale surface water management issues. However, runoff from small, low-frequency storms dominates long-term hydrologic impacts, and existing hydrologic models are usually of limited use in assessing the long-term impacts of land-use change. A long-term hydrologic impact assessment (L-THIA) model has been developed using the curve number (CN) method. Long-term climatic records are used in combination with soils and land-use information to calculate average annual runoff and NPS pollution at a watershed scale. The model is linked to a geographic information system (GIS) for convenient generation and management of model input and output data, and advanced visualization of model results. The L-THIA/NPS GIS model was applied to the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana, USA. Historical land-use scenarios for 1973, 1984, and 1991 were analyzed to track land-use change in the watershed and to assess impacts on annual average runoff and NPS pollution from the watershed and its five subbasins. For the entire watershed between 1973 and 1991, an 18% increase in urban or impervious areas resulted in an estimated 80% increase in annual average runoff volume and estimated increases of more than 50% in annual average loads for lead, copper, and zinc. Estimated nutrient (nitrogen and phosphorus) loads decreased by 15% mainly because of loss of agricultural areas. The L-THIA/NPS GIS model is a powerful tool for identifying environmentally sensitive areas in terms of NPS pollution potential and for evaluating alternative land use scenarios for NPS pollution management.     

RUNOFF VOLUME ESTIMATION USING GIS TECHNIQUES1

ABSTRACT: Rainfall runoff of six watersheds was modeled via the Soil Conservation Service runoff curve number model in two ways: conventionally (manually) and via a geographic information system (GIS). Input data (elevation, soils, and landcover) were digital for the latter method. In contrast to previous studies, the GIS was ised for all phases of the modeling process, including watershed delineation and routing of runoff. A comparison between the two methods was consistent with results reported by others and indicates that the use of a GIS is an acceptable alternative to the conventional method for watersheds lacking relatively flat terrain. Given this limitation, the GIS method may prove advantageous over manual methods when study areas are large or numerous, runoff is modeled repetitively, alternative landcover scenarios are explored, or a digital database already exists for the study area.     

SUSCEPTIBILITY OF INDIANA WATERSHEDS TO HERBICIDE CONTAMINATION1

ABSTRACT: An index of watershed susceptibility to surface water contamination by herbicides could be used to improve source water assessments for public drinking water supplies, prioritize watershed restoration projects, and direct funding and educational efforts to areas where the greatest environmental benefit can be realized. The goal of this study is to use streamflow and herbicide concentration data to develop and evaluate a method for estimating comparative watershed susceptibility to herbicide loss. United States Geological Survey (USGS) concentration data for five relatively water soluble herbicides (alachlor, atrazine, cyanazine, metolachlor, and simazine) were analyzed for 16 Indiana watersheds. Correlation was assessed between observed herbicide losses and: (1) a herbicide runoff index using GIS‐based land use, soil type, SCS runoff curve number, tillage practice, herbicide use estimates, and combinations of these factors; and (2) predicted herbicide losses from a non‐point source pollution model (NAPRA‐Web, an Internet‐based interface for GLEAMS). The highest adjusted R²value was found between herbicide concentration and the runoff curve number alone, ranging from 0.25 to 0.56. Predictions from the simulation model showed a poorer correlation with observed herbicide loss. This indicates potential for using the runoff curve number as a simple herbicide contamination susceptibility index.     

A FAST ALGORITHM FOR AUTOMATICALLY COMPUTING STRAHLER STREAM ORDER1

ABSTRACT: An efficient algorithm was developed to determine Strahler stream order for segments of stream networks represented in a Geographic Information System (GIS). The algorithm correctly assigns Strahler stream order in topologically complex situations such as braided streams and multiple drainage outlets. Execution time varies nearly linearly with the number of stream segments in the network. This technique is expected to be particularly useful for studying the topology of dense stream networks derived from digital elevation model data.     

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.     

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.     

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.     

WATERSHED WEIGHTING OF EXPORT COEFFICIENTS TO MAP CRITICAL PHOSPHOROUS LOADING AREAS1

ABSTRACT: The Export Coefficient model (ECM) is capable of generating reasonable estimates of annual phosphorous loading simply from a watershed's land cover data and export coefficient values (ECVs). In its current form, the ECM assumes that ECVs are homogeneous within each land cover type, yet basic nutrient runoff and hydrological theory suggests that runoff rates have spatial patterns controlled by loading and filtering along the flow paths from the upslope contributing area and downslope dispersal area. Using a geographic information system (GIS) raster, or pixel, modeling format, these contributing area and dispersal area (CADA) controls were derived from the perspective of each individual watershed pixel to weight the otherwise homogeneous ECVs for phosphorous. Although the CADA‐ECM predicts export coefficient spatial variation for a single land use type, the lumped basin load is unaffected by weighting. After CADA weighting, a map of the new ECVs addressed the three fundamental criteria for targeting critical pollutant loading areas: (1) the presence of the pollutant, (2) the likelihood for runoff to carry the pollutant offsite, and (3) the likelihood that buffers will trap nutrients prior to their runoff into the receiving water body. These spatially distributed maps of the most important pollutant management areas were used within New York's West Branch Delaware River watershed to demonstrate how the CADA‐ECM could be applied in targeting phosphorous critical loading areas.     

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.     

Identifying hydrologically sensitive areas: bridging the gap between science and application

Researchers have noted that current water quality protection strategies, like nutrient management plans, lack a sound hydrological underpinning for pollutant transport processes. This is especially true for areas like the northeastern U.S. where copious research has shown that variable source area hydrology largely governs runoff generation. The goal of this study was to develop a scientifically justified method to identify the locations that generate overland flow. Furthermore, this methodology must be computationally simple enough that it can be utilized or incorporated into nutrient management plans and other established water quality tools. We specifically tested the reliability of the 'distance from a stream,'D(s), and the 'topographic index,'lambda, to predict areas with a high propensity for generating overland flow, i.e. hydrologically sensitive areas (HSA). HSAs were defined by their probability of generating runoff, P(sat), based on 30 year simulations using a physically based hydrological model. Using GIS, each location's P(sat) was correlated with D(s) and lambda. We used three Delaware Co., NY watersheds in the New York City watershed system with areas varying in size from 1.6 to 37 km2 and with forested and agricultural land uses. The topographic index gave stronger, more regionally consistent correlations with P(sat) than did D(s). Equations correlating lambda and P(sat) for each month are presented and can be used to estimate hydrological sensitivity in the region surrounding our study watersheds, i.e. in Delaware Co. This work is currently being incorporated into an Internet Mapping System to facilitate user-friendly, on-line identification of HSAs.     

A Water Quality Model for Regional Stream Assessment and Conservation Strategy Development

Non-point-source (NPS) pollution remains the primary source of stream impairment in the United States. Many problems such as eutrophication, sedimentation, and hypoxia are linked with NPS pollution which reduces the water quality for aquatic and terrestrial organisms. Increasingly, NPS pollution models have been used for landscape-scale pollution assessment and conservation strategy development. Our modeling approach functions at a scale between simple landscape-level assessments and complex, data-intensive modeling by providing a rapid, landscape-scale geographic information system (GIS) model with minimal data requirements and widespread applicability. Our model relies on curve numbers, literature-derived pollution concentrations, and land status to evaluate total phosphorus (TP), total nitrogen (TN), and suspended solids (SS) at the reach scale. Model testing in the Chesapeake Bay watershed indicated that predicted distributions of water quality classes were realistic at the reach scale, but precise estimates of pollution concentrations at the local scale can have errors. Application of our model in the tributary watersheds along Lake Ontario suggested that it is useful to managers in watershed planning by rapidly providing important information about NPS pollution conditions in areas where large data gaps exist, comparisons among stream reaches across numerous watersheds are required, or regional assessments are sought.