The entrainment and suspension of sand and gravel are important for the evolution of rivers, deltas, coastal areas, and submarine fans. The prediction of a vertical profile of suspended sediment concentration typically consists of assessing (1) the concentration near the bed using an entrainment relation and (2) the upward vertical distribution of sediment in the water column. Considerable uncertainty exists in regard to both of these steps, especially the near-bed concentration. Most entrainment relations have been tested against limited grain-size-specific data, and no relations have been evaluated for gravel suspension, which can be important in bedrock and mountain rivers. To address these issues, we compiled a database with suspended sediment data from natural rivers and flume experiments, taking advantage of the increasing availability of high-resolution grain size measurements. We evaluated 12 dimensionless parameters that may determine entrainment and suspension relations and applied multivariate regression analysis. A best-fit two-parameter equation (r² = 0.79) shows that near-bed entrainment, evaluated at 10 % of the flow depth, decreases with the ratio of settling velocity to skin-friction shear velocity (w_(si)/u_(∗ skin)), as in previous relations, and increases with Froude number (Fr), possibly due to its role in determining bedload-layer concentrations. We used the Rouse equation to predict concentration upward from the reference level and evaluated the coefficient β_i, which accounts for differences in the turbulent diffusivity of sediment from the parabolic eddy viscosity model used in the Rouse derivation. The best-fit relation for β_i (r² = 0.40) indicates greater relative sediment diffusivities for rivers with greater flow resistance, possibly due to bedform-induced turbulence, and larger w_(si)/u_(∗ skin); the latter dependence is nonlinear and therefore different from standard Rouse theory. In addition, we used empirical relations for gravel saltation to show that our relation for near-bed concentration also provides good predictions for coarse-grained sediment. The new relations extend the calibrated parameter space over a wider range in sediment sizes and flow conditions compared to previous work and result in 95 % of concentration data throughout the water column predicted within a factor of 9
