Soulsby (1997) proposed the following equation for the total load sediment transport rate

1) q_{bk}^{*}=\left\{\begin{array} 0.005Uh\left(\frac{U-U_{crk}}{\sqrt{R_{k}gd_{k}}}\right)^{2.4}\left(\frac{d_{k}}{h}\right)^{1.2}\,\,\,\,\,\mathrm{for}\,\,U>U_{crk}\,\,\\ 0\,\,\,\,\,\,\,\text{otherwise}\,\, \end{array}\right.
2) q_{sk}^{*}=\left\{\begin{array} 0.012Uh\left(\frac{U-U_{crk}}{\sqrt{R_{k}gd_{k}}}\right)^{2.4}\left(\frac{d_{k}}{h}\right)d_{*k}^{-0.6}\,\,\,\,\mathrm{for}\,\,U>U_{crk}\\ 0\,\,\,\,\,\,\text{otherwise}\,\, \end{array}\right.

where

{q_{bk}}^* = fractional bed-load sediment transport potential [L2/T]

{q_{sk}}^* = fractional suspended-load sediment transport potential [L2/T]

R_{k} = ρsk/ρw − 1 = submerged specific gravity of a particle [-]

ρ_{sk} = sediment density [M/L3]

ρ_{w} = water density [M/L3]

U = effective depth-averaged current velocity [m/s]

U_{crk} = critical depth-averaged velocity for incipient motion [m/s]

The Soulsby-van Rijn formula was developed by calibrating the above equations to the van Rijn (1993) sediment transport model. The formulas were originally proposed for well-sorted sediments. The formulas have been modified here for nonuniform sediments by replacing the median grain size with the grain class diameter and multiplying the critical depth-averaged current velocity with a hiding and exposure correction factor. Here the Wu et al. (2000) hiding and exposure correction factor is utilized but in principle others may also be used. van Rijn (1984a,b; 2007a,b) computed the critical depth-averaged current velocity using the van Rijn formula.