The South Atlantic Anticyclone as a key player for the representation of the tropical Atlantic climate in coupled climate models

The key role of the South Atlantic Anticyclone (SAA) on the seasonal cycle of the tropical Atlantic is investigated with a regionally coupled atmosphere–ocean model for two different coupled domains. Both domains include the equatorial Atlantic and a large portion of the northern tropical Atlantic, but one extends southward, and the other northwestward. The SAA is simulated as internal model variability in the former, and is prescribed as external forcing in the latter. In the first case, the model shows significant warm biases in sea surface temperature (SST) in the Angola-Benguela front zone. If the SAA is externally prescribed, these biases are substantially reduced. The biases are both of oceanic and atmospheric origin, and are influenced by ocean–atmosphere interactions in coupled runs. The strong SST austral summer biases are associated with a weaker SAA, which weakens the winds over the southeastern tropical Atlantic, deepens the thermocline and prevents the local coastal upwelling of colder water. The biases in the basins interior in this season could be related to the advection and eddy transport of the coastal warm anomalies. In winter, the deeper thermocline and atmospheric fluxes are probably the main biases sources. Biases in incoming solar radiation and thus cloudiness seem to be a secondary effect only observed in austral winter. We conclude that the external prescription of the SAA south of 20°S improves the simulation of the seasonal cycle over the tropical Atlantic, revealing the fundamental role of this anticyclone in shaping the climate over this region.


Introduction
The tropical oceans are one of the key areas of weather and climate variability. Coupled atmosphere-ocean general circulation models (AOGCMs) typically show warm biases in sea surface temperature (SST) over both the tropical Pacific and Atlantic oceans, particularly in their eastern parts (Mechoso et al. 1995). Richter et al. (2012b) reported that the tropical Atlantic biases persist in most AOGCMs participating in the Coupled Model Intercomparison Project (CMIP) Phase 3 and 5 (Meehl et al. 2007;Taylor et al. 2012). Several studies have explored the causes and structure of these biases: for example, Li and Xie (2012) suggested that the biases could be both of oceanic or atmospheric origin. They proposed the existence of two Abstract The key role of the South Atlantic Anticyclone (SAA) on the seasonal cycle of the tropical Atlantic is investigated with a regionally coupled atmosphere-ocean model for two different coupled domains. Both domains include the equatorial Atlantic and a large portion of the northern tropical Atlantic, but one extends southward, and the other northwestward. The SAA is simulated as internal model variability in the former, and is prescribed as external forcing in the latter. In the first case, the model shows significant warm biases in sea surface temperature (SST) in the Angola-Benguela front zone. If the SAA is externally prescribed, these biases are substantially reduced. The biases are both of oceanic and atmospheric origin, and are 1 3 predominant types of tropical-wide errors in CMIP models. The first can be traced back to biases in atmospheric simulation of cloud cover and shows a broad meridional structure of the same sign across all basins, with models with enhanced cloud cover featuring lower biases in tropical-wide SST. The second originates from differences in the simulation of the thermocline depth, with models with a deeper thermocline typically showing a warmer equatorial cold tongue. Over the tropical Atlantic, the oceanic and atmospheric processes are strongly coupled and this coupling is a major determinant of climate variability over the adjacent tropical continents. For example, the year-toyear variability of the seasonal development of Atlantic cold tongue has strong impacts on rainfall variability over West Africa in general, particularly on the onset of the West African monsoon (Kushnir et al. 2006;Brandt et al. 2011). Thus, the warm SST biases in coupled models severely inhibits the ability of these models to correctly reproduce both the observed SST variability in the equatorial Atlantic and the dynamics governing that variability (Ding et al. 2015).
The present paper focuses on the southeastern tropical Atlantic (SETA; Fig. 1a). Various mechanisms have been proposed for the development of warm biases in this region. Some of the mechanisms are related and can be linked to a specific season. A possible cause involves the models' inability to correctly simulate the wind stress in the SETA. Colberg and Reason (2006) demonstrate the sensitivity of the Angola Benguela Frontal Zone position to the wind stress forcing as its curl controls the southeastward flow of the South Equatorial Counter Current towards the Angolan coast, and the opposing northward flow of cold water is mainly the result of geostrophic adjustment to coastal upwelling driven by the meridional wind stress. These local wind stress inaccuracies can be related to the large-scale atmospheric circulation due to the incorrect simulation of the strength of the South Atlantic Anticyclone (SAA, e.g., Lübbecke et al. 2010). Another possible source of errors arises from a poor simulation of convection over the continents adjacent to the tropical Atlantic during boreal spring: a deficient (excessive) precipitation over tropical South America (Africa) could cause an anomalous zonal sea-level pressure gradient in the equatorial Atlantic that produces weaker easterlies in the southern hemispheric tropical Atlantic Ocean. This bias is present in atmospheric GCMs (AGCM) forced by observed SST (Richter and Xie 2008;Wahl et al. 2011).
The models' excessive radiative forcing due to errors in reproducing the observed low, marine stratus cloud cover over the southeastern ocean can also lead to an anomalous warming (Tanimoto and Xie 2002;Huang et al. 2007). Other possible mechanisms that lead to the generation of warm biases in the SETA are the modification of the upper-ocean stratification and barrier layers by incorrectly simulated rainfall (e.g., Liu et al. 2009) and the ocean GCMs (OGCM) inability to correctly simulate the sharp thermocline characteristic of the equatorial and southeastern tropical Atlantic (Xu et al. 2014;Richter 2015). Additionally, an inadequate simulation of the seasonal shifts of the intertropical convergence zone (ITCZ) in timing and location contributes to errors in the modeled wind stress and low marine clouds.
As mentioned-above, the variability of the South Atlantic Anticyclone (SAA) affects the southern tropical Atlantic through the wind stress variability in the region, which plays a leading role in the generation of SST anomalies in the SETA (e.g., Lübbecke et al. 2010). The seasonal cycle of the SAA is mainly controlled by the convective heat sources over South America and southern Africa, and by the interaction of orography with the trade easterlies and the midlatitude westerlies (Rodwell and Hoskins 2001;Miyasaka and Nakamura 2010). It is also influenced by remote, cross-hemispheric effects of heat sources located in the Northern Hemisphere (Ji et al. 2014). These fundamental mechanisms lead to a notably stronger SAA in austral winter than during austral summer. In the latter season, the SAA is mainly a response to heating over the adjacent Southern Hemisphere continental land masses, while the northern hemispheric heat sources are responsible for the strengthening of the anticyclone in austral winter Lee et al. 2013;Ji et al. 2014).
According to Seager et al. (2003), air-sea interactions could provide a stronger contribution to the SAA development during the austral summer. The authors propose a positive feedback mechanism for this season, in which a strengthening of the SAA is associated with decreased SST in the east and with warmer conditions in the west due to poleward advection of warmer air in the lower troposphere, thereby enhancing the subtropical zonal SST gradient. The strengthened SST gradient reinforces in turn the anticyclonic circulation by stabilizing the atmosphere in the east and making it more unstable and susceptible to deep convection in the west. Under such conditions, ocean upwelling and advection drive the seasonal cycle of the SST asymmetry and control the strengthening of the SAA without greatly influencing their seasonal cycle. However, based on experiments with an uncoupled AGCM,  found that the regional zonal gradients of SST have only a minor impact on the SAA, probably due to the absence of air-sea interactions in their uncoupled model. Another factor that could influence this result is that their analysis focuses on the austral winter season, when the strength of the SAA is mostly determined by monsoonal heating in the northern hemisphere (Lee et al. 2013). Finally, Miyasaka and Nakamura (2010) proposed an additional mechanism, in which cool and warm anomalies over ocean and land surfaces, respectively act as sources of anticyclonic and cyclonic potential vorticity, creating along shore wind anomalies towards the equator.
The influence of horizontal resolution has also been explored with a limited area model. Seo et al. (2006) showed that resolution could play a role in the generation of SETA biases, while Patricola et al. (2012) suggested that simply increasing spatial resolution does not fully remove the biases. Moreover, Doi et al. (2012) concluded that greater resolution does not suffice to reduce the warm SST biases in global models. On the other hand, Harlaß et al. (2015) recently provided evidence that higher horizontal resolution accompanied by an enhanced vertical resolution does substantially improve the simulation of TA SST, precipitation and the representation of the upper ocean thermal structure, especially at the equator. In recent years, the use of regional coupled atmosphere ocean models has increased, as they enable a more detailed analysis of regional climate conditions and variability (e.g., Döscher et al. 2002;Mikolajewicz et al. 2005;Ratnam et al. 2008;Sein et al. 2014). Additionally to the higher horizontal and vertical resolution, regional climate models imply a sharper focus on selected domain location (Sein et al. 2015). The freedom of choice regarding the domain allows not only simulating the coupled interactions between components of the climate system, but also to either include or exclude some of the key sources of variability for the region of interest. This feature is particularly useful for studies in regions where complex feedbacks between the various components of the regional and global climate system occur (e.g., Sein et al. 2014)-like the SETA study region.
In the present paper, the high-resolution regionally coupled atmosphere ocean model ROM (Sein et al. 2015) is used to assess (1) the impact of the SAA location and strength on the simulation of the southern tropical Atlantic climate and (2) the role of atmosphere-ocean coupling in its seasonal cycle. The model encompasses a region where a limited area atmospheric model is coupled to a global ocean with uncoupled ocean dynamics outside the coupled region. This allows for a better understanding of local versus remote forcings of the SAA and to better resolve low-level cloud and deep convective dynamics, as well as air-sea interactions.
Our strategy is based on the comparison of ROM simulations that differ by the coupled domain placement. In the first setup, the coupled model domain includes all atmospheric features relevant for the regional southern tropical Atlantic dynamics, namely the ITCZ and the SAA. In the second setup, the coupled domain is shifted towards the Northern Hemisphere to explicitly exclude the core and southern part of the SAA. With this experimental design, the role of the SAA as an internal versus external forcing is assessed, as well as the importance of a correct simulation of SAA for the southern tropical Atlantic climate. Additional stand-alone simulations with the ocean and atmospheric components of our regional coupled model provide further insight on the impact of ocean-atmosphere interaction processes for the simulated SAA and SST. The models and experiments are briefly described in Sect. 2 and the results are discussed in Sects. 3-6. Finally, our conclusions are given in Sect. 7.

3 2 Experimental design
The regionally coupled model ROM (Sein et al. 2015) consists of a global ocean-sea ice model (MPIOM) coupled to an atmospheric regional model (REMO) and global terrestrial hydrology model (HD). The global ocean is divided into a coupled subdomain, where the ocean and the atmosphere interact, and an uncoupled region where the ocean model is driven by prescribed "perfect" atmospheric forcing and runs without any feedbacks. This allows including or excluding regions that potentially act as sources of forcing for the southern tropical Atlantic from the area of active coupling. Moreover, it permits to explore the impact of model resolution in a computationally affordable way. The coupling interval between the ocean and the atmosphere is 3 h. Lateral atmospheric and upper oceanic boundary conditions outside the coupled domain were prescribed using ERA40 reanalysis (Uppala et al. 2005) data. Following the Davies relaxation scheme (Davies 1983), no nudging is applied outside the eight points used as boundary conditions. The simulation period is 1958-2001. The MPIOM configuration used for all experiments features the grid poles over North America and Northeastern Africa (Fig. 1b). The horizontal resolution reaches 5 km in the north tropical Atlantic and gradually diminishes to 100 km in the southern oceans. The model has 40 vertical levels, with increasing level thickness towards the ocean bottom.
To discern the impact of the SAA on the simulation of the seasonal cycle in the tropical Atlantic, we use two different REMO domains with constant grid spacing in rotated grids (Fig. 1b). Both domains include the southern tropical Atlantic. The AFR domain has 50 km resolution and covers the whole SAA region and includes Africa, parts of the Indian Ocean, and the Mediterranean region. The NAT domain has also 50 km resolution and only includes the northern part of the SAA, but extends to the North Atlantic, the eastern tropical Pacific, and the Mediterranean Sea regions.
The impact of coupling on the atmosphere and ocean is explored with the help of additional uncoupled runs for MPIOM and AFR REMO configurations. The MPIOM standalone ocean simulation (hereafter TR04) is forced with ERA40 Reanalysis data. The list of experiments analyzed here is summarized in Table 1. The sensitivity of the model results to changes in resolution and the parameterization of large-scale precipitation has also been explored. As the results of these sensitivity experiments do not change our main conclusions, we do not show them here in detail. However, we refer to the results of these simulations where it is relevant for the discussion.
The outcome of our simulations is compared in the paper against different observational products. Mean sea level pressure (MSLP) and 10-m winds are taken from the ERA40 reanalysis. Higher resolution (25 km) wind data over the ocean from the National Aeronautics and Space Agency's (NASA) Quick Scatterometer (QuikSCAT) satellite are also considered. As an observational reference for the SST, we resort to the National Oceanic and Atmospheric Administration's (NOAA) Optimum Interpolation SST V2 dataset (OISSTV2; Reynolds et al. 2007), and also to the Global Ocean Data Assimilation System (GODAS; Behringer and Xue 2004) from the National Centers for Environmental prediction (NCEP). The latter SST data are considered in conjunction with winds from the NCEP reanalysis (Kalnay et al. 1996). Sub-surface ocean temperature (for the computation of thermocline depth) is taken from the Polar Science Center Hydrographic Climatology (PHC 3.0; Steele et al. 2001). Information on Sea Surface Height is retrieved from the NCEP Climate Forecasting System Reanalysis (CFSR, Saha et al. 2010) and, finally, incoming solar radiation at the Earth's surface data used here belong to the Surface Solar Radiation Dataset-Heliosat (SARAH; Müller et al. 2015), from the Satellite Application Facility on Climate Monitoring (CM SAF) consortium.

General features of the circulation in the southern tropical Atlantic
The general characteristics of the circulation in the southern tropical Atlantic are analyzed as represented by the different experiments (Table 1). The DJF mean sea level pressure (MSLP) for the period 1980-1999 is depicted in Fig. 2 for ERA40, together with the corresponding biases of AFR coupled and uncoupled simulations, as well as NAT coupled simulations. The AFR domains simulate an austral summer SAA with its northern, westward flowing branch systematically weakened compared to both ERA40 and the NAT domain. This leads to a cyclonic anomaly centered east of the southern tropical Atlantic, which weakens the equatorward winds in the SETA region (see Fig. 1c). This feature of the coupled AFR simulation does not depend on the parameterizations used (not shown) and is particularly clear when the AFR (Fig. 2b, d) and NAT domains ( Fig. 2c) are compared. The weaker MSLP values for the simulated SAA are associated with stronger biases in the SSTs for the AFR domain compared to NAT (Fig. 3c, d). These SST biases are generally higher in both of the coupled simulations than in the standalone TR04 run, but there is quite some similarity in structure. The small negative biases in the equatorial strip, already present in TR04 simulation, become stronger with the coupling in both setups, especially in the eastern and western parts of the domain. On the other hand, the strengthening of the cold bias north of the equator in the region of stronger variability (box Atl3 in Fig. 1c) is compensated by a warmer SST south of the equator, reducing the mean bias here. While the structure of the equatorial biases is not similar to those that are characteristic of global coupled models (see e.g., Li and Xie 2012), it does correspond to a weakening of the zonal SST gradient. The biases become more positive (negative) south (north) of the equatorial line. The SST biases in the Benguela Upwelling System in the AFR simulation corresponds to the typical bias pattern often found in AOGCMs (e.g., Richter et al. 2012a). This simulation shows the strongest warm bias along the coast of southern Africa near the Angola Benguela Front Zone, and along the Benguela Upwelling System expanding through the SETA. In the NAT domain, the SST bias in DJF is significantly reduced, as illustrated in the corresponding panel of Fig. 3. Our sensitivity studies with different atmospheric and ocean parameterizations indicate that although SETA SST bias can change, its large-scale pattern is rather robust and can be related to the errors in the simulation of the SAA (not shown).
The coupling also affects the precipitation patterns ( Fig. 4). The uncoupled REMO simulation shows positive precipitation biases both east and west along the equator, off the coasts of Brazil and Africa, extending southward to around 10°S, while negative biases are mainly located in the Gulf of Guinea north of the equator. In the coupled simulations, air-sea interactions cause anomalous convergence south of the equator and divergence at the equator and to the north (not shown). The precipitation bias is associated with northerly wind anomalies and cold SST biases at the equator and warm anomalies to the south. This suggests that a mechanism similar to the Wind-Evaporation-SST feedback like the one operating in the meridional mode may act (cf. Richter and Xie 2008). Consistently, the ITCZ is effectively displaced southward, or more precisely, a "double ITCZ" is created. While the magnitude and details of this effect depend on the set-up of the coupled model, this erroneous pattern is present in all coupled simulations.
The differences in simulated MSLP are smaller in austral winter (Fig. 5), when the SAA, together with the North  Atlantic Anticyclone, reach their maximum strength. For instance, the bias patterns in both AFR and NAT sets of simulations are more similar in the equatorial region, and the AFR simulations do not show the cyclonic bias in the SETA region identified for DJF. The pattern of the negative anomalies of MSLP to the south of the Gulf of Guinea for the coupled runs resemble those of the strong positive SST bias, and are weaker for the NAT simulations, with more moderate SST biases (cp. Fig. 6c with d). Our sensitivity studies with different parameterizations further suggest that the improvement of the MSLP biases in AFR50C is associated with a somewhat better simulation of SSTs in the Gulf of Guinea (not shown).
The SST biases can be associated with a particular feature of the atmospheric or oceanic component of our coupled system. To explore these possible connections, coupled and uncoupled simulations can be compared. A particular uncoupled simulation (ERA-40 driven) can show biases inherent to the specific model formulation or parameterizations. These biases can be reduced or enhanced by the coupling, revealing the possible mechanisms depending on how they affect the simulated SAA and SSTs.
In DJF, the uncoupled atmospheric model in the AFR domain shows MSLP and surface wind biases that broadly resemble those found in the coupled setup-compare Figs. 2b versus d and 4b versus d. This bias is in fact amplified by the coupling, which further weakens the MSLP in the SAA region. On the other hand, the coupling increases the MSLP values over land, and near the equator.
In JJA, the simulated SAA reaches its maximum strength and the ITCZ its northernmost position (Fig. 5). During this season, the atmospheric uncoupled run in the AFR domain presents a weaker than ERA40 SAA core and small MSLP biases, which are negative in the Atlantic and mostly positive over the African continent (Fig. 5b). The coupling increases the MSLP over the whole domain, reducing the negative biases over the Atlantic at the expense of the growth of positive biases over Africa and the western equatorial Atlantic (Fig. 5d). The northeasterly trade winds are underestimated both in the NAT and AFR simulations, and this likely reduces the contribution of the northern monsoons to the SAA strengthening in austral winter. In our model, this pattern is independent both from parameterizations and resolution.
In this season, the uncoupled TR04 simulation exhibits a zonally stretching cold bias in equatorial Atlantic Ocean. Ocean-atmosphere interaction causes the bias to extend to the north of North West Africa and to the south in the western equatorial Atlantic, both in AFR and NAT coupled simulations (Fig. 6). In SETA, TR04 shows warm biases that 1 3 stretch from the Benguela upwelling system up to the equator, which remain mainly close to the coast. These biases, that in the uncoupled run could be caused by errors in the prescribed momentum fluxes, are strengthened by the coupling. They do not extend as far into the basin's interior as in DJF but expand to the interior of the basin as approaching to the equator, especially in NAT.
The opposite occurs for the warm biases in the northern tropical Atlantic that extends further into the ocean interior in JJA. The SST warm biases in the northern tropical Atlantic show similar structure in NAT and AFR, but with larger values in the latter.
It is worth noting the existence of a clear equatorial westerly wind bias in the coupled runs in JJA that is not evident in the uncoupled AFR (Fig. 7). Some authors (Patricola et al. 2012;Richter et al. 2012b;Xu et al. 2014) postulate that westerly wind bias over the equatorial Atlantic contributes to the insufficient vertical temperature stratification and warm SST biases in the Benguela upwelling system. In AFR50C, precipitation shows dry biases over Africa and dry biases just north of the equator and wet biases further north (Fig. 7). The positive precipitation bias can be related to the positive SST anomaly off northern Africa. This SST bias is also present in the standalone ocean simulation TR04 and is enhanced in the coupled simulations.
The precipitation bias north of the equator in boreal winter is reflected in the northward ITCZ shift in the coupled runs (Fig. 8), which is located too far north from boreal winter through summer, and retracts southwards too early in autumn (especially in AFR50C). This suggests that the ocean-atmospheric feedback involves SST, evaporation and wind. Over the Gulf of Guinea, the biases could be related to a mechanism similar to the Bjerkness feedback: the excessive precipitation over the Gulf of Guinea is apparently the main cause of the equatorial westerly wind bias in the central and western Atlantic.

Analysis of SST and ocean circulation
The ocean response to the coupling in terms of SST anomalies is linked to changes in the regional ocean circulation and present strong seasonal variability. The role of the ocean heat transport in the TA climate has been explored e.g., by Foltz and McPhaden (2006), Joyce et al. (2004), andSeager et al. (2001). They found that the most important effect of horizontal advection is the damping of SST anomalies in the subtropical Atlantic, although Seager et al. (2001) also stresses the importance of the advection by the mean ocean circulation for the propagation of SST signals. The uncoupled TR04 simulation shows a cold bias at the equator and a moderate warm bias at the Angola Benguela Front Zone. According to Xu et al. (2014), an overshooting of the Angola Current (AC), a feature common to many ocean models, can cause this warm bias. To further analyze the interplay of the SSTs and the ocean circulation, we explore the seasonal cycles in the uncoupled TR04 and the coupled NAT50C and AFR50C configurations. The uncoupled ocean model TR04 simulates the tropical Atlantic SST seasonal cycle reasonably well (Fig. 9, upper row), with biases within half a degree for most of the tropical Atlantic for all seasons. Larger positive biases are found in the northeastern tropical Atlantic, with two local maxima in MAM at the Guinea Dome and south of Cape Blanc (cf. Fig. 1 for locations).
These positive northeastern tropical Atlantic biases weaken in JJA and almost disappear in SON, but reinforce again in DJF, which points to deficiencies in the simulation of the mixed layer depth in the area of the Guinea Dome (Doi et al. 2010). Negative biases occur near the Equator, largest during JJA and SON, while they almost vanish in DJF and MAM.
The most remarkable differences with Reynolds SST (Reynolds et al. 2007) are found on the southeastern tropical Atlantic. The larger departures are close to the Angola Benguela Front Zone with a year-round local maximum off Cape Frio (cf. Fig. 1c), indicative of a too southward position of the Angola Benguela Front in TR04, as reported for the ensemble of CMIP5 models by Xu et al. (2014). In DJF, the maximum SST bias is restricted to a small area by the African coast by Cape Frio. During MAM and JJA, the 1 3 maximum SST bias extends as a narrow coastal strip southward along the Benguela Upwelling System down to 25°S, while in SON this coastal warming weakens and the warm anomaly is shifted towards the ocean interior.
In the western part of the southern tropical Atlantic, the NAT50C SSTs are slightly colder than in TR04 but is generally warmer in the central part throughout the year. NAT50C SSTs are clearly cooler along the Equator and along the coasts of the Angola Basin in DJF, but warmer in the central-eastern equatorial Atlantic and along the Benguela Upwelling System. These differences decrease in JJA and SON, when the SST in the Gulf of Guinea is about 1 °C warmer in NAT50C than in TR04: this offsets the cold TR04 bias in the eastern equatorial ocean (5°N-5°S) but increases the warm bias in the southeastern low latitude tropical ocean (5°S-10°S). The impact of the coupling in the AFR domain is more prominent, as most of the southern tropical Atlantic shows consistently warmer SSTs than the uncoupled TR04. The warm bias exhibits the maximum westward and southward extension for DJF and SON, reaching up to 4 °C in the Benguela Upwelling System and being significant throughout the year, except MAM, in SETA. As it could be expected, differences between NAT and TR04 outside of the coupling domain (roughly south of 15°) tend to be lesser due to the bulk heat fluxes.
These SST biases lead to a modification of the zonal SST gradients in the coupled simulations, as seen in the departure of the zonal mean depicted in Fig. 10. The DJF Reynolds climatology shows a clear negative zonal SST gradient through the whole southern tropical Atlantic with the exception of the easternmost region between the Angola Benguela Front and the Gulf of Guinea [15°S-0°], which presents a weak positive zonal SST gradient.
During austral winter months (JJA), the negative zonal SST gradient is enhanced, particularly at more equatorial latitudes. As a consequence, the area of positive SST gradients shrinks to a small band along the equator within the Gulf of Guinea as a consequence of the seasonal equatorial Atlantic cold tongue appearance (Weingartner and Weisberg 1991). The standalone TR04 does not show big differences against the climatology, although the austral winter reinforcement of the SST gradient is not strong due to an insufficient cold tongue development. An evident effect of coupling is that the SST zonal gradient between the Angola Benguela Front and the Gulf of Guinea does not change sign from DJF to JJA (see Fig. 10, AFR50C). AFR also presents a weaker zonal SST gradient in the whole SETA for JJA (Fig. 10).
The implications to the surface ocean circulation are depicted in Fig. 11. As suggested by Xu et al. (2014) based on numerical results, our AFR runs also provide evidence of a clear connection between Angola Benguela Front warm bias and the strength and southward extension of the south flowing Angola current. Geostrophic currents in CFSR show the northward Benguela current that meets with the Angola current at the position of the Angola Benguela Front (between 16°S and 17°S). The northward current is in geostrophic balance with an onshore pressure gradient likely caused by coastal upwelling, which reduces the sea surface elevation locally   (Philander and Yoon 1982), while the Angola current seems to be related to cyclonic wind stress curl (Colberg and Reason 2006), which is prevalent north of 15°S in the atmospheric component of the CFSR reanalysis (not shown) and experiences a strong seasonality with a maximum in DJF. The uncoupled TR04 shows a wider and stronger Angola current in DJF, but its southern limits are close to CFSR.
A different situation is found in the coupled runs: both setups show a clear overshooting of the Angola current, reaching as far south as 26°S in AFR in DJF, while in JJA it is NAT that reaches the southernmost position (Fig. 11). This causes a southward shift of the Angola Benguela Front and larger warm SST bias in that area. This strengthening of the southward flow in AFR simulations can be related to reduced upwelling-favoring winds associated with a weaker SAA, and a stronger negative wind stress curl north of the Angola Benguela Front (Xu et al. 2014). Richter (2015) hypothesizes that westerly wind biases on the equator may also be partly responsible for the overshooting of the Angola current in GCMs because they lead to higher than observed sea level at the eastern boundary contributing to the strength of the current. In turn, warmer water is transported to the interior of the southern tropical Atlantic by advection and eddy transport related to the SEC, spreading the warm bias to this region (see Fig. 9).
In austral winter (JJA), when the SAA is of similar magnitude in both domains, the impact of the SAA on the ocean circulation is harder to discern. The spatial pattern of the Sea Surface Height and the associated geostrophic currents in both sets of simulations for JJA are similar to CFSR (Fig. 11). Therefore, the contribution of the transport out to the central basin to the biases in JJA is apparently not as important for JJA as for the austral summer.
Another proposed mechanism explaining the warm SST biases in the eastern subtropical oceans is related to the poor skill of the GCMs to reproduce the observed stratocumulus cloud fraction, leading to a modification of the surface radiative balance (e.g., Ma et al. 1996). In the AFR runs, the bias in downwelling surface radiation is negative over almost the entire Atlantic Ocean basin except off the southwest African coast in JJA (Fig. 12). The negative bias is due to an overestimation of the cloud cover. Additional experiments with the coupled model suite were carried out with increased thresholds for sub-grid cloud formation and reduced altitudes for cloud top before the cloud may rain. The large negative bias in radiation was removed due to a cloud cover decrease, but at the expense of a significant overestimation of rainfall in the ITCZ (not shown). In these experiments with increased and more realistic solar irradiance over the Atlantic Ocean basin, the patterns of the SST warm bias remained the same, thus corroborating the hypothesis of a dominant role of the response of the ocean dynamics to changes in surface wind stress. However, the increased local radiative forcing enhanced the positive SST anomalies in amplitude (not shown). Still, note that a positive radiation bias is observed west of the Angolan in JJA in the AFR simulations (Fig. 12), implying that an impact from anomalously low cloudiness in this marine stratocumulus zone to the warm bias is very likely. However, its exact contribution would require a thorough study of surface and ocean mixed layer energy budgets which is beyond the present study.

Relevance of the thermocline
The alongshore biases north of the Angola Benguela Front detected in the uncoupled simulations can be related to ocean model deficiencies in simulating the thermocline along the Angola coast, and could be remotely linked to biases in the equatorial thermocline (Lübbecke et al. 2010). They may also arise from inaccuracies in the data used to drive the ocean model. While the overshooting of the Angola current contributes to the biases south of the Angola Benguela Front, similar biases in the SETA SSTs are also found in experiments with the NCEP Coupled Forecast System (CFSR; Huang et al. 2007) and these have been attributed to a failure in simulating the sharp thermocline within the Benguela upwelling system (Xu et al. 2014). In order to gain more insight into the role of the upper ocean in the simulated SSTs, we compare the seasonal cycle of model SST, and the thermocline (represented by the depth of the 23 °C isotherm) with the corresponding data from Reynolds OISSTV2, PHC and the GODAS reanalysis for the two key regions (see Fig. 1c): the Angola Benguela Area (ABA: 20°S-10°S, 8°E-15°E) to represent the Angola Benguela Front (and the center of the Benguela Niños), and the Atlantic3 region (ATL3: 3°S-3°N, 20°W-0°E) to represent the equatorial Atlantic cold tongue. Additionally, we define the South Eastern Equatorial Atlantic (SEEA: 15°S-5°S, 10°W-8°E) to represent the region where the uncoupled TR04 bias off the African coasts expands in the coupled models.
The standalone TR04 provides a good simulation of the seasonal cycle of SST in the ATL3 region (Fig. 13). The model captures well the evolution of the seasonal cycle, reaching its maximum in March-April and featuring a minimum during July-August. Still, it is consistently colder than both GODAS and OISSTV2 through the year. This colder SST is associated with a shallower thermocline (Fig. 13). Both the AFR and NAT setups also simulate well the seasonal cycle of SST. In general, the coupling increases the SST in ATL3, especially in the second half of the year. This warming is similar for both setups, up to 2 °C in JJA. Therefore, the simulated ATL3 SST is closer to both observations and reanalysis in NAT and AFR runs than in the standalone TR04. This warmer SST in the coupled runs is associated with a deeper thermocline in these simulations compared to the standalone TR04 (Fig. 13) for the MAM and JJA seasons. A possible explanation for this in JJA is as follows: during these seasons, a coupled mode is excited, similar to the delayed oscillator operating interanually in the Pacific (Ding et al. 2009;Burls et al. 2011) and a deeper mean thermocline leads to a slower SST response. Moreover, the effect of the equatorial wind bias should be considered for the equatorial upwelling (which is reduced given a westerly wind bias). Additionally, the deeper thermocline in the coupled setups could weaken the Ekman feedback, turning the SST less sensible to the cooling by wind forcing. During austral winter and spring, with the shoaling of the thermocline due to the seasonal strengthening of the winds, local forcing by the equatorial zonal wind component and together with vertical advection of colder subsurface water across the thermocline from the ocean interior driven by Bjerknes-like ocean-atmosphere feedbacks play the main role in the middle and western Atlantic (Burls et al. 2011). For this season, two other mechanisms for easterly wind acceleration over the equatorial Atlantic are identified: (1) the pressure gradients induced by monsoon heating over western Africa and the northern tropical Atlantic is limited to the eastern Atlantic south of West Africa; and (2) the equatorward momentum advection is a significant contributor across the basin. The accelerated easterlies contribute to the equatorial cooling by inducing upwelling and by shallowing the thermocline in the east (Okumura and Xie 2006). Richter et al. (2012b) showed that JJA SST biases in the eastern Equatorial Atlantic in the CCSM global coupled model can be associated with a weakened east-west tilt of the thermocline due to a deeper thermocline in the eastern half of the equatorial Atlantic, which develops in response to a reduced strengthening of the MAM zonal wind stress.
In JJA, the SST is warmer in the east and colder in the west compared to the uncoupled run for both AFR and NAT domains. These improvements in the simulated surface temperature are associated with a shallower thermocline in the west and to a deeper thermocline in the east. This behavior is apparently associated with a weaker simulated easterly wind stress in MAM compared to ERA40 (Fig. 14). In both coupled simulations, the weaker easterlies in the equatorial band during MAM cause a shallower thermocline in the western part of this region, and a deeper thermocline in the center and the east. The ocean-atmosphere feedbacks active during the JJA season lead to a further weakening of the thermocline tilt and a warmer surface temperature in the eastern equator.
In the ABA region, AFR simulations show a much higher SST all across the year that is also associated with a deeper thermocline (Fig. 13). In the SEEA region, warm SST is associated with a deeper mean thermocline from July to December (causing a decrease of the interannual variability). As in the equatorial region, the remote forcing is also 1 3 less direct for the ABA region in JJA than in DJF. For JJA, the main source of warm biases is the too deep thermocline, that can be caused remotely as well as locally. This feature is also present in the uncoupled model. Because of the deeper thermocline, the mixing caused by the stronger alongshore winds brings warmer water to the surface layer. As the thermocline in the coupled setups is deeper than in TR04, the alongshore biases are more prominent for the coupled setups, this effect being stronger for AFR. The warmer subsurface temperature and deeper thermocline can be related to a warmer and stronger Angola Current observed in all the runs, but to a larger extent in AFR. According to Lass et al. (2000), when the AC encounters with the Benguela current at the Angola Benguela Front, it is forced to flow further poleward below the surface.
All simulations exhibit a deeper than observed thermocline in the region between 10S and 10N, this feature being stronger for AFR (Fig. 15a). The thermocline in the AFR coupled simulations also show a weaker tilt, with the western thermocline shallower in AFR compared to both TR04 and NAT simulations. The weaker SAA in AFR thus appears to have a strong effect in the vertical structure of the upper ocean layers across the subtropics. For example, the importance of a realistic simulation of the upwelling along the South African coast for a correct representation of the equatorial mixed layer temperature has been recently highlighted by Vizy and Cook (2013) using a coupled regional climate model/mixed layer ocean model. They found that the Atlantic cold tongue is much improved due to the horizontal advection of cold upwelled waters from the coastal South Atlantic with some contribution from vertical entrainment of cold water from below the mixed layer. Using the heat balance in the mixed layer, Hazeleger and Haarsma (2005) showed that a negative bias in the equatorial cold tongue can be attributed to entrainment flux at the base of the oceanic mixed layer. An enhanced entrainment deepens the mixed layer and causes strong reduction in the upper ocean divergence in the central equatorial Atlantic.

A coupling mechanism for the southern tropical Atlantic
The advection of the relatively cold water by the eastern Benguela current seems to play an important role in the thermal balance of the southern tropical Atlantic, and misrepresentations in the simulation of this transport could contribute to the warm bias there. Moreover, an inadequate representation of the transport can affect the representation of the whole SETA climate. Seager et al. (2003) suggested that a strong subtropical east-west SST gradient is important for the large-scale manifestation of the oceanatmosphere feedback responsible for the SAA strengthening. This east-west SST gradient is composed in the western part by a warm patch located between 20°S and 10°S due mainly to the transport of warm water by the Brazil current and a band of cold water in the east mainly due to the transport by the Benguela current and associated upwelling. The extension and position of this patch of cold SST is largely determined by the location of the Angola Benguela Front. In our simulations, a weaker DJF SAA is associated to a strong southward flowing Angola current with an insufficient equatorward advection of colder water by the Benguela current and as a consequence, a southward displacement of the Angola Benguela Front. This is documented in Fig. 11, which depicts the geostrophic velocity in the SETA for the CFSR reanalysis and our simulations. The resultant warming leads to an east-west gradient that is slightly weaker in NAT and much weaker in AFR than the gradient for standalone TR04 (Fig. 10).
In AFR this is true not only in DJF but in all seasons, and is related to much stronger warm anomalies in the eastern part of SETA that extend further west compared to the TR04 and the NAT simulations. Figure 9 illustrates the seasonality of these differences: for all seasons, the SST is higher for AFR than for TR04 and NAT all across the southern tropical Atlantic, except for the western subtropics. In the atmosphere, the weaker AFR SAA and the associated anomalous cyclonic circulation leads to weaker atmospheric cooling in the eastern ocean by subsidence over the warmer SST, and to slightly less warming in the southwestern Atlantic, in the region of the South Atlantic Convergence Zone by a weaker poleward advection of warm air from the equatorial region (Seager et al. 2003). The warmer SST acts in turn as a source of anomalous positive potential vorticity, further weakening the SAA (Miyasaka and Nakamura 2010).
As pointed out before, the austral summer SST for the AFR setup is higher in most of the south eastern tropical Atlantic and these biases may be associated with local (Benguela Upwelling System) and basin-wide feedbacks. In this season the AFR domains show a weaker land-sea thermal contrast that causes weaker alongshore winds in this region compared to ERA40. According to Miyasaka and Nakamura (2010), this could lead to a weaker South Atlantic high pressure system. The weaker local alongshore winds cause a strong southward flowing Angola current and a weaker northward flux of cool water along the South African coast and less upwelling and mixing in the ocean and by the South Atlantic current, warming further the water in the upwelling region. The remote forcing of the westerly bias in the equatorial winds may also contribute to modify the Angola current. The transport of this anomalous warm water off the Angola-Namibian coast by the South Equatorial current results in a net warming of the ocean interior and a warmer equatorial countercurrent. This anomalous warming is partially compensated by the higher heat loss from the ocean to the atmosphere over the cooler regions. The latent heat loss is particularly enhanced over the Gulf of Guinea due to the reduced subsidence. The warm bias is associated with enhanced atmospheric water content transport to this region in the AFR setup. These two factors lead to stronger precipitation over the Gulf of Guinea for the AFR setup. In turn, the atmosphere responds to warm SST bias by creating cyclonic potential vorticity. This potential vorticity contributes to the weakening of the alongshore winds, and as a consequence to the weakening of the SAA. Figure 15b illustrates the possible importance of this feedback for the simulation of the SAA: a weaker SAA is associated with a warmer SST.
The above results provide evidence that the seasonal biases in the southern tropical Atlantic SST are stronger when the SAA is weaker. In turn, the simulated SAA could be influenced by the SST biases, through a weaker eastwest zonal asymmetry (see Fig. 10). However, the experiments of  with uncoupled atmospheric models show that imposing an artificial east-west SST asymmetry in the subtropical regions does not change significantly the anticyclone in July. They argue that this fact could indicate that the response for standalone atmospheric models is largely due to local thermodynamic effects. We hypothesize that the coupled feedbacks enhance the atmospheric response to the SST east-west gradient through changes in the South Equatorial Current. Another important factor that can influence their results is that their experiments take place in the austral winter, while the importance of coupling seems to be manifested primarily in austral summer.
On the other hand, our NAT simulations display a region of stronger positive SST deviations from the zonal mean than the uncoupled simulations forced by ERA40 fluxes outside the coupled region. This could indicate that changes in the ocean circulation can be attributed to the coupling, and that the coupled feedbacks seem to make the atmospheric response to changes in the east-west SST gradient stronger and non-local through changes in the south Equatorial currents.
The importance of the described coupled feedbacks is clear in austral summer. The coupling also influences the JJA biases, albeit indirectly, through changes in the simulated thermocline depth. Figures 13 and 15a show that the AFR simulations have a deeper thermocline in the Angola Benguela Front Zone. In JJA, when the upwelling is stronger, the deeper thermocline leads to the upwelling of comparatively warmer water, contributing to the strong positive bias.

Conclusions
We assess the impact of a more realistic simulation of the SAA on the southern tropical Atlantic seasonal cycle with the regionally coupled model ROM. Our two coupled setups (AFR and NAT in Fig. 1b) include the southern tropical Atlantic region. The first domain covers the whole SAA, while the second domain only includes the northern part of the SAA but extends northeastward. With these two different coupled domains, we study the impact of the SAA on the tropical Atlantic seasonal cycle. Our main findings are as follows: • In DJF the AFR setup shows stronger SST biases off the African coast associated with a weaker SAA compared to both ERA40 and NAT, leading to a cyclonic surface wind bias. The atmospheric biases already present in standalone runs are modified in the coupled setup likely following a mechanism similar to the Wind-Evaporation-SST feedback (cf. Richter and Xie 2008). • The mechanisms responsible for the identified circulation biases in the tropical Atlantic are different for the DJF and JJA seasons. The differences in simulated MSLP and SST are smaller in JJA, when the SAA reaches its maximum strength. • In DJF, both local (along the Angola Benguela Front Zone) and basin wide ocean-atmosphere feedbacks (that could involve the ocean transport and the SAA) seem to influence the SST biases. A weaker SAA modifies the advective regime along the southern African coast, with a too strong Angola current bringing warm waters too far south, and a reduced northward flux of cold waters by the South Atlantic current. Together with a decrease in upwelling and mixing, and with possible remote influences from the equatorial westerly wind bias, this yields an anomalous warming within the upwelling region that spreads into the ocean interior by means of the South Equatorial current. These SST biases feedback onto the atmosphere through the generation of cyclonic potential vorticity that lessen the winds and, therefore, the SAA • In JJA, a too deep thermocline seems to be the most important factor for the SST biases. A deeper thermocline in the Angola Benguela Front Zone region in the AFR simulations leads to the upwelling of warmer water that contributes to the strong positive SST bias. These modeled ocean surface conditions impact on the West African boreal summer monsoon, leading to enhanced continental precipitation and reduced rainfall in the northern Gulf of Guinea. Concomitant subsidence may have enhanced modeled JJA rainfall over the Congo Basin in a Walker-cell type response. • While our analysis suggests that the model errors in downwelling radiation are more important in JJA than in DJF, they seem to be less important for the warm bias than ocean dynamics when seen over all seasons. However, a definite answer to this question would require a full consideration of atmospheric and ocean heat budgets. • The specification of the SAA south of the Angola Benguela Front Zone as external forcing improves the simulation of the seasonal cycle in the tropical Atlantic, revealing the fundamental role of this anticyclone in shaping the climatology of the region. This conclusion represents an instance of the utility of regional coupled models and the domain-shifting technique implemented here for the identification of the key contributors to the climate of a given area of interest.
We can conclude that coupled models aiming at reproducing the natural climate variability in this region of the globe should include an adequate representation of the SAA, which in the light of our results demands the nontrivial task of capturing the intricate loop of feedbacks involving the subtropical zonal SST gradient, land-sea heating contrasts, and ocean dynamics. The methodology developed in this paper would allow us to study the impact of the SAA on the interannual variability in the SETA region. In such a study, inter-basin influences, like the atmospheric bridge between the Pacific and Atlantic Oceans described in Alexander et al. (2002), would need to be taken into consideration. The atmospheric bridge paradigm relates SST warming in the tropical Atlantic Ocean to enhanced subsidence and reduced cloudiness as a consequence of warmer eastern Pacific waters. On the other hand, Atlantic Niños seem to enhance Pacific Niño events in recent decades (Rodríguez-Fonseca et al. 2009;Ding et al. 2012). Future work should therefore also analyze in detail the interaction of the northern and southern hemispheric subtropical anticyclones in the Atlantic Ocean with both the Atlantic and Pacific Hadley and Walker cell systems.