An assessment of a new settling velocity parameterisation for cohesive sediment transport modeling

An important element within the Defra funded Estuary Process Research project “EstProc” was the implementation of the new or refined algorithms, produced under EstProc, into cohesive sediment numerical models. The implementation stage was important as any extension in the understanding of estuarine processes from EstProc was required to be suitable for dissemination into the wider research community with a level of robustness for general applications demonstrated. This report describes work undertaken to implement the new Manning Floc Settling Velocity Model, developed during EstProc. All Manning component algorithms could be combined to provide estimates of mass settling flux. The algorithms are initially assessed in a number of 1-D scenarios, where the Manning model output is compared against both real observations and the output from alternative settling parameterisations. The Manning model is then implemented into a fully 3-D computational model (TELEMAC3D) of estuarine hydraulics and sediment transport of the Lower Thames estuary. The 3-D model results with the Manning algorithm included were compared to runs with a constant settling velocity of 0.5 mm s⁻¹ and settling velocity based on a simple linear multiplier of concentration and with the above mentioned observations of suspended concentration. The findings of the 1-D case studies found the Manning empirical settling model could reproduce 93% of the total mass settling flux observed over a spring tidal cycle. The floc model fit was even better within the turbidity maximum (TM) zone. A constant 0.5 mm s⁻¹ only estimated 15% of the TM mass flux, whereas the fixed 5 mm s⁻¹ settling rate over-predicted the TM mass flux by 47%. Both settling velocity as a simple linear function of concentration, and van Leussen's method, did not fare much better estimating less than half the observed flux during the various tidal and sub-tidal cycle periods. When the Manning-settling model was applied to a layer with suspended concentrations approaching 6 g l⁻¹, it calculated 96% of the observed mass flux. The main conclusions of the implementation exercise were that it was feasible to implement a complex relationship between settling velocity and concentration in a 3-D computational model of estuarine hydraulics, without producing any significant increase in model run times or reducing model stability. The use of the Manning algorithm greatly improved the reproduction of the observed distribution of suspended concentration both in the vertical and horizontal directions compared to the other simulations. During the 1-D assessments, the Manning-settling model demonstrated flexibility in adapting to a wide range of estuarine environmental conditions (i.e. shear stress and concentration), specifically for applied modelling purposes.