Skip to main content
Geosciences LibreTexts

6.8.2: Critical shear stress and settling velocity

  • Page ID
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    Sand, silt and clay particles are continuously entrained, transported and deposited in the coastal system by currents and waves. Due to the difference in grain size and cohesiveness, the sediment transport characteristics of sand and mud are different. In general, sand is transported as bed load (at the bed-water interface) and as suspended load (higher in the water column). For low sediment concentrations, silt and clay are only transported in suspension. For high concentrations, silt and clay can also be transported as fluid mud in a layer near the bed. In that case, the fluid mud will act as a viscous layer on top of the bed, blocking exchange of sand between the bed and the water column.

    We have seen that the fall velocity is an important parameter determining the vertical distribution of suspended sediment (see for instance Eq. The fall velocity depends on the size, shape and specific density of the particles and of the water temperature and the sediment concentration. As a reference, an individual silt particle of \(10\ \mu m\) has a fall velocity in still water of about \(0.1\ mm/s\). For sand, the particle properties are relatively independent of time (because the quartz crystals are chemically hard and stable, it takes years to change the size of a sand grain). But, as we have seen in Sect. 2.6.2, the size and shape of mud can change quite easily. As a result of flocculation, the size of individual mud particles increases. With this increase, particles become heavier, which will speed up the settling of the particles. On the other hand, due to the increase in size, also the drag force, the resistance the particle ‘feels’ when settling, will increase. This will slow down the particle and can eventually cause deflocculation.

    截屏2021-11-03 下午10.52.44.png
    Figure 6.17: Effect of mud content on the critical erosion shear stress for two different sand-mud mixtures (data from Torfs, 1995). Mud content on horizontal axis represents % fines smaller than \(63\ \mu m\), hence including silt and clay fractions.

    Another factor influencing the sediment transport is the exchange between the bottom layer and the water column. The upward sediment flux from the bed into the water column strongly depends on the bed composition. From laboratory experiments it was concluded that two regimes can be distinguished: a non-cohesive and a cohesive regime (Van Ledden & Wang, 2001). The clay content (the weight percentage of particles with a grain size smaller than \(4\ \mu m\)) of the bed material is the governing parameter for the transition between both regimes. For clay content of less than 5% to 10% the bed behaves more or less non-cohesively. Sand and mud are eroded more or less independently, so the ‘normal’ formulations for the bed erosion can be used. For clay content exceeding 5% to 10% the bed behaves cohesively. The sand and mud particles are eroded simultaneously. Figure 6.17 gives the critical velocities as a function of mud (clay and silt) content in the bottom material.

    In sediment transport modelling, the exchange processes between water column and bottom layer are often assumed to be governed by the critical bed shear stress of the bottom material and the actual occurring bed shear stresses (Sect. 6.3).

    For actual bed shear stresses larger than the critical bed shear stress of the bed material, the grains can be eroded, whereas for actual bed shear stresses smaller than the critical bed shear stress of the bed material, the grains can settle:

    \[\tau_{b, \text{actual}} > \tau_{b, cr} \to \text{erosion}\]

    \[\tau_{b, \text{actual}} < \tau_{b, cr} \to \text{sedimentation}\]

    For non-cohesive material the critical bed shear stress \(\tau_{b, cr}\) for erosion is equal to the critical bed shear stress for sedimentation. For cohesive material, where there is a strong binding force between the particles, this is not the case. The critical bed shear stress for erosion is larger than the critical bed shear stress for sedimentation. This implies that for a range of actual bed shear stresses, there is no exchange with the bottom layer in the case of cohesive sediments:

    \[\tau_{b, \text{actual}} > \tau_{b, cr, er} \ \ \ \ \to \text{erosion}\]

    \[\tau_{b, cr, sed} < \tau_{b, \text{actual}} < \tau_{b, cr, er} \ \ \to \text{no exchange with bottom layer}\]

    \[\tau_{b, \text{actual}} < \tau_{b, cr, sed}\ \ \ \ \to \text{sedimentation}\]

    The exact formulations for sediment transport of mixtures of sand and mud are still a subject of research.

    This page titled 6.8.2: Critical shear stress and settling velocity is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Judith Bosboom & Marcel J.F. Stive (TU Delft Open) via source content that was edited to the style and standards of the LibreTexts platform.