We present in this paper a new three-dimensional (3-D) model for bed-load sediment transport, based on a Lagrangian description. We analyze generalized sub-models for the velocities after collision and the representation of the bed-roughness. The free-flight sub-model includes the effect of several forces, such as buoyancy, drag, virtual mass, lift, Basset and Magnus, and also addresses the particle rotation. A recent methodology for saving computational time in the Basset force is also employed. The sub-models for the post-collision velocity and rotation are based on the conservation of linear and angular momentum during the collision with the bed. We develop a new 3-D representation for the bed roughness by using geometric considerations. In order to address the interaction of particles with the turbulent flow, we tracked the particles through a computed turbulent velocity field for a smooth flat plate. This velocity field was used as a surrogate of the 3-D turbulent conditions close to the bed in streams. We first checked that the basic turbulence statistics for this velocity field could be used to approximate those in an open-channel flow. We then analyzed the interaction of the sediment and the turbulence for a single and multiple particles. We compared numerical results with experimental data obtained by Niño and García (1998b). We show that model predictions are in good agreement with existing data, in the sand size range. © 2011 ASCE.