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.     

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.     

SENSITIVITY CONSIDERATIONS WHEN MODELING HYDROLOGIC PROCESSES WITH DIGITAL ELEVATION MODEL1

ABSTRACT: The purpose of this paper is to investigate the sensitivity of a hydrologic models to the type of DEM used. This was done while modeling basin water quality with 1:24,000 and 1:250,000 U.S. Geological Survey DEMs as input to model hydro‐logical processes. The manner in which the model results were sensitive to the choice of raster cell size (scale) is investigated in this study. The Broadhead watershed, located in New Jersey, USA, was chosen as a study area. Curve numbers were estimated by a trial and error to match simulated and observed total discharge. Monthly runoff for the watershed was used in the calibration process. Higher runoff volumes were simulated by the model when the 1:24,000 DEM were used as input data, probably due to the finer resolution which simulated increased average slope and hence higher estimated runoff from the watershed. As the simulated slope of the watershed is flatten with the 1:250,000 DEM, the response of stream flow was delayed and simulated less runoff volume.     

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

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

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.     

Development of a SWAT Patch for Better Estimation of Sediment Yield in Steep Sloping Watersheds1

Abstract: The watershed scale Soil and Water Assessment Tool (SWAT) model divides watersheds into smaller subwatersheds for simulation of rainfall‐runoff and sediment loading at the field level and routing through stream networks. Typically, the SWAT model first needs to be calibrated and validated for accurate estimation through adjustment of sensitive input parameters (i.e., Curve Number values, USLE P, slope and slope‐length, and so on). However, in some instances, SWAT‐simulated results are greatly affected by the watershed delineation and Digital Elevation Models (DEM) cell size. In this study, the SWAT ArcView GIS Patch II was developed for steep sloping watersheds, and its performance was evaluated for various threshold values and DEM cell size scenarios when delineating subwatersheds using the SWAT model. The SWAT ArcView GIS Patch II was developed using the ArcView GIS Avenue program and Spatial Analyst libraries. The SWAT ArcView GIS Patch II improves upon the SWAT ArcView GIS Patch I because it reflects the topographic factor in calculating the field slope‐length of Hydrologic Response Units in the SWAT model. The simulated sediment value for 321 subwatersheds (watershed delineation threshold value of 25 ha) is greater than that for 43 subwatersheds (watershed delineation threshold value of 200 ha) by 201% without applying the SWAT ArcView GIS Patch II. However, when the SWAT ArcView GIS Patch II was applied, the difference in simulated sediment yield decreases for the same scenario (i.e., difference in simulated sediment with 321 subwatersheds and 43 subwatersheds) was 12%. The simulated sediment value for DEM cell size of 50 m is greater than that for DEM cell size of 10 m by 19.8% without the SWAT ArcView GIS Patch II. However, the difference becomes smaller (3.4% difference) between 50 and 10 m with the SWAT ArcView GIS Patch II for the DEM scenarios. As shown in this study, the SWAT ArcView GIS Patch II can reduce differences in simulated sediment values for various watershed delineation and DEM cell size scenarios. Without the SWAT ArcView GIS Patch II, variations in the SWAT‐simulated results using various watershed delineation and DEM cell size scenarios could be greater than those from input parameter calibration. Thus, the results obtained in this study show that the SWAT ArcView GIS Patch II should be used when simulating hydrology and sediment yield for steep sloping watersheds (especially if average slope of the subwatershed is >25%) for more accurate simulation of hydrology and sediment using the SWAT model. The SWAT ArcView GIS Patch II is available at http://www.EnvSys.co.kr/~swat for free download.     

STREAM DISCHARGE PREDICTION VIA A GRID BASED SOIL WATER ROUTING WITH PADDY FIELDS1

ABSTRACT: A grid based daily hydrologic model for a watershed with paddy fields was developed to predict the stream discharge. ASCII formatted elevation, soil, and land use data supported by the GRASS Geographic Information System are used to generate distributed results such as surface runoff and subsurface flow, soil water content, and evapotranspiration. The model uses a single flow path algorithm and simulates a water balance at each grid element. A linear reservoir assumption was used to predict subsurface runoff components. The model was applied to a 75.6 km² watershed located in the middle of South Korea, and observed stream flow hydrographs from 1995 and 1996 were compared to model predictions. The stream flow predictions of 1995 and 1996 generally agreed with the observed flow, resulting in a Nash‐Sutcliffe efficiency R² of 0.60 and 0.62, respectively. The hydraulic conductivity for percolating water through the saturated layer affected baseflow generation. The levee height of the paddy influenced the time and magnitude of the surface runoff, depending on irrigation management. The model will be used for making low flow management decisions by evaluating the role of each land use to stream flow, especially in case of paddy decrease by gradual urbanization of a watershed.     

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.     

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.     

Assessing Surface Water Quality and Its Relation with Urban Land Cover Changes in the Lake Calumet Area, Greater Chicago

Urban land use and land cover change significantly affect spatial and temporal patterns of runoff, which in turn impacts surface water quality. With the exponential growth in urban areas over the past three decades, changes in land use and land cover to cater for the growth of cities has been a conspicuous spectacle in urban spaces. The main goal of this study was to assess the impacts of land cover change on runoff and surface water quality using a partial area hydrology framework. The study employed ArcHydro GIS extension and a modified version of Long-Term Hydrologic and Nonpoint Source Pollution model (L-THIA-NPS) in estimating runoff and nonpoint source pollutant concentration around Lake Calumet between 1992 and 2001. Data employed include National Land Cover Data set, rainfall data, digital elevation model (DEM), Soil Survey Geographic (SSURGO) data, and The United States Environmental Protection Agency’s STORET (storage and retrieval) water quality data. The model was able to predict surface water quality reasonably well over the study period. Sensitivity analysis facilitated a manual calibration of the model. Model validation was executed by comparing simulated results following calibration and observed water quality data for the study area. The study demonstrates that the level of concentration of nonpoint source pollutants in surface water within an urban watershed heavily depends on the spatiotemporal variations in areas that contribute towards runoff compared to the spatial extent of change in major land use/land cover.     

DEVELOPMENT OF WATERSIIED MODELS FOR TWO SIERRA NEVADA BASINS USING A GEOGRAPHIC INFORMATION SYSTEM1

ABSTRACT: Techniques were developed using vector and raster data in a geographic information system (GIS) to define the spatial variability of watershed characteristics in the north-central Sierra Nevada of California and Nevada and to assist in computing model input parameters. The U.S. Geological Survey's Precipitation-Runoff Modeling System, a physically based, distributed-parameter watershed model, simulates runoff for a basin by partitioning a watershed into areas that each have a homogeneous hydrologic response to precipitation or snowmelt. These land units, known as hydrologic-response units (HRU's), are characterized according to physical properties, such as altitude, slope, aspect, land cover, soils, and geology, and climate patterns. Digital data were used to develop a GIS data base and HRIJ classification for the American River and Carson River basins. The following criteria are used in delineating HRU's: (1) Data layers are hydrologically significant and have a resolution appropriate to the watershed's natural spatial variability, (2) the technique for delineating HRU's accommodates different classification criteria and is reproducible, and (3) HRU's are not limited by hydrographic-subbasin boundaries. HRU's so defined are spatially noncontiguous. The result is an objective, efficient methodology for characterizing a watershed and for delineating HRU's. Also, digital data can be analyzed and transformed to assist in defining parameters and in calibrating the model.     

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.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

Classification and spatial mapping of riparian habitat with applications toward management of streams impacted by nonpoint source pollution

Management of riparian habitats has been recognized for its importance in reducing instream effects of agricultural nonpoint source pollution. By serving as a buffer, well structured riparian habitats can reduce nonpoint source impacts by filtering surface runoff from field to stream. A system has been developed where key characteristics of riparian habitat, vegetation type, height, width, riparian and shoreline bank slope, and land use are classified as discrete categorical units. This classification system recognizes seven riparian vegetation types, which are determined by dominant plant type. Riparian and shoreline bank slope, in addition to riparian width and height, each consist of five categories. Classification by discrete units allows for ready digitizing of information for production of spatial maps using a geographic information system (GIS). The classification system was tested for field efficiency on Tom Beall Creek watershed, an agriculturally impacted third-order stream in the Clearwater River drainage, Nez Perce County, Idaho, USA. The classification system was simple to use during field applications and provided a good inventory of riparian habitat. After successful field tests, spatial maps were produced for each component using the Professional Map Analysis Package (pMAP), a GIS program. With pMAP, a map describing general riparian habitat condition was produced by combining the maps of components of riparian habitat, and the condition map was integrated with a map of soil erosion potential in order to determine areas along the stream that are susceptible to nonpoint source pollution inputs. Integration of spatial maps of riparian classification and watershed characteristics has great potential as a tool for aiding in making management decisions for mitigating off-site impacts of agricultural nonpoint source pollution.     

APPLICATIONS OF A GIS FOR MODELING THE SENSITIVITY OF WATER RESOURCES TO ALTERATIONS IN CLIMATE IN THE GUNNISON RIVER BASIN,COLORADO1

ABSTRACT: The Gunnison River drains a mountainous basin in western Colorado, and is a large contributor of water to the Colorado River. As part of a study to assess water resource sensitivity to alterations in climate in the Gunnison River basin, climatic and hydrologic processes are being modeled. A geographic information system (GIS) is being used in this study as a link between data and modelers - serving as a common data base for project personnel with differing specialties, providing a means to investigate the effects of scale on model results, and providing a framework for the transfer of parameter values among models. Specific applications presented include: (1) developing elevation grids for a precipitation model from digital elevation model (DEM) point-elevation values, and visualizing the effects of grid resolution on model results; (2) using a GIS to facilitate the definition and parameterization of a distributed-parameters, watershed model in multiple basins; and (3) nesting atmospheric and hydrologic models to produce possible scenarios of climate change.     

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.     

MODELING FRAMEWORK FOR MANAGING COPPER RUNOFF IN URBAN WATERSHEDS1

ABSTRACT: A modeling framework was developed for managing copper runoff in urban watersheds that incorporates water quality characterization, watershed land use areas, hydrologic data, a statistical simulator, a biotic ligand binding model to characterize acute toxicity, and a statistical method for setting a watershed specific copper loading. The modeling framework is driven by export coefficients derived from water quality parameters and hydrologic inputs measured in an urban watershed's storm water system. This framework was applied to a watershed containing a copper roof built in 1992. A series of simulations was run to predict the change in receiving stream water chemistry caused by roof aging and to determine the maximum copper loading (at the 99 percent confidence level) a watershed could accept without causing acute toxicity in the receiving stream. Forecasting the amount of copper flux responsible for exceeding the assimilation capacity of a watershed can be directly related to maximum copper loadings responsible for causing toxicity in the receiving streams. The framework developed in this study can be used to evaluate copper utilization in urban watersheds.     

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.     

SPATIAL DISTRIBUTION OF POINT SOIL MOISTURE ESTIMATES USING LANDSAT TM DATA AND FUZZY-C CLASSIFICATION1

ABSTRACT: Many hydrologic models have input data requirements that are difficult to satisfy for all but a few well-instrumented, experimental watersheds. In this study, point soil moisture in a mountain watershed with various types of vegetative cover was modeled using a generalized regression model. Information on sur-ficial characteristics of the watershed was obtained by applying fuzzy set theory to a database consisting of only satellite and a digital elevation model (DEM). The fuzzy-c algorithm separated the watershed into distinguishable classes and provided regression coefficients for each ground pixel. The regression model used the coefficients to estimate distributed soil moisture over the entire watershed. A soil moisture accounting model was used to resolve temporal differences between measurements at prototypical measurement sites and validation sites. The results were reasonably accurate for all classes in the watershed. The spatial distribution of soil moisture estimates corresponded accurately with soil moisture measurements at validation sites on the watershed. It was concluded that use of the regression model to distribute soil moisture from a specified number of points can be combined with satellite and DEM information to provide a reasonable estimation of the spatial distribution of soil moisture for a watershed.     

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.