Drainage Impacts on Surficial Water Retention Capacity of a Prairie Pothole Watershed

Wetland restoration has been proposed as a tool to mitigate excess runoff and associated nonpoint source pollution in the Upper Midwestern United States. This study quantified the surficial water retention capacity of existing and drained wetlands for the Greater Blue Earth River Basin (GBERB), an intensively drained agricultural watershed. Using airborne light detection and ranging, the historic depressional storage was determined to be 152 mm. Individual depression analysis suggested that the restoration of most drained areas would have little impact on the storage capacity of the GBERB because the majority (53%) of retention capacity was in large depressions (>40 ha) which comprised only a small proportion (<1.0) of the observed depressions. Accounting for change in storage and the difference in annual evapotranspiration (ET) between wetlands and the croplands that replaced them, restoration of all depressions in the Minnesota portion of GBERB would provide a maximum of 131 mm additional capacity over and above the modern day capacity (193 mm; 56 mm depressional storage; 60 mm wetland ET; and 77 mm cropland ET). Considering that depressional depths in smaller areas are within the range of uncertainty of the lidar digital elevation models and larger depressions have the most storage, we conclude that efforts to increase the surficial water‐holding capacity of the GBERB would be best served in the restoration of large (>40 ha) depressions.     

Depressional runoff cascade networks of the Des Moines Lobe of Iowa

In this study, we develop a general mathematical framework and algorithm for routing cumulative precipitation excess through depressional fill–spill cascade networks in a landscape using only information about depression morphology, local contributing areas, and potential overland flow pathways. The framework also allows for the classification of depressions according to their landscape position within a network, and calculation of precipitation- and non-precipitation-dependent network properties, including measures of network complexity and runoff connectivity. To demonstrate its use, we applied our framework to the 167,287 drained depressions of the Des Moines Lobe of Iowa, a sub-region of the Prairie Pothole Region of North America, over a large range of historically observed precipitation amounts for scenarios both neglecting and incorporating infiltration in runoff generation. Our results show that 85.3% of depressions in this region form 18,851 unique depressional runoff cascade networks, with the remainder being disjunct features. Most of the properties of the region's networks appear to conform to either a truncated power-law or lognormal distribution. For a given rainfall amount, surface runoff connectivity between depressions within networks, and between networks and off-network areas, is controlled primarily by available aggregate depressional volumetric storage and contributing area, and to a lesser degree, network complexity.     

Rainfall Runoff Relationships for a Cloud Forest Watershed in Central America: Implications for Water Resource Engineering1

Caballero, Luis A., Alon Rimmer, Zachary M. Easton, and Tammo S. Steenhuis, 2012. Rainfall Runoff Relationships for a Cloud Forest Watershed in Central America: Implications for Water Resource Engineering. Journal of the American Water Resources Association (JAWRA) 48(5): 1022‐1031. DOI: 10.1111/j.1752‐1688.2012.00668.x Abstract: Understanding the basic relationships between rainfall and runoff is vital for effective management and utilization of scarce water resources. Especially, this is important in Central America with widespread potable water shortage during the dry months of the monsoon. Potential good water sources are cloud forests, but little information concerning its potential is available to water supply engineers. Our objective is to define rainfall‐runoff‐base flow relationships for a cloud forest catchment. Flumes were installed for measuring river flow in four subwatersheds in La Tigra National Park, Honduras. One of the four watersheds was a 636‐ha subwatershed (WS1) with 60% cloud forest coverage. Precipitation averaged 1,130 mm/yr over the entire basin. About half of the total rainfall became runoff for the cloud forest watershed whereas, for the adjacent undisturbed forested watershed, the total discharge was <20% of the amount of precipitation. Infiltration rates were generally greater than rainfall rates. Therefore, most rainfall infiltrated into the soil, especially in the upper, steep, and well‐drained portions of the watershed. Direct runoff was generated from saturated areas near the river and exposed bedrock. This research provides compelling evidence that base flow is the primary contributor to streamflow during both wet and dry seasons in cloud forest catchments. Protecting these flow processes over time is critical for the sustained provision of potable water.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

Integrating Temporal Inequality into Conservation Planning to Improve Practice Design and Efficacy

In contrast to spatial inequality, there are currently no methods for leveraging information on temporal inequality to improve conservation efficacy. The objective of this study was to use Lorenz curves to quantify temporal inequality in surface runoff and tile drainage, identify controls on nutrient loading in these flowpaths, and develop design flows for structural conservation practices. Surface runoff (n = 94 site‐years) and tile drainage (n = 90 site‐years) were monitored on 40 fields in Ohio. Results showed, on average, 80% of nitrate‐nitrogen, soluble reactive phosphorus (P), and total P loads occurred between 7 and 12 days per year in surface runoff and between 32 and 58 days per year in tile drainage. Similar temporal inequality between discharge and load provided evidence that loading was transport‐limited and highlighted the critical role hydrologic connectivity plays in nutrient delivery from tile‐drained fields. Design flow criterion for sizing structural practices based on load reduction goals was developed by combining Lorenz curves and flow duration curves. Comparing temporal inequality between fields and the Maumee River, the largest tributary to the western Lake Erie Basin, revealed challenges associated with achieving watershed load reduction goals with field‐scale conservation. In‐field (i.e., improved nutrient and water management), edge‐of‐field (i.e., structural practices), and instream practices will all be required to meet nutrient reduction goals from tile‐drained watersheds.     

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.     

APPLICATION OF RORB MODEL TO A CATCHMENT IN SINGAPORE1

ABSTRACT: Runoff Routing model (RORB) is a general model applicable to both rural and urban catchments. The performance of the model is illustrated through its simulation of flood runoff hydrographs in an urban catchment in Singapore. The essential feature of the model is the routing of rainfall excesses on subareas through some arrangement of concentrated storage elements, which represent the distribution of temporary storage of flood runoff on the watershed. This nonlinear routing procedure of the storage elements has two common parameters, k_c and m. With the limited data available, these two parameter values were determined through calibration runs. The same set of values of k_c and m were then used in the model to determine the runoff hydrographs of five other storms selected from the rainfall events between 1979 and 1981. It was found that the simulated runoff hydrographs matched reasonably well with the recorded hydrographs.     

A New Open‐Access HUC‐8 Based Downscaled CMIP‐5 Climate Model Forecast Dataset for the Conterminous United States

Watershed‐scale hydrologic simulation models generally require climate data inputs including precipitation and temperature. These climate inputs can be derived from downscaled global climate simulations which have the potential to drive runoff forecasts at the scale of local watersheds. While a simulation designed to drive a local watershed model would ideally be constructed at an appropriate scale, global climate simulations are, by definition, arbitrarily determined large rectangular spatial grids. This paper addresses the technical challenge of making climate simulation model results readily available in the form of downscaled datasets that can be used for watershed scale models. Specifically, we present the development and deployment of a new Coupled Model Intercomparison Project phase 5 (CMIP5) based database which has been prepared through a scaling and weighted averaging process for use at the level of U.S. Geological Survey (USGS) Hydrologic Unit Code (HUC)‐8 watersheds. The resulting dataset includes 2,106 virtual observation sites (watershed centroids) each with 698 associated time series datasets representing average monthly temperature and precipitation between 1950 and 2099 based on 234 unique climate model simulations. The new dataset is deployed on a HydroServer and distributed using WaterOneFlow web services in the WaterML format. These methods can be adapted for downscaled General Circulation Model (GCM) results for specific drainage areas smaller than HUC‐8. Two example use cases for the dataset also are presented.     

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.     

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

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

A Hydrologic Data Screening Procedure for Exploring Monotonic Trends and Shifts in Rainfall and Runoff Patterns

A thorough understanding of past and present hydrologic responses to changes in precipitation patterns is crucial for predicting future conditions. The main objectives of this study were to determine temporal changes in rainfall‐runoff relationship and to identify significant trends and abrupt shifts in rainfall and runoff time series. Ninety‐year rainfall and runoff time series datasets from the Gasconade and Meramec watersheds in east‐central Missouri were used to develop data screening procedure to assess changes in the rainfall and runoff temporal patterns. A statistically significant change in mean and variance was detected in 1980 in the rainfall and runoff time series within both watersheds. In addition, both the rainfall and runoff time series indicated the presence of nonstationary attributes such as statistically significant monotonic trends and/or change in mean and variance, which should be taken into consideration when using the time series to predict future scenarios. The annual peak runoff and the annual low flow in the Meramec watershed showed significant temporal changes compared to that in the Gasconade watershed. Water loss in both watersheds was found to be significantly increasing which is potentially due to the increase in groundwater pumping for water supply purposes.     

Impacts of Climate Variability and Land Use Alterations on Frequency Distributions of Terrestrial Runoff Loading to Coastal Waters in Southern California1

Abstract: The transport of water, sediment, dissolved and particulate chemicals, and bacteria from coastal watersheds affects the nearshore marine and estuarine waters. In southern California, coastal watersheds deliver water and associated constituents to the nearshore system in discrete pulses. To better understand the pulsed nature of these watersheds, frequency distributions of simulated runoff events are presented for: (1) three land use conditions (1929, 1998, 2050); (2) three time periods (all water years 1989‐2002), only El Nino years (1992, 1993, 1995, 1998); and only non‐El Nino years; and (3) three regions (watershed, uplands, and lowlands). At the watershed scale, there was a significant increase (>200%) in mean event runoff from 1929 to 2050 (0.4‐1.3 cm) due to localized urbanization, which shifted the dominant sources of runoff from the mountains in 1929 (78% of watershed runoff) to the coastal plane for 2050 conditions (51% of watershed runoff). Inter‐annual climate variability was strong in the rainfall and runoff frequency distributions, with mean event rainfall and runoff 66 and 60% larger in El Nino relative to non‐El Nino years. Combining urbanization and climate variability, 2050 land conditions resulted in El Nino years being five times more likely to produce large (>3.0 cm) runoff events relative to non‐El Nino years. Combining frequency distributions of event runoff with regional nutrient export relationships, we show that in El Nino years, one in five events produced runoff ≥2.5 cm and temporary nearshore nitrate and phosphate concentrations of 12 and 1.4 μM, respectively, or approximately 5‐10 times above ambient conditions.     

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.     

Impact of land drainage on peatland hydrology

There is a long history of drainage of blanket peat but few studies of the long-term hydrological impact of drainage. This paper aims to test differences in runoff production processes between intact and drained blanket peat catchments and determine whether there have been any long-term changes in stream flow since drainage occurred. Hillslope runoff processes and stream discharge were measured in four blanket peat catchments. Two catchments were drained with open-cut ditches in the 1950s. Ditching originally resulted in shorter lag times and flashier storm hydrographs but no change in the annual catchment runoff efficiency. In the period between 2002 and 2004, the hydrographs in the drained catchments, while still flashy, were less sensitive to rainfall than in the 1950s and the runoff efficiency had significantly increased. Drains resulted in a distinctive spatial pattern of runoff production across the slopes. Overland flow was significantly lower in the drained catchments where throughflow was more dominant. In the intact peatlands, matrix throughflow produced by peat layers below 10 cm was rare and produced <1% of the runoff. However, in drained peatlands, matrix throughflow in deeper peat layers was common and provided around 23% of the runoff from gauged plots. Macropore flow, the density of soil piping, and pipeflow were significantly greater in drained peatlands than in intact basins. Gradual changes to peat structure could explain the long-term changes in river flow, which are in addition to those occurring in the immediate aftermath of peatland drainage.     

An Application of the T‐TEL Assessment Method to Evaluate Connectivity in a Lake‐Dominated Watershed after Drought

Lakes are landscape features that influence connectivity of mass and energy by being foci for the reception, mixing, and provision of water and material. Where lake fractions are high, they influence hydrological connectivity. This behavior was exemplified in the Baker Creek watershed in Canada's Northwest Territories during a two‐year drought in which many lake levels declined below outlet elevations. This study evaluated how lakes controlled surface runoff connectivity reestablishment following the drought using a new assessment method, T‐TEL (time scales — thresholds, excesses, losses). Analysis of daily data showed that during a summer period following the drought, connectivity occurred between 0% and 41% of the time. The size of run‐of‐the‐river lakes relative to their upstream watershed area, and the upstream lake fraction, are two factors for connectivity. These terms represent a lake's ability to control the size of storage deficits relative to rainfall, and evaporation and storage losses along pathways. The connectivity magnitude–duration curve only aligned with the watershed flow duration curve during high‐water conditions, implying lakes functioned as individuals rather than as part of a perennial watercourse during much of the study. The T‐TEL method can be used to quantify consistent metrics of hydrologic connectivity that can be used for regionalization exercises and understanding hydrologic controls on material transport.     

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.     

Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain1

Epps, Thomas H., Daniel R. Hitchcock, Anand D. Jayakaran, Drake R. Loflin, Thomas M. Williams, and Devendra M. Amatya, 2012. Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12000 Abstract: Hydrologic monitoring was conducted in two first‐order lower coastal plain watersheds in South Carolina, United States, a region with increasing growth and land use change. Storm events over a three‐year period were analyzed for direct runoff coefficients (ROC) and the total storm response (TSR) as percent rainfall. ROC calculations utilized an empirical hydrograph separation method that partitioned total streamflow into sustained base flow and direct runoff components. ROC ratios ranged from 0 to 0.32 on the Upper Debidue Creek (UDC) watershed and 0 to 0.57 on Watershed 80 (WS80); TSR results ranged from 0 to 0.93 at UDC and 0.01 to 0.74 at WS80. Variability in event runoff generation was attributed to seasonal trends in water table elevation fluctuation as regulated by evapotranspiration. Groundwater elevation breakpoints for each watershed were identified based on antecedent water table elevation, streamflow, ROCs, and TSRs. These thresholds represent the groundwater elevation above which event runoff generation increased sharply in response to rainfall. For effective coastal land use decision making, baseline watershed hydrology must be understood to serve as a benchmark for management goals, based on both seasonal and event‐based surface and groundwater interactions.     

Nitrate Removal in a Restored Spring‐Fed Wetland,Pennsylvania, USA1

ABSTRACT:  In 2001, the 1.04‐ha Hornbaker wetland in south‐central Pennsylvania was restored by blocking an artificial drainage ditch to increase water storage and hydraulic retention time (HRT). A primary goal was to diminish downstream delivery of nitrate that enters the wetland from a limestone spring, its main source of inflow. Wetland inflow and outflow were monitored weekly for two years to assess nitrate flux, water temperature, pH, and specific conductivity. In Year 2, spring discharge was measured weekly to allow calculation of nitrate loads and hydraulic retention time. Surface runoff was confirmed to be a small fraction of wetland inflows via rainfall‐runoff modeling with TR‐55. The full dataset (n = 102) was screened to remove 13 weeks in which spring discharge constituted < 85% of total inflows because of high precipitation and surface runoff. Over two years (n = 89), mean nitrate‐nitrogen concentrations were 7.89 mg/l in inflow and 3.68 mg/l in outflow, with a mean nitrate removal of 4.19 mg/l. During Year 2 (n = 47), for which nitrate load data were available, the wetland removed an average of 2.32 kg N/day, 65% of the load. Nitrate removal was significantly correlated with HRT, water temperature, and the concentration of nitrate in inflow and was significantly greater during the growing season (5.36 mg/l, 64%) than during the non‐growing season (3.23 mg/l, 43%). This study indicates that hydrologic restoration of formerly drained wetlands can provide substantial water quality benefits and that the hydrologic characteristics of spring‐fed wetlands, in particular, support effective nitrogen removal.     

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

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

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.