<figure><img src="image_1.png" /></figure>

HAL
open science

# Phase diagrams of dune shape and orientation depending on sand availability

Xin Gao, Clément Narteau, Olivier Rozier, Sylvain Courrech Du Pont

► To cite this version:

Xin Gao, Clément Narteau, Olivier Rozier, Sylvain Courrech Du Pont. Phase diagrams of dune shape and orientation depending on sand availability. Scientific Reports, 2015, 5, pp. 585-597. 10.1038/srep14677. insu-03579557

HAL Id: insu-03579557

https://insu.hal.science/insu-03579557

Submitted on 18 Feb 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Phase diagrams of dune shape and orientation depending on sand availability

Received: 19 May 2015

Accepted: 07 September 2015

Published: 30 September 2015

Xin Gao<fnref n="1" />, Clément Narteau<fnref n="1" />, Olivier Rozier<fnref n="1" /> & Sylvain Courrech du Pont²

New evidence indicates that sand availability does not only control dune type but also the underlying dune growth mechanism and the subsequent dune orientation. Here we numerically investigate the development of bedforms in bidirectional wind regimes for two different conditions of sand availability: an erodible sand bed or a localized sand source on a non-erodible ground. These two conditions of sand availability are associated with two independent dune growth mechanisms and, for both of them, we present the complete phase diagrams of dune shape and orientation. On an erodible sand bed, linear dunes are observed over the entire parameter space. Then, the divergence angle and the transport ratio between the two winds control dune orientation and dynamics. For a localized sand source, different dune morphologies are observed depending on the wind regime. There are systematic transitions in dune shape from barchans to linear dunes extending away from the localized sand source, and vice-versa. These transitions are captured fairly by a new dimensionless parameter, which compares the ability of winds to build the dune topography in the two modes of dune orientation.

Dunes have been primarily classified according to their shape considering both sand availability and wind directional variability<fnref n="1" />. The wind directional variability is usually measured by the RDP/DP dimensionless parameter. It is the ratio between the resultant drift potential (RDP, the norm of the mean sand flux vector) and the drift potential (DP, the mean of the norms of the individual sand flux vectors). In zones of low sand availability under nearly unidirectional winds (i.e., for high, close to one, RDP/DP-values), crescentic barchan dunes propagate on a non-erodible ground. Under the same wind conditions, transverse linear dunes form where sand availability increases. Star dunes are observed in major depositional centers exposed to wide multidirectional wind regimes (i.e., where the RDP/DP-value tends to zero). For moderate wind directional variability, linear dunes develop. Given the seasonal changes in wind direction on Earth, these linear dunes are widespread in terrestrial sand seas². Dunes have also been classified according to their orientation. The crest lines of barchan dunes are approximately perpendicular to the resultant transport direction, whereas, by definition, the interlaced arms of star dunes have no unique orientation³. For linear dunes, Hunter et al.⁴ used the angle $\phi$ between their crestlines and the resultant sand transport direction to describe them as longitudinal ($\phi \le 15°$), oblique ($15° < \phi < 75°$) or transverse ($\phi \ge 75°$).

A significant breakthrough in understanding the relationships between bedform orientation and flow regime was the work of Rubin and Hunter⁵. When considering a bidirectional flow regime, they experimentally showed that sediment bedforms select the orientation for which the sum of the normal to crest components of the two transport vectors reaches its maximum value. Thus, they introduced the gross bedform-normal transport rule (from now the GBNR), which can be generalized to multidirectional

<footnote n="1">¹Equipe de Dynamique des Fluides Géologiques, Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, UMR 7154 CNRS, 1 rue Jussieu, 75238 Paris Cedex 05, France. ²Laboratoire Matière et Systèmes Complexes, Sorbonne Paris Cité, Université Paris Diderot, UMR 7057 CNRS, Bâtiment Condorcet, 10 rue Alice Domon et Léonie Duquet, 75205 Paris Cedex 13, France. Correspondence and requests for materials should be addressed to X.G. (email: gao@ipgp.fr)</footnote>
flows⁶,⁷. Using this rule together with the dune type classification, a given wind regime should therefore be associated with a specific dune shape and a single orientation⁸,⁹. This single line of reasoning has been recently challenged on both experimental and field observations by *Courrech du Pont et al.*¹⁰. They have demonstrated that there are two modes for dune orientation depending on sand availability. Each mode is associated with a specific dune growth mechanism:

1. On an erodible sand bed, dunes grow in height and migrate selecting an orientation that maximizes the normal to crest components of transport, a result consistent with Rubin and Hunter⁵. This is the bed instability mode, which refers to the ability for a flat sand bed to organize in periodic bedforms as soon as there is sediment transport¹¹.

2. On a non-erodible ground with a localized source of sediment, dunes grow by extension away from the source in the direction of the resultant sand flux at the crest. This is the fingering mode. Classical geomorphologic sand sources are depositional areas related to coastal or fluvial systems, pre-existing dunes or topographic obstacles.

Hence, the same wind regime can produce two dune orientations according to the sand availability. In addition, these two modes may locally coexist as a result of changes in sand availability or due to the development of superimposed bedforms.

*Courrech du Pont et al.*¹⁰ have developed a theoretical model to predict dune orientation as a function of a multidirectional wind regime and the prevailing growth mechanism. Model predictions have been confronted against a few experimental data and field examples. Although these comparisons support the validity of the approach, the model was not precisely checked against a wide range of wind conditions. Here we use numerical modeling to describe the diversity of dune patterns that can be produced under bidirectional wind regime using two different conditions of sand availability: an erodible sand bed or a non-erodible ground with a localized sand source. Over the entire parameter space of the bidirectional wind regime, the numerical results are compared to the prediction of the model of *Courrech du Pont et al.*¹⁰ using only one fitting parameter. This parameter is the wind speed-up, a key parameter in the physics of dunes, which introduces the effect of dune topography on the flow and on the subsequent sand flux over dune flanks. Furthermore, we study the phase diagram of linear dunes for the two dune growth mechanisms. Whereas linear dunes systematically occur in the bed instability mode, the dune morphology in the fingering mode exhibits transitions from barchans to linear dunes, and vice versa. We find that these changes in dune morphology cannot be associated to changes in the RDP/DP-value. Instead, these transitions can be captured by a dimensionless parameter that compares the dune height growth rates of the two growth mechanisms. Finally, we show a field example where this parameter map and the dune morphology map correspond.

## Results

In a cellular automaton dune model (see Methods), we set two types of numerical experiments to investigate the effect of sand availability. In a first set of experiments dedicated to conditions of high sand availability, a flat sand bed with a thickness larger than the flow depth ensures that bedforms never reach the non-erodible ground at the bottom of the cellular space. In addition, boundary conditions are periodic. The vertical edges of the cubic lattice are pasted together (see Fig. 1a at $t=0$). In a second set of experiments dedicated to conditions of low sand availability, a circular source of sediment of diameter $20l_0$ is located upstream of a non-erodible ground (see Fig. 1b at $t=0$). The input flux is not fixed, but there is always at least one layer of sedimentary cells on the circular source. Thus, the input flux adapts to the experimental conditions. In these simulations, boundary conditions are open to remove all sedimentary cells that reach the downstream border of the cellular space. For these two conditions of sand availability, we systematically investigate the dune morphology and orientation for periodic bidirectional wind regimes. The wind transport capacity is kept constant but the direction and the duration of the secondary wind are changed from one simulation to another (see Supplementary Note 1 and Table S1). Then, the periodic bidirectional wind regime is fully characterized by the divergence angle $\theta$ and the transport ratio $N$ between the two winds. This transport ratio is simply the ratio between the time spent in the primary and secondary winds over a period of wind reorientation (i.e., $N \ge 1$).

**Sand availability control on dune shape and orientation in bidirectional wind regimes.** Figure 1a shows the development of bedforms starting from a flat sediment layer. Over short times, the bed instability mechanism is responsible for the development of dunes transverse to each individual wind (i.e., two distinct dune orientations may coexist). After a few cycles of wind reorientation, dunes are well established and exhibit a single alignment as they continue to grow in height by pattern coarsening (see Supplementary Note 2 and Movie S1). As shown in Fig. 1b, a completely different linear dune growth scenario occurs when sand availability is restricted to a localized source: a straight finger-like structure develops and elongates with a constant orientation. This finger extends but does not grow in height and width except at the source location. This shape is remarkably stable in space and time despite a permanent flux of sediment parallel to the crest on both sides of the dune (see Supplementary Note 2 and Movie S2). Taken together, Fig. 1a,b reveal the two dune growth mechanisms with obvious differences in orientation according to conditions of sand availability.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Formation and evolution of dune fields for $\theta=130^\circ$, $N=2$ and two conditions of sand availability. (a) Dunes grow in height and wavelength from a flat bed with no restriction in sand availability; (b) A finger dune extends on a non-erodable ground from a localized sand source. The red and blue arrows show the sand flux vectors of the dominant and secondary winds, respectively. The non-erodable ground is shown in gray. The cellular space has a square basis of side $L=600l_0$. As shown by the orientation of the superimposed bedforms in (a) and the orientation of the finger tip in (b), images are taken after the secondary wind at the end of the cycle of wind reorientation.</figcaption></figure>

Figure 2a,b show the dune fields once they have reached a steady-state in morphology and orientation in the parameter space {$\theta, N$} of bidirectional wind regimes under the same two conditions of sand availability as in Fig. 1. When bedforms grow from a flat sand bed (Fig. 2a), a periodic dune alignment is systematically observed. For $\theta < 90^\circ$, the two winds blow always from the same side of the dune and both participate to the formation of linear transverse dunes with gentle upstream slopes (~10°) and slip faces in the lee¹². For $\theta > 90^\circ$, the two winds blow alternatively from both sides of the dune. As a consequence, the slip face switches alternatively from one side to the other. Dunes are more symmetric in shape with slopes of approximately 20°. Their orientations are always more perpendicular than parallel to the dominant wind direction. Hence, the dune aspect ratio seen by the secondary wind is smaller than the dune aspect ratio seen by the primary wind. It explains why superimposed bedforms perpendicular to the secondary wind are observed more easily than those perpendicular to the dominant wind (e.g., Fig. 2a for $\theta=110^\circ$). However, they are transient dune features, which are rapidly blown away by the next wind (see Supplementary Note 2 and Movie S1). All these observed dune orientations are in good agreement with the GBNR.

When bedforms develop from a localized sand source, the steady-state dunes exhibit different shapes. Straight fingers are observed for a range of divergence angles, which depends on the transport ratio (e.g., $65^\circ \lesssim \theta \lesssim 155^\circ$ for $N=1.5$). Outside of this range, finger-like structures extending away from a sand source are unstable and eventually break up into asymmetric barchans as already observed for {$\theta=25^\circ, N=1$} by Reffet et al.¹³ in underwater experiments. Trains of asymmetric barchans are always observed for $\theta < 60^\circ$ or transport ratio $N \ge 3.5$ (Fig. 2b; Supplementary Note 2 and Movie S3). For a given $N$-value smaller than 3.5, two transitions are observed when $\theta$ increases. For small $\theta$-values, trains of barchans with small horns are observed (e.g., $\theta < 65^\circ$ for $N=2$). For intermediate $\theta$-values, finger dunes are observed (e.g., $70^\circ \lesssim \theta \lesssim 150^\circ$ for $N=2$). For large $\theta$-values (e.g., $\theta > 150^\circ$ for $N=2$), asymmetric barchans are observed (see also Taniguchi et al.¹⁴ and Parteli et al.¹⁵). Unlike finger dunes, which are straight and seem to extend indefinitely, these asymmetric barchans have a curved crest line and arms with a limited extension. Their shape can be understood as the result of the interaction of the two growth mechanisms. The barchanoid base is transverse and stretched by two asymmetric arms. The arm that is upwind of the secondary flow has a smaller extension than the other one, whose tip tends to align in the mean transport direction (e.g., $\theta=155^\circ, N=2$ in Fig. 2b). Only one arm can elongate due to its position relative to the barchanoid base, which is a mobile source of sediment. In the numerical simulations with a localized sand source, the transition from a finger dune to asymmetric barchans occurs for a $\theta$-value that depends on the transport ratio $N$, but it always larger than 90°. This is consistent with the flume experiments of Taniguchi et al.¹⁴ and the numerical simulations of Parteli et al.¹⁵ who studied the phase diagram of asymmetric barchans. The migration of asymmetric barchans modifies the orientation of the curvy elongated arm. However, in the entire parameter space {$\theta, N$}, a clear orientation is given by either the finger trend, the alignment of the train of barchans or the elongated horn of asymmetric barchans
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Steady-state dune field morphologies in the parameter space {θ, N} of bidirectional wind regimes for two conditions of sand availability. (a) Dunes grow in height and wavelength from a flat bed with no restriction in sand availability; (b) Dunes extend or propagate on a non-erodable ground from a localized sand source. The red and blue arrows show the sand flux vectors of the dominant and secondary winds, respectively. In all cases, the dominant wind blows from left to right. The non-erodable ground is shown in gray. The cellular space has a square basis of side L = 600 l₀. As shown by the orientation of the superimposed bedforms in (a) and the orientation of the finger tip in (b), images are taken after the secondary wind at the end of the cycle of wind reorientation.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Dune field orientation with respect to the angle $\theta$ between the two winds for different transport ratio *N* and two conditions of sand availability. red *N* = 1, blue *N* = 1.5, green *N* = 2, black *N* = 5. (a) Dunes grow in height and wavelength from a flat bed with no restriction in sand availability; (b) Dunes extend or propagate on a non-erodable ground from a localized sand source. All orientations are defined with respect to the dominant wind direction using 2D spatial autocorrelation. Symbols differentiate finger dunes (circles), trains of barchan dunes (triangles) and coexistence of barchans and finger dunes (squares). For linear dunes, symbols also differentiate between transverse (full circles), oblique (double circles) and longitudinal dunes (open circles). Solid lines show the predicted dune orientation derived from equations (1) and (2) using $\gamma$ = 1.6, the best-fit flux-up ratio. Error bars show the standard deviation for 10 realizations with different diffusion rates to evaluate the role of defects (see Supplementary Note 4). Note the good agreement between the prediction of the numerical model and the analytical solutions.</figcaption></figure>

(see Supplementary Note 3). Contrary to dunes that develop from an erodable ground (Fig. 2a), this orientation is always more parallel than perpendicular to the dominant wind direction.

Figure 2 as a whole shows the diversity of dune patterns and orientations that can be produced under two different conditions of sand availability. In the bed instability mode, note the abrupt change in dune orientation when the $\theta$-value crosses 90°. Such a change does not occur for linear dunes in the fingering mode because finger dunes are longitudinal bedforms built by two winds blowing on either side of their crest whatever the divergence angle $\theta$. Then, finger dunes can be sustained for $\theta$-values smaller than 90° (see Supplementary Note 2 and Movie S4).

**Dune orientation.** Using the numerical model output, we now turn to quantitative estimates of dune trends and sediment transport properties in the parameter space {$\theta, N$} of bidirectional wind regimes. Symbols in Fig. 3a show the dune orientation $\alpha_I$ in the bed instability mode when dunes grow in height
from the underlying sediment layer. For the same {$\theta$, $N$} -values, symbols in Fig. 3b show the dune trend $\alpha_F$ in the fingering mode when dunes grow by extension away from a localized sand source. These two dune trends $\alpha_I$ and $\alpha_F$ measured in the numerical experiments compare well to the solutions of the theoretical model proposed by Courrech du Pont et al.¹⁰.

Let us consider a linear dune of height $H$, width $W$ and orientation $\alpha$ with respect to $\vec{i}$ in the standard basis ($\vec{i}$, $\vec{j}$). On the upstream slope of dunes, the convergence of flow lines results in an acceleration of flow speed and an increase of the bed shear stress¹⁶. Following Courrech du Pont et al.¹⁰ (see Supplementary Note 5), the sand flux at the crest is

$$ \vec{Q}_s(t) = Q_0 (1 + \gamma |\sin(\psi(t) - \alpha)|) (\cos(\psi(t))\vec{i} + \sin(\psi(t))\vec{j}), \quad (1) $$

where $Q_0$ is the saturated sand flux over a flat sand bed, $\psi$ the wind direction and $\gamma = \beta H/W$ the flux-up ratio. $\beta$ is a dimensionless coefficient, which accounts for all the physical ingredients, e.g., roughness, that may affect the speed-up in addition to the apparent dune aspect ratio¹⁷⁻¹⁹. In our numerical simulations of dunes exposed to bidirectional wind regimes, only the time duration changes between the two winds, which ensures that $Q_0$ is constant. For dunes that are large enough to integrate the entire cycle $\Delta T$ of wind reorientation, the dune height growth rate $\sigma$ is assumed to scale like the characteristic flux divergence averaged over a cycle:

$$ \sigma \propto \frac{1}{\Delta T} \int_{\Delta T} \frac{\|\vec{Q}_s(t)\| \sin(\psi(t) - \alpha)}{HW} dt. \quad (2) $$

As in the GBNR⁵, the selected orientation $\alpha_I$ for the bed instability mode is the one with the largest growth rate $\sigma_I$, such that $d\sigma/d\alpha = 0$. For the fingering mode, the selected orientation $\alpha_F$ is the direction of the resultant sand flux at the crest of the dune, i.e., the orientation for which the normal to crest components of transport cancel each other. Obviously, these two orientations depend on the $\gamma$-value. In bidirectional wind regimes, Courrech du Pont et al.¹⁰ determined the ranges of dune trend that may be observed from $\gamma = 0$ to $\gamma \to +\infty$.

Inversely, given some specific dune orientations, we can estimate the $\gamma$-value that fits the best our data. We obtain $\gamma = 1.6$ using a least-squares method fitting simultaneously all the dune alignments measured in the numerical simulations. As shown in Fig. 3, both numerical and analytical results exhibit the same behavior across the entire parameter space {$\theta$, $N$}. Such a good agreement supports that the two proposed dune growth mechanisms are at work in the numerical model, and that equations (1) and (2) capture the essential features of sediment transport at the crest of the simulated dunes. A flux-up ratio of 1.6 is consistent with the average values documented in experiments and aeolian dune fields¹⁹⁻²⁵. This flux-up ratio corresponds to a fractional speed-up ratio smaller than 0.6 when considering a reasonable range of wind speed at a height of 10 m (e.g., $u_{10} < 25 \text{ m s}^{-1}$) and classical aeolian transport laws with a threshold wind velocity of $4.4 \text{ m s}^{-1}$ (see Supplementary Notes 6 and 7).

**Phase diagram of finger linear dunes.** We now address the transition from finger dunes to trains of barchans in the numerical experiments (Fig. 2b). For different transport ratios $N$, the black lines in Fig. 4 show the ranges of $\theta$-values over which finger dunes are observed. This range decreases with an increasing transport ratio. Close to the transitions, it is common to observe fingers with a finite extension that emit barchans from their tip (e.g., for {60°, 1.5} or {155°, 2} in Fig. 2b). This shows that the two dune growth mechanisms coexist; superimposed dunes in the bed instability mode can develop on the top of dunes in the fingering mode to ultimately break them up into sets of barchans. Note that without a localized sand supply but starting from an initial sand pile, the source of sediment that potentially feeds the extension of a finger dune becomes mobile. In this case, a wider variety of dune shapes would be observed across the entire parameter space {$\theta$, $N$}. This variety includes barchans, asymmetric barchans with an elongated arm, tear drop and chestnut-like dunes¹³⁻¹⁵.

Although a finger dune is longitudinal and aligned with the mean sand flux direction, its slopes are built and sustained by the two winds blowing from either side. Equivalently to the height growth rate $\sigma_I$, one can derive the height growth rate $\sigma_F$, which is the typical rate to shape a cross-section of a finger dune. In practice, the $\sigma_F$-value is computed by injecting $\alpha = \alpha_F$ in equation (2). Then, when the dune develops from a localized sand source, the relative contribution of the two dune growth mechanisms to the final dune shape can be quantified by the ratio $\sigma_F/\sigma_I$ between the height growth rates of the fingering and the bed instability modes. By definition $\alpha_I$ is the orientation for which the growth rate is maximum so that $0 < \sigma_F/\sigma_I < 1$.

Within the parameter space {$\theta$, $N$}, Fig. 4a shows that transitions in dune shape are captured fairly well by the ratio $\sigma_F/\sigma_I$. Finger dunes are observed when $\sigma_F/\sigma_I \gtrsim 0.6$. For a given $N$-value and an increasing $\theta$-value, the transition from barchan to finger dunes occurs for a $\sigma_F/\sigma_I$-value which slightly increases with $N$, from $0.46 \pm 0.03$ when $N=1$ to $0.57 \pm 0.02$ when $N=3$ (see black dots in Fig. 5). For the same transport ratio $N$ but larger angles of divergence $\theta$, the transition from a finger dune to asymmetric barchans occurs for a $\sigma_F/\sigma_I$-value which slightly decreases with $N$, from $0.73 \pm 0.05$ when $N=1.5$ to $0.6 \pm 0.01$ when $N=3$ (see red dots in Fig. 5). For $N \ge 3.5$, no stable finger dunes are observed. The
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

**Figure 4.** Transitions in dune shape in the parameter space {θ, N} when the dune develops from a localized sand source. (a) Growth rate ratio σ<sub>F</sub>/σ<sub>I</sub> and (b) RDP/DP-value in the parameter space {θ, N} of bidirectional wind regimes. The growth rates {σ<sub>P</sub>, σ<sub>F</sub>} are computed with equation (2) using the corresponding dune orientations {α<sub>P</sub>, α<sub>F</sub>}, γ = 1.6 and the same {H, W}-values (i.e., the same dune shape). The σ<sub>I</sub>/σ<sub>F</sub>-value depends only on the wind regime. For numerical simulations with low sand availability and localized sand source (Fig. 2b), black and white lines show the zones of the parameter space {θ, N} in which fingers and barchans are observed, respectively. Gray lines indicate the zones for transition in dune shape.

<figure><img src="image_9.png" /></figure>

**Figure 5.** The σ<sub>F</sub>/σ<sub>I</sub>-values for transitions in dune shape as a function of the transport ratio N. For a given N-value and an increasing θ-value, dune morphology changes from barchans to finger dune (black symbols) and from finger dune to asymmetric barchans (red symbols). For a given N-value, error bars give the range of σ<sub>F</sub>/σ<sub>I</sub>-values over which the transition is observed.
$\sigma_F/\sigma_I$-value is then always smaller than 0.53. The parameter $\sigma_F/\sigma_I$ is also consistent with the finger domain extending up to 180° when $N=1$. Note also that for a given $N$-value close to 1 and $\gamma=1.6$, the $\sigma_F/\sigma_I$-value is maximum for divergence angles $\theta$ which are significantly larger than 90°. For larger $N$-values, this $\theta$-value for which the $\sigma_F/\sigma_I$-value is maximum asymptotically tends to 90°. This is a consequence of the speed-up, the $\sigma_F/\sigma_I$-value is maximum for $\theta=90°$ when $\gamma=0$.

Despite the transition does not correspond to a sharp constant $\sigma_F/\sigma_I$-value ($\sigma_F/\sigma_I \approx 0.60 \pm 0.15$, Fig. 5), this simple parameter quite separates the finger dune and barchan domains. The ratio $||\vec{Q}_F||/||\vec{Q}_I||$ between the norms of sand fluxes associated to each mode, which also makes sense, captures as well the transition. Yet it is highly dependent of the $\gamma$-value (see Supplementary Note 5). It is worth mentioning that $\sigma_F/\sigma_I$ does not consider that the two modes may interact and is computed with a constant flux-up ratio ($\gamma=1.6$) with no assumption on dune shape. It depends on the sole wind regime as the widely used RDP/DP parameter.

As the $\sigma_F/\sigma_I$-value, the RDP/DP-value depends on the divergence angle $\theta$ and the transport ratio $N$ between the two winds. But this classic measurement of the wind variability fails to capture all the transitions in dune shape (Fig. 4b). Although the transitions from barchan to finger dunes are observed for $0.7 < \text{RDP/DP} < 0.9$, the transitions from finger to asymmetric barchans correspond to a wide range of RDP/DP-values, from 0.6 for $N=3$ to 0 for $N \to 1$. Basically, the RDP/DP cannot distinguish dune shape in complex wind regime if different sand flux components cancel each other.

Hence, we infer that, in zone of low sand availability, $\sigma_F/\sigma_I$ can be used to predict the dune field morphology. For large $\sigma_F/\sigma_I$-values, the fingering mode should prevail and linear dunes aligned parallel to the resultant sand flux should be observed. For small $\sigma_F/\sigma_I$-values, finger dunes are unstable, the bed instability should prevail and migrating barchans should be observed.

## Discussion

To test these predictions against natural conditions in modern terrestrial sand seas, we use field images supplied by Google Earth and wind data supplied by the ERA-Interim dataset$^{2,27}$. This global atmospheric reanalysis provides updated wind speed and orientation at 10 m above ground since 1979, with a 0.25° horizontal resolution and a 6 h time resolution. Undisturbed saturated sand fluxes $Q_0$ are computed with the transport law of *Ungar and Haff*$^{28}$, which scales as the wind velocity square minus the transport onset velocity square. The transport onset friction velocity is chosen to be $19 \text{ cm s}^{-1}$, consistently with the measurements of *Iversen and Rasmussen*$^{29}$. {$\alpha_I, \alpha_F$} and growth rates {$\sigma_I, \sigma_F$} are then computed using $\gamma=1.6$ (see Supplementary Note 6).

Here, we concentrate on the ergs around the Tibesti Massif in east central Sahara, where *Wilson*$^{30}$ and *Mainguet and Callot*$^{31}$ have identified major changes in flow orientations and a variety of dune shapes. In this area of low sand availability where the inter-dune area is free of sand or composed of coarse-grain sediment, we identify and locate the different dune types from recent satellite images. Then, we compare the spatial distribution of barchans and linear finger dunes to the $\sigma_F/\sigma_I$-value map derived from the wind data (Fig. 6a). Barchans occur in zones where $\sigma_F/\sigma_I < 0.2$ (Fig. 6b), whereas linear finger dunes are always observed where $\sigma_F/\sigma_I > 0.4$ (Fig. 6d). In most cases, zones of intermediate values of $\sigma_F/\sigma_I$ exhibit both barchans and finger dunes (Fig. 6c). Although the transition is not observed for the same $\sigma_F/\sigma_I$-value in the simulations and in the field, these observations are consistent with our numerical results and predictions. The quantitative discrepancy between the transition values of $\sigma_F/\sigma_I$ could be ascribed to non fully accurate ERA-Interim wind data in this mountainous area, which is distant of any meteorological station.

In the ergs around the Tibesti Massif where linear dunes occur, dunes always roughly align with the resultant sand flux. It indicates that these linear dunes recognized as seifs by *Bagnold*$^{32}$ or silks by *Mainguet and Callot*$^{31}$ are finger dunes that grow by extension from localized depositional areas. Note that zones of transition in this field coincide to places where the wind tends to be more unidirectional. It corresponds to the transition from barchans to finger dunes when the $\theta$-value increases for a given $N$-value, which is also captured by the RDP/DP parameter (Fig. 4b). Nevertheless, barchan dunes are not restricted to large RDP/DP values. In agreement with our predictions, barchans can occur in zones of moderate RDP/DP-value, providing that $\sigma_F/\sigma_I$ is small enough. This is for example the case in sand corridors between the Badain Jaran and the Tengger deserts where two wind components cancel each other (see Fig. 7).

Our numerical experiments show that sand availability acts on both dune shape and orientation. For a given condition of sand availability, the dune shape and orientation are uniquely determined by the wind regime. For periodic bidirectional wind regimes and two different conditions of sand availability, we document all the possible dune morphologies and orientations in the numerical model. In quantitative agreement with the model of *Courrech du Pont et al.*$^{10}$, these orientations can be derived from the sand flux at the crest and two independent dune growth mechanisms. When dunes develop from an erodible sand bed, linear dunes are observed. Their orientation maximizes the dune height growth rate in agreement with the GBNR$^5$. Depending on the wind regime, their orientation can be transverse, oblique or longitudinal to the direction of the resultant transport direction. When dunes develop from a localized source of sediment, two major dune types are observed: finger linear dunes and barchans. Fingers grow by extension as a result of the sediment flux parallel to the dune crest while the dune topography is built and sustained by the sediment flux transverse to the dune crest. The orientation of finger linear dunes
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Transitions in dune shape around the Tibesti Massif (east central Sahara), a zone of low sand availability. (a) Map of σ<sub>F</sub>/σ<sub>I</sub>-value derived from the wind data. The continuous map is generated by linear interpolation based on the horizontal spatial resolution of 0.25° × 0.25° of the wind data. Local dune fields exhibit (b) barchans (image credit: Google Earth), (c) fingers breaking up into barchans (image credit: Google Earth) and (d) fingers elongating in the direction of the resultant sand flux (image credit: Google Earth). All these different dune types are reported in (a) using black symbols. Insets in (b–d) show the local flux roses. The map used in figure (a) is generated by GMT (The general mapping tools).</figcaption></figure>

is the one of the sand flux resultant at the dune crest. For strongly asymmetric wind regime ($\theta \to 0^\circ$ or $N > 5$) or if transport components cancel each other ($\theta \to 180^\circ$), fingers break up into trains of barchans because superimposed dunes in the bed instability mode develop and propagate parallel to the main dune crest. In zones of low sand availability, the dune morphology (i.e., barchan or finger dune) can be inferred from the growth rate ratio $\sigma_F/\sigma_I$, which captures fairly well the transitions in dune types both in the simulations and in a field example. Although the chosen field example supports the validity of the approach for a multidirectional wind regime, the robustness of the parameter $\sigma_F/\sigma_I$ to predict the dune morphology should be tested in complex wind regimes. The transition in dune shape would also deserve to be investigated further by exploring more precisely the dynamics of superimposed bedforms and the role of cohesion<sup>33,34</sup>. In addition, it would be interesting to systematically test the stability and the behavior of dune fields to a spatial or a temporal change of wind regime or sand supply.
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Barchan dunes in a zone of low RDP/DP and $\sigma_F/\sigma_I$-values. These dunes are located in a sand corridor between the Badain Jaran and the Tengger deserts (China, 39°23′25″ N, 102°59′26″ E, image credit: Google Earth). Inset shows the local sand flux roses. In this area, the mean RDP/DP-value and the flux roses derived from the ERA-Interim data are in agreement with the direct measurements of Zhang *et al.*³⁸, ten kilometers away from the dunes. The map used in this figure is generated by GMT (The general mapping tools).</figcaption></figure>

## Methods

The numerical model combines a cellular automaton of sediment transport with a lattice-gas cellular automaton for high Reynolds-flow simulations¹¹,³⁵,³⁶. This model is a fully coupled system in which there is a complete feedback mechanism between the flow and the evolving bed elevation profile.

In the cellular automaton of sediment transport, the lattice is made of cubic cells of elementary length scale $l_0$. A cell is defined by its state: fluid, mobile and immobile sediment. Erosion, transport, deposition, diffusion and avalanching are simulated using different sets of transitions between pairs of nearest neighbor cells. Meanwhile, the particles of the lattice-gas cellular automaton can propagate and collide within the fluid state of the model of sediment transport. Local wind speed and direction are computed by averaging the velocity vectors of the fluid particles in space and time. There is a permanent effect of topography on the flow because the fluid particles rebound back along their initial trajectories when they impact on sedimentary cells. Reciprocally, the flow continuously affects the topography through sediment transport in zones where the basal shear stress is strong enough to mobilize sedimentary particles. In order to determine the effect of flow velocity on the erosion rate $\Lambda_e$, we locally compute the spatial variation of the tangential velocity profile. In the referential of the topography, it writes

$$ \tau_s = \tau_0 l_0 \frac{\partial V_x}{\partial z}, \qquad (3) $$

where $\tau_0$ is the stress scale of the model and $V_x$ is dimensionless and expressed as a number of fluid particles. Then, we consider that the erosion rate $\Lambda_e$ is linearly related to the bed shear stress $\tau_s$ beyond a threshold value $\tau_1$
$$ \Lambda_e = \begin{cases} 0 & \text{for } \tau_s \le \tau_1, \\ \Lambda_0 \frac{\tau_s - \tau_1}{\tau_2 - \tau_1} & \text{for } \tau_1 \le \tau_s \le \tau_2, \\ \Lambda_0 & \text{else,} \end{cases} \quad (4) $$

where $\Lambda_0$ is a constant rate, $\tau_1$ is the threshold for motion inception and $\tau_2$ is a parameter to adjust the slope of the linear relationship. By definition, $(\tau_s - \tau_1)$ is the excess shear stress from which we can account for the feedback mechanism of the bed shear stress on topography. This erosion rate and the rates associated to each set of transition translate into the model the characteristic time scale of the underlying physical processes. The characteristic length and time scales {$l_0, t_0$} of the model are entirely defined with respect to the most unstable wavelength for the formation of dunes and the saturated sediment flux¹¹. Therefore, the model can be used to quantitatively address dune morphodynamics in various physical environments (see Table 3 in Zhang et al.³⁷).

The lattice gas cellular automaton and the model of sediment transport have two independent time scales. As in nature, the evolution of dune topography in the numerical simulations is many orders of magnitude slower than the evolution of the surrounding turbulent flow. Practically, the temporal coupling is controlled by the dynamic topography and we set the frequency $f_c$ at which the propagation cycles of the flow model occur within the cellular space of the model of sediment transport. Thus, for high $f_c$-values, expressed in units of frequency $1/t_0$, the flow always adapts to the current state of the topography. When the overall wind conditions are changing (e.g., after rotation), a stabilization of the flow is systematically performed before starting the processes of sediment transport.

A rotation of the cellular space simulates a change in wind orientation. After each rotation, the model of sediment transport is restarted once the flow has reached a steady-state. We analyze periodic bidirectional wind regimes with two winds {$W_N, W_1$} of equal strength but of different durations {$\Delta T_N, \Delta T_1$}. The saturated flux over a flat sand bed has a constant value $Q_0$ over the entire period of wind reorientation $\Delta T = \Delta T_N + \Delta T_1$. The bidirectional wind regime is defined by the divergence angle $\theta \in [0; 180^\circ]$ and the transport ratio $N = \Delta T_N/\Delta T_1$ between the two winds. $W_N$ is the dominant wind so that $N \ge 1$. More details are given in Supplementary Notes.

## References

1. Wasson, R. & Hyde, R. Factors determining desert dune types. *Nature* **304**, 337–339 (1983).

2. Pye, K. & Tsoar, H. *Aeolian Sand and Sand Dunes* (Unwin Hyman, London, 1990).

3. Zhang, D., Narteau, C., Rozier, O. & Courrech du Pont, S. Morphology and dynamics of star dunes from numerical modelling. *Nat. Geosci.* **5**, 463–467 (2012).

4. Hunter, R. E., Richmond, B. M. & Alpha, T. R. Storm-controlled oblique dunes of the Oregon coast. *Geol. Soc. Am. Bull.* **94**, 1450–1465 (1983).

5. Rubin, D. M. & Hunter R. E., Bedform alignment in directionally varying flows. *Science* **237**, 276–278 (1987).

6. Werner, B. & Kocurek, G. Bed-form dynamics: Does the tail wag the dog? *Geology* **25**, 771–774 (1997).

7. Ping, L., Narteau, C., Dong, Z., Zhang, Z. & Courrech du Pont, S. Emergence of oblique dunes in a landscape-scale experiment. *Nat. Geosci.* **7**, 99–103 (2014).

8. Kocurek, G. & Ewing, R. Aeolian dune field self-organization-implications for the formation of simple versus complex dune-field patterns. *Geomorphology* **72**, 94–105 (2005).

9. Fenton, L. K., Michaels, T. I. & Beyer, R. A. Inverse maximum gross bedform-normal transport 1: How to determine a dune-constructing wind regime using only imagery. *Icarus* **230**, 5–14 (2014).

10. Courrech du Pont, S., Narteau, C. & Gao, X. Two modes for dune orientation. *Geology* **42**, 743–746 (2014).

11. Narteau, C., Zhang, D., Rozier, O. & Claudin, P. Setting the length and time scales of a cellular automaton dune model from the analysis of superimposed bed forms. *J. Geophys. Res.* **114**, F03006 (2009).

12. Zhang, D., Narteau, C. & Rozier, O. Morphodynamics of barchan and transverse dunes using a cellular automaton model. *J. Geophys. Res.* **115**, F03041 (2010).

13. Reffet, E., Courrech du Pont, S., Hersen, P. & Douady, S. Formation and stability of transverse and longitudinal sand dunes. *Geology* **38**, 491–494 (2010).

14. Taniguchi, K., Endo, N. & Sekiguchi, H. The effect of periodic changes in wind direction on the deformation and morphology of isolated sand dunes based on flume experiments and field data from the Western Sahara. *Geomorphology* **179**, 286–299 (2012).

15. Parteli, E. J. et al. Origins of barchan dune asymmetry: insights from numerical simulations. *Aeolian Res.* **12**, 121–133 (2014).

16. Jackson, P. & Hunt, J. Turbulent wind flow over a low hill. *Q. J. Roy. Meteorol. Soc.* **101**, 929–955 (1975).

17. Britter, R., Hunt, J. & Richards, K. Air flow over a two-dimensional hill: Studies of velocity speed-up, roughness effects and turbulence. *Q. J. Roy. Meteorol. Soc.* **107**, 91–110 (1981).

18. Fourrière, A., Claudin, P. & Andreotti, B. Bedforms in a turbulent stream: formation of ripples by primary linear instability and of dunes by nonlinear pattern coarsening. *J. Fluid Mech.* **649**, 287 (2010).

19. Courrech du Pont, S. Dune morphodynamics. *Comptes Rendus Physique* **16**, 118–138 (2015).

20. Lancaster, N., Nickling, W., Neuman, C. & Wyatt, V. Sediment flux and airflow on the stoss slope of a barchan dune. *Geomorphology* **17**, 55–62 (1996).

21. Wiggs, G. F., Livingstone, I. & Warren, A. The role of streamline curvature in sand dune dynamics: evidence from field and wind tunnel measurements. *Geomorphology* **17**, 29–46 (1996).

22. Neuman, C. M., Lancaster, N. & Nickling, W. Relations between dune morphology, air flow, and sediment flux on reversing dunes, Silver Peak, Nevada. *Sedimentology* **44**, 1103–1111 (1997).

23. Walker, I. J. & Nickling, W. G. Simulation and measurement of surface shear stress over isolated and closely spaced transverse dunes in a wind tunnel. *Earth Surf. Processes Landforms* **28**, 1111–1124 (2003).

24. Wu, X., Zou, X., Zheng, Z. & Zhang, C. Field measurement and scaled-down wind-tunnel model measurement of airflow field over a barchan dune. *J. Arid. Environ.* **75**, 438–445 (2011).
25. Claudin, P., Wiggs, G. & Andreotti, B. Field evidence for the upwind velocity shift at the crest of low dunes. *Boundary-Layer Meteorol.* **148**, 195–206 (2013).

26. Uppala, S. M. et al. The ERA-40 re-analysis. Q. J. Roy. Meteorol. Soc. **131**, 2961–3012 (2005).

27. Dee, D. et al. The ERA-interim reanalysis: Configuration and performance of the data assimilation system. Q. J. Roy. Meteorol. Soc. **137**, 553–597 (2011).

28. Ungar, J. & Haff, P. Steady state saltation in air. Sedimentology **34**, 289–299 (1987).

29. Iversen, J. D. & Rasmussen, K. R. The effect of wind speed and bed slope on sand transport. Sedimentology **46**, 723–731 (1999).

30. Wilson, I. G. Desert sandflow basins and a model for the development of ergs. Geogr. J. **137**, 180–199 (1971).

31. Mainguet, M. & Callot, Y. *L'Erg de Fachi-Bilma (Tchad-Niger): Contribution à la connaissance de la dynamique des ergs et des dunes des zones arides chaudes* (Éd. du Centre national de la recherche scientifique, 1978).

32. Bagnold, R. A. *The Physics of Blown Sand and Desert Dunes* (Chapman and Hall, London, 1941).

33. Rubin, D. M. & Hesp, P. Multiple origins of linear dunes on Earth and Titan. Nat. Geosci. **2**, 653–658 (2009).

34. Araújo, A. D., Parteli, E. J., Pöschel, T., Andrade, J. S. & Herrmann, H. J. Numerical modeling of the wind flow over a transverse dune. Sci. Rep. **3**, 2858 (2013).

35. Narteau, C., Le Mouël, J., Poirier, J., Sepúlveda, E. & Shnirman, M. On a small-scale roughness of the core-mantle boundary. Earth Planet. Sci. Lett. **191**, 49–60 (2001).

36. Rozier, O. & Narteau, C. A real-space cellular automaton laboratory. Earth Surf. Processes Landforms **39**, 98–109 (2014).

37. Zhang, D., Yang, X., Rozier, O. & Narteau, C. Mean sediment residence time in barchan dunes. J. Geophys. Res. Earth Surf. **119**, 451–463 (2014).

38. Zhang, Z., Dong, Z. & Li, C. Wind regime and sand transport in China's Badain Jaran Desert. Aeolian Res. **17**, 1–13 (2015).

## Acknowledgements

The data in this paper are the result of numerical simulations using the Real-Space Cellular Automaton Laboratory (ReSCAL), a free software under GNU general public license. The source codes can be downloaded from http://ipgp.fr/rescal. The wind data used in Figs. 6 and 7 are available at http://apps.ecmwf.int/datasets/, the website of the European Centre for Medium-Range Weather Forecasts. Images of Figs. 6b–d and 7 are courtesy of Google Earth. We acknowledge financial support from the UnivEarthS LabEx program of Sorbonne Paris Cité (ANR-10-LABX-0023 and ANR-11-IDEX-0005-02) and the French National Research Agency (ANR-09-RISK-004/GESTRANS and ANR-12-BS05-001-03/EXO-DUNES).

## Author Contributions

X.G. carried out all the numerical simulations and all the statistical data analysis. O.R. developed the real-space cellular automaton, a free software under GNU general public license. X.G., C.N. and S.C.P. designed the study and wrote the manuscript. All authors shared ideas and discussed the results.

## Additional Information

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Gao, X. et al. Phase diagrams of dune shape and orientation depending on sand availability. Sci. Rep. **5**, 14677; doi: 10.1038/srep14677 (2015).

<figure><img src="image_12.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/