Numerical Morphological Modeling of Open-Check Dams

Open-check dams are built in mountain streams to control sediment transport during a flood. Sediment passes through them at the lowest discharges, whereas deposition occurs during the highest discharges. Open-check dams are currently designed based mainly on construction experience. Modeling of hydraulics and bed morphology in check dams involves mixed flows (supercritical and subcritical) as well as discontinuities such as hydraulic jumps. In this paper an unsteady coupled numerical mobile-bed model that can tackle rapid varying flows and discontinuities is applied. The numerical technique is based on the classical staggered grids and implicit integration schemes, together with a proper mass and momentum balance. The 1D numerical model is successfully verified with experimental data of slit-check dams. The applicability of the model in the design of open-check dams is also illustrated.     

Dynamic Modeling of Storm Surge and Inland Flooding in a Texas Coastal Floodplain

The purpose of this study is to analyze the combined effects of storm surge and inland rainfall on the floodplain of a coastal bayou in the Houston area by using dynamic hydraulic modeling. Most existing floodplains in the Houston area are defined using only rainfall as an input into steady-state hydraulic models and do not consider the impact of hurricane-induced storm surge on the floodplain. HEC-RAS, a one-dimensional flow model, was run for both steady- and unsteady-states to analyze the additional effect storm surge has on the coastal floodplain. Storm surge and rainfall data from Hurricane Ike were utilized to run an unsteady hydraulic model on Horsepen Bayou near Galveston Bay. The dynamic model generated a good match between the modeled hydrograph and measured data in the watershed. Additionally, a timing sensitivity analysis was completed by shifting the timing of the storm surge both earlier and later in time. The dynamic model revealed that the timing of both rainfall and storm surge play a significant role in the magnitude of inland flooding.     

Upwind Conservative Scheme for the Saint Venant Equations

An upwind conservative scheme with a weighted average water-surface-gradient approach is proposed to compute one-dimensional open channel flows. The numerical scheme is based on the control volume method. The intercell flux is computed by the one-sided upwind method. The water surface gradient is evaluated by the weighted average of both upwind and downwind gradients. The scheme is tested with various examples, including dam-break problems in channels with rectangular and triangular cross-sections, hydraulic jump, partial dam-break problem, overtopping flow, a steady flow over bump with hydraulic jump, and a dam-break flood case in a natural river valley. Comparisons between numerical and exact solutions or experimental data demonstrated that the proposed scheme is capable of accurately reproducing various open channel flows, including subcritical, supercritical, and transcritical flows. The scheme is inherently robust, stable, and monotone. The scheme does not require any special treatment, such as artificial viscosity or front tracking technique, to capture steep gradients or discontinuities in the solution.     

Quantifying Rainfall and Flooding Impacts on Groundwater Levels in Irrigation Areas: GIS Approach

Landscapes continuously irrigated without proper drainage for a long period of time frequently experience a rise in water-table levels. Waterlogging and salinization of irrigated areas are immediate impacts of this situation in arid areas, especially when groundwater salinity is high. Flooding and heavy rainfall further recharge groundwater and accelerate these impacts. An understanding of regional groundwater dynamics is required to implement land and water management strategies. The purpose of this study is to quantify the impact of flood and rain events on spatial scales using a geographic information system (GIS). This paper presents a case study of shallow water-table levels and salinity problems in the Wakool irrigation district located in the Murray irrigation area with groundwater average electrical conductivity greater than 25,000?μS/cm. This area has experienced several large flood events during the past several decades. Piezometric data are interpolated to generate a water-table surface for each event by applying the Kriging method of spatial interpolation using the linear variogram model. Spatial and temporal analysis of major flood events over the last four decades is conducted using calculated water-table surfaces to quantify the change in groundwater storage and shallow water-table levels. The drainage impact of a subsurface drainage scheme partially covering the area has also been quantified in this paper. The results show that flooding and local rainfall have a significant impact on shallow groundwater. The study also found that postflood climatic conditions (evaporation and rainfall) play a significant role in the groundwater dynamics of the area. The spatial net average groundwater recharge during the flooding events ranges from 0.19 to 0.52?ML/ha. The GIS-based techniques described in this paper can be used for net recharge estimation in semiarid regions where it is important to quantify net recharge impacts of regional flooding and local rainfall. The spatial visualization of the net recharge in a GIS environment can help prioritize management actions by local communities.     

Modeling Soluble Phosphorus Desorption Kinetics in Tile Drainage

This paper describes a simple model for the desorption and transport of soluble reactive phosphorus (SRP) to subsurface drains. The model assumes first-order kinetically rate-limited desorption in a soil surface mixing layer over a soil profile layer that rests on an underlying, shallow restricting layer. Input data include precipitation, soil hydraulic properties, drain outflow, free water surface fluctuation, sorbed P concentrations for the mixing layer and profile, desorption rate and equilibrium soil-SRP partitioning. Model results are compared to data on flow and SRP concentrations in drain outflow collected during natural rainfall events under field conditions. The concentration time series simulated follow the sharp rise, peak, and gradual recession of the observed field data. Predicted event mass loads resulting from observed and simulated tile discharges differ from the observed load by +8.2% and ?9.7%, respectively. Sensitivity analysis indicate that equilibrium assumptions would not provide satisfactory results and that mass transfer limits SRP release to the tile drain.     

Case Study: Dam Safety during Construction, Lessons of the Overtopping Diversion Works at Aguamilpa Dam

A study on risk analysis by overtopping of the diversion works of Aguamilpa Dam in Mexico was carried out to examine the overtopping event of January, 1992. Specifically, this study focuses on the upstream water surface elevation during the flood considering the tunnel discharge actual width and roughness. During this maximum flow event, the upstream cofferdam was overtopped and the discharge tunnels exceeded their maximum hydraulic capacity. The risk analysis was made considering the original hydrologic data (until 1992), the final constructed conditions, to establish a performance function that compares the risk of original deterministic analysis with a probabilistic analysis made in 1992. The study also focuses on practical improvements in the construction stage of tunnels; like lining the floor with hydraulic concrete and shotcrete in vaults and walls, that improves the safety of dams during construction and increases the real return period. The work here is presented as a case study because of the unique large-scale flow conditions that are registered in diversion works of dams.     

Prediction of Rainfall for Short Term Irrigation Planning and Scheduling—Case Study in Victoria, Australia

In the short-term planning (7–14?days) and operation of complex irrigation systems, an estimate of irrigation water demand (IWD) is of a fundamental concern. To predict the IWD, a reliable estimate of the expected rainfall during any irrigation period is of fundamental importance. Rainfall is generally predicted with a certain probability of exceedance. However, the standard flood flow-frequency distributions cannot be used for prediction of rainfall of such short durations because these rainfall series in general consist of zero values. Two methods, the total probability theorem (TPT) with three normalizing transformations (i.e., power, log, and square root), and the leaky law (LL) were used to predict the rainfall of short durations (7–14?days, depending upon the number of irrigations per season) in the Goulburn irrigation area (GIA) of Victoria, Australia. Investigations indicated that the TPT using the power transformation (TPTP) was more effective in modeling the short-term data series than the log (TPTL) and square-root transformations (TPTS). Although the overall fitting of the short-term rainfall data series by the LL method was significantly (99%) better than the TPT method, some series could not be fitted by the LL method. This revealed that the LL method could not model all short-term rainfall data series. Results showed that although both the TPT and LL were quite satisfactory in predicting short-term seasonal rainfall of short durations in the study area, none of them was individually able to model all the short-term rainfall series. Hence, the joint use of TPT and LL methods was recommended for short-term rainfall prediction of the study area.     

Root Zone Moisture Routing and Water Demand Calculations in the Context of Integrated Hydrology

An important issue that integrated hydrologic models (IHMs) for river basins can address is the management of water resources in heavily inhabited and cultivated basins. To address this issue, these models need to simulate water demands and root zone flows in a basin. Irrigation scheduling models (ISMs) have been widely used by professionals to compute farm level water demands and root zone flows. Available ISMs are neither suitable for use at basin scale nor can they be easily linked to IHMs. This paper describes a new model that utilizes methods used by ISMs to compute root zone flows and water demands in river basins and can be linked to IHMs. The model was applied to a basin in California, and the simulated water demands were compared with data compiled for the basin. The differences in the results were attributed to differences in input potential evapotranspiration rates. The paper demonstrates that simulated water demands for rice are very sensitive to saturated soil hydraulic conductivity, whereas demands for other crops are sensitive to the pore size distribution index.     

Experimental Parametric Study of Suffusion and Backward Erosion

Within hydraulic earth structures (dikes, levees, or dams), internal seepage flows can generate the entrainment of the soil grains. Grain transportation affects both particle size distributions and porosity, and changes the mechanical and hydraulic characteristics of the earth’s structure. The occurrence of failures in new earth structures due to internal erosion demonstrates the urgency of improving our knowledge of these phenomena of erosion. With this intention, a new experimental device has been developed that can apply hydraulic stresses to reconstituted consolidated cohesive soils without cracks in order to characterize the erosion evolution processes that might be present. A parametric study was conducted to examine the influence of three critical parameters on clay and sand erosion mechanisms. When the hydraulic gradient was low, it was concluded that the erosion of the structure’s clay fraction was due to suffusion. When the hydraulic gradient increased, it was concluded that the sand fraction erosion initiation was due to backward erosion. The extent of the erosion was dependent on the clay content. The study underlines the complexity of confinement stress effects on both erosion phenomena.     

Hydraulic Models of the Flow Distribution in a Four Branch Open Channel Junction with Supercritical Flow

Intense rainfall on urban areas can generate severe flooding in the city, and if the conditions are right, the flow in the streets can be supercritical. The redistribution of the flow in street intersections determines the flow rates and water levels in the street network. We have investigated the flow that occurs when two supercritical flows collide in a 90° junction formed by streets of identical cross section. Several flow configurations within the intersection are possible, depending on the position of the hydraulic jumps that form in and upstream of the intersection. Previous work has identified three flow types, with Type II flows being further classified into three subregimes. Hydraulic models have been developed, based on the principles of the conservation of flow and momentum flux in the intersection, which predict the angles at which the jumps will form. These models can be used to determine the flow type that will occur. Moreover, additional models have been developed for computing the outflow discharge distribution. For Type I flows, it has not been possible to develop such a hydraulic model for the discharge distribution, but some data are provided for one configuration to indicate the influence of different parameters. For Type II and Type III flows, such models are developed, and their predictions agree with data obtained from the channel intersection facility at the Laboratory of Fluid Mechanics and Acoustics in Lyon.     

Embankment Dam Breach Parameters and Their Uncertainties

Potential flood hazards that would be created by breached embankment dams need to be evaluated to select spillway design floods and to prepare emergency action plans. The breaches are often modeled simply, usually in the shape of a trapezoid that is defined by its final height, base width or average width, and side slopes, along with the time needed for the opening to form completely. Data collected from 74 embankment dam failures were used to develop mathematical expressions for the expected values of the final width and side slope of a trapezoidal breach along with its formation time. Information is provided that allows variances of the predicted quantities to be calculated as well. The findings of the statistical analysis were then applied in a Monte Carlo simulation to estimate the degree of uncertainty of predicted peak flows and water levels downstream from breached embankment dams.     

Hydraulic Jet Control for River Junction Design of Yuen Long Bypass Floodway, Hong Kong

The Yuen Long Bypass Floodway (YLBF) was designed to collect flows from the Sham Chung River (SCR) and the San Hui Nullah (SHN) and to serve as a diversion channel of the Yuen Long Main Nullah (YLMN). Under a 200-year return period design condition, the floodway was designed (1) to divert a flow of approximately 38?m3/s from the supercritical YLMN flow and (2) to convey a total combined flow of 278?m3/s to downstream within acceptable flood levels. The success of the design depends critically on complicated junction flow interactions that cannot be resolved by 1D unsteady flow models. These features include the supercritical-subcritical flow transition at the San Hui-Floodway (SHN-YLBF) junction and the diversion of part of the supercritical flow from the Main Nullah (YLMN). A laboratory Froude scale physical model was constructed to study water stages and flow characteristics in the floodway and to investigate optimal design arrangements at channel junctions and transitions. This paper summarizes the main features of the unique river junction network, in particular the use of the hydraulic jet principle at the SHN-YLBF junction to lower flood levels. In addition, a numerical flow model is employed to study flow details at the river junctions. The model is based on the general 2D shallow water equations in strong conservation form. The equations are discretized using the total variation diminishing finite-volume method which captures the discontinuity in hydraulic jumps. The numerical model predictions are well supported by the laboratory data, and the theoretical and experimental results offer useful insights for the design of urban flood control schemes under tight space constraints.     

Experimental Study of a Right-Angled End Junction between a Pipe and an Open Channel

This paper presents an experimental study of a junction between a closed conduit and an open channel. This study was undertaken to explore hydraulic properties of outlets of subsurface drainage or sewage networks into an open air stream during flood events. Experiments were conducted in a laboratory flume, with a main rectangular channel joined at right angle to a lateral circular pipe. Both branches were supplied with independent flow rates and downstream water level was controlled by an adjustable weir. Several flow patterns were identified, combining free-surface and pressurized flows. Transitions between these flow patterns, as well as changes in water level or energy, in response to the modifications of experimental variables, were studied and could be linked to known properties of single channels, single pipes, and homogeneous junctions. Transitions between free-surface and pressurized pipe flow appeared to be strongly dependent on the whole set of experimental variables and the pipe longitudinal slope. This work contributes to a better knowledge of hydraulic and hydrologic key processes for point source discharging.     

Scour Vulnerability of River Bridge Piers

A simple procedure is proposed to assess the vulnerability of bridge piers in rivers, taking into account the phenomena governing fluvial dynamics during flood events. The procedure requires an estimation of the maximum scour depth of the soil surrounding both the pier and the foundation as well as an analysis of the bearing capacity of the pier–foundation–soil geotechnical system. The scour depth is determined in terms of the physical and mechanical properties of the streambed soil, the shape of the pier foundation and the destabilizing effects induced by hydrodynamic forces. The coupling of both the hydraulic and geotechnical analyses enables to identify the most significant factors characterizing scour depth and affecting pier vulnerability. Two levels (low, medium) of allowable vulnerability, bounded by an extreme condition of high vulnerability, are defined and analytically determined in function of the maximum scour depth and the foundation depth. Specific diagrams corresponding to each category of foreseen actions allow a quick evaluation of the vulnerability of a bridge pier.     

Rockfill Modulus and Settlement of Concrete Face Rockfill Dams

A method is presented for estimating the modulus of compacted rockfill in dams based on the particle size, unconfined compressive strength of the rock, compaction layer thickness, compactive effort, and the applied vertical stress. Also presented are methods for predicting the crest settlement and face slab deformation of concrete face rockfill dams during first filling and in the long term. It is demonstrated that the modulus is stress dependent and guidance is provided on how to assess this, as well as effects of construction in narrow valleys where arching may affect the stresses in the dam. These methods are based on analysis of the 35 dams with good quality information on construction materials and placement methods, and good quality internal and surface settlement monitoring records.     

Effects of Hysteresis on Steady-State Infiltration in Unsaturated Slopes

Hysteresis is a common feature exhibited in hydraulic properties of an unsaturated soil. For a specific matric suction, water content or coefficient of permeability on a wetting curve is always lower than those found on a drying curve. This paper focuses on hysteresis observed in steady-state infiltration tests in a laboratory slope model. The slope model consisted of a 400 mm thick fine sand layer overlying a 200 mm thick gravelly sand layer at a slope angle of 30°. The slope model was subjected to artificial rainfalls of different intensities. The slope model was instrumented to continuously measure the changes in pore-water pressure or matric suction, volumetric water content, and water balance during an experiment. Two experiments with similar applied precipitation intensities were conducted on soils that experienced adsorption and desorption processes. For the adsorption process, the slope model was first subjected to an antecedent steady-state rainfall with an intensity lower than the intensity of the incident steady-state rainfall. In the adsorption process, the water content of the soils increased during the incident rainfall prior to achieving the steady-state condition. For the desorption process, the slope model was first subjected to an antecedent steady-state rainfall with an intensity higher than the intensity of the incident steady-state rainfall. In the desorption process, the water content of the soils actually decreased during the incident rainfall prior to achieving the steady-state condition. The results indicate that the matric suction distributions in soils experiencing the desorption process were higher than those observed in soils experiencing the adsorption process. The matric suctions within the slope during a steady-state infiltration were affected by the initial water content of the soil prior to the infiltration process. Numerical analyses, employing both drying and wetting hydraulic properties of the soils, were performed to study the difference in matric suctions as observed in the experiments. The results suggest that the hysteretic behavior of the soil affects the matric suction distribution within the slope at steady-state conditions. The appropriate hydraulic properties of the soils (i.e., drying or wetting) should be used in accordance with the process that the soils actually experience (i.e., desorption process or adsorption process) even though the slope is under a steady-state rainfall condition.     

Is it worth a dam?

Once a sign of modernization and growth, dams are often seen today as symbols of environmental and social devastation. Over 800,000 dams have been built worldwide to provide drinking water, flood control, hydropower, irrigation, navigation, and water storage. Dams do indeed provide these things,but at the cost of several adverse, unexpected effects: disruption of ecosystems, decline of fish stocks, forced human and animal resettlements, and diseases such as malaria, which are borne by vectors that thrive in quiet waters.     

Riparian Vegetation Mapping for Hydraulic Roughness Estimation Using Very High Resolution Remote Sensing Data Fusion

For detailed hydraulic modeling, accurate spatial information of riparian vegetation patterns needs to be derived in automatic fashion. We propose a supervised classification for heterogeneous riparian corridors with a low number of spectrally separate classes using data fusion of a Quickbird image and LIDAR data. The approach considers nine land cover classes including three woody riparian species, brush, cultivated areas, grassland, urban infrastructures, bare soil and water. The classical “stacked vector” approach is adopted for data fusion, while the nonparametric weighted feature-extraction method and the pixel-oriented maximum likelihood algorithm are used for feature-reduction and classification purposes, respectively. We test the approach over a 14-km stretch of the Sieve River (Tuscany Region, Italy). A one-dimensional river modeling is applied over the study reach comparing the results of a classification-derived hydraulic roughness map and a traditional ground-based approach. Despite the complex study reach, the classification method produced encouraging accuracies (OKS = 0.77) and represents a useful tool to delineate application domains of flow resistance models suited to different hydrodynamic patterns (e.g., stiff/flexible vegetation). Hydraulic modeling results showed that the remotely derived floodplain roughness parameterization captures the equivalent Manning coefficient over 20 test cross sections with uncertainty distributions described by low mean and standard deviation values.     

Solution for the Closed-Loop Problem in Pressurized Multipipe Systems

The fundamental background of the solution for the steady-state flow in pressurized water closed-loop pipe systems, without any reservoirs in-between, is presented in this paper. The use of the steady-state mass balance and energy equations to calculate discharges and heads in this type of hydraulic system leads to an undetermined problem. The way to solve this indeterminacy is to consider an additional continuity equation associated with the difference between initial and final conditions, taking into account fluid compressibility and pipe-wall deformability. A complete formulation is derived considering pressure and temperature changes in the hydraulic system. Simplified formulae are presented for isothermal flows in simple systems and multiple closed-loops with pipes in series and in parallel. This problem can also be solved by a pseudotransient analysis technique applied to steady-state conditions. Proposed solutions for this problem are applied to steady-state flows and tested for different system configurations.