An original mathematical formulation for suspended sediments in a two-dimensional wave boundary layer is presented. The model accounts for non-immediate adaptation of sediments to the hydrodinamic conditions, and allows to include the effect of sorting of the different diameters considered. The mathematical model is numerically solved through a finite difference scheme. It is suitable that results compare favourably with experimental data by Staub et alii.