1. Introduction
The sand dune system is important, unique and involves complex physical (e.g., [
1,
2,
3]), ecological (e.g., [
4]), and geomorphological (e.g., [
5]) processes. Vast desert areas (∼20%) comprise dunes [
5,
6], and they can degrade human habitations, agricultural fields, roads, etc. [
7,
8]. However, stable (fixed) dunes may have economic value; for example, the stable dunes in the Sahel and southern Africa have pastoral and agricultural uses. Sand dunes can be classified into two main categories: active (such as barchans, transverse and seif dunes) and fixed (mainly parabolic and vegetated linear dunes), in which the dune’s mobility state depends on several factors including vegetation and biogenic crust cover, wind power, and precipitation. Sand dunes may be affected by climate changes (e.g., [
9,
10,
11]). The wind is the driving force of a sand dune; yet, precipitation affects the vegetation and biogenic crust covers, which can stabilize the sand, thus playing an important role in dune mobility. Previous studies mainly associated dune formation with arid conditions and dune fixation with humid conditions (e.g., [
12,
13,
14,
15,
16,
17]). Droughts can affect dune mobility by suppressing the growth of vegetation that stabilizes the dunes and can lead to dune reactivation as, for example, occurred in North America, in the Mojave Desert [
18], in the Great Sand Dunes [
19], and in the Canadian Prairies [
20].
The relationship between vegetation and the earth’s climate has been intensively studied in the past. A well-known conceptual model regarding the effect of vegetation on the earth’s climate is the “Daisyworld” model of [
21]; in this model, black and white daisies regulate the earth’s temperature to be roughly constant, in spite of the drastic changes in solar insolation. In 1974, Otterman and his colleagues [
13,
22] analyzed satellite images of the Israel-Egypt border. These studies reported a clear difference between the Israeli and the Egyptian sides in which the former was much darker. This difference was attributed to vegetation and soil characteristics—the Egyptian side experienced massive overgrazing, wood gathering etc., which led to the destruction of vegetation, while the Israeli side was protected from such activities, leaving the vegetation in its natural distribution. Later, it was suggested that the trampling of the biogenic crust cover by Bedouin-tended herds of goats and camels (a practice that was prohibited on the Israeli side) was the primary cause of the brighter (yellowish) sand on the Egyptian side in comparison to the grayish color of the biogenic crust on the Israeli side [
4,
23,
24,
25,
26]. Otterman further suggested that the differences in the soil color led to a difference in the radiative temperature through changes in the surface albedo—the brighter color of the bare sand resulted in a greater reflection of solar radiation, leading to colder temperatures. Furthermore, a positive feedback was suggested in which colder temperatures caused a decrease in precipitation and thus a further decrease in vegetation cover, which led to further cooling, etc. Charney [
12] used a relatively simple atmospheric model to explain the severe drought occurring in the Sahel at that time. In his model, areas with a higher albedo experienced radiative cooling that was connected with other atmospheric processes (such as the Hadley circulation). Charney’s hypothesis was verified using a general circulation model (GCM, [
27]).
Since the work of Charney [
12] and Otterman et al. [
13] positive feedbacks associated with vegetation cover were studied extensively. Such mechanisms were used to explain the severe and prolonged droughts in the Sahel [
28], and the differences between the climate across borders, such as the USA-Mexico [
29,
30] and the Israel-Egypt borders [
31,
32]. More recent studies (e.g., [
33]) have suggested that the severe drought in the Sahel was more closely linked to changes in precipitation rather than anthropogenic activities [
34]. Changes in vegetation cover were linked to paleoclimatic phenomena, including the mid-Holocene warming (∼6000 years ago, [
35]), Taiga-Tundra dynamics [
36], the amplification of glacial-interglacial temperature change [
37], the north African wet phase (the ”green Sahara”) during the early and mid-Holocene, Saharan dynamics [
38], and long-term Saharan change [
35,
39,
40].
Vegetation on dunes is affected by the wind (e.g., [
4,
10]). The main effect occurs via sand transport that either covers vegetation (sand deposition) or exposes roots (sand erosion), leading to vegetation suppression. In addition, the biogenic soil crust, which stabilizes dunes, may also be affected by wind through sand transport that covers the crust, blocking light and eventually leading to crust suppression; wind can also erode crust on the interface between the biogenic crust area and the bare sand area. As sand dunes cover a substantial area of arid and semi-arid regions, it is important to study the mutual feedbacks between climate and dune dynamics. Many studies have examined the effects of vegetation on climate through albedo, temperature, and precipitation–to our knowledge, none have examined the feedbacks between wind and vegetation (and biogenic crust) cover on sand dunes.
Here we propose and examine the following feedback between local climate and vegetation: we start from two neighboring sand dune regions, one with fixed dunes that are covered by vegetation (and biogenic crust) and another with bare and mobile dunes. This distinction results in an albedo difference between the two regions, with warmer temperatures on the vegetated fixed dunes due to their larger absorption of solar radiation. This temperature difference leads to a pressure difference and to air flow from the bare side to the vegetated side, with increased winds on the vegetated side. This enhanced wind activity decreases the vegetation cover, which increases the albedo, etc. In principle, this feedback should eventually lead to the same albedo in the two regions. However, due to the dynamics of vegetation cover on sand dunes, this is not always the case. We studied the proposed feedback in two stages. First, using a state-of-the-art atmospheric general circulation model, the Weather Research and Forecasting (WRF) model [
41], we examined the changes in wind magnitude and pattern when prescribing the albedo in neighboring regions (i.e., we changed the albedo such that the albedo difference between the neighboring regions varied). We found that the wind drift potential (DP) increased approximately linearly in the vegetated side as the albedo difference between the regions increased. We then used this result in a simple model for vegetation cover on sand dunes [
42,
43,
44] to study the effect of albedo on dune vegetation cover (and hence dune mobility), which then yielded an updated wind DP.
The rest of the manuscript is organized as follows. We first (
Section 2) briefly introduce the atmospheric model and the setups we used. We then present the results obtained by this model (
Section 3). Then, we introduce a vegetation cover on sand model that includes the sand albedo feedback and discuss the results we obtained using this model (
Section 4). A summary of the results concludes the paper (
Section 5).
2. The Atmospheric Model, WRF
We used the WRF 3.5 model [
41] in our study. The WRF is a mesoscale numerical weather forecasting model that is used both to study atmospheric phenomena and for operational forecasting. The WRF was used to investigate a wide spectrum of meteorological problems and can handle a wide range of spatial scales, from thousands of kilometers down to tens of meters. It can perform simulations using real observational data or can be run under idealized conditions. The WRF is a Euler non-hydrostatic, fully compressible model. We used the WRF single-moment 6-class micro-physics scheme. The planetary boundary layer scheme was that of Yonsei University and for the cumulus parameterization, we used the Kain-Fritsch scheme.
The first region we studied is centered in the northwestern Negev Desert, along the Israel-Egypt border. In this region, there is a noticeable albedo difference between the active dune Egyptian side and the fixed dune Israeli side (
Figure 1a). The coarser grid resolution was ∼15.4 km, with 70 × 70 zonal and meridional grid points. The resolution of the second grid was ∼5.1 km and had 60 × 60 grid points in the zonal and meridional directions while the innermost grid had a resolution of 1.7 km with 50 × 50 grid points in the zonal and meridional directions. We studied two years, one with relatively weak winds (September 2004 to August 2005) and one with strong winds (September 2002 to August 2003).
The second region we studied is the Sahel/Sahara region (
Figure 1b). We used similar grid resolutions as for the first region with ∼17.5 km, ∼5.8 km, and 1.9 km resolutions and 70 × 70, 60 × 60, and 50 × 50 grid points in the zonal and meridional directions. Here we simulated a 10-day episode, during a storm (17–27 January 1990). Ten days were sufficient to achieve convergence.
To study the effect of albedo difference across two neighboring regions, we ran the simulation using various prescribed albedo values in the regions. We assume that the mean albedo of the region remains the same while the albedo contrast between the two regions is changing; the mean albedo was chosen to be the average of vegetation albedo of 0.1 and bare sand albedo of 0.4, see below. Thus, the mean albedo of the two regions was
where the albedo in the neighboring regions was
where
. We investigate both negative and positive
as the background synoptic activity and nonlinear processes associated with it may break the expected symmetry around
. This is a wide range of albedo and is used to study the effect of albedo difference under extreme albedo conditions. For the Israel-Egypt example, the country border line was chosen to be the albedo border line, and the albedo was changed only in the innermost grid (
Figure 1a and
Figure 2a). For the Sahel/Sahara case, the albedo border line was chosen to be 14.7
N, where here we repeated the numerical experiments when the albedo was changed, first only within the innermost grid and then in all three grids.
3. Results
We show the results of a particular numerical experiment from the Israel-Egypt border, where the albedo is low (0.13) on the Israeli side and high (0.37) on the Egyptian side. These values are in range of the albedo of vegetation (and biogenic crust), which, in most cases, is higher than 0.1, and the albedo of dry sand rarely exceeds 0.4 (see, e.g., [
45,
46]). Below we examine a range of albedo differences. In the current example, the albedo in the innermost grid (within which the simulations are the most realistic) is prescribed while the albedo values of the other grids remain unchanged (
Figure 1a and
Figure 2a). We show the results of the winter season (three winter months), from December 2002 to February 2003, during which the winds were relatively strong.
Figure 2b depicts the mean surface temperature (colors) and the mean wind field (arrows). The surface temperature over the vegetated/crusted Israeli (northern) region is higher by several degrees than the Egyptian (southern) bare dune region. This is a direct consequence of the lower albedo on the Israeli side that leads to greater absorption of solar radiation. In general, the weather during winter in this region is often affected by the Cyprus low system (e.g., [
47,
48]), which usually leads to southwesterly winds. The temperature contrast across the border enhances the winds over the vegetated and warmer Israeli side in comparison to the bare and colder Egyptian side (
Figure 2b). This temperature difference leads to a pressure difference and hence to enhanced flow over the warmer Israeli region in which the air convects upward.
To visualize the effect of the albedo difference, we compared the results discussed above to results obtained when the albedo on both sides of the Israel-Egypt border is the same and equal to 0.25 (the mean albedo of the example discussed above). In
Figure 2c, we present the difference between the surface temperature and the wind field shown in
Figure 2b and the temperature and wind field of the run with the same albedo on both sides of the border. Clearly, the vegetated/crusted Israeli side is warmer than the Egyptian side by a few
C, and the wind is stronger over the warmer Israeli side. These are the winter mean values; during the winter, the sand activity is strongest. Since the motivation of the present study is the dune activity, we also present the difference (relative to the run in which the albedo is the same on both sides of the Israel-Egypt border) in wind shear velocity (
Figure 2d)—wind shear velocity is closely related to sand transport and is calculated by WRF using a logarithmic wind speed profile. The wind shear velocity is enhanced over the warmer, vegetated Israeli side and is suppressed over the colder, bare Egyptian side, indicating enhanced sand activity on the vegetated side. The shear velocity difference is small since we present the winter mean values—as mentioned above, the wind activity during winter in this region is dominated by a few storms while otherwise the wind is weak, leading to a small mean shear velocity. Yet, the winter storms can lead to enhanced sand activity; see the results below regarding the wind drift potential. As expected, the wind shear velocity difference over the sea is about zero.
We next systematically examined the effect of the albedo difference on the temperature and wind.
Figure 3a,b depicts the surface (skin) temperature on the Israeli and Egyptian sides, from September to November and December to February in the years 2002–2003 and 2004–2005, as a function of the albedo difference between the Egyptian and Israeli sides; a specific example of this temperature is depicted in
Figure 2b. In
Figure 3c,d we plot surface air (2 m) temperature and, as expected, the albedo difference effect on the 2 m air temperature is smaller than the effect on the surface skin temperature. It is apparent that, in all cases, both the surface and air temperature increase, almost linearly, as the surface albedo decreases. Similar to the temperature, we plot in
Figure 3e,f the mean wind (10 m height) speed, calculated by WRF, as a function of the Egypt-Israel albedo difference. Also here the relation is linear and the mean wind speed increases as the albedo decreases. This is true in both seasons and during years with both relatively strong and weak winds.
A common measure for the wind sand drift potential (DP) was suggested by [
49] and is defined as
where
U is the wind speed (in knots: 1 knot = 0.514 m/s) measured at a standard height of 10 m, and
is the minimal threshold velocity (=12 knots) necessary for sand transport. The “
” represents the temporal mean. The units of DP are vector units (VU) and, based on DP the wind energy can be classified into three categories: (i) low energy with
; (ii) intermediate energy with
; and (iii) high energy with
. Due to the cubic dependence of DP on the wind speed, extreme wind events are much more dominant, and the dependence of DP on the albedo difference
is not necessarily similar to the dependence of the wind speed
U on
. As for the temperature and wind speed, DP also depends linearly on the albedo difference, with increased wind activity over the low albedo region (
Figure 4a,b). The effect of the albedo difference on DP can be drastic, in some cases doubling the DP. Note that the DP was much larger during the winter months due to winter storms. In addition, storms were stronger in the 2002–2003 winter season than in the 2004–2005 winter season, while during the fall (September to November), the wind was stronger in 2004. Also note that DP can be larger than zero, even if the seasonal mean wind speed (shown in
Figure 3e,f) is smaller than the threshold velocity
, as the DP is based on the mean of the cube of the wind speed and not on the cube of the mean wind speed.
The DP ignores the direction variability of the wind and the Resultant Drift Potential (RDP) together with the Resultant Drift Direction (RDD) provide some information about the wind direction. RDP is defined in a similar way to the DP, where
,
where
is the wind direction and
. RDD is defined as the ratio between RDP and DP,
such that when the wind direction is constant in time
(
) while otherwise it is less than 1. Smaller RDD indicates greater variability in wind direction.
Figure 4c,d depict the RDP, and, as expected, is smaller than the DP (shown in
Figure 4a,b). There is a significant difference between September to November, 2002 and September to November, 2004, where for the former the wind power (expressed by DP) is lower as well as the RDP. In addition, wind direction variability was larger during 2002. As for the winter months, the wind power during December 2002 to February 2003 was significantly larger than that of December 2004 to February 2005. Moreover, wind direction is relatively uniform during the winter months.
In the above example, the albedo changed only within the innermost grid that spans a region of approximately 70 × 70 km; we refrain from altering the albedo over all grids as the coarser ones are not appropriate for albedo changes (especially the northern regions that lacks sand areas). However, winds can intensify when the region is larger. To validate the results obtained for the Israel-Egypt border case and to study the albedo difference effect in a large domain size, we considered the example of the Sahara/Sahel for which we prescribed the albedo in the northern and southern parts of the domains (grids), first in the inner grid only and then in all three grids, changing the land albedo over the entire domain. We considered a 10-day episode (17–27 January 1990) during which strong winds were observed. In
Figure 5, we show the DP as a function of the difference in albedo between the southern and the northern parts. When the albedo was changed only in the inner grid, the DP depended linearly on the albedo difference and increased (decreased) in the lower (higher) albedo region. When the albedo was changed in all three grids, the DP was up to twice as large as the results obtained when the albedo was prescribed only in the inner grid. Here, DP varied linearly in the southern part but, surprisingly, remained roughly constant in the northern region. It is possible that a simulation over the entire season may yield different, more linear relations.
In summary, albedo changes across neighboring regions led to large changes in DP, where DP increased linearly in the lower albedo (vegetated/crusted) region and decreased in the higher albedo (bare) region as a function of the albedo difference.
4. The Effect of Albedo Difference across Regions on Vegetation Cover
After establishing the linear relation between the DP and the albedo difference across two neighboring regions, we then studied the ways that this increased/decreased DP affects the vegetation cover, which, in turn, affects the albedo and hence the DP. For this purpose, we used the vegetation cover model of Yizhaq et al [
42,
43]. For the sake of clarity and simplicity, we used the simplest model version, i.e., the spatially independent model with an ordinary differential equation for vegetation cover, ignoring the effect of biogenic crust [
50,
51] and spatial dynamics [
44].
Here we briefly describe the model and its coupling to winds. We assumed (i) two neighboring regions with specific albedo values and that (ii) the local wind DP is composed of background wind (associated with, for example, synoptic scale activity) plus a linear deviation in DP due to changes in albedo. That is,
where
is the drift potential attributed to background wind,
is the albedo difference between the two regions, and
is a parameter that quantifies the effect of the albedo difference on DP. The albedo in each region is the weighted mean of the relative cover of sand and vegetation within the region and is
for
(the index of the two regions) and where
is the known albedo of vegetation and sand. Thus,
and
where the upper “+” sign refers to the first region (“1”) and the lower “-” sign refers to the second region.
Once the
of each region was determined (as a function of the difference in albedo or the difference in relative vegetation cover), we used the model for relative vegetation cover on sand dunes suggested in [
42,
43]:
where
if
; otherwise
. In addition,
. For more details regarding the model’s construction and parameters, see [
43].
The coupling between the relative vegetation cover of the two areas,
and
, occurs through DP
that depends on the vegetation cover difference,
. Thus, when the vegetation cover in the two regions is equal (
), one receives the steady state solution of one region. This solution is shown in
Figure 6 (green curves) and exhibits bi-stability and hysteresis. Yet, other steady states do exist, and, as expected, these depend more on
(see Equation (
4)), the parameter that quantifies the effect of the albedo difference. The precipitation rate
p affects the steady states less drastically (at least in the parameter regimes we examined). Basically, the coupling of the vegetation with the atmosphere (winds) “widens” the regime of hysteresis behavior, i.e., the range of bi-stability is larger; the vegetation cover of the fixed-dune (vegetated) state becomes larger while the bare-dune state becomes smaller. These new states are associated with the wind field that is different from the one obtained when the albedo difference feedback is ignored. They can also affect the surface temperature of these regions, altering the local climate.
With respect to the Egypt/Israel border case discussed in the previous section (shown in
Figure 3 and
Figure 4), the DP value are low (
Figure 4a,b) and thus, under the current climate conditions, it is likely that the dunes there are associated with the fixed dune (more vegetated) state (see also, e.g., [
43]).
5. Summary
Many places around the globe show regions of active bare dunes that border a region of fixed vegetated dunes (see, e.g., [
42,
43]). This difference in vegetation cover can lead to a difference in surface albedo that further enhances the wind activity on the vegetated side and weakens the wind activity on the bare side. This, in turn, can affect the vegetation cover, as vegetation can be suppressed (directly and indirectly) by the wind [
42,
43]. We studied this hypothesis in two steps. We first used a state-of-the-art atmospheric general circulation model, the WRF, to investigate how an albedo difference across regions affects local climate variables, such as surface wind and temperature; we found that surface temperature, wind speed, and wind drift potential (DP) in the lower albedo region increased linearly as a function of the albedo difference, and decreased linearly in the higher albedo region. We also found that this dependence intensified when the region in which the albedo changed was larger.
We then used this linear relation to study how the altered wind field affects the vegetation cover on sand dunes. For this purpose, we used a simple model for vegetation cover on sand dunes, and we found that the steady states of the model became more extreme (i.e., the vegetation cover of the vegetated state becomes larger while the vegetation cover of the bare state becomes smaller) and that the bi-stability and the hysteresis curve regime were extended. The altered vegetation cover states were associated with a new wind field.
We also examined the change in precipitation due to the albedo difference and found that it is not significant. In addition, for the sake of clarity, we restricted ourselves to a simple vegetation cover model and ignored spatial effects and biogenic crust dynamics. The inclusion of these effects may enrich and alter the dynamics of the front between vegetated and bare sand dune areas.
Local wind has a time scale that is much faster than the time scale of sand dunes. More specifically, winds converge to a statistical steady state (after altering the surface albedo) in a few days. Sand dunes have a much longer adjustment time, more than tens of years. This time scale separation allows us to easily couple the vegetation cover on sand dune model to the atmospheric winds, as in terms of the sand dune time scale, changes in dune cover instantly translated to changes in wind power (DP).
The coupled model we propose exhibits multiple steady states of active and fixed dunes. These were reported previously (e.g., [
42]), but here, we link this bi-stability of sand dunes to the bi-stability of the local atmosphere. This atmospheric bi-stability influences and is influenced by the sand dune states, through the feedback mechanism we propose. We conjecture that such atmospheric bi-stability may have underlain past rapid changes in climate.
The mechanism we proposed here may be applicable to other climate systems. For example, a key component of the climate system is the high latitude sea ice, and many mechanisms in past and present day climate change were linked to sea ice dynamics (e.g., [
52,
53,
54,
55]). It is possible that the contrast between the open sea’s low albedo and the sea ice’s high albedo enhances the winds towards the sea, dragging the sea ice farther into the ocean. This mechanism may accelerate the rapid growth of the sea ice.