Erosion of Cohesive Sediments: Resuspension, Bed Load, and Erosion Patterns from Field Experiments

New field data on cohesive sediment erosion is presented and discussed, with particular focus on partitioning the total erosion into resuspension and bed load. The data were obtained using a recently developed in situ flume of the National Institute of Water and Atmospheric Research, New Zealand. The erosion rate is estimated from direct measurements of bed surface elevations by acoustic sensors, whereas resuspension rate is obtained using data on sediment concentrations measured by optical backscatter sensors. The bed- load contribution to the total erosion rate is evaluated from the conservation equation for sediments. To test repeatability, the data from the in situ flume are compared with those from a previous version of the flume. The results show that comparative studies of in situ flumes and standardized deployment procedures enable direct comparison of experimental data on cohesive sediment erosion. Overall, the data show that a commonly used assumption that the erosion rate is equal to the resuspension rate is not always valid as bed load plays a significant role in cohesive sediment erosion. The data also highlight the importance of clay content and other sediment physical characteristics in the sediment mixture.     

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.     

Straight Benthic Flow-Through Flume for In Situ Measurement of Cohesive Sediment Dynamics

The design of a straight benthic flow-through flume for in situ studies of cohesive sediment dynamics is described including the flume structure and probes installed for routine measurements of suspended sediments and flow velocity. The flume was calibrated for two roughness types covering the range of possible cohesive bed roughnesses. The calibration included a set of three-dimensional velocity measurements using acoustic Doppler velocimeter. These measurements were used to develop calibration relationships between the bed shear stress (which is difficult to measure directly in routine deployments) and the flume centerline flow velocity, which is routinely measured. An example of a successful deployment of the flume is presented. The limitations and potential for further improvements are also briefly discussed.     

Cohesive Sediment Characterization by Combined Sedimentation and Rheological Measurements

A laboratory characterization of cohesive sediment has been carried out in which data obtained from standard sedimentation and rheological measurements were combined in a determination of the critical solid concentration for the detection of elasticity in a weakly cohesive suspension. The corresponding storage modulus and shear stress are very critical in any in situ rheometry of sediments, especially in the study of mud-water surface erosion in a flume. Sedimentation results showed that particle size distribution rather than surface treatment controlled the rheological behavior of the suspension while the critical solid concentration for the appearance of three-dimensional space-filling network, showing some measurable elasticity in the suspension, occurred in the region of 0.015. This parallel between the consolidation behavior and shear rheology development for the flocculating system has been established. This technique could be an adjunct to the laboratory characterization of cohesive sediments for the estimation of critical shear stress for surface erosion, especially in a typical flume experiment under water wave pressure.     

Flume Measurements of Sediment Erodibility in Boston Harbor

To obtain in situ measurements of sediment erodibility in defined bottom shear stress environments, a portable, straight flume was built, tested, and deployed in the field for six experiments at three locations in Quincy Bay of Boston Harbor, Mass. The flume had a 1.0-m-long inlet section, which included a boundary-layer trip and a roughened, plexiglass bottom; this design prevented erosion of the sediment bed in the boundary-layer-development region. Downstream of the inlet section was a 1.2-m-long sediment test section, which had a laboratory-verified, uniform bottom stress. In the absence of algal mats, our flume experiments on sites exhibiting a range of bed properties indicated quite uniform erodibility, with a critical shear stress τc of 0.10 ± 0.04 Pa and an erosion rate constant M of 3.2 ± 0.2 × 10?3 kg m?2 s?1 Pa?1 (R2 = 0.92, N = 17, where N is the total number of erosion rate measurements made in the absence of algal mats). The measured rates were consistent with those of many other in situ studies. We observed markedly reduced erodibility in early October 1995 when the sediment was covered by a benthic diatom mat, and measured erosion rates were lessened by 50–80%. The possibility of depth-dependent sediment erodibility in near surface (top 3 mm) was investigated by calculating a set of depth-dependent erosion parameters. The parameters obtained suggested that both the critical shear stress and the erosion rate constant were depth-sensitive (both doubling by 1 mm into the sediment).     

Sediment Erosion Characteristics in the Anacostia River

Four in situ experiments on sediment erosion characteristics were conducted at the Anacostia River that runs through Washington, D.C. Supplemental erosion rate data were also obtained by carrying out five laboratory experiments using sediment samples collected at the field. In laboratory experiments, the sediment samples were mixed with tap water and placed in the flume to form beds for finding the difference in terms of erosion characteristics caused by different sediment composition among the five samples. This approach enables the finding of erosion characteristics for the entire tidal Anacostia River with limited resources. The in situ measured critical bed-shear stresses τcr for erosion at the water-sediment interface z = 0 varies from 0.03 to 0.08?Pa. Field results indicated that τcr(z) increases with the depth z and becomes more than 0.6 to 0.7?Pa with an erosion thickness of less than 1?cm. Sediment beds prepared at a laboratory appear having an upper limit on how much τcr(z) can be developed.     

Flume Test Section Length and Sediment Erodibility

This paper analyzes the effect of flume test section length on sediment erodibility measurements. A modular flume was constructed and experiments were conducted with two test section lengths: 0.15 and 1.10?m. The internal height and width of the flume were 0.11 and 0.13?m, respectively. A fine (7?μm) commercially available quartz sediment was used for the tests. The expectation was that the shorter flume test section would experience a significantly higher erosion rate (per unit surface area) due to its greater sensitivity to edge effects (i.e., scour) at the entrance and exit of the flume test section. However, the measured erosion rates at comparable bottom stresses were only 35% greater in the short test-section tests. These results were consistent with the lack of significant scour development at the entrance or exits of the test sections. Hence, flume test section length alone does not appear to significantly affect erodibility measurements provided edge effects (i.e., scour) are minor.     

Case Study: Two-Dimensional Model Simulation of Channel Migration Processes in West Jordan River, Utah

The paper presents the application of a two-dimensional depth-averaged numerical model to simulate the lateral migration processes of a meandering reach in the West Jordan River in the state of Utah. A new bank erosion model was developed and then integrated with a depth-averaged two-dimensional hydrodynamic model. The rate of bank erosion is determined by bed degradation, lateral erosion, and bank failure. Because bank material in the West Jordan River is stratified with layers of cohesive and noncohesive materials, a specific bank erosion model was developed to consider stratified layers in the bank surface. This bank erosion model distinguishes itself from other models by relating bank erosion rate with not only flow but also sediment transport near the bank. Additionally, bank height, slope, and thickness of two layers in the bank surface were considered when calculating the rate of bank erosion. The developed model was then applied to simulate the processes of meandering migration in the study reach from 1981 to 1992. Historical real-time hydrographic data, as well as field survey data of channel geometry and bed and bank materials, were used as the input data. Simulated cross-sectional geometries after this 12-year period agreed with field measurements, and the R2 value for predicting thalweg elevation and bank shift are 0.881 and 0.706, respectively.     

Test of a Method to Calculate Near-Bank Velocity and Boundary Shear Stress

Detailed knowledge of the flow and boundary shear stress fields near the banks of natural channels is essential for making accurate calculations of rates of near-bank sediment transport and geomorphic adjustment. This paper presents a high-resolution laboratory data set of velocity and boundary shear stress measurements and uses it to test a relatively simple, fully predictive, numerical method for determining these distributions across the cross-section of a straight channel. The measurements are made in a flume with a fairly complex cross-section that includes a simulated, cobble-roughened floodplain. The method tested is that reported by Kean and Smith in Riparian Vegetation and Fluvial Geomorphology in 2004, which is modified here to include the effects of drag on clasts located in the channel. The calculated patterns of velocity and boundary shear stress are shown to be in reasonable agreement with the measurements. The principal differences between the measured and calculated fields are the result of secondary circulations, which are not included in the calculation. Better agreement with the structure of the measured streamwise velocity field is obtained by distorting the calculated flow field with the measured secondary flow. Calculations for a variety of narrower and wider configurations of the original flume geometry are used to show how the width-to-depth ratio affects the distribution of velocity and boundary shear stress across the channel.     

Turbulent Effects on the Settling Velocity of Suspended Sediment

The mean settling velocities of suspended sediments in turbulence have been examined. The settling velocities in a flume are directly measured by using an acoustic Doppler velocimeter. The results indicate the same trend as previous work in homogeneous isotropic turbulence. In addition to the flume experiment, the numerical experiments were conducted in the velocity field of homogeneous isotropic turbulence simulated by Kraichnan’s technique. The experimental and numerical results show that at high turbulence intensity the relative settling velocity increases with the increasing relative turbulence intensity regardless of the Stokes number. At intermediate turbulence intensity, it seems that the settling data bifurcate, i.e., the particles at the large Stokes number tend to be slowed, whereas the settling velocity of particles is increased at the small Stokes number.     

Experimental Study of the Flow Field over Bottom Intake Racks

Bottom racks made by longitudinal bars are hydraulic structures widely adopted for engineering purposes. In the present paper we revisit the problem of their hydraulic design, analyzing the data obtained from a systematic series of experiments carried out in a laboratory flume. For each run we measured the diverted discharge, the water surface longitudinal profile, and using a two-dimensional backscatter laser Doppler anemometer, we measured the velocity field over the rack and in the slit between two adjacent bars. The latter measurements, in particular, allow us to obtain the along-rack distributions of the discharge coefficient to be used to determine the rate of change of the diverted discharge. We use such distributions to derive a physically based relationship relating the overall diverted discharge to the length of the rack, the void ratio, the discharge coefficient measured under static conditions, the specific head of the stream approaching the rack, and a modified Froude number. The robustness of the proposed relationship is confirmed by the comparison between the discharges calculated through the proposed relationship and those measured in an extensive series of experiments available in the literature, characterized by ranges of the relevant flow parameters much larger than those investigated in the present contribution.     

Flow Measurement Using Flying ADV Probes

An innovative measurement system of “flying” acoustic Doppler velocimeters was designed in order to allow rapid velocity measurement over a large flow field. Such measurements are necessary, for example, when measuring over a temporally varying and locally nonuniform rough bed. The measurement technique was verified by comparison with measurements taken in the same flows using a traditional stationary probe technique. Comparison showed that the flying-probe approach performs similarly to stationary measurements in capturing the mean flow field and turbulent fluctuations. The data obtained from flying probe experiments can be used to describe the flow in terms of double-averaged hydrodynamic variables, obtained by averaging in time and spatial domains within a thin slab parallel to the mean bed. Examples are presented of flow measurements over a fixed flat bed, a fixed dune bed, and over mobile developing bed forms. It is shown that near-bed measurements suffer from boundary reflection interference, though affected data can be filtered out based on the ADV-measured correlation coefficient. Measurement below roughness tops is possible, with in-bed records being detectable by spikes in measured signal-to-noise-ratio and by comparison with measured bed topography.     

Field Testing of a Simple Flume (SMBF) for Flow Measurement in Open Channels

In this paper the stage–discharge relationship of a new flume named SMBF (Samani, Magallanex, Baiamonte, Ferro), originally proposed by Samani and Magallanez and tested by Baiamonte and Ferro, for measuring flow discharge in open channels is reviewed. The flume is obtained inserting two semicylinders in a rectangular cross section. The results of some experimental runs carried out using horizontal flumes characterized by different values of the contraction ratio (ranging from 0.17 to 0.81) are used for determining the two coefficients of the power stage–discharge equation. The stage–discharge equation is tested using flow measurements carried out in the period between December 2004 and March 2006 in the Sicilian experimental SPA1 basin. Field testing of the SMBF flume is developed using discharge measurements carried out by a Khafagi–Venturi flume placed in the field measurement channel.     

Automated Sediment Erosion Testing System Using Digital Imaging

Measurement of vertical profiles of the critical shear stress, τc, and the erosion rate, E, from the same undisturbed sediment core is crucial for modeling the resuspension of fine-grained natural sediments. The automated sediment erosion testing system (ASETS) was developed to determine profiles of τc and E with centimeter spatial (vertical) resolution in an undisturbed (Shelby tube) sediment core, whose surface was eroded by steady turbulent flow through a flume. The unique feature of ASETS is that it is a real-time imaging method that accurately determines the position of the core surface during erosion for both calculating the vertical profile of E and controlling a motor-driver system that automatically pushes up the core to maintain its surface flush with the flume bottom. Undisturbed, field cores were tested over a range of flow (average bed shear stress, τb) conditions. The amount of eroded sediment from both optical backscattering measurements and the imaging method were in good agreement, which validated ASETS. Measured vertical profiles of τc and E were similar to those reported in literature. E correlated well with (τb?τc)2, which agrees with previous results in literature.     

Critical Bed-Shear Stress for Cohesive Sediment Deposition under Steady Flows

Results of two laboratory experiments on the cohesive sediment deposition behavior are presented. Data indicate that an annular flume, either with or without the channel bed and the top ring rotating in opposite directions to minimize the secondary circulation, can be used to study the deposition behavior of cohesive sediments. Direct observations on when and where the bed is formed suggest that deposition only occurs when the local bed-shear stress (τb) is less than a critical value. Secondary circulation in the flume, which produced upward current near the inner corner and downward current near the outer corner, did not prevent deposition near the inner corner (because of the small τb) nor promote deposition near the outer corner (because of the relatively large τb).     

New Method for Estimation of Discharge

A new technique for drawing isovel patterns in an open or closed channel is presented. It is assumed that the velocity at each arbitrary point in the conduit is affected by the hydraulic characteristics of the boundary. While any velocity profile can be applied to the model, a power-law formula is used here. In addition to the isovels patterns, the energy and momentum correction factors (α and β), the ratio of mean to maximum velocity (V/umax), and the position of the maximum velocity are calculated. To examine the results obtained, the model was applied to a pipe with a circular cross section. A comparison between the profiles of the proposed model and the available power-law profile indicated that the two profiles were coincident with each other over the majority of the cross section. Furthermore, the predicted isovels were compared with velocity measurements in the main flow direction obtained along the centerline and lateral direction of a rectangular flume. The estimated discharge, based on measured points on the upper half of the flow depth away from the boundaries was within ±7% of the measured and much better in comparison to the prediction of one- and two-point methods. The prediction of the depth-averaged velocity values for the River Severn in the United Kingdom shows a good agreement with the measured data and the best analytical results obtained by the depth-averaged Navier–Stokes equations.     

Statistically Derived Bedload Formula for Any Fraction of Nonuniform Sediment

Based on a method of combining stochastic processes with mechanics, a new bedload formula for the arbitrary kth size fraction of nonuniform sediment is theoretically developed by using a stochastic model of sediment exchange and the probabilistic distribution of fractional bedload transport rates. The relations, proposed recently by Sun, for the probability of fractional incipient motion and for the average velocity of particle motion are introduced to bedload formula. Plenty of experimental data for the bedload transport rate of uniform sediment are used to determine two constants. The theoretical bedload formula for any fraction of nonuniform sediment possesses several advantages, including a clear physical concept, a strict mathematical derivation, and a self-adaptability to uniform sediment. The formula is verified with natural data expressing the transport of nonuniform sediment under full motion in laboratory flume. The result shows that the experimental observations agree well with the predicted fractional bedload transport rates. Comparison of the theory with field data finds that the proposed formula still applies to partial transport of bedload in gravel-bed streams as long as the immobile percentage of bed material is taken into account.     

SRICOS: Prediction of Scour Rate in Cohesive Soils at Bridge Piers

A new method called SRICOS is proposed to predict the scour depth z versus time t around a cylindrical bridge pier of diameter D founded in clay. The steps involved are: (1) taking samples at the bridge pier site; (2) testing them in an erosion function apparatus to obtain the scour rate ? versus the hydraulic shear stress applied τ; (3) predicting the maximum shear stress τmax, which will be induced around the pier by the water flowing at vo before the scour hole starts to develop; (4) using the measured ? versus τ curve to obtain the initial scour rate ?i corresponding to τmax; (5) predicting the maximum depth of scour zmax for the pier; (6) using ?i and zmax to develop the hyperbolic function describing the scour depth z versus time t curve; and (7) reading the z versus t curve at a time corresponding to the duration of the flood to find the scour depth that will develop around that pier. A new apparatus is developed to measure the ? versus t curve of step 2, a series of advanced numerical simulations are performed to develop an equation for the τmax value of step 3, and a series of flume tests are performed to develop an equation for the zmax value of step 5. The method is evaluated by comparing predictions and measurements in 42 flume experiments.     

Implicit Bidiagonal Scheme for Depth-Averaged Free-Surface Flow Equations

A general fast implicit bidiagonal numerical scheme, based on the MacCormack's predictor-corrector technique requiring the inversion of only block bidiagonal matrices, has been developed and subsequently applied for subcritical and supercritical free-surface flow calculations. The model has been applied to depth-averaged steady flows. There are two main advantages of the proposed method: the technique has fast convergence and utilizes a body fitted nonorthogonal local coordinate system to simulate irregular geometry flows. The model is used to analyze a wide variety of hydraulic engineering problems including flows in a converging-diverging subcritical flume, supercritical expansions at various Froude numbers, and supercritical converging chutes. For each of these test cases, the calculated results are compared with experimental data. The comparisons with measurements as well as with other numerical solutions show that the proposed method is comparatively accurate, fast, and reliable.