We consider a model for sand surface dynamics equivalent to an evolutionary variational inequality, discuss its possible generalizations and relation to some other continuous models and to the Monge-Kantorovich problem. Numerical approximations are derived for both primal and dual formulations of this variational inequality.