Principles of spontaneous bedform generation

Bedforms often form spontaneously when water flows over a sandy sediment bed. Processes responsible for the instability of flat sediment beds are discussed in several Coastal Wiki articles listed below under #Related articles. There are several types of seafloor instability, which give rise to the great diversity of bedforms observed in the field. The common denominator of bedform generation processes is the adjustment of the flow field to an initial random disturbance of the seafloor, such that gradients in the corresponding adjusted sediment flux amplify the initial disturbance. A mathematical derivation of bedform generation processes is generally based on the identification of leading nonlinearities in the equations describing the flow adaptation to an infinitesimally small disturbance of the seafloor. This article aims to provide a qualitative and more intuitive understanding of bedform generation processes. The focus is on bedforms with a rhythmic character with crest and trough lines perpendicular to the undisturbed flow, such as bed ripples, subaqueous dunes and tidal sand waves.

Bedload transport maximum

The principle of spontaneous bedform generation is schematically illustrated in a series of figures. Fig. 1 shows a vertical longitudinal section of a smooth 2-dimensional bed perturbation in the form of a small ridge and a shallow trench. A steady flow from left to right adapts to the bed perturbation as shown by the streamlines. The fluid layer near the bottom where the streamlines are deflected by the bed perturbation is called the boundary layer. The bending over the ridge is opposed by the inertia of the flow (inertia is the momentum carried by the flow, a nonlinear effect). The streamlines following the bottom therefore almost collide with the upstream ridge flank. At the location where they come closest to the ridge flank, the sediment experiences the strongest shear stress. This is the location where the capacity of the flow to move sediment along the bottom is greatest.

Fig. 1. Steady flow streamlines over an emerging bed ripple in a vertical longitudinal section. The streamlines in the boundary layer contract on the upstream side of the ripple and diverge at the downstream side. The velocity shear is maximum where the contraction of the streamlines is greatest.

Fig. 2. The emerging ripple is eroded where the shear stress increases due to contraction of the streamline toward the sediment bed. Accretion occurs where streamlines diverge from the sediment bed. The erosion-accretion boundary is shifted by the bedload saturation length [math]l_{sat}[/math], the distance over which the sediment flux adjusts to a change in shear stress. The ripple moves downstream as its height increases.

Sediment transport along the bottom by rolling or jumping of sediment grains is called bedload transport. However, it takes a while before the full bedload transport capacity is realized. The location of maximum bedload transport is therefore located a short distance downstream of the location of maximum bed shear stress (Fig. 2) . Before the maximum is reached, the bedload transport increases on the upstream flank of the ridge, by moving an increasing amount of sediment grains. The ridge flank upstream of the location of maximum bedload transport therefore loses sediment and erodes. Downstream of the location of maximum bedload transport, the capacity of the flow to move sediment grains decreases. Sediment grains are therefore left behind, which means accretion of the ridge downstream of the location of maximum bedload transport - the ridge migrates in a downstream direction. If the location of maximum bedload transport is located upstream of the ridge crest, the crest will accrete. In this case, the bed perturbation grows as it migrates downstream. Figures 1 and 2 illustrate the essential principle underlying the spontaneous formation of bedforms in steady flow.

Residual circulation cells

Fig. 3. The velocity profiles at the ripple trough and the ripple crest corresponding to the contraction-divergence of the streamlines.

Figure 3 shows the adjusted velocity profiles at the ridge crest and trough. The boundary layer is thinner above the ridge crest than above the trough, implying flow acceleration-deceleration. At the seabed (below the bedload transport layer) the velocity is zero. The fluid mass in the boundary layer is thus bending over, which means that is has a rotational component called vorticity. The near-bed clockwise vorticity is clearly greater at the ridge crest than at the trough. Figure 4 shows the transport of vorticity at the ridge crest and the trough. The boundary layer region from trough to crest experiences loss of clockwise vorticity, while the boundary layer region from crest to trough experiences a gain (a nonlinear second-order effect). The loss of clockwise vorticity in the trough-to-crest region drives a small counterclockwise rotational flow, superimposed on the main flow shown in Fig. 5. The gain of clockwise vorticity in the crest-to-trough region likewise drives a clockwise rotational flow component. These superimposed circulations contribute to moving sediment grains from the trench to the ridge crest.

Fig. 4. Transport of vorticity (rotational component of velocity) at the ripple trough and ripple crest. The vorticity balance on the upstream side of the ripple is negative (more clockwise vorticity exits than enters), while the vorticity balance is positive on the downstream side.

Fig. 5. The negative vorticity balance on the upstream side of the ripple drives a counterclockwise circulation, while the positive balance on the downstream side drives a clockwise circulation. This circulation pattern, superimposed on the steady main flow, contributes to moving sediment particles from the ripple trough to the ripple crest.

If the main steady flow is reversed (Fig. 6), the resulting adjusted flow pattern is the exact mirror of Fig 5. All flow directions are reversed except the superimposed circulation pattern resulting from the vorticity transport balance. Figure 7 shows the net flow pattern that results from considering successive periods of alternating steady flow. In this case no bedform migration occurs. However, the bed perturbation can grow because the residual circulation pattern in the boundary layer contributes to moving sediment grains from the trench to the ridge.

Equating oscillating flow (e.g. horizontal wave orbital flow, or tidal flow) with successive periods of alternating steady flow is not entirely correct. The boundary layer structure in oscillating flow differs from steady flow, primarily because of the time required for a turbulent boundary layer to form after flow reversal. For surface waves, the thickness of the turbulent boundary layer is therefore much smaller than for steady flow. However, mathematical models taking the temporal dynamics of the boundary layer into account yield a residual near-bed circulation pattern similar to Fig. 7.

Fig. 6. Flow pattern adjusted to a bed perturbation similar as in Fig. 5, but with the flow direction reversed.

Fig. 7. If oscillating flow is assumed to be equivalent to successive periods of steady flow from the right and from the left, then the net resulting flow pattern consists of circulation cells that move sediment grains from the troughs to the crests.

Ripple coarsening in steady flow

Ripples are the first bedforms to arise spontaneously when water flows over a sandy sediment bed. They form from a small perturbation of a flat sediment bed when the flow velocity is just above the critical value for setting sediment grains in motion. The bed shear stress maximum is located on the upstream flank of the perturbation at a distance of the order of [math] \lambda / 8[/math] from the crest (see Appendix; [math]\lambda[/math] is the wavelength of the perturbation). Assuming that the bedload saturation length [math]l_{sat}[/math] is 10-100 times the median sediment grain diameter [math]d_{50}[/math] , then the fastest growing perturbations have a wavelength of 100-1000 [math]d_{50}[/math] (see Appendix). These emergent ripples (wavelets) will grow in height, but also in wavelength. The migration speed of the ripples depends on the height and on the wavelength. As not all the ripples have the same height and wavelength, ripples tend to coalesce into larger ripples with greater wavelengths. This process is called ripple coarsening.

The seabed is then no longer smooth; the flow along the bottom becomes turbulent and the boundary layer extends to the free water surface. At this development stage, the streamlines will not follow the bed everywhere, but they separate at the ripple crest, to which a recirculating vortex is attached at the downstream side. Sediment suspension by the vortex greatly increases the saturation length. The location of maximum sediment transport may be shifted downstream of the ripple crest, resulting in ripple decay. Wavelength growth is also inhibited when the ripple wavelength [math]\lambda = 2 \pi / k[/math] reaches values at which surface water waves become resonant with the bed perturbation. When the current velocity [math]u[/math] equals or exceeds the wave celerity [math]c= \sqrt{g \, \tanh(kh)/k}[/math], the flow becomes supercritical and antidunes appear. However, the development of subaqueous dunes from bed perturbations with wavelengths [math]\lambda[/math] much larger than the water depth [math]h[/math] is possible for rough (rippled) sediment beds. The dune wavelength scales with the water depth; the fastest growing perturbations have wavelengths in the range [math]\lambda = 6-60 \; h[/math]. The development of dunes with smaller wavelengths is suppressed by gravity-driven transport from crest to trough.

Wavelength selection in oscillating flow

Spontaneous bedform generation in oscillating flow is generally assumed to be driven by the residual circulations shown in Fig. 7. These residual circulation cells exist for bed perturbations of any wavelength. However, periodic bedforms in field situations only occur at certain wavelength ranges. Wave ripples typically have a wavelength of 5 to 50 cm, while tidal sand waves (submarine dunes in tidal seas) have a wavelength of 100 to 1000 m. The strength of the residual circulation cells differs for different wavelengths, but this is not the only explanation of the wavelength ranges. Bedload transport is generally assumed to have a downslope component that counteracts the contribution of the residual circulation cells to the upslope movement of sediment grains. The downslope transport component is stronger for bed perturbations with a small wavelength (because they have steeper slopes). Thus, the growth of bed perturbations with a small wavelength will be slow or even impossible. On the other hand, if the wavelength of the bed perturbation is large, the residual circulation cells are weak and the distance for sediment grains to move from trough to crest is large. So, bed perturbations will start developing ripples with intermediate wavelengths.

Fig. 8. Sand transport over a vortex bed ripple in the Aberdeen oscillatory flow tunnel. The orbital velocity to the right creates a sand-laden vortex on the lee-side of the ripple crest, which is swept back over the ripple right after flow reversal. Courtesy: Jebbe van der Werf.

When the ripple height grows, flow separation at the ripple crest will generate vortexes laden with sediment that are ejected into the flow (Fig. 8). The resulting modification (inversion) of the sediment fluxes creates zones of erosion and deposition, allowing for the creation and development of new ripples. Merging occurs when two ripples are situated between two other ripples that generate vortices of a larger magnitude. By this mechanism, a field of wave ripples can arise from a single initial seabed perturbation (Fig. 9).

Fig. 9. Numerical simulation of the development of a ripple field from of a single sandy Gaussian hump which interacts with an oscillatory flow. The dotted lines represent the location of ripple crests at successive times, while evolving towards the final ripple field (solid lines). Adapted from Marieu et al. (2008).

Wave ripples are induced by the horizontal component of the wave orbital velocity near the sediment bed. Ripple wavelengths can therefore not be much larger than the horizontal orbital excursion given by [math]U_0 T / \pi[/math], where [math]U_0[/math] is the maximum wave orbital velocity along the seabed and [math]T[/math] the wave period. Sheet flow transport prevails when wave orbital currents are strong ([math]U_0 \ge 0.5-1 \; m/s[/math]); bed ripples are wiped out. This explains why wave ripple lengths larger than 50 cm seldom occur. Megaripples with wavelength of several meters and complex 3-dimensional shapes are also observed in wave-dominated environments, but their generation mechanism is not yet well understood.

Sand waves in shelf seas are induced by tidal currents through essentially the same residual circulation process (Fig. 7) as wave ripples. The wavelength of sand waves is much smaller than the tidal velocity excursion, which can exceed 10 km in tidal seas. The weakness of the residual circulation cells associated with bed perturbations of great length is the main reason of the relatively small length of sand waves compared to the tidal excursion.

Related articles

Stability models

Rhythmic shoreline features

Wave ripples

Wave ripple formation

Sand ridges in shelf seas

Beach Cusps

Nearshore sandbars

Appendix: Ripple wavelength in steady flow (current ripples)

Fig. A1. The phase lead angle [math]\phi[/math] of the bed shear stress relative to the bed perturbation as a function of [math]k \, d_{50}[/math]. The blue dots are from electrochemical measurements. The solid red line is from a full numerical model. Figure adapted from Charru et al. (2013).

We use the following symbols:

[math]x = [/math] the coordinate along the main flow

[math]t =[/math] time

[math]\lambda =[/math] the wavelength of the bed perturbation, [math]k = 2 \pi /\lambda =[/math] the corresponding wavenumber

[math]z(x,t) = z_b e^{i k (x – ct) + \sigma t}[/math] the complex bed perturbation (only the real part counts), which is assumed periodic, with growth rate [math]\sigma[/math] and migration rate [math]c[/math]

[math]l_{sat} =[/math] the saturation length, of the order of 10-100 median grain diameters [math]d_{50}[/math]

[math]q_0 =[/math] the sediment flux over the flat bottom

[math]q_{sat} = C_{sat} \, k \, z(x,t) \, e^{i \phi } =[/math] the adjustment of the saturated sediment flux to the perturbed bottom ([math]C_{sat}=[/math] constant)

[math]\phi =[/math] the phase lead of the bed shear stress relative to the bed perturbation

[math]q(x,t) = q_0 + q’(x,t) = [/math] the true sediment flux, which adjusts to the saturated sediment flux according to the relaxation equation.

[math]\Large\frac{\partial q}{\partial x}\normalsize = i k q’ = \Large\frac{1}{l_{sat}}\normalsize \, (q_{sat} -q’) , \qquad[/math] giving [math]\qquad q’ = \Large\frac{q_{sat}}{1 + ik l_{sat}}\normalsize . \qquad (1)[/math]

The seabed evolution equation, [math]\quad \Large\frac{\partial z}{\partial t}\normalsize = - \Large\frac{\partial q}{\partial x}\normalsize \qquad[/math] gives [math]\qquad (\sigma - ikc) \, z = - ikq’ = - \Large\frac{ikq_{sat}}{1 + ikl_{sat}}\normalsize = - \Large\frac{ik^2 C_{sat} z e^{i \phi}}{1 + ikl_{sat}}\normalsize. \qquad (2)[/math]

Taking the real part of equation (2) gives the growth rate

[math]\sigma = \Large\frac{C_{sat}}{l_{sat}^2}\frac{y^2 \sin \phi \, – \, y^3 \cos \phi}{1 + y^2}\normalsize . \qquad (3)[/math]

Here [math]y = k l_{sat}[/math]. The growth rate is positive if [math]y \lt \tan \phi[/math].

The phase [math]\phi[/math] is a slowly varying function of the wavenumber [math]k \;[/math] for values of [math]k \, d_{50}[/math] in the range 10^-3-10^-2 (Fig. A1). We will ignore this variation when deriving an expression of the fastest growing wavelength. The fastest growing wavelength [math]k[/math] is then found by solving the equation [math]d \sigma /dk =0[/math]. The result is

[math]k = \Large\frac{1}{l_{sat}}\normalsize \Big( Y - \Large\frac{1}{Y}\normalsize \Big) , \quad Y = \Big( \Large\frac{1 + \sin \phi}{\cos \phi}\normalsize \Big)^{1/3} . \qquad (4)[/math]

Ripple development starts in laminar flow for perturbations with small wavelengths (large [math]k \, d_{50}[/math]), for which the phase [math]\phi \approx 45^o \;[/math]. Using [math]l_{sat} = 10-100 \: d_{50}[/math], the wavelength [math]\lambda[/math] of the fastest growing perturbation [math]\lambda = 100-1000 \: d_{50}[/math]. If downslope sediment transport (from crest to trough) is taken into account, a slope-dependent parameter must be subtracted from the phase [math]\phi[/math]. This implies that the wavelength of the fastest growing perturbations will be longer.

References

The main author of this article is Job Dronkers Please note that others may also have edited the contents of this article. Citation: Job Dronkers (2024): Principles of spontaneous bedform generation. Available from http://www.coastalwiki.org/wiki/Principles_of_spontaneous_bedform_generation [accessed on 14-09-2024] For other articles by this author see Category:Articles by Job Dronkers For an overview of contributions by this author see Special:Contributions/Dronkers J