Search tips
Search criteria 


Logo of ncommsLink to Publisher's site
Nat Commun. 2017; 8: 14966.
Published online 2017 May 25. doi:  10.1038/ncomms14966
PMCID: PMC5529672

Skilful prediction of Sahel summer rainfall on inter-annual and multi-year timescales


Summer rainfall in the Sahel region of Africa exhibits one of the largest signals of climatic variability and with a population reliant on agricultural productivity, the Sahel is particularly vulnerable to major droughts such as occurred in the 1970s and 1980s. Rainfall levels have subsequently recovered, but future projections remain uncertain. Here we show that Sahel rainfall is skilfully predicted on inter-annual and multi-year (that is, >5 years) timescales and use these predictions to better understand the driving mechanisms. Moisture budget analysis indicates that on multi-year timescales, a warmer north Atlantic and Mediterranean enhance Sahel rainfall through increased meridional convergence of low-level, externally sourced moisture. In contrast, year-to-year rainfall levels are largely determined by the recycling rate of local moisture, regulated by planetary circulation patterns associated with the El Niño-Southern Oscillation. Our findings aid improved understanding and forecasting of Sahel drought, paramount for successful adaptation strategies in a changing climate.

The Sahel, which we define as the region between 16° W to 36° E and 10° N to 20° N, extends across the African continent, transitioning between two distinct climatological regimes: the Sahara desert to the north and the lush tropical rainforests to the south. Dry conditions pervade the Sahel throughout most of the year, with a wetter period between July and September associated with a northward shift of the main tropical rainfall band1,2. The Sahel summer rainy season exhibits one of the largest signals of global climatic variability3, reflected in the fragility of agricultural productivity in north Africa. Prolonged drought ravished the region in the early 1980s, shattering crop and livestock farming and resulting in an estimated 100,000 deaths in rural communities from starvation, malnutrition and disease4. Sahelian summer rainfall has since recovered to the levels of the 1960s3, but future projections are uncertain5. Superimposed on this multi-decadal change are substantial swings in inter-annual rainfall levels2. Understanding the physical drivers of both inter-annual and multi-year Sahelian rainfall variability is crucial to gain confidence in future predictions.

Climate models vary widely in their ability to capture Sahel rainfall anomalies. There is some statistical evidence of multi-year predictability, especially using multi-model means, for timescales greater than about 4 years and in the Western Sahel (Supplementary Table 1)6,7,8,9. No previous studies have documented predicting summer Sahelian rainfall from the previous winter, although seasonal predictions with lead times of 0–3 months have had considerable promise10,11,12,13,14,15,16,17. Model initialization is reported to improve both the amplitude and predictability of Sahel summer rainfall anomalies, although the results are model dependent7,9,18,19. Improved skill is also noted for models, which better simulate teleconnections between Sahel rainfall and sea-surface temperatures (SSTs)20,21. However, to our knowledge, no studies have clearly demonstrated skilfull predictions of the mechanisms by which global SSTs influence the moisture budget of the Sahel. Here we present improved understanding and prediction of Sahel rainfall variability on both multi-year and inter-annual timescales for lead times >8 months, using a comprehensive set of retrospective forecasts (hereafter hindcasts) from the latest Met Office Decadal Climate Prediction system, DePreSys3 (ref. 22). Covering the period since 1960, DePreSys3 uses a much higher resolution (~60 km in the atmosphere and 0.25° in the ocean) than previous versions23,24, enabling better representation of the key physical processes25 (Methods).


Skilfull prediction of Sahel summer rainfall

The DePreSys3 hindcasts skilfully predict post-1960 summer (that is, the mean of July, August and September) rainfall levels across the Sahel on both inter-annual and longer timescales (Fig. 1). For linearly de-trended mean summer rainfall averaged over the 2–5 year period after initialization (Year 2–5), correlations with observations are 0.73, 0.75 and 0.63 for the whole, west and east Sahel regions, respectively, and all are significant to the 95% level (see Methods). The west and east Sahel are delineated by 10° E. At shorter, 8 month leads, the skill of DePreSys3 for predicting the upcoming summers (Year 1) is more modest, yet still significant: r=0.48, 0.50 and 0.33 for the whole, west and east Sahel regions, respectively. The drought period of the 1980s and the subsequent recovery are particularly well captured by the model, alongside some of the major inter-annual peaks and troughs in rainfall: for example in 1997–98 and 2012 (Fig. 1c,d). Accounting for ensemble size, DePreSys3 shows high skill compared with other initialized models and even the Coupled Model Intercomparison Project (CMIP5) multi-model mean hindcasts, in particular on inter-annual timescales (Supplementary Table 1). However, understanding the physical mechanisms that modulate Sahel rainfall patterns is key for gaining further confidence in future predictions. We now assess these in detail using DePreSys3.

Figure 1
Skilful predictions of summer rainfall in the Sahel.

Circulation patterns associated with Sahel rainfall change

The Sahel sits at the confluence of several atmospheric circulation features that are illustrated schematically in Fig. 2. Examination of the summer climatological meridional flow reveals two areas of peak convergence and ascent (Fig. 3a): a deep circulation cell at 10°N, co-located with the Tropical Easterly Jet (TEJ) and the core of the tropical rain belt, and a shallower meridional overturning cell further to the north (~17°) that coincides with the surface manifestation of the Inter-tropical Convergence Zone (ITCZ)2,26. Zonally, Sahelian surface ascent peaks at the West African coast, just east of where the easterly Atlantic trade winds morph into the onshore West African Westerly Jet (WAWJ)2,27 (Fig. 3b). This low-level flow, alongside a band of mid-atmospheric upwelling and upper level easterly flow spans much of the African continent. These atmospheric climatological patterns are consistent with several re-analysis products, although the deep circulation cell does not extend quite as far to the north in DePreSys3 (Supplementary Note 1 and Supplementary Fig. 3).

Figure 2
Schematic of climatological atmospheric features that modulate Sahel summer rainfall.
Figure 3
Atmospheric circulation associated with Sahel rainfall.

To distinguish the character of Sahelian rainfall drivers over different timescales, we analyse a large set of DePreSys3 1 year lead time hindcasts, separated into multi-year (that is, >5 years) and inter-annual timescales (Methods). Resultant timeseries are divided into wet and dry Sahel composites (Methods). Differences between wet and dry Sahel summer circulation patterns in DePreSys3 for both multi-year and inter-annual variability are shown in Fig. 3c–f. On multi-year timescales the deep convection cell is strengthened slightly and shifted northward by ~3.5° latitude during wet periods, displacing the main rainbelt further into the Sahel28. On inter-annual timescales, wet summers are modulated by a strengthened (as opposed to displaced) deep convection cell, coincident with a strongly enhanced zonal Walker-type circulation29,30. A poleward shift in the shallower circulation cell (that is, at 17°N) is apparent with increased rainfall on both timescales. As the lower cell meridional shifts however remain within the Sahel in DePreSys3, they will not impact domain-averaged rainfall levels. Similar analysis for re-analysis data are discussed in the Supplementary Note 1.

Sahelian summer atmospheric moisture fluxes

We next investigate changes in the moisture fluxed into the Sahel region at different atmospheric levels, that is, the summer mean specific humidity, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m1.jpg multiplied by the wind, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m2.jpg. Wet-minus-dry composites are shown in Fig. 4. Moisture flux changes are dominated by atmospheric levels below 400 hPa, reflecting the higher humidity content here. On both inter-annual and multi-year timescales, wetter summers are associated with increased moisture from the following: stronger northward winds between 600 and 850 hPa along the southern Sahel boundary; increased humidity carried by the Mediterranean sourced surface winds (even though they actually weaken slightly)—this is particularly important in the east of the region31 Supplementary Fig. 7; and a stronger WAWJ. Moisture flux anomalies associated with the mid-altitude African Easterly Jet32 shows a much stronger modulation on multi-year timescales and mostly counteract those from the lower level westerlies (Supplementary Note 2).

Figure 4
Changes in moisture between wet and dry Sahel summers in DePreSys3.

Analysis of column integrated moisture fluxes, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m3.jpg, highlights the role of meridional circulations in modulating Sahel rainfall. The humidity of southward flowing surface winds from the Mediterranean and the strength of the northward monsoon flow dominate the overall moisture convergence in the Sahel, whereas little convergence is associated with the zonal WAWJ and African Easterly Jet components (Fig. 4). Changes in moisture convergence between wet and dry periods are similar on multi-annual and inter-annual timescales (0.17 and 0.16 mm per day, respectively, Fig. 4). However, the total rainfall change is higher on inter-annual timescales (0.28 mm per day compared with 0.23 mm per day for multi-year timescales). This difference is accounted for by a stronger contribution from evaporation (0.12 mm per day compared with 0.06 mm per day for multi-year timescales).

External moisture supply versus local moisture recycling

The relative roles of moisture advected into the region Pa versus the recycling of local moisture Pm are further assessed by computing timeseries of recycling ratios ρ=Pm/P, where P=Pa+Pm is the total rainfall (see Methods)33. Sahel recycling ratios in DePreSys3 range from 0.3 to 0.4, suggesting that recycled water typically contributes just over a third of the total summer Sahel precipitation. However, we find a distinct contrast between the character of rainfall variability on different timescales: multi-year precipitation rate anomalies are weakly anti-correlated with the recycling ratio, whereas inter-annual fluctuations show a strong positive relationship: r=0.77(0.62, 0.88), where numbers in brackets represent the 5–95% confidence intervals. (Fig. 4). These relationships suggests the dominant role for dynamically driven local moisture recycling on inter-annual timescales, with external moisture supply being more important on multi-year timescales.

To further investigate the apparently different mechanisms driving Sahel rainfall change on different timescales, we compute the moist static energy (MSE), which enables the relative roles of moisture and temperature changes to atmospheric stability (and hence ascent) to be assessed. The MSE represents the sensible, latent and geopotential energy, such that MSE=cpT+Lq+gz, where cp is the specific heat of air at constant pressure, T the temperature, L the latent heat of vapourization of water, q the specific humidity, z the height and g the acceleration due to gravity34. An increasing MSE with altitude denotes a stable atmosphere. On multi-year timescales, increased moisture content (Lq) during wet years reduces the vertical stability throughout the atmosphere, but most intensely below 500 hPa (Fig. 5a). By comparison, on inter-annual timescales, it is the upper atmosphere (above 500 hPa), which is de-stabilized by cooler upper atmospheric temperatures (cpT) and therefore ascent in the deep convection cell enhanced, whereas at lower levels cooler temperatures counteract the increased moisture content (Fig. 5b). These results are consistent with the apparent changes in circulation patterns and moisture fluxes between wet and dry Sahel summers on different timescales (Figs 3 and and4):4): meridional shifts in the circulation and moisture transport govern rainfall change on multi-year timescales, whereas upper-atmosphere temperature changes impact upper-atmosphere moisture ascent on inter-annual timescales.

Figure 5
Drivers of atmospheric stability in DePreSys3.

SSTs and teleconnections

To gain further confidence in the skill of the forecasts, we now assess the global drivers of moisture flux and local recycling that were found to be important on the two timescales. Several studies discuss the high sensitivity of the Sahel rainy season to the tropical and north Atlantic, the tropical Pacific, the Indian Ocean Dipole, the Mediterranean and the Saharan Low3,18,20,31,35,36,37,38,39,40,41,42. Similar results are found in DePreSys3: on multi-year timescales, a strong relationship exists between north Atlantic and Mediterranean SSTs and Sahelian summer rainfall (Fig. 6a), whereas the eastern tropical Pacific, eastern equatorial Atlantic and Indian Ocean are significantly correlated with Sahel precipitation inter-annually (Fig. 6b). Considering these teleconnections, we now return to consider the source of the skill in the year 2–5 and year 1 lead-time hindcast sets (Fig. 1). DePreSys3 is able to skilfully predict north Atlantic and Mediterranean SST on multi-year timescales, and that of the tropical Pacific (that is, the El Niño-Southern Oscillation (ENSO)) and the western Indian Ocean (but not Indian Ocean Dipole) for year one hindcasts (Fig. 6c,d). However, we do not find significant skill in predicting the Saharan Low. Thus, we conclude that the north Atlantic and the Mediterranean on multi-year timescales, and the equatorial east Pacific and west Indian ocean on inter-annual timescales, are the fundamental sources of Sahel rainfall skill in DePreSys3 and evaluate the key mechanisms in more detail below.

Figure 6
Relationship between global surface temperature patterns and Sahel rainfall.

On multi-year timescales, a comparatively warm north Atlantic acts to displace the marine ITCZ northward24,43,44, shifting the tropical rainbelt and the associated meridional moisture convergence into the Sahel region (Figs 3c, ,4a4a and and7).7). Warmer northern hemisphere temperatures also lead to increased evaporation over the Mediterranean and the Atlantic Ocean, enhancing the specific humidity of air advected into the Sahel from these regions (Fig. 4a)31. In addition to elevating moisture convergence, this increased low-level moisture content over the Sahel acts to destabilize the atmosphere below 400 hPa, promoting convection and precipitation, particularly in the lower atmosphere (Fig. 5a). Similar patterns in multi-year circulation and MSE anomalies are seen for composites based on north Atlantic SST (25–60° N, 7–70° E; Supplementary Figs 8c,d and 9). In support of these processes underpinning the skilful prediction of multi-year rainfall change, we find that DePreSys3 year 2–5 hindcasts are able to forecast the precipitation, the low-level humidity and the wind anomalies associated with the sustained droughts of the 1970s and 1980s (Fig. 8). The equatorward migration of the marine ITCZ during the drought period is reflected in the north–south dipole pattern of rainfall anomaly maps (Fig. 8a,b)38 and anomalous northerly winds fluxing moist air out of the Sahel across its southern boundary (Fig. 8c,d). Regions of reduced surface humidity associated with the drought are concentrated where winds originating from anomalously low humidity regions in the Mediterranean and Atlantic Ocean converge (Fig. 8e,f). Furthermore, we find that meridional shifts in the marine ITCZ, assessed as the latitude where 925 hPa meridional winds averaged between 50° W and 16° W reach zero magnitude, are skilfully predicted at the 95% level for year 2–5 hindcasts compared with re-analysis (r=0.57). Low-level (that is, 850 hPa) specific humidity changes over the Sahel are also skilful in the year 2–5 hindcasts, with a correlation of 0.89, significant at the 95% level.

Figure 7
Anomalous circulation patterns of the marine ITCZ associated with Sahel rainfall change.
Figure 8
DePreSys3 predictions of the prolonged droughts of the 1970s and 1980s and the underlying mechanisms.

The influence of ENSO and the western Indian Ocean on Sahelian rainfall is largely confined to inter-annual timescales (Fig. 6b). La Nina events are known to strengthen upper level easterlies (that is, the TEJ) across the Sahel through intensified surface temperature gradients, both zonally across the Pacific/Indian oceans and between the Tibetan Plateau and India30,42,45 (Figs 3e,f and 9a,b). These upper level wind anomalies act to dynamically strengthen the deep convective upwelling cell in the Sahel through enhanced north–south wind divergence2. Furthermore, La Niña events and a cooler west Indian ocean destabilize the troposphere over the Sahel: anomalous downwelling and reduced latent heat release over the tropical Pacific and Indian oceans cool the global tropical upper atmosphere46,47,48, resulting in reduced atmospheric stability between 200 and 500 hPa (Fig. 5b). Note also the cooler surface temperature across the Sahel associated with La Niña events (Fig. 9a,b), which acts to stabilize the lower atmosphere (Fig. 5b). Composites based on an ENSO index (SST in −5 to 5° N, −170 to −120° W) or west Indian Ocean SSTs (−10 to 10° N, 50 to 70° W), produce similar patterns to Figs 3e,f and and5b5b (Supplementary Figs 8e,f and 9b). These two mechanisms, which are both skilfully predicted by DePreSys3 (Fig. 9c,d and e), amplify deep convection in the Sahel, leading to increased moisture recycling on inter-annual timescales (Fig. 4b).

Figure 9
Skilful predictions of the mechanisms which modulate Sahel rainfall on inter-annual timescales.


We have shown that Sahel rainfall and the associated driving mechanisms are predictable on both multi-year and inter-annual timescales. However, we note that climatological rainfall levels and the amplitude of variability in DePreSys3 are too low, particularly poleward of 10° N (Fig. 10). Other modelling studies commonly report similar issues12,49,50. We also note that the magnitude of moisture flux anomalies in the re-analysis data are generally greater than those in DePreSys3 (Fig. 5 and Supplementary Fig. 4).

Figure 10
The climatology and variability of Sahel rainfall in DePreSys3 compared with observations.

We investigate the lack of rainfall in DePreSys3 by analysing the first empirical orthogonal function (EOF1) for Sahel precipitation, which accounts for ~30% of the variance in both the observations and DePreSys3. This analysis highlights that peak precipitation variability in DePreSys3 is about 2.5 (2) degrees too far south for hindcasts with year 2–5 (year 1) lead times (Fig. 10d,f). These EOF1 patterns also explain the sharp drop-off in correlation to the south of the Sahel (Fig. 1), which corresponds to opposite signs in the Global Precipitation Climatology Centre (GPCC) and DePreSys3 EOF1 patterns. To investigate whether the southward displacement of rainfall variability in DePreSys3 could account for the low mean and variance in modelled rainfall levels in the Sahel region, we re-compute rainfall timeseries for a box shifted by the appropriate degrees latitude, as evaluated using EOF1 patterns (Fig. 10d,f, orange boxes). We find that this shift accounts for the low climatological Sahel rainfall levels in DePreSys3, but does not wholly account for the low modelled variability. For example, for lead times of 2–5 years, the mean of each ensemble member standard deviation, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m4.jpg, is 0.16 for the shifted box compared to 0.31 for the Sahel GPCC timeseries. Furthermore, the ratio of predictable components51 in observations relative to the model is greater than one (1.18) for DePreSys3, in agreement with results from CMIP5 multi-model ensemble51. This result suggests that Sahel rainfall is constrained too weakly by the predictable signal relative to the noise. Similar issues are seen for seasonal forecasts of the North Atlantic Oscillation52. Thus, the discrepancy in climatological rainfall levels is potentially explained by the main rainbelt in DePreSys3 not extending far enough north. However, even after correcting for this, the variability and signal to noise ratio of the model remains too low. We speculate that other factors, such as the Saharan Low, act to amplify Sahelian rainfall36 in reality but are not represented properly in DePreSys3. Resolving these issues is beyond the scope of the this study, but could potentially lead to better skill in future.

In summary, this study demonstrates the ability of initialized climate models to skilfully predict and better understand both periods of sustained multi-year Sahel drought and year-to-year summer rainfall variability. Skilfull predictions of the underlying driving mechanisms increases our confidence that the high statistical skill represents physically realistic processes. On multi-year timescales, meridional moisture convergence plays a leading role in regulating Sahelian rainfall levels. Prolonged droughts are the result of equator-ward migration of the tropical rainbelt, weakened southerly monsoon winds and reduced low-level moisture content, in response to a cooler north Atlantic and Mediterranean Sea. On inter-annual timescales, the imprint of ENSO and western Indian ocean SSTs on tropospheric temperatures and zonal circulation modulates the ascent and recycling of local Sahel moisture through changes in upper-atmospheric stability and upper troposphere meridional wind divergence. These findings may also have implications for how the character of the Sahelian summer rains respond to global climate patterns: since tropical Pacific temperatures predominantly regulate the vertical stability and rate of local moisture recycling or threshold for deep convection, their impact may be largely manifest in the frequency of rainy days in the Sahel; north Atlantic and Mediterranean SSTs, on the other hand, may be expected to affect the median intensity of rainfall events through their modulation of moisture supply. In fact, data from Senegal indicates that the persistent droughts of the 1970s and 1980s, and the subsequent recovery, are associated with changes in the median intensity of daily rainfall events as opposed to number of rainy days29,53, consistent with our expectations from multi-year changes driven by north Atlantic and Mediterranean SSTs. Our findings contribute to a deeper understanding of Sahel rainfall change and show that skilful predictions are possible months to years ahead. Such results are paramount to building resilience for adaptation to climate variability and change, improving future food security and economic stability in this vulnerable region.



DePreSys3 is based on the Hadley Centre Global Environment Model version 3, HadGEM3-GC2 (ref. 54) and is initialized by relaxing towards observed analyses of: ocean temperature and salinity from global covariance analysis55 (nudged monthly with a 10-day relaxation timescales); sea-ice concentrations from the HadISST data set56 (monthly with 24 h relaxation); and ERA-interim57 atmospheric temperature and winds (nudged 6 h with 6 h relaxation). Data are assimilated from the surface and ramped up to 1 km altitude. The atmospheric resolution is 60 km with 85 quasi-horizontal atmospheric levels. The oceanic resolution is 0.25° with 75 levels. Hindcasts are fully influenced by external forcings (for example, greenhouse gases, aerosols, ozone, solar radiation and volcanoes) as per the CMIP5 protocol58. Ten ensemble members are run up to 16 months from initiation, for every year since 1980, and roughly every 2–3 years between 1960 and 1981 (1960, 1962, 1965, 1968, 1970, 1972, 1975, 1978 and 1980): 42 start dates in total. A smaller subset of these ten ensemble member model runs were continued to 5 years from intialization (1960, 1962, 1965, 1968, 1970, 1972, 1975, 1978, 1980, 1982, 1985, 1988, 1990, 1992, 1995, 1998, 2000, 2002, 2005, 2008 and 2009): 21 start dates. All members are initiated on 1 November. Time series are detrended by subtracting a linear trend line to focus on multi-year and inter-annual variability. We note that DePreSys3 precipitation estimates in the Sahel show deficiencies in capturing multi-decadal trends. A posteriori adjustment of trends is therefore required to make real-time predictions59,60. No other pre-processing was performed.

Multi-year and inter-annual variability

The multi-year component of year one lead-time hindcasts is extracted by smoothing time series with a 5-year running mean. The residual is used as the inter-annual signal. The time series are shown in Supplementary Fig. 2. Decomposing the Year 1 hindcasts in this way provides a fair comparison, because both multi-year and inter-annual time series have the same lead time and number of data points (this would not be the case if we compared the inter-annual component of the year 1 lead-time data with the hindcasts averaged 2–5 years from initialization). We do however note that moisture budget analysis performed on hindcasts averaged 2–5 years from initialization, is largely consistent with that from the multi-year component of the Year 1 hindcast data (Supplementary Fig. 10). Resultant timeseries are divided into wet and dry Sahel composites by grouping wetter/dryer than average years and using a weighted average with weights determined by the absolute values of Sahel precipitation. In this way, all data points are included in wet and dry composites but very wet or dry years contribute more.

Observations and re-analysis products

Precipitation data were obtained from the monthly GPCC land-based dataset (1901–present)61. We note that results were similar when using the Climate Research Unit62 precipitation data set and a subset of years with precipitation estimates from the Combined Precipitation Data Set (GPCP, available 1979–2015)63. The HadCRUT4 data set was used for global surface temperature estimates64. For winds, specific humidity and temperatures at all atmospheric levels, we use the NCEP/NCAR Reanalysis 1 project data provided by the NOAA National Center for Environmental Prediction65. Monthly means of daily means of variables were used. We do not include analysis, which requires the use of re-analysis precipitation and evaporation estimates, due to the large to uncertainties associated with these variables (Supplementary Fig. 3)66. For example, we found that the moisture budget in the Sahel did not close in the re-analysis data (ΔP−ΔE≠ΔDIV), stressing the limitations of diagnostic variables relating to the hydrological cycle in re-analysis products and the utility of model data in understanding Sahelian summer rainfall character.

Moisture budget analysis and recycling ratio calculations

Monthly means of specific humidity, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m5.jpg, and wind, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m6.jpg, at different atmospheric levels are output by DePreSys3 (angled brackets represent the monthly mean). To evaluate the moisture flux contribution to the total Sahel region from each side of the domain as a function of atmospheric pressure, we average An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m7.jpg across the domain edge for each pressure level and divide by the domain length along the trajectory of the moisture flux33. In addition, moisture fluxes plotted as a function of pressure in Fig. 4, have been divided by ρwg, where g is the acceleration due to gravity, and ρw is the density of water. The result is to give units of mm per day per hPa. We note that estimates of An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m8.jpg represent the component of the total moisture flux, An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m9.jpg, due the mean as opposed to the transient flow67.

Exact column integrated moisture fluxes are computed during the model integration as An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m10.jpg, where An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m11.jpg represents the monthly mean, g is the acceleration due to gravity, ρw is the density of water and ps is the surface pressure. Output of this integral is in units of m2 s−1, that is, it represents the flux of total moisture depth across each box edge. We average integrated fluxes across each domain edge and scale them by the length along the flux trajectory33, such that units are in mm per day and are directly comparable to domain averaged values of precipitation and evaporation. We denote these scaled column integrated moisture fluxes as An external file that holds a picture, illustration, etc.
Object name is ncomms14966-m12.jpg. The total moisture flux into, Fin, and out of, Fout, the Sahel are evaluated by summing the appropriate combinations of column integrated moisture fluxes, depending on their sign (that is, whether they represent flux advection into or out of the domain).

Recycling ratios, ρ, are evaluated as Pm/P=E/(P+2F), where Pm is the contribution of precipitation due to the recycling of local moisture, P the total precipitation, E the evaporation and F=0.5(Fin+Fout) is the average horizontal moisture flux per unit area due to the total flow. The recycling ratio is derived from moisture conservation using PE=FinFout (ref. 33). Smaller values of ρ indicate a dominance of the advective contribution to precipitate over the rate of recycling of local moisture.

Assessing the significance of correlations

To account for a finite ensemble size, Nens, and finite number of points in each timeseries, Nt, on skill estimates of DePreSys3 ensemble mean outputs, we create an additional 1,000 hindcast time series as follows: Step 1: we randomly sample with replacement Nens ensemble members, from which we create an ensemble mean. Step 2: we resample and shuffle the resulting ensemble mean timeseries using block bootstrapping with a length of 4 years, being careful to account for irregular hindcast start times. The corresponding observational data points are reshuffled in the same way. Step 3: we compute the correlation between the re-shuffled ensemble mean and observational time series. Step 4: steps 1–3 are repeated 1,000 times. Correlation error bars for precipitation and surface temperatures as quoted in the text are computed as the 5–95% confidence interval limits from the resultant probability distribution of 1,000 correlation values. Stippling on maps represent grid points where >95% of the sampled hindcasts skill exceed zero68. For Sahel precipitation correlations shown in Fig. 1, we also evaluate whether the skill is significantly better than that which could be obtained by chance: we check that the skill exceeds the 95th percentile of a correlation distribution computed using the 1,000 hindcast sub-samples (as above) and the original unshuffled observation time series.

Data availability

The data that support the findings of this study are available from the corresponding authors upon request.

Additional information

How to cite this article: Sheen, K. L. et al. Skilful prediction of Sahel summer rainfall on inter-annual and multi-year timescales. Nat. Commun. 8, 14966 doi: 10.1038/ncomms14966 (2017).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:

Supplementary Figures, Supplementary Notes and Supplementary References

Peer Review File:


This work was supported by the Joint DECC/Defra Met Office Hadley Centre Climate Programme (GA01101) and the EU FP7 SPECS project. The contribution of D.P.R. has received funding from the NERC/DFID Future Climate for Africa programme under the AMMA-2050 project, grant number NE/M019977/1. We are also grateful for some discussions with A. Scaife and C. Pomposi. We also thank two reviewers for their helpful suggestions.


The authors declare no competing financial interests.

Author contributions K.L.S. led the data analysis and writing, with suggestions and comments from all authors. K.L.S., D.M.S. and N.J.D. conceived the ideas for the work. N.J.D. set-up and ran the DePreSys3 hindcasts. R.E. analysed the Sahel precipitation skill of other CMIP5 models.


  • Thorncroft C. D., Nguyan H., Zhang C. & Peyrille P. Annual cycle of the West African Monsoon: regional circulations and associated water vapour transport. Q. J. R. Meteorol. Soc. 137, 129–147 (2011).
  • Nicholson S. E. The west African Sahel: a review of recent studies on the rainfall regime and its interannual variability. ISRN Meteorol. 2013, 453521 (2012).
  • Martin E. R. & Thorncroft C. D. The impact of the AMO on the West African monsoon annual cycle. Q. J. R. Meteorol. Soc. 140, 31–46 (2014).
  • Kandji S. T., Verchot L. & Mackensen J. Climate change and variability in the Sahel region: impacts and adaptation strategies in the agricultural sector. World Agroforestry Centre (ICRAF) and United Nations Environment Programme (UNEP) (2006).
  • Rowell D. P., Senior C. A., Vellinga M. & Graham R. Can climate projection uncertainty be constrained over Africa using metrics of contemporary performance? Clim. Change 134, 621–633 (2016).
  • Otero N., Mohino E. & Gaetani M. Decadal prediction of Sahel rainfall using dynamics-based indices. Clim. Dyn. 47, 3415–3431 (2016).
  • Martin E. R. & Thorncroft C. Sahel rainfall in multimodel CIMP5 decadal hindcasts. Geophys. Res. Lett. 41, 2169–2175 (2014).
  • Garcia-Serrano J., Doblas-Reyes F. J., Haarsma R. J. & Polo I. Decadal prediction of the dominant West African monsoon rainfall modes. J. Geophys. Res. 118, 5260–5279 (2013).
  • Gaetani M. & Mohino E. Decadal prediction of the Sahelian Precipitation in CMIP5 simulations. J. Clim. 26, 7708–7719 (2013).
  • Folland C. K., Owen J., Ward M. N. & Colman A. W. Prediction of seasonal rainfall in the Sahel region using empirical and dynamical methods. J. Forecast. 10, 21–56 (1991).
  • Ward M. N. et al. . Experimental seasonal forecasting of tropical rainfall at the Meteorological Office. Predict. Inter. Clim. Variations NATO ASI Ser. 10, 197–216 (1993).
  • Rowell D. P., Folland C. K., Maskell K. & Ward M. N. Variability of summer rainfall over tropical north Africa (1906-92): observations and modelling. Q. J. R. Meteorol. Soc. 121, 669–704 (1995).
  • Fontaine B. & Camberlin P. An improvement of June–September rainfall forecasting in the Sahel based upon region April–May moist static energy content (1968–1997). Geophys. Res. Lett. 26, 2041–2044 (1999).
  • Garlic G., Douville H. & Deque M. Prospects for improved seasonal predictions of monsoon precipitation over Sahel. Int. J. Climatol. 22, 331–345 (2002).
  • Ndiaye O., Goddard L. & Ward M. N. Using regional wind fields to improve general circulation model forecasts of July–September Sahel rainfall. Int. J. Climatol. 29, 1262–1275 (2009).
  • Ndiaye O., Ward M. N. & Thiam W. M. Predictability of seasonal Sahel rainfall using GCMs and lead-time improvements through the use of a coupled model. J. Clim. 24, 1931–1949 (2011).
  • Badr H. S., Zaitchik B. F. & Guikema S. D. Application of statistical models to the prediction of seasonal rainfall anomalies over the Sahel. J. Appl. Meteor. Clim. 53, 614–636 (2013).
  • Mohino E. et al. . Changes in the interannual SST-forced signals on West Africa rainfall. AGCM intercomparison. Clim. Dyn. 37, 1707–1725 (2011).
  • Doblas-Reyes F. J. et al. . Initialized near-term regional climate change prediction. Nat. Commun. 302, 1027–1030 (2003).
  • Rowell D. P. Simulating SST Teleconnections to Africa: what is the state of the art? J. Clim. 26, 5397–5418 (2013).
  • Martin E. R., Thorncroft C. D. & Booth B. B. B. The multidecadal Atlantic SST-Sahel rainfall teleconnection in CMIP5 simulations. J. Clim. 27, 784–806 (2014).
  • Dunstone N. J. et al. . Skilful predictions of the winter North Atlantic Oscillation one year ahead. Nat. Geosci. 9, 809–814 (2016).
  • Smith D. M. et al. . Improved surface prediction for the coming decade from a global climate model. Science 317, 796–799 (2007). [PubMed]
  • Smith D. M. et al. . Skilful multi-year predictions of Atlantic hurricane frequency. Nat. Geosci. 3, 846–849 (2010).
  • Vellinga M. et al. . Sahel decadal rainfall variability and the role of model horizonal resolution. Geophys. Res. Lett. 43, 326–333 (2016).
  • Nicholson S. E. A revised picture of the structure of the ‘monsoon’ and land ITCZ over West Africa. Clim. Dyn. 32, 1155–1171 (2009).
  • Pu B. & Cook K. H. Dynamics of the west African westerly jet. J. Clim. 23, 6263–6276 (2010).
  • Nicholson S. E. The intensity, location and structure of the tropical rainbelt over west Africa as factors in interannual variability. Int. J. Climatol. 28, 1775–1785 (2008).
  • Grist J. E. & Nicholson E. A study of the dynamic factors influencing the rainfall variability in the West African Sahel. J. Clim. 14, 1337–1359 (2001).
  • Pomposi C. P., Giannini A., Kushnir Y. & Lee D. E. Understanding Pacific Ocean influence on interannual precipitation variability in the Sahel. Geophys. Res. Lett. 43, 9234–9242 (2016).
  • Rowell D. P. The impact of Mediterranean SSTs on the Sahelian rainfall season. J. Clim. 16, 849–862 (2003).
  • Cook K. Generation of the African easterly jet and its role in determining west African precipitation. J. Clim. 12, 1165–1184 (1999).
  • Trenberth K. E. Atmospheric moisture recycling: Role of advection and local evaporation. J. Clim. 12, 1368–1381 (1999).
  • Pu B. & Cook K. H. Role of the west African westerly jet in Sahel rainfall variations. J. Clim. 25, 2880–2896 (2012).
  • Folland C. K., Palmer T. N. & Parker D. E. Sahel rainfall and worldwide sea temperatures, 1901-85. Nature 320, 602–607 (1986).
  • Biasutti M., Sobel A. H. & Camargo S. J. The Role of the Sahara low in summertime Sahel rainfall variability and change in the CMIP3 models. J. Clim. 22, 5755–5771 (2009).
  • Park J., Bader J. & Matei D. Anthropogenic Mediterranean warming essential driver for present and future Sahel rainfall. Nat. Clim. Change 6, 941–945 (2016).
  • Giannini A., Saravanan R. & Chang P. Oceanic forcing of Sahel rainfall on interannual to interdecadal time scales. Science 302, 1027–1030 (2003). [PubMed]
  • Knight J. R., Folland C. K. & Scaife A. A. Climate impacts of the Atlantic multidecadal oscillation. Geosphys. Res. Lett. 33, L17706 (2006).
  • Zhang R. & Delworth T. L. Impact of Atlantic mutlidecadal oscillations on Indian/Sahel rainfall and Atlantic hurricanes. Geophys. Res. Lett. 33, L17712 (2006).
  • Rodríguez-Fonseca B. et al. . Variability and predictability of West African droughts: a review on the role of seas surface temperature anomalies. J. Climate 28, 4034–4060 (2015).
  • Rowell D. P. Teleconnections between the tropical Pacific and the Sahel. Q. T. R. Meteorol. Soc. 127, 1683–1706 (2001).
  • Kang S. M., Held I. M., Frierson D. M. W. & Zhao M. The response of the ITCZ to extratropical thermal forcing: Idealised slab-ocean experiments with a GCM. J. Clim. 21, 3521–3532 (2008).
  • Dunstone N. J., Smith D. M. & Eade R. Multiyear predictability of the tropical Atlantic atmosphere driven by the high latitude North Atlantic Ocean. Geophys. Res. Lett. 38, L14701 (2011).
  • Hulme M. & Tosdevin N. The tropical easterly jet and Sudan rainfall: a review. Theor. Appl. Climatol. 39, 179–187 (1988).
  • Chiang J. C. H. & Sobel A. H. Tropical tropospheric temperature variations caused by ENSO and their influence on the remote tropical climate. J. Clim. 15, 2616–2631 (2002).
  • Bader J. & Latif M. The impact of decadal-scale Indian Ocean sea surface temperature anomalies on Sahelian rainfall and the North Atlantic Oscillation. Geophys. Res. Lett. 30, 2169 (2003).
  • Parhi P., Giannini A., Gentine P. & Lall U. Resolving contrasting regional rainfall responses to El Niño over tropical Africa. J. Clim. 29, 1461–1476 (2016).
  • Scaife A. A. et al. . The CLIVAR C20C project: selected twentieth century climate events. Clim. Dyn. 33, 603–614 (2009).
  • Biasutti M. Forced Sahel rainfall trends in the CMIP5 archive. Geophys. Res. Lett. 118, 1613–1623 (2013).
  • Eade R. et al. . Do seasonal-to-decadal climate predictions underestimate the predictability of the real world? Geophys. Res. Lett. 41, 5620–5628 (2014). [PMC free article] [PubMed]
  • Scaife A. A. et al. . Skillful long-range prediction of European and North American winters. Geophys. Res. Lett. 41, 2514–2519 (2014).
  • Giannini A. et al. . A unifying view of climate change in the Sahel linking intra-seasonal, interannual and longer time scales. Environ. Res. Lett 8, 024010 (2013).
  • Willams K. et al. . The Met Office Global Coupled Model 2.9 (GC2) configuration. J. Geosci. Model Dev. 88, 1509–1524 (2015).
  • Smith D. M. & Murphy J. M. An objective temperature and salinity analysis using covariances from a global climate model. J. Geophys. Res. 112, C02022 (2007).
  • Rayner N. A. et al. . Global analyses of sea surface temperature, sea ice and night marine air temperature since the late nineteenth century. J. Geophys. Res. 108, 4407 (2003).
  • Dee D. P. The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Quart. J. Roy. Meteor. Soc. 137, 553–607 (2011).
  • Taylor K. E., Stouffer R. J. & Meehl G. A. An overview of CMIP5 and the experimental design. Bull. Am. Meteorol. Soc. 93, 485–498 (2012).
  • Kharin V. V., Boer G. J., Merryfield W. J., Scinocca J. F. & Lee W. S. Statistical adjustment of decadal predictions in a changing climate. Geophys. Res. Lett. 39, L19705 (2012).
  • Fučkar N. S., Volpi D., Guemas V. & Doblas-Reyes F. J. A posteriori adjustment of near-term climate predictions: accounting for the drift dependence on the initial conditions. Geophys.Res. Lett. 41, 5200–5207 (2014).
  • Schneider U. et al. . GPCC’s new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle. Theor. Appl. Climatol. 115, 15–40 (2014).
  • Harris I., Jones P. D., Osborn T. J. & Lister D. H. Updated high-resolution grids of monthly climatic observations the CRU TS3.10 Dataset. Int. J. Climatol. 34, 623–642 (2014).
  • Xie P. et al. . GPCP Pentad precipitation analyses: an experimental dataset based on gauge observations and satellite estimates. J. Clim 16, 2197–2214 (2003).
  • Morice C. P., Kennedy J. J., Rayner N. A. & Jones P. D. Quantifying uncertainties in global and regional temperature change using an ensemble of observational estimates: the HadCRUT4 data set. J. Geophys. Res. 117, D08101 (2012).
  • Kalnay et al. . The NCEP/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77, 437–470 (1996).
  • Seager R. & Henderson N. Diagnostic computation of moisture budgets in the ERA-Interim reanalysis with reference to analysis of CMIP-archived atmospheric model data. J. Clim. 26, 7876–7901 (2013).
  • Pomposi C. P., Kushnir Y. & Giannini A. Moisture budget analysis of SST-driven decadal Sahel precipitation variability in the twentieth century. Clim. Dyn. 44, 3303–3321 (2014).
  • Knight J. R. et al. . Predictions of climate several years ahead using an improved decadal prediction system. J. Clim. 27, 7550–7567 (2014).
  • Trenberth K. E. The Definition of El Niño. Bull. Am. Meteorol. Soc. 78, 2771–2777 (1997).

Articles from Nature Communications are provided here courtesy of Nature Publishing Group