Future evolution of the Sahel precipitation zonal contrast in CESM1

The main focus of this study is the zonal contrast of the Sahel precipitation shown in the CMIP5 climate projections: precipitation decreases over the western Sahel (i.e., Senegal and western Mali) and increases over the central Sahel (i.e., eastern Mali, Burkina Faso and Niger). This zonal contrast in future precipitation change is a robust model response to climate change but suffers from a lack of an explanation. To this aim, we study the impact of current and future climate change on Sahel precipitation by using the Large Ensemble of the Community Earth System Model version 1 (CESM1). In CESM1, global warming leads to a strengthening of the zonal contrast, as shown by the difference between the 2060–2099 period (under a high emission scenario) and the 1960–1999 period (under the historical forcing). The zonal contrast is associated with dynamic shifts in the atmospheric circulation. We show that, in absence of a forced response, that is, when only accounting for internal climate variability, the zonal contrast is associated with the Pacific and the tropical Atlantic oceans variability. However, future patterns in sea surface temperature (SST) anomalies are not necessary to explaining the projected strengthening of the zonal contrast. The mechanisms underlying the simulated changes are elucidated by analysing a set of CMIP5 idealised simulations. We show the increase in precipitation over the central Sahel to be mostly associated with the surface warming over northern Africa, which favour the displacement of the monsoon cell northwards. Over the western Sahel, the decrease in Sahel precipitation is associated with a southward shift of the monsoon circulation, and is mostly due to the warming of the SST. These two mechanisms allow explaining the zonal contrast in precipitation change.


Introduction
The third and fifth phases of the Climate Model Intercomparison Projects (CMIP; see Taylor et al. 2012) have allowed the analysis of a large variety of climate simulations. Results have highlighted that the impact of climate change on Sahel precipitation, as a whole coherent region, is strongly modeldependent, ranging from a decrease to an increase in future precipitation (Druyan 2011;Monerie et al. 2016;Yan et al. 2018), hence leading to a rather low signal-to-noise ratio (Caminade and Terray 2010). Therefore, no robust conclusions were reached from climate projections in the fourth and fifth assessment reports (Solomon et al. 2007; Barros et al. 2014).
Yet, numerous studies have been devoted to understand and quantify the impact of climate change on Sahel precipitation. Several mechanisms have then been pointed out. First, it is commonly accepted that the increase in global mean surface temperature leads to an increase in specific humidity following the Clausius-Clapeyron relation, allowing an increase in moisture convergence over the Sahel (i.e., thermodynamic changes) (Held and Soden 2006;Giannini 2010;Kitoh et al. 2013). Second, the warming of the northern Hemisphere is associated with a northward displacement of the monsoon cell, and with an increase in Sahel precipitation (Haarsma et al. 2005;Skinner et al. 2012;Chadwick 1 3 et al. 2019). Besides, the increase in greenhouse gas concentration and temperature is associated with an increase in evaporation (Biasutti 2013) and with changes in moist static energy gradients, which are key factors for the development of the monsoon circulation (Bordoni and Schneider 2008). Third, changes in Sahel precipitation have been associated with changes in sea surface temperature (SST). The warming of the global ocean heats the troposphere and imposes stability, which may decrease moisture transport and convection over land (Held et al. 2005;Caminade and Terray 2010;Giannini 2010;Biasutti 2013;Gaetani et al. 2017;Hill et al. 2018). Patterns in SST change, such as the cross-equatorial Atlantic temperature gradient (Hoerling et al. 2006) and the northern-hemispheric differential warming (Park et al. 2015) have been shown to be one of the causes of Sahel precipitation change and model uncertainties. Therefore, the effects of climate change on Sahel precipitation are the result of complex interactions. Fontaine et al. (2011b) have shown that a robust pattern emerges over the Sahel in response to global warming, with the majority of CMIP3 climate projections showing an increase in central Sahel precipitation and a decrease in western Sahel precipitation. This has also been shown within the CMIP5 ensemble (Monerie et al. 2012(Monerie et al. , 2013(Monerie et al. , 2016Biasutti 2013;Roehrig et al. 2013;James et al. 2015;Diallo et al. 2016;Gaetani et al. 2017;Akinsanola and Zhou 2018;Dunning et al. 2018) and the CMIP6 ensemble (Almazroui et al. 2020;Monerie et al. 2020). The multi-model ensemble also projects a modulation of the seasonal cycle of Sahel precipitation, with a decrease over the western Sahel occurring mainly during the early Sahel rainy season (i.e., May-June-July-August) and an increase over the central Sahel occurring mainly during the last months of the Sahel rainy season (i.e., August-September-October) (Biasutti 2013;Seth et al. 2013;Monerie et al. 2016;Dunning et al. 2018).
The central/western Sahel zonal contrast is found to be robust in climate projections, because it is simulated by a majority of independent climate realizations, but appears paradoxical in a region that is, at first approximation, considered as homogeneous (Nicholson 2013). However, several studies have shown that the recent observed recovery in Sahel precipitation, which has followed the 1960s-80 s drought (Lebel and Ali 2009;Sanogo et al. 2015;Panthou et al. 2018), has a similar behaviour, i.e., stronger over the eastern than over the western Sahel (Lebel and Ali 2009;Panthou et al. 2018). It has also been suggested to correspond to an observed mode of climate variability, as highlighted by Nicholson and Palao (1993), who have defined two spatially consistent sub-domains, the Sahel and a West Coast domain.
The presence of this zonal contrast in the climate projections could have strong implications for Sahel populations because regionally-adapted policies would have to be done to respond to heterogeneous climate change across the Sahel. For instance, reginal climate models also simulate a zonal contrast in precipitation change at the end of the twentyfirst century, highlighting strong effects of climate change on extreme events, with strong local future modulations on the number of dry days and extreme precipitation events (Diallo et al. 2016;Akinsanola and Zhou 2019). However, the zonal contrast is, as highlighted by Panthou et al. (2018), not fully understood yet. In this study we aim at documenting and understanding the zonal contrast to help providing reliable and usable information to decision makers.
Based on previous studies, we raised several hypotheses that can explain an amplification of the zonal contrast in Sahel precipitation towards the end of the twenty-first century: [1] Several authors have shown that Sahel precipitation anomalies are strongly associated with the magnitude and patterns of the SSTs anomalies (Folland et al. 1986;Palmer 1986;Rowell et al. 1992;Knight et al. 2006;Mohino et al. 2011a;Martin and Thorncroft 2014b). Warming of the SST and changes in SST anomaly patterns could thus be responsible for the zonal contrast in Sahel precipitation.
[2] The zonal contrast could be directly associated with the warming over land. The warming is associated with an increase in evaporation and with changes in moist static energy gradients over land, both having an impact on the monsoon circulation. Besides, a strengthened SHL may favour moisture divergence on its western edge (advecting the relatively dry northerlies) and moisture convergence on its eastern edge (strengthening the westerlies) (Roehrig et al. 2011;Lavaysse et al. 2016), then leading to the zonal contrast. However, Shekhar and Boos (2017) demonstrated that the focus has to be made on the shift of the SHL location, and that the impact of a deepened SHL on Sahel precipitation is unclear.
[3] The zonal contrast could arise from the competitive roles of increases in land surface temperature and changes in SST, combining both hypothesis [1] and [2]. Both can lead to opposite effects but can act together to shape Sahel precipitation changes, as shown in Biasutti (2013) and Gaetani et al. (2017).
[4] Finally the emergence of a strong zonal contrast could be an artefact due to the inability of climate models to simulate the West African monsoon dynamic, as proposed in James et al. (2015) for the western Sahel.
Our goal is to revisit these mechanisms to shed light on their roles on the zonal dipole in Sahel precipitation change. We investigate the evolution of the zonal contrast over both the 20th and the twenty-first century by using a climate model. We separate the total response into a forced response [i.e., the response of a climate model to external forcing like the Greenhouse Gases (GHGs)] and a variation that is associated with the internal climate variability, i.e., in absence of external forcing.
We thus explore the drivers of the zonal contrast with the presence and in the absence of the forced response. On that purpose we use the Large Ensemble of the Community Earth System Model (CESM) (Kay et al. 2015).
The structure of the paper is summarized as follows: Sect. 2 describes the simulations and the methodologies used. In Sect. 3, we revisit the different mechanisms explaining the establishment and evolution of a zonal contrast in Sahel precipitation. A discussion is given in Sect. 4 and a conclusion in Sect. 5.

CESM1
We use the 40 members of the Community Earth System Model (CESM) of the Large Ensemble Project (LENS, Kay et al. 2015). The model is hereafter noted CESM1. CESM1 is a coupled climate model, which includes ocean, land, atmosphere and sea ice components. The atmospheric model is the Community Atmosphere Model (CAM5.2; Hurrell et al. 2013). The land surface component is the Community Land Model version 4 (Lawrence et al. 2011), the ocean model is the Parallel Ocean Program version 2 (Smith et al. 2010) and the sea ice component is the Community Ice Code version 4 (Hunke et al. 2010). All components were run at ~ 1° horizontal resolution (Kay et al. 2015). Each member shares the same radiative forcing scenario (historical evolution in GHG concentration, volcanic forcing and anthropogenic aerosols) for the period 1920-2005 and the Radiative Concentration Pathway 8.5 (hereafter noted RCP8.5) for the period 2006-2100. The RCP8.5 emission scenario describes a future trajectory in anthropogenic emissions eventually resulting in 8.5 W m −2 radiative forcing in 2100 (Meinshausen et al. 2011). To build the large ensemble, the 40 members only differ by their initial atmospheric conditions. Details of the simulations are found in Kay et al. (2015). Model's outputs are publicly and freely available on the NCAR website (www. cesm.ucar.edu/proje cts/commu nity-proje cts/LENS/datasets.html). CESM1 simulates an increase in Sahel precipitation over the central Sahel and a decrease in precipitation over the western Sahel (Monerie et al. 2017) (and Fig. 1), as in the CMIP5 ensemble, and is therefore particularly suitable for this study.

AMIP simulations
The simulations of the Atmospheric Model Intercomparison Project (AMIP), within the CMIP5 exercise (Taylor et al. 2012), are also considered. The idealised AMIP experiments are complementary to the simulations of the CESM1 large ensemble, and are used to identify the effects of the heterogeneous global SST warming, to the homogeneous SST warming, and to the land warming only. In the AMIP experiments, the models have been run in an atmospheric-only configuration. For the reference simulation (hereafter noted amip), the atmosphere is forced by prescribing observed SST and sea ice conditions for the period 1979-2008, and by the observed evolution in the atmospheric composition (e.g., CO 2 atmospheric composition and aerosol emission), and the land use and land cover changes. Three other AMIP-type simulations are analysed. In the amip4K, a uniform increase of 4 K is added to the global SST forcing; in the amipFuture simulation, an SST pattern of climate change is scaled to have a global mean increase of + 4 K. The pattern is taken from the CMIP3 ensemble mean change in the 1%/yr CO 2 increase experiments at quadrupling (Bony et al. 2011). In the amip4xCO2, there is a quadrupling of the CO 2 atmospheric concentration with prescribed SST. All other boundary conditions remain the same as in the AMIP control simulation (i.e., amip). We used a set of 8 models whose outputs were available at the time of the study (MRI-CGCM3, IPSL-CM5A-LR, IPSL-CM5B-LR, CNRM-CM5, HadGEM2-A, MIROC5, MPI-ESM-LR and MPI-ESM-MR).

Observations
To discuss the ability of CESM1 to simulate the zonal contrast in Sahel precipitation we use several observations. Observed precipitation data are given by the Global Precipitation Climatology Center (GPCC) version v7 (Schneider et al. 2014), available from 1901 to present on a 0.5° × 0.5° longitude on a global grid. Data of the Climate Research Unit (CRU; version 4.01) allow a good estimation of continental rainfall at a 0.5° horizontal resolution from 1901 to 2016 (Harris et al. 2014), Data from the University of Delaware (UDEL; version 4.01; Willmott et al. 2001) are at the 0.5° horizontal resolution from 1901 to 2014. The Precipitation Reconstruction over Land (PREC/L) is at the 1.0° horizontal resolution and is available from 1948 to present (Chen et al. 2002). The twentieth century reanalysis is produced by assimilating only surface pressure observations of synoptic pressure, monthly sea surface temperature and sea ice distribution (Compo et al. 2006). The twentieth century reanalysis is available from 1851 to present and the ensemble mean of the 56 generated members is use to evaluate CESM1 simulations of the Sahel precipitation and teleconnections. For surface air temperature we use the Berkeley Earth Surface Temperature dataset (BEST), which is an interpolated dataset of surface air temperature anomalies (with respect to the period 1951-1980) over land and ocean built from temperature records by using statistical techniques (Rohde et al. 2013). These data are available from 1850 to present at a 1° × 1° horizontal resolution.

Estimating the forced response and internal climate variability
A large ensemble (LENS) allows providing a broad statistical sampling of coupled internal climate variability (IV hereinafter). A climate model simulation can exhibit internal fluctuations of climate, but these oscillations are more likely to be out of phase among the 40 members, when not mostly driven by external forcing. The forced response can be defined as the common characteristic to all members of the ensemble. Hence, the forced response ( Δ v ) is estimated by the ensemble mean of the total response ( Δ mv ) for a given variable ( v ) by averaging over all members (m) (Deser et al. 2014). Internal climate variability ( IV mv ) is obtained by subtracting the forced response from each individual member, i.e., We use the IV component to define how CESM1 simulates anomalies in a variable v, in absence of external forcing. We assess the impact of climate change using the ensemble mean. Held and Soden (2006) have shown that climate change leads to a general weakening of the tropical circulation, which they explained through decomposing tropical precipitation, by making the assumption that precipitation can be approximated by,

Decomposing precipitation change
where P is precipitation, M* is a proxy for convective massflux from the boundary layer to the free troposphere (with M* = P/q) and q is near surface specific humidity (see also   Kent et al. 2015). Following Chadwick et al. (2016), changes in precipitation can then be decomposed, as Precipitation change is then reformulated in terms of thermodynamic ( ΔP therm ), dynamic ( ΔP dyn ) and a cross nonlinear ( ΔP cross ) components.
ΔP therm is the thermodynamic component of ΔP , that is the change in precipitation associated with changes in specific humidity, at constant M*. ΔP dyn is the dynamic component and represents the contribution to ∆P due to changes in circulation (M*) at constant q. ΔP cross is the change in precipitation, explained by simultaneous changes in both circulation and specific humidity (Chadwick et al. 2013).
In a warmer climate, q is expected to increase at approximately the Clausius-Clapeyron rate of 7% K −1 (see Held and Soden 2006), suggesting a strengthening of the contribution of the thermodynamic term for future climate conditions. However, precipitation increases at a weaker rate than specific humidity (Mitchell et al. 2007;Stephens and Ellis 2008;Chadwick et al. 2013) in climate models, due to a limited increase in tropospheric radiative cooling (Stephens and Ellis 2008), constraining the dynamic term to decrease with warming (i.e., M* = P/q) (Held and Soden 2006;Vecchi and Soden 2007a;Chadwick et al. 2013). According to Vecchi and Soden (2007) the weakening of the tropical mean circulation is associated with less precipitation over land, globally, and with a strong shift of the circulation over the Pacific Ocean and in the Indian Ocean (not shown).
Further decomposition of ΔP dyn allows to document changes that are due to the weakening of the tropical mean circulation ( ΔP weak ) and to a shift in the pattern of the circulation ( ΔP shift ), as With, where, α = − (tropical mean ∆M*/tropical mean M*). α is scaled by the strength of the mean tropical circulation, suggesting that changes in the tropical mean circulation are inversely proportional to the climatology mass-flux field, as shown in Chadwick et al. (2013).
The thermodynamic term can be partitioned into a component that corresponds to an increased maritime moisture advection onto land, and a residual that is associated with changes in circulation, evaporation and vertical mixing (Chadwick et al. 2016;Rowell and Chadwick 2018), by decomposing the change in near-surface specific humidity. A nonuniform circulation weakening could manifest itself through an increase in convection in some regions and a decrease in convection in other regions, and will be seen, through the decomposition, as a shift in the circulation. We made the assumption that this approximation is reasonable, although the interpretation of mechanisms would need to be reviewed.
The change in land specific humidity, with no change in circulation, and due to an increase in advection from the Ocean to land, in response to the change in specific humidity over the Ocean, is obtained by scaling historical mean moisture over land by the zonal mean fractional change in moisture over the Oceans (Chadwick et al. 2016).
where, q ocean_future q ocean_historical is the zonal mean fractional change in moisture over the ocean, at the same latitude, with no change in the circulation. This approximation implies that changes in maritime advection mainly have a zonal source, as shown by Lélé et al. (2015). The change in moisture advection from the ocean onto land grid points is then obtained as,

And the thermodynamic advection term is
The residual thermodynamic is obtained as the change in precipitation associated with the change in specific humidity that is not due to the increase in maritime moisture advection ( Δq adv ), over land, at constant dynamic, i.e.,

And
The decomposition is performed at the monthly time step (Chadwick et al. 2016;Rowell and Chadwick 2018) prior to compute the seasonal means and the area-weighted averages over a certain region. Changes in precipitation are evaluated by comparing precipitation of the period 2060-2099, under the RCP8.5 emission scenario, to the reference period 1960-1999, under the historical emission scenario. Timeevolutions of the anomalies are also assessed relative to the year 1920. We used the 925 hPa specific humidity for the decomposition because the 2-m specific humidity was only provided over the land areas in CESM1, but checked that over land this does not affect the results (focusing on ΔP therm , ΔP dyn and ΔP cross ) (not shown).

The African rainfall pattern index
The African Rainfall Pattern Index (ARPI; Monerie et al., 2012) is a proxy of the zonal dipole and allows analysing the evolution of the Sahel precipitation zonal contrast. We first compute two Sahelian indices: the first is defined over the western Sahel (WS; 10°-20°N; 20°W-5°W) and the second is defined over the central Sahel (CS; 10°-20°N; 5°W-20°E).
The central and western Sahelian indices are standardized (stan) to avoid western Sahel precipitation to dominate variations of the ARPI metric. In the standard ARPI, precipitations are averaged over the July-August-September (JAS) season.

Computing the strength of the Saharan heat low
We compute an index to account for the strength of the SHL. We followed Dixon et al. (2016) and compared Sahara 925 hPa geopotential height ( ZG Sahara ) with tropical mean (20°S-20°N) 925 hPa geopotential height ( ZG tropical_mean ) . As is Dunning et al. (2018) we compute the weighted area average over the Sahara between 15°W and 30°E, and 20°N and 30°N.
We highlight that negative (positive) values of the index correspond to intense (weak) SHL. Using reanalysis datasets, Shekhar and Boos (2017) have shown that the impact of a strengthened SHL on the Sahel precipitation is unclear, whereas the effect of the meridional shift of the SHL is suggested to be robust. Therefore, we also analyse meridional shifts in the location of the SHL, by computing the zonal average of ZG-ZG tropical_mean , between 15°W and 30°E, following Shekhar and Boos (2017). The location of the SHL is identified by selecting the latitude of the minimum value of ZG-ZG tropical_mean , after cubic splines interpolation.
In addition, we account for the 3D structure of the SHL by using another SHL definition, computing the difference between geopotential height at 700 hPa and at 925 hPa (i.e., ZG 700 -ZG 925 ), averaged over the Sahara, hereafter called low-level atmospheric thickness (LLAT; Lavaysse et al. 2009), after Shekhar and Boos (2017).

Sahel precipitation forced response to global warming
The spatial pattern of the precipitation change simulated by CESM1 is similar to the one of the CMIP5 multi-model ensemble-mean (Akinsanola and Zhou 2018; among others): an increase in precipitation over the central Sahel and a decrease in precipitation over the western Sahel at the end of the twenty-first century, relative to the end of the twentieth century (Fig. 1a). Over the central Sahel, both the increase in precipitation, north of 10°N, and the decrease in precipitation, south of 10°N, denote a northward shift of the West African monsoon. Over the western Sahel, monsoon precipitation tends to move southward. Although at a weaker intensity, similar changes have been obtained at a near-term time horizon (i.e., over the period 2010-2049) (Monerie et al., 2017) and with the RCP4.5 emission scenario (Monerie et al. 2012;Akinsanola and Zhou 2018), and the CMIP6 Shared Socioeconomic Pathways (Almazroui et al. 2020;Monerie et al. 2020).
It is worth noting that precipitation change simulated by CESM1 is particularly robust, as it is reproduced by most of the members (see the stippling on Fig. 1a). It has been argued that uncertainties in simulating tropical precipitation change at the end of the twenty-first century in CMIP5 models are mostly due to differences in climate models parameterizations rather than to internal climate variability (Hawkins and Sutton 2011). This would explain why the spread among the CESM1 members, is much reduced, in comparison to the multi-model CMIP5 ensemble spread, as also shown in Monerie et al. (2017).
In CESM1 the decrease in western Sahel precipitation is particularly strong from June to August (hereafter called JJA). The increase in central Sahel precipitation occurs at the end of the rainy season (August-October; hereafter called ASO) (Fig. 1b, c). These seasonal changes have also been shown within CMIP5 simulations (Biasutti 2013;Seth et al. 2013; among others), and have been associated with an anomalously late monsoon onset over the western Sahel and late monsoon cessation date over the central Sahel (Dunning et al. 2018). Therefore, we also assess impact of climate change on the longer JJASO period, and find similar results than for JAS, for both sign and magnitudes of precipitation change (Fig. S1).

Evolution of the Sahel zonal dipole
To measure the evolution of the strength of the precipitation zonal contrast we compute the ARPI values following Eq. 15 (Sect. 2.2.3). Note that the definition of the central and western Sahel is slightly different from that in Monerie et al. (2012) to better account for the zonal precipitation contrast of CESM1. In addition, a refined index (ARPI refined ) has also been used, selecting two different seasons (i.e., JJA for the western Sahel and ASO for the central Sahel), instead of JAS for both areas, to better account for the temporal inhomogeneity. ARPI anomalies are computed relative to the year 1920.
The ARPI ensemble mean shows that the zonal contrast is projected to strengthen in the future and in JAS, relative to the twentieth century (Fig. 2a), as shown within CMIP3 (Monerie et al. 2012) and CMIP5 (Biasutti 2013; among others) simulations. The forced ARPI variations are highly correlated with variations in global mean surface temperature with a correlation coefficient of 0.89, suggesting that forced APRI variations might be explained through changes in global mean surface temperature. Besides, it is worth noting that impact of climate change emerges from internal climate variability after 1980, highlighting a strong impact of climate change on the zonal contrast ( Fig. 2a grey shading).
The refined ARPI index shows a larger increase than the original ARPI index, maximizing the impacts of climate change on the zonal precipitation dipole. However, the time evolution of the refined index shows a very similar behaviour to that of the original ARPI. We thus focus on only the ARPI index, which can be more easily used than the refined ARPI index, with other models for which the timing in precipitation change could be different from that in CESM1. Results obtained with the refined index are presented in the supplementary material (Fig. S2).

The ARPI index and internal climate variability
We begin by exploring the hypothesis that the zonal dipole contrast could be caused by changes in SSTs in areas key for Sahel rainfall. Changes in SSTs in such areas would, in turn, be shaped by the forced response. Therefore, we show the SST-related factors controlling the zonal dipole in the absence of a forced response (i.e., the internal climate variability component of the ARPI index), before showing how the forced response could impact these factors.
We regress the time-series corresponding to the IV component of the ARPI time series onto the IV component of surface air temperature, over the ocean (Fig. 2b) [see equation (1)]. Note that for computing the ARPI-IV time-series the ensemble-mean shown in Fig. 2a is removed from the ensemble members and this results in time-series without a long temporal trend. The regression pattern obtained in Fig. 2b i.e., a cooling over the eastern Pacific Ocean and warming over the maritime continent, is very similar to the negative phase of the El-Niño Southern Oscillation (ENSO) pattern in the tropical Pacific and suggests a possible contribution of ENSO in shaping the zonal dipole in the Sahel precipitation (Rowell et al. 1992;Joly et al. 2007;Fontaine et al. 2011a;Mohino et al. 2011b). In addition, the regression pattern also shows an Atlantic cross-equatorial gradient of SSTs, suggesting a role of this structure in the Sahel zonal contrast in the CESM1 model. To further understand the (a) (b) Fig. 2 a Ensemble mean ARPI (red; i.e., the difference between the standardized JAS central Sahel precipitation minus the JAS western Sahel precipitation) and refined ARPI (blue; i.e., the difference between the standardized ASO central Sahel precipitation minus the JJA western Sahel precipitation) change, relative to the year 1920. Blue and red shadings indicate the ensemble member spread, as defined by 2 times the standard error ( ensemble−member √ n , with n = 40 ). Intermember spread in the internal variability component of ARPI (grey shading; i.e. obtained by inter-member standard deviation) and of the refined ARPI (black lines). b Ensemble-mean of internal component of ARPI regressed onto surface temperature internal variability, over the period 1920-2100, as first computed for each member. Stippling indicates that at least 80% of the members are simulating a regression that is significantly different to zero, according to a Spearman's rank test, at the 95% confidence level relationship between the Atlantic cross-equatorial SST gradient, ENSO, and Sahel precipitation, we analyse changes in precipitation associated with SSTs over the tropical Pacific and Atlantic oceanic basins in the next subsections, separately for the western and central Sahel. Firstly, we verified that ENSO and Atlantic cross-equatorial gradients are independent, with a correlation coefficient value of r = − 0.13 between both indices. Secondly, we assess the relationship between western/central Sahel precipitation and surface temperature in the absence of external forcing (i.e., using the IV component), that is, how natural climate variability shapes the zonal dipole. This will help to elucidate how the external forcing could impact the Sahel zonal dipole via the SST-rainfall relationship.

Internal climate variability
The internal component of western Sahel precipitation (IV pr ) is regressed onto the internal component of ocean surface air temperature (IV tas ) over the long period 1920-2100 for each member individually. Changes in precipitation are strong in June-July-August (JJA) (Fig. 1c) so we analysed western Sahel precipitation over the JJA period in this section. The average of the 40 regression patterns is displayed in Fig. 3a [see equation (1)]. Results reveal an important role of the Atlantic cross-equatorial temperature gradient on the western Sahel precipitation IV component (IV pr ;  Fig. 3a). We only found significant correlations over the Atlantic Ocean, showing that the Atlantic cross-equatorial temperature gradient is key for the western Sahel (Fig. S3). The Tropical South Atlantic/Tropical North Atlantic (hereafter noted TSA/TNA) gradient is computed through the difference between the area-weighted temperature average over the TSA and TNA regions (see boxes on Fig. 3a). The regression of the IV component of the TSA/TNA gradient onto the IV component of the precipitation field highlights the connection between the tropical Atlantic Ocean and the meridional migration of the monsoon precipitation (Fig. 3b). This is consistent with the known observed modes of variability, for which an anomalously warm eastern equatorial Atlantic Ocean (or Gulf of Guinea) is associated with a decrease in Sahel precipitation (Nicholson 2005;Joly et al. 2007;Losada et al. 2012;Fontaine et al. 2011a). Analysis also reveals that the TSA/TNA impact is mostly mediated by changes in atmospheric circulation ( ΔP dyn ) (Fig. 3c), while changes in the thermodynamic component are much more moderate (Fig. S4). This analysis is repeated by using the raw western Sahel precipitation (i.e., internal + forced components) indices and the raw fields of surface temperature, which also shows a strong impact of the cross-equatorial Atlantic temperature gradient on western Sahel precipitation, over the historical period 1920-2005 (Fig. 3a). This suggest that the TSA/ TNA gradient could have a role in the total response of the zonal precipitation contrast over the Sahel, with the forced response driving changes in the TSA/TNA gradient. This will be addressed in the following section.

Evolution of the tropical Atlantic thermal gradient
Changes in precipitation can be strongly impacted by changes in the patterns of the SSTs. For instance, Hoerling et al. (2006) have shown that long-term changes in the temperature contrast, between the North and the South Atlantic gradient, can lead to changes in Sahel precipitation. An increase in the TSA/TNA temperature gradient, due to global warming, could then lead to long-term changes in western Sahel precipitation.
In JJA, the TSA/TNA gradient increases strongly over the 1920-2000 period (Fig. 3d). The positive trend is statistically different to zero according to a Spearman's rank correlation test (at the 95% confidence level). Over the 1920-2000 period, western Sahel precipitation decreases with time (Fig. 4a), in agreement with a strengthening of the TSA/TNA gradient. This suggests that the cross-equatorial Atlantic gradient is one of the drivers of the ARPI evolution over the twentieth century.
Impact of climate change on the TSA/TNA thermal gradient is unclear, and the gradient is not stronger over the 2060-2099 period than over the 1960-1999 period. Unlike the TSA/TNA thermal gradient, western precipitation decreases over the twenty-first century, showing that the TSA/TNA thermal gradient is unlikely to be one of the main causes of the projected strengthening of the zonal contrast.

Internal climate variability
We regress the internal component of central Sahel precipitation (IV pr ) onto the internal component of ocean surface air temperature (IV tas ) over the long period 1920-2100. As changes in central Sahel precipitation are stronger in August-September-October (ASO) than throughout the year (Fig. 1c), we focus the analysis on the change in central Sahel precipitation on the ASO period. ASO central Sahel precipitation is strongly related to ENSO in CESM1 (Fig. 5a), highlighting the strong impact of the Tropical Pacific Ocean on the Sahel precipitation zonal contrast in this model. A zonal dipole appears between the eastern Pacific Ocean (i.e., Niño3.4 box) and the SSTs around the maritime continent (MAR hereinafter). We document ENSO evolution over the twentieth and twenty-first century by computing the difference between the Niño3.4 and MAR indices, and we estimate the impact of the ENSO by regressing the internal component variability of the Pacific zonal temperature gradient onto the internal component variability of the precipitation. The regression pattern shows that a negative (positive) ENSO phase is associated with a precipitation increase (decrease) over the central Sahel (Fig. 5b). ENSO is thus able to force a shift of the rain belt over Africa between 0° and 30°E. The impact of the Tropical Pacific Ocean is mostly mediated by changes in atmospheric dynamic (Fig. 5c) rather than changes in the thermodynamic (Fig. S5). The analysis performed by computing the regression between the en3.4 index and Sahel precipitation yields the same pattern as in Fig. 5b (not shown), showing that the EN34-MAR difference provides a robust estimate of the impacts of ENSO over Sahel precipitation. Similar results are obtained with the raw time series over the historical period (green contour in Fig. 5a), so in the following we explore the role of the EN34/MAR gradient in the forced response.

Evolution of the tropical Pacific thermal gradient
The EN34/MAR gradient is projected to increase with time ( Fig. 5d) in ASO. As a consequence, the impact of the Pacific thermal gradient on central Sahel precipitation increases over the 20 th and the twenty-first century. However, our results indicate that the impact of the Pacific thermal gradient in CESM1 is moderate on central Sahel precipitation and does not impact western Sahel precipitation (Fig. 5b). Moreover, central Sahel precipitation is projected to increase (Figs. 1c and 4b) while a strengthening of the EN34/MAR gradient leads to a decrease in central Sahel precipitation, showing that this thermal gradient is unlikely to explain the projected changes in central Sahel precipitation.

Decomposing changes in Sahel precipitation
We have shown that, when not influenced by the anthropogenic forcing, the Sahel zonal contrast is associated with SST patterns of natural climate variability. However, results shown in the previous sections suggest that the long-term change in the Sahel zonal contrast cannot be explained by long-term changes in surface temperature patterns over the ocean, as defined in Figs. 2b and 3b i.e. the long-term changes in the TSA/TNA and EN34/MAR temperature gradients. Therefore, we take a step forward at better understanding the zonal contrast in this section by decomposing the changes in Sahel precipitation (Sect. 2).

Western Sahel precipitation change
Over the western Sahel, precipitation decreases continuously from 1920 to 2100 (Fig. 4a). It is worth noting that the decrease in precipitation is occurring at different rates, slower during the twenty-first century than during the twentieth century (see the dotted line in Fig. 4a). In spite of the local increase in surface air temperature (Fig. S6), evaporation decreases over the western Sahel. This is mostly due to the fact that, in the semi-arid region, changes in evaporation mostly depend on the changes in soil moisture, rather than in temperature (Cook et al. 2014). The change in moisture flux convergence also undergoes a slight decrease in its changing rate, weaker during the 21st than during the twentieth century. Moreover, the moisture flux convergence change is stronger than the local precipitation recycling change, showing that the decrease in precipitation is mostly due to a decrease in moisture transport by the trade winds. Changes in moisture flux can be due to either or both changes in moisture air content and atmospheric circulation. Therefore, we decompose precipitation change into its dynamic, thermodynamic and cross nonlinear components.
By using the precipitation decomposition described in Sect. 2, we show that over the western Sahel the decrease in precipitation is mostly associated with the ΔP dyn component (Fig. 6a), i.e., dynamic changes. The early decrease in ∆ P is only due to changes in ΔP dyn while ΔP therm and ΔP cross impacts are only distinguishable during the twentyfirst century (Fig. 6a). Over the western Sahel the decrease in ΔP dyn is linear through time and the change in the ∆ P rate is explained by ΔP therm that allows damping the decrease in precipitation, over the twenty-first century.
The 1920-1980 change in precipitation is only associated with changes in the pattern of the atmospheric circulation (i.e.ΔP shift ) (Fig. 6c)  , with n = 40) Fig. 6 Ensemble mean change in a precipitation (∆ P; red), ΔP dyn (orange), ΔP therm (green) and ΔP cross (purple), in c ΔP dyn (orange), ΔP weak (dark orange) and ΔP shift (pink), e ΔP therm (green), ΔP qadv (blue) and ΔP qres (brown), and in g ∆ P (red), ΔP qadvw (orange), ΔP shiftres (light blue) and ΔP cross (purple), for the Western Sahel in JJA. b-d, f-h same as in a-c, e-g but for the central Sahel and in ASO. Anomalies are given relative to the year 1920. Units are in mm.day −1 with a strengthening of the moisture advection from the ocean to land ( ΔP qadv ) (Fig. 6e). The warming of the tropical SSTs is associated with both an enhanced moisture advection from the oceans ( ΔP qadv ), which increases precipitation over Sahel (Fig. 6e), over the ascent regions, and with a weakening of the tropical circulation ( ΔP weak ) that reduces precipitation (Ma et al. 2018) (Fig. 6c). Thus, both ΔP qadv and ΔP weak are strongly negatively correlated and there is a physical cancelation between ΔP qadv and ΔP weak (Kent et al. 2015;Rowell and Chadwick 2018). Therefore, we combined both terms, with ΔP qadvw = ΔP weak + ΔP qadv . ΔP res is a residual that is associated with changes in circulation, evaporation and vertical mixing (Chadwick et al. 2016). We combined ΔP res and ΔP shift that are associated with changes in the dynamic: ΔP shiftres = ΔP shift + ΔP res . ΔP weak and ΔP qadv partially cancel out and ΔP qadvw is moderate and only increases after 2010 (Fig. 6g). Finally, the change in western Sahel precipitation is mostly associated with a shift in the circulation.

Central Sahel precipitation change
Central Sahel precipitation tends to decrease over the period 1920-1960, to increase over the period 1980-2080 and to be stabilized at the end of the twenty-first century (Fig. 4b).
Over the central Sahel, the change in precipitation is apparently explained by ΔP therm , because ΔP dyn and ΔP cross are moderate (Fig. 6b). However, the change in the atmospheric dynamic is moderate because its two components, ΔP shift and ΔP weak partially cancel each other (Fig. 6d). Moreover, the pattern of ΔP dyn denotes a northward shift of the West African Monsoon east of 0°W (Fig. S7), with positive (negative) anomalies north (south) of 15°N, leading to a rather low change when averaged between 10°N and 20°N. ΔP weak decreases after the 1980s while ΔP qadv increases after the 1980s, and ΔP qadvw is then small over the central Sahel (Fig. 6h) suggesting that the change in central Sahel precipitation is mostly due to shifts in the circulation (Fig. 6h).
ΔP qadv and ΔP weak are relatively homogeneous over West Africa, but with strong values over the Guinean coast and the Cameroon highs (Fig. S7), where climatological precipitation is the highest. ΔP qadvw is then homogeneous over the Sahel. ΔP qadvw is slightly positive because ΔP qadv is stronger than ΔP weak (Fig. 7a , d). Both ΔP shift and ΔP res have a similar pattern (not shown), and ΔP shiftres is associated with a zonal dipole in precipitation change over the Sahel (Fig. 7b, e). As in Chadwick et al. (2016), the shift in the tropical circulation mostly explains thus the pattern in precipitation change. The shift in the circulation is stronger in JJA than in ASO over the western Sahel, explaining why the western Sahel precipitation decrease is stronger in the early rainy season. Changes in the cross term ( ΔP cross ) are not negligible and also show a dipole (i.e., a decrease over the western Sahel and an increase over the central Sahel) (Fig. 7c, f). However, the pattern of ΔP cross exhibits some dissimilitude with the pattern of ΔP and its intensity is lower than the dynamic term, over the Sahel.

Role of regional changes
Shifts in the circulation are strong over the Sahel. Changes in the circulation can be associated with both changes in sea surface temperature, both locally and globally, and with the warming of the Saharan desert that is associated with a strengthening of the monsoon circulation (Chadwick et al. 2019). To assess the respective roles of the SST warming and of the warming over land we use the AMIP simulations (Taylor et al. 2012), which are described in Sect. 2.

Impacts of sea surface temperature change and CO 2 concentration increase
The net impact of a homogeneous warming of the SSTs is a southward shift of the West African monsoon and a decrease in Sahel precipitation (Fig. 8a). Gaetani et al. (2017) found that a warming of the SSTs leads to a decrease in Sahel precipitation through large-scale changes of the atmospheric circulation, resulting in increased subsidence over the Sahel. They showed that the large-scale circulation is partly driven by the warming of the tropical and equatorial Atlantic Ocean, where ascents are associated with largescale changes over the Indian Ocean, in a Walker-type circulation. In addition, Hill et al. (2017Hill et al. ( , 2018 have shown that the decrease in Sahel precipitation is associated with an increase in meridional moist static energy gradient spanning the Sahel and northerly dry air advection from the Sahara desert to the Sahel. The amipFuture simulation shows that a heterogeneous increase of the SSTs (Fig. S8) is also associated with a decrease in precipitation over the Sahel (Fig. 8b). Though the difference between the amip4k and amipFuture is not negligible, two-thirds of the precipitation anomalies are due only to the homogeneous warming of the global ocean, shedding light on the prominent importance of the tropical mean warming of the SSTs in comparison to a change in SST gradients. For both homogeneous and heterogeneous SST warming, changes in precipitation are mostly due to changes in the location of the atmospheric circulation, with the warming of the SST leading to a southward shift of the monsoon (Fig. 9a and Fig. 9b; and Fig S8). Besides, a warming of the SST is associated with a strengthening of the maritime moisture advection, which is partly offset by the weakening of the tropical mean circulation (i.e., ΔP qadvw is close to zero) (Fig. S10). The impacts of both homogeneous and heterogeneous warming of the SSTs on ΔP pattern are explained by ΔP shift (Fig. 9a, b) and ΔP shiftres (Fig. S9). Interestingly, a homogeneous warming of the SST does not only lead to a decrease in Sahel precipitation, but also to a southward shift of the monsoon.
The amip4xCO2 simulations show a relatively large warming over land in JAS over the northern Hemisphere, including the Saharan desert, north of the Sahel (Fig. S8). The maritime advection (i.e.ΔP qadv ) and the change in the tropical mean circulation (i.e., ΔP weak ) are low (not shown) because SSTs are prescribed and do not increase in the  S8). The monsoon circulation shifts northward and Sahel precipitation increases (Figs. 8c and 9c). This is consistent with the results of Gaetani et al. (2017), who have shown that the moist static energy content increases in the lower troposphere, in response to an increase in land temperature, driving the monsoon system to shift northward. Changes in SSTs are associated with a strong decrease in precipitation over the western Sahel, and with a southward shift of the ITCZ over the Atlantic Ocean (Fig. 8a,  b). The increase in greenhouse gases concentration is associated with an increase in precipitation over the central Sahel, while changes over the western Sahel are not robust (Fig. 8c). A zonal dipole in precipitation over the Sahel does not clearly appear due to the warming of the SST nor to the warming of northern Africa. The addition of the impacts of a warming of the SSTs (amip4K), that is the slow response to climate change, and of the increase in greenhouse gases (amip4xCO2), that is the fast response to climate change, shows a zonal dipole with a northward shift of the monsoon over the central Sahel and a southward shift of the monsoon over the western Sahel (Figs. 8d and 9d). This suggests that the zonal dipole is neither due to changes in SSTs nor to the warming of northern Africa, but to a combination of both.

Role of the strength and location of the Saharan Heat low
No clear changes of the SHL index [Eq. 11] are obtained during the early twentieth century, and the decrease in western Sahel precipitation is thus not associated with changes of the Heat Low strength and location between the year 1920 and the year 1980 (Fig. 10a, c). After 1980, the decrease in geopotential height field at 925 hPa (ZG925) is robust as reproduced by all members. The post-1980s deepening of the Heat Low (Fig. 10a) is associated with an increase in surface air-temperature and a decrease in sea level pressure over northern Africa (Fig. S6). The regression pattern of the IV component of SHL index onto IV component of precipitation shows a relationship between the strength of the SHL and central Sahel precipitation (Fig. 10b). A strengthening of the SHL (i.e.,negative anomalies in ZG925) is associated with an increase in precipitation over the central Sahel and over West Africa, South of 10°N (Fig. 10b). To confirm this result, we also use a second metric, accounting for the 3D structure of the shallow circulation, with the LLAT index (i.e., ZG700-ZG925). LLAT increases with time, also showing a strengthening of the SHL in a warmer climate (Fig.  S11). However, the regression between the LLAT index and Sahel precipitation does not show significant impacts due to a strengthened heat low over Sahel (Fig. S11). This is consistent with Shekhar and Boos (2017), which have shown that impact a change in the strength of the SHL is unclear on Sahel precipitation. Shekhar and Boos (2017) have shown that the impact of a shift in the SHL location is associated with robust changes over the Sahel in reanalysis. In a warmer climate, SHL is projected to move northward (Dunning et al. 2018) (Fig. 10c). A northward shift of the SHL is also shown with the LLAT index (Fig. S11). A northward shift of the Saharan low is associated with an increase in precipitation over the central Sahel in CESM1 (Fig. 10d  and Fig. S11). The projected northward shift of the SHL could therefore lead to an increase in precipitation over the Sahel in a warmer climate. However, the impact of a northward shifted heat low is not significant over the western Sahel (Fig. 10d). Therefore, we suggest that changes in the location of the SHL contributes to the zonal contrast because impacting more strongly the central than the western Sahel, causing heterogeneous changes in precipitation (Fig. 10d). This is also (a) (b) (c) (d) Fig. 9 As in Fig. 8 but for ΔP shift consistent with the results of the 4xCO2 simulations, in which an increase in northern African surface temperature is associated with a strengthening of the SHL (as shown in Gaetani et al. 2017), and with an increase in precipitation over the central Sahel only (Fig. 8c). However, impacts on the SHL cannot explain, by themselves, the establishment and future evolution of the zonal contrast in precipitation over the Sahel.

Discussion
The results obtained with the CMIP5 amip simulations might not be representative of the response of CESM1. For instance, the SST pattern used to force the amipFuture simulations is characterised by a future strengthening of the TSA/TNA temperature gradient (Fig. S9) that is not found in CESM1. Within the amip ensemble, we analysed the outputs of the CCSM4 model, a model developed at NCAR, and that produces the same pattern in precipitation , with n = 40 ). b Ensemble mean of precipi-tation regressed onto the SHL index (mm.day −1 /m) in JAS over the period 1920-2100 for the internal variability component of both precipitation and SHL index, as first computed for each member. c as in a but for the location, in latitude, of the SHL. d as in b but with precipitation regressed onto the latitudinal location of the SHL (mm. day −1 /°latitude). Stippling indicates that at least 80% of the members are simulating a regression that is significantly different to zero, according to a Spearman's rank test at the 95% confidence level change as in CESM1, under the RCP8.5 emission scenario. The warming of the SSTs (amip4K and amipFuture) also produces in CCSM4 a southward shift of the monsoon circulation and a decrease in Sahel precipitation (Fig. S12). These results bring confidence on the fact that the results obtained within the CMIP5 ensemble could be similar to CESM1, that a warming of the SST (increase in GHG concentration) is associated with a decrease (an increase) in Sahel precipitation in CESM1. Meanwhile, we might have underestimated the importance of SST gradients on Sahel precipitation with the AMIP simulations. The forcing that is imposed in amipFuture is defined from an ensemble mean of the CMIP3 simulations, and anomalies are expected to be of weaker magnitude than in each model. By analysing the anomalies in western Sahel precipitation, for each month of the year, we also show that the decrease in precipitation, due to the warming of the SST is the strongest in July-August, during the first half of the rainy season, while the land warming is associated with an increase in precipitation, later in the year, from August to October (Fig. S13). The sign (negative) and timing (during the first half of the rain seasons) of western Sahel precipitation change, could thus be due to the warming of the SSTs.
We assess whether the zonal contrast in precipitation over the Sahel in model simulations is also a mode of variability in observations. We performed an Empirical Orthogonal Function (EOF) on observed JAS Sahel precipitation (20°W-20°E; 10°-20°N), using several datasets, i.e., CRU, GPCC, PRECL, UDEL and 20CR. Results indicate that the zonal contrast is a mode of variability in observations, but only explains between 6 and 9% of the total precipitation variability (Fig. S14). The same analysis is repeated using the twentieth century reanalysis, for which the zonal contrast has a stronger importance, explaining 20% of the total precipitation variability. This may suggest that the zonal contrast could then be overestimated in climate models. The regression of the observed ARPI index onto the surface air temperature shows contrasted results between CESM1 and observations (Fig. 2, Fig. S15). Differences between CESM1 and observations at simulating the relationship between SST and ARPI could be due to the inability of CESM1 to simulate the Sahel precipitation variability (Monerie et al. 2017). Therefore we cannot exclude that simulated changes in Sahel zonal contrast could be due to biases in climate models, as suggested by James et al. (2015) for western Sahel precipitation.
The ARPI could be used as a metric to evaluate and select models that are able or unable to simulate Sahel precipitation variability, as done in Cook and Vizy (2006) with the Sahel/Guinean coast Sahel precipitation dipole. However, we show here that the zonal dipole is associated with different mechanisms when accounting or not for the forced response. Unlike in a warmer climate, the simulated zonal contrast, over the historical period, is mainly associated with the thermal gradients over different oceanic basins (i.e., variations in the tropical Atlantic and Pacific Ocean). This could limit the possibility to select a set of the best models based on their ability to simulate the historical observed climate variability since climate variability could be associated with mechanisms of different magnitudes in both historical and warmer climate.

Conclusion
We used the Community Earth System Model CESM1 (Kay et al. 2015) Large Ensemble to investigate the zonal dipole in Sahel precipitation. The large ensemble consists of a set of 40 simulations and allows analysing the simulated climate over the twentieth century (under the historical scenario) and the twenty-first century (under the RCP8.5 emission scenario). Moreover, with 40 members we can robustly estimate the forced response and internal climate variability. We confirm that a zonal contrast (i.e., western/central) in Sahel precipitation is a robust feature of the CESM1 projections.
We clarified here that: Over the central Sahel, precipitation increases due to a northward shift of the monsoon circulation ( ΔP shiftres ). The moistening of the low-level atmosphere through moisture advection (the thermodynamic ΔP qadvw component) compensate the decrease in precipitation that is associated with the weakening of the tropical mean circulation ( ΔP weak ), leading to a moderate increase in precipitation ( ΔP qadvw is positive) (Fig. 6). Western Sahel precipitation decreases strongly between June and August over both the 20th and the twentyfirst century. These changes are mostly contributed by a change in the atmospheric dynamics, while thermodynamic changes damp the decrease in precipitation (Fig. 6). During the historical period the decrease in western Sahel precipitation is associated with the tropical South Atlantic/tropical North Atlantic thermal gradient, which pulls the monsoon cell southward. However, the cross-equatorial Atlantic gradient does not strengthen with time, and is therefore not able to explain the magnitude of the decrease in western Sahel precipitation over the twenty-first century in CESM1. We propose that the global warming induces a zonal contrast in precipitation due to the competitive effects of the warming over land, that is due to an increase in greenhouse gases concentration, and that is leading to an increase in precipitation over the central Sahel, and to the warming of the SST, which is associated with a southward shift of the monsoon circulation, that is par-ticularly strong over the western Sahel. This is consistent with Gaetani et al. (2017), which have shown that the Sahel precipitation zonal contrast is contributed by both increases in SST and CO 2 concentration. We show that changes in SST contrasts (i.e., the TNA/ TSA thermal gradient and ENSO) do not play a prominent role on the zonal precipitation contrast over the twenty-first century. The Saharan Heat Low is often pointed out to explaining the East/West difference in Sahel precipitation change (Biasutti et al. 2009;Roehrig et al. 2011;Lavaysse et al. 2016). However, we only find a relationship between the location of the Heat Low and central Sahel precipitation.
In summary, we found that the Sahel precipitation zonal contrast in CESM1 is associated with two mechanisms: (1) the warming of northern Hemisphere is associated with a northward shift of the monsoon and an increase in central Sahel precipitation; (2) the warming of the SSTs is associated with a southward shift of the circulation.
We speculate that results could be strengthened by local feedbacks, through a modification of the location and strength of the Sahara Heat Low and, of the content in soil moisture through land-atmosphere feedbacks (Koster et al. 2004;Wang and Alo 2012;Mariotti et al. 2014).
A caveat of this study relies on the limited ability of CESM1 to simulate the observed Sahel precipitation variability (Monerie et al. 2017) and the observed SST-ARPI relationship (Fig. S13). Further work is then needed to better document changes in observed ARPI. A multi-model analysis should also be performed to assess the relationship between Atlantic tropical SSTs with the West African monsoon, as done with the Atlantic Multidecadal Variability (Martin and Thorncroft 2014a), the Tropical Easterly Jet (Whittleston et al. 2017) and the Heat Low (Dixon et al. 2016). Another caveat of the study is that we used monthly means while start and demise dates of the monsoon are projected to be delayed by a few days over Sahel (Dunning et al. 2018). Results might be slightly different using simulated rainy season period, especially when accounting for seasonal changes.