Research ArticleCLIMATOLOGY

On the decadal scale correlation between African dust and Sahel rainfall: The role of Saharan heat low–forced winds

See allHide authors and affiliations

Science Advances  09 Oct 2015:
Vol. 1, no. 9, e1500646
DOI: 10.1126/sciadv.1500646


A large body of work has shown that year-to-year variations in North African dust emission are inversely proportional to previous-year monsoon rainfall in the Sahel, implying that African dust emission is highly sensitive to vegetation changes in this narrow transitional zone. However, such a theory is not supported by field observations or modeling studies, as both suggest that interannual variability in dust is due to changes in wind speeds over the major emitting regions, which lie to the north of the Sahelian vegetated zone. We reconcile this contradiction showing that interannual variability in Sahelian rainfall and surface wind speeds over the Sahara are the result of changes in lower tropospheric air temperatures over the Saharan heat low (SHL). As the SHL warms, an anomalous tropospheric circulation develops that reduces wind speeds over the Sahara and displaces the monsoonal rainfall northward, thus simultaneously increasing Sahelian rainfall and reducing dust emission from the major dust “hotspots” in the Sahara. Our results shed light on why climate models are, to date, unable to reproduce observed historical variability in dust emission and transport from this region.

  • Dust emission
  • the Saharan Heat Low
  • Sahelian rainfall
  • wind speed


By mass, aeolian dust is the most pervasive aerosol in the atmosphere (1). The presence of suspended dust affects the local energy balance through direct (2) and indirect effects (3), modifies the hydrological cycle via radiative forcing (4) and modification of cloud microphysical properties (5), and is a vehicle for the long-range transport of nutrients to global oceans (6) and terrestrial land surfaces (7). However, despite the importance of dust to the Earth system, coupled climate models are not able to reproduce historical variability in dust emission and transport (8), and thus, there is little confidence in projections of how atmospheric dust concentrations will change in the future. Hence, there is a need to understand the controls on historical changes in dust using observational data.

North Africa is the world’s largest dust source, accounting for more than 55% of global dust emissions (9, 10), and nearly 60% of the dust transported off the continent is deposited in the Atlantic Ocean (1). Long-term in situ observations at Barbados (11) and coral reef proxy and satellite data (12, 13) have shown that dust emission from North Africa peaked during the mid-1980s and has followed a downward trend through at least the late 2000s. Studies have found that dust cover over the Atlantic is anticorrelated with previous-year Sahelian precipitation (11), and a causal relationship between the decline in dust emission and the increase in Sahelian rainfall is plausible because increasing soil moisture would strengthen soil cohesion forces (14), and more dense vegetation cover would increase surface roughness (15), both acting to limit erosion and the mobilization of particles from the surface into the atmosphere. However, there are several aspects of such a theory that are not consistent. For example, satellite imagery suggests that most dust-emitting regions lie to the north of the vegetated region of the Sahel (10, 16), and in situ observations have shown that a more northward propagation of the monsoon, which results in greater Sahelian rainfall, should increase the occurrence of dust emission from sources in the Sahel via haboobs (17). Additionally, a recent study demonstrated that the downward trend in dust that started in the mid-1980s could be reproduced in an aerosol emission and transport model when forced only by historical wind fields from reanalysis and when holding the Sahelian vegetation density constant (18).

Therefore, on the one hand, it would appear as though changes in surface winds over the hyperarid Sahara Desert should be sufficient to describe year-to-year changes in summertime dust emission, yet on the other hand, there is clear evidence showing that dust transported across the Atlantic is inversely proportional to rainfall, as well as vegetation density, in the Sahel. Here, we reconcile this contradiction by showing that both Saharan surface wind fields over the major regional dust-emitting regions and the northward propagation of the monsoon flow, and thus Sahel rainfall, are forced by the thermodynamic state of a meteorological feature termed the Saharan heat low (SHL).


Dust and precipitation

The Sahelian region of North Africa is an arid to semiarid region lying between the Sahara and the more tropical regions near the Guinea coast that is characterized by average summertime precipitation rates from 2 to 6 mm day−1. Numerous studies have suggested that Atlantic dust cover is sensitive to rainfall in the Sahel, given the negative correlations between dust and monsoon season precipitation (11, 17). Indeed, the time series of summertime, June to August, the season when dust production peaks (19), with tropical North Atlantic dust aerosol optical depth (τd) estimated from satellites (12), shows that τd has been decreasing since the mid-1980s, whereas Sahelian rainfall, from the Global Precipitation Climatology Project (20), has been increasing over the same time period (Fig. 1). The time series of monsoon season Sahelian rainfall and satellite-retrieved τd, over the period 1982–2010, are statistically significantly correlated at an r value of −0.53 (P = 0.014). Using dust concentration measurements from Barbados (11) and the same rainfall data, but over the period 1965–2009, the correlation between dust and precipitation is −0.52.

Fig. 1 Seasonal (June to August) time series of Zhl (blue), precipitation (P, green), and τd (red).

The monthly mean precipitation series is averaged over 5°W to 20°E, 10° to 20°N, and τd is averaged over 10° to 65° W and 0° to 30°N. The thin lines are the seasonal means, and the thick lines are the five-season smoothed time series (via a 1-4-7-4-1 filter).

The basis for the theory that Sahelian rainfall affects dust emission is that although the correlation between concurrent-year rainfall and dust is statistically significant and negative, the magnitude of the correlation between previous-year rainfall and dust is larger in magnitude (11) (−0.68 using satellite data and −0.66 using the Barbados data). However, on the basis of the results from a Fisher’s z test to determine the significance of the difference in correlation coefficients, we cannot rule out the null hypothesis that the lagged correlations are statistically indistinguishable from the nonlagged correlations. That is, the fact that the correlation magnitude peaks when Sahelian rainfall leads dust by 1 year may simply be a statistical artifact.

The Saharan heat low

We next examine the concurrent roles of the SHL in dust emission and rainfall in the summer. The SHL is a region of a local surface and lower tropospheric temperature maximum and surface pressure minimum that is present in the western Sahara Desert during the summer (June to August) months, and is characterized by a broad low-level cyclonic circulation and mid-level anticyclonic circulation, about the SHL center (21) (Fig. 2A). The magnitude and extent of these two circulation features can be described by the atmospheric thickness between the 700- and 925-hPa levels, averaged spatially over the SHL center (Zhl), using geopotential height from ERA-Interim reanalysis (22). Larger values of Zhl indicate a strengthening of the characteristic SHL circulation features (21). The SHL center is detected in a manner identical to the Lagrangian SHL tracking method described in (21). The temperature and thickness of the SHL change on time scales ranging from synoptic (23) to interannual (24), and over the last three decades, Zhl has increased (Fig. 1), a trend that is corroborated by in situ observations (24).

Fig. 2 Climatology and regressed precipitation and winds.

(A) Shown is a map of long-term mean dust emission (shading), seasonal mean precipitation (magenta contour lines, mm day−1), and seasonal mean 925-hPa climatological wind vectors (June to August). Dust emission rates are in nondimensional units of % of total annual mean of North Africa dust emission rates. The red circle indicates the location of the mean seasonal (June to August) position of the SHL. The brown boxes are the major dust hotspots of the Bodélé depression (14° to 24°E, 14° to 21°N), the depression in the lee of the Aïr and Adrar Mountains (2° to 10°E, 16° to 24°N), and the Mauritania and Western Sahara source region (12° to 2°W, 18° to 27°N, dashed lines). (B) Map of the coefficients of the regression of 925-hPa meridional and zonal winds (vectors and shading) and precipitation onto Zhl, for June to August and 1979 through 2012. Vectors are all unit length to indicate anomalous wind direction, and the shaded regions are the wind speed coefficients, where only coefficients statistically significant at the 90% level are shown.

As the SHL warms and Zhl increases, there is an intensification of the low-level cyclonic circulation that has its center of action about the SHL (23), and this strengthened low-level circulation increases Sahelian rainfall during the monsoon season via a northward displacement of the monsoon flow (24). Covariability of Zhl and precipitation is evidenced by a statistically significant correlation of their seasonal time series (r = 0.54, P = 0.01) (Fig. 1). This covariability is also seen in Fig. 2B, which shows coefficients of the regression of monthly precipitation anomalies onto monthly mean Zhl anomalies, for the summer months (June to August) and from 1979 through 2012. Here, given an anomalously warm SHL, there is an increase in Sahelian precipitation and a reduction in rainfall along the Gulf of Guinea (Fig. 2B), indicating an anomalously northward position of the monsoon rain band, which is consistent with other studies (24). The dependency of the Sahelian rainfall on the intensity of the SHL is also evident in the regression of 925-hPa winds onto Zhl, also for the months of the monsoon season, which show an anomalous low-level southwesterly flow between the Gulf of Guinea and the Sahel (Fig. 2B) that intensifies the climatological southerly monsoon flow (Fig. 2A) and thus increases the northward transport of moisture to the Sahel (24).

Simultaneous with the increase in Sahel rainfall that is associated with an anomalously warm SHL, there is also a reduction in wind speeds over the three most important dust-emitting regions of the Sahara: the Bodélé depression, a depression in the lee of the Aïr and Adrar Mountains, and western Mauritania and the Western Sahara (10, 16, 24, 25) (Fig. 2A). We note that although the dust emission map from (25) does not show high (relative) emission rates over Mauritania and the Western Sahara, another independent estimate of North African dust emission (10) suggests that this region has the highest emission rates for the entire region.

For all three major dust source regions, the climatological low-level winds are northerly to northeasterly (Fig. 2A), but the anomalous wind fields associated with an increase in Zhl are southwesterly (Fig. 2B), a result of the anomalous cyclonic rotation about the SHL center of action. Consequently, as Zhl increases and the SHL circulation strengthens, there is a net reduction in the scalar 925-mbar wind speeds over these important source regions, as evidenced by the statistically significant coefficients of the regression of wind speed onto Zhl (shading, Fig. 2B). Thus, the regressions of monsoon rainfall and winds onto Zhl (Fig 2B) demonstrate that changes in the low-level circulation associated with the SHL can simultaneously increase Sahelian rainfall and reduce North African dust emission, via the well-established relationship between surface wind speeds and emission [for example, (14)].

Surface wind speeds during strong and weak SHL phases

To more clearly link the dynamics of the SHL to dust emission over the three major source regions (Fig. 3A), we created histograms of reanalysis (22) of six hourly 10-m wind speed for the 300 largest and 300 smallest daily values of Zhl (Fig. 3, A to C). The surface wind speed threshold for dust emission to occur in the Sahara is about 5 to 7 m s−1 (15, 26, 27), and thus only changes in the wind speed distribution within and above this threshold range have an effect on dust emission. The climatological 10-m wind speed histograms for the Bodélé (Fig. 3A), Aïr and Adrar Mountains lee region (Fig. 3B), and Mauritania and Western Sahara (Fig. 3C) source regions all peak at speeds below the threshold for emission. Hence, besides changes in the mean wind speed, changes in tails of the distribution are important for understanding the winds’ influence on dust emission.

Fig. 3 Wind speed histograms.

(A to C) Histogram counts (per 0.2 m s−1) of 10-m wind speed from ERA-Interim reanalysis for the 300 strongest (red) and weakest (blue) phases of the SHL in the Bodélé Depression (A), a depression in the lee of the Aïr and Adrar Mountains (B), and the Mauritania and Western Sahara (C) in the summertime from 1979 through 2012. The mean values of the wind speed for the strong and weak SHL are statistically significantly different in (A), (B), and (C) based on the t score from a two-tailed t test. (D to F) Wind speed histogram differences for the warm minus cool SHL cases. The gray-shaded region represents wind speed threshold for dust emission. A qualitatively identical shift in the wind speed distribution is also observed during the month of May (not shown).

In all three cases, the difference between the two wind speed histograms (large minus small Zhl) demonstrates that there is a shift in the wind speed distributions over the source regions toward lower wind speeds for larger Zhl (Fig. 3, D to F). Quantitatively, there is a 22.5% (Bodélé), 6.5% (Aïr and Adrar), and 4.1% (Mauritania and Western Sahara) increase in the frequency of occurrence of wind speeds below the 5 m s−1 threshold when comparing large and small Zhl. We find a qualitatively identical shift in the wind speed distribution during the month of May. This reduction in dust emission with an anomalously warm SHL and strong SHL circulation is reflected in the time series of τd and Zhl (Fig. 1), which are statistically significantly correlated (r = −0.73, P < 0.01).


Our results suggest that the long-noted negative correlation between dust and rainfall is by association and that both features are simultaneously responding to large-scale forcing by the SHL. Such an interpretation of the data is consistent with previous work regarding the influence of the SHL on the physical characteristics of the West African monsoon (21, 23, 24), and with the Ridley et al. model (18) and observational studies (10, 16, 25) of North African dust sources. We cannot exclude the possibility that dust emission in Sahel is causally correlated to rainfall via changes in vegetation cover and soil moisture [for example, (15)]; rather we suggest that such effects are to first-order negligible when compared to the coincident changes in dust emission over the vastly larger Sahara.

Finally, when forced with observed sea surface temperatures, climate models are unable to reproduce the observed year-to-year variability of dust over the North Atlantic (8). In light of this work, it is plausible that an underlying cause of this deficiency is an inaccurate representation of the location of the major dust-emitting regions (8). Additionally, initial work (4) suggests that climate models may not sufficiently represent the SHL and its influence on the regional climate, which we have shown is requisite to represent the dust cycle in this part of the world, although more work is needed to comprehensively evaluate the representation of the SHL in climate models.



Monthly mean τd, retrieved from satellite radiance measurements from the Advanced Very High Resolution Radiometer, is used to construct the time series of τd over the northern tropical Atlantic for the summer months (June to August) from 1982 through 2009 (12). In Fig. 1, τd is averaged over 10° to 65°W and 0° to 30°N. We use monthly mean precipitation data from the Global Precipitation Climatology Project Version 2 to construct the time series of the Sahel precipitation (21). The precipitation series is averaged over 5°W to 20°E, 10° to 20°N in the summer months (June to August) from 1979 through 2012. We use geopotential height at 925 and 700 hPa at 0600 UT from ERA-Interim daily product to calculate the Zhl, which is the difference in the geopotential heights at these two pressure levels over North Africa (22). Dust emission rates (16) were derived from Spinning Enhanced Visible and Infrared Imager infrared dust index images from March 2006 to February 2008. Last, 925 hPa and 10-m meridional and zonal winds at 0600 UT are from the ERA-Interim daily product from 1979 through 2012. Although wind speeds in the data-sparse Sahara are subject to error, a recent comparison of reanalysis and observed surface winds within the Sahel demonstrated that the ERI-Interim is best suited for analysis of dust emission processes there (28).

Statistical analysis

Correlation coefficients and significance of correlations (using a two-tailed t test) for the time series in Fig. 1 were calculated using the seasonal (June to August) means of each variable and over the time period 1982–2009. We used a Fisher’s z transformation to test the statistical significance of the difference between the correlation coefficients of previous- and concurrent-year Sahelian rainfall and τd. More specifically, this test provides a method of obtaining uncertainty estimates on correlation coefficients, based only on the degrees of freedom in the data. For the zero-lagged correlation between Sahel rainfall and dust, using the 43-year-long Barbados dust data, the correlation coefficient and its 95% confidence interval is 0.53 ± 0.22, and the same quantity for the 1-year lagged correlation is 0.66 ± 0.18.

In Fig. 2B, the anomalous wind field was obtained by separate regressions of the meridional and zonal 925-mbar 0600 UT winds onto daily Zhl, using data from the months June to August and the period 1979–2012. Notably, the results in Figs. 2B and 3 (D to F), which were based on daily data, were obtained without first removing the annual cycle in the various data sets examined. However, we were able to reproduce all results after first removing the annual cycles in these data (not shown).

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We would like to thank two anonymous reviewers for their helpful comments on an earlier version of this manuscript. Funding: This work benefitted from the support of the Agence Nationale de la Recherche (ANR) grant ANR-10-LABX-18-01 of the national Programme Investissements d’Avenir. Funding for this work was also provided by Laboratoire d’excellence Institut Pierre Simon Laplace (L-IPSL). Author contributions: W.W. and A.T.E. designed the study. W.W. carried out the main analysis and wrote the manuscript. C.L. provided historical heat low data. All authors discussed the main findings and contributed to the manuscript writeup. Competing interests: The authors declare that they have no competing interests. Data and materials availability: Reanalysis data from the European Centre for Medium Range Weather Forecasts was downloaded from Global Precipitation Climatology Project data were downloaded from the National Climatic Data Center Web site (
View Abstract

Stay Connected to Science Advances

Navigate This Article