Abstract. Entrainment and suspension of sand and gravel is 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, and especially the near-bed concentration. Most entrainment theories 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, as well as powerful turbidity currents. 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 14 dimensionless parameters that may determine entrainment and suspension relations, and applied multivariate regression analysis. A best-fit two-parameter equation (r2 = 0.79) shows that near-bed entrainment, evaluated at 10 % of the flow depth, increases with the ratio of skin-friction shear velocity to settling velocity (u*skin / wsi), as in previous relations, and 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 between turbulent diffusivities of sediment and momentum. The best-fit relation for βi (r2 = 0.40) indicates greater relative sediment diffusivities for rivers with greater flow resistance, possibly due to bed-form induced turbulence, and smaller u*skin / wsi; the latter effect makes the dependence of Rouse number on u*skin / wsi nonlinear, and therefore different from standard Rousean 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 are a significant improvement compared to previous work, extend the calibrated parameter space over a wider range in sediment sizes and flow conditions, and result in 95 % of concentration data predicted within a factor of nine.