Dune mobility and vegetation cover in the Southwest Kalahari desert

As part of a wider project investigating the palaeoenvironmental significance of partially vegetated linear dunes in the southwest Kalahari, data collected in the latter part of 1992 concerning dune movement and vegetation cover suggest that sediment transport is occurring on some dune surfaces, and that the majority of surface activity occurs on the crests and upper slopes of the dunes. The data suggest that the limiting variables on surface sediment movement vary on different parts of a dune. On interdunes and lower dune slopes the primary limiting variable is available wind energy, while on dune crests and upper slopes it is vegetation cover. Ground cover by litter has much greater importance in protecting the surface sediment from erosion than rooted vegetation. From individual data points, no evidence is found to support a threshold vegetation cover below which sediment movement occurs. Rather, a gradient of activity is suggested whereby a reduction in vegetation cover increases the potential for sediment movement and surface change. However, dunes with differing amounts of mean vegetation cover display differing degrees of surface activity, and at this scale, a vegetation cover threshold in the region of 14 per cent may be recognized.     

Estimation of surface roughness (z0) over a stabilizing coastal dune field based on vegetation and topography

Accurate knowledge of the surface roughness and the resultant wind speed are important for many applications, such as climatic models, wind power meteorology, agriculture and erosion hazards, especially on sand dunes in arid and semi‐arid environments, where vegetation cover is scarce. In this study we aimed at quantifying the effects of vegetation cover and topography on surface roughness over a stabilizing dune field on the southern coast of Israel. Forty‐six wind measurements were made at various distances from the coastline, ranging from 10 to 2800 m, and z₀ values were calculated from the wind measurements based on the ratio between the wind gust and the average wind speed. We estimated vegetation cover using the soil adjusted vegetation index (SAVI) from Landsat satellite images for the upwind sector at various lengths, ranging from 15 to 400 m, and based on digital elevation models and differential GPS field measurements we calculated the topographic variable of the relative heights of the stations. z₀ values were positively correlated with the winter SAVI values (r = 0·87 at an upwind length of 200 m) and negatively correlated with the relative height (r = ?0·68 at an upwind length of 200–400 m for the inland dune stations). Using these variables we were able to create a map of estimated z₀ values having an accuracy of over 64%. Such maps provide a better understanding of the spatial variability in both wind speed and sand movement over coastal dune areas. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluating morphological estimates of the aerodynamic roughness of debris covered glacier ice

Aerodynamic roughness length (z₀), the height above the ground surface at which the extrapolated horizontal wind velocity profile drops to zero, is one of the most poorly parameterised elements of the glacier surface energy balance equation. Microtopographic methods for estimating z₀ have become prominent in the literature in recent years, but are rarely validated against independent measures and are yet to be comprehensively analysed for scale or data resolution dependency. Here, we present the results of a field investigation conducted on the debris covered Khumbu Glacier during the post‐monsoon season of 2015. We focus on two sites. The first is characterised by gravels and cobbles supported by a fine sandy matrix. The second comprises cobbles and boulders separated by voids. Vertical profiles of wind speed recorded by a tower comprising five cup anemometers and deployed over both sites enable us to derive measurements of aerodynamic roughness that reflect their observed surface characteristics (0.0184 m and 0.0243 m, respectively). At the second site, z₀ also varied through time following snowfall (0.0055 m) and during its subsequent melt (0.0129 m), showing the importance of fine resolution topography for near‐surface airflow. To compare the wind profile data with microtopographic methods, we conducted structure from motion multi‐view stereo (SfM‐MVS) surveys across each patch and calculated z₀ using three previously published approaches. The fully three‐dimensional cloud‐based approach is shown to be most stable across different scales and these z₀ values are most correct in relative order when compared with the wind tower data. Popular profile‐based methods perform less well providing highly variable values across different scales and when using data of differing resolution. These findings hold relevance for all studies using microtopographic methods to estimate aerodynamic roughness lengths, including those in non‐glacial settings. Copyright © 2017 John Wiley & Sons, Ltd.     

Aeolian dynamics over a coastal foredune,Prince Edward Island,Canada

Near‐surface airflow over a morphologically simple, vegetated, 8 m high foredune with a small wave‐cut scarp was measured for onshore to oblique‐onshore conditions during a low‐moderate (5–6 m s^‐1 ) wind event and a high velocity (11–18 m s^‐1) sand‐transporting gale event. Flow across the foredune was characterized by significant flow compression and acceleration up and across the foredune during both events. During the gale, a pronounced jet (speed bulge) developed at the foredune crest, which increased in magnitude with increasing wind speed. The vertical (W) velocity component of the 3D flow field was positive (upwards) across the stoss slope under low wind conditions but negative (downwards) during gale wind conditions, with upslope acceleration. During the low velocity event, there was speed‐down within the vegetation canopy, as would be expected for a porous roughness cover. During the strong wind event there was speed‐up in the lower portion of the vegetation canopy, and this was found up the entire stoss slope. Sediment transport during the gale force event was substantial across the beach and foredune despite the moderate vegetation cover and minimum fetch. Aeolian suspension was evident in the lee of the dune crest. The observations presented herein show that strong storm winds are an effective mechanism for translating sediment landwards across a high vegetated foredune, contributing sediment to the stoss slope, crest and leeward slopes of the foredune and backing dunes. Copyright © 2013 John Wiley & Sons, Ltd.     

Controls on the aerodynamic roughness length and the grain‐size dependence of aeolian sediment transport

Estimates of the wind shear stress exerted on Earth's surface using the fully rough form of the law‐of‐the‐wall are a function of the aerodynamic roughness length, z₀. Accurate prediction of aeolian sediment transport rates, therefore, often requires accurate estimates of z₀. The value of z₀ is determined by the surface roughness and the saltation intensity, both of which can be highly dynamic. Here we report field measurements of z₀ values derived from velocity profiles measured over an evolving topography (i.e. sand ripples). The topography was measured by terrestrial laser scanning and the saltation intensity was measured using a disdrometer. By measuring the topographic evolution and saltation intensity simultaneously and using available formulae to estimate the topographic contribution to z₀, we isolated the contribution of saltation intensity to z₀ and document that this component dominates over the topographic component for all but the lowest shear velocities. Our measurements indicate that the increase in z₀ during periods of saltation is approximately one to two orders of magnitude greater than the increase attributed to microtopography (i.e. evolving sand ripples). Our results also reveal differences in transport as a function of grain size. Each grain‐size fraction exhibited a different dependence on shear velocity, with the saltation intensity of fine particles (diameters ranging from 0.125 to 0.25 mm) saturating and eventually decreasing at high shear velocities, which we interpret to be the result of a limitation in the supply of fine particles from the bed at high shear velocities due to bed armoring. Our findings improve knowledge of the controls on the aerodynamic roughness length and the grain‐size dependence of aeolian sediment transport. The results should contribute to the development of improved sediment transport and dust emission models. © 2018 John Wiley & Sons, Ltd.     

A field study of mean and turbulent flow characteristics upwind,over and downwind of barchan dunes

Field‐measured patterns of mean velocity and turbulent airflow are reported for isolated barchan dunes. Turbulence was sampled using a high frequency sonic anemometer, deriving near‐surface Reynolds shear and normal stresses. Measurements upwind of and over a crest‐brink separated barchan indicated that shear stress was sustained despite a velocity reduction at the dune toe. The mapped streamline angles and enhanced turbulent intensities suggest the effects of positive streamline curvature are responsible for this maintenance of shear stress. This field evidence supports an existing model for dune morphodynamics based on wind tunnel turbulence measurements. Downwind, the effect of different dune profiles on flow re‐attachment and recovery was apparent. With transverse incident flow, a re‐attachment length between 2·3 and 5·0h (h is dune brink height) existed for a crest‐brink separated dune and 6·5 to 8·6h for a crest‐brink coincident dune. The lee side shear layer produced elevated turbulent stresses immediately downwind of both dunes, and a decrease in turbulence with distance characterized flow recovery. Recovery of mean velocity for the crest‐brink separated dune occurred over a distance 6·5h shorter than the crest‐brink coincident form. As the application of sonic anemometers in aeolian geomorphology is relatively new, there is debate concerning the suitability of processing their data in relation to dune surface and streamline angle. This paper demonstrates the effect on Reynolds stresses of mathematically correcting data to the local streamline over varying dune slope. Where the streamline angle was closely related to the surface (windward slope), time‐averaged shear stress agreed best with previous wind tunnel findings when data were rotated along streamlines. In the close lee, however, the angle of downwardly projected (separated) flow was not aligned with the flat ground surface. Here, shear stress appeared to be underestimated by streamline correction, and corrected shear stress values were less than half of those uncorrected. Copyright © 2011 John Wiley & Sons, Ltd.     

Wind tunnel tests of the dynamic processes that control wind erosion of a sand bed

Aeolian sand transport is a complicated process that is affected by many factors (e.g. wind velocity, sand particle size, surface microtopography). Under different experimental conditions, erosion processes will therefore produce different results. In this study, we conducted a series of wind tunnel experiments across a range of wind velocities capable of entraining sand particles (8.0, 10.0, 12.0, and 14.0 m s^-1) to study the dynamic changes of the shear velocity, aerodynamic roughness length, and sand transport. We found that the shear velocity and aerodynamic roughness length are not constant; rather, they change dynamically over time, and the rules that describe their changes depend on the free-stream air velocity. For wind tunnel experiments without feeding sand into the airflow, the sand bed elevation decreases with increasing erosion time, and this change significantly affected the values of shear velocity and aerodynamic roughness length. A Gaussian distribution function described the relationships between the sand transport rate (q_T) and the duration of wind erosion (T). It is therefore necessary for modelers to consider both deflation of the bed and the time scale used when calculating sand transport or erosion rates. © 2018 John Wiley & Sons, Ltd.     

Dunes on the skeleton coast,Namibia (South West Africa): Geomorphology and grain size relationships

Simple, and locally compound, transverse and barchanoid dunes dominate the 2000 km² Skeleton Coast dunefield in northwestern Namibia/South West Africa. Dune height and spacing are closely correlated (r = 0-89) and decrease across the dunefield from southwest to northeast, with an accompanying change from transverse to barchanoid ridges and ultimately barchans. The dunes are aligned transverse to the dominant strong south and south southwest onshore winds. Alignment patterns indicate that surface roughness changes between coastal plain and dunes cause dune-forming winds to swing to the right over the dunes, but resume their original direction beyond. Grain size and sorting vary at three scales: the dune, the dune landscape and through the dunefield. Overall the sands, derived from three localities by deflation from beaches supplied by vigorous longshore drift, become progressively finer and better sorted across the dunefield paralleling changes in dune height and spacing. A statistically significant relationship (r = ?0?65) was established between dune spacing and the phi grain size of the coarser fraction of the dune sands, demonstrating the importance of the protective effects of coarse grains, and suggesting that the morphometry of simple transverse dunes may be controlled by the scale of turbulence associated with the threshold wind speed required to move the coarsest fraction of the dune sand.     

A test of granulometric control of desert dune geometry

I. G. Wilson's hypothesis of coarse-tail grain-size control of dune spacing was derived from three dunefields in the Sahara and, although it was supplemented by measurements of spacing from other parts of the world, the grain-size data come only from North Africa. In this paper the hypothesis is tested in the Australian dunefields. Australian dunes do not form separate categories on a P₂₀/s (twentieth percentile/spacing) plot and, when placed on Wilson's diagram, the Australian data form a continuum between dunes and draas. Ripples maintain their identity, suggesting that the average saltation length of sand controls ripple wavelength while dunes and Wilson's draas are formed by secondary flow of some kind. In Australia the spread of data on the P₂₀/s diagram indicates that grain-size is not the prime control on s. Data from Australia and the Sahara indicate that direct linear relationships between s and h (dune spacing and height means respectively for blocks of dunes) occur but they have different slopes in different areas. These differences possibly reflect variations in vegetation and substrate as well as differences in wind regime. The separation of dunes from draas in the Sahara reflects the greater role of grain-size in an area where coarse grains frequently occur in dune crests. The relative paucity of coarse grains in Australian dune crests may reflect the fine-grained alluvium from which the dunes are derived. These differences may be ascribed to differing topography in Australia and the Sahara.     

Relations of sand trapping efficiency and migration speed of transverse dunes to wind velocity

A self‐consistent model which describes transverse dune migration in equilibrium is introduced. It shows that an equilibrium expression for dune migration speed (c _d) must take into account sand trapping efficiency (T _E), and that T _E is strongly related to the wind speedup over the windward surface. An expression for sand trapping efficiency (T _E) is analytically derived from a microscale analysis of sand grain deposition on the slip face. Sand trapping efficiency (T _E) is mainly determined by shear velocity on a level surface (u_*(−∞)), and rapidly decreases as u_*(−∞) increases. For each dune height (H), dune migration speed (c _d) first increases, and then decreases monotonically after reaching the maximum, as the shear velocity on a level surface (u_*(−∞)) increases. Dune migration speed (c _d) is not inversely proportional to dune height (H). For low dunes, small sand trapping efficiency (T _E) suppresses c _d, whereas for high dunes, wind speedup and large T _E resist the decrease of c _d. Some field data show the same tendency. The dune‐to‐plane‐bed transition observed in subaqueous and venusian bedforms could be associated with the decrease of sand trapping efficiency (T _E). Copyright © 2000 John Wiley & Sons, Ltd.     

Decadal and seasonal development of embryo dunes on an accreting macrotidal beach: North Lincolnshire,UK

Embryo dunes are often ephemeral, but can develop to become established coastal foredunes. In 2001 a patch of embryo dunes 13.11 m² appeared on a beach in north Lincolnshire, UK and had expanded to over 3600 m² by 2011. The rate of expansion is linked to storm occurrence, where expansion is slowed during years with a higher incidence of storm surges. From July 2009–October 2010 seasonal changes in dune field topography were determined using terrestrial laser scanning (TLS) data. Vegetation is important in the development of embryo dunes, but can cause errors in TLS data. Tests evaluating the impact of vegetation on the TLS data suggest the minimum elevation value from the TLS point cloud within a 0.05 m grid cell gives a good approximation of the ground surface. Digital elevation models (DEMs) of the dunes constructed using filtered data showed the embryo dunes underwent a classic seasonal cycle of erosion during the winter and accretion during the summer. For example from October 2009 to April 2010 over 375 m³ of sediment was eroded from the dunes whereas during spring and summer 2010 the dune field gained over 600 m³ of sand. The overall magnitude of change in dune height and volume from season to season exceeded the errors associated with the construction of the DEM from the TLS data and the vegetation filtering process, which suggests TLS can be useful for documenting topographic change in vegetated dunes. After 10 years, the patch of embryo dunes is still expanding but has not yet merged with more established foredunes to landward. Aeolian process measurements indicate that, at present, the embryo dunes do not prevent sand from reaching the foredunes, however the rate of foredune progradation has slowed concurrently with the expansion of the embryo dune field. Copyright © 2013 John Wiley & Sons, Ltd.     

Sand flux and wind profiles in the saltation layer above a rounded dune top

The near-bed airflow and the movement of sand dune sediments by wind are fundamental dune geomorphological processes.This research measured the wind profiles and sand mass flux on the rounded top of a transverse dune at the southern edge of the Tengger Desert to examine how to best predict the vertical profile of sand flux.This work also tested the accuracy of previously developed models in predicting the apparent roughness length during saltation.Results show that mass flux vertical distribution over the dune top is underestimated by an exponential function,overestimated by a power function,but closely matches the predictions made using the LgstcDoseRsp function.Given suitable values ofα,βandγaccording to the grain size composition,S?rensen equation with the peaked shape of the mass transport curve will well predict the dimensionless mass flux qg/ρu*3against dimensionless shear velocity u*/u*t.The modified Charnock model works best of the previously published models tested,with an R2of 0.783 in predicting the enhanced roughness over the moving sand surface,as opposed to an R2of0.758 for the Owen model and an R2of 0.547 for the Raupach model.For the rounded dune top in this study,C m=0.446±0.016.     

A process‐based hypothesis for the barchan–parabolic transformation and implications for dune activity modelling

The introduction of vegetation to bare barchan dunes can result in a morphological transformation to vegetated parabolic dunes. Models can mimic this planform inversion, but little is known about the specific processes and mechanisms responsible. Here we outline a minimalist, quantitative, and process‐based hypothesis to explain the barchan–parabolic transformation. The process is described in terms of variations in the stabilization of wind‐parallel cross‐sectional dune slices. We hypothesize that stabilization of individual ‘dune slices’ is the predictable result of feedbacks initiated from colonization of vegetation on the slipface, which can only occur when slipface deposition rates are less than the deposition tolerance of vegetation. Under a constant vegetation growth regime the transformation of a barchan dune into a parabolic dune is a geometric response to spanwise gradients in deposition rates. Initial vegetation colonization of barchan horns causes shear between the anchored sides and the advancing centre of the dune, which rotates the planform brinkline angle from concave‐ to convex‐downwind. This reduces slipface deposition rate and allows vegetation to expand inward from the arms to the dune centre. The planform inversion of bare barchans dunes into vegetated parabolic dunes ultimately leads to complete stabilization. Our hypothesis raises several important questions for future study: (i) are parabolic dunes transitional landforms between active and vegetation‐stabilized dune states? (ii) should stabilization modelling of parabolic dune fields be treated differently than linear dunes? and (iii) are stabilized parabolic dune fields ‘armoured’ against re‐activation? Copyright © 2012 John Wiley & Sons, Ltd.     

Locally modified surface winds on linear dunes as derived from directional raingauges

The study deals with local airflow patterns induced by 12–16 m high, partly vegetated dunes, under various prevailing wind conditions. Detailed empirical data were derived by means of an unconventional use of directional raingauges. The utility and accuracy of the method have been previously tested in a combined empirical and numerical study over ridges on a similar scale. Results are presented on the nature, extent and intensity of locally modified patterns of surface flow and on their dependence on wind conditions and slope inclination. Results may serve to verify or refine existing perceptions or models of the flow affecting sand movement over dunes. For example, deflected surface wind flow and even stationary vortices were found to form under winds with resultant angles of incidence of only 5–10° relative to the axes of the dunes, and on slopes inclined only 12–15°, i.e. lower angles than often assumed in previous studies of dune morphology. Various effects of the resulting surface flow on the dune ecosystem are also referred to, e.g. the distribution of settling dust and seed dispersal. Copyright © 2002 John Wiley & Sons, Ltd.     

Early Holocene dune field development in Dalarna,central Sweden: A geomorphological and geophysical case study

Bonäsheden, Sweden's largest continuous dune field, situated in the county of Dalarna, central Sweden, has been investigated using LiDAR (light detection and ranging) remote sensing, ground penetrating radar as well as by field observations and luminescence dating. The use of LiDAR in conjunction with geographic information system (GIS) software proved to be efficient in mapping the inactive dune field and classifying the dune morphology, especially when slope raster images were used. The dunes have formed mostly by winds from the northwest (NW) and are of a transverse type. Still other dune types, such as parabolic dunes, and transverse dunes with a deviating orientation are present. Also, there seems to be different generations of dunes, suggesting a complex palaeowind environment with a change from predominantly north‐westerly winds to more westerly winds. Luminescence dating finally allows us to have an absolute chronology of the development of the Bonäsheden dune field, revealing formation of the dune field closely following the de‐glaciation of this part of Sweden (c. 10.5 ka). The well preserved transverse shape of the majority of the dunes suggests rapid stabilization by vegetation, although sand drift still seems to have been active on a noticeable scale for at least 1500 years and also, occasionally and patchy, as coversand deposition during the Late Holocene. A simple model is proposed for the dune field development of Bonäsheden based on our findings. This model is a useful addition since the majority of present day dune field models focus on the formation of parabolic dunes or large unvegetated dune fields. Our results suggest that most models cannot adequately simulate the formation of such small dune fields as that of Bonäsheden, with apparently rapidly fixated transverse dunes in a previously glaciated, now vegetated area. Copyright © 2017 John Wiley & Sons, Ltd.     

Aeolian sediment flux decay: Non-linear behaviour on developing deflation lag surfaces

Wind tunnel simulations of the effect of non-erodible roughness elements on sediment transport show that the flux ratio q/q_s, shear velocity U*, and roughness density λ are co-dependent variables. Initially, the sediment flux is enhanced by kinetic energy retention in relatively elastic collisions that occur at the roughness element surfaces, but at the same time, the rising surface coverage of the immobile elements reduces the probability of grain ejection. A zone of strong shearing stress develops within 0·03 to 0·04 m of the rough bed because of a relative straightening of velocity profiles which are normally convex with saltation drag. This positive influence on fluid entrainment is opposed by declining shear stress partitioned to the sand bed. Similarly, because the free stream velocity U_f is fixed while U_* increases, velocity at height z and particle momentum gain from the airstream decline, leading eventually to lower numbers of particles ejected on average at each impact. When the ratio of the element basal area to frontal area σ is approximately equal to 3·5, secondary flow effects appear to become significant, so that the dimensionless aerodynamic roughness parameter Z₀/h and shear stress on the exposed sand bed T_s decrease. It is at this point that grain supply to the airstream and saltation drag appear to be significantly reduced, thereby intensifying the reduction in U_*. The zone of strong fluid shear near the bed dissipates.     

Effects of atmospheric conditions on wind profiles and aeolian sand transport with an example from white sands national monument

Quantifying aeolian sand transport rates relies upon the computation of the near-surface shear velocity (u*) determined from velocity profiles of the wind. While it has been recognized that various conditions, such as saltation, surface roughness, surface slope and atmospheric conditions, have an effect on the velocity profile, it is commonly assumed that measurements made above the surface will be representative of the near-surface shear velocity. Airflow and temperature data collected over a flat substrate at White Sands National Monument in New Mexico, however, show the significant effects that atmospheric conditions have on velocity profiles. During the day, when solar insolation is heating the surface, atmospheric conditions in the lowest several metres become unstable, resulting in enhanced convection and vertical mixing so that the velocity gradient changes little with height. As a result, the shear stress in this region of vertical mixing lessens, while the near-surface shear stress is increased because the higher wind speeds are now nearer the surface. At night, the near-surface atmospheric conditions are stable, thereby reducing convection and vertical mixing, resulting in stratified airflow and increased shear velocity away from the surface. Unless this atmospheric effect is accounted for, estimates of sediment transport rates may be in error by as much as a factor of 15 times when wind speeds are near threshold velocity. At wind speeds approaching 10 ms¹, at 5m above the surface, this error in computing sediment transport is reduced to a factor of only two to three times, and may be within the range of measurement error.     

Growth by extension,and reworking,of a south‐western Kalahari linear dune

Optically stimulated luminescence (OSL) dating studies of linear (longitudinal) dunes have been used extensively to elucidate late Quaternary environments and climates in arid or formerly arid regions, yet understanding of the development of such dunes is incomplete. In particular, conflicting opinions have been presented regarding the propensity of linear dunes to migrate laterally, the degree to which they rework their own sediment during accumulation and whether they form primarily by extension, as opposed to lateral sand movement from adjacent interdunes. This study focuses on this last point, although the importance of the other controversies is discussed in context. A simple linear dune in the south‐western Kalahari, which has a prominent termination on a pan (playa) surface, provides an opportunity to directly test hypotheses of dune extension. Chronostratigraphy along a ~600 m transect along the crest of the dune, constrained by 42 OSL ages, reveals that the dune grew by extension on occasions in the late Pleistocene and early Holocene, but has also been subject to reworking along its length, which has continued until recent times. Dune development by extensional growth is suggested to operate under environmental conditions more conducive to net accumulation, whereas reworking is largely independent of conditions throughout the last ~18 ka, and may represent seasonal fluctuations in the position of the dune crest. The relative significance of these two modes of development is suggested to be a key control on the efficacy of linear dunes as archives of environmental and climatic change. Copyright © 2011 John Wiley & Sons, Ltd.     

Runoff generation in a sandy area—the nizzana sands,Western Negev,Israel

This paper presents the results of sprinkling experiments conducted over the lower part of vegetated and crusted linear dunes as well as over flat playa surfaces that appear in the interdune corridors. Data obtained show that these two surface units respond quickly to rainstorms. Runoff generation can be expected for any storm exceeding 2-3 mm and runoff coefficients are high. When the topsoil algal crust, 1 to 2 mm thick, is removed from the surface of the dune, infiltration increases drastically and eliminates any possibility of runoff generation under present-day rainfall conditions. This data may be of great help in the understanding of the geomorphology and sedimentary sequence of the corridors separting linear dunes.