Stream Bank Erosion: In Situ Flume Tests

A new technique for testing the erodibility of cohesive stream banks using an in situ flume is presented. The erosion rate is estimated from direct measurements of bed surface elevations by acoustic sensors. The sediment resuspension rate is obtained using data on sediment concentrations measured by optical backscatter sensors and from water samples. The bed-load contribution to the total erosion rate is evaluated from the conservation equation for sediments. Temporal patterns of erosion and resuspension rates are studied employing stepwise increments of bed-shear stress. The data show that bed load plays a significant role in cohesive bank erosion. The data analysis suggests that erosion and resuspension thresholds observed in experiments were very low or equal to zero. The data support the power type equation for the erosion and resuspension rates with bed-shear stress as the key factor. The data also highlights the potential importance of mud content and water content on erosion.     

Finite Element Modeling of the Time-Dependent Deformation of a Slope Bounding a Hydroelectric Reservoir

The south slope of the Vernago hydroelectric reservoir in northern Italy is subjected to a slow time-dependent movement that produced some structural damages in a service shaft located within it. With their continuation, these movements could also damage the diversion tunnel of the reservoir, leading to possible operational problems for the power plant. Two nonlinear, plane strain, finite element models of the slope were developed to evaluate the effectiveness of a berm to be constructed at the slope toe to reduce its present rate of deformation. They are based on different constitutive laws (namely: Viscoplastic and elastoplastic softening) the parameters of which were calibrated on the basis of the available geotechnical information and through the back analysis of inclinometer measurements. The critical evaluation of the results of the two numerical models singled out the inherent limits of one of them, in spite of the satisfactory outcome of the mentioned back analysis. The conclusion that can be reached from this example is that a successful back analysis does not necessarily guarantee a successful forward analysis.     

Coupling Bank Stability and Bed Deformation Models to Predict Equilibrium Bed Topography in River Bends

Attempts to include river bank erosion predictors within morphological models are becoming increasingly common, but uncertainty surrounds the procedures used to couple bed deformation and bank erosion submodels in a way that maintains the mass continuity of eroded bank sediment. Herein we present a coupling procedure that comprises two discrete elements. First, immediately following bank failure, the slumped debris comes to rest at the bank toe as a planar surface inclined at an “angle of repose.” Second, we use a fractional transport model to simulate the subsequent erosion and deposition of the failed bank material debris. The method is demonstrated with an example in which an equilibrium bed topography model is combined with a river bank erosion model to predict the morphological response of a river bend to a large flow event.     

Case Study of Fluvial Modeling of River Responses to Dam Removal

A plan was made to remove Matilija Dam on the Ventura River. With dam removal, the delta in the reservoir and the downstream channel were expected to undergo major changes in morphology. The FLUVIAL-12 model was employed to simulate reservoir and river channel responses after dam removal. As a first step, the model was calibrated using the Ventura River data to establish its validity. In calibration, the model was used to simulate the fluvial processes starting from the time of dam completion. The simulated sediment deposition above the Matilija Dam matches closely with the deposits measured by the U.S. Bureau of Reclamation. A large amount of sediment was stored in the reservoir; some of the stored sediment would be transported downstream after dam removal. An important consequence of dam removal is the major increase of sediment release to the river channel downstream. The sediment supply to the downstream reach is not only from the stored sediment in the reservoir but also from natural sediment inflow from the upstream watershed. Therefore, sediment supply to the downstream reach will exceed the natural sediment flow before the dam presence. This situation tends to overload the downstream reach with sediment, resulting in excessive deposition. The amount of sediment release from the area above the removed dam is closely related to the changes in reservoir morphology. It is necessary to model changes in the channel bed profile and channel width during erosion in order to determine the amount of sediment removal. The amount of sediment release may not be simulated using an erodible-bed model but it may be determined using an erodible boundary model.     

Reservoir Routing using Steady and Unsteady Flow through Rockfill Dams

In this paper a two-dimensional (2D) model for flow through rockfill dams is presented and its results have been compared to 1D model. The model is based on the 2D continuity equation. In the model, an exponential relationship between Reynolds number (R) and Darcy-Weisbach coefficient (f) is suggested and combined with the continuity equation. Coefficients of this relationship are estimated by using real data and a nonlinear optimization technique. Introducing inflow hydrograph to the reservoir and rockfill characteristics as input data and utilizing the above model the outflow hydrograph can be determined. The model has been calibrated and verified using real data. The results of the numerical solution have been shown to be more reliable than the 1D model. To demonstrate the model sensitivity to different parameters, a parametric sensitivity analysis has been conducted. Finally, a comparison between the steady- and unsteady-state results is introduced.     

Modeling Reservoir Irrigation in Uncertain Hydrologic Environment

In a detailed model for reservoir irrigation taking into account the soil moisture dynamics in the root zone of the crops, the data set for reservoir inflow and rainfall in the command will usually be of sufficient length to enable their variations to be described by probability distributions. However, the potential evapotranspiration of the crop itself depends on the characteristics of the crop and the reference evaporation, the quantification of both being associated with a high degree of uncertainty. The main purpose of this paper is to propose a mathematical programming model to determine the annual relative yield of crops and to determine its reliability, for a single reservoir meant for irrigation of multiple crops, incorporating variations in inflow, rainfall in the command area, and crop consumptive use. The inflow to the reservoir and rainfall in the reservoir command area are treated as random variables, whereas potential evapotranspiration is modeled as a fuzzy set. The model’s application is illustrated with reference to an existing single-reservoir system in Southern India.     

Analytical Approach to Calculate Rate of Bank Erosion

Bank erosion consists of two processes: basal erosion due to fluvial hydraulic force and bank failure under the influence of gravity. Because bank resistance force varies with the degree of saturation of bank material, the probability of bank failure is the probability of the driving force of bank failure being greater than the bank resistance force. The degree of saturation of bank material increases with river stage; therefore, the frequency of bank failure is correlated to the frequency of flooding. Consequently, the rate of bank erosion is due to both basal erosion and bank failure, and bank failure is a probabilistic phenomenon. In this paper, for cohesive bank material experiencing planar bank failure, a deterministic approach was adopted for basal erosion analysis, whereas the probability of bank failure was included in the analysis of bank failure. A method for calculating the rate of bank erosion was derived that integrates both basal erosion and bank failure processes, and accounts for the effects of hydraulic force, bank geometry, bank material properties, and probability of bank failure.     

Modeling Turbidity in a Water Supply Reservoir: Advancements and Issues

The development and testing of a “turbidity” model is documented for a water supply reservoir, Schoharie Reservoir, NY, where inorganic terrigenous particles received during runoff events in turbid density currents from the primary tributary cause distinct periodic degradation. The model state variables are fractions (two or three) of the beam attenuation coefficient at 660?nm (c660), a surrogate optical metric of turbidity. The fractions of c660 correspond to slow and rapidly settling components; the latter implicitly accommodates particle aggregation. The transport framework is a two-dimensional (laterally averaged), independently tested, hydrodynamic model. Model testing is supported by detailed measurements of the dynamics of tributary and meteorological drivers and c660 within the reservoir, during and following twelve runoff events. The model is demonstrated to meet the demanding temporal and spatial predictive needs of water supply lakes and reservoirs, by performing well in simulating the timing and magnitude of c660 peaks, the vertical and longitudinal patterns of c660, diminishment following runoff events, and the dependence of impact on magnitude of a runoff event. Further advancements in turbidity modeling, including multiple particle size classes as state variables and explicit representation of particle aggregation and resuspension inputs, are considered.     

Time-Domain Hydrodynamic Forces on Rigid Dams with Reservoir Bottom Absorption of Energy

In this investigation, a two-dimensional time-domain closed-form mathematical model for the hydrodynamic forces on the upstream vertical face of a given rigid dam subjected to a specified horizontal ground motion accelerogram was developed. The model includes the absorption of energy at the elastic reservoir bottom, characterized by the impedance ratio of the sub-bottom materials with respect to water (α). The formulated boundary-value problem is solved in Laplace’s domain and subsequently transformed back to the time domain. Response spectra for the hydrodynamic base shear force and overturning moment are constructed for extreme values of the parameter α. It is found that, frequently, including the solid-foundation elasticity in the reservoir model attenuates the resultant hydrodynamic forces on a rigid barrier, as compared to the results for the case of a rigid reservoir foundation. In this case, the elasticity of the sub-bottom materials constitutes an effective energy dissipating mechanism (radiation damping). By contrast, for sub-bottom materials with less-than-water impedance, amplification of the effective earthquake forces is obtained, as compared to the results for the case of a rigid reservoir foundation.     

Modeling Resuspension in a Dynamic Water Supply Reservoir

Enhancements to the two-dimensional lake and reservoir water quality model W2Tn to simulate the effects of currents and waves on sediment resuspension and turbidity are described. Bed stress attributable to currents was computed by the hydrothermal component of W2Tn, whereas a surface wave component was added to W2Tn to determine bed stress owing to waves. Resuspension flux is computed from bed stress and is included as a source of turbidity to the water column. The model is tested through application to Schoharie Reservoir, a drinking water supply that experiences episodes of elevated turbidity caused by runoff events and exacerbated by drawdown. Model predictions of bed stress attributed to currents are validated by using measurements obtained from acoustic Doppler instrumentation. The surface wave component of the model is established on a framework that has been previously validated for Schoharie Reservoir. Testing of the enhanced turbidity component of W2Tn was completed for a 3.5-year period of historical observations, which included a number of runoff events covering a range of severity and variations in reservoir drawdown. The enhanced model performed well in simulating observed conditions in the water column. The resuspension mechanism made a significant contribution to the predicted turbidity during periods of reservoir drawdown and during a severe runoff event. The model also performed well in simulating the observed turbidity of the drinking water withdrawal. Resuspension of particles contributing to turbidity was largely attributable to reservoir currents with surface wave-induced resuspension playing a smaller role. The potential application of this model to other water bodies and water quality issues is discussed.     

Modeling of Subsidence and Reservoir Compaction under Waterflood Operations

A full-field, three-dimensional (3D) computer model has been developed to numerically simulate reservoir compaction and surface subsidence for a weak, water-sensitive, hydrocarbon reservoir during field-wide water-injection operations. The developed model was used for modeling the compaction and subsidence processes under waterflood operations at the Ekofisk Field in the North Sea. The model was thoroughly validated through the comparison of model results to extensive field measurements with good agreement being achieved. The validated model has been successfully employed as a tool to forecast subsidence and to assist in the development of a subsidence risk assessment. For practical field applications, important quantitative information, that includes reservoir compaction, seafloor subsidence, and seafloor horizontal movement, may be generated from the full-field, 3D model and is presented in this paper.     

Capillary Fringe and Unsaturated Flow in a Porous Reservoir Bank

Steady, 2D Darcian seepage from a zero-depth reservoir into a homogeneous porous bank is studied analytically. Using the Green-Ampt assumption for hydraulic conductivity as a function of pressure and the Vedernikov model for the tension-saturated zone, a free boundary problem with the capillary fringe spread along the bank surface is solved. Accurate direct calculations are made for the extent of the capillary rise along a vertical and horizontal bank. Unsaturated flow from a vertical phreatic surface into a bank with either an impermeable or isobaric vertical soil slope is analyzed in terms of the Philip model. Explicit expressions for the Darcian velocity components, stream function, and/or Kirchhoff potential are presented. The flow topology and characteristics are shown to depend strongly on the boundary condition at the soil surface.     

Cross-Well Radar. II: Comparison and Experimental Validation of Modeling Channel Transfer Function

Close agreement between theory and experiment is critical for adequate understanding and implementation of the cross-well radar (otherwise known as cross-borehole ground penetrating radar) technique, mentioned in a previous paper by the authors. Comparison of experimental results to simulation using a half-space dyadic Green’s function in the frequency domain requires development of transfer functions to transform the experimental data into a compatible form. A channel transfer function (CTF) was developed to avoid having to model the transmitting and receiving characteristics of the antennas. The CTF considers electromagnetic wave propagation through the intervening media only (soil in this case) and hence corresponds to the simulation results that assume ideal sources and receivers. The CTF is based on assuming the transmitting antenna, soil, and receiving antenna as a cascade of three two-port microwave junctions between the input and output ports of the vector network analyzer used in the experimental measurements. Experimentally determined CTF results are then compared with computational model simulations for cases of relatively dry and saturated sandy soil backgrounds. The results demonstrate a reasonable agreement, supporting both the model and CTF formulation.     

Three-Dimensional Modeling for Estimation of Hydraulic Retention Time in a Reservoir

A three-dimensional computational fluid dynamics model is used to estimate the hydraulic residence time for a portion of the Wachusett Reservoir in central Massachusetts. The basin under consideration has several major inflows and exhibits complex flow patterns. The basin is modeled using the FLUENT software package with particles used to track travel time in a steady-state flow field. A tetrahedral mesh with over 1.6 million cells is used with accurate depiction of basin bathymetry and inlet and outlet geometries. Modeling is performed to simulate behavior during a period when conditions are isothermal. It is determined that mean hydraulic residence time is 3–4?days; approximately half of what would be expected assuming strictly plug flow. The presence of a primary flow path, large scale eddies and stagnation zones contribute to the faster travel times. Reductions in inflow rates produce increased residence times and significant changes in flow patterns.     

Contact Erosion at the Interface between Granular Coarse Soil and Various Base Soils under Tangential Flow Condition

Under embankment dams and dykes, horizontal groundwater seepage prevails. If the subsoil is layered, and if some coarse layers are not appropriate filters for finer layers, there can be contact erosion at the interface between fine and coarse soils. In order to study contact erosion threshold, some base-soil and coarse-soil combinations were submitted to a flow parallel to the interface between the coarse soil and the base soil. Critical velocities and critical hydraulic gradients were measured for various base soils. Using effective base-soil grain diameter, an empirical expression for critical velocity was proposed that is well adapted for silts or sand/clay mixtures as well as for sands. The mass of eroded soil was measured relative to the flow velocity for each base-soil/coarse-soil setup. The shear stress applied to the interface between base soil and coarse soil was derived from the hydraulic gradient. Using an empirical relationship between applied shear stress and measured eroded mass, erosion rate was estimated for each base-soil/coarse-soil setup.     

Life of Maithon Reservoir on Ground of Sedimentation: Case Study in India

This case study discusses the life of the Maithon Reservoir studied based upon the past reservoir sedimentation trends and downstream demands for which the reservoir was built, where increased or projected demands were not considered during the study. Various studies have been done to observe declining storage in different zones of the reservoir, including reliability of the conservation storage, rise of reservoir bed level, reduction of reservoir trap efficiency, etc., at different projected years. An attempt has also been made to discuss the full life of the reservoir and to provide a short discussion of different strategies for its prolonged life because sediment is filling the reservoir much faster than its designed life. An additional life of 58?years for the Maithon Reservoir may be obtained if the Balpahari Reservoir is constructed as a siltation trap at just the upstream of the Maithon Reservoir within the year 2005. Actually, this case study provides a guideline for dealing with a practical sedimentation problem along with the serviceability of a reservoir based upon Murthy’s recommendations.     

Modeling the Evolution of Incised Streams. II: Streambank Erosion

Incision and ensuing widening of alluvial stream channels is widespread in the midsouth and midwestern United States and represents an important form of channel adjustment. Streambanks have been found to contribute as much as 80% of the total suspended load. The location, timing, and magnitude of streambank erosion are difficult to predict. Results from field studies to characterize the resistance of fine-grained materials to hydraulic and geotechnical erosion, the impact of pore-water pressures on failure dimensions and shearing resistance, and the role of riparian vegetation on matric suction, streambank permeability, and shearing resistance are used to enhance the channel evolution model CONCEPTS (conservational channel evolution and pollutant transport system). This paper discusses the conceptualization of the above-mentioned physical processes, and demonstrates the ability of the derived model to simulate streambank-failure processes. The model is tested against observed streambank erosion of a bendway on Goodwin Creek, Miss. between March 1996 and March 2001, where it accurately predicts the rate of retreat of the outside bank of the bendway. The observed change in average channel width within the central section of the bendway is 2.96?m over the simulation period, whereas a retreat of 3.18?m (7.4% larger) is simulated. The observed top-bank retreat within the central section of the bendway is 3.54?m over the simulation period, whereas a retreat of 3.01?m (15% smaller) is simulated.     

Relating Damping to Soil Permeability

Published comparisons of complex moduli in dry and saturated soils have shown that viscous behavior is only evident when a sufficiently massive viscous fluid (like water) is present. That is, the loss tangent is frequency dependent for water saturated specimens, but nearly frequency independent for dry samples. While the Kelvin–Voigt (KV) representation of a soil captures the general viscous behavior using a dashpot, it fails to account for the possibly separate motions of the fluid and frame (there is only a single mass element). An alternative representation which separates the two masses, water and frame, is presented here. This Kelvin–Voigt–Maxwell–Biot (KVMB) model draws on elements of the long standing linear viscoelastic models in a way that connects the viscous damping to permeability and inertial mass coupling. A mathematical mapping between the KV and KVMB representations is derived and permits continued use of the KV model, while retaining an understanding of the separate mass motions.     

Forced Resonant Oscillations as a Response to Periodic Winds in a Stratified Reservoir

The response of the Sau reservoir to a wind field characterized by having marked periodicities of 12?h and 24?h has been studied. Measurements of temperature, with a thermistor string, and currents, with an acoustic Doppler current meter, show that the reservoir also responds with the same water periodicities. During certain times of the stratified period, some of the natural oscillation modes of the reservoir are close to these forcing wind periods. In particular, in mid-July the vertical Mode V2 is close to 12?h and in mid- to end of September the vertical Mode V3 is close to 24?h. In these situations, these modes are selected out of the spectrum of possible internal waves and the reservoir behaves as a forced oscillator in resonance with the wind. The structure and the period of these vertical modes have been elucidated by using the 3D model ELCOM. Both modes are affected by the Earth’s rotation at the widest part of the reservoir.     

Improved Approach to Construct Constitutive Surfaces for Stable-Structured Soils Covering Both Saturated and Unsaturated Conditions

Constitutive surfaces are indispensable for investigation of the behavior of soils. Saturated and unsaturated soils coexist in most engineering problems and it is meaningful to develop constitutive surfaces covering both saturated and unsaturated conditions which help to investigate the behavior for both saturated and unsaturated soils in a unified way. At present, the methodologies used for saturated and unsaturated soils are different and few researchers consider the constitutive surfaces for saturated soils. For unsaturated soils, the suction-controlled triaxial tests are usually laborious, time consuming, costly, and may not justify routine engineering projects. This paper discusses the role of constitutive surfaces in soil mechanics and presents an improved approach over existing interpolation methods to construct the constitutive surfaces for saturated and unsaturated conditions for a stable-structured soil using simple laboratory tests.