Impacts of the Atlantic Equatorial Mode in a warmer climate

The main source of sea surface temperature (SST) variability in the Tropical Atlantic at interannual time scales is the Equatorial Mode or Atlantic El Niño. It has been shown to affect the adjacent continents and also remote regions, leading to a weakened Indian Monsoon and promoting La Niña-type anomalies over the Pacific. However, its effects in a warmer climate are unknown. This work analyses the impact of the Equatorial Mode at the end of the twenty first century by means of sensitivity experiments with an atmosphere general circulation model. The prescribed boundary conditions for the future climate are based on the outputs from models participating in the coupled model intercomparison project-phase V. Our results suggest that even if the characteristics of the Equatorial Mode at the end of the twenty first century remained equal to those of the twentieth century, there will be an eastward shift of the main rainfall positive anomalies in the Tropical Atlantic and a weakening of the negative rainfall anomalies over the Asian monsoon due to the change in climatological SSTs. We also show that extratropical surface temperature anomalies over land related to the mode will change in regions like Southwestern Europe, East Australia, Asia or North America due to the eastward shift of the sea level pressure systems and related surface winds.

section_content: Interannual variability in the tropical Atlantic basin is dominated by the Equatorial Mode (Zebiak 1993). Its signature on sea surface temperature (SST) anomalies is a warming of the tropical eastern Atlantic basin related to a deepening of the thermocline (Carton and Huang 1994; Carton et al. 1996). The mode starts in Angola/Benguela region and propagates west and north via Rossby waves (Polo et al. 2008). 

The Equatorial Mode is also known as the Atlantic El Niño due to its similarities with its Pacific counterpart. The underlying mechanism for this mode is the Bjerknes feedback (Chang et al. 2006; Keenlyside and Latif 2007; Jansen et al. 2009), by which a warm anomaly in the eastern equatorial basin weakens the trade winds, which, in turn, reduce upwelling and deepen the thermocline in the east causing further warming (Bjerknes 1969). Both El Niño events are well described by the recharge oscillator model (Jansen et al. 2009 ). However, the Atlantic Equatorial Mode is further damped than the Pacific El Niño (Zebiak 1993; Lübbecke and McPhaden 2013). The Atlantic Equatorial Mode is tightly tied to the seasonal cycle but, unlike the Pacific El Niño, it shows maximum SST anomalies in boreal summer. 

The Atlantic Equatorial Mode is known to affect interannual climate variability over adjacent and remote regions. A positive phase of the mode leads to increases of rainfall over the Gulf of Guinea and reductions over the Sahel (Janicot et al. 1998; Vizy and Cook 2002; Giannini et al. 2005; Losada et al. 2010a; Mohino et al. 2011, among others). This dipole of rainfall over West Africa is established through a reduction of the surface temperature gradient that weakens the monsoon flow and shifts it southwards (Janicot et al. 1998; Losada et al. 2010a ). On the other side of the basin, rainfall over the brazilian Nordeste region is increased by the Equatorial Mode through a southward shift in the location of the Intertropical Convergence Zone (Giannini et al. 2004). Garcia-Serrano et al. (2011) also report impacts over central Europe due to poleward Rossby wave propagation. In turn, Losada et al. (2012) suggest that it can also affect Mediterranean rainfall during boreal summer, although differently depending on the years taken into account. 

Away from the Atlantic basin, Kucharski et al. (2009) and Losada et al. (2010b) argue that the Equatorial Mode promotes anomalous rising over the Equatorial Atlantic and anomalous subsidence over the rest of the tropics. This can lead to a weakening of the Asian monsoon (Kucharski et al. 2008; Losada et al. 2010b ) and also explain the reduction in the anticorrelation between ENSO and the Indian monsoon rainfall (Kucharski et al. 2007 ). On the Pacific basin the anomalous subsidence has been shown to favour La Niña development in the last part of the twentieth century (Rodriguez-Fonseca et al. 2009; Ding et al. 2012), which can have impacts for ENSO prediction (Keenlyside et al. 2013). However, to our knowledge, there is no study on the possible changes of these impacts in a future warmer climate, which is the main aim of this work. 

The usual tool to explore climate change for the next century are general circulation models (GCMs), which couple several components of the climate system, including the atmosphere and the ocean. These GCMs respond to prescribed time-varying concentrations of various atmospheric constituents like greenhouse gases (Taylor et al. 2011). However, GCMs are known to have serious systematic biases which are common to many models (Wang et al. 2014). Ashfaq et al. (2011) show that these SST biases can substantially affect the simulated tropical rainfall response to climate change. One of the regions with severe biases is the tropical southeastern Atlantic, which shows a warm bias with a magnitude comparable to the projected climate change of SSTs in the area (Fig. 1 ). This bias has been shown to hinder the predictability skill in the region (Stockdale et al. 2006). Richter et al. (2014) showed that there has been very modest improvement in reducing this bias in the last years. In both coupled model intercomparison projects phase 3 (CMIP3) and phase 5 (CMIP5) most models show a zonal gradient of equatorial SSTs in boreal summer that is reversed with respect to the observations (Richter et al. 2014). The bias seems related to weaker than observed easterlies in the previous boreal spring, which deepen the thermocline in the east and lead to its shoaling in the west (Richter et al. 2012 (Richter et al., 2014)). Part of the weak easterlies seem to come from the atmospheric part of the model through less rainfall than observed over South America and more rainfall than observed over Equatorial Africa and a southward location of the ITCZ (Richter and Xie 2008; Richter et al. 2012 ). In addition, very few current GCMs are able to simulate a realistic Atlantic Equatorial Mode (Richter et al. 2014 ) and coupled models fail to simulate many of the observed teleconnections of the Tropical Atlantic SSTs, possibly due to model climatological biases, especially over the Tropical Atlantic (Barimalala et al. 2012). 

In this work we aim at diagnosing the possible changes in the impacts of the Atlantic Equatorial Mode in a warmer climate. On the one hand, these changes could be due to modifications in the characteristics of the mode. On the other, changes in the background state of the climate system could alter the impacts of the mode through the modification of climatological features relevant for the teleconnections with remote regions. The presence of severe biases in the Tropical Atlantic and the difficulty for GCMs to simulate a realistic Equatorial Mode and its impact hinder the study of changes in the characteristics of the mode. However, GCM projected SST changes can be used to evaluate the change of the impact of the mode due to a modification of the background state. This is the approach we use in this work. We analyze the change of the impact of the Atlantic Equatorial Mode by means of sensitivity experiments performed with the SPEEDY Atmospheric General Circulation Model (AGCM) that use as boundary conditions SSTs 2 ). b Average summer (July to September) change in SST (°C, shaded) and sea ice concentration (hatched) between the period 2071-2100 and the period 1971-2000 simulated by the CMIP5 models used. For sea ice concentration regions with a reduction between 15 and 30 points are simply hatched, while those with a reduction above 30 points are double hatched from the HadISST1 reconstruction (Rayner et al. 2003) and projected SST changes given by 11 GCMs participating in CMIP5 (Taylor et al. 2011 ). 

The study is organized in four sections. In Sect. 2 we describe the data and model used, as well as the experimental design followed. The main results of the work are presented in Sect. 3. Finally, Sect. 4 provides a summary of the results and the main conclusions of the study.
section_name: Introduction

section_name: Data and methods

section_content: HadISST1 (Rayner et al. 2003) dataset is used in this work. It is a reconstructed dataset of SST and sea ice concentration with a horizontal resolution of 1° in latitude and longitude and global coverage. It spans the period 1870 to present with a time resolution of 1 month. We also use monthly mean precipitation estimates over land from the University of East Anglia Climate Research Unit (CRU), version TS3. 1 (Mitchel and Jones 2005).
section_name: Observations

section_content: Associated with the fifth report of the Intergovernmental Panel on Climate Change (IPCC 2013), there is a coordinated project known as the fifth phase of the Coupled Model Intercomparison Project (CMIP5) that involves some of the most important modelling groups in the world. This project provides a framework for coordinated modelling experiments with state-of-the-art general circulation models (GCM) (Taylor et al. 2011). Two sets of long-term simulations were used from CMIP5: historical and future projection RCP8. 5. The former is forced by observed atmospheric changes and typically covers the period from the mid-nineteenth century to present. The RCP8. 5 future projection used is representative of a high emission scenario, with an approximative radiative forcing of 8. 5 Wm -2 at the end of the twenty first century. It prescribes a CO 2 and CH 4 concentration of close to three times greater (2. 9 CO 2 and 2. 7 CH 4 ) at the end of the twenty first century than at the end of the twentieth. 

In this work we use the SPEEDY AGCM model version 40 run with a spectral truncation at a total wavenumber 30 (approximate resolution of 3. 8° in latitude and longitude) and eight vertical levels. Details for this model can be found in Molteni (2003) and Kucharski et al. (2007). Despite some biases in the representation of the summer Asian monsoon, the model presents a realistic distribution of boreal summer rainfall over the East Pacific and the Tropical Atlantic (Molteni 2003). The model has been extensively used in climate studies regarding the teleconnections of the Tropical Atlantic (e. g. Martín-Rey et al. 2012; Barimalala et al. 2012; Kucharski et al. 2007 Kucharski et al., 2008 Kucharski et al., 2009 Kucharski et al., 2011;; Losada et al. 2012; Rodriguez-Fonseca et al. 2009). Rauscher et al. (2011) and Herceg Bulic et al. (2012) also argue that SPEEDY can be used in climate change impact studies as it is capable of simulating similar anomalies than those from complex coupled models given the appropriate boundary conditions.
section_name: Models

section_name: Experimental design

section_content: In this work we analyze the impact of the Equatorial Mode in two different periods: the period 1971-2000, which we name hereafter the "present climate" and the period 2071-2100, which we name hereafter "future climate". To simulate the present climate, we use SSTs and sea ice concentrations averaged in the period 1971-2000 from the HadISST1 database (Rayner et al. 2003). This is labeled as 1971-2000 in Table 1. In the present climate simulations, the absorption coefficient in the CO 2 band is fixed to be consistent with late-twentieth century CO 2 concentrations. 

To simulate the future climate we rely on simulations performed by some of the models participating in CMIP5 historical and RCP8. 5 experiments. However, Ashfaq et al. (2011) showed that GCM SST biases can greatly affect projections of climate change. To avoid the effect of these biases we use a similar approach to the one proposed in Patricola and Cook (2011), Rauscher et al. (2011) and Chung et al. (2014) : for both, SST and sea ice concentration, we estimate the future climatologies as the addition of the changes simulated by the models in each variable to the present climatologies derived from HadISST1 dataset. This is labeled as 2071-2100 in Table 1. The SST and sea ice concentration changes are calculated as the multi-model average difference between the climatologies simulated by each model for the present and future climates. Figure 1b illustrates the average change in climatological summer SSTs simulated by the CMIP5 models. For consistency reasons, the absorption coefficient in the CO 2 band in the future climate simulations is fixed as 2. 9 times the one used in the present climate simulations. In total, 11 models were selected based on the availability of their historical and RCP8. 5 simulations at the time of analysis (Table 2 ). We also discarded different versions of a same model.
section_name: Climatology

section_content: To define the SST anomalies corresponding to the Equatorial Mode we calculate the empirical orthogonal functions (EOFs) of the June-July-August-September (JJAS) anomalies of SST for the Atlantic ocean covering the latitudes 30°S to 20°N and for the period 1971-2000. SST data is detrended previously to the EOF calculation. Then, we select those years in which the principal component (PC)

of the first EOF is larger (smaller) than one (-1) standard deviation for the 4 months. This gives us 1973, 1984, 1987, 1998 as positive years and 1976, 1982, 1992 as negative years. We calculate the anomalous pattern of the Equatorial Mode for the model experiments as the difference between the positive and negative years. Finally, to avoid discontinuities in the SST fields, we apply a linear decrease of the SST anomalies in a 10° band polewards from 30°S and 20°N. 

In Fig. 2a we show the summer (July to September) mean SST anomalies corresponding to the Equatorial Mode that we use in the AGCM experiments. The warm anomalies cover the whole equatorial Atlantic from March to December, with the maximum warming during boreal summer, especially during the months of July and August (not shown). The analysis presented in this paper focuses on the boreal summer season, defined hereafter as the average in July, August and September.
section_name: Definition of the Atlantic Equatorial Mode

section_content: Table 1 summarizes the simulations performed to evaluate the impact of the Equatorial Mode in the present and future climates. Each simulation consists of a 60-year run started from an atmosphere at rest. The first 10 years of simulation are discarded. We perform two sets of simulations, one for the present climate and another one for the future one. Each set is composed by three simulations, one for the positive phase of the mode, one for the negative and one for the 1 ). In this work we focus on the linear part of the impact of the Equatorial Mode and estimate it as the difference between the positive and negative simulations averaged over the last 50 years of the simulations for a given set. To test whether this difference is statistically significant, a two-tailed t test is used at a significance level of alpha 0. 05. The main conclusions of the work hold when defining the impact in terms of the positive minus control runs or the negative minus control ones (not shown). 

To evaluate the change of impact, we subtract the difference between the positive and negative simulations in the present climate set to the same difference obtained in the future climate set. To test the significance of this change, we use the Wilcoxon-Mann-Whitney rank-sum test. This is a non parametric test that evaluates the difference in location of two independent samples. Details of the test can be found in Wilks (2005). Again, we use a significance level of alpha 0. 05.
section_name: AGCM experiments

section_name: Results

section_content: We begin our analysis by focusing on the tropical response to the Equatorial Mode. In Fig. 3 we show the impact of the mode under present conditions, future conditions and the change between these impacts. We focus on the boreal summer season (July to September), when the amplitude of the Equatorial Mode is highest. 

Under the present climate, the model shows a rainfall response to the mode consistent with previous studies. In the Tropical Atlantic a dipole of rainfall is present with an anomalous increase below 10°N and an anomalous decrease above this latitude. Over the surrounding continents, an increase of rainfall over the Guinea Gulf and Northern Brazil is simulated, while significant rainfall deficits are present over East Africa, in accordance with Regions where differences are statistically significant (alpha = 0. 05) are marked with crosses Giannini et al. (2004) and Losada et al. (2010a). Over Central America, rainfall is enhanced while it is reduced over the Florida Peninsula. Such anomalies are consistent with the precipitation response shown by Wang et al. (2008) in their study on the atmospheric impact of the Atlantic Warm Pool and could, thus, be related to the anomalous warming of the southern part of the Atlantic Warm Pool that accompanies the Equatorial Mode in our simulations (Fig. 2a ). In the Indian basin there are negative rainfall anomalies over the Asian monsoon and positive ones over the Equatorial Indian Ocean (Fig. 3a ). These anomalies are related together by a Hadley-type cell (not shown), as suggested by Losada et al. (2010b). Over the Tropical Pacific, there is reduced rainfall over the dateline in connection to a reinforcement of the trades (not shown), which is consistent with previous studies (Rodriguez-Fonseca et al. 2009; Losada et al. 2010b ). 

The simulated tropical rainfall anomalies over land show consistency with observations, especially over Africa and Central and South America (Figs. 2b, 3a). However, the model does not simulate the positive anomalies of rainfall over the Maritime continent shown in observations. Such discrepancy is related to the SST anomalies over the Pacific basin that were not used in our experiments (Fig. 2a ) and which reduce rainfall over the central Pacific and enhance it over the Maritime Continent (not shown). 

In the future climate the impact of the Equatorial Mode is increased in the Tropical Atlantic basin, with enhanced rainfall over the Amazonian basin and over a zonal band between the Equator and 10°N spanning the Atlantic ocean east of 35°W, and West Africa (Fig. 3b, c ). This is consistent with increased convergence at low-levels in those areas (Fig. 4c ) and a more negative anomaly of the velocity potential at 200 hPa over the Tropical North Atlantic off the South America coast (not shown). The maximum anomalous rainfall and low-level divergence are shifted eastwards from approximately 40°W in the present climate to 30°W in the future one (Figs. 3, 4 ). This change in the impact of the Equatorial Mode can be related to the increase in SSTs over the Tropical Atlantic from the present to the future. In Fig. 5a, b we show that the area where the SST approximate threshold of 28 °C for deep convection (Zhang 1993; Sud et al. 1999 ) is reached is amplified in the Tropical Atlantic basin in a future climate, including much of its eastern part (Fig. 5b ). This change of climatological SSTs could promote the increase and shift towards the north and east of the climatological low-level convergence (Fig. 4d, e ) and associated rainfall over the North Tropical Atlantic (not shown). In this future scenario, an additional SST anomaly over the Tropical Atlantic related to the Atlantic Equatorial Mode has an increased local impact that is shifted to the east of the basin. 

However, the main differences in rainfall between the present and future simulations are found in the tropical Indian Ocean. The negative anomalies of rainfall located over the Bay of Bengal and the Indian Continent weaken, while those located over the western Indian Ocean are reinforced (Fig. 3 ). This change contradicts the results of previous works (Kucharski et al. 2008; Losada et al. 2010b ); following those previous results one would expect that an increase in low level convergence over the Tropical Atlantic and West Africa would lead to an increase in low level divergence over the Indian Peninsula and Bay of Bengal. Nevertheless, in our results the impact of the Equatorial Mode in a future climate shows a weakening in the low level divergence over that region (Fig. 4c ). The explanation for this apparent inconsistency can be related to the change in climatological divergence (Fig. 4d, e ): In the present simulation there is a region of low level convergence over the Bay of Bengal related to the Asian Monsoon. In the future climatology this convergence is weakened, while the convergence over the north of India is intensified. To explain such a change in the climatological divergence, we look into the characteristics of the climatological SSTs over the Indian Ocean, for the present and future climates. In addition to a general increase in the temperature of the whole basin, the main gradients of SST do not change in the western part of the Ocean, while they do in its eastern part (Fig. 6 ). In the present conditions, there is a relative maximum (minimum) of SST north (south) of 10°N over the Bay of Bengal (Fig. 6a ). Such maximum would be responsible for the local maximum of low level climatological convergence (Fig. 4d ), which is also related to a maximum of climatological precipitation there (not shown). Conversely, in the future SSTs such maximum disappears (Fig. 6b ), and so does the local maximum of climatological low level convergence (Fig. 4e ) and precipitation (not shown). That is, due to changes in the climatological SSTs, in the future simulations the main centers of climatological divergence and precipitation in the Indian basin are shifted to the west. Because the impact of any remote source of anomalous divergence would be directed to the regions with maximum climatological divergences, the main subsidence over the Indian Ocean related to the Atlantic Equatorial Mode in the future climate is shifted to the west of the basin and the observed link between the Indian Monsoon and the Equatorial Atlantic mode will be weakened.
section_name: Tropics

section_content: In response to the Equatorial Mode, the model simulates rainfall increases over the continental areas surrounding the North Atlantic basin centered around 40°N (Fig. 3a ). 

In particular, there are significant increases in the present climate over the eastern North America, and over Western Europe, above the north coast of the Mediterranean Sea, in correspondence with Losada et al. (2012). The positive phase of the Equatorial Mode also enhances rainfall over the South Atlantic, in a wide area centered around 35°S and 10°W in the present climate. These South Atlantic rainfall anomalies shift approximately 10° to the east in a warmer climate (Fig. 3b ). Such shift could be connected to the eastward shift of the maximum local rainfall anomalies over the Tropical Atlantic Ocean that were discussed previously. 

In the following, we show that this eastward shift in positive rainfall anomalies (and related divergence) over the equatorial Atlantic has consequences for the teleconnections of the Equatorial Mode to the extratropics. 

In addition to rainfall anomalies, the Atlantic Equatorial Mode impacts extratropical continental surface temperatures. In Fig. 5a we show the anomalous surface temperatures associated with a positive Atlantic Equatorial Mode event in the present climate. There are regions of increased surface temperature, like the Western North America, Central southern Europe, South Africa and a broad region that spans from Iran to Kazakhstan, and also regions of negative anomalies, like central Canada to the south of the Hudson Bay, Mongolia and Southern East Russia and the western Australia. The advection of mean temperatures by the anomalous meridional wind can explain most of these surface anomalies (Fig. 5a ). These extratropical anomalous winds are, in turn, related to the anomalous sea level pressure (Fig. 5d ), which shows a correspondence with anomalies of geopotential height at high levels (Fig. 7a ), suggesting an extratropical barotropic response to the Atlantic Equatorial Mode. This could be related to stationary Rossby waves, as proposed by García-Serrano et al. ( 2011 ), who noted that Rossby waves could propagate to central Europe in response to the mature phase of the Atlantic Equatorial Mode. 

To further study the extratropical teleconnections due to stationary Rossby Waves we show in Fig. 7b the stationary total wave number (K s ) calculated from the climatological zonal wind at 200 hPa as:

where β M is cos (latitude) times the meridional gradient of the absolute vorticy on the sphere, R t is the Earth's radius and ν is the climatological mean of the relative rotational rate of the atmosphere (Hoskins and Karoly 1981; Hoskins and Ambrizzi 1993). Rossby waves can propagate provided their zonal wavenumber is smaller or equal to K s and local maxima of K s provides Rossby waveguides. In the Southern Hemisphere there are two clear waveguides centered around 30°S and 50°S, consistent with findings of Ambrizzi et al. (1995) for the boreal summer atmosphere. The anomalous cyclones and anticyclones found in the Southern Hemisphere are thus consistent with stationary Rossby waves trapped into these two waveguides. In the North Hemisphere there is also a waveguide related to the North Atlantic jet and to the North African-Asian jet with a total wavenumber of approximately five, which is consistent with the apparent wavenumber-5 structure in Fig. 7a. 

(1) To look for the sources of these Rossby waves the anomalous rossby wave sources (RWS) have been calculated using Eq. ( 6 ) in the appendix and are shown in Fig. 7d. Following Rodwell and Jung (2008) we have chosen to integrate vertically the Rossby wave sources between 300 and 100 hPa to avoid biasing the diagnostic due to changes in the level of convective outflow in the tropics. In our case, the second term in Eq. ( 6 ), related to the anomalous divergence, dominates the anomalous Rossby wave sources. This can be seen when comparing the anomalous Rossby wave sources in Fig. 7d with the anomalous divergence in those same levels (Fig. 7c ). There is a clear correspondence between the main anomalies of divergence and the total anomalous RWS. In addition to the tropical sources, there are extratropical Rossby wave sources (Fig. 7d ). The most notable of them is the positive source located in the south Atlantic, with a center in approximately 40°S, 10°W. The anomalous divergence associated with this RWS could be related to negative divergence anomalies at approximately 15°S through a meridional overturning circulation (Fig. 7c ); the anomalous convergence in 15°S, in turn, would be related to the main tropical divergence at roughly ), e first (--→ ν ψ ′∇ξ ) and f second (--→ ν ψ ∇ξ ′ ) response terms (10 -11 s -2 ) to the anomalous RWS averaged over 100-300 hPa pressure levels. (See details in the text). For (d, e and f) contour intervals are marked at ±1, ±2, ±4, ±6, ±8, ±10, ±12 and ±14. Negative values are marked with dashed lines 10°N through a Hadley-type anomalous cell. This extratropical RWS is responsible for much of the Rossby wave response seen in the South Atlantic and Indian Oceans (Figs. 5d, 7a), which is more clearly seen when looking at the stationary response to the RWS. 

To analyze the stationary response to the RWS we show in Fig. 7e, f the last two terms in Eq. ( 7 ) in the appendix, which represent the advection of mean vorticity by the anomalous rotational wind and the advection of anomalous vorticity by the mean rotational wind, respectively. In both of them, the propagation of Rossby waves at three different latitudes can be clearly seen: one in the northern Hemisphere, between 30° and 60° N, and two in the Southern Hemisphere, centered at 30°S and 50°S, respectively. These correspond very well to the waveguides highlighted in Fig. 7b that were previously commented. There is a nice correspondence between both response terms in the south Indian Ocean over the most southern waveguide, where the RWS is very small (Fig. 7d-f ). This model response can, thus, be understood as a Rossby wave being forced by the anomalous RWS in the South Atlantic, which is related to the local anomalous divergence (Fig. 7c, d ). Such anomalous RWS is also in part responsible for the Rossby waves propagating at 30°S, though the RWS at other locations in this waveguide seem to be also contribution to them. In particular, the ones related to anomalous divergence at upper levels over South Africa, and to anomalous convergence over the south Indian Ocean off coast of Australia and over the South Pacific. In the Northern Hemisphere the response terms in Fig. 7e, f seem to be excited by different RWS, as the ones over south France, the Caspian sea, North America below the great Lakes and over the North Atlantic offshore of the Gulf of St. Lawrence, which are all related to anomalous upper level divergence (Fig. 7c ) and statistically significant positive rainfall anomalies (Fig. 3a ). These extratropical RWS are ultimately forced by the tropical forcing. 

There is also a correspondence between the response terms in Fig. 7e, f and the geopotential anomalies in upper levels in Fig. 7a. The anomalous rotational winds related to an anomalous low are directed southward to the west of the low and northward to the east in the Northern Hemisphere. Taking into account the poleward gradient of absolute mean vorticity and the negative sign in the definition of the second term in Eq. ( 7 ) in the appendix, the centers of the extratropical anomalous lows in the geopotential height are expected to be located over the regions where there is a zonal change between positive and negative values of this second term in Eq. ( 7 ) in the appendix (advection of mean vorticity by the anomalous rotational wind) in the Northern Hemisphere, and the reverse in the Southern Hemisphere, which can be seen when comparing Fig. 7a, e. 

Regarding the future climate, the analysis of the stationary total wave number K s (Fig. 8b ) reveals that the location of the waveguides in the southern and northern hemisphere does not change with respect to the present simulations. Thus, we expect that the changes in the extratropical Rossby waves in the future simulations will be due to the changes in the response to the Equatorial Mode, and not to changes in the mean flow. 

As we explained in the previous section, in a warmer climate the impact of the Atlantic Equatorial Mode changes in the tropics. There is a shift to the east in the rainfall maximum anomalies in the equatorial Atlantic related to a shift in the convergence at low levels (Figs. 3, 4). At high levels the maximum in tropical divergence is also shifted to the east, and there is an increased response over West Africa (Fig. 8c ). This eastward shift in the tropical divergence over the Atlantic is accompanied by a shift to the east in the pole of extratropical divergence located over the South Atlantic of approximately 10° and the reduction of the divergence pole over the North Atlantic, offshore of the Gulf of St. Lawrence. The divergence at high levels over southern Europe is also shifted to the east from the south of France in the present simulations to the north of the Adriatic sea in the future ones. Over the Indian Ocean anomalous rainfall and divergence at low and high levels change in the tropics: in a warmer climate the impact of the Atlantic Equatorial Mode in rainfall is shifted to the west of the Indian basin and it is greatly reduced (Fig. 3 ). Consistently, the extratropical divergence at high levels is also reduced over the South Indian Ocean and over the Asian continent (Fig. 8c ). 

These changes in the divergent field have their implication in the RWS (Fig. 8d ): the main RWS in the South Atlantic is displaced eastwards to the Greenwich meridian in the future simulation. This, in turn, displaces to the east the other Rossby wave sources in the southern Hemisphere, like the one in the eastern South Pacific. Some of the poles of RWS in the Northern Hemisphere are also shifted to the east, as the ones over the North Atlantic and over South Europe. The former one is also weakened in a future climate in correspondence with a weaker divergence at high levels (Fig. 8c ). In addition, most of the RWS over Asia are weakened. 

These changes in the location and strength of the RWS modify the location and strength of the response and, thus, of the extratropical geopotential height anomalies (Fig. 8a, e, f ). The displacement to the east of the stationary Rossby waves in the southern Hemisphere is also seen in the mean sea level pressure (Fig. 5e ). In agreement with the displacement of the Rossby wave response over Australia (Figs. 7e, f, 8e, f), the negative anomaly in geopotential height at 200 hPa located over east Australia in the present simulation is displaced to the dateline in the future simulation (Figs. 7a, 8a). This structure is barotropic and advects cold surface mean temperatures through the anomalous wind in the present simulation (Fig. 5a ). The eastward shift in the future simulation reduces meridional wind anomalies and, thus, surface temperature anomalies. For this reason, over Australia the impact of the Atlantic Equatorial Mode on surface temperatures significantly changes in the future (Fig. 5b ). Following the shift of the local RWS over Southwestern Europe, the eastward shift of the cyclonic barotropic anomaly from the present to the future simulation also implies a change in surface anomalous winds over the western Mediterranean from southerly in the present climate (Fig. 5a ) to northerly in the future climate (Fig. 5b ). This, in turn, changes anomalous surface temperatures over Western Europe from warm anomalies in the present simulation to cold ones in the future simulation. Over North America, the eastward shift and weakening of the barotropic cyclonic anomaly located over the Labrador Peninsula (Figs. 7a, 8a, 5d, e ) weakens the anomalous northerly winds at the surface over the Great Lakes and thus the cold surface temperature anomalies in the future simulation (Fig. 5a, b ). The weakening of the RWS over Asia (Fig. 8d ) is translated into a weakening of the Rossby wave response in the region (Fig. 8e ) and, thus, of the anomalies in geopotential height (Figs. 8a, 5e ). At the surface, the northerly surface wind anomalies associated with the barotropic cyclone located over the Caspian sea in the present simulation are weakened in the future one (Fig. 5f ) and so is the temperature surface anomaly (Fig. 5b ). 

In summary, the eastward shift in the Tropical Atlantic divergence and the weakening of the convergence over the tropical Indian Ocean in the future simulation with respect to the present one lead to an eastward shift of the divergence over the extratropical Atlantic and a weakening of the divergent anomalies over the Indian Ocean, respectively. These changes are connected through the RWS to the strength and location of the stationary Rossby Wave response which marks the location of anomalies in geopotential height in the extratropics. These extratropical geopotential height anomalies can, in turn, explain the surface temperature anomalies in the extratropics through the anomalous advection of mean surface temperature. In this way, in a future climate the Atlantic Equatorial Mode changes the impact in surface temperatures over regions like Southwestern Europe, East Australia, Asia or North America. Previous works have posed that the tropical diabatic heating in the Tropical Atlantic (that could be, in turn, produced by anomalous SSTs associated to the Equatorial Mode) is related to temperature conditions and summer heat waves in western Europe, through the excitation of extratropical Rossby waves (Cassou et al. 2005; Della-Marta et al. 2007 ). Our results imply a possible change in the Rossby wave response to an equal warming in the Tropical Atlantic, and thus a change of the impact in the European heat waves in summer.
section_name: Extratropics

section_content: Up to now we have analyzed the changes in the impact of the Atlantic Equatorial Mode in a future warmer climate. We have fixed the mode and we have changed the climatological SSTs, sea ice concentration and CO 2 concentration in the simulations. In this section we want to draw light into the origin of the changes in the impacts. For this sake we have performed additional simulations changing the prescribed CO 2 concentration (FutuCO2 simulations), the climatological sea ice concentration (futuSIC simulations) or the SST fields (futuSST simulations), leaving the other two factors as in the present climate. As for the future and present simulations, each set of simulations is composed by three simulations (positive and negative phase of the Equatorial Mode and climatology). 

A description of this new sets is displayed in Table 3. In Fig. 9 we show the comparison between each additional set of simulations and the present ones by means of the changes in the impact of the Equatorial Mode. The change of sea ice concentration to future conditions does not alter significantly the impact of the Atlantic Equatorial Mode in rainfall (Fig. 9a ). It influences the response of surface temperatures in some extratropical locations as northwestern Russia (Fig. 9f ) by modifying the position of mean sea level pressure anomalies in the Northern Hemisphere (not shown). However, these anomalies do not agree with the ones obtained when comparing the future simulations with the present ones (Fig. 5c ). 

Though altering only CO 2 concentration cannot account for the full response of future simulations, it has a clear impact on tropical rainfall (Fig. 9b ). It can explain the increase of rainfall related to the Equatorial Mode over northern South America and contributes to rainfall changes over the equatorial Atlantic, western equatorial Africa and South Asia and eastern Indian Ocean. Conversely, the changes in surface temperature related to the modification of CO 2 are weak (Fig. 9g ). The changes with respect to the present simulations obtained with futuSST are very similar to the ones obtained in the future simulation (Fig. 9c, h ). This suggests that it is the change in climatological SSTs which matter the most for the changes in the impacts of the Equatorial Mode on rainfall and surface temperatures that we have highlighted in the previous sections. The change in climatological SSTs show a world wide increase of temperatures (Fig. 1b ). The average (70°S to 70°N) difference in SST from the present to the future CMIP5 simulations is 2. 8 °C. However, some regions warm more than others. For instance, the tropical oceans, notably the eastern Pacific, eastern Atlantic and western Indian Oceans warm further than the surrounding areas. In particular, the westward enhancement of climatological SSTs over the tropical Indian Ocean (Figs. 1b, 6 ) could be responsible for the change in climatological divergence at low levels in this basin (Fig. 4c, d ), which, in turn, seems to be the reason for the change in the impact of the Atlantic Equatorial Mode in this basin (Fig. 3 ). Chiang and Lintner (2005) argue that more important than the overall anomaly in SST is its gradient for establishing the tropical rainfall anomalies. To evaluate if the main changes in the impact of the Equatorial Mode are related to the overall warming or to the change in local SST gradients (due to the uneven tropical warming, Fig. 1b ) we have performed two sets of simulations, FutuGrads and FutuTemp. In the first set, FutuGrads, we have changed the SST climatological field to the future one (2071-2100, using RCP85 scenario) minus 2. 8 °C at each grid point. In this way, we retain the future SST gradients but reduce the average SST to the present conditions. In the second set, FutuTemp, we have added 2. 8 °C to each grid point in the present climatological SST field. In this way, we retain the future average increase in SSTs but leaving SST gradients as in the present simulations (see Table 3 for more details). 

The difference between the impacts of the Equatorial Mode in FutuGrads set of simulations and the present one is weak (Fig. 9d, i ). Its main effect is to contribute to the modification of rainfall anomalies over the Indian Ocean, thus confirming the results shown in Sect. 3. 1. Conversely, the results from FutuTemp set of simulations (Fig. 9e, j ) are very similar to the ones obtained with the SST future climatology (Fig. 9c, h ) and also to the future set (Figs. 3c, 5c). This is, most of the change in the impact of the Equatorial Mode from present to future conditions that was highlighted in the previous section can be explained by the average increase of 2. 8 °C in SST climatology. This result diverges from the one obtained by Chung et al. (2014) regarding the changes of the Pacific El Niño impacts from the present to a future climate. They showed that, unlike our case, the patterns of the precipitation response were more strongly influenced by the spatial structure of the change in SST climatologies. However, in our experiments it is the mean increase in climatological SSTs that drives eastward the main rainfall increase associated with the Equatorial Mode over the Tropical Atlantic, which, in turn, modifies the location of the divergence at high levels over the Atlantic and, thus, the Rossby Wave sources (not shown). These, in turn, shift eastward the location of extratropical pressure anomalies and, thus, modify surface temperature anomalies through the anomalous advection of climatological temperatures (not shown).
section_name: Origin of the changes

section_content: We have analyzed the local and remote impacts of the Equatorial Mode under present and future climate conditions. We followed an AGCM-based approach in which the prescribed SSTs were calculated as a superposition of the same anomalous pattern associated to the Equatorial Mode upon two different climatologies, representative of the present and future climates, respectively. The SST anomalous pattern is consistent with the Atlantic Equatorial Mode pattern observed during the period 1971-2000. The present climatological SSTs correspond to the observed SSTs averaged over the 1971-2000 period. For the future climate, the SSTs were calculated as the present climatological SSTs plus the changes in SSTs simulated by 11 models from the CMIP5 experiment (Table 2 ). With this approach we try to avoid errors derived from the strong SST biases present in current GCMs, especially over the Tropical Atlantic Ocean. For consistency reasons, the sea ice concentration and the absorption coefficient in the CO 2 band were changed accordingly in the future climate simulations. 

Our results show that the impact of the Equatorial Mode will change under future climate conditions in both the tropics and the extratropics. The changes will be necessarily linked to changes in the climatological SST, sea ice concentration or CO 2 absorption coefficient, as they are the only three factors that vary between present and future simulations. 

In the tropics, the main centers of anomalous divergence located over the Tropical Atlantic in response to the Equatorial Mode will be shifted northeasterly. Over the Indian Ocean, the impact of the Atlantic Equatorial Mode over the Indian Monsoon precipitation will be weakened. We show that the former is linked to the increase in the climatological SST over the Tropical Atlantic: In the future simulations the convection threshold of 28 °C is reached over the whole Tropical Atlantic basin, north of the equator and up to 15°N, approximately. Such future climatology allows tropical convection to take place over a wider area than in the present climate. In turn, the change in the response of the Indian Monsoon to the Equatorial Mode is linked to changes in the location of the low level climatological convergence and precipitation in the Bay of Bengal. We have shown that such changes can also be related to changes in the climatological SSTs over the Indian Ocean. 

In the extratropics, the anomalous response to the Equatorial Mode in both present and future climates has the form of extratropical Rossby waves that project onto surface temperature anomalies due to the anomalous advection of climatological surface temperature. However, the main centers of those waves appear shifted to the east in the future climate with respect to the present one, which, in turn, produces a change in the extratropical surface temperature response to the Equatorial Mode. We show that such a shift simulated by the AGCM is mainly related to the eastward shift of the anomalous divergence in the Tropical Atlantic. 

We have also analyzed the contribution of climatological SSTs, sea ice concentration and CO 2 to the changes in the impact of the Equatorial Mode described above. We have shown that the change in sea ice concentration is the least important factor in the AGCM response, while the direct change in CO 2 produces a bigger impact, especially on tropical precipitation. However, the dominant factor are climatological SSTs. Regarding the latter, we have investigated separately the role of the world-wide mean SST increase and the change in the SST spatial structure from the present to the future climate. Our results show that, although the changes in the SST spatial structure and concentration of CO 2 have a non-negligible impact in the response over the Indian sector, the key factor controlling the changes in the AGCM response to the Equatorial Mode is the spatial averaged increase of 2. 8 °C in climatological SSTs. (5)

Fig. 9 Changes in the impact of the Equatorial Mode in rainfall (mm/day) (left column) and land surface temperatures (°C, right column) between the simulations FutuSIC (first row), FutuCO2 (second row), FutuSST (third row), Futugrads (fourth row), Fututemp (fifth row) and the present simulations. Statistically significant regions (alpha = 0. 05) are marked by a thick gray contour. Changes are calculated as the differences between the impact of the Equatorial mode (defined as the positive minus negative simulations) in each set (FutuSIC, FutuCO2, FutuSST, Futugrads and Fututemp, see Table 3 ) minus the impact in the present set (Table 1)
section_name: Summary and conclusions

section_name: Acknowledgments

section_content: The research leading to these results has received funding from the European Union Seventh Framework Programme ( FP7/2007(FP7/ -2013 ) ) under grant agreement n° 603521. This work was also supported by the Spanish project CGL2012-38923-C02-01. We acknowledge the World Climate Research Programme 's Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 2 of this paper) for producing and making available their model output. For CMIP the U. S. Department of Energy's Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
section_name: Acknowledgments

section_name: Appendix

section_content: To analyze Rossby wave sources (RWS) and their stationary response we start with the vertical vorticity equation. Following Rodwell and Jung (2008) we use the shallow-atmosphere approximation and neglect the terms related to vertical advection, tilting and friction. In this case, the vertical vorticity equation reads:

where ξ is the vertical component of the absolute vorti- city, t is the time and ν is the horizontal wind. If we sepa- rate the wind into its divergent and rotational components (� ν = � ν χ + � ν ψ, where � ν χ = ∇χ is the divergent wind and � ν ψ = � k × ∇ψ is the rotational wind, χ is the velocity potential, ψ the streamfunction and k the local unit vertical vector), Eq. ( 2 ) can be written as:

If we move all the terms that are related to the divergent flow to the right-hand side of the equation and leave the ones dependent on the rotational flow in the left-hand side, the vertical vorticity equation can be expressed as:

The right-hand side term of Eq. ( 4 ) is known as the RWS (Sardeshmukh and Hoskins 1988). The anomalous RWS can be calculated following Qin and Robinson (1993) as:

where () and () ′ are the climatological mean and anomaly, respectively. Equation ( 5 ) can be further decomposed in four terms:

Neglecting the seasonal mean of the anomalous local tendency of vorticity, Eq. ( 4 ) can be re-formulated in terms of the anomalous RWS and their stationary response as:

The second and third terms in Eq. ( 7 ) can be regarded as the response to the anomalous RWS (Rodwell and Jung 2008). These two terms should cancel each other in those regions where the anomalous RWS is very small.
section_name: Appendix