This article analyzes SST remote forcing on the interannual variability of Sahel summer (June–September) moderate (below 75th percentile) and heavy (above 75th percentile) daily precipitation events during the period 1981–2016. Evidence is given that interannual variability of these events is markedly different. The occurrence of moderate daily rainfall events appears to be enhanced by positive SST anomalies over the tropical North Atlantic and Mediterranean, which act to increase low-level moisture advection toward the Sahel from the equatorial and north tropical Atlantic (the opposite holds for negative SSTs anomalies). In contrast, heavy and extreme daily rainfall events seem to be linked to El Niño–Southern Oscillation (ENSO) and Mediterranean variability. Under La Niña conditions and a warmer Mediterranean, vertical atmospheric instability is increased over the Sahel and low-level moisture supply from the equatorial Atlantic is enhanced over the area (the reverse is found for opposite-sign SST anomalies). Further evidence suggests that interannual variability of Sahel rainfall is mainly dominated by the extreme events. These results have implications for seasonal forecasting of Sahel moderate and heavy precipitation events based on SST predictors, as significant predictability is found from 1 to 4 months in advance.
The Sahel is facing an increase of extreme rainfall occurrence (Ly et al. 2013). Recently, Taylor et al. (2017) have shown that extreme thunderstorms over the Sahel are nowadays 3 times more frequent than they were 35 years ago. These extremes constitute the primary impact of climate change on society (Katz and Brown 1992), more than changes in the mean climate (Mitchell et al. 1990). Due to the vulnerability of the Sahel region, precipitation extremes often directly (and indirectly) impact numerous economy sectors and are associated with elevated death tolls (New et al. 2006; Lobell et al. 2011; Anyamba et al. 2014; Sané et al. 2015, 2016). The predictability of these extreme events is a big challenge for the climate science community in the Sahel area.
Among processes controlling Sahel rainfall variability, oceanic forcing has been found to be the dominant driver (Giannini et al. 2003; Rowell 2013; Rodríguez-Fonseca et al. 2011, 2015; Gómara et al. 2018). Many works have shown that seasonal rainfall variability over the Sahel is related to sea surface temperature (SST) anomalies over the tropical Atlantic, Pacific, and the Mediterranean Sea (Rodríguez-Fonseca et al. 2015, and references therein). In particular, a warming of Pacific and Atlantic tropical waters has been found to promote a decrease in seasonal Sahelian rainfall (e.g., Rowell 2001; Janicot et al. 2001; Mohino et al. 2011a,b; Losada et al. 2010), while a warmer Mediterranean leads to an increase (Rowell 2013; Fontaine et al. 2010). These relations are somewhat linear, as opposite responses in rainfall take place when colder-than-usual SSTs arise over the same marine regions.
Regarding Sahel daily rainfall intensity/frequency, Giannini et al. (2013) found that the recovery during the last two decades was accompanied by the relative predominance of an increase in median rainfall intensity. Such an increase was shown to coincide with positive SST anomalies over the tropical North Atlantic (TNA), potentially associated with a positive phase of the Atlantic multidecadal oscillation (AMO). Additionally, it was also shown that extreme dry spell occurrence over West Africa is significantly linked to anomalous SSTs in the Indian and South Atlantic oceans (Bader and Latif 2011; Salack et al. 2014). Parhi et al. (2016) also revealed that El Niño–Southern Oscillation (ENSO) significantly influences the interannual variability of Sahel rainy days frequency during boreal summer.
In contrast, the possible influence of SST forcing on Sahel precipitation, stratified and isolated in terms of event intensity (e.g., moderate, extreme), has received little attention in the literature so far. Ta et al. (2016) recently endeavored this subject, but their analysis focused on tropical Atlantic SST forcing over the Sahelian and Guinean regions, and they looked for oceanic forcing associated with trends on extreme rainfall. Results from this study pointed out the existence of a significant link between extreme rainfall trends over West Africa and SST anomalies in the equatorial Atlantic.
Distinctively, the present paper focuses on the interannual variability of Sahel rainfall events of different intensities (e.g., low, moderate, heavy, and extreme), searching for oceanic sources of predictability worldwide to provide skillful seasonal forecasts. The article is structured as follows. Section 2 describes the data and methodology used. Section 3 is centered in the analysis of Sahel moderate and heavy precipitation indices, their associated teleconnections and mechanisms, and their seasonal predictability. A summary of the main findings concludes this article.
2. Data and methodology
A daily rainfall dataset, the Climate Hazards Group Infrared Precipitation with Station Data (CHIRPS), is utilized. The CHIRPS dataset spans from 1981 to present day at 0.25° × 0.25° longitude–latitude horizontal resolution and is based upon satellite infrared and rain gauge measurements over West Africa (Funk et al. 2014, 2015). Daily June–September (JJAS) precipitation measurements are considered to characterize Sahel precipitation indices over the area (20°W–10°E, 10°–20°N). Note that similar results are obtained for daily July–September (JAS) measurements. Depending on percentile thresholds, indices are stratified in terms of daily rainfall intensity: low (<50th), moderate (<75th), heavy (>75th), and extreme (>95th) events. Therefore, the moderate category also includes low precipitation events. Likewise, the heavy category contains extreme daily events. Thresholds are calculated over each Sahel grid point separately, with rainy days defined as days in which the rainfall amount is greater than 1 mm. Subsequently, for each precipitation category, the gridpoint precipitation values fulfilling the intensity criterion are retained and averaged over the total Sahel area for each JJAS season. Then, the seasonal anomalies of these indices are calculated by removing the seasonal mean. The obtained time series are standardized and high-pass filtered (Lanczos; 8-yr cutoff period) to keep only with interannual variability (one data point per year). The resultant time series are considered as our rainfall indices along the whole study.
The SST dataset from the Met Office Hadley Centre Sea Ice and Sea Surface Temperature dataset (HadISST) is used and regressed on precipitation indices to search for possible oceanic teleconnections. The SST data have a horizontal grid of 1° × 1° in longitude–latitude (Rayner et al. 2003). The indices of ENSO, TNA, and the Mediterranean Sea Index (MSI) are computed based on high-pass-filtered SST anomalies, averaged over the regions 180°–80°W, 5°S–5°N; 60°–20°W, 5°–30°N; and 5°W–30°E, 30°–45°N, respectively.
To analyze the underlying physics of the identified teleconnection mechanisms, atmospheric data from the European Center for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis (period 1979–2017) are also utilized. ERA-Interim is provided at 0.75° × 0.75° horizontal resolution with 60 vertical atmospheric levels, from the surface up to 0.1 hPa (Dee et al. 2011). The dynamical fields analyzed to compute JJAS seasonal anomalies are (i) specific humidity at 850 hPa; (ii) moisture flux at 850 hPa; and (iii) vertical instability, which is the difference of wind divergence between the 200- and 850-hPa pressure levels (hereafter DIV200/850; Gómara et al. 2018). Positive DIV200/850 values thus indicate vertical destabilization and upward motion of air (the reverse is true for negative values). The sea level pressure (SLP) is also used (instead of the DIV200/850 parameter) for the MSI, as the influence of the Mediterranean Sea in Sahelian rainfall has been found to take place through anomalous moisture advection at surface levels, mediated by the SLP gradient between the Gulf of Guinea and the Sahara (Rowell 2013; Fontaine et al. 2010; Gómara et al. 2018).
To identify covariability between fields, regression maps are calculated by projecting a time-evolving field onto the time series associated with a particular index. For each grid point (i, j), a linear regression is done between an index I, which is standardized (e.g., rainfall or TNA), and the value of anomalous of the variable Y at that grid point:
To further analyze the atmospheric stability associated with moderate/weak and heavy/extreme precipitation events over the Sahel, moist static energy (MSE) is calculated:
where T is the temperature, CP is the specific heat at constant pressure, z is the height, g is the gravitational acceleration, Lυ is the latent heat of vaporization at 0°C, and q is the specific humidity. MSE thus combines the contributions from latent heat (LH), sensible heat (SH), and potential energy (PE) terms. The MSE is a useful tool to investigate the relative roles of moisture and temperature changes to this stability (Sheen et al. 2017; Hill et al. 2017). An increasing MSE with altitude denotes a stable atmosphere.
Due to the time limitations posed by the observational datasets, the time interval considered for the study is 1981–2016, a period in which a strong Sahel–ENSO teleconnection has been documented (Janicot et al. 2001; Mohino et al. 2011a; Rodríguez-Fonseca et al. 2011, 2015; Gómara et al. 2018). To analyze the predictability of Sahel daily precipitation events based on SST anomalies, a stepwise multilinear regression statistical model is utilized for the same period (cf. section 3d). Finally, a two-tailed t test that accounts for autocorrelation of the time series (Bretherton et al. 1999) is used throughout the paper for testing statistical significance of the results (95% confidence interval).
a. Interannual variability of Sahel moderate and heavy daily rainfall
Total Sahelian JJAS rainfall variability is separated into the contributions of low, moderate, heavy, and extreme daily rainfall events. A first result arises when performing this separation, as moderate and heavy daily rainfall variability indices show very different interannual time evolutions, appearing to be caused by distinctive physical mechanisms. Figure 1 shows that, except for some periods and years like 1993–96, 1998–2002, 2014, and 2015, the time evolution of these indices is decoupled. This is especially evident before 1993 and after 2005.
Correlation values between low, moderate, heavy, and extreme indices are presented in Table 1. Quantitatively, correlation values between moderate and heavy indices are found to be positive (0.24), but not statistically significant, therefore supporting the conclusion that variability of moderate and heavy rainfall events over the Sahel is decoupled.
To quantify the imprint of each of the rainfall intensity categories on the total rainfall index, their correlation values are provided in Table 2, together with their standard deviations. The correlation value between moderate and total precipitation indices is 0.46 and that between heavy and total precipitation is 0.97. Consequently, interannual variability of total Sahel rainfall is found to be mainly dominated by heavy daily precipitation events. However, regarding the contribution to yearly precipitation amount, both heavy and moderate rainfall are important (51% and 49% of the total annual precipitation, respectively). For the remaining categories, the contribution of low daily rainfall events to the total precipitation is found to be more important in terms of yearly rainfall accumulation (23%) than in terms of explained interannual variability (cf. Table 2, right column).
b. SST teleconnection patterns
Consistent with previous studies (e.g., Rowell 2001; Janicot et al. 2001; Mohino et al. 2011b; Fontaine et al. 2010; Rodríguez-Fonseca et al. 2011), an increase of JJAS Sahel total rainfall is associated with a La Niña–like pattern in the equatorial Pacific and warmer-than-usual SSTs over the Mediterranean and subtropical North Atlantic (Fig. 2a). By linear construction, the opposite is also true for a decrease in total seasonal rainfall. Figures 2b–d show JJAS SST forcing patterns associated with moderate, heavy, and extreme rainfall events, respectively. In this context, please note that similar SST patterns have been found considering slight variations in the total study period (1981–2011, 1981–2012, or 1981–2013). The enhanced occurrence of moderate daily Sahelian rainfall events is related to an SST warming of the tropical/subtropical North Atlantic (Fig. 2b), while a La Niña–like signal appears over the equatorial Pacific, associated with heavy and extreme daily rainfall events (Figs. 2c,d). Over the North Atlantic a horseshoe pattern, typically associated with the negative phase of the North Atlantic Oscillation (NAO), is also observed in Fig. 2b. Indeed, under a negative NAO phase, the weakening of the trade winds decreases the latent heat flux and sea level pressure over subtropical latitudes, warming the oceanic waters beneath through the wind–evaporation–SST (WES) feedback process (Okumura et al. 2001; Czaja and Frankignoul 1999). The Mediterranean (MSI) also shows a significant link with both moderate and heavy/extreme events (correlation values of 0.67, 0.62, and 0.31 for moderate, heavy, and extreme events, respectively; cf. Table 3). The patterns associated with heavy (Fig. 2c) and extreme daily precipitation events (Fig. 2d) mostly resemble those associated with total rainfall (Fig. 2a). Nevertheless, the intensity/extent of the SST anomalies over the Pacific appears enhanced for the case of heavy/extreme daily precipitation events, pointing out to a tropical Pacific forcing as a main driver of intense rainfall, compared to the secondary role of extratropical patterns (North Atlantic and Mediterranean; Figs. 2c,d).
c. Mechanisms of SST teleconnections on Sahel moderate and heavy daily rainfall
In Fig. 3 the underlying physical mechanisms of the teleconnections identified in Fig. 2 are analyzed. Figure 3a shows that Sahel total rainfall variability is associated with anomalous southwesterly low-level moisture flux, which transports humidity from the north tropical and equatorial Atlantic toward the Sahel (arrows). The moisture advection enhances the available low-level humidity over this area (shadings). Positive vertical instability (DIV200/850) anomalies are also observed over the same region for total rainfall, reflecting low-level wind convergence and upper-level divergence associated with vertical air movement and, potentially, convection (contours).
A similar picture is observed for heavy and extreme daily precipitation events (Figs. 3c,d). Organized deep convective systems are characterized by a baroclinic structure with opposite divergent circulation patterns at upper levels, a feature denoted in contours. The maximum anomalies of low-level humidity are located over the central and eastern Sahel (Figs. 3c,d). This flow is driven by the large-scale circulation and probably associated with La Niña, which acts to promote vertical air destabilization over the Sahel (Joly and Voldoire 2009).
Conversely, for moderate daily rainfall events (Fig. 3b), local negative DIV200/850 anomalies are present over the Sahel region, suggesting that local convection is not enhanced for these events. There is, however, enhanced low-level humidity stemming from the TNA region, which is consistent with a warming in the tropical/subtropical North Atlantic waters (Fig. 2b). Such warming acts to increase evaporation over these areas and therefore moisture advection toward the Sahel.
These results indicate that heavy and extreme daily rainfall events seem to be mainly associated with ENSO, while moderate events seem to be linked to SST anomalies from the TNA. SST variability over the Mediterranean appears to be important for both moderate and heavy/extreme rainfall (cf. Figs. 2b–d). However, regarding the mechanism, Figs. 3b and 3c show no clear significant low-level humidity flux from the Mediterranean toward the Sahel.
To test the robustness of the main dynamical mechanisms proposed, a reverse methodological approach is carried out in Fig. 4. First, worldwide SST anomalies are regressed onto ENSO, TNA, and MSI indices (Figs. 4a–c). Second, atmospheric variables from Fig. 3 are also regressed on TNA and ENSO indices (Figs. 4d,e). The same analysis is done for the MSI index by just replacing the DIV200/850 parameter by SLP, as explained before (Fig. 4f).
For a TNA positive phase (as for moderate rainfall events), the regression map shows enhanced low-level specific humidity over the Sahel region, transported from the tropical/subtropical North Atlantic (Fig. 4d). In this area, a horseshoe pattern (Fig. 4a), together with a negative phase of the NAO, may act to reduce the trade winds, thus increasing humidity advection over the Sahel (from a warming of the TNA).
For an ENSO negative phase (as for Sahel heavy/extreme events), results show enhanced northeastward low-level moisture flux, transporting humidity from the northern tropical and equatorial Atlantic toward the Sahel. Such anomalous flow is again consistent with the vertical destabilization (stabilization) that La Niña (El Niño) induces over the Sahel, promoting wind convergence (divergence) at 850 hPa and divergence (convergence) at 200 hPa (contours in Fig. 4e). Figure 4b also reveals that an equatorial Pacific cooling (warming) is associated with a concomitant equatorial Atlantic warming (cooling). The simultaneous appearance of these anomalies for the period considered (1981–2016) has been reported in previous studies (Rodríguez-Fonseca et al. 2009).
Finally, for MSI, Fig. 4f shows that a warmer Mediterranean enhances evaporation and strengthens the Sahara heat low, thus increasing the pressure gradient between the Sahara and Gulf of Guinea in agreement with previous studies (Rowell 2013; Peyrillé et al. 2007; Fontaine et al. 2010; Gaetani et al. 2010; Gómara et al. 2018). The increased meridional pressure gradient promotes a further penetration of the monsoonal flow from the Gulf of Guinea (anomalous southwesterly positive moisture flux). Figure 4f also shows that this meridional circulation strengthens with changes in the heat low, combined with a large change of westerly flow from the tropical North Atlantic, which also increases moisture and heat transport into the Sahel (Martin and Thorncroft 2014). In addition, as it will be shown later, the inclusion of the MSI leads to enhanced predictive skill (cf. section 3d).
To further analyze the origin of rainfall variability and its relation to intensity, an analysis of MSE [Eq. (2)] is subsequently performed in Fig. 5. For the variability of moderate and extreme daily precipitation events, the contribution of the PE component (Fig. 5d) is negligible compared to the LH and SH terms (Figs. 5b,c). Attending to moderate and extreme daily precipitation vertical profiles, the stronger differences arise in the low tropospheric levels (Figs. 5a,b).
MSE associated with low to moderate daily precipitation decreases strongly with altitude compared to that associated with extremes (Fig. 5a), which suggests instability in the low levels (Sheen et al. 2017). Figure 5c shows that this strong decrease is mainly explained by the LH component. However, the LH component decreases with a similar rate for extremes too. Surface evaporation is important over land where soil moisture is available, a feature that is very important in moderate rainfall events, where MSE is controlled by LH. In contrast, during extreme daily precipitation events, the contribution of the LH component (Fig. 5c) is balanced at low tropospheric levels by strong negative values of SH (Fig. 5b). These strong negative SH values could depict the equilibrium response once convection has overturned the vertically integrated temperature column (Wei et al. 2014). At medium and upper levels (500–200 hPa), the MSE profile is more unstable for extreme than for moderate precipitation events (Fig. 5a). This can be related to both LH and SH components. On the one hand, the upper troposphere (400–200 hPa) is colder in the extreme case, consistent with a tropical cooling due to La Niña (Lintner and Chiang 2007), which affects the SH component. On the other hand, the middle troposphere is more humid in the extreme case, which would be related to a stronger vertical rise of moisture-rich air (Fig. 3d).
In summary, moderate events seem to be controlled by instability in lower levels associated with changes in latent heat fluxes due to surface evaporation. The influence of the subtropical North Atlantic and Mediterranean SSTs is felt through surface processes and moisture flux advection from the warm seas. In contrast, extreme rainfall events are controlled by upper-level instability and convective precipitation systems, a mechanism consistent with ENSO remote influence, as the teleconnection with ENSO is related to changes in the upper-level divergence, favoring (La Niña) or inhibiting (El Niño) the instability in the upper levels.
This hypothesis is confirmed in Fig. 6, where air temperature anomalies from the surface up to 200 hPa are regressed on moderate and extreme precipitation indices and averaged over the corridor 0°–17°N. Similar results (but with a weaker signal) are obtained when limiting the latitudinal band to that of Sahel (10°–20°N). The extended band has therefore been chosen in order capture a better ENSO signature, as the teleconnection associated with the latter is stronger near the equator.
During moderate daily rainfall events, as shown by the SH anomalies profile (Fig. 5b), atmospheric temperature anomalies are significantly warmer over the tropical North Atlantic (Fig. 6a). These results are consistent with the presence of warmer waters over the near subtropical North Atlantic, which heat up the air immediately above. During extreme daily rainfall events, in addition to weak upper-level cooling over the Sahel, associated with La Niña influence (the opposite takes place for El Niño), the temperature anomalies exhibit a strong and significant low-level cooling in the lower levels (Fig. 6b). This cooling could be associated with the temperature equilibrium response, once convection has overturned the column (Fig. 5b). Consequently, MSE is quasi constant with altitude for extreme daily precipitation at lower levels, due to the balance of SH and LH terms. In this context, it is well known that the evaporation (sublimation) of rain (ice) falling droplets produces a strong cooling in convective systems, which subsequently strengthens the downward current itself and cold pools. The latter are considered as fundamental ingredients for deep convection (Torri et al. 2015). However, due to their mesoscale characteristics, their direct causality on seasonal and spatial averaged rainfall indices from this study is extremely hard to detect (Wei et al. 2014).
The results from this section are entirely consistent with the same analysis applied to nonfiltered data (linear trends removed), as the SST patterns considered are mainly interannual (ENSO, TNA, Mediterranean).
d. Predictability of Sahel moderate and heavy/extreme daily rainfall at interannual time scales
As described in the previous section, the key SST indices associated with JJAS Sahel rainfall variability at interannual time scales are the TNA, ENSO, and MSI.
To assess the potential predictability of rainfall events over the Sahel, different predictors are selected based on these SST indices. Potential predictors for each Sahel rainfall index are identified by applying a stepwise multilinear regression. This method permits us to automatically identify the relevant SST predictors for each rainfall index and disregard those which do not increase model skill (forward selection and backward elimination; von Storch and Zwiers 1999). Subsequently, the selected predictors for each index are utilized to perform a leave-one-out cross-validated hindcast at lag 0 (predictors/predictands based on the same monthly period of August). In this way, the linear regression model is built for each year to be predicted, calculating the coefficients of the model with all the years in our database except the one that is predicted (ter Braak and Juggins 1993; Birks 1981). The hindcast is performed and correlated with the omitted observations. The same analysis is extended by utilizing the same SST predictors from the previous months to assess predictability of Sahel moderate and heavy daily rainfall events. Unlike in the previous analyses, all considered indices and anomalies in this section are calculated on a monthly basis. As a result, all selected rainfall predictands (precipitation indices) are computed for the month of August. SST predictors considered for lag 0 are thus simultaneous SST anomalies. Lag-1 forecasts take into account SST anomalies from the previous month of the same year (July), and so on back to one year in advance.
In Table 3, the main SST predictors utilized to predict rainfall indices at lag 0 are provided. As expected, the regression model considers the TNA as a main predictor of moderate Sahel rainfall, together with the MSI. For heavy/extreme precipitation, both ENSO and MSI are considered as skillful predictors. This outcome is fully consistent with the findings of Figs. 2–6.
Based on these predictors, the cross-validation hindcast for each rainfall index at lag 0 is provided (Fig. 7). The model is able to reproduce the temporal evolution of Sahel rainfall indices with high significant skill (0.73, 0.69, and 0.50 correlation coefficients for moderate, heavy, and extreme rainfall indices, respectively; 95% confidence interval).
In Fig. 8 the same analysis is extended to consider SST predictors from lag 0 (August) to −12 (August of year −1). In the figure, lagged correlations between rainfall predictands and SST predictors are provided with bars, whereas skill scores from the multiregression model are provided with markers. For moderate rainfall events, significant forecast skill is obtained from MSI and TNA predictors 2 months ahead (Fig. 8a; 95% confidence level). In other words, the model is able to predict August moderate rainfall based on MSI and TNA SST predictors from June. This outcome is also consistent with lagged correlation values, which are statistically significant (same confidence interval) for both MSI (thin bars) and TNA (wide bars) back to lag −2.
For heavy and extreme rainfall indices the forecast window (based on MSI and ENSO SST predictors) is even increased to 3 and 4 months in advance (May/April), respectively (Figs. 8b,c). Also, the lagged correlations with the MSI index decrease more rapidly backward in time than for ENSO. Correlation values between ENSO and the rainfall indices remain significant from April on (lag −4), whereas those of MSI just reach the previous July (lag −1). The fact that the Mediterranean variability is stronger during the summer months may be behind this behavior (Gaetani et al. 2010). Finally, no correlation is found between August heavy/extreme precipitation anomalies and ENSO during the previous winter (December/January), suggesting that the observed teleconnection takes place during the onset phase of La Niña/El Niño.
4. Conclusions and discussion
The contributions, variability, and teleconnections of different types of precipitation, stratified in terms of intensity, have been analyzed over the Sahel. For this purpose, a precipitation dataset (CHIRPS), a sea surface temperature database (HadISST), and a reanalysis product (ERA-Interim) have been analyzed during the summer boreal season (June–September) over 36 consecutive years (1981–2016). Rainfall categories have been built upon gridpoint percentile thresholds during rainy days (days with precipitation amount greater than 1 mm) as follows: low (<50th), moderate (<75th), heavy (>75th), and extreme (>95th) daily rainfall events. For each category, precipitation values have been retained and averaged over the area (20°W–10°E, 10°–20°N). Last, seasonal precipitation anomalies have been calculated, and a high-pass filter has been applied to the time series to keep only the interannual variability (allowed oscillation periods from 1 to 8 years).
Evidence has been presented that total Sahelian rainfall variability is dominated by the occurrence/absence of heavy and extreme daily rainfall events, rather than moderate or low rainfall episodes, whose variability and teleconnections are markedly different (Tables 1, 2; Fig. 1).
On the one hand, Sahel heavy and extreme precipitation variability is mainly linked to ENSO and, to a lower extent, Mediterranean SST anomalies (MSI). During La Niña episodes (Figs. 2c,d), the induced anomalous circulation over the Sahel enhances vertical destabilization, promoting air convergence at lower levels and divergence at upper levels (Figs. 3c,d and 5). Such conditions, together with the penetration of moisture-rich air from the Gulf of Guinea toward the Sahel (Figs. 3c,d), create a very favorable environment for convection and thus heavy and extreme rainfall occurrence (Figs. 5, 6b). The reverse holds during El Niño and negative MSI episodes.
On the other hand, moderate rainfall variability is primarily associated with TNA and also Mediterranean SST anomalies (Fig. 2b). When these regions’ SSTs are warmer than average, enhanced moisture transport from the equatorial and north tropical Atlantic toward the Sahel is triggered (Sheen et al. 2017), producing an increase of humidity availability over the region (Fig. 3b). In this case, the vertical atmospheric stability is reduced over the Sahel (Fig. 3b) and environmental conditions may be favorable for the formation of horizontal stratiform clouds (Figs. 5, 6a), potentially responsible for the moderate precipitation events. As before, these conditions are reversed for negative TNA and Mediterranean SST anomalies.
The Atlantic influence is associated with enhanced available humidity in the lower and medium levels of the troposphere. The main difference with the Pacific forcing is not the moisture advection, but the vertical destabilization of the atmosphere. For the latter, air uplift is promoted and, therefore, so are heavy precipitation events from deep convection, including organized mesoscale or synoptic systems such as mesoscale convective systems or African easterly waves (AEWs). According to Okonkwo (2014), a La Niña event over the Pacific increases the frequency of AEWs. One can therefore expect a stronger link between extreme rainfall events and AEWs. These AEWs are a major source of synoptic-scale rainfall variability throughout West Africa (Skinner and Diffenbaugh 2013; Martin and Thorncroft 2015). Following the approach from Diedhiou et al. (1998) for the identification of AEWs, in our case the number of AEWs associated with heavy rainfall (22%) is 2 times higher than that associated with moderate ones (11%; cf. Table 4).
In terms of yearly precipitation amounts, moderate and heavy daily rainfall contributions are of the same order of magnitude (~50%; Table 2). This is an important result when tackling, for instance, climate services and impacts in agriculture, as heavy daily precipitation events have more impact than extreme ones and contribute more on interannual variability (cf. Table 2, first column). Some signals over the Mediterranean have been found in relation to heavy daily rainfall events (Fig. 2c), thus suggesting that these results might complement those relating recent trends in Mediterranean warming on Sahelian recovery (Park et al. 2016; Biasutti 2016).
Finally, a stepwise multilinear regression model has been applied to identify the most relevant SST predictors of each rainfall index at different forecast windows. For this purpose, indices of moderate, heavy, and extreme daily precipitation have been computed based on August monthly anomalies (predictands). Monthly SST anomaly indices (ENSO, TNA, and MSI) from the same month (lag 0), and the previous ones (back to lag −12), have been considered as predictors. Results indicate that moderate, heavy, and extreme Sahel rainfall anomalies can be predicted with high confidence utilizing SST predictors several months in advance. TNA and MSI predictors from 2 months in advance, together, are able to provide model skill for the prediction of moderate August Sahelian rainfall. Furthermore, Sahel heavy and extreme precipitation anomalies can be predicted up to 3 and 4 months ahead, respectively, considering ENSO and MSI as predictors. The findings provided in this article may have implications on seasonal forecasting of Sahel moderate, heavy, and extreme rainfall daily events. And in this sense, it must be noted that, due to the relative short time period considered (36 years), SST–rainfall links obtained from this study are assumed to remain stable in time. To analyze the impact of nonstationarity in these relationships (Janicot et al. 2001; Mohino et al. 2011a; Rodríguez-Fonseca et al. 2011, 2015; Suárez-Moreno et al. 2018), a longer-period rainfall dataset should be required. An analogous analysis based on regional model data (e.g., CORDEX-Africa), which considers potential nonlinearities in the SST–precipitation relationships, could also be another interesting path for future research.
The research leading to these results has received funding from the CSIC ICOOP project, ICOOPB20204, the UCM Cooperation Project 7 from XIV Call; the NERC/DFID Future Climate for Africa programme under the AMMA-2050 project, Grant NE/M020428/1; the Spanish Project CGL2017-86415-R; and the EU/FP7 PREFACE project under Grant Agreement 603521. Iñigo Gómara was supported by the Spanish Ministry of Economy and Competitiveness (Juan de la Cierva-Formación contract FJCI-2015-23874) and Universidad Politécnica de Madrid (Programa Propio - Retención de Talento Doctor). The Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) are available from the Climate Hazards Group (http://chg.geog.ucsb.edu/data/chirps/). The Hadley Centre Sea Ice and Sea Surface Temperature dataset (HadISST) is available from the Met Office (https://www.metoffice.gov.uk/hadobs/hadisst/). The European Center for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis (period 1979–2017) is available from the ECMWF MARS server (https://www.ecmwf.int/en/research/climate-reanalysis/era-interim). We would also like to acknowledge Julián Villamayor and Chris Taylor for constructive remarks. Finally, we thank the three anonymous reviewers for their helpful comments and suggestions, which have contributed to improve this manuscript.