FLOOD ROUTING THROUGH A FLAT,COMPLEX FLOODPLAIN USING A ONE-DIMENSIONAL UNSTEADY FLOW COMPUTER PROGRAM1

ABSTRACT: The routing of flood waves through the Central Basin of the Passaic River in New Jersey is complex because of flat gradients and flow reversals. The one-dimensional unsteady flow program DWOPER, developed by the National Weather Service, was used to simulate flood wave movement through the Basin. A historical event was used for calibration and two synthetic events were simulated. Boundary conditions consisted of discharge hydrographs at inflow points to the study area, local flow hydrographs at interior points, and a stage discharge relation for flow over the crest of a diversion dam at the basin outlet. Manning's n values were adjusted based on stage and discharge data for the historical event; however, verification data were not available for events comparable in magnitude to the synthetic events. Aspects of the investigation reported include techniques for characterizing the flow system, model calibration, techniques for representing a tunnel diversion, and simulation results.     

A DISTRIBUTED DYNAMIC WATERSHED MODEL1

ABSTRACT: A distributed watershed model was developed to mathematically simulate overland and channel flow for a single-event storm. The modeled watersheds in the study were subdivided into rectangular grid elements. All hydrologically significant parameters, such as land slope, rainfall and precipitation excess, were assumed to be uniform within each element. The Green-Ampt method was adopted to generate precipitation excess for each element during the simulation period. A two-dimensional diffusion wave model was used for overland flow routing and an iterative Alternative Direction Implicit scheme was used to solve the simultaneous overland flow equations. Once the overland flow became inflow to the channel, a one-dimensional dynamic wave flood routing technique, based on a four-point, implicit, non-linear finite difference solution of the St. Venant equation of unsteady flow, was applied. A limited number of comparisons were made between simulated and observed hydrographs for areas of about one square mile. Given the appropriate parameters, the model was able to accurately simulate runoff for single-event storms. This paper describes a distributed watershed model developed to simulate overland and channel flow. Comparisons were made between simulated and observed hydrographs for three watersheds. The model was able to accurately simulate the runoff for single-event storms using 61-m by 61-m (200 ft by 200 ft) watershed grid elements.     

Comparison of Unsteady and Quasi‐Unsteady Flow Models in Simulating Sediment Transport in an Ephemeral Arizona Stream1

Hummel, Ryan, Jennifer G. Duan, and Shiyan Zhang, 2012. Comparison of Unsteady and Quasi‐Unsteady Flow Models in Simulating Sediment Transport in an Ephemeral Arizona Stream. Journal of the American Water Resources Association (JAWRA) 48(5): 987‐998. DOI: 10.1111/j.1752‐1688.2012.00663.x Abstract: Hydrodynamic and sediment transport models are useful engineering tools for predicting unsteady flood flow and sediment transport. Many models such as HEC‐RAS, HEC‐6, and IALLUVIAL apply quasi‐unsteady flow model, whereas others apply the unsteady flow model. It remains unknown if a quasi‐unsteady flow model is sufficiently accurate for simulating sediment transport in rapidly varied unsteady flood events, especially in ephemeral rivers in arid and semiarid regions. This study compared the quasi‐unsteady HEC‐RAS 4.1 model with one‐dimensional (1D) Finite Volume Method (FVM) based model in simulating flood flow and sediment transport in the Pantano Wash, a dryland river in the state of Arizona. The objective is to determine which sediment transport method is appropriate in predicting bed elevation changes in an ephemeral stream, Pantano Wash, and if an unsteady model is more accurate than a quasi‐unsteady flow model in predicting sediment transport. Results showed that the quasi‐unsteady HEC‐RAS model and the 1D FVM yielded similar results of bed degradation and aggradation for this dryland stream, although the FVM model predicted better flood hydrographs. Among the seven sediment transport formulas embedded in HEC‐RAS, Yang’s and Engelund‐Hansen’s equations gave the best matches with the field measurements for this particular case study.     

去学水电站生态流量泄放措施优化研究

水电站的建设会造成坝址下游河段减水,对水环境和水生生物造成不利影响,因此电站生态流量的确定和泄放措施的可行性对于维护河流生态环境至关重要。以硕曲河去学水电站为例,基于可行性研究阶段提出的生态流量泄放措施方案,采用数学模型和理论分析方法,从生态流量取水水温、泄洪洞事故条件下生态流量保证率以及水资源的合理高效利用等方面进一步对生态流量泄放方案进行了优化比选。提出的优化方案提高了生态流量取水口的布置高程,在泄洪洞之外单独布置生态泄放系统,并增设生态流量发电机组等。对比分析表明,该优化方案在改善生态流量的下泄水温、提高生态流量泄放保证率、合理高效利用水能资源等方面具有较大优势。本研究为水电工程中生态流量泄放措施保障等方面提供了科学依据和技术参考。     

A GIS‐ASSISTED DISTRIBUTED WATERSHED MODEL FOR SIMULATING FLOODING AND INUNDATION1

ABSTRACT: A distributed watershed model combining kinematic wave routing, 1‐D dynamic channel‐flow routing, and 2‐D diffusive overland‐flow routing has been developed to simulate flooding and inundation levels of large watersheds. The study watershed was linked to a GIS database and was divided into an upstream mountainous area and a downstream alluvial plain. A kinematic wave routing was adopted at the mountainous area to compute the discharge flowing into the alluvial plain. A 1‐D dynamic channel routing solving the St. Venant equations by the Preissmann method was performed for the main channel of the alluvial plain, whereas a 2‐D overland‐flow routing solving the diffusion wave equation with the Alternating Direction Explicit scheme was used for floodplains. The above two routings were connected by weir‐link discharge formula. The parameters in the model were calibrated and independently verified by single‐event storms. An example application of flooding/inundation analysis was conducted for the Taichung station and the Woozi depot (Taiwan High Speed Rail). Suggested inundation‐proofing measures ‐ including raising ground surface elevation of the station and depot and building a waterproofing exterior wall and their combination ‐ were investigated. It was concluded that building the waterproofing exterior wall had a strong tendency to decrease peak inundation depth.     

Diffusive wave models for operational forecasting of channel routing at continental scale

Operational forecast models require robust, computationally efficient, and reliable algorithms. We desire accurate forecasts within the limits of the uncertainties in channel geometry and roughness because the output from these algorithms leads to flood warnings and a variety of water management decisions. The current operational Water Model uses the Muskingum-Cunge method, which does not account for key hydraulic conditions such as flow hysteresis and backwater effects, limiting its ability in situations with pronounced backwater effects. This situation most commonly occurs in low-gradient rivers, near confluences and channel constrictions, coastal regions where the combined actions of tides, storm surges, and wind can cause adverse flow. These situations necessitate a more rigorous flow routing approach such as dynamic or diffusive wave approximation to simulate flow hydraulics accurately. Avoiding the dynamic wave routing due to its extreme computational cost, this work presents two diffusive wave approaches to simulate flow routing in a complex river network. This study reports a comparison of two different diffusive wave models that both use a finite difference solution solved using an implicit Crank–Nicolson (CN) scheme with second-order accuracy in both time and space. The first model applies the CN scheme over three spatial nodes and is referred to as Crank–Nicolson over Space (CNS). The second model uses the CN scheme over three temporal nodes and is referred to as Crank–Nicolson over Time (CNT). Both models can properly account for complex cross-section geometry and variable computational points spacing along the channel length. The models were tested in different watersheds representing a mixture of steep and flat topographies. Comparing model outputs against observations of discharges and water levels indicated that the models accurately predict the peak discharge, peak water level, and flooding duration. Both models are accurate and computationally stable over a broad range of hydraulic regimes. The CNS model is dependent on the Courant criteria, making it less computational efficient where short channel segments are present. The CNT model does not suffer from that constraint and is, thus, highly computationally efficient and could be more useful for operational forecast models.     

TWO-DIMENSIONAL DIFFUSION-PROBABILISTIC MODEL OF A SLOW DAM BREAK1

ABSTRACT: A two-dimensional model of a dam-break flood wave is developed by simplifying the St. Venant equations to eliminate local acceleration and inertial terms and combining the simplified equations with continuity to form a diffusion type partial differential equation. This model is cascaded with a two point probability estimate scheme to account for uncertainty in the dam break flood hydrograph and channel roughness. The development and application of the probabilistic model is the main contribution of this paper. The approach is applied to a hypothetical dam break of Long Valley Dam on the Owens River above Bishop, California.     

Regional scale flood modeling using NEXRAD rainfall, GIS, and HEC-HMS/RAS: a case study for the San Antonio River Basin Summer 2002 storm event

This paper develops a framework for regional scale flood modeling that integrates NEXRAD Level III rainfall, GIS, and a hydrological model (HEC-HMS/RAS). The San Antonio River Basin (about 4000 square miles, 10,000 km2) in Central Texas, USA, is the domain of the study because it is a region subject to frequent occurrences of severe flash flooding. A major flood in the summer of 2002 is chosen as a case to examine the modeling framework. The model consists of a rainfall-runoff model (HEC-HMS) that converts precipitation excess to overland flow and channel runoff, as well as a hydraulic model (HEC-RAS) that models unsteady state flow through the river channel network based on the HEC-HMS-derived hydrographs. HEC-HMS is run on a 4 x 4 km grid in the domain, a resolution consistent with the resolution of NEXRAD rainfall taken from the local river authority. Watershed parameters are calibrated manually to produce a good simulation of discharge at 12 subbasins. With the calibrated discharge, HEC-RAS is capable of producing floodplain polygons that are comparable to the satellite imagery. The modeling framework presented in this study incorporates a portion of the recently developed GIS tool named Map to Map that has been created on a local scale and extends it to a regional scale. The results of this research will benefit future modeling efforts by providing a tool for hydrological forecasts of flooding on a regional scale. While designed for the San Antonio River Basin, this regional scale model may be used as a prototype for model applications in other areas of the country.     

EFFECT OF RAINFALL EXCESS CALCULATIONS ON MODELED HYDROGRAPH ACCURACY AND UNIT-HYDROGRAPH PARAMETERS1

ABSTRACT: Two methods of computing rainfall excess in the U.S. Army Corps of Engineers’flood hydrograph package (HEC-1), the Initial and Uniform method and the Exponential method, are compared to evaluate the effects on modeled hydrograph accuracy. Two computed unit-hydrograph parameters, time of concentration and storage coefficient, were also compared. Rainfall and runoff data from 209 storms in 32 gaged basins in Illinois were used to calibrate the HEC-1 model. Three hydrograph characteristics - sum of incremental flows, peak discharge, and time of peak discharge - were used to evaluate modeled hydrograph accuracy. Mean percent error for each basin and hydrograph characteristic was computed. An evaluation of the mean errors indicates that, although some bias in modeled hydrograph accuracy is evident, rainfall excess computed using either method results in a computed hydrograph accuracy that is within generally accepted limits. Application of a linear-regression model shows no significant differences in computed values of unit-hydrograph parameters.     

MANAGEMENT OF FLOOD CONTROL SUMPS AND POLLUTANT TRANSPORT1

ABSTRACT: Levee sump systems are used by many riverine communities for temporary storage of urban wet weather flows. The hydrologic performance and transport of stormwater pollutants in sump systems, however, have not been systematically studied. The objective of this paper is to present a case study to demonstrate development and application of a procedure for assessing the hydraulic performance of flood control sumps in an urban watershed. Two sumps of highly variable physical and hydraulic characteristics were selected for analysis. A hydrologic modeling package was used to estimate the flow hydrograph for each outfall as part of the flow balance for the sump. To validate these results, a water balance was used to estimate the total runoff using sump operational data. The hydrologic model calculations provide a satisfactory estimate of the total runoff and its time‐distribution to the sump. The model was then used to estimate pollutant loads to the sump and to the river. Although flow of stormwater through a sump system is regulated solely by flood‐control requirements, these sumps may function as sedimentation basins that provide purification of stormwater. A sample calculation of removals of several conventional pollutants in the target sumps using a mass balance approach is presented.     

COMPUTATION OF UNSTEADY FLOWS IN THE ALABAMA RIVER1

ABSTRACT: An application is described of the branch-network flow model, BRANCH, to the upper Alabama River system in central Alabama. The model is used to simulate one-dimensional unsteady flows and water surface elevations in approximately 60 river miles of the Alabama River system. Preliminary calibration was made using 72 hours of observed data. Simulated discharges are about 10 percent lower than observed discharges at higher discharge rates and computed flows lag observed flows by about 30 minutes.     

PREDICTING TEMPORAL AND SPATIAL FLOOD DYNAMICS USING A PRECALIBRATED MODEL1

ABSTRACT: A large storm in December 1990 allowed the evaluation of flood predictions from a hydrologic model (TOPMODEL) that had been previously calibrated on the West Fork of Walker Branch Watershed, a gauged 37.5 ha catchment near Oak Ridge, Tennessee. The model predicts both hydrograph dynamics and the spatial distribution of overland flow using an index based on topography. Maximum extent of overland flow during the storm was determined from patterns of leaf litter removal from valley bottoms. Both the flood hydrograph and the extent of overland flow were accurately predicted using model parameters obtained from a three-month period of normal flow conditions during 1983.     

Managing the cooling capacity of the upper rhine: A case study

In most industrialized countries, environmental standards exist which prescribe the maximum allowable man-made increase in water temperature of a river. Together with flowrate and weather conditions, these standards determine the rate at which waste heat may be discharged into a river at any moment. Power generating stations with variable cooling systems can adjust their heat discharge into the river in compliance with environmental standards and by doing so exert an influence on power generation capacity. In this paper, a scheme is developed that allows a chain of power stations discharging into the same river to operate their cooling systems such that the output of total electricity is maximized and water temperature standards are accomodated. The optimum balance between stations is determined through dynamic programming. From the results of a simulation model using historical data, simple decision rules for day-to-day operation are abstracted. These rules are based solely on the river flow rates at each power station.     

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.     

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.     

NONLINEAR MODELING AND PREDICTION OF A RIVER FLOW SYSTEM1

ABSTRACT: Model estimation and prediction of a river flow system are investigated using nonlinear system identification techniques. We demonstrate how the dynamics of the system, rainfall, and river flow can be modeled using NARMAX (Nonlinear Autoregressive Moving Average with eXogenuous input) models. The parameters of the model are estimated using an orthogonal least squares algorithm with intelligent structure detection. The identification of the nonlinear model is described to represent the relationship between local rainfall and river flow at Enoree station (inputs) and river flow at Whitmire (output) for a river flow system in South Carolina.     

APPLICATION OF THE U.S.G.S. DIFFUSION HYDRODYNAMIC MODEL FOR URBAN FLOODPLAIN ANAYSIS1

ABSTRACT: The two-dimensional Diffusion Hydrodynamic Model, DHM, is applied to the evaluation of floodplain depths resulting from an overflow of a leveed river. The environmental concerns of flood protection and high flow velocities can be better studied with the help of the two-dimensional DHM flow model than by use of the one-dimensional modeling techniques. In the test case, some of the predicted flood depth differences between the DHM and the one-dimensional approach (i.e., HEC-2) are found to be significant. Although the DHM generates considerable information, it is easy to use and does not require expertise beyond that required for use of the one-dimensional approaches.     

MEASURED AND PREDICTED VELOCITY AND LONGITUDINAL DISPERSION AT STEAI)Y AND UNSTEADY FLOW,COLORADO RIVER,GLEN CANYON DAM TO LAKE MEAD1

ABSTRACT: The effect of unsteadiness of dam releases on velocity and longitudinal dispersion of flow was evaluated by injecting a fluorescent dye into the Colorado River below Glen Canyon Dam and sampling for dye concentration at selected sites downstream. Measurements of a 26-kilometer reach of Glen Canyon, just below Glen Canyon Dam, were made at nearly steady dam releases of 139, 425, and 651 cubic meters per second. Measurements of a 380-kilometer reach of Grand Canyon were made at steady releases of 425 cubic meters per second and at unsteady releases with a daily mean of about 425 cubic meters per second. In Glen Canyon, average flow velocity through the study reach increased directly with discharge, but dispersion was greatest at the lowest of the three flows measured. In Grand Canyon, average flow velocity varied slightly from subreach to subreach at both steady and unsteady flow but was not significantly different at steady and unsteady flow over the entire study reach. Also, longitudinal dispersion was not significantly different during steady and unsteady flow. Long tails on the time-concentration curves at a site, characteristic of most rivers but not predicted by the one-dimensional theory, were not found in this study. Absence of tails on the curves shows that, at the measured flows, the eddies that are characteristic of the Grand Canyon reach do not trap water for a significant length of time. Data from the measurements were used to calibrate a one-dimensional flow model and a solute-transport model. The combined set of calibrated flow and solute-transport models was then used to predict velocity and dispersion at potential dam-release patterns.     

CHAOTIC FORECASTING OF DISCHARGE TIME SERIES: A CASE STUDY1

ABSTRACT: This paper considers the problem of forecasting the discharge time series of a river by means of a chaotic approach. To this aim, we first check for some evidence of chaotic behavior in the dynamic by considering a set of different procedures, namely, the phase portrait of the attractor, the correlation dimension, and the largest Lyapunov exponent. Their joint application seems to confirm the presence of a nonlinear deterministic dynamic of chaotic type. Second, we consider the so‐called nearest neighbors predictor and we compare it with a classical linear model. By comparing these two predictors, it seems that nonlinear river flow modeling, and in particular chaotic modeling, is an effective method to improve predictions.     

A RATING-CURVE METHOD FOR HYDRODYNAMIC SIMULATIONS IN CONTAMINANT TRANSPORT MODELING1

ABSTRACT: Accurate prediction of hydrodynamics is of great importance to modeling contaminant transport and water quality in a river. Flow conditions are needed in estimating potential exposure contamination levels and the recovery time for a no-action alternative in contaminated sediments remediation. Considering highly meandering characteristics of the Buffalo River, New York, a three-dimensional hydrodynamic model was selected to route upstream flows through the 8-km river section with limited existing information based on the model's fully predictive capability and process-oriented feature. The model was employed to simulate changes in water depth and flow velocity with space and time in response to variation in flow rate and/or water surface elevation at boundaries for given bottom morphometry and initial conditions. Flow conditions of the river reach where historical flow data are not available were computed. A rating-curve approach was developed to meet continuous and event contaminant modeling needs. Rating curves (depth-discharge and velocity-discharge relationships) were constructed at selected stations from the 3-D hydrodynamic simulations of individual flow events. The curves were obtained as steady solutions to an unsteady problem. The rating-curve approach serves to link flow information provided by the hydrodynamic model to a contaminant transport model. With the approach, the linking problem resulting from incompatible model dimensions and grid sizes can be solved. The curves will be used to simulate sediment movement and to predict contaminant fate and transport in the river.