Dynamics of mud transport
Mud in coastal areas is mainly found in intertidal deposits .
Migniot  was probably the first to present a comprehensive overview of all the processes involved in mud dynamics.
In order to understand the dynamics of mud in coastal environments, it is necessary to properly define mud and its properties, in contrast to sand and other non-cohesive particles.
Mud is defined as a mixture of mainly fine-grained sediments (clays, silt and sand), organic matter and water , where the cohesive properties of the clay fraction, enhanced by the properties of the organic matter, dominate the overall behaviour. Studies on erosion behaviour of sand-mud mixtures indicate that the bed exhibits cohesive behaviour for clay contents above 15-20% .
In daily language "mud" refers to the deposited state of mud particles. In this state mud can occur as a fluid-like of soil-like entity. The dynamics of mud then refers to the formation, deformation and erosion of such layers.
A key feature of mud particles is their cohesive nature that distinguishes them from non-cohesive solid particles such as sand.
Clay particles in an aquatic environment tend to stick together (coagulate) as the result of the van der Waals forces into aggregates or flocs. This process is enhanced by slimes (EPS) and mucus produced by micro-benthos and bacteria (that feed on decaying organic matter).
The size, structure and density of flocs are determined by the forces the aggregate-particles undergo. These forces comprise: hydrodynamic forces (especially shear), collisions between particles and electrochemical forces (determined by the composition of solid particles and dissolved ions in the ambient water). The latter explains also why mud particles in fresh and saline water have a different structure.
The basic building blocs of flocs are the primary particles and/or flocculi (compact aggregates, O(10 µm), which rarely break down into primary particles), which form micro-flocs (silt size), macro-flocs (fine-sand size) and organic-rich mega-flocs (coarse-sand size, linked to seasonal biological events, such as algae-bloom). The simultaneous occurrence of micro- and macro-flocs is attributed to the tidal dynamics, and/or sometimes to the mixing of sea- and river-born aggregates. Further research is necessary to understand better the dynamics of the different floc populations.
Despite the trend to characterize the floc structure by a fractal number , the structure is not self-similar. In general, the fractal dimension decreases with increasing floc size, which implies that the floc structures becomes more and more open and its strength decreases.
|Schematic representation of mud floc dynamics.|
© Lee et al. (2012) .
Since size, structure and density of a mud particle vary dynamically, the same applies to its settling velocity. The value of the parameter “settling velocity” actually should correspond to the averaged value for the entire local floc population, i.e. the settling flux per unit concentration. It has been demonstrated that the sediment flux can much more accurately be calculated when considering the major floc populations individually .
For hindered settling of mud particles, it is important to express the hindrance correction factor in terms of the effective volumetric concentration occupied by the aggregates, including the immobilized water captured in the porous structure of the floc.
Mud Bed Formation
|Fluid mud at the Suriname Coast. |
© Hydraulics Laboratory, KULeuven
The critical volumetric concentration at which soil formation starts is called the gel point. In principle it could be predicted if the effective volume of the flocs is known. For North-Sea mud, the gel point is found for soil bulk densities of about 1100 kg/m³, which implies a solids concentration of the order of only 8%. This explains the apparent fluid-like behaviour of freshly deposited mud, since its structure is easily disturbed. This is better understood in the light of the rheological characteristics of fluid mud.
Due to flocculation and hindered settling, instant bed formation is retarded and a near-bottom high-concentrated (HC) suspension layer (also known as high-concentrated benthic suspension, HCBS) is formed.
Mud Bed Destruction
Considering the fact that the soil consists of aggregated particles, which are bound by electrochemical forces of varying origin and strength, these bonds may be broken by various mechanical forces, either shear forces (at the surface or internally) or excess pore pressures that exceed the effective stress (usually due to wave action, but also by e.g. earthquakes). Local micro-cracks may grow to larger cracks and eventually create failure planes in the bed. Erosion is the process where the structure is broken to such degree that the loose parts may be picked-up by the flow and transported. A distinction is furthermore made into the following three erosion modes:
- Surface erosion = particle-by-particle erosion at the surface;
- Mass or Bulk erosion = erosion of a patch of mud above a failure plane;
- Cliff erosion = break-up of solid clumps of over consolidated mud by (boat) waves, which after mechanical rolling erosion may become rounded mud pellets.
|Cliff erosion of an intertidal mudflat in the Western Scheldt Estuary.|
© Michiel Taal / Deltares
The deformation and flow of fluid mud requires a rheological description, i.e. mathematical expressions (closure laws) relating stress to deformation and/or deformation rate.
Deformation under wave action has traditionally been described by visco-elastic or plastic models, which allow analytical solutions under (over)simplified (1D) conditions.
However, considering the link between strength and bed structure, a description by a thixotropic characterization seems more promising. This approach shows (logical) analogies to the structural kinetics approach to flocculation .
Suspended Load Transport of Mud
Traditionally, it has been assumed that cohesive sediments have such low settling velocities that the dominant transport mode is by suspended load. Usually, only dilute suspension transport is considered.The importance of high-concentrated (HC) suspension transport in the inner layer above the bed (often named “fluid mud”, but this term is more consistently restricted to another state > see Fluid Mud) is often underestimated or ignored. However, the amount of sediment transported in this layer can be very significant. Research on this topic is still ongoing.
The thickness of the HC suspension layer above the bottom can be significantly larger than in the case of sand. A relatively sharp interface, a lutocline, can be found between this layer and the dilute layer above. Instabilities can be observed along this lutocline in the form of internal waves.
Contrary to HC sand suspension layers, HC mud suspension layers usually exhibit strong turbulence damping (or even laminarization) and drag reduction. A well known example is that of the Yellow River (China), where roughness values have to be taken corresponding to smoother than a smooth glass plate, in order to predict the hydrodynamic resistance correctly.
Due to interaction of heavier saline sea water and lighter fresh river water, a zone with high suspended loads of cohesive sediments, an estuarine turbidity maximum (ETM), is formed in estuaries. The turbulence due to the interaction keeps particles in suspension, while large fresh water flocs may alter into more compact flocs with a lower settling rate due to the increased salinity. The ETM location moves up- and downstream with the tide, depending on the strength of the tide and the river discharge. In some estuaries, more than one ETM can occur.
Bedload Transport of Mud
Despite the traditional approach of only considering suspended transport for mud, two observations suggest that bedload transport of mud should not too easily be discarded.
- High-concentrated near-bottom transport of mud should not be treated as suspended load, since the theory for the latter requires fully-developed turbulent flow, while the former typically shows the properties of low-Reynolds-number flow, requiring another mathematical treatment (which is still under development). Because of the high amounts of sediment transported in this layer, and the physical similarities with bedload transport of non-cohesive sediments, it seems justified to consider this mode as bedload transport of mud.
- Large clumps of overconsolidated may break off from mud layers that are exposed to the air for longer time such that they dry and show cracks. Once the water line reaches these areas in combination with significant wave action (either wind waves or boat traffic induced waves), large clumps may break off, resulting in cliff erosion. Pulled up and down by the tidal currents, they role over the bottom, they may break up into smaller clumps, erode more or less, and eventually become rounded mud pebbles which can role even better and be transported over longer distances or accumulate in deeper areas. For example, mud pebbles are formed along the low water line of the intertidal flats of the Yzer Mouth (Nieuwpoort, Belgium) and found back in the mud dredged from the navigation channel. – Laboratory observations suggest that bulk erosion may also generate submarine mud pebbles. As far as known, no evidence from the field on the latter is available.
|Mud pebble formation after cliff erosion along the Yzer mouth intertidal mudflat.|
© Erik Toorman / KULeuven
Fluid Mud Flow & Density Currents
Mud bottoms in underconsolidated state are prone to fluidization and/or liquefaction. The resulting fluid mud can deform under wave action and, when there is no barrier, even flow driven by gravity.
As long as the density remains below the gel point, the effective viscosity is so high that the flow behaviour will remain laminar. When accelerating, the shear and resulting instabilities at the interface will cause entrainment in the two directions (i.e. water into the fluid mud and mud particles into the water column), resulting in dilution of the fluid mud below the gel point, and the water-sediment mixture should now be considered a (highly-concentrated) suspension, which can flow turbulently.
When depositing on a sloping bottom (e.g. river banks on dredged channel slopes), self-weight induced shear may avoid consolidation and keep the deposit liquefied, such that the deposit may flow as a gravity current to the lowest point.
This knowledge is used in low-cost agitation dredging, where a mud deposit is mechanically disturbed and liquefied in low-energy locations, such as lock entrances, with the purpose to accumulate the fluid mud in deeper areas where the mud, if necessary, can be removed during maintenance dredging.
Mud gravity currents (or avalanches) may also be induced due to liquefaction by earth-quakes.
Mud gravity currents may carry large amounts of mud to the deeper ocean where rivers end in a canyon. Example: It has been estimated that about 50% of the mud from the Amazon flows down the continental slope in front of its mouth. About 20% is transported to the west by wave action as migrating mud banks along the Guyanas coast and eventually flows down the Orinoco river canyon into the deeper ocean.
|Mud density currents generated in laboratory flume.|
© Erik Toorman / KULeuven
An important property of fluid mud layers is their capacity to absorb energy from surface waves in the overflowing water layer. Famous examples are the Guyanas coast  and the Lousiana coast  , where this has been studied.
It has been demonstrated that wave damping can well be simulated using a thixotropic closure .
The well known nearshore spectral wave model SWAN has also been extended with a semi-empirical sink term to account for wave damping by fluid mud .
Mud transport modelling
Cohesive sediment transport models for short to medium term predictions of morphological (trend) studies combine the hydrodynamic equations for the sediment-water mixture with a sediment mass balance (or transport) equation. The latter is an advection-diffusion equation with sink (deposition) and source (erosion) terms to allow exchange of sediment fluxes with the bottom.
The most famous surface erosion law for cohesive sediments has been named after Partheniades’ experimental work . The erosion flux is the product of an erosion rate multiplied with a probability factor as a function of the shear stress in excess of a critical erosion shear stress. Subsequent research has found that modifications had to allowed in the case of soft, freshly deposited mud . Partheniades’ law has recently been extended to account for the statistical turbulent fluctuations on the hydrodynamic shear stress .
Deposition laws (2D vs. 3D)There still is much confusion on applying the correct deposition law for (cohesive) sediments. The exact bottom boundary condition at the bottom is simply imposing the sedimentation flux (the product of sediment concentration and settling velocity) at the bottom. This has to be used in 3D and 1- and 2DV models where the vertical water column is resolved.
The famous empirical deposition law from Krone  may only be used in the case of depth-averaged (2DH) modelling. Here the bottom flux is approximated by the product of depth-averaged sediment concentration with the settling velocity of a single particle, multiplied with a probability p of the fraction of the total suspended load which may deposit. i.e., the complementary probability actually represents the amount of sediment which can be held in suspension by the turbulent energy in the flow. Introducing this probability compensates for the fact the no vertical energy balance is computed in 2DH models.
Since the settling velocity of the dynamically varying mud aggregate particles changes as the floc structure changes, many attempts have been made to account for flocculation. The first simple models were empirical relationship correcting the settling velocity with a factor as a function of instantaneous concentration and turbulent dissipation rate. The next generation introduced a kinetic equation describing the aggregation and break-up of flocs . Following, more elaborate population balance equation modelling has been introduced, applied to multiple floc size fractions. In order to reduce the computational effort for large scale applications, new strategies have been proposed to work with a continuous particle size distribution  and to distinguish at least two floc populations  .
Sand-mud mixture modelling
In most coastal environments mud and sand coexist and influence each others’ behaviour. Because cohesive sediment was thought to be transported predominantly as suspended load and sand as bed load, the two were often treated independently. Current understanding of sediment transport strongly suggests that the two should be solved in a coupled way. The major reason is the interaction in the inner layer above the bottom, where there are high concentrations and four-way coupling fluid-particle interactions. Furthermore, the sand-mud ratio in the bed determine the availability of each fraction and the erosion behaviour of the surface.
To date, the most advanced modelling approaches have been proposed by van Ledden  and Waeles .
In order to account for dynamical erosion properties of the bed surface, several models have been proposed to predict the erosion resistance. The most detailed models rely on the geotechnical self-weight consolidation theory of Gibson .
Because of the computational effort and the high uncertainty on the model parameters, Sanford  proposed a much simpler relaxation model.
The latter has been investigated to be extended with a fluidisation model, to account for the counter action of wave induced pore pressure build up.
In all these models the bed is subdived into layers. Exchange with the water column is allow with the surface layer only (the “active” layer). A practical problem remains the book-keeping of the thickness of each layer and of the sand-mud ratio, when the model takes into account sediment mixtures.
The importance of organic matter in mud has been acknowledged for a long time. They change time scales of certain processes (such as consolidation) and strength parameter by one order of magnitude.
Unfortunately, the characterization and quantification of the effects of organic matter is very complex because the many biological agents and seasonal effects.
Benthic organisms tend to modify cohesive sediment characteristics . Individual particles are ingested and form large, low density faecal pellets with settling velocities much higher than the individual particles contained in them. Estuarine sediment transport conditions are altered e.g. by change of bed roughness due to deposition of fecal pellets or by biological activity in the surface layers causing either binding or destabilizing of the bed surface .
Density, size and settling velocity of fecal pellets are important in describing sediment transport processes . Faecal pellets are low density faeces from suspension-feeding and mud-ingesting marine and estuarine benthic animals including oysters, worms, and barnacles. They consist of large amounts of fine sediments mainly single particles in the 1-5 μm range  and flocs ingested by the feeding animal, compacted in the gut and excreted as large aggregates glued together by mucus. Up to 25 % of the deposited sediment may in some areas originate from biodeposition.
Decay of organic matter generates methane gas in the pores of the soil skeleton. This results in changes in the mechanical properties of the mud layer .
Major application fields
Morphodynamics & Sediment transportErosion and/or deposition change the boundaries between land and water. Human interventions into the natural environment disturb the natural trend of streams to reach their regime. The impact of anthropogenic changes on the environment needs to be assessed and, if possible, quantified through model studies. The interest goes to the prediction of sediment fluxes, sediment balances and morphological trends, necessary to assess their impact on economy (navigability), safety (against flooding) and ecology.
Because of the large complexity and the multi-scale processes involved, the required modelling technique will depend very much on the purpose of the study and the spatial and temporal scales involved.
Dredging & Nautical Depth
Navigation channels, lock entrances and harbour docks are prone to accumulation of mud on the bottom. In order to maintain navigation, these areas must frequently be dredged.
The design of dredging equipment, the estimation of dredging operation costs and the environmental impact assessment of dredging activities requires the understanding of hydraulic transport of mud and sand-mud mixtures.
The least studied problem related to the occurrence of mud in navigation channels, is the determination of a practical and meaningful criterion to determine when a fluid mud bottom becomes harmful to navigation. This is of importance to optimize the economical cost of the maintenance of navigation ways and harbours.
At present, the criterion is still based on a critical bulk density of fluid mud, since it can relatively easily be determined and there is a clear relationship between bulk density and rheological parameters. The traditional density criterion, however, overlooks the time-dependence of the rheological parameters (i.e. thixotropy). For this reason, further research is required to establish a better criterion which accounts the latter. For this purpose a research programme has started at Flanders Hydraulics (Antwerp, Belgium) .
General books on cohesive sediments (mud):
- Partheniades, E. (2009). “Cohesive sediment in open channels”. Butterworth-Heinemann, Burlington (MA), xv + 358 pp.
- Whitehouse, R., R. Soulsby, W. Roberts & H. Mitchener (2000). “Dynamics of estuarine muds”. HR Wallingford & Thomas Telford, London (UK), xxii + 210 pp.
- Winterwerp, J.C. & W. van Kesteren (2004). “Introduction to the physics of cohesive sediment in the marine environment”, Elsevier, Amsterdam (NL), xiii + 466 pp. + annexes.
- Coast erosion
- Types and background of coastal erosion
- Protection against coastal erosion
- Manual Sediment Transport Measurements in Rivers, Estuaries and Coastal Seas
- Eisma, D. et al. (1997). Intertidal Deposits: River Mouths, Tidal Flats and Coastal Lagoons. CRC Press, Boca Raton (FL), 525 pp.
- Migniot, C. (1968). Etude des propriétés physiques de différents sediments très fins et leur comportement sous des actions hydrodynamiques. La Houille Blanche, 1968(No.7):591-620 (in French).
- Berlamont, J., Ockenden, M., Toorman, E. & Winterwerp, J. (1993). The characterisation of cohesive sediment properties. Coastal Engineering, 21:105-128.
- Mitchener, H. & Torfs, H. (1996). Erosion of mud/sand mixtures. J. Coastal Engineering, 29: 1-25.
- Kranenburg, C. (1994). On the fractal structure of cohesive sediment aggregates. Estuarine, Coastal and Shelf Science, 39:451-460.
- Lee, B.J., Fettweis, M., Toorman, E., Moltz, F. (2012). Multimodality of a particle size distribution of cohesive suspended particulate matters in a coastal zone. J. Geophysical Research, 117, C03014, 17pp. (doi:10.1029/2011JC007552)
- Toorman, E.A. (1997). Modelling the thixotropic behaviour of dense cohesive sediment suspensions. Rheologica Acta Vol.36 (No.1):56-65.
- Wells, J.T. & Coleman, J.M. (1981). Physical processes and fine-grained sediment dynamics, coast of Surinam, South-America. J. Sedimentary Petrology, 51(4):1053-1068.
- Sheremet, A., Jaramillo, S., Su, S.-F., Allison, M.A. & Holland, K.T. (2011). Wave-mud interaction over the muddy Atchafalaya subaqueous clinoform, Louisiana, United States: wave processes. J. Geophysical Research, 116, C06005, 14 pp. (doi:10.1029/2010JC006644).
- Toorman, E.A. (2008). An investigation into the thixotropic wave dissipation potential of fluid mud. AGU Chapman Conference on Physics of Wave-Mud Interaction (Amelia Island, Florida, November 2008). Book of Abstracts, p.24.
- Kranenburg, W., Winterwerp, J., de Boer, G., Cornelisse, J., & Zijlema, M. (2011). SWAN-Mud: engineering model for mud-induced wave damping. J. Hydraul. Eng., 137(9), 959–975.
- Partheniades, E. (1965). Erosion and deposition of cohesive soils. J. Hydraulic Division ASCE, 91(HY1) :105-139.
- Parchure, T.M. & Mehta, A.J. (1985). Erosion of soft cohesive sediment deposits. J. Hydraulic Engineering, 111 (10) :1308-1326.
- Van Prooijen, B.C. & Winterwerp, J.C. (2010). A stochastic formulation for erosion of cohesive sediments. J. Geophysical Research, 115, C01005, 15 pp. (doi:10.1029/2008JC005189)
- Krone, R.B. (1962). Flume studies of the transport of sediment in estuarial shoaling processes. Final Report; Hydraulic Engineering Laboratory and Sanitary Research Engineering Laboratory, University of California, Berkeley (CA).
- Winterwerp, J.C. (2002). On the flocculation and settling velocity of estuarine mud. Continental Shelf Research, 22:1339-1360.
- Maerz, J., Verney, R. , Wirtz, K. & Feudel, U. (2011). Modeling flocculation processes: Intercomparison of a size class-based model and a distribution-based model. Continental Shelf Research, 31:S84–S93.
- Lee, B.J., E.A. Toorman, F. Moltz & J. Wang (2011). A two-class population balance equation yielding bimodal flocculation of marine or estuarine sediments. Water Research, 45:2131-2145.
- Van Ledden, M., van Kesteren, W.G.M. & Winterwerp, J.C. (2004). A conceptual framework for the erosion behaviour of sand–mud mixtures. Continental Shelf Research, 24:1–11.
- Le Hir, P., Cayocca, F. & Waeles, B. (2011). Dynamics of sand and mud mixtures: a multiprocess-based modelling strategy. Continental Shelf Research, 31:S135–S149.
- Gibson, R. E., England, G. L., & Hussey, M. J. L. (1967). The theory of one-dimensional consolidation of saturated clays: I. Finite nonlinear consolidation of thin homogeneous layers. Geotechnique, 17:261–273. &
Gibson, R. E., Schiffman, R. L., & Cargill, K. W. (1981). The theory of one-dimensional consolidation of saturated clays: II. Finite non-linear consolidation of thick homogeneous layers. Can. Geotech. J., 18:280–293.
- Sanford, L. (2008). modeling a dynamically varying mixed sediment bed with erosion, deposition, bioturbation, consolidation and armoring. Computers and Geosciences, 34:1263-1283.
- Edelvang, K. 1995. "The significance of aggregation in an estuarine environment." PhD thesis printed in Geographica Hafniensia A5 ISBN 87-87945-18-5
- Nowell ARM, Jumars PA & Eckman JE (1981) Effects of biological activity on the entrainment of marine sediments. Mar Geol 42, p 133-153
- Taghon GL, Nowell ARM & Jumars PA (1984) Transport and breakdown of faecal pellets: Biological and sedimentological consequences. Limnol Oceanogr Vol 29(1), p 64-72
- Haven DS & Morales-Alamo R (1972) Biodeposition as a factor in sedimentation of fine suspended solids in estuaries. Geol Soc Mem, Vol 133, p 121-130
- van Kesteren, W. & van Kessel Th. (2002). Gas bubble nucleation and growth in cohesive sediments. Proceedings in Marine Science, Vol.5:329–341, Elsevier, Amsterdam (NL).
- Claeys, S. (2011). Rheology as a survey tool. Hydro International, 15:3
Please note that others may also have edited the contents of this article.