A method for estimating a suspended sediment distribution in an ocean region. The method includes receiving time-dependent variables including wave height, wave period, wave propagation direction, steady current magnitude, steady current direction, and water depth, and an estimate of bottom sediment distribution, calculating a new bottom sediment distribution using a sediment suspension and advection module, calculating a concentration of sediment by size at a plurality of depths using an ocean bottom boundary layer hydrodynamic module. A central module uses C++ language to transfer calculations from the sediment suspension and advection module and the ocean bottom boundary layer hydrodynamic module.