Numerical modelling of aeolian erosion over rough surfaces

The presence of non‐erodible roughness elements on erodible surfaces has the effect of absorbing part of the wind shear stress and thus protecting the erodible surface from wind erosion. This paper examines the shear stress distribution over roughness arrays of varying density, representing the progress of erosion on a bed of erodible and non‐erodible particles. Three‐dimensional numerical simulations, simulating wind flow over a bed of particles covered by roughness elements, were conducted in order to investigate the effect of roughness elements on the shear stress near the surface. The results of these simulations confirm that the erosion of soil by wind is strongly attenuated by the presence of roughness elements on the surface and depends on the geometric properties of the roughness elements. Based on the new numerical results obtained, a refinement of existing theoretical approaches is developed to describe the dependence of the friction velocity upon roughness frontal area and real exposed cover rate. The new formulation proposed will allow a more accurate evaluation of shear stress partitioning as a function of topographic changes. Copyright © 2010 John Wiley & Sons, Ltd.