JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 117, C03014, doi:10.1029/2011JC007552, 2012
Multimodality of a particle size distribution of cohesive suspended particulate matters in a coastal zone Byung Joon Lee,1 Michael Fettweis,2 Erik Toorman,1 and Fred J. Molz3 Received 27 August 2011; revised 20 December 2011; accepted 10 January 2012; published 9 March 2012.
[1] Particle size distributions (PSDs) of suspended particulate matters in a coastal zone are lognormal and multimodal in general. The multimodal PSD, which is caused by the mixing of multiple particle and aggregate size groups under flocculation and erosion/resuspension, is a record of the particle and aggregate dynamics in a coastal zone. Curve-fitting software was used to decompose the multimodal PSD into subordinate lognormal PSDs of primary particles, flocculi, microflocs, and macroflocs. The curve-fitting analysis for a time series of multimodal PSDs in the Belgian coastal zone showed the dependency of the multimodality on (1) shear-dependent flocculation in a flood and ebb tide, (2) breakage-resistant flocculation in the spring season, and (3) silt-sized particle erosion and advection in a storm surge. Also, for modeling and simulation purposes, the curve-fitting analysis and the settling flux estimation for the multimodal PSDs showed the possibility of using discrete groups of primary particles, flocculi, microflocs, and macroflocs as an approximation of a continuous multimodal PSD. Citation: Lee, B. J., M. Fettweis, E. Toorman, and F. J. Molz (2012), Multimodality of a particle size distribution of cohesive suspended particulate matters in a coastal zone, J. Geophys. Res., 117, C03014, doi:10.1029/2011JC007552.
1. Introduction [2] Particle size distributions (PSDs) of suspended particulate matters in a coastal zone are often multimodal with superposition of the subordinate lognormal PSDs [Christiansen et al., 2006; Mikkelsen et al., 2006]. The log normality describes a more or less skewed distribution toward small size and the multimodality describes a distribution consisting of multiple modal peaks. The multimodal PSD is common in coastal water because of mixing of multiple particle and aggregate size groups. Suspended particulate matter in coastal water is composed of various constituents including clay minerals, organic matter, salts, trace metals and so on [Berlamont et al., 1993; Maggi, 2009]. When such various constituents are subject to flocculation (i.e., aggregation and breakage) and transport (i.e., in and out flux) in highly dynamic tidal currents of a coastal zone, they will never reach an equilibrium state but rather undergo a transient state while alternating the rise and sink of the modal peaks [Gibbs et al., 1989; Li et al., 1993; Yuan et al., 2009]. [3] Erosion of large sand particles and/or sand-sized clay aggregates from the seabed often causes the mixing of different particle and aggregate size groups and develops the 1 Hydraulics Laboratory, Department of Civil Engineering, Katholieke Universiteit Leuven, Heverlee, Belgium. 2 Management Unit of the North Sea Mathematical Models, Royal Belgian Institute of Natural Sciences, Brussels, Belgium. 3 Department of Environmental Engineering and Earth Sciences, Clemson University, Anderson, South Carolina, USA.
Copyright 2012 by the American Geophysical Union. 0148-0227/12/2011JC007552
multimodal PSDs in a coastal zone [Christiansen et al., 2006; Yuan et al., 2009]. For example, Yuan et al. [2009] reported that the mixing of advected clays/silts from the upstream and eroded sands from the local seabed causes a bimodal PSD of suspended sediments. Individual clay/silt particles, which are presumably noncohesive considering no reported flocculation at the measuring site, become predominant in a mild tide without erosion and resuspension of local sands and compose a unimodal PSD. However, sands start eroding and resuspending from the seabed at the critical shear stress for erosion (t c), and the mixture of “background” clays/silts and “eroded/resuspended” sands develop a bimodal PSD. When the tidal current becomes mild again, slow-settling clays/silts remain in suspension but fast-settling sands deposit back to the seabed. The mixture of clays/silts and sands will continuously change their relative mass fractions and alternately develop a unimodal, bimodal or multimodal PSD in the flowvarying tidal cycle of a coastal zone. Often sediments in coastal areas, such as the Belgian nearshore zone, consist of fine mud and coarse sand. The mud–sand ratio influences the transition between cohesive and noncohesive behavior, and therefore affects erosion, SPM concentration and composition as well as on the PSD of the suspended material [Fettweis et al., 2012; van Ledden et al., 2004; Manning et al., 2010]. If the ratio of cohesive mud to noncohesive sand is high enough to prevent the seabed from erosion and resuspension in a flood and ebb tide, a multimodal PSD could be caused by processes other than erosion and resuspension. [4] Flocculation of fine-grained cohesive sediments may also develop a multimodal PSD in a coastal zone. Flocculation is a dynamic process combining the counteracting processes of aggregation and breakage in a fluid shear field [Eisma, 1986; van Leussen, 1994; Maggi, 2005; Mietta,
C03014
1 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 1. A multimodal PSD and schematic diagrams of the discrete aggregate groups of primary particles, flocculi, microflocs, and macroflocs. 2010; Winterwerp, 1998]. Whereas aggregation assembles aggregates with building blocks (i.e., fine-grained cohesive sediments), breakage disassembles aggregates down to building blocks. Flocculation continuously changes the relative volume fractions of building blocks and aggregates and develops distinct peaks in a PSD, while balancing aggregation and breakage in a flow-varying tidal cycle [Verney et al., 2011]. Furthermore, flocculation of marine or coastal clay minerals has been reported to develop a four-level ordered conceptual structure, namely, primary particles, flocculi, microflocs and macroflocs (Figure 1) [van Leussen, 1994]. Primary particles consist mostly of fine particles with a wide size range of 0.25 2.5 mm, including fragmented clay minerals, organic and calcareous particles, picophytoplankton and heterotrophic bacteria [Andrews et al., 2010]. Flocculi, a compound word of floc and nuclei, consist of strongly bound clay minerals with a size range of 10 20 mm. Flocculi are the smallest clay-based aggregates which are seldom broken down to the lowest-level primary particles (e.g., clay minerals) even at the highest turbulent shear of a coastal zone, and therefore they may be regarded as the basic and major building blocks of aggregates besides primary particles [van Leussen, 1994]. Microflocs consist of flocculi and primary particles with a size range of 50 200 mm. Macroflocs are built from microflocs, primary particles and flocculi, and have a size range of hundreds to thousands of micrometers. Macroflocs are sometimes referred to as marine snow because of their large and fluffy structure [Alldredge and Silver, 1988]. In the flow-varying tidal cycle of a coastal zone, primary particles, flocculi, microflocs, and macroflocs are highly interactive, resulting in a transient, multimodal PSD [Manning and Bass, 2006; Mikkelsen et al., 2006]. [5] A time-varying multimodal PSD constitutes a scientific record of flocculation and transport of the constituent particles and aggregates. One can use a time series of such PSDs to understand the possible causes of a multimodal PSD and investigate the particle and aggregate dynamics. For example,
C03014
the occurrence of an additional peak at sand size besides an existing peak at clay/silt size can be evidence of erosion of sands or aggregation of fine-grained cohesive sediments. However, most current researchers have used observed multimodal PSDs in a qualitative and descriptive manner for investigating the particle and aggregate dynamics in a coastal zones or estuaries [Chen et al., 2005; Manning and Bass, 2006; Mikkelsen et al., 2006]. Thus, a quantitative and systematic statistical approach may provide new information for interpreting the particle and aggregate dynamics underlying multimodal PSDs. [6] A quantitative statistical method is also required for practical purposes in modeling and simulation. The recent mathematical models with multiple and discrete aggregate size groups are more accurate for simulating flocculation and transport of multimodal PSDs than models with a single size group under unimodal approximation [Lee et al., 2011; Maerz et al., 2011]. Therefore, characterization of each modal peak of a multimodal PSD and logical simplification of such a modal peak as a discrete aggregate group are required for better modeling and simulation. A statistical method, which had been developed and applied to investigate aerosol dynamics [Banerjee et al., 2009; Hussein et al., 2005; Whitey, 2007], is adopted herein to analyze and interpret a coastal zone multimodal PSD in a quantitative and systematic way. The statistical method with automated curve-fitting software (DistFit™, Chimera Technologies Inc., USA) helped us to decompose a multimodal PSD to the subordinate lognormal PSDs of primary particles, flocculi, microflocs and macroflocs, to quantify the mean diameter, standard deviation and volume fraction of the subordinate PSDs, and to investigate the statistics of a large and complicated PSD time series in a quantitative and systematic way. For modeling and simulation purposes, this paper will apply the selected statistical method to analyze the time series of the multimodal PSDs in terms of the characteristics of the subordinate lognormal PSDs.
2. Materials and Methods 2.1. Regional Settings [7] The study area is situated in the Belgian coastal zone facing the southern North Sea. The coastal zone is characterized by a mean tidal range between 4.3 and 2.8 m at spring and neap tide, respectively, and by maximum current velocities of more than 1 m/s. The wind from the southwest dominates for 33% of the year, followed by the wind from the northeast. The maximum wind speed coincides with the wind from the southwest, but the highest wave action occurs under the wind from the northeast [Fettweis et al., 2012]. The area is of interest due to the occurrence of a highly turbid coastal zone. The SPM concentration ranges between a minimum of 20–70 mg/L and a maximum of 100–1000 mg/L. SPM concentrations decrease during periods of southwest wind due to the advection of less turbid English Channel water, but increase during northeast winds due to an increased SPM outflow from the Westerschelde estuary [Baeye et al., 2011]. The nearshore area is characterized by bottom sediments varying from pure sand to pure mud [Verfaillie et al., 2006]. The sediment underneath fine-grained cohesive sediment suspension at the measuring site consists of silts and fine
2 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
sands with a median grain size of about 150 mm. Ephemeral fluffy layers of fine-grained cohesive sediment have been frequently observed on the top of the silty and sandy bed [Baeye et al., 2011]. The suspended matters in the same area consist mostly of clay (41%) and silt (57%) sized particles with a small amount of sand (2%); the mean primary particle size is 1.1 3.7 mm. 2.2. Instrumentation [8] A tripod was deployed at about 1 km from the shore line in a water depth of 5.5 m below MLLWS (mean lowest low water at spring tide), and current velocity, turbulent shear, salinity, and suspended sediment data were collected for the two field campaigns from 28 January 2008 to 11 February 2008 and from 15 April 2008 to 30 April 2008. The measuring instruments on the tripod were a 5 MHz SonTekR Acoustic Doppler Velocimeter (ADV Ocean), a 3 MHz SonTekR Acoustic Doppler Profiler (ADP), two D&AR optical backscatter sensors (OBS), a Sea-BirdR SBE37 CT system, a Sequoia ScientificR Laser In Situ Scattering and Transmissometer 100C (LISST-100C), and two SonTekR Hydra systems for data storage and batteries. The measured values from OBS, LISST, and ADV were averaged over each 10 min period and recorded. Two OBSs measured SPM concentrations at 0.2 and 2 m above the bed (mab) and the LISST was attached at 2 mab. The ADV measured a time series of flow velocity and turbulent intensity at 0.2 mab. The high-frequency ADV (25 Hz) measurements permitted to decompose the velocity in terms of a mean and a fluctuating part and to calculate the turbulent kinetic energy (TKE). TKE is proportional to bottom shear stress [Pope et al., 2006; Thompson et al., 2003; Verney et al., 2007]. However, in the presence of waves the bottom shear stress was estimated with the inertial dissipation method adopting the spectrum of the vertical velocity component and a correction for the advection by waves [Fettweis et al., 2010; Sherwood et al., 2006; Trowbridge and Elgar, 2001]. The LISST-100C measured particle size distributions (PSDs) in 32 logarithmically spaced size groups over the range of 2.5–500 mm [Agrawal and Pottsmith, 2000]. The volume concentration of each size group was estimated with an empirical volume calibration constant, which was obtained under a presumed sphericity of particles. It is important to note that uncertainties using LISST-100C detectors may arise due to nonspherical flocs, floc sizes exceeding the instrument range, a too high SPM concentration, or stratification of the water column [Agrawal and Pottsmith, 2000; Andrews et al., 2010; Fettweis, 2008; Mikkelsen et al., 2007]. Air bubbles can make errors of the LISST measurements, especially due to breaking waves during the storm period. However, the water depth was high enough to prevent breaking waves and trapping air bubbles during the measuring campaign [U.S. Army Corps of Engineers, 2003; Fettweis et al., 2012]. 2.3. Curve Fitting and Decomposition of a Multimodal PSD [9] The observed multimodal (4-peaked) PSD was assumed to be formed by overlapping four unimodal lognormal PSDs of primary particles, flocculi, microflocs, and macroflocs (Figure 1), and these data served as input to the selected software which performed the curve fitting analysis and the decomposition of the observed PSD into the subordinate
C03014
lognormal PSDs (equation (1)) [Hussein et al., 2005; Mäkelä et al., 2000; Whitey, 2007]. " # 4 iÞ 2 Vi dV X 1 lnðD=D pffiffiffiffiffiffi ¼ exp dD 2 lnðsi Þ 2p lnðsi Þ i¼1
ð1Þ
[10] In equation (1), V and D are the volumetric concentration and diameter of each size interval of the LISST-100C measured PSD. dV/dD is the normalized volumetric fraction by the width of the size interval and used for curve fitting to i is the geometric a lognormal distribution [Hinds, 1999]. D mean diameter (Dg;i ), which is equivalent to the median diameter (D50,i) for a lognormal distribution (E. W. Weisstein, Log normal distribution, from MathWorld— A Wolfram Web Resources, 2006, available at http://mathworld.wolfram.com/LogNormalDistribution.html), si is the geometrical standard deviation, and V i is the volumetric concentration of an ith unimodal PSD. Equation (1) has i, si and V i of primary partitwelve fitting parameters for D cles, flocculi, microflocs and macroflocs. Based on observation, the mean diameters of microflocs and macroflocs were set to change as curve-fitting parameters in the ranges of 20 200 and 200 500 mm, respectively. However, the mean diameters of primary particles and flocculi were fixed at 3 and 15 mm because the bands of those mean diameters were found to be narrow within about 1 and 5 mm, respectively. These approximations reduce the number of curve-fitting parameters and the complexity of the curve-fitting procedure. The standard deviation of a subordinate PSD (si) was limited under 3 to prevent unrealistically wide PSDs comprising the best quality PSD. [11] DistFit™ (Chimera Technologies, USA), which has been widely used in analyzing aerosol particles, was used to find the best quality fit to a measured PSD. The best quality PSD was defined by the minimum error between the simulated and measured PSDs [Hussein et al., 2005; Whitey, 2007]. Noteworthy is that the recent version of DistFit™ has an automated fitting algorithm to handle a large number of PSDs. Therefore, DistFit™ can free users from the laborious task for finding the best quality PSD curves for a large number of measured PSDs. (In the present application, data were collected at 10 min intervals for two to three weeks.) i, si and The best quality PSD provided parameter values (D V i) for the subordinate lognormal PSDs of primary particles, flocculi, microflocs and macroflocs. The quality of the curve fitting analysis was monitored with the sum of errors between simulations and experimental data (c) (e.g., Table 1). The resulting parameter values were then used as indicators of the particle and aggregate dynamics. 2.4. Estimation of the Settling Flux of a Multimodal PSD [12] The reference sediment settling flux of a PSD (SFref) was calculated by summing the settling fluxes of all the 32 aggregate size intervals of the LISST instrument. Each settling flux was calculated by multiplying the settling velocity and volume concentration of an aggregate size interval of the LISST instrument. The settling velocity was calculated by the modified Stokes’ equation with floc fractal dimensions (equation (2)) [Winterwerp, 1998]. The reference settling flux was then compared to the estimated settling fluxes,
3 of 17
4 of 17
1546 427 2530 532 1392 104 278 1224 391
143 512 188 96 175 260 194
119.80 119.93 120.08 31.21 31.44 31.59 31.69
VT
111.96 112.13 112.22 112.41 39.06 39.19 39.30 39.47 119.65
Day
0.130 0.375 0.160 0.031 0.046 0.080 0.052
0.399 0.358 0.753 0.313 0.288 0.090 0.158 1.086 0.256
SFref
119 98 109 43 37 44 36
26 104 31 80 26 126 75 126 79
D50
0.087 0.258 0.106 0.022 0.034 0.060 0.036
0.211 0.230 0.410 0.220 0.191 0.067 0.109 0.793 0.161
SFD50
33.1 31.2 33.8 29.0 26.1 25.0 30.8
47.1 35.8 45.6 29.7 33.7 25.6 31.0 27.0 37.1
ErD50
187 149 173 63 51 60 53
52 172 59 118 41 175 114 181 133
Dwt
0.135 0.388 0.165 0.032 0.046 0.081 0.053
0.417 0.371 0.783 0.322 0.296 0.092 0.163 1.119 0.266
SFDwt
3.8 3.5 3.1 3.2 0.0 1.3 1.9
4.5 3.6 4.0 2.9 2.8 2.2 3.2 3.0 3.9
ErDwt
142 501 175 11195 173 254 188
1496 420 2414 504 1347 100 269 1192 388
V1
150 110 154 46 40 48 40
27 136 35 95 29 145 86 150 86
Dg,1
3.0 2.5 2.4 2.4 2.1 2.1 2.4
3.0 3.0 3.0 2.2 2.5 2.3 2.4 2.4 3.0
sg,1
0.201 0.088 0.166 0.036 0.040 0.042 0.056
0.070 0.158 0.085 0.071 0.063 0.148 0.057 0.119 0.089
c1M
0.191 0.434 0.199 0.034 0.048 0.083 0.057
0.379 0.512 0.796 0.335 0.310 0.102 0.174 1.286 0.310
SF1M
Curve Fitting: One Mode
46.9 15.7 24.4 9.7 4.3 3.8 9.6
5.0 43.0 5.7 7.0 7.6 13.3 10.1 18.4 21.1
Er1M
29 55 39 21 40 45 59
408 85 1053 80 596 10 39 124 55
V1
15 15 15 15 15 15 15
15 15 15 15 15 15 15 15 15
Dg,1
2.4 2.4 2.3 2.3 2.2 2.1 2.4
2.6 2.5 2.8 2.3 2.6 2.1 2.4 2.3 2.5
sg,1
113 453 149 75 135 214 134
1095 340 1391 451 772 92 238 1099 332
V2
195 120 180 57 48 55 53
33 170 56 104 42 157 96 160 102
Dg,2
2.0 2.2 1.9 2.0 1.8 1.8 1.9
3.0 2.1 2.3 2.0 1.9 2.0 2.1 2.1 2.5
sg,2
0.139 0.067 0.083 0.009 0.009 0.006 0.019
0.069 0.106 0.077 0.019 0.046 0.112 0.022 0.087 0.078
c2M
Curve Fitting: Two Modes
0.141 0.391 0.170 0.031 0.045 0.078 0.053
0.403 0.392 0.737 0.311 0.294 0.094 0.160 1.174 0.270
SF2M
8.5 4.3 6.3 0.0 2.2 2.5 1.9
1.0 9.5 2.1 0.6 2.1 4.4 1.3 8.1 5.5
Er2M
a VT and Vi represent the volumetric concentration of a multimodal PSD and a decomposed PSD, respectively. SF represents the settling flux of a PSD. Er represents the error between the reference and estimated fluxes for a PSD. Here c is the sum of errors between measurement and simulation for a PSD, representing the quality of simulation.
Storm surge
Shear-resistant flocculation
Neap tide
Spring tide
Group
Measurement
Table 1. A Summary of Median and Mean Diameters, Best Fit Parameters From Curve Fitting Analysis With One to Four Particle and Aggregate Groups, and Estimated Settling Fluxes for the PSDs in Figures 4, 6, 8, and 10a
C03014 LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS C03014
5 of 17
Storm surge
Shear-resistant flocculation
Neap tide
Spring tide
Group
15 15 15 15 15 15 15 15 15
15 15 15 15 15 15 15
24 61 38 21 41 46 59
Dg,1
803 82 1135 82 615 9 41 127 64
V1
Table 1. (continued)
2.2 2.3 2.3 2.3 2.1 2.1 2.4
2.6 2.5 2.7 2.3 2.5 2.0 2.3 2.2 2.4
sg,1
71 393 114 75 132 212 134
629 252 1251 444 741 65 230 899 271
V2
118 106 142 57 48 55 53
43 126 57 104 42 112 94 132 85
Dg,2
1.9 2.0 1.8 2.0 1.8 1.8 1.9
2.2 1.9 2.0 1.9 1.9 1.9 2.0 1.9 2.0
sg,2
48 58 36 0 1 1 0
75 92 73 7 11 29 7 197 55
V3
315 368 325 213 240 277 224
262 337 347 423 225 288 370 354 352
Dg,3
Curve Fitting: Three Modes
1.3 1.2 1.2 2.9 1.4 1.4 2.9
1.5 1.3 1.3 1.1 1.5 1.3 1.2 1.2 1.2
sg,3
0.039 0.030 0.037 0.009 0.006 0.005 0.019
0.065 0.043 0.070 0.014 0.045 0.029 0.015 0.045 0.033
c3M
0.131 0.379 0.160 0.031 0.046 0.079 0.053
0.408 0.363 0.769 0.315 0.298 0.089 0.161 1.101 0.255
SF3M
0.8 1.1 0.0 0.0 0.0 1.3 1.9
2.3 1.4 2.1 0.6 3.5 1.1 1.9 1.4 0.4
Er3M
0 0 0 1 0 0 4
104 6 180 1 70 0 1 2 0
V1
3 3 3 3 3 3 3
3 3 3 3 3 3 3 3 3
Dg,1
1.1 1.6 1.6 1.1 2.4 1.1 1.1
1.1 1.1 1.1 1.1 1.1 1.5 1.1 1.1 2.2
sg,1
23 61 38 21 41 45 53
682 72 901 79 541 9 39 122 64
V2
15 15 15 15 15 15 15
15 15 15 15 15 15 15 15 15
Dg,2
2.2 2.3 2.3 2.2 2.1 2.1 2.1
2.1 2.2 2.1 2.2 2.0 2.0 2.2 2.2 2.4
sg,2
71 393 114 75 132 212 136
676 256 1352 442 758 65 231 901 271
V3
117 106 142 57 48 55 53
47 123 57 103 45 112 94 132 85
Dg,3
1.9 2.0 1.8 2.0 1.8 1.8 1.9
2.1 1.9 1.9 1.9 1.8 1.9 2.0 1.9 2.0
sg,3
48 58 36 0 1 1 0
77 93 86 9 19 29 8 198 55
c3M
Curve Fitting: Four Modes
314 368 325 200 240 266 429
262 336 329 388 200 288 370 354 352
SF3M
1.3 1.2 1.2 2.0 1.4 1.4 1.1
1.5 1.3 1.3 1.1 1.5 1.3 1.2 1.2 1.2
Er3M
0.039 0.030 0.037 0.007 0.008 0.005 0.009
0.024 0.038 0.025 0.014 0.015 0.029 0.014 0.044 0.033
V1
0.8 1.1 0.0 0.0 0.0 1.3 0.0
0.8 1.1 1.1 0.6 0.3 1.1 1.9 1.4 0.4
Dg,1
0.3 0.2 0.6 1.8 1.3 1.9 1.7
3.2 0.2 3.0 0.7 1.9 1.8 0.1 0.5 1.4
sg,1
C03014 LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS C03014
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
which were calculated with the median diameter (D50), the volume-weighted mean diameter (Dwt) and the mean dia i) obtained from meters of the discrete aggregate groups (D the curve-fitting analysis. The purpose of the comparison between the reference and estimated settling fluxes was to check if the discrete aggregate groups of primary particles, flocculi, microflocs and macroflocs could approximate a continuous multimodal PSD in the settling flux calculation. SF ¼
X X 1 ðrs rw Þ g 3ni Dni i 1 DP Vi ws;i ¼ Vi 18 m 1 þ 0:15Rei0:687 i i
ð2Þ
[13] In equation (2), the volumetric concentration of an aggregate size interval (Vi) was obtained from the LISST instrument or the curve fitting analysis. The densities of primary particle and seawater (rs and rw) were fixed at 2600 and 1030 kg/m3. The diameter of primary particles (DP) was approximated at 2 mm from the previous research [Fettweis, 2008]. The gravitational acceleration and the fluid viscosity (g and m) were 9.81 m/s2 and 0.0001 kg/m/s. Rei represents the Reynolds number of an aggregate. The two stepwise fractal dimensions (ni) of 2.5 and 2.0 were used to represent the two distinct structures of less porous and hard building blocks (primary particles and flocculi; Di < 20 mm) and highly porous and fluffy flocs (microflocs and macroflocs; 20 mm < Di < 500 mm) [Winterwerp and van Kesteren, 2004].
3. Results and Discussion [14] Each particle size distribution (PSD) alternately developed one to four modal peaks in a flood/ebb tidal cycle. The first and second modal peaks of the PSD appeared consistently at about 3 mm and 15 mm for primary particles and flocculi, whereas the third and fourth modal peaks varied in the size ranges of 20 200 mm and > 200 mm for microflocs and macroflocs (Figure 1). In the entire measuring period, the median diameters (D50) of the PSDs seemed dependent on turbulent shear and temperature (Figure 2). High turbulent shear decreased the median diameter down to about 50 mm. Low turbulent shear and/or temperature rise concurred with increasing the median diameters up to about 100 mm. [15] The hydrodynamic condition (i.e., flow velocity and turbulent intensity) of the measuring site varied with storm, neap, and spring tidal periods (Figure 2). A storm period was characterized by more disordered flow vectors and five to ten times higher wave height and turbulent shear than the other regular neap and spring tidal periods (Figure 2a). Both the neap and spring tidal periods were found to have ordered flow vectors, which direct toward or away from the shore line for the flood and ebb currents, respectively (Figures 2a and 2b). However, the neap tide had two to three times lower flow velocity and turbulent shear than the spring tide, and it had a more asymmetric flood and ebb tidal current (i.e., high flow velocity for a flood current but low flow velocity for an ebb current) (Figure 2a). [16] The following paragraphs will explain in more detail the unique behaviors of the PSDs under flocculation and transport for the periods of (1) a spring tide with a symmetric
C03014
tidal current, (2) a neap tide with a flood-dominant asymmetric tidal current, (3) breakage-resistant flocculation in the spring season and (4) a storm surge. 3.1. Spring Tide With a Symmetric Tidal Current [17] The spring tide in Figure 2b had a symmetric pattern of flow velocity and turbulent shear for the flood and ebb currents. The highest turbulent shear coincided with the flood and ebb peak flow, and the lowest shear coincided with the flood-to-ebb and ebb-to-flood slack waters (i.e., zero flow velocity) (Figure 3a). The spring tide with a symmetric tidal current also caused a symmetric behavior of a multimodal PSD for the flood and ebb currents. The highest turbulent shear at both the flood and ebb peak flows increased the volume fractions of building blocks of primary particles and flocculi (i.e., the section between the bottom and the line of PP + Flocculi in Figure 3b) but decreased the volume fractions of aggregates of microflocs and macroflocs (i.e., the section between the line of PP + Flocculi and the top in Figure 3b) and skewed the PSD toward small size with a rising tail at the smallest size group (Figures 4a and 4c). In contrast, the lowest shear at both the slack waters increased the volumetric fractions of aggregates but decreased the volumetric fractions of building blocks (Figure 3b), and skewed the PSD toward large size with a distinct peak or hump of macroflocs (Figures 4b and 4d). [18] The multimodality of a PSD appears influenced by shear-dependent flocculation, in which aggregation and breakage counteract each other in a flow- and shear-varying tidal cycle [Fettweis et al., 2006; Winterwerp and van Kesteren, 2004; Winterwerp et al., 2006]. The highest shear of the peak flow would be expected to enhance breakage of aggregates to building blocks, while the lowest shear at slack water would enhance aggregation of building blocks to aggregates [Lee et al., 2011; Verney et al., 2009, 2011]. In addition to shear-dependent flocculation, one may suspect that turbulence-mediated erosion and resuspension would increase the influx of large particles and aggregates from the seabed into the water column [Orange et al., 2005; Yuan et al., 2009]. However, the peak flows, which are the best condition for erosion and resuspension of larger particles, instead decreased the volume fraction of large particles and aggregates, pointing thus to mainly floc breakup. Only a small amount of sands, less than 2% of the total suspended matters (see also section 2.1) was responsible for the tiny peaks of the macrofloc-sized class at the peak flows (Figures 4a and 4c). Therefore, erosion and resuspension might neither cause the influx of large particles and aggregates nor affect the multimodality in the nonstormy regular tidal cycle of the measuring site at 2 mab or higher in the water column. 3.2. Neap Tide With a Flood-Dominant Asymmetric Current [19] A neap tide at the measuring site had a flood-dominant asymmetric current, which consists of a shorter flood cycle with a higher flow velocity but a longer ebb cycle with a lower flow velocity. A sharp peak of flow velocity and turbulent shear was found in a flood cycle, whereas a stretched mound was found in an ebb cycle (Figure 5a). The geographical characteristics of a coastal zone and/or the periodic
6 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 2. A time series of the flow velocity vector with x and y direction components (U and V), median diameter (D50), turbulent shear (t), and temperature (°C) in the measuring periods of (a) 28 January 2008 to 11 February 2008 and (b) 15 April 2008 to 5 June 2008. The dotted lines show the smoothed curves with the moving averages of 70 data points. The shaded area represents a tidal cycle which has a specific flow condition and flocculation behavior.
7 of 17
C03014
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 3. (a) A time series of time-averaged flow velocities in x and y direction (U and V) and turbulent shear (t) and (b) a time series of the cumulative volumetric fractions of primary particles (PP), flocculi, microflocs (Micro), and macroflocs (Macro), for the spring tide with a symmetric tidal current. The symbols represent the measured data, and the lines show the smoothed curves with the moving averages of three data points. In Figure 3b, each section between two lines represents the individual volume fraction of primary particles, flocculi, microflocs, and macroflocs. The four vertical lines, from the left to right, indicate the flood peak flow, the flood-to-ebb slack water, the ebb peak flow, and the ebb-to-flood slack water, of which PSDs are plotted in Figure 4 for comparison. change of oceanic currents have been reported to cause such an asymmetric tidal current [Brown and Davies, 2010; Latteux, 1995; van der Molen et al., 2004]. The flooddominant asymmetric current caused asymmetric behavior of the resulting multimodal PSD for the flood and ebb currents. For example, macroflocs appeared in the entire ebb cycle even at the ebb peak flow but disappeared in the flood cycle. Building blocks (i.e., primary particles and flocculi) developed a sharp rise of the volumetric fraction at the flood peak flow but did not at the ebb peak flow (Figure 5b). The PSDs in the ebb cycle and the slack water skewed toward large size with a peak or hump of macroflocs (Figures 6b, 6c, and 6d), but the PSDs in the flood cycle skewed toward small size (Figure 6a).
C03014
[20] The neap tide with a flood-dominant asymmetric current generally enhanced the overall flocculation rate (i.e., increased aggregation and decreased breakage) more than the spring tide with a symmetric current, considering that the neap tide had a smaller volume of building blocks but a larger volume of aggregates especially during the ebb cycle. Why did the neap tide enhance the flocculation rate? First, the lower turbulent shear in the neap tide, which was about one third of the shear in the spring tide, necessarily enhanced the flocculation rate (Figures 3a and 5a). Flood/ebb flow imbalance further decreased turbulent shear and enhanced the flocculation rate during the ebb cycle by keeping the turbulent shear under a threshold of floc breakage. Second, stratification of sediment concentration in the water column might also enhance the flocculation rate, especially in the mild flow and shear condition of the ebb current. Winterwerp [2011] reported that the lower turbulent shear of the ebb current could reduce vertical dispersion and stratify sediment concentration (i.e., the deeper water depth, the higher sediment concentration) in the cohesive sediment suspension in the Ems-Dollard estuary, Germany/Netherlands. Similarly, the low-turbulence ebb current at the measuring site might enhance vertical stratification of the cohesive sediment suspension. The highly concentrated sediments near the bottom could enhance the flocculation rate by increasing the frequency of aggregate collision and attachment. In contrast, the higher turbulent shear of the flood current could increase vertical dispersion and homogenize sediment concentration in the water column in a breakage-dominant condition. In flocculation theory, turbulent shear (G) and sediment concentration (C) control the flocculation rate [van Leussen, 1994; Winterwerp and van Kesteren, 2004], we therefore conclude that substantial turbulent shear reduction and possible sediment stratification during the ebb cycle are the major causes of the enhanced flocculation rate in the neap tide, compared to the spring tide. 3.3. Breakage-Resistant Flocculation [21] The pattern of flow velocity during the observed breakage-resistant period was similar to that in the spring tide with a symmetric tidal current. Turbulent shear remained two to three times lower than the shear in the spring tide, but was still higher than the shear in the neap tide (Figures 3a and 7a). However, the flocculation rate in this period was higher than in any other periods, i.e., less building blocks of primary particles and flocculi but more aggregates of microfloc and macroflocs due to increasing aggregation but decreasing breakage. For example, the volume fraction of macroflocs (i.e., the section between the line of PP + Flocculi + Micro and the top in Figure 7b) reached up to 40% at slack water and still remained high at 10% even at the peak flow. It is more important to note that the volume fraction of building blocks (i.e., the section between the line of PP + Flocculi and the bottom in Figure 7b) remained low (under 20%) without substantial floc breakage in the entire flood and ebb cycle. Therefore, all the PSDs skewed toward large size with a large volume fraction of aggregates, even at the highest turbulent shear of the peak flow (Figure 8). [22] This period had calm weather and low tidal amplitude without storms, which could increase the flocculation rate due to an aggregation-dominant condition with low turbulent shear. However, flocculation in this period did not follow
8 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 4. (a) A PSD in the flood peak flow, (b) a PSD in the flood-to-ebb slack water, (c) a PSD in the ebb peak flow, and (d) a PSD in the ebb-to-flood slack water for the spring tide with a symmetric tidal current. Meas represents a PSD measured with the LISST instrument. PP, Flocculi, Micro, Macro, and Sum represent the decomposed PSDs of primary particles, flocculi, microflocs and macroflocs, and the superposition of the decomposed PSDs.
9 of 17
C03014
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
C03014
would be expected to stimulate growth of biological species (e.g., phytoplankton, bacterioplankton and benthic fauna) in the coastal zone [Droppo et al., 2005]. Simultaneously, the algae bloom and associated biological activity might enhance biologically mediated flocculation by increasing the production of particle-binding microbial exudates such as Extracellular Polymeric Substances (EPS) or Transparent Exopolymeric Particles (TEP) [Chen et al., 2005; Droppo, 2001; Droppo and Amos, 2001]. Such biologically mediated flocculation might increase not only floc size and but also breakageresistance at the measuring site, and therefore prevented floc breakage into building blocks. This may explain the fact that the volume fractions of building blocks remained low under 20% irrespective of turbulent shear (Figure 7b).
Figure 5. (a) A time series of time-averaged flow velocities in x and y direction (U and V) and turbulent shear (t) and (b) a time series of the cumulative volumetric fractions of primary particles (PP), flocculi, microflocs (Micro), and macroflocs (Macro), for the neap tide with a flood-dominant asymmetric tidal current. The symbols represent the measured data, and the lines show the smoothed curves with the moving averages of three data points. In Figure 5b, each section between two lines represents the individual volume fraction of primary particles, flocculi, microflocs, and macroflocs. The four vertical lines, from the left to right, indicate the flood peak flow, the flood-to-ebb slack water, the ebb peak flow, and the ebb-to-flood slack water, of which PSDs are plotted in Figure 6 for comparison. the typical pattern of shear-dependent flocculation (see sections 3.1 and 3.2). For example, the volume fractions of building blocks remained low without developing a significant peak even at the highest turbulent shear of the flood or ebb peak flow (Figure 7b). Aggregates seemed armored against breakage. This might indicate the involvement of other breakage-resistant flocculation mechanisms besides the shear-dependent mechanism. Biologically mediated flocculation could be one of the most probable causes of breakageresistant flocculation [Maggi, 2009]. [23] An algae bloom was observed in the vicinity of the measuring site, concurring with temperature rise and high nutrient influx in the manuring season (Figure 2b). This
3.4. Storm Period [24] Turbulent shear in the storm period was about ten times higher than in the other nonstorm periods, and so did not favor flocculation (Figure 9a). Under such strong turbulent shear, each volumetric fraction of primary particles, flocculi, microflocs and macroflocs was rather consistent. The volumetric fractions of primary particles and flocculi were consistent at about 0 and 20%, respectively. The volumetric fraction of macroflocs remained low, mostly under 10% (i.e., the section between the line of PP + Flocculi + Micro and the top in Figure 9b). Therefore, all the PSDs in the storm period developed a similar structure, having a dominant modal peak at about 50 mm with a small tail of building blocks (Figure 10). [25] The unique PSD during the storm period might be caused by influx of sand- and silt-sized particles but not by flocculation of suspended cohesive sediments. Advection of oceanic water masses by wind-induced currents with lower SPM concentration resulted in disappearance of high concentrated mud suspension near the bed and consequently in erosion of the previously covered local sand- and silt-sized particles at the measuring site. Changes in PSD distribution and median diameter during these events have been reported [Fettweis et al., 2012]. Results from the PSD analysis in this research also support this finding. A PSD is composed of breakage-resistant primary sediments, when it is subject to the highest turbulent shear in a storm surge or a peak flood/ ebb flow, and therefore the PSD gives information about the origin and characteristics of the lowest-level sediment composition. The PSD in the storm period consisted mostly of silt-sized particles (80% volume fraction; 50 60 mm diameter) (e.g., Figure 10a), whereas the PSD at the peak flow in the regular tidal period consisted of primary particles (5 10% volume fraction; 2.5 mm diameter), flocculi (40 50% volume fraction; 15 mm diameter) and small-size microflocs (30 40% volume fraction; 40 50 mm diameter) (e.g., Figure 4a). The observed difference of the sediment composition proves an influx of eroded and advected silt-sized particles replacing suspended cohesive sediments during the storm period. [26] Generally, high flow velocity and turbulent shear made particles/flocs small because of decreasing the flocculation rate (i.e., decreasing aggregation but increasing breakage). However, the shear-dependent flocculation was not the only factor affecting the PSD. Other physicochemical and biological factors also changed the particle/floc size and the multimodality of a PSD. For example, flocculation
10 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 6. (a) A PSD in the flood peak flow, (b) a PSD in the flood-to-ebb slack water, (c) a PSD in the ebb peak flow, and (d) a PSD in the ebb-to-flood slack water for the neap tide with a flood-dominant asymmetric tidal current. Meas represents a PSD measured with the LISST instrument. PP, Flocculi, Micro, Macro, and Sum represent the decomposed PSDs of primary particles, flocculi, microflocs and macroflocs, and the superposition of the decomposed PSDs.
11 of 17
C03014
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 7. (a) A time series of time-averaged flow velocities in x and y direction (U and V) and turbulent shear (t) and (b) a time series of the cumulative volumetric fractions of primary particles (PP), flocculi, microflocs (Micro), and macroflocs (Macro), for the period of breakage-resistant flocculation. The symbols represent the measured data, and the lines show the smoothed curves with the moving averages of three data points. In Figure 7b, each section between two lines represents the individual volume fraction of primary particles, flocculi, microflocs, and macroflocs. The four vertical lines, from the left to right, indicate the flood peak flow, the flood-to-ebb slack water, the ebb peak flow, and the ebb-to-flood slack water, of which PSDs are plotted in Figure 8 for comparison. became more breakage resistant in the spring season with decreasing floc breakage against turbulent shear and increasing particle aggregation. Besides flocculation, the previously covered silts and find sands might be eroded and advected to the measuring site by a strong wind-driven oceanic current and developed a unique PSD consisting mostly of silt-sized particles in the storm period [Baeye et al., 2011; Fettweis et al., 2012]. 3.5. Characterization and Discretization of the Multimodal PSD [27] Particle size distributions (PSDs) developed a multimodality in a flow-varying tidal cycle at the measuring site.
C03014
Based on the decomposition of the measured PSDs into a superposition of lognormal PSDs, it was concluded that the multimodal PSDs consisted of the four unimodal PSDs of primary particles, flocculi, microflocs and macroflocs that varied their dynamic behavior in a flow-varying tidal cycle. This led us to the idea that the continuous curve of a multimodal PSD might be discretized and simplified to the multiple particle/aggregate groups for practical purposes in modeling and simulation [Lee et al., 2011]. Therefore, it is necessary to review the characteristics of each discrete particle/aggregate group for reasonable discretization and simplification in modeling and simulation. The following paragraphs will review the characteristics and approximate nature of the discrete particle/aggregate groups as a replacement for a continuous multimodal PSD. [28] Primary particles have been reported to consist of various organic and inorganic particles (e.g., fragmented clay minerals, organic and calcareous particles, picophytoplankton and heterotrophic bacteria). A substantial amount of such primary particles are out-of-range fine particles with a wide size range of 0.25 2.5 mm, although they appear as a single peak at the smallest measurable size bin of the LISST instrument [Andrews et al., 2010; Fettweis et al., 2012]. Andrews et al. [2010] reported that such out-of-range particles sharply increased the volume concentration of the smallest and second smallest size bins, and irrationally outof-range particles decreased the volume concentrations of the nearby size bins up to 7 mm but increased the volume concentration of the largest size bin. This bias might further enhance the separation between the two peaks of primary particles and flocculi and develop a small peak of macroflocs in the breakage-dominant condition during the peak flows (e.g., Figures 4a and 4c). Considering the analytical error of the LISST instrument and the reported composition of out-ofrange fine particles, the actual volume concentration peak of primary particles, including out-of-range fine particles with a wide size range of 0.25 2.5 mm, might be a long extended tail of the peak of flocculi. [29] Clay minerals have been reported to aggregate into breakage-resistant flocculi and sequentially into microflocs due to the physicochemical binding mechanisms of (1) Coulombic attraction and (2) charge dilution [Mietta, 2010; Mpofu et al., 2003; Tombacz and Szekeres, 2004]. In spite of the net negative surface charge, the sporadic positive sites available on the edges of clay minerals can induce strong Coulombic attraction between clay minerals (e.g., face-to-edge attachment). Otherwise, the charge dilution by abundant salt ions in brackish and marine water, which, in aquatic chemistry, is called electric double layer compression [Stumm and Morgan, 1996], can reduce the net negative surface charge and electrostatic repulsion of clay minerals and enhance aggregation. Based on the observed difference between flocculi and microflocs with respect to structure and shear resistance, we conclude that the Coulombic attraction is likely to make strongly bound flocculi but the charge dilution favors loosely bound microflocs (Figure 1). [30] Macroflocs developed an additional peak in the PSD beyond the major peak of microflocs (e.g., Figure 8). The gluing capability of microbial exudates (EPS or TEP) is known to enhance the building of organic-rich macroflocs [Chen et al., 2005; Droppo, 2001; Droppo and Amos, 2001]. Floc sizes of marine sediments have been reported to range
12 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 8. (a) A PSD in the flood peak flow, (b) a PSD in the flood-to-ebb slack water, (c) a PSD in the ebb peak flow, and (d) a PSD in the ebb-to-flood slack water for the period of biologically mediated flocculation. Meas represents a PSD measured with the LISST instrument. PP, Flocculi, Micro, Macro, and Sum represent the decomposed PSDs of primary particles, flocculi, microflocs and macroflocs, and the superposition of the decomposed PSDs.
13 of 17
C03014
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 9. (a) A time series of time-averaged flow velocities in x and y direction (U and V) and turbulent shear (t) and (b) a time series of the cumulative volumetric fractions of primary particles (PP), flocculi, microflocs (Micro), and macroflocs (Macro), for the storm period. The symbols represent the measured data, and the lines show the smoothed curves with the moving averages of three data points. In Figure 9b, each section between two lines represents the individual volume fraction of primary particles, flocculi, microflocs, and macroflocs. The four vertical lines, from the left to right, indicate the flood peak flow, the flood-to-ebb slack water, the ebb peak flow, and the ebb-to-flood slack water, of which PSDs are plotted in Figure 10 for comparison. from hundreds to thousands of micrometers in organic-rich conditions, but only from tens to a few hundred micrometers in organic-scarce conditions [Fettweis et al., 2006]. The reported size ranges of the two floc groups in organic-scarce and organic-rich conditions agrees with the measured size ranges of microflocs and macroflocs in this research (20 200 mm and > 200 mm). This suggests that our observed microflocs and macroflocs were built by the respective physicochemical and biological flocculation mechanisms in the organic-rich seawater. In fact, the subordinate aggregate groups of primary particles, flocculi, microflocs and macroflocs have their unique assemblage mechanisms and physicochemical and biological characteristics. However, the interactions between the particle and aggregate groups and
C03014
their kinetics remain unknown and need further investigation, because they involve the complexity and heterogeneity of floc binding mechanism, composition, morphology, etc. [e.g., Maggi, 2009]. [31] We further investigated if the discrete groups of primary particles, flocculi, microflocs and macroflocs could approximate a continuous multimodal PSD in the settling flux calculation. The estimated settling fluxes of the discrete aggregate groups approached the reference settling flux of a multimodal PSD, as the number of the discrete aggregate groups increased (Figure 11 and Table 1). For example, the estimated settling fluxes with a single median diameter (D50) had the largest error of 25 50% against the reference settling flux, but the estimated settling fluxes with a volumeweighted mean diameter (Dwt) had smaller errors up to 4.5%, because the volume-weighted mean diameter takes into account the volumes of all 32 aggregate size groups of the LISST instrument. Although one could accurately estimate the sediment settling flux with a volume-weighted mean diameter, it would be very difficult to track all the 32 groups in large-scale multidimensional simulation. Therefore, decomposing a continuous multimodal PSD to the discrete aggregate groups of primary particles, flocculi, microflocs and macroflocs appears to be a practical strategy for simulating flocculation and transport in a coastal zone. [32] When calculating the estimated settling flux with the discrete aggregate groups of primary particles, flocculi, microflocs and macroflocs, we used the expected mean mean;i) that is statistically the most probable value, diameter (D g;i) obtained from instead of the geometric mean diameters (D the curve fitting analysis (equation (3)). The latter is equivalent to the median diameters (D50,i) and thus subject to erroneous estimation of the settling flux, as discussed in the g;i ) were earlier paragraph. The geometric mean diameters (D mean;i ) with converted to the expected mean diameter (D equation (3) (Weisstein, online reference, 2006). mean;i ¼ exp ln D g;i þ 1 ð lnsi Þ2 D 2
ð3Þ
[33] The estimated settling fluxes with a single discrete aggregate group had up to 45% errors against the reference settling flux of a continuous multimodal PSD. However, as the number of discrete aggregate groups increased, the estimated settling fluxes approached the reference values. The settling fluxes with two discrete groups of building blocks (i.e., primary particles and flocculi) and aggregates (i.e., microflocs and macroflocs) produced up to 9.5% errors. The settling flux with three discrete groups of building blocks, microflocs and macroflocs has up to 3.5% errors. The settling flux with four discrete groups of primary particles, flocculi, microflocs and macroflocs had up to 3.2% errors. The smallest error between the estimated and reference settling fluxes indicates that the multiple discrete aggregate groups could be a reasonable approximation of a continuous multimodal PSD.
4. Conclusion [34] The new curve fitting analysis allowed us to decompose multimodal PSDs to the subordinate unimodal lognormal
14 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Figure 10. (a) A PSD in the flood peak flow, (b) a PSD in the flood-to-ebb slack water, (c) a PSD in the ebb peak flow, and (d) a PSD in the ebb-to-flood slack water, for the storm period. Meas represents a PSD measured with the LISST instrument. PP, Flocculi, Micro, Macro, and Sum represent the decomposed PSDs of primary particles, flocculi, microflocs and macroflocs, and the superposition of the decomposed PSDs.
15 of 17
C03014
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
C03014
aggregate groups, such as well-controlled laboratory experiments, microscopic analyses for floc morphology and chemical analyses for organic or chlorophyll-A content of aggregates. The new curve fitting analysis for a large data pool of multimodal PSDs, combined with direct physicochemical and biological analyses, will allow us to investigate closely the highly complicated flocculation mechanisms in a coastal zone. [36] Acknowledgments. The authors would like to acknowledge the Flemish Science Foundation (FWO Vlaanderen) for funding the FWO project G.0263.08. This study was partly funded by the Maritime Access Division of the Ministry of the Flemish Community (MOMO project). We wish to acknowledge F. Francken and D. Van den Eynde for providing the bed shear stress data.
References
Figure 11. A comparative plot of the reference settling flux (x axis) versus the estimated settling flux (y axis) of a specific PSD. D50 and Dwt represent the estimated settling fluxes, which were calculated with a median diameter and a volume-weighted mean diameter. The 4, 3, 2, and 1 modes represent the estimated settling fluxes, which were calculated with mean diameters obtained from curve fitting analysis with 4, 3, 2, and 1 discrete aggregate group(s). PSDs and to investigate a time series of the multimodal PSDs in a qualitative and quantitative way. Results from the curve fitting analysis for a time series of multimodal PSDs agree with a general flocculation theory in that shear-dependent flocculation was a main contributor to the changing multimodality of a PSD of suspended particulate matters in a coastal zone. The PSDs at low turbulent shear skewed toward large size with a large volume fraction of aggregates in an aggregation-dominant condition, whereas the PSDs at high turbulent shear skewed toward small size with a large volume fraction of building blocks in a breakage-dominant condition. The results further revealed an important finding that flocculation became breakage-resistant against turbulent shear in the spring season. Breakage-resistant flocculation, which concurred with a temperature rise and high nutrient influx, was presumably caused by biologically mediated flocculation. This supports the findings of the recent numerical study on biological flocculation in a nutrient-enriched aquatic environment [Maggi, 2009]. However, the PSD and the constructed PSD in the storm period were a function of eroded and advected matters, irrespective of shear-dependent flocculation in a flood/ebb tidal current, as substantiated in the observations. Furthermore, the curve-fitting analysis as well as the settling flux estimation for the multimodal PSDs showed the possibility of the discrete groups of primary particles, flocculi, microflocs and macroflocs being used as an approximation of a continuous multimodal PSD in modeling and simulation. [35] However, the assemblage mechanisms and properties of primary particles, flocculi, microflocs and macroflocs were rationalized with rather indirect evidences from this research and the previous references. Therefore, we suggest for future research direct measurements of the physicochemical and biological properties of the subordinate
Agrawal, Y., and H. Pottsmith (2000), Instruments for particle size and settling velocity observations in sediment transport, Mar. Geol., 168, 89–114, doi:10.1016/S0025-3227(00)00044-X. Alldredge, A., and M. Silver (1988), Characteristics, dynamics and significance of marine snow, Prog. Oceanogr., 20, 41–82, doi:10.1016/00796611(88)90053-5. Andrews, S., D. Nover, and S. Schladow (2010), Using laser diffraction data to obtain accurate particle size distributions: The role of particle composition, Limnol. Oceanogr. Methods, 8, 507–526, doi:10.4319/ lom.2010.8.507. Baeye, M., M. Fettweis, G. Voulgaris, and V. Van Lancker (2011), Sediment mobility in response to tidal and wind-driven flows along the Belgian inner shelf, southern North Sea, Ocean Dyn., 61, 611–622, doi:10.1007/s10236-010-0370-7. Banerjee, S., B. Menkiel, and L. Ganippa (2009), Multi-lognormal soot particle size distribution for time-resolved laser induced incandescence in diesel engines, Appl. Phys. B, 96, 571–579, doi:10.1007/s00340-0093589-6. Berlamont, J., M. Ockenden, E. Toorman, and J. Winterwerp (1993), The characterization of cohesive sediment properties, Coastal Eng., 21, 105–128, doi:10.1016/0378-3839(93)90047-C. Brown, J., and A. Davies (2010), Flood/ebb tidal asymmetry in a shallow sandy estuary and the impact on net sand transport, Geomorphology, 114, 431–439, doi:10.1016/j.geomorph.2009.08.006. Chen, M., S. Wartel, and S. Temmerman (2005), Seasonal variation of floc characteristics on tidal flats, the Scheldt estuary, Hydrobiologia, 540, 181–195, doi:10.1007/s10750-004-7143-6. Christiansen, C., J. Bartholdy, and C. Sorensen (2006), Composition and size distributions of local and advected sediment trapped over a tidal flat during moderate and storm conditions, Dan. J. Geogr., 106(1), 1–11. Droppo, I. (2001), Rethinking what constitutes suspended sediment, Hydrol. Processes, 15, 1551–1564, doi:10.1002/hyp.228. Droppo, I., and C. Amos (2001), Structure, stability, and transformation of contaminated lacustrine surface fine-grained laminae, J. Sediment. Res., 71(5), 717–726, doi:10.1306/2DC40963-0E47-11D7-8643000102C1865D. Droppo, I., G. Leppard, S. Liss, and T. Milligan (2005), Flocculation in Natural and Engineered Environmental Systems, CRC Press, Boca Raton, Fla. Eisma, D. (1986), Flocculation and de-flocculation of suspended matter in estuaries, Neth. J. Sea Res., 20, 183–199, doi:10.1016/0077-7579(86) 90041-4. Fettweis, M. (2008), Uncertainty of excess density and settling velocity of mud flocs derived from in situ measurements, Estuarine Coastal Shelf Sci., 78, 426–436, doi:10.1016/j.ecss.2008.01.007. Fettweis, M., F. Francken, V. Pison, and D. Van den Eynde (2006), Suspended particulate matter dynamics and aggregate sizes in a high turbidity area, Mar. Geol., 235, 63–74, doi:10.1016/j.margeo.2006.10.005. Fettweis, M., F. Francken, D. Van den Eynde, T. Verwaest, J. Janssens, and V. Van Lancker (2010), Storm influence on SPM concentrations in a coastal turbidity maximum area with high anthropogenic impact (southern North Sea), Cont. Shelf Res., 30, 1417–1427, doi:10.1016/j.csr.2010.05. 001. Fettweis, M., M. Baeye, B. J. Lee, P. Chen, and J. C. S. Yu (2012), Hydrometeorological influences and multimodal suspended particle size distributions in the Belgian nearshore area, Geo-Mar. Lett., doi:10.1007/ s00367-011-0266-7, in press. Gibbs, R., D. Tshudy, L. Konwar, and J. Martin (1989), Coagulation and transport of sediments in the Gironde estuary, Sedimentology, 36, 987–999, doi:10.1111/j.1365-3091.1989.tb01536.x.
16 of 17
C03014
LEE ET AL.: A MULTIMODAL PSD OF COHESIVE SEDIMENTS
Hinds, W. (1999), Aerosol Technology: Properties, Behavior, and Measurement of Airborne Particles, 2nd ed., John Wiley, New York. Hussein, T., M. Dal Maso, T. Petaja, I. Koponen, P. Paatero, P. Aalto, K. Hameri, and M. Kulmala (2005), Evaluation of an automatic algorithm for fitting the particle number size distributions, Boreal Environ. Res., 10, 337–355. Latteux, B. (1995), Techniques for long-term morphological simulation under tidal action, Mar. Geol., 126, 129–141, doi:10.1016/00253227(95)00069-B. Lee, B., E. Toorman, F. Molz, and J. Wang (2011), A two-class population balance equation yielding bimodal flocculation of marine or estuarine sediments, Water Res., 45(5), 2131–2145, doi:10.1016/j.watres.2010.12.028. Li, Y., E. Wolanski, and Q. Xie (1993), Coagulation and settling of suspended sediment in the Jiaojiang River estuary, China, J. Coastal Res., 9, 390–402. Maerz, J., R. Verney, K. Wirtz, and U. Feudel (2011), Modeling flocculation processes: Intercomparison of a size class-based model and a distributionbased model, Cont. Shelf Res., 31, S84–S93, doi:10.1016/j.csr.2010.05.011. Maggi, F. (2005), Flocculation dynamics of cohesive sediment, PhD dissertation, Delft Univ. of Technol., Delft, Netherlands. Maggi, F. (2009), Biological flocculation of suspended particles in nutrientrich aqueous ecosystems, J. Hydrol., 376, 116–125, doi:10.1016/j.jhydrol. 2009.07.040. Mäkelä, J. M., I. K. Koponen, P. Aalto, and M. Kulmala (2000), One-year data of submicron size modes of tropospheric background aerosol in southern Finland, J. Aerosol Sci., 31(5), 595–611, doi:10.1016/S00218502(99)00545-5. Manning, A., and S. Bass (2006), Variability in cohesive sediment settling fluxes: Observations under different estuarine tidal conditions, Mar. Geol., 235, 177–192, doi:10.1016/j.margeo.2006.10.013. Manning, A., J. Baugh, R. Spearman, and R. Whitehouse (2010), Flocculation settling characteristics of mud:sand mixtures, Ocean Dyn., 60, 237–253, doi:10.1007/s10236-009-0251-0. Mietta, F. (2010), Evolution of the floc size distribution of cohesive sediments, PhD dissertation, Delft Univ. of Technol., Delft, Netherlands. Mikkelsen, O., P. Hill, and T. Milligan (2006), Single-grain, microfloc and macrofloc volume variations observed with a LISST-100 and a digital floc camera, J. Sea Res., 55, 87–102, doi:10.1016/j.seares.2005.09.003. Mikkelsen, O., K. Curran, P. Hill, and T. Milligan (2007), Entropy analysis of in situ particle size spectra, Estuarine Coastal Shelf Sci., 72, 615–625, doi:10.1016/j.ecss.2006.11.027. Mpofu, P., J. Addai-Mensah, and J. Ralston (2003), Influence of hydrolyzable metal ions on the interfacial chemistry, particle interactions, and dewatering behavior of kaolinite dispersions, J. Colloid Interface Sci., 261, 349–359, doi:10.1016/S0021-9797(03)00113-9. Orange, D., A. Garcia-Garcia, T. Lorenson, C. Nittrouer, T. Milligan, S. Miserocchi, L. Langone, A. Correggiari, and F. Trincardi (2005), Shallow gas and flood deposition on the Po Delta, Mar. Geol., 222–223, 159–177, doi:10.1016/j.margeo.2005.06.040. Pope, N., J. Widdows, and M. Brinsley (2006), Estimation of bed shear stress using the turbulent kinetic energy approach: A comparison of annular flume and field data, Cont. Shelf Res., 26, 959–970, doi:10.1016/j.csr.2006. 02.010. Sherwood, C., J. Lacy, and G. Voulgaris (2006), Shear velocity estimates on the inner shelf off Grays Harbor, Washington, USA, Cont. Shelf Res., 26, 1995–2018, doi:10.1016/j.csr.2006.07.025. Stumm, W., and J. Morgan (1996), Aquatic Chemistry: Chemical Equilibria and Rates in Natural Waters, 3rd ed., John Wiley, New York. Thompson, C., C. Amos, T. Jones, and J. Chaplin (2003), The manifestation of fluid-transmitted bed shear stress in a smooth annular flume: A comparison of methods, J. Coastal Res., 19, 1094–1103.
C03014
Tombacz, E., and M. Szekeres (2004), Surface charge of kaolinite in aqueous suspension in comparison with montmorillonite, Appl. Clay Sci., 27, 75–94. Trowbridge, J., and S. Elgar (2001), Turbulence measurements in the surf zone, J. Phys. Oceanogr., 31, 2403–2417, doi:10.1175/1520-0485 (2001)0312.0.CO;2. U.S. Army Corps of Engineers (2003), Surf zone hydrodynamics, in Coastal Engineering Manual, chap. 4, pp. II-4-1–II-4-40, Washington, D. C. van der Molen, J., J. Gerrits, and H. de Swart (2004), Modelling the morphodynamics of a tidal shelf sea, Cont. Shelf Res., 24, 483–507, doi:10.1016/j.csr.2003.12.001. van Ledden, M., W. van Kesteren, and J. Winterwerp (2004), A conceptual framework for the erosion behaviour of sand-mud mixtures, Cont. Shelf Res., 24, 1–11, doi:10.1016/j.csr.2003.09.002. van Leussen, W. (1994), Estuarine macroflocs: Their role in finegrained sediment transport, PhD dissertation, Utrecht Univ., Utrecht, Netherlands. Verfaillie, E., M. Van Meirvenne, and V. Van Lancker (2006), Multivariate geostatistics for the predictive modelling of the surficial sand distribution in shelf seas, Cont. Shelf Res., 26, 2454–2468, doi:10.1016/ j.csr.2006.07.028. Verney, R., J. Deloffre, J. Brun-Cottan, and R. Lafite (2007), The effect of wave-induced turbulence on intertidal mudflats: Impact of boat traffic and wind, Cont. Shelf Res., 27, 594–612, doi:10.1016/j.csr.2006.10.005. Verney, R., R. Lafite, and J. Brun-Cottan (2009), Flocculation potential of estuarine particles: The importance of environmental factors and of the spatial and seasonal variability of suspended particulate matter, Estuaries Coasts, 32, 678–693, doi:10.1007/s12237-009-9160-1. Verney, R., R. Lafite, J. C. Brun-Cottan, and P. Le Hir (2011), Behaviour of a floc population during a tidal cycle: Laboratory experiments and numerical modelling, Cont. Shelf Res., 31, S64–S83, doi:10.1016/j.csr.2010.02.005. Whitey, K. (2007), The physical characteristics of sulfur aerosols, Atmos. Environ., 41, supplement, 25–49, doi:10.1016/j.atmosenv.2007.10.057. Winterwerp, J. (1998), A simple model for turbulence induced flocculation of cohesive sediment, J. Hydraul. Eng., 36(3), 309–326, doi:10.1080/ 00221689809498621. Winterwerp, J. (2011), Fine sediment transport by tidal asymmetry in the high-concentrated Ems River: Indications for a regime shift in response to channel deepening, Ocean Dyn., 61, 203–215, doi:10.1007/s10236010-0332-0. Winterwerp, J., and W. van Kesteren (2004), Introduction to the Physics of Cohesive Sediment in the Marine Environment, Elsevier, Amsterdam. Winterwerp, J., A. Manning, C. Martens, T. de Mulder, and J. Vanlede (2006), A heuristic formula for turbulence-induced flocculation of cohesive sediment, Estuarine Coastal Shelf Sci., 68, 195–207, doi:10.1016/ j.ecss.2006.02.003. Yuan, Y., H. Wei, L. Zhao, and Y. Cao (2009), Implications of intermittent turbulent bursts for sediment resuspension in a coastal bottom boundary layer: A field study in the western Yellow Sea, China, Mar. Geol., 263, 87–96, doi:10.1016/j.margeo.2009.03.023. M. Fettweis, Management Unit of the North Sea Mathematical Models, Royal Belgian Institute of Natural Sciences, Gulledelle 100, B-1200 Brussels, Belgium. B. J. Lee and E. Toorman, Hydraulics Laboratory, Department of Civil Engineering, Katholieke Universiteit Leuven, Kasteelpark Arenberg 40, B-3001 Heverlee, Belgium. (
[email protected]) F. J. Molz, Department of Environmental Engineering and Earth Sciences, Clemson University, 342 Computer Ct., Anderson, SC 29625, USA.
17 of 17