Sand Dune Formation: Physics of Desert Landforms

1. Aeolian Transport: How Wind Moves Sand

Aeolian transport (from Aeolus, Greek god of wind) is the wind-driven movement of sediment. For sand grains — typically 0.1 to 0.5 mm in diameter — the key parameter is the fluid threshold: the wind speed at which aerodynamic drag and lift first overcome gravity and inter-grain friction to set a grain in motion.

The German geologist Ralph Bagnold, studying North African deserts during World War II, established the foundational relationship between wind friction velocity and sand transport. He showed that the threshold friction velocity u*_t scales with grain diameter:

Critically, once grains are already bouncing there is a lower impact threshold (roughly 80% of the fluid threshold). Impacting grains knock stationary ones into motion, so transport can sustain itself in winds weaker than those required to initiate it. This hysteresis means dune fields can persist under conditions that would never have started them from a flat surface.

Bagnold's sand transport equation — that mass flux Q scales with the cube of friction velocity — remains the cornerstone of aeolian sediment budgeting:

2. Saltation and Reptation in Detail

Wind-driven sand moves by three distinct mechanisms, each operating at different scales and contributing differently to dune construction.

Saltation

Saltation (from the Latin saltare, to jump) is the primary transport mode for medium sand. A grain lifted from the surface follows a ballistic arc — accelerated forward by drag, pulled down by gravity — and lands 10 to 100 grain diameters downwind. A typical saltation trajectory has a low liftoff angle (around 10–20 degrees) but a steep impact angle (50–80 degrees). The asymmetry is caused by the grain quickly accelerating to a forward speed close to the wind but gaining only modest vertical velocity.

On impact, a saltating grain does two things: it may bounce and resume its own trajectory with partial energy recovery, and it transfers momentum to the bed, ejecting several smaller grains into short hops. The quantity of grains ejected per impact (the splash function) determines how quickly transport ramps up downwind of a bare patch — the saturation length, typically 1–10 metres for desert sand and a key parameter in dune formation models.

Reptation

Reptation (creeping motion) describes grains nudged by the impact of saltating grains but not given enough energy to become airborne themselves. They roll or slide a few grain diameters downwind. Reptation accounts for roughly 20–25% of total sand mass flux and is the dominant supplier of sand to the avalanching slip face of a dune.

Suspension

Very fine particles (below ~0.07 mm) can be lifted into suspension by turbulent eddies and carried hundreds or thousands of kilometres. Suspension builds loess sheets and hazes rather than dunes; classic dune-forming sand is too coarse to stay aloft for long.

3. Dune Initiation: From Flat Bed to Ripple to Dune

A perfectly flat sand sheet is unstable in a sustained wind. Any small surface irregularity deflects the airflow, creating a sheltered wake on its lee side where transport flux drops, sand accumulates, and the bump grows. This is the seed of a dune. The instability is formalised as a linear stability analysis of the coupled system of airflow and sediment flux.

The key insight from such analyses (pioneered by Werner 1995 and extended by Andreotti, Claudin and others in the 2000s) is that two length scales compete:

• The saturation length L_sat: the distance sand flux needs to adjust to changes in the wind. It is related to the inertia of saltating grains and is typically 1–20 m. Below this scale, flux cannot respond quickly enough to amplify bumps, so very small perturbations decay. This imposes a minimum dune size — dunes cannot be arbitrarily small.
• The hydrodynamic drag length: set by the height at which the wind velocity profile feels the bump on the surface. This controls which wavelengths of bump receive maximum wind speedup on their windward side and therefore grow fastest.

The fastest-growing wavelength from linear theory is roughly 20 times L_sat, which for typical desert sand comes out to tens of metres — consistent with the observed spacing of incipient dunes (embryo dunes). As bumps grow taller, nonlinear effects take over: the lee slope steepens until it exceeds the angle of repose (about 33–34 degrees for dry quartz sand), at which point avalanching creates the characteristic sharp slip face.

Sand ripples, with wavelengths of centimetres, form by a different mechanism (shadow- zone instability driven by the ballistic impact pattern of saltating grains) and are distinct from dunes in their dynamics; they often cover the windward face of a dune.

4. Anatomy of a Barchan Dune

The barchan (from the Kazakh word for a crescent-shaped dune) is the most mathematically tractable and visually striking desert landform. It forms in environments with a unidirectional wind, a limited sand supply, and a hard underlying surface (pavement). Its crescent shape is not arbitrary — it is a stable attractor of the dune dynamics.

Windward (Stoss) Slope

The windward face is a gently inclined ramp, typically 10–15 degrees. Wind accelerates as it climbs this ramp (speed-up effect), eroding sand from the flanks and base and transporting it toward the crest. The rate of sand removal from the stoss side determines how fast the dune migrates.

Crest and Brink

Sand reaches the crest and cascades over the sharp brink — the top edge of the slip face. The brink is where the wind separates from the surface, creating a recirculation zone (separation bubble) on the lee side. Inside the bubble, wind speeds drop to near zero and sand deposited there is sheltered from further transport. This is the deposition zone that builds the slip face.

Slip Face

The slip face is the steep avalanche face on the lee side, maintained at the angle of repose (~33 degrees). Sand grains tumble down in thin avalanche sheets when the local slope exceeds this threshold. The slip face advances the dune forward: as the windward side is eroded and the slip face is built up, the entire dune translates downwind while preserving its shape — a kinematic wave of sand.

Horns

The two tapering horns of a barchan extend downwind, flanking the slip face. They form because sand supply is lower at the edges of the dune — grains that reach the horn tip are quickly transported away by the wind, keeping the horns thin and fast-moving. The horns always point downwind. The asymmetry between the two horns reveals any obliquity in the wind direction.

5. Dune Morphology and Wind Regime

Barchans are just one member of a family of dune types, each associated with a different combination of wind regime, sand supply, and vegetation.

Transverse Dunes

Where sand supply is abundant, individual barchans merge laterally into continuous transverse ridges that run perpendicular to the dominant wind. The slip face spans the entire ridge width. Saudi Arabia's Rub' al Khali and the Namib Sand Sea contain vast transverse dune fields.

Longitudinal (Seif) Dunes

Under bimodal winds blowing from two directions at an acute angle, sand forms long ridges aligned roughly parallel to the resultant wind vector. Seif dunes (from the Arabic word for sword) can extend hundreds of kilometres and are the dominant dune type in the Australian interior and large portions of the Sahara.

Star Dunes

When winds blow from three or more directions without a strongly dominant vector, sand builds up into star (or pyramid) dunes with multiple radiating arms and a central peak. Star dunes are the tallest dune type — some in the Sahara and the Erg Chebbi exceed 300 m — and they are remarkably stationary because net sand flux over any given period is nearly zero.

Parabolic Dunes

Where sparse vegetation partially anchors the sand, the dune shape inverts compared to a barchan: the horns trail upwind, anchored by plant roots, while the central nose advances downwind as a blowout. Parabolic dunes are common in coastal and semi-arid environments where grass colonises the edges.

6. Dune Migration and Speed

A migrating dune is a mass-conserving wave: every grain eroded from the windward slope is deposited on the slip face, moving the dune bodily downwind without loss (in the absence of external sand inputs or outputs). The migration speed follows directly from sand flux conservation.

This inverse relationship between height and speed has a striking consequence for dune interactions. When a faster-moving small barchan catches up with a slower large one, the two can collide and exchange sand. Several outcomes are possible depending on their size ratio: the small dune may pass through the large one by tunneling sand via reptation, they may merge into a single larger dune, or the collision may produce two new dunes of different sizes. These interactions were first catalogued in numerical simulations and later confirmed in time-lapse aerial photography of the Moroccan and Peruvian deserts.

Avalanche Dynamics and Self-Organised Criticality

The slip face is maintained in a state of self-organised criticality: the slope fluctuates around the angle of repose, with intermittent avalanche sheets triggered whenever a sufficient mass of reptating grains accumulates near the brink. The size distribution of these avalanche events follows a power law over several decades, analogous to the Gutenberg-Richter law for earthquakes. This universality reflects the same underlying mathematics: a driven dissipative system maintaining itself at a critical point without external tuning.

7. Key Takeaways