Physics of Cohesive Sediment Flocculation and Transport: State-of-the-Art Experimental and Numerical Techniques

Due to climate change, sea level rise and anthropogenic development, coastal communities have been facing increasing threats from flooding, land loss, and deterioration of water quality, to name just a few. Most of these pressing problems are directly or indirectly associated with the transport of cohesive fine-grained sediments that form porous aggregates of particles, called flocs. Through their complex structures, flocs are vehicles for the transport of organic carbon, nutrients, and contaminants. Most coastal/estuarine models neglect the flocculation process, which poses a considerable limitation of their predictive capability. We describe a set of experimental and numerical tools that represent the state-of-the-art and can, if combined properly, yield answers to many of the aforementioned issues. In particular, we cover floc measurement techniques and strategies for grain-resolving simulations that can be used as an accurate and efficient means to generate highly-resolved data under idealized conditions. These data feed into continuum models in terms of population balance equations to describe the temporal evolution of flocs. The combined approach allows for a comprehensive investigation across the scales of individual particles, turbulence and the bottom boundary layer to gain a better understanding of the fundamental dynamics of flocculation and their impact on fine-grained sediment transport

Current- and wave- generated bedforms on mixed sand-clay intertidal flats : a new bedform phase diagram and implications for bed roughness and preservation potential

Funding: This work was supported by the United Kingdom’s Natural Environment Research Council (NERC) under Grant NE/I027223/1 (COHBED). JM, Julie Hope, and Daniel Parsons were partially funded by a Horizon 2020 European Research Council Consolidator Award (GEOSTICK, Grant 725955). The GEOSTICK project also kindly contributed the article processing fees. Andrew Manning’s contribution toward this research was partly supported by the National Science Foundation grants OCE-1924532 and OCE-1736668, TKI-MUSA project 11204950-000-ZKS-0002, and the HR Wallingford company research project FineScale (Grant Nos. ACK3013_62). DP received funding from the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland) and their support is gratefully acknowledged. MASTS is funded by the Scottish Funding Council (Grant reference HR09011) and contributing institutions.The effect of bedforms on frictional roughness felt by the overlying flow is crucial to the regional modelling of estuaries and coastal seas. Bedforms are also a key marker of palaeoenvironments. Experiments have shown that even modest biotic and abiotic cohesion in sand inhibits bedform formation, modifies bedform size, and slows bedform development, but this has rarely been tested in nature. The present study used a comprehensive dataset recorded over a complete spring–neap cycle on an intertidal flat to investigate bedform dynamics controlled by a wide range of wave and current conditions, including the effects of wave–current angle and bed cohesion. A detailed picture of different bedform types and their relationship to the flow, be they equilibrium, non-equilibrium, or relict, was produced, and captured in a phase diagram that integrates wave-dominated, current-dominated, and combined wave–current bedforms. This bedform phase diagram incorporates a substantially wider range of flow conditions than previous phase diagrams, including bedforms related to near-orthogonal wave–current angles, such as ladderback ripples. Comparison with laboratory-derived bedform phase diagrams indicates that washed-out ripples, lunate interference ripples and upper-stage plane beds replace the subaqueous dune field; such bedform distributions may be a key characteristic of intertidal flats. The field data also provide a means of predicting the dimensions of these bedforms, which can be transferred to other areas and grain sizes. We show that an equation for the prediction of equilibrium bedform size is sufficient to predict the roughness, even though the bedforms are highly variable in character and only in equilibrium with the flow for approximately half the time. Whilst the effect of cohesive clay is limited under more active spring conditions, clay does play a role in reducing the bedform dimensions under more quiescent neap conditions. We also investigated which combinations of waves, currents, and bed clay contents in the intertidal zone have the highest potential for bedform preservation in the geological record. This shows that combined wave–current bedforms have the lowest preservation potential and equilibrium current ripples have the highest preservation potential, even in the presence of moderate and storm waves. Hence, the absence of wave ripples and combined-flow bedforms and their primary stratification in sedimentary successions cannot be taken as evidence that waves were absent at the time of deposition.Publisher PDFPeer reviewe

Current- and Wave-Generated Bedforms on Mixed Sand–Clay Intertidal Flats:A New Bedform Phase Diagram and Implications for Bed Roughness and Preservation Potential

The effect of bedforms on frictional roughness felt by the overlying flow is crucial to the regional modelling of estuaries and coastal seas. Bedforms are also a key marker of palaeoenvironments. Experiments have shown that even modest biotic and abiotic cohesion in sand inhibits bedform formation, modifies bedform size, and slows bedform development, but this has rarely been tested in nature. The present study used a comprehensive dataset recorded over a complete spring–neap cycle on an intertidal flat to investigate bedform dynamics controlled by a wide range of wave and current conditions, including the effects of wave–current angle and bed cohesion. A detailed picture of different bedform types and their relationship to the flow, be they equilibrium, non-equilibrium, or relict, was produced, and captured in a phase diagram that integrates wave-dominated, current-dominated, and combined wave–current bedforms. This bedform phase diagram incorporates a substantially wider range of flow conditions than previous phase diagrams, including bedforms related to near-orthogonal wave–current angles, such as ladderback ripples. Comparison with laboratory-derived bedform phase diagrams indicates that washed-out ripples, lunate interference ripples and upper-stage plane beds replace the subaqueous dune field; such bedform distributions may be a key characteristic of intertidal flats. The field data also provide a means of predicting the dimensions of these bedforms, which can be transferred to other areas and grain sizes. We show that an equation for the prediction of equilibrium bedform size is sufficient to predict the roughness, even though the bedforms are highly variable in character and only in equilibrium with the flow for approximately half the time. Whilst the effect of cohesive clay is limited under more active spring conditions, clay does play a role in reducing the bedform dimensions under more quiescent neap conditions. We also investigated which combinations of waves, currents, and bed clay contents in the intertidal zone have the highest potential for bedform preservation in the geological record. This shows that combined wave–current bedforms have the lowest preservation potential and equilibrium current ripples have the highest preservation potential, even in the presence of moderate and storm waves. Hence, the absence of wave ripples and combined-flow bedforms and their primary stratification in sedimentary successions cannot be taken as evidence that waves were absent at the time of deposition

The role of biophysical cohesion on subaqueous bed form size

Biologically active, fine-grained sediment forms abundant sedimentary deposits on Earth's surface, and mixed mud-sand dominates many coasts, deltas, and estuaries. Our predictions of sediment transport and bed roughness in these environments presently rely on empirically based bed form predictors that are based exclusively on biologically inactive cohesionless silt, sand, and gravel. This approach underpins many paleoenvironmental reconstructions of sedimentary successions, which rely on analysis of cross-stratification and bounding surfaces produced by migrating bed forms. Here we present controlled laboratory experiments that identify and quantify the influence of physical and biological cohesion on equilibrium bed form morphology. The results show the profound influence of biological cohesion on bed form size and identify how cohesive bonding mechanisms in different sediment mixtures govern the relationships. The findings highlight that existing bed form predictors require reformulation for combined biophysical cohesive effects in order to improve morphodynamic model predictions and to enhance the interpretations of these environments in the geological record

Bedform migration in a mixed sand and cohesive clay intertidal environment and implications for bed material transport predictions

Many coastal and estuarine environments are dominated by mixtures of non-cohesive sand and cohesive mud. The migration rate of bedforms, such as ripples and dunes, in these environments is important in determining bed material transport rates to inform and assess numerical models of sediment transport and geomorphology. However, these models tend to ignore parameters describing the physical and biological cohesion (resulting from clay and extracellular polymeric substances, EPS) in natural mixed sediment, largely because of a scarcity of relevant laboratory and field data. To address this gap in knowledge, data were collected on intertidal flats over a spring-neap cycle to determine the bed material transport rates of bedforms in biologically-active mixed sand-mud. Bed cohesive composition changed from below 2 vol% up to 5.4 vol% cohesive clay, as the tide progressed from spring towards neap. The amount of EPS in the bed sediment was found to vary linearly with the clay content. Using multiple linear regression, the transport rate was found to depend on the Shields stress parameter and the bed cohesive clay content. The transport rates decreased with increasing cohesive clay and EPS content, when these contents were below 2.8 vol% and 0.05 wt%, respectively. Above these limits, bedform migration and bed material transport was not detectable by the instruments in the study area. These limits are consistent with recently conducted sand-clay and sand-EPS laboratory experiments on bedform development. This work has important implications for the circumstances under which existing sand-only bedform migration transport formulae may be applied in a mixed sand-clay environment, particularly as 2.8 vol% cohesive clay is well within the commonly adopted definition of “clean sand”

Magnetic and Geochemical Evidence of Yellow and Yangtze River Influence on Tidal Flat Deposits in Northern Jiangsu Plain, China

The formation of a broad tidal flat along the coast of the northern Jiangsu Province of China depends largely on the sediment supply from the Yellow and Yangtze Rivers although the relative contributions from each of these two large rivers remain uncertain. Knowledge of sediment sources to the tidal flat is critical for understanding the evolution of this muddy coast impacted by the large rivers. The present study focuses on tracking the sediment source of the tidal flat deposits based on sediment magnetic properties and geochemical analyses as well as statistical analysis. The study shows that sediments to the north of Yangkougang (site 9) have lower values of saturation isothermal remanence magnetization (SIRM), magnetic susceptibility (χ), demagnetization parameter S -100 ratios, lower Fe/Al values and higher Ca concentrations, while the opposite is true for samples south to Haozhigang (site 18). In addition, SIRM values of the \u3c16μm and \u3e63μm fractions generally display increasing trend from north to south. These results suggest that sediment sources rather than particle size variation are the dominant factors influencing the bulk magnetic properties. In light of the sediment composition comparison between the Yellow River and Yangtze River, factor analysis is used to identify sediment source of the tidal flat sediments. It is indicated that the Yellow River in its former course is the dominant supplier for sites to the north of Yangkougang (site 9), while the Yangtze River is the dominant supplier for sites to the south of Haozhigang (site 18). The coast between Yangkougang (site 9) and Haozhigang (site 18) is a transition zone influenced by both rivers. The inferred provenance contrasts are consistent with the pattern of coastal geomorphological evolution. Our data also suggest that sediments eroded from the former Yellow River delta have contributed to the evolution of the Yangtze Estuary in historical time. This study demonstrates that a combined magnetic and geochemical fingerprinting techniques in couple with statistical analysis is valuable for identifying sediment sources of tidal flats and deltas influenced by large rivers in the world

Oil-Mineral Flocculation and Settling Dynamics

In recent decades, oil spill contamination has tended to occur more commonly in deltaic and estuarial systems. The management of oil spillages has been a major challenge in the surrounding deltas due to the highly sensitivity nature of deltaic ecosystems. Many deltas have an abundance of clay minerals that can flocculate, and these play an important role in determining the transport of spilled oil contamination and its eventual fate, particularly given that suspended sediment and microbial activities are often prevalent and diverse in natural environments. The primary work presented here focuses on laboratory experimental studies that help develop improved parameterizations of flocculation processes for oil-sediment-biogeochemical modeling. Oil-mineral flocs (OMA) have been successfully created from a series of laboratory flocculation experiments. A floc video instrument LabSFLOC-2 has been adopted for the first time to study the settling dynamics of OMAs. Experimental results reveal OMAs can easily form in any oil, cohesive sediment, and seawater mixtures. However, Kaolin and Bentonite forms dramatically different OMA structures, which leads to their variable characteristics. In the Bentonite clay cases, the oil flocs tend to be much larger and with higher densities than those in Kaolin clay cases, resulting in significant variability of flocs settling velocities

Integrating field and laboratory approaches for ripple development in mixed sand–clay–EPS

The shape and size of sedimentary bedforms play a key role in the reconstruction of sedimentary processes in modern and ancient environments. Recent laboratory experiments have shown that bedforms in mixed sand–clay develop at a slower rate and often have smaller heights and wavelengths than equivalent bedforms in pure sand. This effect is generally attributed to cohesive forces that can be of physical origin, caused by electrostatic forces of attraction between clay minerals, and of biological origin, caused by ‘sticky’ extracellular polymeric substances (EPS) produced by micro-organisms, such as microalgae (microphytobenthos) and bacteria. The present study demonstrates, for the first time, that these laboratory experiments are a suitable analogue for current ripples formed by tidal currents on a natural mixed sand–mud–EPS intertidal flat in a macrotidal estuary. Integrated hydrodynamic and bed morphological measurements, collected during a spring tide under weak wave conditions near Hilbre Island (Dee Estuary, north-west England, UK), reveal a statistically significant decrease in current ripple wavelength for progressively higher bed mud and EPS contents, and a concurrent change from three-dimensional linguoid to two-dimensional straight-crested ripple planform morphology. These results agree well with observations in laboratory flumes, but the rate of decrease of ripple wavelength as mud content increased was found to be substantially greater for the field than the laboratory. Since the formation of ripples under natural conditions is inherently more complex than in the laboratory, four additional factors that might affect current ripple development in estuaries, but which were not accounted for in laboratory experiments, were explored. These were current forcing, clay type, pore water salinity and bed EPS content. These data illustrate that clay type alone cannot explain the difference in the rate of decrease in ripple wavelength, because the bed clay contents were too low for clay type to have had a measurable effect on bedform development. Accounting for the difference in current forcing between the field and experiments, and therefore the relative stage of development with respect to equilibrium ripples, increases the difference between the ripple wavelengths. The presence of strongly cohesive EPS in the current ripples on the natural intertidal flat might explain the majority of the difference in the rate of decrease in ripple wavelength between the field and the laboratory. The effect of pore water salinity on the rate of bedform development cannot be quantified at present, but salinity is postulated herein to have had a smaller influence on the ripple wavelength than bed EPS content. The common presence of clay and EPS in many aqueous sedimentary environments implies that a re-assessment of the role of current ripples and their primary current lamination in predicting and reconstructing flow regimes is necessary, and that models that are valid for pure sand are an inappropriate descriptor for more complex mixed sediment. This study proposes that this re-assessment is necessary at all bed clay contents above 3%