Modeling flocculation processes of fine-grained particles using a size-resolved method: Comparison with published laboratory experiments

The transport of fine-grained particles in estuarine and coastal waters is influenced by flocculation processes (aggregation and floc breakup). As a consequence, the particle size varies with time in the water column, and can be orders of magnitude larger than those of primary particles. In this study the variations in floc size is simulated using a size-resolved method, which approximates the real size distribution of particles by a range of size bins and solves a mass balance equation for each bin. To predict the size distribution both aggregation and breakup processes are included. The conventional rectilinear aggregation kernel is used which considers both turbulent shear and differential settling. The breakup kernel accounts for the fractal dimension of the flocs. A flocculation simulation is compared to the settling column lab experiments of Winterwerp [1998. A simple model for turbulence induced flocculation of cohesive sediment, Journal of Hydraulic Research, 36, 309–326], and a one-dimensional sediment transport model is verified with the observed variations in floc size and concentration over tidal cycles in a laboratory flume experiment of Bale et al. [2002. Direct observation of the formation and break-up of aggregates in an annular flume using laser reflectance particle sizing. In: Winterwerp, J.C., Kranenburg, C. (Eds.), Fine Sediment Dynamics in the Marine Environment. Elsevier, pp. 189–201]. The numerical simulations compare qualitatively and quantitatively well with the laboratory measurements, and the analysis of the two simulation results indicates that the median floc size can be correlated to the sediment concentration and Kolmogorov microscale. Sensitivity studies are conducted to explore the role of settling velocity and erosion rate. The results are not sensitive towards the formulation of settling velocity, but the parameterization of erosion flux is important. The studies show that for predicting the sediment deposition flux it is crucial to include flocculation processes.     

An idealized model study of flocculation on sediment trapping in an estuarine turbidity maximum

Flocculation has an important impact on particle trapping in estuarine turbidity maximum (ETM) through associated increases in particle settling velocity. To quantify the importance of the flocculation processes, a size-resolved flocculation model is implemented into an ocean circulation model to simulate fine-grained particle trapping in an ETM. The model resolves the particle size from robust small flocs, about 30 μm, to very large flocs, over 1000 μm. An idealized two-dimensional model study is performed to simulate along-channel variations of suspended sediment concentrations driven by gravitational circulation and tidal currents. The results indicate that the flocculation processes play a key role in generating strong tidal asymmetrical variations in suspended sediment concentration and particle trapping. Comparison with observations suggests that the flocculation model produces realistic characteristics of an ETM.     

A simple numerical model for inlet sedimentation at intermittently open–closed coastal lagoons

This contribution aims to model the dominant processes that control sedimentation within the ocean inlet to intermittently open–closed coastal lagoons; focussing on the role of infilling, backfilling and morphodynamic feedback. The key elements that have been included in the present model are: (1) the delivery of sediment to the mouth of the estuary by littoral processes; (2) sediment transport processes within the inlet due to non-linear tidal flow; and (3) the down-slope diffusion of sediment. The model discussed here includes a simple one-dimensional (1-D), non-linear flow parameterisation that predicts the spatial variability in the magnitude and non-linearity of the tidal flow. The predicted third and fourth velocity moments are used to drive a morphodynamic module. Down-slope diffusion of sediment is dealt with in a separate diffusion term in the bed-evolution equation. Feedback between the evolving morphology and the flow field are integral to the model. Numerical simulations are used to investigate different modes of evolution for this type of inlet system. Inlet closure due to infilling is critically controlled by the balance between sediment addition at the mouth due to littoral processes and the removal of sediment by non-linear tidal flow. Rapid widening or deepening of the inlet at its landward margin leads to the deposition of a flood shoal. Under conditions of high sedimentation (and low diffusion) the flood shoal can become sufficiently well developed to present a physical barrier to sediment entering the lagoon. Under these circumstances backfilling can become significant. The infilling and backfilling processes are ameliorated by efficient down-slope diffusion which is controlled in the present model by a diffusion parameter, D. High diffusion coefficients slow inlet closure and allow sediment to propagate further into the lagoon.     

Freshwater flocculation of suspended sediments in the Yangtze River, China

This study focuses on suspended sediments and in situ flocculation in the Yangtze River, with the goal of improving our understanding of the relationship between freshwater and estuarine flocculation. A field survey with state-of-the-art instruments was carried out in January 2008 in the reach from downstream of the Three Gorges Dam to the estuary. The data show that in situ floc mean diameters range from 22 to 182???m in the river, whereas the median dispersed grain sizes are 4.4?C11.4???m. This demonstrates that flocculation is an important process during the transport of suspended sediments along the river. The flocculation characteristics, suspended sediment concentration and dispersed grain sizes all vary longitudinally in the main stream of the Yangtze River. Biochemical factors are likely be more significant in the freshwater flocculation than in the estuary, where hydrodynamics and biochemical factors are both important. Flocculation is found in the freshwater river, in the estuary and in coastal waters, which indicates that dynamic break-up/reflocculation processes take place during the suspended sediment transport. The freshwater flocs may behave as parent flocs to the estuarine flocculation. This study enhances our understanding of flocculation from estuarine and coastal areas to fresh river systems and provides insights into the effects of input of riverine flocs to the estuarine flocculation and into the sources and fate of flocs.     

Floc size and settling velocity within a Spartina anglica canopy

There is increasing interest in tidal wetlands as mechanisms for sustainable and long-term coastal defence. The complexities of the interaction between the deposition of suspended particulate matter (SPM) and submerged vegetation, however, is to a large extent poorly understood. Consequently, accurate parameterisation of cohesive sediment settling fluxes in these environments is a crucial requirement for the development of high-resolution numerical models of wetland morphodynamics. A novel laboratory experiment is described in which the turbulent flow structure within a canopy of the halophytic macrophyte Spartina anglica is examined, and floc characteristics quantified using a unique floc camera configuration able to measure directly the full spectral floc size (D) and settling velocity (W_s). We provide the first quantitative observations of floc characteristics from shallow (h<0.5 m), vegetated flows and investigate the potential influence that variations in vegetative density may have on flocculation, and thus depositional fluxes, in comparison to unvegetated flows.     

Lattice Boltzmann simulation of turbulence-induced flocculation of cohesive sediment

Both the floc formation and floc breakup of cohesive sediment are affected by turbulent shear which is recognized as one of the most important parameters, and thus, on the settling and transport of cohesive sediment. In this study, the development of floc characteristics at early stage and steady-state of flocculation were investigated via a three-dimensional lattice Boltzmann numerical model for turbulence-induced flocculation. Simulations for collision and aggregation of various size particles, floc growth, and breakup in isotropic and homogenous turbulent flows with different shear stresses were conducted. Model results for the temporal evolution of floc size distribution show that the normalized floc size distributions is time-independent during early stage of flocculation, and at steady-state, shear rate has no effect on the shape of normalized floc size distribution. Furthermore, the size, settling velocity, and effective density of flocs at the non-equilibrium flocculation stage do not change significantly for shear stresses in the range 0–0.4 N m^?2. The relationships between floc size and settling velocity established during floc growth stages and that during steady-states are different.     

Controls on effective settling velocities of mud flocs in the Changjiang Estuary,China

Theoretical and experimental studies were undertaken to gain insight into physical parameters controlling the flocculation and settling properties of mud flocs in the Changjiang Estuary, China. The Rouse equation is applied to vertical profiles of suspended sediment concentration to determine the bulk mean settling velocity (w_s) of sediment suspended in the Changjiang Estuary. Both in situ point‐sampled and acoustically measured profiles of suspended mud concentrations were fit selectively. The calculated settling velocities w_s mainly ranged from 0·4 to 4·1 mm s^?1 for the point‐sampled data set, and from 1·0 to 3·0 mm s^?1 for the acoustically measured data set. Furthermore, the settling velocities of mud flocs increased with mean concentration (C?) of mud flocs in suspension and were proportional to increasing bottom shear stress (τ_b) of tidal flow. The best equation for the field settling velocity of mud flocs in the Changjiang Estuary can be expressed by the power law: w_s = mC?ⁿ (m, 1·14–2·37; n, 0·84–1·03). It is suggested that C? and τ_b were the dominant physical parameters controlling the flocculation and w_s of mud flocs in suspension. Copyright © 2004 John Wiley & Sons, Ltd.     

Modeling time dependence in a tidal boundary layer

Although the rate of acceleration and deceleration in a depth-limited tidal flow can be quite large, it has been found that quasi-steady logarithmic models or steady eddy diffusivitics can explain the turbulent field quite well (Gross andNowell, 1983;Long, 1981). A simple, time-dependent model was devised based on kinetic energy closure (Bradshaw et al., 1967) to simulate a depth-limited tidal boundary layer without the assumption of a time-independent eddy diffusivity. Computer solutions revealed that the parameterS = u_s/δΩ (u_s is a turbulent velocity scale, δ is depth, Ω is tidal frequency) governs the phasing of the velocity to the forcing pressure gradient. However, for values ofS reasonable for depth-limited estuarine tidal flow the ratio of the turbulent energy production-dissipation balance to rate of change of turbulent energy is always quite large, indicating that the turbulent adjustment rate is fast enough for the flow to be modeled with quasi-steady techniques.     

Flocculation and the loss of sediment from the Po River plume

In October 2000, a 100-year flood event in the Po River resulted in the formation of a fine-grained sediment deposit extending up to 10 km from the river mouth. Soon after this event, and for a subsequent period of 2 years, box cores were collected on a grid of stations off of the Po Delta to observe the evolution of the flood deposit. Using a process-based parameterization of the disaggregated inorganic grain size distribution, the evolution of the surficial sediment on the Po shelf since the 2000 flood has been interpreted in the context of particle flocculation dynamics. This method produces estimates of floc limit, the diameter at which the flux of single grains to the seabed equals the flux of flocs, and floc mass fraction, the amount of material deposited as flocs. Floc limit depends on the extent of flocculation in suspension, and floc mass fraction describes the extent of flocculation in the sediment. Immediately after the flood, these two parameters were high at stations located beneath the path of the flood plume as observed in satellite images. The occurrence of a highly flocculated deposit below the path of the plume leads to two hypotheses: (1) the high sediment concentration in the river flood plume produced extensive sediment flocculation in the plume and (2) post-depositional remobilization of sediment delivered from the plume to the seabed was limited in the relatively low-energy environment of the Po prodelta. Floc limit and floc mass fraction estimated from bottom sediment sampling 3 and 10 months after the initial sampling were lower, indicating that during normal discharge, flocculated fine-grained sediment from the Po River settles close to the mouth, leaving only a small amount of material in suspension in the plume for direct deposition onto the prodelta. These findings are consistent with laboratory studies of suspended sediment that show that sediment concentration and turbulent energy exert dominant control on the extent of flocculation and the loss of sediment from suspension.     

Fine sediment transport by tidal asymmetry in the high-concentrated Ems River: indications for a regime shift in response to channel deepening

This paper describes an analysis of the observed up-river transport of fine sediments in the Ems River, Germany/Netherlands, using a 1DV POINT MODEL, accounting for turbulence-induced flocculation and sediment-induced buoyancy destruction. From this analysis, it is inferred that the net up-river transport is mainly due to an asymmetry in vertical mixing, often referred to as internal tidal asymmetry. It is argued that the large stratification observed during ebb should be attributed to a profound interaction between turbulence-induced flocculation and sediment-induced buoyancy destruction, as a result of which the river became an efficient trap for fine suspended sediment. Moreover, an asymmetry in flocculation processes was found, such that during flood relative large flocs are transported at relative large flow velocity high in the water column, whereas during ebb, the larger flocs are transported at smaller velocities close to the bed??this asymmetry contributes to the large trapping mentioned above. The internal tidal asymmetry and asymmetry in flocculation processes are both driven by the pronounced asymmetry in flow velocities, with flood velocities almost twice the ebb values. It is further argued that this efficient trapping is the result of a continuous deepening of the river, and occurs when concentrations in the river become typically a few hundred mg/l; this was the case during the 1990 survey analyzed in this paper. We also speculate that a second regime shift did occur in the river when fluid mud layers become so thick that net transport rates are directly related to the asymmetry in flow velocity itself, probably still in conjunction with internal asymmetry as well. This would yield an efficient mechanism to transport large amounts of fine sediment far up-river, as currently observed.     

Flocculation: a key process in the sediment flux of a large,glacier‐fed lake

This study investigates the consequences of flocculation for sediment flux in glacier‐fed Lillooet Lake, British Columbia based on density, fractal dimension, in situ profiles of sediment concentration and size distribution, and settling velocity equations presented in the literature. Sediment flux attributed to macroflocs during the late spring and summer accounts for a significant portion of sediment flux in the lake, equivalent to at least one‐quarter of the average annual sediment flux. Fine sediment is reaching the lake floor faster in flocs than occurs if settling as individual grains. This flux varies both spatially and temporally over the observation period, suggesting a link between deposition via flocculation and the properties of bottom sediments. Macrofloc flux increased through June, reached a peak during July, and then declined into August. Macrofloc flux was greatest in the distal end of the first basin, approximately 10 km from the point of inflow. Relatively high excess densities (～0·1 g cm^–3 at 500 µm) for flocs in situ are consistent with a composition dominated by inorganic primary particles. Microlaminations within Lillooet Lake varves have been linked by earlier workers to discharge events, and the action of turbidity currents, emanating from the Lillooet River. While turbidity currents undoubtedly occur in Lillooet Lake, these results demonstrate flocculation as an adjunct process linking discharge, lake level, macrofloc flux, bulk density and microlaminations. In situ measurements of sediment settling velocity in glacier‐fed lakes are required to better constrain flux rates, and permit comparison between flocculation in lacustrine environments with existing studies of estuarine, marine and fluvial flocculation. Copyright © 2009 John Wiley & Sons, Ltd.     

Modeling flocculation in a hypertidal estuary

When fine particles are involved, cohesive properties of sediment can result in flocculation and significantly complicate sediment process studies. We combine data from field observations and state-of-the-art modeling to investigate and predict flocculation processes within a hypertidal estuary. The study site is the Welsh Channel located at the entrance of the Dee Estuary in Liverpool Bay. Field data consist of measurements from a fixed site deployment during 12–22 February 2008. Grain size, suspended sediment volume concentration, and current velocity were obtained hourly from moored instruments at 1.5 m above bed. Near-bottom water samples taken every hour from a research vessel are used to convert volume concentrations to mass concentrations for the moored measurements. We use the hydrodynamic model Proudman Oceanographic Laboratory Coastal Ocean Modelling System (POLCOMS) coupled with the turbulence model General Ocean Turbulence Model (GOTM) and a sediment module to obtain three-dimensional distributions of suspended particulate matter (SPM). Flocculation is identified by changes in grain size. Small flocs were found during flood and ebb periods—and correlate with strong currents—due to breakup, while coarse flocs were present during slack waters because of aggregation. A fractal number of 2.4 is found for the study site. Turbulent stresses and particle settling velocities are estimated and are found to be related via an exponential function. The result is a simple semiempirical formulation for the fall velocity of the particles solely depending on turbulent stresses. The formula is implemented in the full three-dimensional model to represent changes in particle size due to flocculation processes. Predictions from the model are in agreement with observations for both settling velocity and SPM. The SPM fortnight variability was reproduced by the model and the concentration peaks are almost in phase with those from field data.     

Characterization of particle size and settling velocity of cohesive sediments affected by a neutral exopolymer

In natural waters,exopolymers or extracellular polymeric substances(EPS) exuded by microorganisms interact with clay particles,resulting in the flocculation of clays and hence alteration to the properties of suspended cohesive sediments.To investigate and further understand how neutral EPS affect cohesive sediment transport and the final sediment yield,an experimental study was conducted on laboratory-prepared clay and guar gum(used as an analog for neutral EPS) suspensions to characterize EPS-induced flocculation and the settling velocity of resultant floes.Four different clays consisting of kaolinite,illite,Ca-montmorillonite,and Na-montmorillonite were studied to examine the influence of different layer charges on clay flocculation induced by neutral EPS.Floc size was determined by a laser particle size analyzer,and settling velocity estimated by analyzing the time-series floc settling images captured by an optical microscope.Results indicate that neutral EPS promote clay-EPS flocculation for all four clays with the particle/floc size significantly increased from～0.1-60μm to as large as～600μm.Clays’ layer charge has a profound influence on the clay-EPS flocculation.With the same floc size,the settling velocity of clay-EPS flocs is typically smaller than that of pure clay flocs,which is attributed to the reduced density of flocs caused by the EPS. However,for flocs of the same composition(e.g.pure clay or hybrid clay-EPS mixture),the settling velocity increases with size.The fractal dimension of these clay-EPS flocs estimated from settling velocity ranges from 1.39 to 1.47,which are smaller than that of pure clay flocs,indicating that these flocs are less compacted than the pure clay flocs.     

A three-dimensional numerical model of particulate transport for coastal waters

A Lagrangian marker particle in Eulerian finite difference cell solution to the three-dimensional incompressible mass transport equation was developed for predicting particulate transport in coastal and estuarine waters. Special features of the solution procedure include a finite difference grid network which translates horizontally and vertically with the mean particle motion and expands with the dispersive growth of the marker particle cloud. The cartesian vertical coordinate of the three-dimensional mass transport equation has been transformed, using instantaneous water column depth to allow adaptation to flow situations with a temporally and spatially varying bottom topography and free surface. Results from this model for turbulent diffusion and advection of a uniform plug flow of sediment in an unbounded uniform flow field with various sediment settling velocities were in excellent agreement with the corresponding analytic solutions. Using current information from a two-dimensional vertically averaged hydrodynamic's model, the model was utilized to predict the long term diffusion and advection of dilute neutrally and negatively buoyant suspended sediment clouds resulting from a hypothetical instantaneous release of dredge spoil waste at Brown's Ledge in Rhode Island Sound.     

Observations of sediment resuspension and settling off the mouth of Jiaozhou Bay,Yellow Sea

We deployed bottom-mounted quadrapod equipped with acoustic Doppler current profiler (ADCP), acoustic Doppler velocimeter (ADV), and optical backscatter sensor (OBS) over two semidiurnal tidal cycles along the western coast of the Yellow Sea, China. In combination with shipboard profiling of CTD and LISST-100, we resolved the temporal and spatial distributions of tidal currents, turbulent kinetic energy (TKE), suspended sediment concentration (SSC) and particle size distributions. During the observations, tidal-induced bottom shear stress was the main stirring factor. However, weak tidal flow during the ebb phase was accompanied by two large SSC and median size events. The interactions of seiche-induced oscillations with weak ebb flow induced multiple flow reversals and provided a source of turbulence production, which stripped up the benthic fluff layers (only several millimeters) around the Jiaozhou Bay mouth. Several different methods for inferring mean suspended sediment settling velocity agreed well under peak currents, including estimates using LISST-based Stokes’ settling law, and ADCP-based Rouse profiles, ADV-based inertial-dissipation balance and Reynolds flux. Suspended particles in the study site can be roughly classified into two types according to settling behavior: a smaller, denser class consistent with silt and clay and a larger, less dense class consistent with loosely aggregated flocs. In the present work, we prove that acoustic approaches are robust in simultaneously and non-intrusively estimating hydrodynamics, SSC and settling velocities, which is especially applicable for studying sediment dynamics in tidal environments with moderate concentration levels.     

Response of tide-dominated mangrove shorelines in Northern Australia to anticipated sea-level rise

The exact response of mangrove shorelines to anticipated sea-level rise will depend upon the balance between sedimentation and sea-level change. Within the Top End of the Northern Territory of Australia there are extensive, relatively unmodified, tide-dominated mangrove forests, where tidal processes redistribute sediment. Harbours, such as Darwin Harbour, and tidal rivers, such as the South Alligator River with its associated coastal and estuarine plains, represent opposite extremes in terms of Holocene sedimentary infill, and will respond differently to sea-level rise. In Darwin Harbour, mangrove assemblages can be recognized in geomorphologically defined habitats. Similar topography within and between creeks implies morphodynamic equilibrium with tidal processes. Tidal reworking of sediment may maintain an equilibrial profile under gradually rising sea level, with resuspension of lower intertidal and subtidal muds and their redeposition within upper intertidal mangrove habitats. In contrast, the plains along the coast and tidal rivers draining into van Diemen Gulf developed during the post-glacial marine transgression, and since sea level stabilized, around 6000 years ago, coastal plains have prograded. These broad plains are presently not extensively influenced by salt water, but are often at elevations close to, or even below, modern high-tide levels. They may, therefore, revert to saline conditions particularly rapidly if the sea rises. The pattern of change may not be directly analogous to marine incursion experienced in the early Holocene, because broad plains have been able to prograde during the last 6000 years of relatively stable sea level.     

Flocculation and particle size analysis of expansive clay sediments affected by biological, chemical, and hydrodynamic factors

Expansive clay sediments are abundant in riverine and estuarine waters and bottom beds, and their particle size distributions (PSD) are important for the analysis of sediment transport. This paper presents an experimental study to evaluate, using a laser particle size analyzer under varying flow conditions, the intrinsic PSD of two expansive clays, a Ca- and a Na-montmorillonite and the influence of biological, chemical, and hydrodynamic factors on their flocculation and PSD. The considered biological factor consists of three extracellular polymeric substances of varying polarity, including xanthan gum, guar gum, and chitosan; the chemical factor is the salinity; and the hydrodynamic factor is the types of flow indicated by the Reynolds number and shear rate. The intrinsic PSD of both clays show a multimodal lognormal distribution with sizes ranging from 0.2 to 50 μm. All three biopolymers, xanthan gum, guar gum, and chitosan, can facilitate flocculation through long-range polymer bridging and short-range ion-dipole interaction, hydrogen bonding, and Coulomb force. The influence of salinity is different for the two clays: the particle size of the Na-montmorillonite increases with salinity, which is caused by flocculation resulting from the suppressed electrical double layer, while that of the Ca-montmorillonite is slightly reduced owing to the decreased basal spacing and cation exchange. For different hydrodynamic conditions, higher shear rate promotes the flocculation of Ca-montmorillonite, but breaks the Na-montmorillonite flocs. The significance of understanding the flocculation and PSD of expansive clays is also discussed in terms of sediment transport under different aquatic environments.     

Turbulence and suspended sediment processes in the Garonne River tidal bore in November 2016

A tidal bore is a water discontinuity at the leading edge of a ood tide wave in estuaries with a large tidal range and funneling topography. New measurements were done in the Garonne River tidal bore on 14 15 November 2016, at a site previously investigated between 2010 and 2015. The data focused on long, continuous, high-frequency records of instantaneous velocity and suspended sediment con- centration (SSC) estimate for several hours during the late ebb, tidal bore passage and ood tide. The bore passage drastically modi ed the ow eld, with very intense turbulent and sediment mixing. This was evidenced with large and rapid uctuations of both velocity and Reynolds stress, as well as large SSCs during the ood tide. Granulometry data indicated larger grain sizes of suspended sediment in water samples compared to sediment bed material, with a broader distribution, shortly after the tidal bore. The tidal bore induced a sudden suspended sediment ux reversal and a large increase in suspended sedi- ment ux magnitude. The time-variations of turbulent velocity and suspended sediment properties indicated large uctuations throughout the entire data set. The ratio of integral time scales of SSC to velocity in the x-direction was on average TE,SSC/TE,x 0.16 during the late ebb tide, compared to TE,SSC/ TE,x 0.09 during the late ood tide. The results imply different time scales between turbulent velocities and suspended sediment concentrations.     

Development of a Lagrangian sediment model to reproduce the bathymetric evolution of the Mersey Estuary

The development and application of a Lagrangian particle-tracking model to simulate sediment transport in the Mersey Estuary (UK) is described. Each of the particles (up to a million in total) is advected horizontally by the 3-D tidal currents. Related vertical movements are: (1) downwards by settlement at a prescribed velocity w _s and (2) both upwards and downwards by a distance related to the vertical eddy diffusivity. In a novel departure from traditional practice, where initial distributions of surficial sediments are specified, all particles are introduced at the seaward boundary of the model. Provenance studies indicate surficial sediments are overwhelmingly of marine origin. For the predominant fine sediments, ‘cyclical convergence’ in suspended sediment concentrations is approximated after about two spring-neap tidal cycles. Comparisons are shown between the suspended sediment concentrations and net deposition rates computed by this model against observed values and earlier computations utilising both 1-D and 3-D Eulerian models. While all of these results are in broad agreement, the flexibility of the Lagrangian approach for simulating flocculation, consolidation and mixed sediments illustrates its future potential.