# COASTAL DUNE DYNAMICS AND PROCESSES

Μηχανική

24 Οκτ 2013 (πριν από 4 χρόνια και 6 μήνες)

88 εμφανίσεις

DESERT DUNE DYNAMICS AND PROCESSES

By LEVENT YILMAZ, (visited Professor at Louisiana State University, Coastal Studies
Institute, Baton Rouge, LA, 70803, USA)

Technical University of Istanbul, Civil Engineering Faculty, Hydraulic Division, 80626,

Masla
k, Istanbul, Turkey, e:mail :
lyilmaz@itu.edu.tr

Abstract

The development of the dunes are governed by the effects of turbulence. Turbulence is a
type of fluid flow that is strongly rotational and apparently cha
otic. Turbulence separates
nearby parcels of air and thus mixed fluid properties . The evolution of sand dunes is
determined by the interactions between the atmosphere, the surface and the transport and
deposition of sand. We are concerned with this physic
al process and its computational
simulation from three perspectives; namely, (1) flow structure; (2) sand transport and
deposition and (3) interactions between flow structure and sand transport
-
deposition,
which determine the dune morphology.

Keywords: Co
astal management, dunes, sediment supply

Introduction

The system of moving bedforms in a flow field can be explained by the sediment
-
continuity equation and the sediment
-
transport equation. It yields

(1)

(2)

where h is the height of the bedforms or topographic height, t is the time,

is the
sediment density, q is the sediment transport capacity in kg.m
-
1
.
s
-
1

,

is t
he shear stress
due to saltation.

Mathematical Model

Equation (2) shows that the sediment transport varies linearly with the shear stress (
Stam, 1994) as

q(x,t)=q
o
+A
1

(x,t)

(3)

where q
o

is a constant basic sediment
-
transport and A
1

is the linearity constant. If
topographic variations are relatively small, it can be generally stated that the shear stress
is formed of a constant basic term (

o
), and a correctio
n term (

1
) that varies with space
and time:

(x,t)=

o
+

1
(x,t)

(4)

If, in a first approximation, this correction is assumed to be linear with the topographic
height:

1
(x,t)=A
2
h(x,t)

(5)

the sediment
-
transport equation becomes:

q(x,t)=q
o
+A
1

o
+A
3
h(x,t)

(6)

where A
3
=A/A
2

Taking the derivative of Equation (6), it yields

(7)

which, substituted in the continuity equation becomes

(8)

This is called the simple wave equation, and it describes the propagation of a wave at
constant velocity. The ratio

(9)

The simple
-
wave equation, in this case, results from
combining a linear shear stress and a
linear sediment
-
transport formula. The ratio in Equation (9) is called the wave velocity. It
concerns the morphodynamics of dunes and this wave velocity is equal to the migration
rate. This wave will advance at a const
ant rate c without changing its shape. The solution
of Equation (8) is given as

h(x,t)=f(x

ct)

(10)

which expresses that at a certain point x at a certain time t the topography will have the
same height as it had at the start (at t=
0) at a point (x

ct). The height of the bedforms
can be any function of ( x

ct ).

Kinematic Wave Approximation

Another formulation of the mathematical model of the topography is given by using
a sediment
-
transport formula instead of a linea
r relationship with the shear stress.
Bagnold’s (1941) model yields

(11)

where C
B

is the Bagnold’s constants with the unit [s
2
m
1/2
kg
1/2
]. If the linear shear stress
are assumed as in the Equations (4) and (5), it is obtained

(12)

and

(13) and (14)

In this mathematical model it can be seen that the migration rate is not a constant but a
more complicated expression that varies with the topographic height, and the wave
velocity is given as

(15)

In this mathematical model the dune shape changes. The migration rate increases with
height, which means that the top of the dune will advance more rapidly than the base. The
lee side of the dune will tend to become steeper and the peak will even
tually overtake the
slipface. In observations, according to this mathematical model, the maximum angle of
repose for sediment will be surpassed and avalanching will occur at the slipface, limiting
this asymmetrical shape. This is a well known type of equat
ion called a breaking
-
wave
equation which shows the breaking
-
wave behaviour at the mathematical solution. This
behaviour of the solution will be similar for any non
-
linear sediment
-
transport equation,
as long as the peak has a higher velocity than the base

which is given in mathematical
formulation as

(16)

The deformation of the bedform in a “dune
-
like shape” is partly due to the non
-
linear
relationship between sediment transport and shear stress.

Figure: 1 The comparison of dimensionless lengths of barchan dunes (Stam, 1994)

Development of an analytical solution for bedform migration and growth

The expression of the shear stress of the velocity in terms of the topography makes
it
adequate for the development of an analytical solution. The analytical solution becomes
simpler if only one wave number is considered, so that the summations in the Fourier
series are reduced to only one term.

Fi
gure:2 Comparison of the Dimensionless lengths of the Dunes (Stam, 1994)

The general line is analogous to the development of the simple and kinematic wave
equations . The continuity equation is given as

(17)

Dimensionl
ess coordinates have been used (indicated by an asterisk *), so that

(18), (19) and (20)

where H is the maximum topographic height [m].

A spatial derivative expression of the sediment transport has to be substitute
d in the
continuity equation. A linearization of Bagnold’s sediment
-
transport formula will be used
(Stam, 1994). For the linearization it has to be considered that the shear stress (

(x
*
,t))
results from Prandtl’s logarithmic profile (

o
) with Jackson and

Hunt’s (1975)
dimensionless first
-
order correction (

1
(x
*
,t))

(21)

Bagnold’s linearized sediment
-
transport formula then becomes

(22)

where:

o

= shear stress from logarithmic profile [Pa
]

1
= first order correction to the shear stress from the logarithmic profile [dimensionless]

= perturbation factor. This is dimensionless number smaller tan 1. (Stam, 1994)

Differentiation of the linearized transport equation to the dimensionless coor
dinate x
*

gives

(23)

Jackson and Hunt’s (1975) expression for the correction of the shear stress can be used
for the developing of this equation. Expressed as a Fourier Transform

1
(k
*
,t) the shear
stress correction is gi
ven as

(24) and (25)

and

=

/4 if k
*
>0 (positive wave number) (25a)

=
-

/4 if k
*
<0 (negative wave number)

where:

l = thickness of th
e inner region [m]

K
o
= Modified Bessel function of the zero
-
order

K
1
= Modified Bessel function of the first

order

k
*
= dimensionless wave number. It should be noted that by introducing the dimensionless
coordinate x
*
=x/L, the maximum wave length (

) has b
ecome equal to one and therefore
the wave number (k=2

/

) has become dimensionless also. For small arguments, the

Figure:3 Comparison of the Dune Dimensions (Yilmaz, 1997)

Conclusions

The evolution of desert dunes is deter
mined by the interactions between the atmosphere,
the surface and the transport and deposition of sand while the morphology and dynamics
of Mediterranean Aeolian sand dunes are governed by sand movement induced by shore
wave shear. In conditions of unidire
ctional constant winds and sand supply, it is well
known that transverse and Mediterranean shore dunes migrate downwind without
changing their shapes in comparison with the desert barchan dunes. Beach ridges or
coastal dunes consist of also compound dunes
made up of two or more dunes of the same
basic type, coalescing or overlapping, and complex dunes in which two or more different
basic types are combined or superimposed.

References

Alfaro, S.C., Gaudichet, A., Gomes, L., and Maille, M. 1997. Modeling th
e size
distribution of a soil aerosol produced by sandblasting. J. Geophys. Res., 102, pp. 11239
-
11249.

Al
-
Hinai, K.G., Moore, J.M. 1987. Monitoring of sand migration in eastern Arabia by
remote sensing . 20
th

International Symposium on Remote Sensing of
Environment,
Proceedings. Environment Research Institute of Michigan: Ann Arbor, Michigan.

Anderson, R.S. 1988. The pattern of grainfall deposition in the lee of Aeolian Dunes.
Sedimentology, 35: pp. 175
-
188.

Anderson, R.S., Hallet, B. 1986. Sediment tra
nsport by wind; toward a general model.
Geological Society of America Bulletin 97; pp. 523
-
535.

Bagnold, R.A. 1941. The Physics of Blown Sand and Desert Dunes . Methuen: London.

CERC (Coastal Engineering Research Center).1998. Shore Protection Manual. U
.S.
Army Corps of Engineer Waterways Experiment Station, Washington, D.C.; U.S.
Government Printing Office.

Coursin, A. 1964. Observations et experiences faites en avril et mai 1956 sur les
barkhanes du Souehel el Abiodh (region est de Port
-
Etienne). Inst
itut Francaise de
l’Afrique Noire, Bulletin, Serie A 26: pp. 989
-
1022.

Fryberger, S.G. 1979.
Dune forms and wind regime. In McKee ED, (ed). A Study of
Global Sand Seas . United States Government Printing Office: Washington.

Fryberger, S.G., Al
-
Sari, A.M.
, Clisham, T.J., Rizvi, S.A.R., Al
-
Hinai, K.G. 1984. Wind
sedimentation in the Jafurah sand sea, Saudi Arabia. Sedimentology 31; pp. 413
-
431.

Gillette, D.A. 1974. On the production of soil wind erosion aerosols having the potential
for long range transpor
t.
Journal de Recherches Atmospheriques, 8, pp. 735
-
744.

Gillette, D. A. and Stockton, P.H. 1989. The effect of nonerodible particles on wind
erosion at erodible surfaces. J. Geophys. Res., 94, pp. 12885
-
12893.

Greeley, R., Iversen, J.D. 1985. Wind as A

Geological Process on Earth, Mars, Venus
and Titan. Cambridge University Press; Cambridge.

Hastenrath, S.L. 1967. The barchans of the Arequipa region, southern Peru.
Zeitschrift
fur Geomorphologie 11; pp. 300
-
311.

Howard, A.D., Morton, J. B., Gad
-
el
-
Hak
, M., Pierce, D. B., 1978. Sand transport model
of barchan dune equilibrium. Sedimentology 25: pp. 307
-
338.

IAHR, 1960. 17
th

Congress of the International Association for Hydraulic Research,
Proceedings: Hydraulic Engineering for improved water management

-
Germany.

Jackson, P.S., Hunt, J.C.R. 1975. Turbulent wind flow over a low hill. Quarterly Journal
of the Royal Meteorological Society, 101; pp. 929
-
955.

Lancaster, N. 1995. Geomorphology of Desert Dunes . Routledge; London.

Long, J.T., S
harp, R.P. 1964. Barchan
-
dune movement in Imperial Valley, California,
Geological Society of America Bulletin 75: pp. 149
-
156.

Lettau, K., Lettau, H.H. 1978.
Experimental and micrometeorological field studies of
dune migration. In Lettau, H.H., Lettau, K
. (eds). Exploring the World’s Driest Climates.
Institute of Environmental Science Report 101. Center for Climatic Research, University
-
147.

Lumley, J.L. 1970. Stochastic Tools in Turbulence, Academic Press, New York.

Luml
ey, J.L., and Panofsky, H.A. 1964. The Structure of Atmospheric Turbulence:
Wiley
-
Interscience, New York.

McKee, E.D. 1979. Introduction to a study of global sand seas . In McKee, E.D., editor,
A Study of Global Sand Seas, number 1052, pp. 1
-
17, Proffessi
onal Paper U.S. Geol.
Surv.

Mulligan, K.R. 1988. Velocity profiles measured on the windward slope of a transverse
dune. Earth Surface Processes and Landforms 13: pp. 573
-
582.

Nickling, W. G., Gillies, J. A. 1989 .
Emission of fine
-
graned particulates fr
om desert
soils. In Leinen, M. and Sarnthein, M., editors, Paleoclimatology and Paleometeorology:
Modern and Past Patterns of Global Atmospheric Transport, pp. 133
-
165, Dordrecht.

Ruelle, D., Takens, F., 1971 . On the nature o
f turbulence. Comm. Math. Phys. 20, pp.
167
-
192, and 23, pp. 343
-
344.

Shao, Y., 2000. Physics and Modelling of Wind Erosion, Kluwer Academic Publishers,
Dordrecht.

Shao, Y., McTainsh, G. H., Leys, J. F., and Raupach, M. R. 1993a. Efficiency of
sediment
samplers for wind erosion measurement.
Aust. J. Soil Res., 31, pp. 519
-
532.

Stam, J. M. T. 1994. Process
-
Based Modelling of Aeolian Bedforms, Ph. Degree Thesis
at the Technical University of Delft, Netherland.

White, B.R., Schulz, J.C. 1977. Magnus effec
t in saltation. Journal of Fluid Mechanics
81; pp. 497
-
512.

Wiggs, G.F.S., Livingstone, I., Warren, A. 1996. The role of streamline curvature in sand
dune dynamics: evidence from field and wind tunnel measurements.
Geomorphology, 17;
pp. 29
-
46.

Yilmaz,
L. 1997. Information on Dunes at the Turkish Mediterranean Caretta
-
Caretta
Coast, edited to the International Project Meeting “ Integration of Groundtruth &
Spaceborne Sensor’s Data with GIS for Coastal Zone Management at the Turkish
Mediterranean Coasts
, February, Istanbul, Turkey.

Zeman, O., Jensen, N.O. 1988.
Progress report on modeling permanent form sand dunes.
Riso National Laboratory Report RISO
-
M
-

2738, Roskilde, Denmark.