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.     

Integrating contributing areas and indexing phosphorus loss from agricultural watersheds

Most states in the USA have adopted P Indexing to guide P-based management of agricultural fields by identifying the relative risk of P loss at farm and watershed scales. To a large extent, this risk is based on hydrologic principles that frequently occurring storms can initiate surface runoff from fields. Once initiated, this hydrological pathway has a high potential to transport P to the stream. In regions where hydrologically active areas of watersheds vary in time and space, surface runoff generation by "saturation excess" has been linked to distance from stream, with larger events resulting in larger contributing distances. Thus, storm-return period and P loss from a 39.5-ha mixed-land-use watershed in Pennsylvania was evaluated to relate return-period thresholds and distances contributing P to streams. Of 248 storm flows between 1997 and 2006, 93% had a return period of 1 yr, contributing 47% of total P (TP) export, while the largest two storms (10-yr return period) accounted for 23% of TP export. Contributing distance thresholds for the watershed were determined (50-150 m) for a range of storm-return periods (1-10 yr) from hydrograph analysis. By modifying storm-return period thresholds in the P Index and thereby contributing distance, it is possible to account for greater risk of P loss during large storms. For instance, increasing return period threshold from 1 (current P indices) to 5 yr, which accounted for 67% of TP export, increased the P-management restricted area from 20 to 58% of the watershed. An increase in impacted area relative to a decreased risk of P loss creates a management-policy dilemma that cannot be ignored.     

COMPARISON OF DESIGN STORM CONCEPTS USING CONTINUOUS SIMULATION WITH SHORT DURATION STORMS1

ABSTRACT: The objectives of this paper were to test the ability of various design storm distributions to simulate the actual rainfall pattern and to compare the runoff rates used in the design of stormwater management devices in the State of Florida using continuous simulation approach. The analyses were performed for four gaged stations to evaluate the applicability of design storm distributions in different parts of the State of Florida. The approach used in this study compared the peak runoff rates from design storms based on the various distributions to those that would result from actual rainfall events. A series of continuous runoff rates were developed through the use of actual fifteen-minute recorded rainfall data, Horton type infiltration decay and recovery rate, and a continuous simulation model. The runoff rates were analyzed using frequency distributions to obtain peak runoff rates associated with different return periods based on the assumption that the continuous simulation approach closely predicts the actual runoff rates from the gaged stations. The results show that the behavior of the design storm distributions varies for different watershed characteristics in different parts of the state. The study also suggests that in general the Florida Department of Transportation and the Suwanne River Water Management (FDOT/ SRWMD) distributions appeared to agree with the continuous simulation results.     

Impacts of Climate Change on Extreme Precipitation Events Over Flamingo Tropicana Watershed1

Abstract: Climate change, particularly the projected changes to precipitation patterns, is likely to affect runoff both regionally and temporally. Extreme rainfall events are expected to become more intense in the future in arid urban areas and this will likely lead to higher streamflow. Through hydrological modeling, this article simulates an urban basin response to the most intense storm under anthropogenic climate change conditions. This study performs an event‐based simulation for shorter duration storms in the Flamingo Tropicana (FT) watershed in Las Vegas, Nevada. An extreme storm, defined as a 100‐year return period storm, is selected from historical records and perturbed to future climatic conditions with respect to multimodel multiscenario (A1B, A2, B1) bias corrected and spatially disaggregated data from the World Climate Research Programme's (WCRP's) database. The cumulative annual precipitation for each 30‐year period shows a continuous decrease from 2011 to 2099; however, the summer convective storms, which are considered as extreme storms for the study area, are expected to be more intense in future. Extreme storm events show larger changes in streamflow under different climate scenarios and time periods. The simulated peak streamflow and total runoff volume shows an increase from 40% to more than 150% (during 2041‐2099) for different climate scenarios. This type of analysis can help evaluate the vulnerability of existing flood control system and flood control policies.     

COMPOSITE VS. DISTRIBUTED CURVE NUMBERS: EFFECTS ON ESTIMATES OF STORM RUNOFF DEPTHS1

ABSTRACT: The U.S. Department of Agriculture Curve Number (CN) method is one of the most common and widely used techniques for estimating surface runoff and has been incorporated into a number of popular hydrologic models. The CN method has traditionally been applied using compositing techniques in which the area weighted average of all curve numbers is calculated for a watershed or a small number of sub-watersheds. CN compositing was originally developed as a time saving procedure, reducing the number of runoff calculations required. However, with the proliferation of high speed computers and geographic information systems, it is now feasible to use distributed CNs when applying the CN method. To determine the effect of using composited versus distributed CNs on runoff estimates, two simulations of idealized watersheds were developed to compare runoff depths using composite and distributed CNs. The results of these simulations were compared to the results of similar analyses performed on an urbanizing watershed located in central Indiana and show that runoff depth estimates using distributed CNs are as much as 100 percent higher than when composited CNs are used. Underestimation of runoff due to CN compositing is a result of the curvilinear relationship between CN and runoff depth and is most severe for wide CN ranges, low CN values, and low precipitation depths. For larger design storms, however, the difference in runoff computed using composite and distributed CNs is minimal.     

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.     

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.     

EVALUATION OF HYDROLOGIC BENEFITS OF INFILTRATION BASED URBAN STORM WATER MANAGEMENT1

ABSTRACT: As watersheds are urbanized, their surfaces are made less pervious and more channelized, which reduces infiltration and speeds up the removal of excess runoff. Traditional storm water management seeks to remove runoff as quickly as possible, gathering excess runoff in detention basins for peak reduction where necessary. In contrast, more recently developed “low impact” alternatives manage rainfall where it falls, through a combination of enhancing infiltration properties of pervious areas and rerouting impervious runoff across pervious areas to allow an opportunity for infiltration. In this paper, we investigate the potential for reducing the hydrologic impacts of urbanization by using infiltration based, low impact storm water management. We describe a group of preliminary experiments using relatively simple engineering tools to compare three basic scenarios of development: an undeveloped landscape; a fully developed landscape using traditional, high impact storm water management; and a fully developed landscape using infiltration based, low impact design. Based on these experiments, it appears that by manipulating the layout of urbanized landscapes, it is possible to reduce impacts on hydrology relative to traditional, fully connected storm water systems. However, the amount of reduction in impact is sensitive to both rainfall event size and soil texture, with greatest reductions being possible for small, relatively frequent rainfall events and more pervious soil textures. Thus, low impact techniques appear to provide a valuable tool for reducing runoff for the events that see the greatest relative increases from urbanization: those generated by the small, relatively frequent rainfall events that are small enough to produce little or no runoff from pervious surfaces, but produce runoff from impervious areas. However, it is clear that there still needs to be measures in place for flood management for larger, more intense, and relatively rarer storm events, which are capable of producing significant runoff even for undeveloped basins.     

HYDROLOGIC BEHAVIOR OF VEGETATED ROOFS1

ABSTRACT: Control of stormwater runoff from impervious surfaces is an important national goal because of disruptions to downstream ecosystems, water users, and property owners caused by increased flows and degraded quality. One method for reducing stormwater is the use of vegetated (green) roofs, which efficiently detain and retain stormwater when compared to conventional (black) roofs. A paired green roof‐black roof test plot was constructed at the University of Georgia and monitored between November 2003 and November 2004 for the green roof's effectiveness in reducing stormwater flows. Stormwater mitigation performance was monitored for 31 precipitation events, which ranged in depth from 0.28 to 8.43 cm. Green roof precipitation retention decreased with precipitation depth; ranging from just under 90 percent for small storms (< 2.54 cm) to slightly less than 50 percent for larger storms (> 7.62 cm). Runoff from the green roof was delayed; average runoff lag times increased from 17.0 minutes for the black roof to 34.9 minutes for the green roof, an average increase of 17.9 minutes. Precipitation and runoff data were used to estimate the green roof curve number, CN = 86. This information can be used in hydrologic models for developing stormwater mitigation programs.     

WATERSHED OPTIMIZATION OF AGRICULTURAL BEST MANAGEMENT PRACTICES: CONTINUOUS SIMULATION VERSUS DESIGN STORMS1

ABSTRACT: Nonpoint source (NPS) models and expert opinions are often used to prescribe best management practices (BMPs) for controlling NPS pollution. An optimization algorithm (e.g., a genetic algorithm, or GA) linked with a NPS model (e.g., Annualized AGricultural Nonpoint Source pollution model, or AnnAGNPS), can be used to more objectively prescribe BMPs and to optimize NPS pollution control measures by maximizing pollutant reduction and net monetary return from a watershed. Pollutant loads from design storms and annual loads from a continuous simulation can both be used for optimizing BMP schemes. However, which strategy results in a better solution (in terms of providing water quality protection) for a watershed is not clear. The specific objective of the study was to determine the differences in watershed pollutant loads, in an experimental watershed in Pennsylvania, resulting from optimization analyses performed using pollutant loads from a series of five 2‐yr 24‐hr storm events, a series of five 5‐yr 24‐hr storm events, and cumulative pollutant loads from a continuous simulation of five years of weather data. For each of these three different event alternatives, 100 near optimal solutions (BMP schemes) were generated. Sediment (Sed), sediment nitrogen (SedN), dissolved N (SolN), sediment organic carbon (SedOC), and sediment phosphorus (SedP) loads from a different five‐year period (an evaluation period) suggest that the optimal BMP schemes resulting from the use of annual cumulative pollutant loads from a continuous simulation of five years of weather data provide smaller cumulative NPS pollutant loads at the watershed outlet.     

THE APPLICATION OF HYDROLOGIC MODELS TO SMALL WATERSHEDS HAVING MILD TOPOGRAPHY1

ABSTRACT: The application of hydrologic models to small watersheds of mild topography is not well documented. This study evaluates the applicability of hydrologic models described by Huggins and the Soil Conservation Service to small watersheds by comparing the simulated and actual hydrograph for both gaged and ungaged situations. The annual maximum rainfall events plus storms exceeding 2.5 inches from 25 years of rainfall and runoff data for two small watersheds were selected for the model evaluations. These storms had a variety of patterns and occurred on many different watershed conditions. Simulated and actual hydrographs were compared using a parameter which contained volume, peak, and shape factors. One-half of the selected storms were used to calibrate the models. For both models, there were no significant differences between the simulated and actual runoff volumes and peak runoff rates. Parameters obtained during the calibration process and relationships developed to estimate antecedent moisture and to modify tabulated runoff curve numbers were used to simulate the runoff hydrograph from the remaining storms. These remaining storms or test storms were simulated only once in order to imitate an ungaged situation. In general, both the Huggins and SCS model performed similarly on the test storms, but the level of model performance was lower than that for the calibration storms. For both models, the two-day antecedent rainfall was more important than the five-day in determining antecedent moisture and modifying tabulated curve numbers. The time of concentration which resulted in good hydrograph simulations was about three times larger than that estimated using published empirical relationships.     

EFFECTS OF STORM PATTERNS ON RUNOFF HYDROGRAPHS1

ABSTRACT: As an alternative to the conventional single-peak design storms commonly used in hydrologic practice, a large number of Southeastern Pennsylvania storm events were selected from hourly U.S. National Oceanographic and Atmospheric Administration (NOAA) records, and their temporal distributions were analyzed. From these recorded events, design storms of a typical distribution were developed for storm durations between 6 and 18 hours. All of these generated design storms have two or more peaks. The conventional single peak as well as the “typical” multi-peak storms were then applied to a simulated watershed. It was found that the multi-peak storms consistently produced more dispersed hydrographs with lower runoff peaks than the conventional single peak storms.     

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.     

RELIABILITY OF THE DESIGN STORM CONCEPT IN EVALUATING RUNOFF PEAK FLOW1

A comparative study was undertaken to evaluate peak runoff flow rates using (1) a continuous series of actual rainfall events and (2) design storms. The ILLUDAS computer model was used to simulate runoff over a catchment within the city of Montreal, Canada. A ten-year period, five-minute increment rainfall data base was used to derive peak flow frequency curves. Two types of design storms were analyzed: one derived from intensity duration frequency curves (Chicago type), the other from averaging actual rainfall patterns (Huff type). Antecedent soil moisture conditions were considered in the analyses. It was found that the probability distribution of runoff peak flow was sensitive to the choice of design storm pattern and to the antecedent soil moisture condition. A symmetrical, Chicago-type design storm with antecedent dry soil moisture produced a flow frequency curve similar to the one obtained from a series of historical rainfall events.     

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.     

SELECTED CHALLENGES IN RUNOFF GENERATION RESEARCH IN FORESTS FROM THE HILLSLOPE TO HEADWATER DRAINAGE BASIN SCALE1

ABSTRACT: When faced with practical forest land management issues such as the impacts of logging or forest conversion to other land uses, planners ideally require a comprehensive understanding of within drainage basin hydrological processes to determine the most vulnerable areas to increased storm runoff and erosion. Land managers in particular need to know the source areas and magnitude of inputs to the storm hydrograph, in terms of water quantity, sediment and solute transport; and the routing of such hydrographs from headwater to larger drainage basins. The latter includes an overall assessment at various scales of the impacts of forest disturbance and conversion on the water balance. This paper will focus on runoff generation in terms of identifying the various pathways and source areas. Such aspects will be linked with the need for a more comprehensive effort towards the field testing of so-called ‘physically based’ models of runoff generation. Some of the controversial issues arising from the difficulties in reconciling results from hydrochemical investigations with complementary hydrometric studies will be highlighted. Subsequently, attention will be given to topographic-wetness models, which have promising applications in forestland management. In addition, alternative simple models for application at the catchment scale will be assessed. The latter is in recognition that at smaller scales, heterogeneity both in time and space of soil hydraulic properties demand a greater number of parameters in modelling. Such considerations can even prove an obstacle in terms of the confident application of ‘physically based’ models.     

CALIBRATION OF STORM LOADS IN THE SOUTH PRONG WATERSHED,FLORIDA, USING BASINS/HSPF1

ABSTRACT: The South Prong watershed is a major tributary system of the Sebastian River and adjacent Indian River Lagoon. Continued urbanization of the Sebastian River drainage basin and other watersheds of the Indian River Lagoon is expected to increase runoff and nonpoint source pollutant loads. The St. Johns River Water Management District developed watershed simulation models to estimate potential impacts on the ecological systems of receiving waters and to assist planners in devising strategies to prevent further degradation of water resources. In the South Prong system, a storm water sampling program was carried out to calibrate the water quality components of the watershed model for total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN). During the period of May to November 1999, water quality and flow data were collected at three locations within the watershed. Two of the sampling stations were located at the downstream end of major watercourses. The third station was located at the watershed outlet. Five storm events were sampled and measured at each station. Sampling was conducted at appropriate intervals to represent the rising limb, peak, and recession limb of each storm event. The simulations were handled by HSPF (Hydrologic Simulation Program‐Fortran). Results include calibration of the hydrology and calibration of the individual storm loads. The hydrologic calibration was continuous over the period 1994 through 1999. Simulated storm runoff, storm loads, and event mean concentrations were compared with their corresponding observed values. The hydrologic calibration showed good results. The outcome of the individual storm calibrations was mixed. Overall, however, the simulated storm loads agreed reasonably well with measured loads for a majority of the storms.     

成都市某区暴雨径流过程模拟分析

为提高城市雨洪管理的效率,最大限度地减少暴雨洪水带来的危害,针对城市防洪排涝的需要,在分析成都市某区降雨径流规律后,建立了该区暴雨径流模拟的数学模型,有助于采取相应措施充分利用雨洪资源。通过对雨洪过程模拟验证表明,模型适合该区域的实际情况,具有一定的合理性和可靠性。     

STORM CHARACTERISTICS OF CONVECTIVE-SCALE PRECIPITATION1

A classification scheme for convective precipitation, having applications in both analysis and modeling of meteorological and hydrological events, is presented. The method is based upon observations of rainfall at the ground, radar scans of storm events, and visible and infrared satellite imagery of larger storm systems. Empirical and theoretical frequency distributions are derived for total storm rainfall, rainfall duration and time between storms for each of the convective categories. This stratification is directly applicable to the experimental design and evaluation of weather modification projects and may be useful for the development and interpretation of meteorological and hydrological models. When atmospheric conditions limit storm development to cells, rainfall was seldom observed. Small clusters also produce small amounts of rainfall but have a longer lifetime than cells and are likely candidates for cloud seeding attempts to encourage their growth to large clusters. Large and nested clusters usually produce large amounts of natural precipitation. A few large storms account for most of a season's rainfall.     

THE EFFECTS OF A MOVING RAINSTORM ON DIRECT RUNOFF PROPERTIES1

ABSTRACT: The effects of a moving rainstorm on flood runoff characteristics were investigated. A flood hydrograph simulation model called “FH-Model” and a natural watershed were used. A hypothetical rainstorm of 50 years recurrence interval, 75 mm depth, and 4 hours duration was used to show the effects of velocity and direction of the moving rainstorm on the runoff characteristics. Compared with an equivalent stationary rainstorm (ESRS), the peak flow caused by a rainstorm moving in a downstream direction with a speed equal to channel velocity, V, was 27.5 percent higher and the peak flow caused by the same rainstorm moving in an upstream direction was 21.7 percent smaller. These percentages reduced to 10.5 percent and 8.6 percent for storms moving downstream and upstream, respectively, at three times the channel velocity, 3V. There were negligible differences in the time of peak, T_p between runoff caused by storms moving downstream and runoff produced by ESRS. However, T_p for a storm moving upstream at V velocity was 82 percent higher than that produced by ESRS, but was reduced to 27 percent higher when the storm velocity was 3V.