A vertical gradient of suspended sediment concentration often exists in estuaries, particularly within the bottom boundary layer where sediment erosion and deposition take place. This results in a vertical density gradient and hence modification of the flow. However, this important effect has often been ignored in past sediment studies. Because of this and because of other empirical assumptions, existing erosion models cannot be used as predictive tools. This paper employs a simplified second-order closure model to simulate the effect of sediment-induced stratification on bottom boundary layer dynamics, and particularly the erosion process. Numerical models, which employ the concept of Richardson number dependent eddy viscosity, have been developed in the past for stratified flows. These models require a large number of data for tuning parameters of the eddy viscosity formula. The simplified second-order closure model used here consists of the dynamic equations of motion for mean variables (velocity, temperature, salinity, density, and suspended sediment concentration) and turbulent quantities of turbulent kinetic energy and turbulence macroscale. Model constants are thus invariant owing to the added physcis. Model simulations of laboratory sediment-laden boundary layers indicate that, as sediments are eroded or resuspended from the bottom, a vertical gradient of suspended sediment concentration can lead to a significant reduction of the turbulent shear stresses and a slowdown of the erosion-resuspension process. Hence significant error may be contained in some empirical sediment erosion rate formulae, which were derived by ignoring flow-sediment interactions. Simulations of wave boundary layer and a thermally stratified boundary layer, obtained with the same numerical model, are also presented. There is good agreement between measured and simulated mean and turbulent quantities in both cases. ¿ American Geophysical Union 1989