Three-dimensional modeling of sediment resuspension in a large shallow lake

We simulated bottom resuspension events in Lake Erie, using a coupled three-dimensional hydrodynamic and water quality model. Key parameters in the model, including critical bottom shear stress (tau(cr)) and resuspension rate (alpha) were calibrated and validated by comparing the model output to observations. These included total suspended solid (TSS) concentrations in the bottom boundary layer (RMSE = 0.74 mg L-1) and water column (RMSE = 0.81 mg L-1), and to time series of acoustic backscatter signal (R-2 > 0.8) and turbidity (R-2 approximate to 0.4) from long-term moorings near the lakebed in 2008-09 and 2013. Signals from phytoplankton, in spring and summer, caused discrepancies between modeled TSS and the observed turbidity data. Although common practice, we show that literature-based or field-observed critical shear stress should not be directly applied in large-scale Reynolds-averaged sediment model as this will likely underestimate resuspension. In agreement with the literature, the model reproduced more frequent and intensive surface-wave driven resuspension in the shallow regions (< similar to 20 m), particularly in the western basin, compared to the deeper central and eastern basins, where bottom stresses induced by mean currents (tau(c)) were comparable with those due to surface waves (tau(w)). However, on the northshore of the eastern basin, tau(c) often predominated over tau(w). We simulated thermocline motion, including up- and down-welling events and swashing of the internal Poincare wave, to contribute to tau(c) in the central basin and form nepheloid layers. (C) 2021 International Association for Great Lakes Research. Published by Elsevier B.V. All rights reserved.