MONTHLY RUNOFF PREDICTIONS BASED ON RAINFALL FORECASTS IN A SMALL OKLAHOMA WATERSHED1

ABSTRACT: Conditions under which monthly rainfall forecasts translate into monthly runoff predictions that could support water resources planning and management activities were investigated on a small watershed in central Oklahoma. Runoff response to rainfall forecasts was simulated using the hydrologic model SWAT. Eighteen scenarios were examined that represented combinations of wet, average, and dry antecedent rainfall conditions, with wet, normal, and dry forecasted rainfall. Results suggest that for the climatic and physiographic conditions under consideration, rainfall forecasts could offer potential application opportunities in surface water resources but only under certain conditions. Pronounced wet and dry antecedent rainfall conditions were shown to have greater impact on runoff than forecasts, particularly in the first month of a forecast period. Large forecast impacts on runoff occurred under wet antecedent conditions, when the fraction of forecasted rainfall contributing to runoff was greatest. Under dry antecedent conditions, most of the forecasted rainfall was absorbed in the soil profile, with little immediate runoff response. Persistent three‐month forecasts produced stronger impacts on runoff than one‐month forecasts due to cumulative effects in the hydrologic system. Runoff response to antecedent conditions and forecasts suggest a highly asymmetric utility function for rainfall forecasts, with greatest decision‐support potential for persistent wet forecasts under wet antecedent conditions when the forecast signal is least dampened by soil‐storage effects. Under average and dry antecedent conditions, rainfall forecasts showed little potential value for practical applications in surface water resources assessments.     

IMPROVED RAINFALL/RUNOFF ESTIMATES USING REMOTELY SENSED SOIL MOISTURE1

ABSTRACT: Remotely sensed soil moisture data measured during the Southern Great Plains 1997 (SGP97) experiment in Oklahoma were used to characterize antecedent soil moisture conditions for the Soil Conservation Service (SCS) curve number method. The precipitation‐adjusted curve number and the soil moisture were strongly related (r²= 0.70). Remotely sensed soil moisture fields were used to adjust the curve numbers and the runoff estimates for five watersheds, in the Little Washita watershed; the results ranged from 2.8 km² to 601.6 km². The soil moisture data were applied at two spatial scales, a finer one (800 m) measuring spatial resolution and a coarser one (28 km). The root mean square error (RMSE) and the mean absolute error (MAE) of the runoff estimated by the standard SCS method was reduced by nearly 50 percent when the 800 m soil moisture data were used to adjust the curve number. The coarser scale soil moisture data also significantly reduced the error in the runoff predictions with 41 percent and 28 percent reductions in MAE and RMSE, respectively. The results suggest that remote sensing of soil moisture, when combined with the SCS method, can improve rainfall runoff predictions at a range of spatial scales.     

Sensitivity of SCS Models to Curve Number Variation1

ABSTRACT: The Soil Conservation Service (SCS) models, including the TR-20 computer program and the simplified methods in TR-55, are widely used in hydrologic design. The runoff curve number (CN), which is an important input parameter to SCS models, is defined in terms of land use tretments, hydrologic, condition, antecedent soil moisture, and soil type. The objective of this study was to evaluate the sensitivity of the SCS models to errors in CN estimates. The results show that the effects of CN variation decrease as the design rainfall depth increases, such as for the larger storm events. The value and use of the sensitivity curves are demonstrated using a comparison of Landsat and conventionally derived curve numbers for three watersheds in Pennsylvania.     

Investigation of the Curve Number Method For Surface Runoff Estimation In Tropical Regions

This study tests the applicability of the curve number (CN) method within the Soil and Water Assessment Tool (SWAT) to estimate surface runoff at the watershed scale in tropical regions. To do this, surface runoff simulated using the CN method was compared with observed runoff in numerous rainfall‐runoff events in three small tropical watersheds located in the Upper Blue Nile basin, Ethiopia. The CN method generally performed well in simulating surface runoff in the studied watersheds (Nash‐Sutcliff efficiency [NSE] > 0.7; percent bias [PBIAS] < 32%). Moreover, there was no difference in the performance of the CN method in simulating surface runoff under low and high antecedent rainfall (PBIAS for both antecedent conditions: ~30%; modified NSE: ~0.4). It was also found that the method accurately estimated surface runoff at high rainfall intensity (e.g., PBIAS < 15%); however, at low rainfall intensity, the CN method repeatedly underestimated surface runoff (e.g., PBIAS > 60%). This was possibly due to low infiltrability and valley bottom saturated areas typical of many tropical soils, indicating that there is scope for further improvements in the parameterization/representation of tropical soils in the CN method for runoff estimation, to capture low rainfall‐intensity events. In this study the retention parameter was linked to the soil moisture content, which seems to be an appropriate approach to account for antecedent wetness conditions in the tropics.     

Decreased Runoff Response to Precipitation,Little Missouri River Basin,Northern Great Plains,USA

High variability in precipitation and streamflow in the semiarid northern Great Plains causes large uncertainty in water availability. This uncertainty is compounded by potential effects of future climate change. We examined historical variability in annual and growing season precipitation, temperature, and streamflow within the Little Missouri River Basin and identified differences in the runoff response to precipitation for the period 1976‐2012 compared to 1939‐1975 (n = 37 years in both cases). Computed mean values for the second half of the record showed little change (<5%) in annual or growing season precipitation, but average annual runoff at the basin outlet decreased by 22%, with 66% of the reduction in flow occurring during the growing season. Our results show a statistically significant (p < 0.10) 27% decrease in the annual runoff response to precipitation (runoff ratio). Surface‐water withdrawals for various uses appear to account for <12% of the reduction in average annual flow volume, and we found no published or reported evidence of substantial flow reduction caused by groundwater pumping in this basin. Results of our analysis suggest that increases in monthly average maximum and minimum temperatures, including >1°C increases in January through March, are the dominant driver of the observed decrease in runoff response to precipitation in the Little Missouri River Basin.     

HYDROLOGY OF A SMALL IMPOUNDMENT AND WATERSHED IN CENTRAL OKLAHOMA1

ABSTRACT: West Bitter Creek floodwater retarding structure site 3 in South Central Oklahoma was instrumented and records obtained and analyzed to obtain information concerning an impoundment water budget that is useful to landowners and designers of these impoundments. On-site loss of water from the impoundment was only 17 percent of the inflow during three years when the annual precipitation averaged 26 inches and the annual inflow averaged 1.4 inches. Runoff from an eroded area with no farm ponds was about 70 percent greater per unit area than from a portion of the watershed where 71 percent of the drainage area was controlled by farm ponds. A previous study indicated, however, that the ponds were reducing runoff only 13 percent. Loss of top soil increases runoff considerably. Only 24 percent of the total runoff into the impoundment was base flow. The flow rate into the impoundment was less than 0.05 cfs 70 percent of the time, and the inflow rate exceeded 10 cfs only 1 percent of the time. SCS runoff curve numbers varied between 57 and 96 for the impoundment watershed with an inverse relation between precipitation amount and curve number apprently caused by partial area runoff from impervious and semi-impervious areas. A comparison of measured event runoff versus event runoff computed by the SCS curve numbers gave an r² of only 0.44. However, the total computed surface runoff for eight years of record was less than 1 percent below the measured runoff which indicated the curve number method was a good tool for predicting long term runoff for the 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.     

Curve Number Derivation for Watersheds Draining Two Headwater Streams in Lower Coastal Plain South Carolina,USA

The objective of this study was to assess curve number (CN) values derived for two forested headwater catchments in the Lower Coastal Plain (LCP) of South Carolina using a three‐year period of storm event rainfall and runoff data in comparison with results obtained from CN method calculations. Derived CNs from rainfall/runoff pairs ranged from 46 to 90 for the Upper Debidue Creek (UDC) watershed and from 42 to 89 for the Watershed 80 (WS80). However, runoff generation from storm events was strongly related to water table elevation, where seasonally variable evapotranspirative wet and dry moisture conditions persist. Seasonal water table fluctuation is independent of, but can be compounded by, wet conditions that occur as a result of prior storm events, further complicating flow prediction. Runoff predictions for LCP first‐order watersheds do not compare closely to measured flow under the average moisture condition normally associated with the CN method. In this study, however, results show improvement in flow predictions using CNs adjusted for antecedent runoff conditions and based on water table position. These results indicate that adaptations of CN model parameters are required for reliable flow predictions for these LCP catchments with shallow water tables. Low gradient topography and shallow water table characteristics of LCP watersheds allow for unique hydrologic conditions that must be assessed and managed differently than higher gradient watersheds.     

EVAPOTRANSPIRATION AND SOIL MOISTURE DYNAMICS ON A SEMIARID PONDEROSA PINE HILLSLOPE1

ABSTRACT: Detailed measurements of soil moisture and ET in semiarid forest environments have not been widely reported in the literature. In this study, soil moisture and water balance components were measured over a four‐year period on a semiarid ponderosa pine hillslope, with evapotranspiration (ET) determined as the residual of measured precipitation, runoff, and change in soil moisture storage. ET accounts for approximately 95 percent of the water budget and has a distinctly bimodal annual pattern, with peaks occurring after spring snowmelt and during the late summer monsoon season, periods that coincide with high soil moisture. Weekly growing season ET rates determined by the hillslope water balance are found to be invariably below calculated potential rates. Normalized ET rates are linearly correlated (r²= 0.62) with soil moisture; therefore, a simple linear relation is proposed. Growing season soil moisture dynamics were modeled based on this relation. Results are in fair agreement (r²= 0.63) with the observed soil moisture data over the four growing seasons; however, for two dry summers with little surface runoff, much better results (r² > 0.90) were obtained.     

Calibration and Verification of SWMM for Low Impact Development

The Storm Water Management Model was used to simulate runoff and nutrient export from a low impact development (LID) watershed and a watershed using traditional runoff controls. Predictions were compared to observed values. Uncalibrated simulations underpredicted weekly runoff volume and average peak flow rates from the multiple subcatchment LID watershed by over 80%; the single subcatchment traditional watershed had better predictions. Saturated hydraulic conductivity, Manning's n for swales, and initial soil moisture deficit were sensitive parameters. After calibration, prediction of total weekly runoff volume for the LID and traditional watersheds improved to within 12 and 5% of observed values, respectively. For the validation period, predicted total weekly runoff volumes for the LID and traditional watersheds were within 6 and 2% of observed values, respectively. Water quality simulation was less successful, Nash–Sutcliffe coefficients >0.5 for both calibration and validation periods were only achieved for prediction of total nitrogen export from the LID watershed. Simulation of a 100‐year, 24‐h storm resulted in a runoff coefficient of 0.46 for the LID watershed and 0.59 for the traditional watershed. Results suggest either calibration is needed to improve predictions for LID watersheds or expanded look‐up tables for Green–Ampt infiltration parameter values that account for compaction of urban soil and antecedent conditions are needed.     

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.     

SAGEBRUSH RANGELAND HYDROLOGY AND EVALUATION OF THE SPUR HYDROLOGY MODEL1

ABSTRACT: An excellent hydrologic record on sagebrush range-land has been developed at the Reynolds Creek Experimental Watershed in southwestern Idaho. The objectives of this paper were two-fold: (1) to analyze and describe the hydrologic record (8–18 years) from four sagebrush watersheds (1–83 ha); and (2) to evaluate the hydrology component of SPUR, a comprehensive rangeland model. The watersheds represent a gradient in elevation (1180–1658 m) and precipitation (240–350 mm/yr). Runoff was a small fraction (> 2 percent) of the total water budget for all of the watersheds. It occurred very infrequently at the three lower elevation watersheds: Summit, Flats, and Nancy Gulch. At Lower Sheep, the highest elevation watershed, runoff occurred most years for a period of 1 to 17 weeks in the winter. Frozen soil combined with rainfall or snowmelt was associated with most of the runoff from Flats and Nancy Gulch. At Summit summertime thunderstorms produced all of the runoff. The average annual sediment yield from all of the watersheds was low (17–950 kg/ha). It was highest from Summit, which had well developed alluvial channels and very steep slopes. SPUR was able to simulate runoff with reasonable accuracy only at Summit, where frozen soils were not a factor. There was poor correlation between predicted and actual annual 8ediment loss. The model tended to overpredict evapotranspiration early in the growing season and underpredict it in the late summer.     

Effects of Watershed Impervious Cover on Dissolved Silica Loading in Storm Flow1

Abstract: Dissolved silica (DSi) availability is a factor that affects the composition of algal populations in aquatic ecosystems. DSi cycling is tightly linked to the hydrological cycle, which is affected by human alterations of the landscape. Development activities that increase impervious cover change watershed hydrology and may increase the discharge of DSi‐poor rainwater and decrease the discharge of DSi‐rich ground water into aquatic ecosystems, possibly shifting algal community composition toward less desirable assemblages. In this study, DSi loadings from two adjacent coastal watersheds with different percent impervious cover were compared during four rain and five nonrain events. Loadings in the more impervious watershed contained a significantly larger proportion of surface runoff than base flow (ground‐water discharge) and had lower [DSi] water during rain events than the less impervious watershed. Application of the Soil Conservation Service Curve Number (CN) method showed that the minimum rainfall height necessary to yield runoff was significantly lower for the more impervious watershed, implying that runoff volumes increase with impervious cover as well as the frequency of runoff‐yielding events. Empirical data collected during this study and estimates derived from the CN method suggest that impervious cover may be responsible for both short‐term DSi limitation during rain events as well as long‐term reduction of DSi inputs into aquatic ecosystems.     

Modeling Postfire Response and Recovery using the Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS)1

Abstract: This paper investigates application of the Army Corps of Engineers’ Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS) to a burned watershed in San Bernardino County, California. We evaluate the HEC‐HMS’ ability to simulate discharge in prefire and postfire conditions in a semi arid watershed and the necessary parameterizations for modeling hydrologic response during the immediate, and subsequent recovery, period after a wildfire. The model is applied to City Creek watershed, which was 90% burned during the Old Fire of October 2003. An optimal spatial resolution for the HEC‐HMS model was chosen based on an initial sensitivity analysis of subbasin configurations and related model performance. Five prefire storms were calibrated for the selected model resolution, defining a set of parameters that reasonably simulate prefire conditions. Six postfire storms, two from each of the following rainy (winter) seasons were then selected to simulate postfire response and evaluate relative changes in parameter values and model behavior. There were clear trends in the postfire parameters [initial abstractions (I_a), curve number (CN), and lag time] that reveal significant (and expected) changes in watershed behavior. CN returns to prefire (baseline) values by the end of Year 2, while I_a approaches baseline by the end of the third rainy season. However, lag time remains significantly lower than prefire values throughout the three‐year study period. Our results indicate that recovery of soil conditions and related runoff response is not entirely evidenced by the end of the study period (three rainy seasons postfire). Understanding the evolution of the land surface and related hydrologic properties during the highly dynamic postfire period, and accounting for these changes in model parameterizations, will allow for more accurate and reliable discharge simulations in both the immediate, and subsequent, rainy seasons following fire.     

Effects of calibration on L-THIA GIS runoff and pollutant estimation

Urbanization can result in alteration of a watershed's hydrologic response and water quality. To simulate hydrologic and water quality impacts of land use changes, the Long-Term Hydrologic Impact Assessment (L-THIA) system has been used. The L-THIA system estimates pollutant loading based on direct runoff quantity and land use based pollutant coefficients. The accurate estimation of direct runoff is important in assessing water quality impacts of land use changes. An automated program was developed to calibrate the L-THIA model using the millions of curve number (CN) combinations associated with land uses and hydrologic soil groups. L-THIA calibration for the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana was performed using land use data for 1991 and daily rainfall data for six months of 1991 (January 1-June 30) to minimize errors associated with use of different temporal land use data and rainfall data. For the calibration period, the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The calibrated CN values were used for validation of the model for the same year (July 1-December 31), and the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The Nash-Sutcliffe coefficient was 0.52 for January 1, 1991 to December 31, 1991 using uncalibrated CN values. As shown in this study, the use of better input parameters for the L-THIA model can improve accuracy. The effects on direct runoff and pollutant estimation of the calibrated CN values in the L-THIA model were investigated for the LEC. Following calibration, the estimated average annual direct runoff for the LEC watershed increased by 34%, total nitrogen by 24%, total phosphorus by 22%, and total lead by 43%. This study demonstrates that the L-THIA model should be calibrated and validated prior to application in a particular watershed to more accurately assess the effects of land use changes on hydrology and water quality.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

Thermal Characteristics of Stormwater Runoff from Asphalt and Sod Surfaces1

Abstract: Urban impervious surfaces absorb and store thermal energy, particularly during warm summer months. During a rainfall/runoff event, thermal energy is transferred from the impervious surface to the runoff, causing it to become warmer. As this higher temperature runoff enters receiving waters, it can be harmful to coldwater habitat. In an urban watershed, impervious asphalt surfaces (roads, parking lots, and driveways) and pervious residential lawns comprise a significant portion of the watershed area. A paired asphalt‐turfgrass sod plot was constructed to compare the thermal runoff characteristics between asphalt and turfgrass sod surfaces, to identify meteorological variables that influence these thermal characteristics, and to evaluate evaporative heat loss for runoff from asphalt surfaces. Rainfall simulations were conducted during the summers of 2004 and 2005 under a range of climatic conditions. Asphalt surface temperatures immediately prior to rainfall simulations averaged 43.6°C and decreased an average of 12.3°C over 60 min as rain cooled the surface. In contrast, presimulation sod surface temperatures averaged only 23.3°C and increased an average of 1.3°C throughout the rainfall events. Heat transferred from the asphalt to the runoff resulted in initial asphalt runoff temperatures averaging 35.0°C that decreased by an average of 4.1°C at the end of the event. Sod runoff temperatures averaged only 25.5°C and remained fairly constant throughout the simulations. Multivariable regression equations were developed to predict (1) average asphalt surface temperature (R² = 0.90) and average asphalt runoff temperature (R² = 0.92) as a function of solar radiation, rain temperature, and wind speed, and (2) average sod surface temperature (R² = 0.85) and average sod runoff temperature (R² = 0.94) as a function of solar radiation, rain temperature, rain intensity, and wind speed. Based on a heat balance analysis, existing evaporation equations developed from studies on lakes were not adequate to predict evaporation from runoff on a heated impervious surface. The combined heat from the asphalt and sod plots was an average of 38% less than the total heat had the total area consisted solely of asphalt.     

Soil nitrate nitrogen dynamics after biosolids application in a tobosagrass desert grassland

Dormant-season application of biosolids increases desert grass production more than growing season application in the first growing season after application. Differential patterns of NO3-N (plant available N) release following seasonal biosolids application may explain this response. Experiments were conducted to determine soil nitrate nitrogen dynamics following application of biosolids during two seasons in a tobosagrass [Hilaria mutica (Buckl.) Benth.] Chihuahuan Desert grassland. Biosolids were applied either in the dormant (early April) or growing (early July) season at 0, 18, or 34 dry Mg ha(-1). A polyester-nylon mulch was also applied to serve as a control that approximated the same physical effects on the soil surface as the biosolids but without any chemical effects. Supplemental irrigation was applied to half of the plots. Soil NO3-N was measured at two depths (0-5 and 5-15 cm) underneath biosolids (or mulch) and in interspace positions relative to surface location of biosolids (or mulch). Dormant-season biosolids application significantly increased soil NO3-N during the first growing season, and also increased soil NO3-N throughout the first growing season compared to growing-season biosolids application in a year of higher-than-average spring precipitation. In a year of lower-than-average spring precipitation, season of application did not affect soil NO3-N. Soil NO3-N was higher at both biosolids rates for both seasons of application than in the control treatment. Biosolids increased soil NO3-N compared to the inert mulch. Irrigation did not significantly affect soil NO3-N. Soil NO3-N was not significantly different underneath biosolids and in interspace positions. Surface soil NO3-N was higher during the first year of biosolids application, and subsurface soil NO3-N increased during the second year. Results showed that biosolids rate and season of application affected soil NO3-N measured during the growing season. Under dry spring-normal summer precipitation conditions, season of application did not affect soil NO3-N; in contrast, dormant season application increased soil NO3-N more than growing season application under wet spring-dry summer conditions.     

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.