RUNOFF MODELING OF A MOUNTAINOUS CATCHMENT USING TOPMODEL: A CASE STUDY1

ABSTRACT: The rainfall‐runoff response of the Tygarts Creek Catchment in eastern Kentucky is studied using TOPMODEL, a hydrologic model that simulates runoff at the catchment outlet based on the concepts of saturation excess overland flow and subsurface flow. Unlike the traditional application of this model to continuous rainfall‐runoff data, the use of TOPMOEL in single event runoff modeling, specifically floods, is explored here. TOPMODEL utilizes a topographic index as an indicator of the likely spatial distribution of rainfall excess generation in the catchment. The topographic index values within the catchment are determined using the digital terrain analysis procedures in conjunction with digital elevation model (DEM) data. Select parameters in TOPMODEL are calibrated using an iterative procedure to obtain the best‐fit runoff hydrograph. The calibrated parameters are the surface transmissivity, T_O, the transmissivity decay parameter, m, and the initial moisture deficit in the root zone, S_r0. These parameters are calibrated using three storm events and verified using three additional storm events. Overall, the calibration results obtained in this study are in general agreement with the results documented from previous studies using TOPMODEL. However, the parameter values did not perform well during the verification phase of this study.     

Variable Source Area Hydrology Modeling with the Water Erosion Prediction Project Model

In nondegraded watersheds of humid climates, subsurface flow patterns determine where the soil saturates and where surface runoff is occurring. Most models necessarily use infiltration‐excess (i.e., Hortonian) runoff for predicting runoff and associated constituents because subsurface flow algorithms are not included in the model. In this article, we modify the Water Erosion Prediction Project (WEPP) model to simulate subsurface flow correctly and to predict the spatial and temporal location of saturation, the associated lateral flow and surface runoff, and the location where the water can re‐infiltrate. The modified model, called WEPP‐UI, correctly simulated the hillslope drainage data from the Coweeta Hydrologic Laboratory hillslope plot. We applied WEPP‐UI to convex, concave, and S‐shaped hillslope profiles, and found that multiple overland flow elements are needed to simulate distributed lateral flow and runoff well. Concave slopes had the greatest runoff, while convex slopes had the least. Our findings concur with observations in watersheds with saturation‐excess overland flow that most surface runoff is generated on lower concave slopes, whereas on convex slopes runoff infiltrates before reaching the stream. Since the WEPP model is capable of simulating both saturation‐excess and infiltration‐excess runoff, we expect that this model will be a powerful tool in the future for managing water quality.     

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.     

Incorporating Variable Source Area Hydrology into a Spatially Distributed Direct Runoff Model1

Buchanan, Brian, Zachary M. Easton, Rebecca Schneider, and M. Todd Walter, 2011. Incorporating Variable Source Area Hydrology Into a Spatially Distributed Direct Runoff Model. Journal of the American Water Resources Association (JAWRA) 48(1): 43‐60. DOI: 10.1111/j.1752‐1688.2011.00594.x Abstract: Few hydrologic models simulate both variable source area (VSA) hydrology, and runoff‐routing at high enough spatial resolutions to capture fine‐scale hydrologic pathways connecting VSA to the stream network. This paper describes a geographic information system‐based operational model that simulates the spatio‐temporal dynamics of VSA runoff generation and distributed runoff‐routing, including through complex artificial drainage networks. The model combines the Natural Resource Conservation Service’s Curve Number (CN) equation for estimating storm runoff with the topographic index concept for predicting the locations of VSA and a runoff‐routing algorithm into a new spatially distributed direct hydrograph (SDDH) model (SDDH‐VSA). Using a small agricultural watershed in central New York, SDDH‐VSA results were compared to those from a SDDH model using the traditional land use assumptions for the CN (SDDH‐CN). The SDDH‐VSA model generally agreed better with observed discharge than the SDDH‐CN model (average, Nash‐Sutcliffe efficiency of 0.69 vs. 0.58, respectively) and resulted in more realistic spatial patterns of runoff‐generating areas. The SDDH approach did not correctly capture the timing of runoff from small storms in dry periods. Despite this type of limitation, SDDH‐VSA extends the applicability of the SDDH technique to VSA conditions, providing a basis for new tools to help identify critical management areas and assess water quality risks due to landscape alterations.     

THE IMPACT OF RUNOFF GENERATION MECHANISMS ON THE LOCATION OF CRITICAL SOURCE AREAS1

ABSTRACT: Identifying phosphorus (P) source areas and transport pathways is a key step in decreasing P loading to natural water systems. This study compared the effects of two modeled runoff generation processes ‐ saturation excess and infiltration excess ‐ on total phosphorus (TP) and soluble reactive phosphorus (SRP) concentrations in 10 catchment streams of a Catskill mountain watershed in southeastern New York. The spatial distribution of runoff from forested land and agricultural land was generated for both runoff processes; results of both distributions were consistent with Soil Conservation Service‐Curve Number (SCS‐CN) theory. These spatial runoff distributions were then used to simulate stream concentrations of TP and SRP through a simple equation derived from an observed relation between P concentration and land use; empirical results indicate that TP and SRP concentrations increased with increasing percentage of agricultural land. Simulated TP and SRP stream concentrations predicted for the 10 catchments were strongly affected by the assumed runoff mechanism. The modeled TP and SRP concentrations produced by saturation excess distribution averaged 31 percent higher and 42 percent higher, respectively, than those produced by the infiltration excess distribution. Misrepresenting the primary runoff mechanism could not only produce erroneous concentrations, it could fail to correctly locate critical source areas for implementation of best management practices. Thus, identification of the primary runoff mechanism is critical in selection of appropriate models in the mitigation of nonpoint source pollution. Correct representation of runoff processes is also critical in the future development of biogeochemical transport models, especially those that address nutrient fluxes.     

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.     

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.     

EFFECTS OF INITIAL ABSTRACTION AND URBANIZATION ON ESTIMATED RUNOFF USING CN TECHNOLOGY1

ABSTRACT: Few studies have been conducted to explore the effects of initial abstraction on estimated direct runoff despite the widespread use of the curve number (CN) method in many hydrologic models to estimate direct runoff. In this study, use of a 5 percent ratio of initial abstraction (I_a) to storage (S) to estimate daily direct runoff with modified CN values for a 5 percent I_a/S value was investigated using the Long‐Term Hydrologic Impact Assessment (L‐THIA) geographic information system (GIS). In addition, the effects on estimated runoff of altering the hydrologic soil group due to urbanization were investigated. The L‐THIA model was applied to the Indiana Little Eagle Creek watershed with 5 percent and 20 percent I_a/S values, considering hydrologic soil group alteration due to urbanization. The results indicate that uses of a 5 percent l_a/S and modified CN values and Hydrologic Soil Group D for urbanized areas in model runs can improve long term direct runoff prediction.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

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.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

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.     

Examining Modeling Approaches for the Rainfall‐Runoff Process in Wildfire‐Affected Watersheds: Using San Dimas Experimental Forest

Wildfire can significantly change watershed hydrological processes resulting in increased risks for flooding, erosion, and debris flow. The goal of this study was to evaluate the predictive capability of hydrological models in estimating post‐fire runoff using data from the San Dimas Experimental Forest (SDEF), San Dimas, California. Four methods were chosen representing different types of post‐fire runoff prediction methods, including a Rule of Thumb, Modified Rational Method (MODRAT), HEC‐HMS Curve Number, and KINematic Runoff and EROSion Model 2 (KINEROS2). Results showed that simple, empirical peak flow models performed acceptably if calibrated correctly. However, these models do not reflect hydrological mechanisms and may not be applicable for predictions outside the area where they were calibrated. For pre‐fire conditions, the Curve Number approach implemented in HEC‐HMS provided more accurate results than KINEROS2, whereas for post‐fire conditions, the opposite was observed. Such a trend may imply fundamental changes from pre‐ to post‐fire hydrology. Analysis suggests that the runoff generation mechanism in the watershed may have temporarily changed due to fire effects from saturation‐excess runoff or subsurface storm dominated complex mechanisms to an infiltration‐excess dominated mechanism. Infiltration modeling using the Hydrus‐1D model supports this inference. Results of this study indicate that physically‐based approaches may better reflect this trend and have the potential to provide consistent and satisfactory prediction.     

Hydrologic Landscape Classification to Estimate Bristol Bay,Alaska Watershed Hydrology

Hydrologic landscapes (HLs) have proven to be a useful tool for broad scale assessment and classification of landscapes across the United States as they help organize larger geographical areas into areas of similar hydrologic characteristics. We developed a HL classification for the Bristol Bay watershed of southwest Alaska that incorporates indices of annual climate and seasonality, terrain, geology, and the influences of large lakes and glaciers. A HL classification is particularly useful in this large watershed because of its hydrologic and landscape variability, important salmon fishery, variety of environmental and potential anthropogenic stressors, and lack of widespread hydrologic data. Following creation of Bristol Bay basin‐wide HL classes, we compared the HL distributions within watersheds grouped by two calculated runoff parameters derived from available long‐term streamflow records and found HL distributions within these groups provided predictive insight on hydrologic behavior. Using these developed runoff groups, we estimated expected hydrologic behavior in watersheds across the larger Bristol Bay watershed that lacked gauged streamflow records. The HL approach provides a scientific basis for estimating the first‐order hydrologic behavior of watersheds and landscapes that lack detailed hydrologic information.     

Comparison of Indexes for Prioritizing Placement of Water Quality Buffers in Agricultural Watersheds1

Dosskey, Michael G. and Zeyuan Qiu, 2011. Comparison of Indexes for Prioritizing Placement of Water Quality Buffers in Agricultural Watersheds. Journal of the American Water Resources Association (JAWRA) 47(4):662‐671. DOI: 10.1111/j.1752‐1688.2011.00532.x Abstract: Five physically based, spatially distributed, empirical indexes were compared for the degree to which they identified the same or different locations in watersheds where vegetative buffers would function better for reducing agricultural nonpoint source pollution. All five indexes were calculated on a 10 m × 10 m digital elevation grid on agricultural land in the 144‐km² Neshanic River watershed in New Jersey. The indexes included the topography‐based Wetness Index (WI) and Topographic Index (TI) and three soil survey‐based indexes (sediment trapping efficiency [STE], water trapping efficiency [WTE], and groundwater interaction [GI]). Results showed that each index associated higher pollution risk and mitigation potential to a different part of the landscape. The WI and TI identified swales and riparian areas where runoff converges, whereas STE and WTE identified upland sites. The STE and WTE lack the fine scale of slope resolution and the accounting for convergent runoff patterns that can be important for properly locating buffers in some watersheds. The GI index indicates the existence of a shallow water table but the correspondence with WI‐ and TI‐identified sites was only modest. For watersheds where pollutant loading is generated by both saturation‐excess (emphasized by TI and WI) and infiltration‐excess processes (emphasized by STE and WTE), the indexes could be complementary. However, techniques would be needed for properly apportioning priority among sites identified by each index.     

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.     

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.     

A Geospatial Methodology to Identify Locations of Concentrated Runoff from Agricultural Fields

A geospatial methodology has been developed that utilizes high resolution lidar‐derived DEMs to help track runoff from agricultural fields and identify areas of potential concentrated flow through vegetated riparian areas in the Coastal Plain of Virginia. Points of concentrated flow are identified across 74 agricultural fields within the Virginia portion of the Chesapeake Bay watershed. On average, 70% of the surface area of the agricultural fields analyzed drains through less than 20 m of the field margin, and on average 81% of the field surface area drains through 1% or less of the field margin. Within the riparian buffer, locations that were predicted by the geospatial model to have high levels of concentrated flow were found to exhibit evidence of channelization. Results indicate that flow concentration and channelized flow through vegetated riparian areas may be common along the margin of agricultural fields, resulting in vegetated riparian areas that are less effective at sediment trapping than assumed. Additional results suggest that the regulations governing the location and width of vegetated riparian may not be sufficient to achieve goals for reducing sediment and nutrient runoff from nonpoint agricultural sources. Combined with the increasing availability of lidar‐derived DEMs, the geospatial model presented has the potential to advance management practices aimed at reducing nonpoint source pollution leaving agricultural fields.     

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.     

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.