The simplest “random flight” models for the paths of heavy particles in turbulence have been tested against previous observations of the deposition of glass beads from an elevated source in the atmospheric surface layer. For the bead sizes examined (diameter 50–100 μm), for which the ratio of particle inertial timescale to turbulence timescale τp/ΓL ≪ 1, it was found sufficient to adapt, as others earlier have done, a well-mixed first-order Lagrangian stochastic (“Langevin”) model of fluid element trajectories, simply by superposing a gravitational settling velocity wg and reducing the velocity autocorrelation timescale along the heavy particle trajectory (Γp) relative to the fluid-Lagrangian timescale (ΓL). That is to say, unless details of the particle distribution very close to ground (where τp/ΓL is not small) are of interest, no advantage other than conceptual clarity can be found in the more faithful approach of explicitly modeling particle acceleration by means of the particle equation of motion. With the timescale reduction parameter β ∼ 2, the Langevin model estimated the location and width of the bead deposit swath very well and fixed the peak deposit density to within (at worst) about 100% error (in very stable stratification), but more generally to within about 20%. In the case where trajectories intersected a tall crop canopy, uncertainties in the treatment of deposition proved more significant than nuances of the trajectory algorithm.