COMPARISON OF METHODS FOR ESTIMATING FLOOD MAGNITUDES ON SMALL STREAMS IN GEORGIA1

ABSTRACT: The U.S. Geological Survey has collected flood data for small, natural streams at many sites throughout Georgia during the past 20 years. Flood-frequency relations were developed for these data using four methods: (1) observed (log-Pearson Type HI analysis) data, (2) rainfall-runoff model, (3) regional regression equations, and (4) map-model combination. The results of the latter three methods were compared to the analyses of the observed data in order to quantify the differences in the methods and determine if the differences are statistically significant. Comparison of regression-estimates with observed discharges for sites having 20 years (1966 to 1985) and 10 years (1976 to 1985) of record at different sites of annual peak record indicate that the regression-estimates are not significantly different from the observed data. Comparison of rainfall-runoff-model simulated estimates with observed discharges for sites having 10 years and 20 years of annual peak record indicated that the model-simulated estimates are significantly and not significantly different, respectively. The biasedness probably is due to a “loss of variance” in the averaging procedures used within the model and the short length of record as indicated in the 10 and 20 years of record. The comparison of map-model simulated estimates with observed discharges for sites having 20 years of annual-peak runoff indicate that the simulated estimates are not significantly different. Comparison of “improved” map-model simulated estimates with observed discharges for sites having 20 years of annual-peak runoff data indicate that the simulated estimates are different. The average adjustment factor suggested by Lichty and Liscum to calculate the “improved” map-model overestimates in Georgia by an average of 20 percent for three recurrence intervals analyzed.     

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).     

VERIFICATION AND APPLICATION OF REGIONAL EQUATIONS FOR THE STORAGE FUNCTION RUNOFF MODEL1

The storage function model is a nonlinear rainfall-runoff model that has been developed for and applied to flood runoff analysis in Japan. This paper extends the model applicability by developing practical equations for estimating model parameters which are appropriate on a regional basis, i.e., so-called regional equations. Previously, the parameters were computed from historical data for a specific basin or from relationships that do not account for land use and topography. To develop the regionalized equations, model parameters were identified for 91 flood events from 22 watersheds in Japan by applying a mathematical optimization technique. Results from 39 of these events were statistically compared and regional relationships were determined as a function of land use, basin area and rainfall intensity. The utility of the estimated equations were tested by computing runoff hydrographs for lumped basins. The estimated parameters were also applied in a distributed watershed model formulation. Both applications showed acceptable results that validate the use of the regionalized relationships.     

THE INFLUENCE OF CLIMATE VARIATION ON THE ESTIMATION OF LOW FLOWS USED TO PROTECT WATER QUALITY: A NATIONWIDE ASSESSMENT1

ABSTRACT: Historical flow records are used to estimate the regulatory low flows that serve a key function in setting discharge permit limits through the National Pollutant Discharge Elimination System, which provides a nationwide mechanism for protecting water quality. Use of historical records creates an implicit connection between water quality protection and climate variability. The longer the record, the more likely the low flow estimate will be based on a broad set of climate conditions, and thus provides adequate water quality protection in the future. Unfortunately, a long record often is not available at a specific location. This analysis examines the connection between climate variability and the variability of biologically based and hydrologically based low flow estimates at 176 sites from the Hydro‐Climatic Data Network, a collection of stream gages identified by the USGS as relatively free of anthropogenic influences. Results show that a record of 10 to 20 years is necessary for satisfactory estimates of regulatory low flows. Although it is possible to estimate a biologically based low flow from a record of less than 10 years, these estimates are highly uncertain and incorporate a bias that undermines water quality protection.     

FLOOD EVENT MODELING - A STUDY OF TWO METHODS1

: Despite the advances in catchment modeling in recent years, engineers still face major problems in estimating flood flows. Application of unit hydrograph and runoff routing models to five United Kingdom catchments shows that either can be tuned to predict, on a test event, the routing effects of a catchment with equal accuracy. The larger remaining problem is the prediction of losses from rainfall and this study shows how alternative ways of describing the within event distribution of these losses can be an important factor controlling the success of the overall model. Other problems include the risks of extrapolation to larger events, baseflow separation methods, rainfall patterns, and inevitable errors in the data.     

Application of Index Procedures to Flood Frequency Analysis in Turkey1

Abstract: This study investigates the regional analysis of annual maximum flood series of 48 stream gauging stations in the basins of the West Mediterranean Region in Turkey. The region is divided into three homogeneous subregions according to both Student‐t test and Dalrymple homogeneity test. The regional relationships of mean annual flood per unit area‐drainage area and coefficient of skew‐coefficient of variation are obtained. Two statistically meaningful relationships of the mean flood per unit area‐drainage area and a unique relationship between skewness and variation coefficients exist. Results show that the index‐flood method may be applicable to each homogenous subregion to estimate flood quantiles in the study area.     

Nonparametric Monte Carlo Simulation for Flood Frequency Curve Derivation: An Application to a Korean Watershed1

Abstract: A mix of causative mechanisms may be responsible for flood at a site. Floods may be caused because of extreme rainfall or rain on other rainfall events. The statistical attributes of these events differ according to the watershed characteristics and the causes. Traditional methods of flood frequency analysis are only adequate for specific situations. Also, to address the uncertainty of flood frequency estimates for hydraulic structures, a series of probabilistic analyses of rainfall‐runoff and flow routing models, and their associated inputs, are used. This is a complex problem in that the probability distributions of multiple independent and derived random variables need to be estimated to evaluate the probability of floods. Therefore, the objectives of this study were to develop a flood frequency curve derivation method driven by multiple random variables and to develop a tool that can consider the uncertainties of design floods. This study focuses on developing a flood frequency curve based on nonparametric statistical methods for the estimation of probabilities of rare floods that are more appropriate in Korea. To derive the frequency curve, rainfall generation using the nonparametric kernel density estimation approach is proposed. Many flood events are simulated by nonparametric Monte Carlo simulations coupled with the center Latin hypercube sampling method to estimate the associated uncertainty. This study applies the methods described to a Korean watershed. The results provide higher physical appropriateness and reasonable estimates of design flood.     

COMPARISON OF ANNUAL AND PARTIAL DURATION SERIES FLOODS ON THE MURRUMBIDGEE RIVER1

ABSTRACT: In current hydrologic practice flood frequency estimates are usually based upon either the annual or the partial duration series of floods. Recurrence intervals generated by each series are not equivalent, however, and conversion of recurrence intervals from one series to the other is usually achieved by reference to a mathematical function developed by Langbein in 1949. Data collected on the Murrumbidgee River in New South Wales suggest, however, that the Langbein conversion function does not always provide a reliable means of comparing recurrence intervals. For discharges more frequent than the three year annual flood the Langbein function understates the discrepancy between the two sets of recurrence interval by approximately 35 percent. Langbein's own North American data appear to be consistent with those collected on the Murrumbidgee River.     

THE PROBABILITY OF CONSECUTIVE RAINLESS DAYS1

ABSTRACT: A major objective of this work was to develop a method applicable to the Intermountain region for estimating the probability of n consecutive dry days, where n ≤ 30 days. One result was a computationally simple method of producing the desired estimates directly from the rainfall record. For a consecutive dry-day period of n days, these estimates are equivalent to those obtained from a Markov model of order n-1. The efficacy of a simple Markov model was also evaluated. In this climatic region, the simple Markov model produces probability estimates of consecutive dry days that are too conservative, especially at long dry-day periods. In this data set, it was found that the longer the dry-day sequence, the more conservative the Markov estimate. The source of these conservative estimates is the strong dry-day persistence, which is characteristic of summer weather in the Intermountain region. In this region, the best estimates of the probability of consecutive dry days are probably those obtained directly from a representative rainfall record and smoothed by the partial sums of a fourier Series.     

FLOOD-RUNOFF FORECASTING WITH HEC1F1

HEC1F is a computer program for making short- to medium-term forecasts of uncontrolled flood runoff. The program employs unit hydrographs and hydrologic routing to simulate runoff from a subdivided basin. Estimates of future rainfall can be accommodated. Runoff parameters for gaged headwater subbasins can be estimated (optimized) in real time. Blending of calculated with observed hydrographs can be performed. HEC1F is a component of an on-line software system that includes capability for data acquisition and processing, precipitation analysis, streamflow forecasting, reservoir system analysis, and graphical display of data and simulation results. The conceptual framework for HEC1F is described, and application of the program is illustrated.     

PARAMETER UNCERTAINTY IN RAINFALL-RUNOFF MODELING FOR POLDER SYSTEM DESIGN1

ABSTRACT: An approach is developed for incorporating the uncertainty of parameters for estimating runoff in the design of polder systems in ungaged watersheds. Monte Carlo Simulation is used to derive a set of realizations of streamflow hydrographs for a given design rainstorm using the U. S. Soil Conservation Service (SCS) unit hydrograph model. The inverse of the SCS curve number, which is a function of the antecedent runoff condition in the SCS model, is the random input in the Monte Carlo Simulation. Monte Carlo realizations of streamfiow hydrographs are used to simulate the performance of a polder flood protection system. From this simulation the probability of occurrence of flood levels for a particular hydraulic design may be used to evaluate its effectiveness. This approach is demonstrated for the Pluit Polder flood protection system for the City of Jakarta, Indonesia. While the results of the application indicate that uncertainty in the antecedent runoff condition is important, the effects of uncertainty in rainfall data, in additional runoff parameters, such as time to peak, in the hydraulic design, and in the rainfall-runoff model selected should also be considered. Although, the SCS model is limited to agricultural conditions, the approach presented herein may be applied to other flood control systems if appropriate storm runoff models are selected.     

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.     

COMPUTATION OF STAGE-DISCHARGE RELATIONSHIPS AFFECTED BY UNSTEADY FLOW1

ABSTRACT: The dynamic relationship between stage and discharge which is unique to a particular flood for a selected station along the river can be determined via a mathematical model based on the complete one-dimensional equations of unsteady flow, i.e., the equations for the conservation of mass and momentum of the flood wave, and the Manning equation which accounts for energy losses. By assuming the bulk of the flood wave moves as a kinematic wave, the need for spatial resolution of the flood can be eliminated, and only the time variation of either the discharge or stage at the selected station is necessary for the computation of the other. The mathematical model can be used in river forecasting to convert the forecast discharge hydrograph into a stage hydrograph which properly reflects the unique dynamic stage-discharge relationship produced by the variable energy slope of the flood discharge. The model can be used also in stream gaging to convert a recorded stage hydrograph into a discharge hydrograph which properly accounts for the effects of unsteady flow. The model is applied to several observed floods at selected stations along the Lower Mississippi, Red, and Atchafalaya Rivers. The root mean square errors between observed and computed discharges are in the range of 3 to 7 percent, values well within the accuracy of the observations. A simple, easily-applied graphical procedure is also provided for estimating the magnitude of the effect of the unsteady flow on stage-discharge ratings. As a general rule, the dynamic effect may be significant if the channel bottom slope is less than 0.001 ft/ft (about 5 ft/mi) when the rate of change of stage is greater than about 0.10 ft/hr.     

ASSESSING FLOOD MITIGATION ALTERNATIVES IN SHIJR AREA IN METROPOLITAN TAIPEI1

The Keelung River Basin in northern Taiwan lies immediately upstream of the Taipei metropolitan area. The Shijr area is in the lower basin and is subject to frequent flooding. This work applies micromanagement and source control, including widely distributed infiltration and detention/ retention runoff retarding measures, in the Wudu watershed above Shijr. A method is also developed that combines a genetic algorithm and a rainfall runoff model to optimize the spatial distribution of runoff retarding facilities. Downstream of Wudu in the Shijr area, five dredging schemes are considered. If 10‐year flood flows cannot be confined in the channel, then a levee embankment that corresponds to the respective runoff retarding scheme will be required. The minimum total cost is considered in the rule to select from the regional flood mitigation alternatives. The results of this study reveal that runoff retarding facilities installed in the upper and middle parts of the watershed are most effective in reducing the flood peak. Moreover, as the cost of acquiring land for the levee embankment increases, installing runoff retarding measures in the upper portion of the watershed becomes more economical.     

APPLICATION OF AN INTEGRATED BASIN-SCALE HYDROLOGIC MODEL TO SIMULATE SURFACE-WATER AND GROUND-WATER INTERACTIONS1

ABSTRACT: Hydrologic models have become an indispensable tool for studying processes and water management in watersheds. A physically-based, distributed-parameter model, Basin-Scale Hydro-logic Model (BSIIM), has been developed to simulate the hydrologic response of large drainage basins. The model formulation is based on equations describing water movement both on the surface and in the subsurface. The model incorporates detailed information on climate, digital elevation, and soil moisture budget, as well as surface-water and ground-water systems. This model has been applied to the Big Darby Creek Watershed, Ohio in a 28-year simulation of rainfall-runoff processes. Unknown coefficients for controlling runoff, storativity, hydraulic conductivity, and streambed permeability are determined by a trial-and-error calibration. The performance of model calibration and predictive capability of the model was evaluated based on the correlation between simulated and observed daily stream discharges. Discrepancies between observed and simulated results exist because of limited precipitation data and simplifying assumptions related to soil, land use, and geology.     

GIS BASED LONG TERM HYDROLOGIC IMPACT EVALUATION FOR WATERSHED URBANIZATION1

ABSTRACT: To adequately manage impacts of ongoing or future land use changes in a watershed, the magnitude of their hydrologic impacts needs to be assessed. A grid based daily streamflow model was calibrated with two years of observed streamflow data, using time periods when land use data are available and verified by comparison of model predictions with observed streamflow data. Streamflow data were separated into direct runoff and baseflow to estimate the impacts of urbanization on each hydrologic component. Analysis of the ratio between direct runoff and total runoff from 30 years of simulation results and the change in these ratios with urbanization shows that estimated annual direct runoff increased from 49.2 percent (1973) to 63.1 percent (1984) and 65.0 percent (1991), indicating the effects of urbanization are greater on direct runoff than on total runoff. The direct runoff ratio also varies with annual rainfall, with dry year ratios larger than those for wet years. This suggests that the impact of urbanization on areas that are sensitive to runoff ratios, such as stream ecosystems, might be more serious during drier years than in wetter years in terms of water quality and water yield. This indicates that sustainable base‐flow is important to maintaining sound stream ecosystems.     

A Framework to Develop Nationwide Flooding Extents Using Climate Models and Assess Forecast Potential for Flood Resilience

The methods used to simulate flood inundation extents can be significantly improved by high‐resolution spatial data captured over a large area. This paper presents a hydraulic analysis methodology and framework to estimate national‐level floodplain changes likely to be generated by climate change. The hydraulic analysis was performed using existing published Federal Emergency Management Agency 100‐year floodplains and estimated 100‐ and 10‐year return period peak flow discharges. The discharges were estimated using climate variables from global climate models for two future growth scenarios: Representative Concentration Pathways 2.6 and 8.5. River channel dimensions were developed based on existing regional United States Geological Survey publications relating bankfull discharges with channel characteristics. Mathematic relationships for channel bankfull topwidth, depth, and side slope to contributing drainage area measured at model cross sections were developed. The proposed framework can be utilized at a national level to identify critical areas for flood risk assessment. Existing hydraulic models at these “hot spots” could be repurposed for near–real‐time flood forecasting operations. Revitalizing these models for use in simulating flood scenarios in near–real time through the use of meteorological forecasts could provide useful information for first responders of flood emergencies.     

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.     

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.     

FLOOD PLAIN MODELING1

Computer simulation of river basin hydrology has rapidly progressed from an interesting academic exercise to a practical engineering procedure of increasing utility. Mathematical models of the many interrelated processes which occur in a basin have been developed along with efficient numerical procedures for their solution. The present paper is concerned with a particular model which has been used to describe the transformation of a temporally and spatially varying rainfall into a time history of stage and discharge on a flood plain. Although developed principally as a model of the physical processes involved, it is envisioned that the model can serve as one component of an information system for flood plain planning and management. The model consists of the following elements: (i) a rainfall simulation which generates stochastic inputs to the model according to specified rainfall statistics, (ii) a catchment-runoff model which converts the rainfall to surface runoff, (iii) a flood stage model which converts the surface runoff to time histories of flood stage and discharge. The model has been used to investigate the effect of various structural flood control measures in a basin and, in particular, to establish frequency-stage information for each of these. Of particular interest in development of the model have been recurring and partially unanswered questions relative to the proper balance among availability of data for use in the model, data requirements of the model and the objectives of the outputs produced by the model.