CLIMATE FACTOR FOR SMALL-BASIN FLOOD FREQUENCY1

ABSTRACT: A climate factor, C_T, (T = 2–, 25-, and 100-year recurrence intervals) that delineates regional trends in small-basin flood frequency was derived using data from 71 long-term rainfall record sites. Values of C_T at these sites were developed by a regression analysis that related rainfall-runoff model estimates of T-year floods to a sample set of 50 model calibrations. C_T was regionalized via kriging to develop maps depicting its geographic variation for a large part of the United States east of the 105th meridian. Kriged estimates of C_T and basin-runoff characteristics were used to compute regionalized T-year floods for 200 small drainage basins. Observed T-year flood estimates also were developed for these sites. Regionalized floods are shown to account for a large percentage of the variability in observed flood estimates with coefficients of determination ranging from 0.89 for 2-year floods to 0.82 for 100-year floods. The relative importance of the factors comprising regionalized flood estimates is evaluated in terms of scale (size of drainage area), basin-runoff characteristics (rainfall. runoff model parameters), and climate (C_T).     

Long‐Term High‐Resolution Radar Rainfall Fields for Urban Hydrology

Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km²) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km², suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km²) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.     

Modeling Stream Channel Characteristics From Drainage‐Enforced DEMs in Puget Sound,Washington, USA1

Abstract: Mapping stream channels and their geomorphic attributes is an important step in many watershed research and management projects. Often insufficient field data exist to map hydromorphologic attributes across entire drainage basins, necessitating the application of hydrologic modeling tools to digital elevation models (DEMs) via a geographic information system (GIS). In this article, we demonstrate methods for deriving synthetic stream networks via GIS across large and diverse basins using drainage‐enforced DEMs, along with techniques for estimating channel widths and gradient on the reach scale. The two‐step drainage enforcement method we used produced synthetic stream networks that displayed a high degree of positional accuracy relative to the input streams. The accuracies of our estimated channel parameters were assessed with field data, and predictions of bankfull width, wetted width and gradient were strongly correlated with measured values (r² = 0.92, r² = 0.95, r² = 0.88, respectively). Classification accuracies of binned channel attributes were also high. Our methodology allows for the relatively rapid mapping of stream channels and associated morphological attributes across large geographic areas. Although initially developed to provide salmon recovery planners with important salmon habitat information, we suggest these methodologies are relevant to a variety of research and management questions.     

Simulated interior floods and the flood experience in urban communities

The Hydrologic Engineering Center (HEC-1) model was used to construct synthetic hydrographs for isolated interior urban floods. Flood peak and lag time were very well preserved in simulated flows. Total volume was not adequately expressed. Lag time varied inversely with both urban development and storm intensity. Peak discharge varied with storm intensity, but this variability was well defined only at very high urbanization levels. An 175% increase in storm intensity produced a change of about 15% in peak discharge. Claims for flood damage correlated well with estimates of peak flow and lag time combined. Other measures of flood experience also correlated with the two features. Within the range of storms utilized, urban development factors consistently outranked storm intensity as a determining factor in flood damage.     

Utility of Gridded Rainfall for IHACRES Daily River Flow Predictions in Southern California Watersheds1

Abstract: Studies to regionalize conceptual hydrologic models generally require rainfall and river flow data from multiple watersheds. Besides the considerable time (cost) to assemble and process rainfall data for many watersheds, investigators often need to choose from a number of candidate gauges, subjectively weighing the relative importance of proximity and elevation to select a representative rainfall dataset. The Unified Raingauge Dataset (URD) is a gridded daily rainfall dataset that covers the conterminous United States at 0.25 × 0.25 degrees spatial resolution and is available from 1948 to present. The objective of this study was to determine whether uncertainty in daily river flow predictions using the conceptual hydrologic model IHACRES in small to moderate size watersheds (50‐400 km²) in southern California would increase if URD gridded rainfall data were used in place of single rain gauge data to calibrate the model. Rain gauge data were obtained from the gauge nearest the watershed centroid and the gauge closest in elevation to the watershed mean elevation. Results from 20 randomly selected watersheds indicated that IHACRES calibration performance was similar using rainfall data from the URD grids and rain gauge data. There was some evidence of greater uncertainties associated with the URD calibrations in areas where topography may affect rainfall amounts. In contrast to the URD data, monthly gridded data produced by the Parameter‐Elevation Regressions on Independent Slopes Model (PRISM) includes adjustments for elevation and produces gridded values at a finer spatial resolution (4 km²). A limited test on two watersheds demonstrated that scaling the URD daily rainfall estimates to match the PRISM monthly values may improve rainfall estimates and model simulation performance.     

A Power Law Model for River Flow Velocity in Iowa Basins

This study explores power law relationships to estimate water flow velocity as a function of discharge and drainage area across river networks. We test the model using empirical data from 214 United States (U.S.) Geological Survey gauging stations distributed over the state of Iowa in the U.S. The empirical data are the measurements of the mean cross‐sectional velocity and concurrent discharge. The data are used to estimate parameters for a state‐wide model and to test for spatial variability for 15 large river basins contained within the state. Spatial differences among the basins are small but some parameters significantly differ from the state‐wide model. Using individual station data, the authors also explore a simpler power law model that disregards dependence on the drainage area. Overall, the study shows that including drainage area improves the model. Our study provides parameter values that can be directly incorporated into a regional scale routing model, and provides a framework for developing flow velocity models for hydraulically similar rivers in the U.S. and the world.     

BASE FLOW RECESSION RATES,LOW FLOWS,AND HYDROLOGIC FEATURES OF SMALL WATERSHEDS IN PENNSYLVANIA,USA1

This paper examines the relationships between measurable watershed hydrologic features, base flow recession rates, and the Q_7,10 low flow statistic (the annual minimum seven‐day average streamflow occurring once every 10 years on average). Base flow recession constants were determined by analyzing hydrograph recession data from 24 small (>130 km²), unregulated watersheds across five major physiographic provinces of Pennsylvania, providing a highly variable dataset. Geomorphic, hydrogeologic, and land use parameters were determined for each watershed. The base flow recession constant was found to be most strongly correlated to drainage density, geologic index, and ruggedness number (watershed slope); however, these three parameters are intercorrelated. Multiple regression models were developed for predicting the recession rate, and it was found that only two parameters, drainage density and hydrologic soil group, were required to obtain good estimates of the recession constant. Equations were also developed to relate the recession rates to Q_7,10 per unit area, and to the Q_7,10/Q₅₀ ratio. Using these equations, estimates of base flow recession rates, Q_7,10, and streamflow reduction under drought conditions can be made for small, ungaged basins across a wide range of physiography.     

AN EVALUATION OF FLOOD FREQUENCY ESTIMATES BASED ON RAINFALL/RUNOFF MODELING1

ABSTRACT: An evaluation of flood frequency estimates simulated from a rainfall/runoff model is based on (1) computation of the equivalent years of record for regional estimating equations based on 50 small stream sites in Oklahoma and (2) computation of the bias for synthetic flood estimates as compared to observed estimates at 97 small stream sites with at least 20 years of record in eight eastern states. Because of the high intercorrelation of synthetic flood estimates between watersheds, little or no regional (spatial) information may be added to the network as a result of the modeling activity. The equivalent years of record for the regional estimating equations based totally on synthetic flood discharges is shown to be considerably less than the length of rainfall record used to simulate the runoff. Furthermore, the flood estimates from the rainfall/runoff model consistently underestimate the flood discharges based on observed record, particularly for the larger floods. Depending on the way bias is computed, the synthetic estimate of the 100-year flood discharge varies from 11 to 29 percent less than the value based on observed record. In addition, the correlation between observed and synthetic flood frequency estimates at the same site is also investigated. The degree of correlation between these estimates appears to vary with recurrence interval. Unless the correlation between these two estimates is known, it is not possible to compute a weighted estimate with minimum variance.     

TOTAL SOLUTES LOADS AND GEOCHEMICAL YIELDS IN THE SASKATCHEWAN RIVER BASIN1

ABSTRACT: River solute loads have seldom been measured in very large, complex drainage basins, nor have the methods of calculating loads been critically examined. For sites in the Saskatchewan River Basin, Canada, rating curves were poor predictors of solute loads because correlations between discharge and total solutes concentration were weak (R² < 0.05 in most cases) and suffered from hysteresis. In contrast, the interval method produced reliable estimates in all seasons and sites tested, and was little affected by sampling schedule. The limit of precision (SE) for estimates of mean annual or seasonal solute load was 10–15 percent of the mean (5 percent in very small basins), reached with 10 years or more of data. Two-thirds or more of total annual solute load was transported during the open-water season, but the proportion carried during winter increased from 8 percent to 34 percent from the upstream to the downstream end of the basin, due to reservoirs retaining and mixing water. Annual loads of total solutes varied from 6.2 × 10⁴ tonnes in foothills tributaries to almost 4.0 × 10⁶ tonnes in the Saskatchewan River near the mouth. But, on an areal basis, the mountain and foothills region was the dominant solute source, producing 43–97 tonnes/km²/yr, compared with only 3–22 tonnes/km²/yr for prairie rivers. This difference is a consequence of greater rainfall and, hence, more rapid erosion in the mountains.     

CAUSES OF PEAK FLOWS IN NORTHWESTERN MONTANA AND NORTHEASTERN IDAHO1

ABSTRACT: Both catchment experiments and a review of hydrologic processes suggest a varying effect of forest harvest on the magnitude of peak flows according to the cause of those peak flows. In northwestern Montana and Northeastern Idaho, annual maximum flows can result from spring snowmelt, rain, mid-winter rain-on-snow, or rain-on-spring-snowmelt. Meteorologic and physical data were used to determine the cause of annual maximum flows in six basins which had the necessary data and were smaller than 150 mi². Rain-on-spring-snowmelt was the most frequent cause of annual maximum flows in all six basins, although there was a strong gradient in the magnitude and cause of peak flows from southwest to northeast. Less frequent mid-winter rain-on-snow events caused the largest flows on record in four basins. Mid-winter rain-on-snow should be distinguished from rain-on-spring-snowmelt because of differences in seasonal timing, the relative contributions of rain vs. snowmelt, and the projected effects of forest harvest. The effects of mixed flood populations on the flood-frequency curve varied from basin to basin. Annual maximum daily flows could not be reliably predicted from rainfall and snowmelt data.     

Uncertainty Models for Estimates of Physical Characteristics of River Segments Over Large Areas

Spatially comprehensive estimates of the physical characteristics of river segments over large areas are required in many large‐scale analyses of river systems and for the management of multiple basins. Remote sensing and modeling are often used to estimate river characteristics over large areas, but the uncertainties associated with these estimates and their dependence on the physical characteristics of the segments and their catchments are seldom quantified. Using test data with varying degrees of independence, we derived analytical models of the uncertainty associated with estimates of upstream catchment area (CA), segment slope, and mean annual discharge for all river segments of a digital representation of the hydrographic network of France. Although there were strong relationships between our test data and estimates at the scale of France, there were also large relative local uncertainties, which varied with the physical characteristics of the segments and their catchments. Discharge and CA were relatively uncertain where discharge was low and catchments were small. Discharge uncertainty also increased in catchments with large rainfall events and low minimum temperature. The uncertainty of segment slope was strongly related to segment length. Our uncertainty models were consistent across large regions of France, suggesting some degree of generality. Their analytical formulation should facilitate their use in large‐scale ecological studies and simulation models.     

Predicting RUSLE (Revised Universal Soil Loss Equation) Monthly Erosivity Index from Readily Available Rainfall Data in Mediterranean Area

Summary Seasonal rainerosivity is important in the structure and dynamics of Mediterranean ecosystems. The present paper contributes to the quantitative assessment of RUSLE's monthly erosion index in a data-scarce Mediterranean region. Therefore, a regionalized relationship for estimating monthly erosion index (EI_30-month) from only three rainfall parameters has been obtained. Knowledge of the seasonal and annual distribution of erosivity index, permit soil and water conservationists to make improved designs for erosion control, water harvesting or small hydraulic structures. Although a few long data sets were used in the analysis, validation with established monthly erosivity index values from other Italian locations, suggest that the model presented (r² = 0.973) is robust. It is recommended to monthly erosivity estimates when experimental data-scarce rainfall become available.     

APPLICATION OF GREY MODEL AND ARTIFICIAL NEURAL NETWORKS TO FLOOD FORECASTING1

The main focus of this study was to compare the Grey model and several artificial neural network (ANN) models for real time flood forecasting, including a comparison of the models for various lead times (ranging from one to six hours). For hydrological applications, the Grey model has the advantage that it can easily be used in forecasting without assuming that forecast storm events exhibit the same stochastic characteristics as the storm events themselves. The major advantage of an ANN in rainfall‐runoff modeling is that there is no requirement for any prior assumptions regarding the processes involved. The Grey model and three ANN models were applied to a 2,509 km² watershed in the Republic of Korea to compare the results for real time flood forecasting with from one to six hours of lead time. The fifth‐order Grey model and the ANN models with the optimal network architectures, represented by ANN1004 (34 input nodes, 21 hidden nodes, and 1 output node), ANN1010 (40 input nodes, 25 hidden nodes, and 1 output node), and ANN1004T (14 input nodes, 21 hidden nodes, and 1 output node), were adopted to evaluate the effects of time lags and differences between area mean and point rainfall. The Grey model and the ANN models, which provided reliable forecasts with one to six hours of lead time, were calibrated and their datasets validated. The results showed that the Grey model and the ANN1010 model achieved the highest level of performance in forecasting runoff for one to six lead hours. The ANN model architectures (ANN1004 and ANN1010) that used point rainfall data performed better than the model that used mean rainfall data (ANN1004T) in the real time forecasting. The selected models thus appear to be a useful tool for flood forecasting in Korea.     

RESPONSE OF GREEN‐AMPT MEIN‐LARSEN SIMULATED RUNOFF VOLUMES TO TEMPORALLY AGGREGATED PRECIPITATION1

ABSTRACT: Data collection frequency in automated systems is user determined and can range from seconds to hours or days. Currently, there is no standard or recommended frequency interval for collecting precipitation data from automated systems for input to event‐based models such as Green‐Ampt Mein‐Larsen (GAML). Data from 47 storm events at seven locations were used to simulate the response of GAML excess rainfall to temporally aggregated precipitation data. No difference in model efficiency was recognized when comparing one‐minute interval data (R²= 1.00) to five‐minute data (R²= 1.00). Very little model efficiency was lost at a 10‐minute (R²= 0.96) interval. After 10‐minutes, decline in efficiency became more rapid with R²= 0.16 at one hour. The combined effect of time interval with respect to drainage area, hydraulic conductivity, maximum 30‐minute intensity, and total precipitation also revealed similar results.     

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.     

MODIFICATION OF NEILL'S EQUATION FOR STORM‐TIDE DESIGN FLOW ANALYSIS1

ABSTRACT: Design of bridges spanning tidal estuaries or bays requires an estimate of peak tidal flow. One common approach to estimating these flows (Neill's method) uses a first‐order approximation of uniform water surface rise in the water body. For larger water bodies, the assumptions of this method are decreasingly valid. This study develops a simple modification that accounts for the spatial variability in the response of tidal waterways to storm surge flows. The peak tidal flow predicted by Neill's equation is compared to the peak flow determined by numerical simulation of estuaries with simple geometries, ranging from 1 to 25 km in length, using the U.S. Army Corps of Engineers one‐dimensional unsteady flow model, UNET. Results indicate that, under certain conditions, it may be appropriate to apply a correction factor to the peak discharge and peak velocity predicted by Neill's method. An algorithm, developed by nonlinear regression, is presented for computing correction factors based on estuary length, shape, mean depth, and storm‐tide characteristics. The results should permit the design of more reliable, cost‐effective structures by providing more realistic estimates of the potential for bridge scour in tidal waterways, especially when a full solution of the unsteady flow equations is impractical.     

Developing Peak Discharges for Future Flood Risk Studies Using IPCC's CMIP5 Climate Model Results and USGS WREG Program

Extreme climate events, floods, and drought, cause huge impact on daily lives. In order to produce society resilient to extreme events, it is necessary to assess the impact of frequent and high intensity storm events on design parameters. This article describes a methodology to develop future peak “design discharges” throughout the United States that can be used as a guidance to map future floodplains. In order to develop a lower and upper limit for anticipated peak flow discharges, two future growth scenarios — Representative Concentration Pathways (RCPs)‐RCP 2.6 and 8.5 were identified as the weak and strong climate scenario respectively based on the output from the global climate models. The Generalized Least Square technique in United States Geological Survey's Weighted Multiple Regression (WREG) program was used to develop regression equations that relate peak discharges to basin and climate parameters of the contributing watershed. The design discharges reflect the most recent climate model results. Number of frost days, heavy rainfall days, high temperature days, and snow depth were found to be the common extreme climate parameters influencing the regression equations. This methodology can be extended to other flood frequency events if rainfall data is available. The future discharges can be utilized in hydraulics models to estimate floodplains that can assist in resilient infrastructure planning and outline climate change adaptation strategies.     

RAINFALL POINT INTENSITIES IN AN AIR MASS THUNDERSTORM ENVIRONMENT: WALNUT GULCH,ARIZONA1

ABSTRACT: Point rainfall intensities for a given return period are often used to formulate design storms for rainfall/runoff models to simulate design floods. These design floods are in turn used to design bridges, culverts, and a variety of drainage and flood control structures. The projected rapid growth in the southwestern United States will require very substantial monetary investments in drainage infrastructure. Accurate estimates of point rainfall intensities are critical to ensure both safe designs while not wasting dollars in overdesign. Rainfall point intensities (accumulated rainfall depth over a specified duration) for 5‐, 15‐, 30‐, and 60‐minute durations for the 2‐, 5‐, 10‐, 25‐, 50‐, and 100‐year return periods were determined for southeast Arizona. Thirty‐five years of rainfall record (1961 to 1995) were used in this study. The records came from 20 stations that were grouped into five sets of four independent stations to extend the rainfall records. The stations are in the USDA‐ARS Walnut Gulch Experimental Watershed (WGEW), which is representative of large portions of the Southwest whose runoff generation is dominated by air‐mass thunderstorms. The 5‐, 15‐, 30‐, and 60‐minute maximum intensities per year followed log‐normal distributions. The mean point rainfall intensities of the five sets of gages are very close (between 0 and 11 percent) to the NOAA values of the 5‐, 15‐, 30‐, and 60‐minute durations for all return periods. Much larger differences between the mean point rainfall intensities for all durations were found when these results were compared to those of a previous study done with a shorter rainfall record (between 14 and 33 percent for the 25‐, 50‐, and 100‐year return‐periods). The difference between the largest and the smallest values of point rainfall intensities recorded by each group, for all durations, usually increases as the return period increases.     

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.