Tidal Variation in Cohesive Sediment Distribution and Sensitivity to Flocculation and Bed Consolidation in An Idealized, Partially Mixed Estuary

: Particle settling velocity and erodibility are key factors that govern the transport of sediment through coastal environments including estuaries. These are di ﬃ cult to parameterize in models that represent mud, whose properties can change in response to many factors, including tidally varying suspended sediment concentration (SSC) and shear stress. Using the COAWST (Coupled Ocean-Atmosphere-Wave-Sediment Transport) model framework, we implemented bed consolidation, sediment-induced stratiﬁcation, and ﬂocculation formulations within an idealized two-dimensional domain that represented the longitudinal dimension of a micro-tidal, muddy, partially mixed estuary. Within the Estuarine Turbidity Maximum (ETM), SSC and median ﬂoc diameter varied by a factor of four over the tidal cycle. Downstream of the ETM, the median ﬂoc size and SSC were several times smaller and showed less tidal variation (~20% or less). The suspended ﬂoc distributions only reached an equilibrium size as a function of SSC and shear in the ETM at peak tidal ﬂow. In general, ﬂocculation increased particle size, which reduced SSC by half in the ETM through increased settling velocity. Consolidation also limited SSC by reduced resuspension, which then limited ﬂoc growth through reduced SSC by half outside of the ETM. Sediment-induced stratiﬁcation had negligible e ﬀ ects in the parameter space examined. E ﬀ orts to lessen the computation cost of the ﬂocculation routine by reducing the number of size classes proved di ﬃ cult; ﬂoc size distribution and SSC were sensitive to speciﬁcation of size classes by factors of 60% and 300%, respectively.


Motiviation
Estuaries are valuable coastal ecosystems that provide habitat and nursery services to many fishery species, including finfish, crustaceans, and mollusks. In estuaries, freshwater interacts with saline water and mixing can be dominated by waves, tides, riverine input or any combination [1,2]. This leads to complex hydrodynamic conditions with a broad range in spatial and temporal variability [3,4], such as those driven by tidal asymmetries, spring-neap cycles, and seasonal fluctuations in river discharge [5]. The Estuarine Turbidity Maximum (ETM) is a key feature of partially mixed estuaries that occurs at the convergence of freshwater and sea water, which can trap sediment leading to a peak in suspended sediment concentration (SSC), and moves with changes in hydrodynmic conditions [6]. Thus, understanding the processes that are critical to maintaining healthy estuarine environments is challenging.

COAWST
The Coupled Ocean-Atmospheric-Wave-Sediment Transport model (COAWST), incorporates hydrodynamics from the Regional Ocean Modeling System (ROMS) and the Community Sediment Transport Modeling System (CSTMS; [73]). ROMS is a terrain-following, free-surface, hydrostatic primative equation numerical model. ROMS simulates hydrodynamics by solving equations for continuity and conservation of mass and momentum using hydrostatic and Boussinesq assumptions with various initial and boundary conditions defined by the user. It has a large user community and has been described in detail elsewhere, see [73,77,78]. The CSTMS calculates erosion and deposition to the sediment bed and suspended transport of sediment. Sediment is removed from the bed at the user defined rate when the applied stress is greater than the critical shear stress for erosion. The transport of suspended sediment in the water column is estimated by solving the same advection-diffusion equations as used by the hydrodynamic model for salinity and temperature but with an added source/sink term to allow exchange with the sediment bed for individual sediment classes. Specifically, each sediment class has predefined hydrodynamic properties including particle density, settling velocity, and critical shear stress for erosion that have not accounted for cohesive processes until recently; for more details see [73]. Recent CSTMS developments have added cohesive modules for flocculation and bed consolidation; and its seabed layering routines now also account for cohesive, non-cohesive or mixed sediment beds [53].
The new flocculation routine (FLOCMOD) is a population balance module, and the bed consolidation routine increases the critical shear stress for erosion with depth into the sediment bed [53]. FLOCMOD requires flocculation growth and breakup parameters, primary particle size, and fractal dimension [49]. Flocculation parameters are typically chosen to follow field or laboratory studies that provide expected ranges for the associated parameters [38,79,80]. Sediment is entrained into the water column from the bed in CSTMS when the applied bed stress exceeds the critical shear stress for erosion [73,81,82]. The cohesive bed consolidation module specifies additional sediment erodibility parameters that change with depth [65]. In the York River estuary, a field study aided in defining the exponential increase of critical shear stress for erosion with depth into the bed due to consolidation [60]. However, model estimates were particularly sensitive to values that are difficult to obtain from field observations, such as the initial and equilibrium critical shear stress profiles, and the consolidation and swelling timescales [45]. The effect of density-induced stratification from gradients in SSC in CSTMS is implemented using the same approach as [73,74].
The addition of these cohesive sediment routines provides a more complete representation of the processes that are important for suspended transport of muds. The flocculation routine provides a means of examining impacts that are difficult to measure in situ. However, as yet few studies have capitalized on the implementation of a floc population balance model such as FLOCMOD within a community sediment transport model to evaluate the role of cohesive processes in muddy estuaries.

Objective and Outline of the Study
The objective of this study is therefore to evaluate the performance of the cohesive sediment routines developed in the COAWST modeling system and examine the impact of these cohesive processes on sediment distribution within an idealized, partially mixed estuarine simulation. The sensitivity of the flocculation module to the number of sediment sizes was also examined, because the computational costs increase with the number of sizes. To that end, the goal of this study is to answer the following questions: i.
Considering sediment-induced stratification, flocculation dynamics, and bed consolidation, how do these processes impact sediment distribution along a partially mixed estuarine model? ii.
How does using fewer sediment sizes constrain our ability to represent sediment dynamics in a cohesive sediment environment?
To address these research questions, an idealized model domain was used to assess the effectiveness, sensitivity to parameterization, and computational costs of the CSTMS flocculation and other cohesive sediment formulations. The simplified, idealized domain was chosen because it carries lower computational costs than the more realistic three-dimensional simulation, allowing us to complete systematic testing in a more feasible time frame. Additionally, idealized models can provide insight on first order impacts and interactions without confounding factors such as complex bathymetry and highly variable forcing functions. Section 2 details the numerical methods and sediment characteristics utilized in this study and the experiments performed. The results and discussion are presented in Section 3, and suggestions for the future handling of flocculation are provided in Section 4. Finally, the study's conclusions are outlined in Section 5.

York River Estuary, Virginia, USA
The York River estuary [17], which inspired the idealized model geometry utilitzed in this study, extends northeast from the lower Chesapeake Bay (where salinities are typically~26 g kg −1 ) to the confluence of the Pamunkey and Mattaponi rivers, with a total length of~60 km ( Figure 1). The depths in the main channel range from~20 m at the mouth to~6 m at the head of salt (~1 g kg −1 isohaline; Figure 2; [83]). The York is a microtidal, partially mixed estuary that becomes relatively more well-mixed as its depth decreases and tidal energy increases towards its head [83]. Under typical conditions, tidal current magnitude is on the order of 50 cm s −1 at 1 meter above the bed (mab) in the ETM region [17]. The sediment bed in the main channel is dominated by mud with percentages often exceeding 80% [84]. Total suspended solids within the ETM of 250 mg L −1 have been measured during slack at 1 mab [83]. For suspended flocs in the York that follow a fractial relationship in size and density distribution, [38] found the median values for the fractal dimension, primary particle size, and primary particle density to be 2.4, 1.7 microns, and 1900 kg m −3 , respectively. for the future handling of flocculation are provided in Section 4. Finally, the study's conclusions are outlined in Section 5.

York River Estuary, Virginia, USA
The York River estuary [17], which inspired the idealized model geometry utilitzed in this study, extends northeast from the lower Chesapeake Bay (where salinities are typically ~26 g kg −1 ) to the confluence of the Pamunkey and Mattaponi rivers, with a total length of ~60 km ( Figure 1). The depths in the main channel range from ~20 m at the mouth to ~6 m at the head of salt (~1 g kg −1 isohaline; Figure 2; [83]). The York is a microtidal, partially mixed estuary that becomes relatively more well-mixed as its depth decreases and tidal energy increases towards its head [83]. Under typical conditions, tidal current magnitude is on the order of 50 cm s −1 at 1 meter above the bed (mab) in the ETM region [17]. The sediment bed in the main channel is dominated by mud with percentages often exceeding 80% [84]. Total suspended solids within the ETM of 250 mg L −1 have been measured during slack at 1 mab [83]. For suspended flocs in the York that follow a fractial relationship in size and density distribution, [38] found the median values for the fractal dimension, primary particle size, and primary particle density to be 2.4, 1.7 microns, and 1900 kg m −3 , respectively.   for the future handling of flocculation are provided in Section 4. Finally, the study's conclusions are outlined in Section 5.

York River Estuary, Virginia, USA
The York River estuary [17], which inspired the idealized model geometry utilitzed in this study, extends northeast from the lower Chesapeake Bay (where salinities are typically ~26 g kg −1 ) to the confluence of the Pamunkey and Mattaponi rivers, with a total length of ~60 km ( Figure 1). The depths in the main channel range from ~20 m at the mouth to ~6 m at the head of salt (~1 g kg −1 isohaline; Figure 2; [83]). The York is a microtidal, partially mixed estuary that becomes relatively more well-mixed as its depth decreases and tidal energy increases towards its head [83]. Under typical conditions, tidal current magnitude is on the order of 50 cm s −1 at 1 meter above the bed (mab) in the ETM region [17]. The sediment bed in the main channel is dominated by mud with percentages often exceeding 80% [84]. Total suspended solids within the ETM of 250 mg L −1 have been measured during slack at 1 mab [83]. For suspended flocs in the York that follow a fractial relationship in size and density distribution, [38] found the median values for the fractal dimension, primary particle size, and primary particle density to be 2.4, 1.7 microns, and 1900 kg m −3 , respectively.

Model Description
Our implementation of COAWST used components of CSTMS [73] and hydrodynamics associated with ROMS [77,85]. Unlike most previous implementations of the CSTMS, ours utilized the flocculation dynamics module, the bed consolidation/swelling routine, and sediment-induced stratification. We neglected the impacts of sand (non-cohesive sediment) because the York River channel is dominantly mud [84]. The flocculation dynamics were represented by a set of differential equations based on population balance equations described in detail elsewhere [48,49,53]. These equations allow the exchange of sediment mass between the defined size classes, depending on the rate of turbulent shear and the sediment concentrations. Bed consolidation was accounted for by increasing the critical shear stress for erosion with depth into the sediment bed, with model parameters based on erosion experiments [60,65,86]. To implement this within the cohesive version of the CSTMS, the user must define initial and equilibrium profiles for critical shear stress for erosion. The effective instantaneous profile of the critical shear stress can be altered by erosion or deposition, but is then nudged toward the equilibrium profile to simulate consolidation or swelling [53,65]. The adjustment for sediment-induced stratification is computed as in [73].

Model Configuration
A model with the following idealized two-dimensional domain was designed to represent the longitudinal dimension of a micro-tidal, partially mixed estuary. The model represented salinity-driven estuarine circulation, but neglected across-channel variation. The primary features of the idealized simulation were similar to the main channel of the York (see Section 2.1), and the sediment characteristics were based on observed values from this system [17,38]. The idealized estuary was 180 km in total length, and the water column was partitioned into 40 vertical layers ( Figure 3). A 60-km-long area of interest was in the center of the grid, with a buffer of about 60 km on either end (please note that in the figures, "0 km" marks the seaward end of the area of interest). The horizontal resolution was 500 m, and the vertical resolution varied with depth and ranged from 0.053-0.79 m, with higher resolution near the surface and the bed. Our analysis focused on the area of interest that encompasses the estuarine system with water depths of~18 m at the seaward end to the head of salt at the landward end with a depth of 6 m. The full model grid included buffer regions to either side, to minimize the effects of boundary conditions within the area of interest. The model applied includes salinity dependence for flocculation. A natural system is likely to have different floc properties in the riverine portion due to an absence of salinity. The "riverine" section of our domain lies within the upstream buffer zone, outside our region of interest.   The tides were driven at the ocean boundary with a period of 12 h and a micro-tidal range of 0.75 m at the mouth that increased to~1 m at the head, as observed in the York River estuary [17], which produced tidal velocities similar to those observed in the York. The salinity at the seaward boundary was based on salinity measurements near the mouth of the York River estuary, which is open to the Chesapeake Bay. It was held constant in time, but varied in the vertical; specifically it was 26 psu near the bed and decreased to~14 psu at the surface. Freshwater entered the grid at the upstream river boundary with a discharge of 70 m 3 s −1 . On average, the salinity ranged between 0 and 25.5 psu in the area of interest. The temperature was held constant at 10 • C. The sediment model included 11 sediment floc size classes, with particle size diameters logarithmically spaced between 1 and 1024 µm, which represented the median diameter for each size class. These particle diameters span the range reported by [38]; from the inferred primary particle size through the largest flocs observed (~500 µm). We added an additional size class with a particle diameter of 1024 µm in order to account for larger flocs near the bed. The settling velocity (w s,i ) and density (ρ f,i ) of each class were calculated using a modified Stokes' settling equation assuming a fractal dimension of 2.4, and a primary particle density and diameter of 2000 kg m −3 and 1 µm, respectively (based on observations from the York River estuary [38]; Equations (1) and (2); Table 1). For definitions of the variables in the equations below, see Table A1 in the Appendix A.
In the equations above, ρ w is the water density, ρ p is the primary particle density, D p is the primary particle diameter, D f,i is the diameter of the floc in size class i, n f is the fractal dimension, g is the acceleration due to gravity, and ν is the kinematic viscosity. No sediment was discharged from the river, and none left the domain through the river boundary, despite the gradient boundary condition. Table 2 lists the open boundary conditions used. Other critical components to the sediment transport model included a logarithmic drag formulation with a constant hydraulic bottom roughness, for simplicity (z 0b = 5 × 10 −5 m [45,75]), and we used the k-ε turbulence closure model [87]. The model setup was designed to approximately reproduce tidal fluctuations in velocity and SSC as observed in the York River estuary [17]. A time step of 30 s and the numerical schemes listed in Table 3 were utilized.   Table 3. Numerical schemes for the idealized estuary model.

Model Experiments
All model simulations presented here used the following procedure. The model was initialized with vertically homogenous salinity contours ranging from 0-26 evenly spaced surrounding the area of interest in the model domain (−22 to 65 km) and the water column was at rest (u = 0 m s −1 ). The model was run for 90 days to allow the water column properties to adjust to the tidal conditions. The velocity and salinity fields output at the end of the 90-day run were used to initialize a 45-day run that included sediment transport and the routines for sediment-induced stratification, bed consolidation and swelling, and flocculation dynamics. In the 45-day run, the water column was initialized with no sediment in suspension, and the sediment bed was initialized with equal proportions by mass of all the floc size classes, with an initial critical shear stress of erosion profile based on data from the York River estuary taken in April 2006 following [45]. The final sediment bed characteristics, suspended sediment concentrations, velocity, and salinity fields from this 45-day run were then used as the initial conditions for a set of 30-day long case runs. Sediment was eroded from the bed in the same size class it was deposited for all runs that included flocculation dynamics. Additionally, when the bed consolidation routine was used, the consolidation timescale was 1 day, the swelling timescale was 100 days, and the equilibrium critical shear stress profile was based on data from the York River estuary taken in September 2006, also following [45]. Cases that neglected the bed consolidation routine used the default (non-cohesive) erosion formulation that calculates an active layer thickness [73]. Table 4 lists the case runs, the cohesive processes incorporated into each, and is summarized here. The reference case run incorporated the three cohesive processes described above, i.e., flocculation, bed consolidation, and sediment-induced stratification. The along estuary velocity (u) and shear rates (G) were used to show changes in hydrodynamic conditions that would transport sediment or influence floc size. The SSC, the floc size distribution in suspension and the sediment bed, and weighted settling velocity by mass were used to evaluate sediment distributions and reasoning behind changes in sediment distributions. The weighted settling velocity by mass (w s_mass ) provided insight on how quickly sediment will settle based on the amount of each floc size that was in suspension. For the reference case, the median floc diameter (D 50 ) by mass was also compared to D 50 estimated by an equilibrium theory to examine the influence of SSC and shear rate on floc size. The impact of individual cohesive processes on sediment distribution was then examined by completing additional case runs, each of which removed one of the processes from the implementation. The average depth-integrated suspended mass and average sediment size distribution from the case runs were used to assess changes in sediment distribution compared to the reference case. To evaluate the relative importance of flocculation and bed consolidation, the weighted settling velocity by mass and the average change in suspended masses for the case runs that neglected these processes were compared to estimates from the reference case. The gradient Richardson number was used to compare the suspended sediment impacts on density-induced stratification relative to salinity alone (Equations (3) and (4)).
Sensitivity to the number of sediment size classes was also evaluated. The number of sediment classes was reduced to five and three sizes from the eleven size classes in the reference case. The classes chosen for the test cases were selected in an attempt to maintain the size range while using fewer size classes. The run time, the depth averaged suspended sediment mass, the average size distribution, median size and mode, as well as the average mass weighted settling velocity were compared to reference case.
In the equations above, N is the buoyancy frequency, u is the along estuary velocity, z is the depth, g is the acceleration due to gravity, s is sediment density divided by the water density, C is SSC, and ρ w is the water density. Two Richardson numbers (Ri grad ) were calculated: one included the contribution of suspended sediment in the buoyancy frequency (N), and the other did not. See Table A1 in the Appendix A for a definition of symbols. For the version of the Ri grad that neglected sediment-induced stratification, only the second term in Equation (4) was used in the calculation of N.

Reference Case
The reference case model that incorporated cohesive sediment processes (flocculation, bed consolidation, and suspended sediment stratification) showed distinct patterns in the suspension and transport patterns linked to changes in hydrodynamics along the estuary. Compared to other locations along the estuary, the up-river boundary of salt intrusion had the strongest vertical salinity gradient, and the near-bed, tidally averaged along-estuary currents converged at this location, as expected for a partially mixed estuary. Also, this convergence trapped and accumulated sediment in the region and created an Estuary Turbidity Maximum (ETM; Figures 4a and 5a). Throughout the estuary, the suspended sediment concentration peaked near the bed and decreased with height above the bed (Figures 4a and 5a). Several processes limited the diffusion of sediment into the upper water column, including reduced mixing across the halocline, and median floc sizes whose relatively fast settling velocities hindered their upward mixing. Shear rates were greatest near the bed and near the ETM region, especially during flood tide (Figures 4b and 5b). At a height of 0.19 mab, maximum tidal shear reached 10 s −1 , which is comparable to the shear rate of 12 s −1 for maximum tidal flow produced in the laboratory for a previous zero-dimensional application of FLOCMOD [49]. Looking at the lower half of the water column where both concentrations and shear were relatively high, the larger sizes formed by aggregation were found near the bed where there was enough shear to keep them in suspension. This was evident from a larger w s_mass near the bed versus smaller settling velocities in the upper layers of the water column (Figures 4c and 5c). The bathymetric change in the upper estuary added some complexity but also relevant realism for this application relative to the CSTMS/COAWST model domains used by [87] and [53]. Specifically, the bathymetric step helped constrain the along-channel location of the maximum in tidal velocity and vertical mixing, along with the resulting salinity front and ETM, to be in the general vicinity of the transition from the deeper estuary to the shallower river, as seen in partially mixed estuaries such as the York, James, and Rappahannock ( [17,90]). Nonetheless, the resulting SSC in the reference simulation was on the same order of magnitude and the ETM was roughly the same distance from the estuary mouth as the model results from [87]. The ETM in [53] typically occurred ~20 km from the mouth and SSC was up to 400 mg L −1 higher than the concentrations shown here. These differences were not surprising, given that the present model used parameterizations chosen for the York River. Specifically, compared to [53], we used a larger fractal dimension, a smaller primary particle size, and lower floc particle densities. This study also used a different equilibrium critical shear stress for erosion, and longer consolidation and swelling timescales for the bed consolidation parameterization.  The median floc sizes, D50, predicted from FLOCMOD within the ETM were interpreted in the context of an equilibrium floc size theory that postulates that a floc distribution at equilibrium will exhibit a linear relationship between D50 and the ratio of SSC by mass (C) to the square root of the shear rate (G) [46,53]. To explore this with our model results, the ETM D50 was calculated as the diameter at which 50% of the suspended flocs were smaller by mass and compared to C/G 1/2 (Figure The bathymetric change in the upper estuary added some complexity but also relevant realism for this application relative to the CSTMS/COAWST model domains used by [87] and [53]. Specifically, the bathymetric step helped constrain the along-channel location of the maximum in tidal velocity and vertical mixing, along with the resulting salinity front and ETM, to be in the general vicinity of the transition from the deeper estuary to the shallower river, as seen in partially mixed estuaries such as the York, James, and Rappahannock ( [17,90]). Nonetheless, the resulting SSC in the reference simulation was on the same order of magnitude and the ETM was roughly the same distance from the estuary mouth as the model results from [87]. The ETM in [53] typically occurred~20 km from the mouth and SSC was up to 400 mg L −1 higher than the concentrations shown here. These differences were not surprising, given that the present model used parameterizations chosen for the York River. Specifically, compared to [53], we used a larger fractal dimension, a smaller primary particle size, and lower floc particle densities. This study also used a different equilibrium critical shear stress for erosion, and longer consolidation and swelling timescales for the bed consolidation parameterization.
The median floc sizes, D 50 , predicted from FLOCMOD within the ETM were interpreted in the context of an equilibrium floc size theory that postulates that a floc distribution at equilibrium will exhibit a linear relationship between D 50 and the ratio of SSC by mass (C) to the square root of the shear rate (G) [46,53]. To explore this with our model results, the ETM D 50 was calculated as the diameter at which 50% of the suspended flocs were smaller by mass and compared to C/G 1/2 ( Figure 6). Hourly data for the last ten days of the model run were analyzed to avoid effects of model spin-up. Throughout much of the water column in the ETM, model results for periods with relatively high shear rates (i.e., between 4-12 s −1 ) did indeed display a linear relationship between D 50 and C/G 1/2 consistent with an equilibrium floc distribution ( Figure 6). Data from 0.70 mab and higher in the water column, during times with shear rates between 4 s −1 and 12 s −1 , and SSC exceeding 100 mg L −1 were used to fit an equilibrium line (inset in Figure 6). These conditions were assumed to be mostly likely at equilibrium.  . The ratio of SSC (C) to the square root of the shear rate (√G) versus the D50 by mass for the last 10 days of the reference case for the ETM. The colors correspond to depth in the water column (blue being near the bed; red being near the water surface). The blue triangles represent the first grid cell above the sediment-water interface (~0.036 mab). The red line represents the best fit during times with shear rates between 4-12 s −1 ; the slope and intercept are provided in the top left corner. The inset in the top right shows the subset of data used for the equilibrium regression.
Next, we explored time-dependent variability within the ETM region for the reference case starting with the hydrodynamic and sediment conditions nearest the seafloor (0.036 mab). The tidal velocities reached around 40 cm s −1 during peak flood, when the near-bed conditions for the shear rates were about 225 s −1 , and SSCs were ~800 mg L −1 (Figure 7a-c). Both the peaks in the shear and SSC corresponded with the velocity maximum (Figure 7a-c) and the larger particles in suspension, with the very largest seen at peak flood (Figure 7d). An examination of floc size distributions ( Figure  7d,g) revealed that during peak tidal flow, the median floc size, at ~190 μm, exceeded the equilibrium and there were larger floc sizes on the sediment bed surface as finer grains had been winnowed (Figure 7f). The flow velocity and shear was again elevated during ebb flow but remained less energetic than during flood. These currents were strong enough to erode sediment from the bed, and At times when the model produced a D 50 that fell well above or below the equilibrium line in Figure 6, the modeled floc size was interpreted as larger or smaller than the equilibrium, which indicated that another process was likely influencing the size of the flocs in suspension. The median floc size in most of the water column alternated between lying near the predicted equilibrium size to being finer-grained than the equilibrium. In contrast, the model output that most consistently fell above the equilibrium line was for the very near-bed grid cell at 0.036 mab (triangles in Figure 6). These deviations from equilibrium were caused by non-local process such as settling from the upper water column and sediment input via seabed resuspension or advection; or local processes that lead to relatively long time scales for aggregation and disaggregation.
Next, we explored time-dependent variability within the ETM region for the reference case starting with the hydrodynamic and sediment conditions nearest the seafloor (0.036 mab). The tidal velocities reached around 40 cm s −1 during peak flood, when the near-bed conditions for the shear rates were about 225 s −1 , and SSCs were~800 mg L −1 (Figure 7a-c). Both the peaks in the shear and SSC corresponded with the velocity maximum (Figure 7a-c) and the larger particles in suspension, with the very largest seen at peak flood (Figure 7d). An examination of floc size distributions (Figure 7d,g) revealed that during peak tidal flow, the median floc size, at~190 µm, exceeded the equilibrium and there were larger floc sizes on the sediment bed surface as finer grains had been winnowed (Figure 7f). The flow velocity and shear was again elevated during ebb flow but remained less energetic than during flood. These currents were strong enough to erode sediment from the bed, and produced an ETM floc size distribution that included larger size classes as during flood flow, but were slightly smaller (not shown). The larger particles in the near-bed suspension during peak flows may be attributed in part to erosion of larger aggregates from the bed, as well as enhanced mixing and elevated SSC that promoted floc growth (Figure 7). However, during slack flow the near-bed D 50 equaled the equilibrium at~70 µm (Figure 7g). As velocity decreased the larger, faster-settling floc sizes settled, leaving smaller sizes in suspension. Also, the sediment bed became enriched in the fractions of smaller floc sizes, as bed stresses decreased and the medium fractions could settle to the bed without being resuspended (Figure 7i). A two-dimensional idealized model utilizing the Chesapeake Bay ETM as a prototype predicted comparable near-bed SSC in the ETM of~200 mg L −1 near slack and 400 mg L -1 during flood flow, and the median floc size during flood was~150 µm, but the largest flocs were predicted during slack flow [91].
Higher in the water column (0.90 mab) at the ETM, the fluctuations in velocity were larger, but the shear and SSCs were smaller compared to near the bed (Figure 7a-c). The floc size distribution had a smaller concentration of the mid-sized sediment classes, which produced a larger D 50 than near the bed. Unlike the near-bed, D 50 at 0.90 mab equaled the predicted equilibrium D 50 of~256 µm at peak flow (Figure 7e), and was smaller at~90 µm than equilibrium at slack water (Figure 7h). At slack flow, the D 50 at 0.90 mab was similar to the value in the near-bed grid cell, as larger floc sizes had settled. The SSC at 0.90 mab was on average~2 times smaller than at 0.036 mab (Figure 6c), whereas G 1/2 at 0.90 mab was 7 times smaller than at 0.036 mab (Figure 7b). Thus, differences in C/G 1/2 and the equilibrium D 50 very near the bed versus higher in the water column were driven more by differences in shear rate than by differences in SSC. Overall, the ETM concentrations and floc sizes simulated in this study were smaller than those predicted at 0.80 mab in the idealized ETM model for the Chesapeake Bay [91].
The conditions in the lower reaches of the estuary differed from the ETM in that the tidal velocities, the square root of the shear, and SSC respectively stayed below 40 cm s −1 , ranged between 2-10 s −1/2 and ranged between 50-90 mg L −1 in the bottom 1.75 m (Figure 8a-c). These SSC values were similar to concentrations predicted in the downstream portion of an idealized estuary in [87]. The tidal fluctuation in velocity was greater at 1.75 mab than directly above the bed, but the shear rate showed greater tidal fluctuation and magnitude near the bed (Figure 8a,b). The SSC, the floc distribution, and D 50 showed minimal tidal variation in the bottom 1.75 m (Figure 8c-e,g,h). Near the bed (~0.078 mab), the floc size was equal or near the predicted equilibrium floc sizes both near slack and peak flood flow. Higher in the water column (1.75 mab), the modeled median diameters were smaller than the equilibrium floc sizes. The SSC to the square root of shear rate ratio was larger higher in the water column, leading to a larger predicted equilibrium floc size than was calculated by the time-dependent model. Case studies from [92] also showed situations where the floc size was over predicted by the [46] model. The floc size within the sediment bed downstream did not change over the tidal cycle (Figure 8i). equilibrium D50 very near the bed versus higher in the water column were driven more by differences in shear rate than by differences in SSC. Overall, the ETM concentrations and floc sizes simulated in this study were smaller than those predicted at 0.80 mab in the idealized ETM model for the Chesapeake Bay [91]. The conditions in the lower reaches of the estuary differed from the ETM in that the tidal velocities, the square root of the shear, and SSC respectively stayed below 40 cm s −1 , ranged between 2-10 s −1/2 and ranged between 50-90 mg L −1 in the bottom 1.75 m (Figure 8a-c). These SSC values were similar to concentrations predicted in the downstream portion of an idealized estuary in [87]. The tidal fluctuation in velocity was greater at 1.75 mab than directly above the bed, but the shear rate showed greater tidal fluctuation and magnitude near the bed (Figure 8a,b). The SSC, the floc distribution, and D50 showed minimal tidal variation in the bottom 1.75 m (Figure 8c-e,g,h). Near the bed (~0.078 mab), the floc size was equal or near the predicted equilibrium floc sizes both near slack and peak flood flow. Higher in the water column (1.75 mab), the modeled median diameters were smaller than the equilibrium floc sizes. The SSC to the square root of shear rate ratio was larger higher in the water column, leading to a larger predicted equilibrium floc size than was calculated by the time-dependent model. Case studies from [92] also showed situations where the floc size was over predicted by the [46] model. The floc size within the sediment bed downstream did not change over the tidal cycle (Figure 8i).  The results from the reference case demonstrate the variability of floc size under simplified tidal conditions, specifically within the ETM region. The results highlight that the assumption of an equilibrium floc size may not be valid here, because suspended floc sizes within the ETM often deviated from the equilibrium values. As also shown, the floc size predicted by the flocculation routine for the water column outside the ETM was fairly constant though it did not equal the equilibrium floc size. Therefore in the lower reaches of the estuary, a simulation with an The results from the reference case demonstrate the variability of floc size under simplified tidal conditions, specifically within the ETM region. The results highlight that the assumption of an equilibrium floc size may not be valid here, because suspended floc sizes within the ETM often deviated from the equilibrium values. As also shown, the floc size predicted by the flocculation routine for the water column outside the ETM was fairly constant though it did not equal the equilibrium floc size. Therefore in the lower reaches of the estuary, a simulation with an appropriately chosen constant floc size and settling velocity could likely yield similar predictions of SSC as the present flocculation routine.

Sensitivity Tests
The impacts of flocculation, bed consolidation, and suspended sediment induced stratification to the distribution of sediment in the idealized estuary was further explored by comparing the results from the reference case to the test cases that neglected each of these processes (Table 4).

Impacts of Flocculation vs. No Flocculation
The role of flocculation was evaluated by comparing results from the reference case, which included flocculation to the "no floc" case (Table 4). Including flocculation dynamics shifted the size distribution of suspended matter to the larger size classes relative to the no flocculation case, as seen in the size distribution for the ETM (Figure 9a,b). The time-averaged D 50 by mass in the ETM size distribution increased to 123 µm with flocculation, compared to 55 µm in the no flocculation case. Overall, the effect of flocculation was most significant in the upstream portion of the estuary, where flocculation in the reference run decreased the depth-integrated suspended mass by~2.5 times in the ETM, compared to the no flocculation run (Figure 10). The size of sediment in suspension increased even though flocculation reduced the SSC, producing a higher w s_mass in the ETM relative to the no flocculation run (Figure 11a,b). Flocculation increased particle size and as a result increased the settling velocities of the suspended sediment in the ETM. This limited the transport of sediment out of the ETM, and increased deposition there. ETM, compared to the no flocculation run ( Figure 10). The size of sediment in suspension increased even though flocculation reduced the SSC, producing a higher ws_mass in the ETM relative to the no flocculation run (Figure 11a,b). Flocculation increased particle size and as a result increased the settling velocities of the suspended sediment in the ETM. This limited the transport of sediment out of the ETM, and increased deposition there.     For individual size classes, the difference in suspended mass between the reference and no flocculation cases showed that in the ETM (found at ~45-55 km), flocculation moved mass from the smaller size classes ≤64 μm into larger classes ≥256 μm (Figure 12). In the lower estuary (≤ ~30 km), including flocculation generally removed mass from the smallest size classes of 1-16 μm in diameter (blue in Figure 12a,b); but increased suspended mass in the median sizes (red in Figure 12c,d). These larger particles with higher settling velocities concentrated closer to the bed, so that flocculation reduced SSC in the surface waters of the lower estuary ( Figure 12). For individual size classes, the difference in suspended mass between the reference and no flocculation cases showed that in the ETM (found at~45-55 km), flocculation moved mass from the smaller size classes ≤64 µm into larger classes ≥256 µm (Figure 12). In the lower estuary (≤~30 km), including flocculation generally removed mass from the smallest size classes of 1-16 µm in diameter (blue in Figure 12a,b); but increased suspended mass in the median sizes (red in Figure 12c,d). These larger particles with higher settling velocities concentrated closer to the bed, so that flocculation reduced SSC in the surface waters of the lower estuary ( Figure 12). flocculation cases showed that in the ETM (found at ~45-55 km), flocculation moved mass from the smaller size classes ≤64 μm into larger classes ≥256 μm (Figure 12). In the lower estuary (≤ ~30 km), including flocculation generally removed mass from the smallest size classes of 1-16 μm in diameter (blue in Figure 12a,b); but increased suspended mass in the median sizes (red in Figure 12c,d). These larger particles with higher settling velocities concentrated closer to the bed, so that flocculation reduced SSC in the surface waters of the lower estuary ( Figure 12).

Impacts of Bed Consolidation vs. No Consolidation
The time-averaged size distribution from the ETM when bed consolidation was included showed the D 50 by mass was reduced to 123 µm, compared to 170 µm when bed consolidation was neglected (Figure 9c). The depth-integrated suspended mass was reduced throughout most of the estuary; ETM concentrations were~1.5 times smaller in the reference case than those calculated in the no consolidation case (Figure 10). The settling velocities (w s_mass ) were lower when bed consolidation was included compared to the case that neglected bed consolidation (Figure 11a,c). The reference run with consolidation showed that relative to no consolidation, suspended mass was removed from the larger size classes (128-1024 µm), and suspended mass was added to the smaller sizes in most of the estuary (Figure 13). Throughout the area of interest, bed consolidation reduced bed erosion, lowering the SSC and the probability of flocculation to occur, which shifted the size distribution to smaller sizes compared to the no consolidation case (Figure 13).
was included compared to the case that neglected bed consolidation (Figure 11a,c). The reference run with consolidation showed that relative to no consolidation, suspended mass was removed from the larger size classes (128-1024 μm), and suspended mass was added to the smaller sizes in most of the estuary (Figure 13). Throughout the area of interest, bed consolidation reduced bed erosion, lowering the SSC and the probability of flocculation to occur, which shifted the size distribution to smaller sizes compared to the no consolidation case (Figure 13).

Impacts of Sediment-Induced Stratification
Density-induced stratification due to gradients in SSC did not significantly impact sediment patterns; the size distribution of suspended sediment was similar to the distribution calculated when sediment-induced stratification was removed, the time averaged D50 was 123 μm for the reference case compared to 116 μm for the case that neglected sediment-induced stratification (Figures 9d). The depth-integrated suspended sediment mass throughout the estuary for the reference case was nearly identical to the no stratification case (Figure 10). These results indicate that the SSC gradients

Impacts of Sediment-Induced Stratification
Density-induced stratification due to gradients in SSC did not significantly impact sediment patterns; the size distribution of suspended sediment was similar to the distribution calculated when sediment-induced stratification was removed, the time averaged D 50 was 123 µm for the reference case compared to 116 µm for the case that neglected sediment-induced stratification (Figure 9d). The depth-integrated suspended sediment mass throughout the estuary for the reference case was nearly identical to the no stratification case ( Figure 10). These results indicate that the SSC gradients estimated in this set of model runs were not sufficient to induce significant stratification. The gradient Richardson number (Ri grad ; Equations (3) and (4); [93]) was calculated for 0.5 and 1.0 mab throughout the area of interest. Suspended sediment on average contributed minimally to the Ri grad . When the Ri grad exceeded the critical value of 0.25, SSC contributed 2.5% to the Ri grad at 0.5 mab, and 1.8% at 1 mab. Please note that the density of the flocs decreased with size, as was observed in the York River estuary (Table 1; [38]). Suspended sediment-induced stratification could be more important in situations where the density of large flocs did not decrease, where higher stresses or higher erodibility led to higher SSC near the bed, or where stronger temperature/salinity stratification more significantly enhanced sediment-induced stratification.
To summarize Section 3.2.1, Section 3.2.2, Section 3.2.3, flocculation dynamics were especially important in the ETM, while bed consolidation impacted suspended sediment throughout the estuary. Based on these results, future modeling studies in cohesive environments should consider using flocculation if variability in SSC and shear rate is comparable or stronger than the ETM modeled here. However, changes in erodibility due to consolidation is important for most studies in regions with muddy or mixed-sediment beds. For this implementation, density-induced stratification due to SSC gradients was less important for suspended sediment dynamics. The SSCs in this study were <1000 mg L −1 and floc densities were ≤2000 kg m −3 ; thus sediment-induced stratification could be neglected in studies with similar sediment dynamics. However, if a user is unsure, it is best to incorporate density-induced stratification from sediment, as this routine has a low computational cost in COAWST and may impact sediment distributions and mixing.

Sensitivity to the Floc Size Distribution
The final set of cases evaluated the sensitivity of model estimates to the number of size classes used to represent the floc population. Results showed that the calculations of SSC were sensitive to the number of size classes used. Reducing the number of sediment sizes from eleven classes in the reference case to five or three increased the amount of suspended sediment throughout most of the estuary (Figure 14). The case using only three size classes produced depth-integrated suspended sediment that exceeded those from the case that neglected flocculation; the peak concentration in the three size class case was approximately 30% higher than the no flocculation case (Figure 14). led to higher SSC near the bed, or where stronger temperature/salinity stratification more significantly enhanced sediment-induced stratification.
To summarize Sections 3.2.1-3.2.3, flocculation dynamics were especially important in the ETM, while bed consolidation impacted suspended sediment throughout the estuary. Based on these results, future modeling studies in cohesive environments should consider using flocculation if variability in SSC and shear rate is comparable or stronger than the ETM modeled here. However, changes in erodibility due to consolidation is important for most studies in regions with muddy or mixed-sediment beds. For this implementation, density-induced stratification due to SSC gradients was less important for suspended sediment dynamics. The SSCs in this study were <1000 mg L −1 and floc densities were ≤2000 kg m −3 ; thus sediment-induced stratification could be neglected in studies with similar sediment dynamics. However, if a user is unsure, it is best to incorporate densityinduced stratification from sediment, as this routine has a low computational cost in COAWST and may impact sediment distributions and mixing.

Sensitivity to the Floc Size Distribution
The final set of cases evaluated the sensitivity of model estimates to the number of size classes used to represent the floc population. Results showed that the calculations of SSC were sensitive to the number of size classes used. Reducing the number of sediment sizes from eleven classes in the reference case to five or three increased the amount of suspended sediment throughout most of the estuary (Figure 14). The case using only three size classes produced depth-integrated suspended sediment that exceeded those from the case that neglected flocculation; the peak concentration in the three size class case was approximately 30% higher than the no flocculation case (Figure 14). The modal and D50 floc size in suspension by mass were sensitive to the resolution of the floc size distribution. When five size classes were used, suspended mass near the bed in the ETM had a modal size of 64 μm, whereas the reference run had a modal size diameter of 128 μm (Figure 15a,b). The D50 by mass of the five-size class run was 74 μm compared to the reference case D50 of 123 μm. The modal and D 50 floc size in suspension by mass were sensitive to the resolution of the floc size distribution. When five size classes were used, suspended mass near the bed in the ETM had a modal size of 64 µm, whereas the reference run had a modal size diameter of 128 µm (Figure 15a,b). The D 50 by mass of the five-size class run was 74 µm compared to the reference case D 50 of 123 µm. For the three-size class run, the floc D 50 near the bed in the ETM was similar to the 5-size class run at 72 µm. However, the modal size was similar to the reference case at 128 µm, but there was much more suspended mass for the three-class case compared to the reference case ( Figure 15). The changes to suspended size distribution likewise reduced the weighted setting velocities calculated for the three and five class cases, compared to the reference case ( Figure 16). The case with only three size classes had the smallest w s_mass of the cases considered (Figure 16c).
For the three-size class run, the floc D50 near the bed in the ETM was similar to the 5-size class run at 72 μm. However, the modal size was similar to the reference case at 128 μm, but there was much more suspended mass for the three-class case compared to the reference case ( Figure 15). The changes to suspended size distribution likewise reduced the weighted setting velocities calculated for the three and five class cases, compared to the reference case ( Figure 16). The case with only three size classes had the smallest ws_mass of the cases considered (Figure 16c).  The computational expense of the flocculation module was related to the number of size classes used. The maximum number of size classes used in this study was eleven because the computational expense to use more was prohibitive. Including flocculation dynamics in the reference model that had eleven sediment classes increased the computational time by a factor of ~130 (run time of ~128 h with flocculation and ~1 h without). Using five or three size classes increased the cost by factors of ~31, and ~13.4, respectively, compared to not using flocculation dynamics. Finding the optimal set of floc sizes may take several test runs to determine the best configuration as shown in previous studies [49]. Size distributions from our test cases indicated that representing the particles that are abundantly found in suspension (i.e., sizes 16-512 μm in our idealized estuary) is more important than maintaining the complete possible size range. This set of sensitivity shows that to apply the For the three-size class run, the floc D50 near the bed in the ETM was similar to the 5-size class run at 72 μm. However, the modal size was similar to the reference case at 128 μm, but there was much more suspended mass for the three-class case compared to the reference case ( Figure 15). The changes to suspended size distribution likewise reduced the weighted setting velocities calculated for the three and five class cases, compared to the reference case ( Figure 16). The case with only three size classes had the smallest ws_mass of the cases considered (Figure 16c).  The computational expense of the flocculation module was related to the number of size classes used. The maximum number of size classes used in this study was eleven because the computational expense to use more was prohibitive. Including flocculation dynamics in the reference model that had eleven sediment classes increased the computational time by a factor of ~130 (run time of ~128 h with flocculation and ~1 h without). Using five or three size classes increased the cost by factors of ~31, and ~13.4, respectively, compared to not using flocculation dynamics. Finding the optimal set of floc sizes may take several test runs to determine the best configuration as shown in previous studies [49]. Size distributions from our test cases indicated that representing the particles that are abundantly found in suspension (i.e., sizes 16-512 μm in our idealized estuary) is more important than maintaining the complete possible size range. This set of sensitivity shows that to apply the The computational expense of the flocculation module was related to the number of size classes used. The maximum number of size classes used in this study was eleven because the computational expense to use more was prohibitive. Including flocculation dynamics in the reference model that had eleven sediment classes increased the computational time by a factor of~130 (run time of~128 h with flocculation and~1 h without). Using five or three size classes increased the cost by factors of~31, and~13.4, respectively, compared to not using flocculation dynamics. Finding the optimal set of floc sizes may take several test runs to determine the best configuration as shown in previous studies [49]. Size distributions from our test cases indicated that representing the particles that are abundantly found in suspension (i.e., sizes 16-512 µm in our idealized estuary) is more important than maintaining the complete possible size range. This set of sensitivity shows that to apply the flocculation model in other environments, using a one-or two-dimensional domain would be useful for designing an appropriate floc distribution and optimal number of sediment classes to represent. Additionally, the availability of observation-based detailed information on the floc size distribution is important for designing the size distributions for FLOCMOD. The flocculation routine was found to be sensitive to the characteristics assigned to size classes, therefore care should be taken in assigning floc properties. For these reasons, it remains challenging to select an optimal size distribution and associated sediment characteristics of sediment size classes that optimize the results while limiting computational costs.

Key Implications and Future Directions
The case presented in this study showed that sediment size distributions vary along a representative muddy, partially mixed estuary. Our results showed that the processes of bed consolidation and flocculation modified the sediment erodibility and settling velocity within the idealized estuary over tidal timescales. However, many sediment transport models hold these sediment characteristics constant when modeling an entire estuarine system, which can lead to under-or overestimates of SSC. The cohesive process sensitivity tests indicated that it was important to include bed consolidation in estimating the SSC along the estuary. Flocculation was important in areas of the estuary with tidal variability, that has a range similar to or greater than 300 mg L −1 and 8 s −1 in SSC and shear rate, respectively, as was the case for 0.90 mab in the ETM of this idealized estuary. However, due to the sensitivity of the flocculation routine to the size class distribution and other sediment characteristics, observational data are still needed to appropriately constrain the resulting SSC. Additional modeling studies are warranted to help determine the optimal sediment parameters for FLOCMOD. Since the flocculation routine is computationally expensive, a set of one-dimensional runs that more systematically vary the sediment size classes and properties may be useful for parameterizing FLOCMOD for natural estuaries.
There are motivations for using FLOCMOD in regional scale domains despite the computational expense. Employing FLOCMOD in regional studies is now relatively straightforward, because the module has been coupled with the CSTMS [53]. The population balance equations fit well within the CSTMS piecewise parabolic settling scheme [73], chosen to reduce limitations by the CFL (Courant Friedrichs Lewy) stability criteria. Coupling less computationally expensive flocculation routines, such as those that track changes in effective floc settling velocity [46,94], would require structural change to the vertical advection scheme used within the CSTMS, which assumes that the settling velocity of each particle class remains constant throughout the water column.
One approach to assessing the impact of flocculation for a location while maintaining a modest computational expense is to apply FLOCMOD within a small one-dimensional vertical domain, with realistic physics and bed dynamics. This could be a reasonable method for evaluating flocculation dynamics for an area of interest such as the York River estuary ETM. In the York River estuary, flocculation dynamics seem more important within regions with elevated SSC, compared to the rest of the estuary. Further work on reducing the computational cost of the flocculation routine by reducing the number of vertical water column layers, or modifying the FLOCMOD convergence criteria, could make a three-dimension model more cost effective and provide added insight on the distribution of flocs under more realistic and complex bathymetric conditions.
The flocculation formulation may also be useful for examining the formation, dispersal, and fate of oil-mineral aggregates. Parameterization of FLOCMOD could be modified to represent the properties of oil-mineral aggregates using data from laboratory experiments. For example, the growth and breakup parameters in Equations (7)-(12) from [49] may be modified to account for hydrocarbon content. This would allow the sediment transport routines within the CSTMS to be used to examine the dispersal of oil from a spill. However, this does not account for the impact of the bacterial communities that are attracted to oil-mineral aggregates and breakdown the oil, and excrete exopolymeric substances (EPS; [95]), which potentially changes the characteristics of the settling aggregates. Linking FLOCMOD to the HydroBioSed CSTMS module that has been used to account for biogeochemistry [96] may provide an avenue for incorporating microbial degradation with the flocculation process. This would be useful for tracking the formation, transformation, and transport of oil-mineral aggregates.

Conclusions
The idealized model developed in this study reproduced the key physical features of an estuary such as estuarine circulation and the formation of an ETM through flow convergence and sediment trapping. The flocculation model adequately represented the dynamics throughout the entire estuary, and the results demonstrated that sediment conditions in the ETM were more variable over the tidal cycle with SSC~3-10 times greater and median diameters~2-8.5 times larger than downstream. In the ETM, the distribution of floc sizes fluctuated with the shear rate and SSC; the largest flocs in suspension were seen during flood tide, and these tended to deposit to the bed during slack flow. An equilibrium floc theory was applied to the modeled conditions that were assumed to be near-equilibrium (ETM, above the bed, and peak flow conditions) to derive a relationship between equilibrium floc diameter, SSC, and the shear rate. The analysis indicated much of the estuary was not at equilibrium but that other processes were impacting the floc size over much of the tidal cycle throughout the estuary. Additional analysis is needed to determine which other processes are influencing floc size and at which stage of the tidal cycle and where in the estuary are these other processes more dominate.
The impact of individual processes and user-defined floc size characteristics were evaluated through a series of case studies that individually excluded bed consolidation, flocculation, or sediment-induced stratification, and varied the number of floc sizes used. Flocculation dynamics had the largest impact on SSC within the ETM region where both suspended sediment and turbulence varied in time and space. Flocculation reduced the average depth-integrated suspended mass by~50% within the ETM. Outside of the ETM, bed consolidation had the largest impact and similarly decreased the average depth-integrated suspended mass by~50%. Sediment-induced density stratification had a much smaller impact on the distribution of sediment than either flocculation or bed consolidation. The vertical gradients in the suspended concentrations within the idealized estuary did not inhibit mixing based on the SSC contributions to the gradient Richardson number. The use of fewer size classes, to reduce computational expense, resulted in estimates of sediment concentration that increased by as much as~3.5 times, and an estimated median floc size that decreased by as much as~1.7 times relative to the reference case. The flocculation routine was sensitive to the number of floc size classes used and is also expected to be sensitive to the characteristics (i.e., settling velocity, critical shear stress) assumed for the floc classes.
The reference run in this study, though it represented an idealized estuarine configuration, was parameterized using observational data from the York River estuary. Other applications of the cohesive sediment model for a different location would require some insight into how to effectively parameterize these routines. The use of both observed SSCs and floc sizes helped us to constrain model parameters because multiple possible model configurations can yield similar SSC. Measurements of the suspended sediment size distribution based on in situ camera systems, laser instruments, and/or laboratory experiments can be especially useful for constraining fractal dimension and primary particle size and density [25,38,40,[97][98][99]. Erodibility estimates obtained under variable environmental conditions, either through erosion experiments or calibrated estimates from a high sample frequency acoustic instrument [60,100] can provide constraints for the bed consolidation parameters.
Our study showed that sediment transport calculations in partially mixed, muddy estuaries may be equally sensitive to incorporating flocculation and bed consolidation. However, the inclusion of flocculation using a population-based floc model added significantly to the computational cost of our coupled sediment transport-hydrodynamic model. This may prove challenging for applying the population-based flocculation model to a regional-scale model domain, but the aggregation and disaggregation processes significantly impacted suspended sediment concentrations, and suspended grain size distribution.