Diapycnal heat flux and mixed layer heat budget within the Atlantic Cold Tongue

Sea surface temperatures (SSTs) in the eastern tropical Atlantic are crucial for climate variability within the tropical belt. Despite this importance, state-of-the-art climate models show a large SST warm bias in this region. Knowledge about the seasonal mixed layer (ML) heat budget is a prerequisite for understanding SST mean state and its variability. Within this study all contributions to the seasonal ML heat budget are estimated at four locations within the Atlantic cold tongue (ACT) that are representative for the western (0°N, 23°W), central (0°N, 10°W) and eastern (0°N, 0°E) equatorial as well as the southern (10°S, 10°W) ACT. To estimate the contribution of the diapycnal heat flux due to turbulence an extensive data set of microstructure observations collected during ten research cruises between 2005 and 2012 is analyzed. The results for the equatorial ACT indicate that with the inclusion of the diapycnal heat flux the seasonal ML heat budget is balanced. Within the equatorial region, the diapycnal heat flux is essential for the development of the ACT. It dominates over all other cooling terms in the central and eastern equatorial ACT, while it is of similar size as the zonal advection in the western equatorial ACT. In contrast, the SST evolution in the southern ACT region can be explained entirely by air-sea heat fluxes.


Introduction
The tropical Atlantic ocean plays a key role for climate variability in sensitive regions of the surrounding continents: Previous studies suggest that inter-annual variability of sea surface temperature (SST) in the equatorial Atlantic impacts rainfall variability over northeast Brazil and the coastal regions surrounding the Gulf of Guinea Kushnir et al. 2006;Brandt et al. 2011b). In particular, inter-annual SST variability within the Atlantic cold tongue (ACT) is known to influence the onset and strength of the West African Monsoon (Marin et al. 2009;Brandt et al. 2011a;Caniaux et al. 2011). Hence, in order to improve predictions of rainfall variability, SST variability and its causes need to be understood.
Sea surface temperature variability in the tropical Atlantic is dominated by an annual signal formed by the seasonal appearance of the ACT (Fig. 1). During the development of the ACT in boreal spring/summer SSTs decrease by about 6 °C. Towards the end of the year SSTs steadily increase completing the annual cycle (Fig. 2a). On the other hand, the net surface heat fluxes within this region are dominated by a semi-annual signal due to the meridional migration of the Inter Tropical Convergence Zone (ITCZ). The discrepancy between the semi-annual atmospheric forcing and the annual cycle of SSTs suggests ocean dynamics to be of fundamental importance within this region. Furthermore, the Abstract Sea surface temperatures (SSTs) in the eastern tropical Atlantic are crucial for climate variability within the tropical belt. Despite this importance, state-of-the-art climate models show a large SST warm bias in this region. Knowledge about the seasonal mixed layer (ML) heat budget is a prerequisite for understanding SST mean state and its variability. Within this study all contributions to the seasonal ML heat budget are estimated at four locations within the Atlantic cold tongue (ACT) that are representative for the western (0°N, 23°W), central (0°N, 10°W) and eastern (0°N, 0°E) equatorial as well as the southern (10°S, 10°W) ACT. To estimate the contribution of the diapycnal heat flux due to turbulence an extensive data set of microstructure observations collected during ten research cruises between 2005 and 2012 is analyzed. The results for the equatorial ACT indicate that with the inclusion of the diapycnal heat flux the seasonal ML heat budget is balanced. Within the equatorial region, the diapycnal heat flux is essential for the development of the ACT. It dominates over all other cooling terms in the central and eastern equatorial ACT, while it is of similar size as the zonal advection in the western equatorial ACT. In contrast, the 1 3 ACT is the oceanic region with strongest net atmospheric heat gain in the Atlantic (Josey et al. 1999). Hence, in order to reduce SST during ACT development, oceanic processes need to redistribute large amounts of heat to explain the observed seasonal decrease of SST. To further understand the seasonal cycle of SSTs within this region, an analysis of all individual contributions to the mixed layer (ML) heat budget is inevitable.
Several investigations of the ML heat budget in the ACT region either rely on model simulations (Peter et al. 2006;Jouanno et al. 2011b) or observational data (Foltz et al. 2003;Wade et al. 2011;Hummels et al. 2013). Both approaches have their advantages and deficiencies. Model approaches can consistently evaluate the ML heat budget including all individual contributions. However, model results might be biased due to their sensitivity to the parameterization of unresolved physics, e.g. horizontal and vertical mixing. On the other hand, observational studies often lack information on individual terms of the budget, which can not be estimated from the available data sets. In addition, some of the estimated terms may lack reliability due to a limited data base in time and/or space.
Previous observational studies investigating the ML heat budget in the ACT region follow different approaches in order to exploit the existing data sets in the best possible way. Nevertheless, many of these studies lack information on the contribution of the diapycnal heat flux due to turbulent mixing at the base of the ML to the heat budget. Instead, this term is assessed as the residual of the budget, which additionally includes all accumulated uncertainties (Foltz et al. 2003(Foltz et al. , 2013Wade et al. 2011). Model studies of the ACT as well as the Pacific cold tongue (PCT) have suggested that diapycnal mixing is an important contributor to the ML heat budget (Chang 1993;Peter et al. 2006;Jouanno et al. 2011b). In an observational study, Gouriou and Reverdin (1992) suggested diapycnal mixing to vary seasonally within the ACT.
It was only recently that measurement programs could indeed resolve seasonal variability of upper ocean turbulence: From multi-year moored microstructure temperature measurements Moum et al. (2013) demonstrated that diapycnal mixing indeed controls the seasonal cooling within the PCT at 0°N, 140°W, while for the equatorial ACT (i.e. within the equatorial belt 2°S-1.5°N, 23°W-2°E), Hummels et al. (2013) highlighted seasonal and regional differences in diapycnal heat flux from the ML into the upper thermocline using a multi-cruise microstructure profiling data set. Hummels et al. (2013) found that the ML heat loss due to diapycnal mixing is of considerable magnitude and amounts to up to 90 W m −2 within the equatorial ACT region.
By solving the ML heat balance on the equator at 10°W, Hummels et al. (2013) showed that the diapycnal ML heat loss is the dominant cooling term for the ML heat budget during ACT development. In fact, they were able to balance the seasonal ML heat budget at this location for the period from May to November when incorporating estimates of the diapycnal heat flux. In the ACT as well as in the PCT region, turbulent mixing could be associated with shear instabilities (Gregg et al. 1985;Hummels et al. 2013). Within the entire equatorial ACT vertical shear of horizontal velocity is significantly elevated due to opposing currents, namely the westward-directed northern branch of the South Equatorial Current (nSEC) at the surface and the Equatorial Undercurrent (EUC) flowing eastward along the thermocline. Accordingly, the diapycnal heat flux estimated from individual cruises carried out at different locations within the equatorial ACT is elevated (Hummels et al. 2013). This suggests that diapycnal mixing could be Fig. 1 Satellite SSTs (TMI) in July 2001 showing the regional extent of the ACT. The four locations of the PIRATA buoys used for this study are marked by red rectangles. The black lines denote the area used for evaluating the average diapycnal heat flux included in the ML budget important for the ML heat budget in the entire equatorial ACT region. The importance of diapycnal mixing for the ML heat budget associated with shear instabilities forced by the nSEC/EUC system within the equatorial ACT was previously diagnosed in numerical simulations (Jouanno et al. 2011a, b).
However, elevated diapycnal heat flux appears to be confined to the equatorial band (2°S to 1°N) only. At the southern tip of the ACT at 10°S, 10°W, Hummels et al. (2013) find diapycnal ML heat loss to be considerably lower (their Fig. 15), which suggests other mechanisms are important to reduce SSTs in the off-equatorial ACT.
This study is motivated by the success of Hummels et al. (2013) in achieving a balance of the ML heat fluxes and tendency at the equator at 10°W between May and November during the times when estimates of the diapycnal heat flux were available and could be included. Our intention here is to clarify whether a closure of the ML heat budgets from different locations within the ACT can be achieved when explicitly incorporating estimates of the diapycnal heat flux from microstructure observations and to assess the relative contribution of the diapycnal ML heat loss to the budgets at the different locations. The diapycnal heat flux observations reported by Hummels et al. (2013) are complemented with some recent cruises and combined with long-term observations from the Prediction and Research moored Array in the Tropical Atlantic (PIRATA; Bourlès et al. 2008) and climatological products to estimate a seasonal climatology of the tendency and the flux terms contributing to the heat budget of the ML at different locations within the ACT. Furthermore, as microstructure data is only sparsely available, a parameterization for mixing based on stratification and shear is developed from the available data.
In Sect. 2 of this study the data set and method to calculate all individual contributions to the ML heat budget will be presented. Section 3 points out the different background setting at the four locations in the western (23°W), central (10°W) and eastern (0°E) equatorial as well as southern (10°S, 10°W) ACT region. It also includes a description of turbulent mixing within the ACT and the results on the ML heat budgets at the four locations. The results are summarized and further discussed in Sect. 4.

Data and methods
In order to accomplish a seasonal description of the ML heat budget within the ACT, several data sets are combined. A limiting factor is the availability of estimates of the diapycnal heat flux due to mixing processes across the base of the ML. This term is derived from estimates of the turbulent dissipation rates inferred from observations made by microstructure profilers during a multi-cruise program (Hummels et al. 2013). During this program several meridional transects were frequently repeated during different stages of ACT development covering the absence (January-April), the development (April-July) and the mature phase (August-December) of the ACT . The highlighted transects were conducted across the equator (2°S-1.5°N) along 23°W, 10°W and 2°E (see Fig. 1, Table 1). At the equatorial position of these transects (or close to the transect in the case of the 2°E transect), atmospheric and oceanic PIRATA observations are available at high temporal resolution (Fig. 1). Hence, all terms contributing to the ML heat budget will be estimated at the four PIRATA buoy sites at 0°N, 23°W; 0°N, 10°W; 0°N, 0°E and 10°S, 10°W. For the determination of the contribution of the heat flux due to advection additional information on surface velocities, SST and mixed layer depth (MLD) as well as their horizontal gradients are required from climatological products.

Shipboard observations
A unique data set of microstructure shear and temperature profiles along with conductivity-temperature-depth (CTD) profiles was collected on ten cruises to the ACT region during different stages of ACT development (Table 1). The multicruise program includes cruises of the French EGEE (Etude de la circulation océanique et des échanges océan-atmosphère dans le Golfe de Guinée) project (EGEE2-6) and the German projects Nordatlantik (M68/2, MSM18/2) and SFB754 (M80/1, GEOMAR4, MSM22). Four cruises (EGEE3, EGEE5, M68/2 and MSM18/2) were scheduled to coincide with the development of the ACT and onset of the West African Monsoon in boreal spring and early summer (end of May to July), while another two cruises (EGEE2 and EGEE6) were carried out during the mature phase of the monsoon (September-October), when the ACT is still well developed. Finally, three cruises (EGEE4, M80/1 and MSM22) were conducted in November while the ACT was warming.
Another cruise was undertaken during the absence of the ACT in March 2008 (GEOMAR4). Due to the latitudinal extent of the EUC, meridional transect data collected in the equatorial belt between 2°S and 1.5°N are representative of the mixing activity at the equator (Hummels et al. 2013). Hence, to improve the statistics of the equatorial estimate of the diapycnal ML heat loss due to turbulence, all available profiles in this latitudinal range are used. The only stations evaluated outside the equatorial belt were collected near 10°S, 10°W during EGEE3. For more details on the cruises the reader is referred to Hummels et al. (2013).
The microstructure shear and temperature data were collected using different microstructure profilers (MSS90L and MSS90D) manufactured by Sea&Sun Technology in cooperation with ISW-Wassermesstechnik. All profilers were equipped with two shear sensors (airfoil), a fast temperature sensor (FP07), an acceleration sensor, tilt sensors and standard CTD sensors. For a detailed description of the instruments the reader is referred to Prandke and Stips (1998). The profilers were adjusted to descend at 0.5-0.6 m s −1 . Noise levels of turbulent dissipation rates from the microstructure profilers are better than 1 × 10 −9 W kg −1 for the MSS90L and better than 5 × 10 −10 W kg −1 for MSS90D (Prandke and Stips 1998;Schafstall et al. 2010).
The sampling strategy pursued for the different cruises differs in comparison to previous microstructure studies conducted in the equatorial Pacific. Instead of sampling at a single location for a period of several days to several weeks (Peters et al. 1988;Moum et al. 1989Moum et al. , 2009Lien et al. 1995;Inoue et al. 2012), profiling was done at several locations during a single cruise usually separated by 0.5° latitude on the meridional transects. This was accomplished by integrating microstructure profiling into the CTD station program during all cruises. A minimum of 3 and up to 20 microstructure profiles were collected at a single station extending from the surface to between 150 and 300 m depth. For this study a total of 615 profiles collected on 178 stations are used.
A few microstructure profiles from a more recent cruise conducted on R/V Maria S. Merian in summer 2011, MSM18/2, were not considered in this study. R/V Maria S. Merian was put into service in 2006 and ship-board microstructure observations had not been performed in the equatorial region prior to this cruise. Although most of the profiles collected during MSM18/2 were in the range of previously collected profiles from the same location and season, a few extreme events in water depths between 15 and 50 m were observed in single profiles. Here, turbulent dissipation rates reach up to 5 × 10 −5 W kg −1 . These intense mixing events that were associated with very short time scales of <30 min were not evident in the other profiles from similar locations during the previous 9 cruises. Similarly, microstructure data collected from a Slocum glider equipped with a microstructure probe (MicroRider, manufactured by Rockland Scientific) during the period of the cruise did not show mixing events of similar magnitude during a 26 day record from the equator at 10°W. We thus assume that the presence of the R/V Maria S. Merian triggered Kelvin-Helmholtz instabilities in the water column locally, but the exact process could not be assessed. Although these extreme events pass through all quality control procedures in post-processing, it was decided to remove the three affected profiles from the data set of this cruise, as the source of these extremely short events remains ambiguous. The remaining MSM18/2 cruise data used here thus consists of 50 profiles from 17 stations.

Microstructure data processing
Processing of microstructure data and further inferring eddy diffusivities (K ρ ) and diapycnal heat fluxes (J h ) out of the ML is explained in detail in Hummels (2012) and Hummels et al. (2013). A summary is given below.
In the Reynolds-averaged equation for turbulent kinetic energy (TKE) the TKE dissipation rate (ε) is given by where u refers to the turbulent velocity fluctuations, x refers to the spatial coordinate system and the standard tensor notation of summing over three components is used (e.g. in Moum et al. 1995); ν is the kinematic viscosity of seawater. Under the assumption of local isotropy, these 12 terms are simply related leaving only one independent component to be estimated: ε = 7.5ν ∂u ∂z 2 (Hinze 1975;Pope 2000). For this study, dissipation rates of turbulent kinetic energy were determined from the airfoil shear data recorded at the sampling rate of 1024 Hz via the variance method. In fact, the shear wave number spectra, E du / dz , are related to the turbulent dissipation rate through integration, Gregg 1998). The shear wave number spectra were calculated from 1-s ensembles (1,024 values) corresponding to a vertical interval of 0.5-0.6 m. After the integration is performed, a variance loss correction due to the limited resolved wavenumber band between k min and k max , which were iteratively adjusted as explained in Schafstall et al. (2010), is applied by fitting to the universal Nasmyth spectrum (Oakey 1982). Turbulent eddy diffusivities for mass were estimated from dissipation rates of turbulent kinetic energy (ε) using the dissipation method of Osborn (1980): Here, Γ is the mixing efficiency and set to the constant value of 0.2 (for details see Hummels et al. 2013) and N −2 is the inverse of the buoyancy frequency. Buoyancy frequency was calculated from least squared fits to salinity and temperature using a depth interval of up to 30 m or less (3 m directly below the ML) to avoid biases in stratification. Furthermore, the diapycnal heat flux is calculated via J h = −ρc p K ρ ∂T ∂z with ρ and c p being the density of seawater and specific heat capacity, respectively, and ∂T ∂z the vertical temperature gradient. Eddy diffusivity and diapycnal heat flux profiles are calculated separately for every station with average station profiles of turbulent dissipation rates obtained from the individual profiles.

PIRATA data
For evaluating the ML heat budgets, incoming solar radiation, subsurface temperature time series, air temperature, relative humidity, wind speed and rainfall data from the PIRATA buoy sites at 23°W, 10°W and 0°E on the equator as well as at 10°S, 10°W were used. To determine the mean seasonal cycle, we used all available daily averages from 1997 to 2012.

Surface velocities
Near-surface velocities are constructed from a combination of the YoMaHa'07 (Lebedev et al. 2007) data set providing surface velocities from ARGO float surface drifts and surface velocities from drifters (www.aoml.noaa. gov/phod/dac/). The combination of data from different platforms enhances the available data base, which is in general limited as the equatorial region is characterized by diverging poleward flow and drifters (and to some extent also ARGO floats) tend to leave the equatorial region shortly after their deployment. Drifters are designed to follow the mean currents of the upper 15 m due to the attachment of a drogue, while ARGO floats drift at the surface. Some of the drifters lost their drogue after some time making their drift behavior similar to that of ARGO floats. The used drifter data set was updated for the new meta data information, where the drogue loss dates were corrected ). The different platforms have different sensitivity to the wind slip that is accounted for as explained below. The ARGO floats are corrected for the wind slip at the surface similar to the corrections applied to the undrogued drifters (Pazan and Niiler 2001) with differing coefficients A: where W x and W y are the zonal and meridional component of wind velocity (in units of m s −1 ) from NCEP/NCAR Reanalysis 6-hourly winds. The coefficient A was set to A1 = 1.14 × 10 −2 for the ARGO floats, A2 = 1.64 × 10 −2 for undrogued drifters and A3 = 7 × 10 −4 for drogued drifters as in Perez et al. (2013) to account for the different response of the observational platforms to the wind induced slip. The derived velocities were then averaged on a 1° × 1° grid of latitude and longitude with an overlap of 2° in both directions to enlarge the amount of observations within the individual grid boxes.
For this study, investigating the contribution of the individual terms to the ML heat budget within the seasonal cycle, average monthly estimates of surface velocities are required. As drifters as well as ARGO floats can be captured within the vortices of tropical instability waves (TIWs), which act on intraseasonal timescales, the possible bias due to these intraseasonal phenomena needs to be eliminated. In order to derive unbiased monthly estimates required for this study the maximum amount of independent data points combining the drifter and ARGO data base is used: data of the individual drifters is available on a 6 hourly grid, which does not represent independent data. Moored observations from Perez et al. (2013) reveal a de-correlation time scale between 7 and 10 days for meridional velocity, which also agrees with the temporal resolution of the ARGO float data. Accordingly, the data of individual drifters were subsampled by averaging over 10-day periods and then combined with the ARGO data base. The use of shorter de-correlation time scales that is associated with a larger weight on drifter data relative to ARGO float data results only in minor changes in the resulting advective heat flux as long as the de-correlation time scale is larger than 1 day.

SST
To accurately determine local horizontal SST gradients required for evaluating heat advection terms (see below) a SST climatology of higher resolution compared to the one available from Reynolds and Smith (1994) (1° × 1°) was needed: Moum et al. (2013) described an underestimation of the meridional heat advection when using a SST climatology of 1° horizontal resolution. Hence here, monthly averages of satellite SSTs from the Tropical Microwave Imager onboard the Tropical Rainfall Measuring Mission (TRMM) satellite (www.ssmi.com/tmi/) were used, which are available at a 0.25° × 0.25° horizontal resolution. A monthly-mean SST climatology was produced by averaging the monthly SST observations between 1998 and 2012 (TMI SSTs are available since December 1997).

Surface heat fluxes
The individual contributions to the net surface heat flux at the ocean's surface from the recently developed Trop-Flux (Praveen Kumar et al. 2012) product are considered. This recent product was especially developed for the tropical oceans and the individual fluxes are available at daily or monthly resolution from 1979 to present on a 1° × 1° grid. Monthly values are used to build an average longterm monthly climatology of the individual surface fluxes using data between 1997 and 2012 coinciding with the period of available PIRATA observations. The resulting climatological surface fluxes of the TropFlux product are then compared to the PIRATA observations after applying the COARE algorithm (Fairall et al. 2003). Additionally, climatological long wave radiation from the daSilva et al. (1994) surface marine atlas available on a 1° × 1° grid was used for comparison with the TropFlux product. The comparison between TropFlux and PIRATA/daSilva derived surface heat fluxes is presented in Appendix 1.

Mixed layer depth
Finally, to determine horizontal gradients of MLD (needed for the calculation of the entrainment term, see below), the ML climatology of de Boyer Montégut et al. (2004) was used for estimating the entrainment term. The long-term monthly means are only available at 2° resolution in both latitude and longitude, but are interpolated on a 1° × 1° grid to match the resolution of the other variables.

Methodology
To assess the relative importance of the diapycnal heat fluxes on the ML heat budget, the individual terms of the budget are estimated. The heat balance equation for the ML in the following form was first introduced by Stevenson and Niiler (1983). Since then, it has been frequently used in observational studies evaluating the individual contributions to the ML heat budget (Moisan and Niiler 1998;Wang and McPhaden 1999;Foltz et al. 2003) and reads: where h represents the thickness of the ML, and v and T are the vertically averaged velocity and temperature in the depth range between the surface and −h. v ′ and T′ are the deviations from the temporal means and v and T are the deviations from the vertical average. The overbar of the third term indicates temporal averaging. w entrain is the entrainment velocity and T −h is the temperature at the base of the ML. q net denotes the net surface heat flux, and q −h the net heat loss through the base of the ML. Here, q −h is determined from a combination of the penetrative shortwave radiation and the diapycnal heat flux at the base of the ML. From left to right the terms represent local heat storage, horizontal advection (divided into a mean and eddy term), entrainment, the vertical temperature/ velocity covariance (e.g. due to baroclinic flow within the ML) and the combination of net atmospheric heating and vertical turbulent diffusion at the base of the ML. The vertical temperature/velocity covariance involves the correlation of the possible departures of horizontal velocity and temperature from their vertically and temporally averaged values in the ML. Swenson and Hansen (1999) estimated the heat flux due to this term considerably smaller than the other contributions and hence this term will be neglected in the following as also described in Foltz et al. (2003). The evaluation of the individual terms closely follows the procedures described by Foltz et al. (2003). Seasonal cycles of different variables are obtained from daily PIRATA data, averaged on every given day of the year and subsequently averaged for each month of the year. The penetrative fraction of the shortwave radiation depends on the MLD and the optical transparency of seawater. Here, it is calculated via an e-folding decay following Wang and McPhaden (1999), who parameterized shortwave radiation absorbed in the ML as Q abs = Q short (1 − 0.45e −h/25m ). Latent and sensible heat fluxes were calculated from the PIRATA buoy data using the COARE algorithm (Fairall et al. 2003) for comparison with the TropFlux product, which uses the same algorithm (see Appendix 1). Mean horizontal advection is determined by multiplying the monthly mean velocities (calculated from the combined ARGO and drifter velocities) with climatological longterm SST gradients from the TMI satellite observations. The eddy term of horizontal advection is estimated indirectly from the residual between mean horizontal advection described above and total horizontal advection estimated as h(v · ∇T ) = h dT dt − ∂T ∂t (Swenson and Hansen 1999). The total time derivative is obtained from drifter SSTs, whereas the local derivative is estimated from the average monthly TMI SSTs. The entrainment velocity can be expressed as w entrain = ∂h ∂t + ∇ · hv (Stevenson and Niiler 1983) representing the local time derivative of the depth of the base of the ML from PIRATA temperature time series and a divergence term of the product of the MLD and surface velocity climatology.
MLDs were calculated from PIRATA subsurface temperature time series and taken from the MLD climatology of de Boyer Montégut et al. (2004). The climatology is built upon the temperature threshold criterion of 0.2 °C. In the study of Foltz et al. (2003) MLDs are calculated with the temperature threshold criterion of 0.5 °C, but using only the profiles between 5.00 and 7.00 local time as a reference for the threshold to avoid the influence of shallow diurnal mixed layers. Different temperature thresholds (0.2 and 0.5 °C) referenced against SSTs between 5.00 and 7.00 local time as well as the temperature gradient criterion described in Lorbacher et al. (2006) for PIRATA subsurface temperatures have been evaluated and compared to the MLDs of the climatology of de Boyer Montégut et al. (2004) at the PIRATA locations. Best agreement was found for the temperature threshold of 0.5 °C, which is therefore further used here to estimate MLDs from PIRATA subsurface temperatures.

Seasonal variability at the PIRATA sites
Among the four considered PIRATA buoy locations, which are representative for different regions within the ACT, large seasonal variability in the background setting is observed. Common to all regions within the ACT is the strong cooling of SSTs starting around April/May (Fig. 2a). The strongest cooling is observed on the equator at 10°W, which is sometimes referred to as the center of the ACT (Jouanno et al. 2011b). The structure of the seasonal variability of SSTs at the western (23°W) and eastern (0°E) edge of the ACT on the equator is very similar to 10°W. However, the southern ACT region exhibits a slower cooling phase followed by a delayed and more rapid warming phase. Towards the end of the year SSTs recover to warm levels at all locations reaching their maximum levels in March/April again.
The seasonal cycle of net surface heat fluxes from the TropFlux product (see Appendix 1 for a comparison between TropFlux and PIRATA) varies considerably within the ACT (Fig. 2b). The seasonal variation at all locations is mainly caused by variations in the incoming solar radiation and the latent heat flux, whereas the sensible heat flux and the outgoing long-wave radiation remain rather constant throughout the year (Foltz et al. 2003;Hummels et al. 2013). However, within the equatorial region, net surface heat flux is typically positive. In the western and central part of the equatorial ACT region (23°W and 10°W) the ML is warmed throughout the year by 50-100 W m −2 .
Note that here the net surface heat flux is calculated with the absorbed shortwave radiation, which is corrected for the amount of heat penetrating through the ML estimated from PIRATA sub-surface temperatures. Within the eastern equatorial ACT (0°E) the superposition of annual and semiannual cycles of the net surface heat fluxes leads to a net heat flux minimum from May to July when SSTs decrease. In the southern ACT region (10°S, 10°W) a strong annual cycle of net surface heat fluxes is observed, which includes a ML warming as well as a cooling phase. At this location the cooling phase of the ML due to a change in atmospheric forcing of nearly 200 W m −2 coincides with a decline of SSTs.
MLDs are in general shallow in the eastern equatorial Atlantic due to the shoaling thermocline. Accordingly, within the equatorial ACT region MLDs decrease from west (23°W) to east (0°E) (Fig. 2c). At all locations within the equatorial belt, MLDs exhibit a seasonal cycle with maximum MLDs of 40-50 m in boreal autumn and 15-20 m during boreal spring. Within the southern ACT region MLDs are generally larger than in the equatorial belt. In the south, the seasonal variation ranges from around 40 m in boreal winter/spring to maximum 80 m during boreal summer/autumn.
Winds in the ACT region are dominated by the westward trades (easterlies) surrounding the ITCZ, which migrates meridionally during the year. For the equatorial ACT region this migration imprints an annual signal on the wind stress magnitude with strongest winds in August, when the ITCZ is at its northernmost position (Fig. 2d). A weaker semiannual signal, which peaks in April/May and September/ October is superimposed on the dominant annual pattern yielding a double peaked structure. The southern ACT region exhibits only an annual variation in wind stress, which peaks in August, increasing the latent heat flux during this period. This increased latent heat flux contributes to the strong cooling via net surface heat fluxes in the southern ACT (Fig. 2b).
Zonal surface velocities in the equatorial ACT region, determined from a combination of ARGO float and surface drifter data, are subject to a relatively strong semi-annual cycle. This appears odd at first glance, as the wind forcing is dominated by the annual cycle with only a weak semiannual imprint. However, the semi-annual cycle in zonal velocity was explained by the presence of a resonantly forced basin mode (Cane and Moore 1981;Thierry et al. 2004;Ding et al. 2009). The current dominating the equatorial region is the northern branch of the South Equatorial Current (nSEC) (Lumpkin and Garraffo 2005;Hummels et al. 2013), (Fig. 2e). During boreal summer maximum westward velocities are observed of about 0.25-0.5 m s −1 depending on the exact location within the equatorial belt. In the central and eastern equatorial Atlantic, zonal velocities even reverse in sign during the seasonal cycle. Due to the vanishing Coriolis parameter at the equator, meridional velocities might be directly forced by the meridional wind component. As winds at the equator have a southerly (from south to north) component and are largest in boreal summer and fall, meridional velocities are directed in the same direction (Rhein et al. 2010). This is also consistent with the results from Perez et al. (2013), where positive (northward) meridional velocities are found at the equator at 23°W and 10°W. However, their magnitude is significantly reduced compared to zonal velocities (Fig. 2f). In the southern ACT region zonal as well as meridional surface velocities are of reduced magnitude compared to the equatorial region and do not show a distinct seasonal variation. Surface velocities within this region are dominated by the Ekman flow. According to the steady trade winds (south-easterlies), the Ekman flow is directed towards the southwest throughout the year (Fig. 2e, f).

Turbulent mixing within the ACT
In this section, the new and unique data set of microstructure observations acquired in the central and eastern equatorial Atlantic (Hummels et al. 2013) is used to estimate the diapycnal ML heat loss directly, rather than relying on residual estimates of this quantity as has been done previously (Wang and McPhaden 1999;Foltz et al. 2003Foltz et al. , 2013Wade et al. 2011). The main findings relevant for this study from Hummels et al. (2013) concerning the regional and seasonal variability of turbulent mixing are summarized in the following: Turbulent dissipation rates (ε) at the equator are significantly increased in the upper thermocline compared to off-equatorial locations (cf. Fig. 3a). Within the equatorial region turbulent dissipation rates in the upper thermocline are: 1. Elevated in the western equatorial ACT region in comparison to the eastern equatorial ACT region (cf. Figure 3a); 2. Most intense mixing is observed in boreal summer in the whole equatorial ACT region; Hummels et al. (2013) described a close correspondence between the seasonal and regional variability of background shear and stratification levels and ensemble mixing intensities, turbulent eddy diffusivities (K ρ ) and diapycnal heat fluxes (J h ): e.g. vertical shear squared of horizontal velocities (S 2 = (∂U/∂z) 2 + (∂V /∂z) 2 ) was found to reduce from 10°W towards the eastern equatorial ACT, while stratification (N 2 ) increased. This reduces the likelihood of shear instabilities to occur in the Gulf of Guinea and indeed turbulent parameters such as dissipation rates as well as diapycnal ML heat losses were observed to decrease in magnitude from the western equatorial region towards the east. Despite stronger zonal subsurface velocities at 23°W compared to 10°W (Brandt et al. 2011a), shear levels are reduced at 23°W compared to 10°W (Jouanno et al. 2011b). This can be explained by the shoaling of the EUC towards the east, which limits the depth range of opposite flowing currents, westwards at the surface and eastwards in the subsurface. The highest shear levels in the central equatorial ACT region correspond to highest mixing activity. Within the southern ACT, shear levels were significantly reduced due to the lack of strong current features and turbulent mixing was found to be low. The relation between background shear and stratification conditions and turbulent mixing activity will be further analyzed later in this study.
According to the described variability in shear levels, average summer dissipation rates below the ML range from 3.5 × 10 −7 W kg −1 in the eastern to 7 × 10 −7 W kg −1 in the western equatorial ACT (Fig. 3a). In the southern ACT  panels (a, b). Note that the vertical temperature gradient is shown in (c) as it is an important constituent in the formulation of the diapycnal heat flux dissipation rates below the ML do not exceed 5 × 10 −8 W kg −1 . At a depth of 50 m below the ML, equatorial dissipation rates have dropped by about one order of magnitude (Fig. 3a). Inferred eddy diffusivities (Sect. 2.1.1) just below the ML range from 1.2 × 10 −4 m 2 s −1 in the eastern to 7 × 10 −4 m 2 s −1 in the western equatorial ACT, while values in the southern ACT only reach 1.6 × 10 −5 m 2 s −1 at maximum. The magnitude and vertical structure of the observed turbulent parameters agrees rather well with those inferred from microstructure measurement programs carried out in the central equatorial Pacific (Gregg et al. 1985;Peters et al. 1988;Moum et al. 1989;Lien et al. 1995). However, in comparison to the central equatorial Pacific, Hummels et al. (2013) report for the equatorial ACT region a reduction in the night time enhancement of turbulence, which is referred to as deep cycle turbulence (Moum and Caldwell 1985). Nevertheless, to avoid possible biases due to an unevenly distributed sampling time of profiles during the day, mean profiles of turbulent parameters are derived here by separately averaging measurements collected during the day (08:00-20:00) and night (20.00-08:00) before calculating mean profiles. The average of the mean day and mean night profiles are then further incorporated into the ML heat budgets.
In order to obtain the turbulent contribution of diapycnal heat fluxes to the ML heat budget, the transition zone between the base of the ML and the stratified region below has to be accurately resolved. In the equatorial Atlantic as well as in the Pacific, profiles of the diapycnal heat flux are highly divergent below the ML. Maximum values are found at the base of the ML that rapidly decrease in deeper layers ( Fig. 3d; Lien et al. 2008). In most of the mean profiles (Fig. 3d), diapycnal heat flux at 20 m below the ML is significantly reduced. Hence, the amount of heat being extracted from the ML into the interior ocean is characterized by the diapycnal heat flux in a rather narrow layer. Here, diapycnal ML heat loss is determined by averaging the diapycnal heat flux profiles in the interval MLD + 5 m to MLD + 15 m. The reason to use this averaging interval is to ensure values from within the ML to be excluded from the estimate as the method from Osborn (1980), which we use here, is only valid in stratified sheared flow. Due to strong variability in stratification within a few meters below the ML, and due to the fact that stratification for the Osborn parameterization needs to be calculated over a larger length scale than turbulent overturns (that can be as much as several meters due to the strong mixing there), we decided to use this depth interval. Due to the elevated vertical divergence of the heat flux profiles, this approach leads to estimates of diapycnal heat loss of the ML being biased low. The error of this approach can be estimated from the average heat flux profiles by extrapolating from MLD + 10 m to the MLD. The extrapolation needs to be performed as the diapycnal heat flux obtained directly below the ML with the Osborn method can not be considered reliable as explained above. Overall, the diapycnal heat flux reduces by about 30 % at MLD + 10 m compared to the value directly below the ML.
Note that the MLD from vertically high resolved CTD profiles was generally calculated using the temperature threshold criterion with ΔT = 0.2 °C compared to the SST. Heat flux profiles and inferred diapycnal ML heat losses were calculated separately for every station. Subsequently, these station averaged profiles of the different sections of an individual cruise between 2°S and 1.5°N were averaged in day and night ensemble (see above) to obtain a single estimate that was taken as representative for diapycnal heat loss of the ML for the month in which the measurements were collected. Uncertainties for each individual estimate were calculated using Gaussian error propagation as described in Hummels et al. (2013).

Mixed layer heat budget
In the following, the contributions to the ML heat budget derived from PIRATA observations, climatological data sets as well as microstructure observations are combined at the four different locations within the ACT. Note that the individual contributions to the ML heat budget at 0°N, 10°W were already analyzed in Hummels et al. (2013). However, to achieve a consistent comparison between heat budgets at the different PIRATA locations, the analysis of the heat budget at 0°N, 10°W is repeated here partly using different data products (concerning the surface velocities and the net surface heat fluxes) compared to Hummels et al. (2013).
Several modeling as well as observational studies have addressed the ML heat budget in the eastern equatorial Atlantic previously (Foltz et al. 2003;Peter et al. 2006;Jouanno et al. 2011b;Wade et al. 2011;Hummels et al. 2013). Among these studies, the definition of high and low frequency advection terms, referred to as mean and eddy advection here (Sect. 2.4), varied. Hence, in order to compare the results amongst the different studies it has to be clarified which oceanic processes are attributed to the different terms. As has been pointed out in the tropical Pacific analysis by Wang and McPhaden (1999) processes reflected in the eddy advection term close to the equator, particularly in the meridional component, significantly depend on the latitudinal averaging interval over which the heat budget analysis is performed. For local heat budgets or budgets performed for small regions around the equator (i.e. ±2° in latitude), the eddy advection warms the ML as the effect of TIWs laterally advecting warm waters are explicitly resolved (Foltz et al. 2003;Jochum and Murtugudde 2006;Peter et al. 2006). When budgets are performed over a larger meridional extent, eddy advection will cool the ML as TIW contributions are averaged out and the Ekman divergence dominates. This study focuses on local heat budgets at the locations of the PIRATA buoys and thus requires evaluating the eddy heat fluxes as local as possible. Opposing this minimum regional extent requirement is the accuracy of velocity and SST gradient data, for which the statistical reliability of the individual variables increases when they are averaged over larger meridional intervals. As a compromise between locality and statistical reliability, a 2° latitudinal and longitudinal interval is used to evaluate mean and eddy advection terms. In doing so, the warming effect of TIWs will dominate the eddy advection term.
When the net surface heat flux is compared to the observed heat storage, large negative residuals are evident in the western (23°W) and central (10°W) equatorial ACT and a reduced residual in the eastern (0°E) equatorial ACT (Figs. 4,5). In contrast, within the southern ACT region, SST variability during ACT development can virtually be explained by the variability in net atmospheric forcing. Hence, the residuals in the equatorial region need to be explained by different oceanic processes, which probably also vary in their relative contribution within the seasonal cycle.
In the following, the respective contributions of atmosphere and ocean processes to the warming and cooling of the ML during the absence, development and mature phase of the ACT will be discussed. As the focus of this study is on the seasonal variability of the individual contributions to the ML heat budget, annual and semi-annual harmonics were fitted to the individual terms (except for the diapycnal ML heat loss) before illustration (Fig. 5).

0°N, 23°W
At the western edge of the ACT (23°W) in the central equatorial Atlantic, the ML is warmed by net atmospheric forcing (Fig. 2b, 5a) and by eddy heat advection. In the central equatorial Atlantic eddy heat advection is predominately controlled by TIWs (Wang and McPhaden 1999;Jochum and Murtugudde 2006;Peter et al. 2006). TIW activity here was reported enhanced in the beginning of the year, in boreal summer and autumn (Bunge et al. 2007;von Schuckmann et al. 2008), which agrees with periods of elevated eddy advection in this analysis (Fig. 5a). Cooling of the ML is achieved by subsurface processes (diapycnal mixing and entrainment) as well as mean heat advection (Fig. 5a).
During the absence of the ACT (January to April), the ML balance is dominated by net atmospheric forcing and eddy heat advection, the latter contributing to a warming of up to 50 W m −2 in January. During March and April, when the tropical Atlantic is uniformly warm and the meridional gradient of SST is very weak, the eddy heat advection reduces to zero.
During the development phase of the ACT (May to August) the net surface heat flux increases, mainly due to increased incoming solar radiation, which counteracts the observed cooling of SSTs. Eddy advection dominated by the effect of TIWs as discussed above contributes to the warming of the ML. However, there is a strong increase in ML cooling from zonal heat advection (60 W m −2 ) and entrainment (25 W m −2 ) during this phase. The elevated zonal advection term is due to the persistent westward flow (Fig. 2e) advecting cooler surface waters from the central ACT towards 23°W. In addition, elevated diapycnal heat loss at and below the ML further contributes to cool the ML: the two independent estimates for June from 2006 and 2011, 58 W m −2 (M68/2) and 54 W m −2 (MSM18/2) respectively, agree well. The cooling dominated by the diapycnal heat flux and zonal heat advection is strong enough to reduce SSTs despite the warming due to net surface heat fluxes and eddy advection.
Within the mature phase of the ACT (August to the end of the year) the net surface heat flux significantly warms the ML. Eddy advection still contributes significantly (50 W m −2 ) to the warming. Cooling provided by zonal heat advection decreases at the beginning of the mature phase, due to the reduction in surface velocities associated Fig. 4 Comparison of the net surface heat flux (absorbed shortwave radiation (corrected for the amount of shortwave radiation penetrating through the ML) + latent heat flux + sensible heat flux + outgoing longwave radiation) in blue and the observed heat storage in black for the four different locations within the ACT with the nSEC (Fig. 2e). Towards the end of the year zonal heat advection re-intensifies in accordance with reintensified zonal velocities associated with the nSEC. The diapycnal heat flux during this period is still large. The two independent estimates for November from 2009 and 2012 yielded 45 W m −2 (MSM22) and 30 W m −2 (M80/1) respectively. The meridional heat advection increases within this phase due to increasing meridional velocities (Fig. 2f) and increasing MLD (Fig. 2c). Entrainment still contributes a cooling of about 15 W m −2 during this phase. However, the gradual overall reduction of the cooling terms lead to a gradual increase in SSTs towards the end of the year.
Comparison of the sum of the individual terms contributing to the ML heat budget and the observed heat storage reveals a large residual of 30-100 W m −2 throughout the year, if the contribution of the diapycnal heat flux is omitted (Fig. 5b). This was also reported in the recent study of Foltz et al. (2013). Similar residuals (around 80 W m −2 ) have been reported in previous observational studies from the central equatorial Pacific (Wang and McPhaden 1999) as well as from the western equatorial ACT , where the diapycnal contribution could not be estimated. Including the resolved seasonal variability of the diapycnal heat flux into the sum of terms reduces the residual in boreal summer and November by more than a factor of 2 and closes the heat budget at least within the uncertainties (Fig. 5b). The diapycnal ML heat loss together with mean zonal advection are identified as the dominant contribution to the cooling of SSTs during ACT development at 0°N, 23°W.

0°N, 10°W
In the center of the ACT (10°W) the ML is warmed by the atmosphere and eddy advection as was observed for 23°W. The ML is cooled by subsurface processes (diapycnal mixing and entrainment) as well as the meridional heat advection (Fig. 5c). Zonal heat advection at this site is significantly reduced in magnitude compared to the western edge of the ACT (23°W). Although westward surface velocities associated with the nSEC are also intensified at this location (Fig. 2e), the low zonal temperature gradient in the center of the ACT and shallow MLDs lead to reduced zonal heat advection compared to 23°W.
During the absence of the ACT (January to April) the ML heat budget at 10°W is dominated by net atmospheric forcing. The largest oceanic contribution is the meridional heat advection cooling the ML, which balances the warming via zonal heat advection and eddy advection. The meridional velocity is as explained previously a direct response to the meridional wind forcing. Although elevated ML cooling due to diapycnal heat fluxes is anticipated during this period, this study lacks observational support for this hypothesis.
Within the development phase of the ACT (May to August), the net heat flux from the atmosphere increases due to the reduction in latent heat flux due to reduced wind speed and an increase in the incoming solar radiation (Figs. 2b, 5c). Additional warming is provided by the eddy advection due to TIWs, similar as discussed for 23°W. However, strong subsurface cooling leads, despite these warming terms, to a cooling of the ML. The relative constitution of the subsurface cooling at this location differs from the observations at 23°W: Zonal heat advection is significantly reduced compared to 23°W, whereas entrainment is of similar magnitude cooling the ML at a rate of about 20 W m −2 , which agrees with the results obtained by Rhein et al. (2010). The striking difference is the clear dominance of the diapycnal ML heat loss of up to 90 W m −2 over the other oceanic cooling contributions. As detailed in Hummels et al. (2013) elevated vertical shear of horizontal velocities increases the occurrence of shear instabilities leading to the elevated diapycnal heat fluxes during this phase.
During the mature phase of the ACT (August towards the end of the year) atmospheric warming stays on a rather high level. The cooling due to the diapycnal heat flux is still of considerable magnitude, but decreasing. The dominance of the diapycnal heat flux within the subsurface cooling terms reduces in favor of the meridional heat advection, which reaches 60 W m −2 towards the end of the year. This increase in meridional heat advection is due to enhanced meridional velocities presumably caused by increased southerly winds, temperature gradients and MLDs (Fig. 2). The warming effect of the eddy advection reduces during the mature phase and is accompanied by the zonal heat advection actually warming the ML towards the end of the year. This is due to the changing sign of the zonal velocity (the reversal of the nSEC between August and October; Fig. 2e) and the subsequent sign change for the zonal temperature gradient occurring in November and December (not shown). The gradual reduction in the total subsurface cooling leads to a gradual increase in SSTs during this phase.

Fig. 5
Left panels: Individual contributions to the ML heat budget at the different locations within the ACT (color code explained in the legend). Vertical red lines denote 95 % confidence limits for the diapycnal ML heat loss. Right panels: Sum of the individual contributions without (black) and with (red) the diapycnal ML heat loss; observed heat storage is in green. Grey shading denotes 95 % confidence limits for the sum of terms excluding the diapycnal ML heat loss, red vertical lines denote 95 % confidence limits including the diapycnal ML heat loss; see Appendix 3 for details ◂ Comparison of the sum of the individual terms of the ML heat budget to the observed heat storage reveals a large residual of up to 110 W m −2 when omitting the contribution of the diapycnal heat flux during boreal summer and autumn (Fig. 5d). This was already reported by Foltz et al. (2003), who performed a similar study at this location. Consideration of the diapycnal heat flux as a contributing term yields in a closure of the budget within the uncertainties from June to November. However, as further described below, ML cooling due to diapycnal mixing is likely to contribute to the ML budget also during the absence of the ACT.

0°N, 0°E
At the eastern edge of the ACT (0°E) the ML is mainly warmed by the net atmospheric forcing and cooled by the diapycnal heat flux (Fig. 5e). The other oceanic contributions do not exceed 20 W m −2 throughout the year. Zonal heat advection at this location is negligible throughout the year due to reduced zonal temperature gradients and MLDs (Figs. 1, 2c), whereas entrainment acts to slightly cool the ML throughout the year, similar to what has been observed at 10°W. However, as pointed out above the net surface heat flux is significantly reduced in the eastern, equatorial ACT region compared to the more western and central locations (23°W, 10°W, Figs. 2b, 4). Accordingly, less subsurface cooling is required to decrease SSTs at this location.
During the absence of the ACT (January to April) in the beginning of the year the ML heat budget is determined by net atmospheric forcing and oceanic contributions do not exceed 10 W m −2 .
Incoming solar radiation reduces during ACT development (May to August) while the latent heat flux slightly increases, which leads to the reduction in net surface heat fluxes during this period. Concurrently, the diapycnal heat flux as inferred from microstructure observations during June 2006 and 2007 increases to its maximum value at this location of 21 and 29 W m −2 respectively and dominates the subsurface cooling (Fig. 5e). Note that at this location the seasonal variability of the diapycnal heat flux was composed from available estimates at 0°N, 0°E as well as around 0°N, 2°E. MLDs are observed extremely shallow at this longitude (Fig. 2c). Hence, the additional subsurface cooling by diapycnal mixing, which superimposes on the reduced net surface heat flux, is sufficient for the strong decrease in SSTs. Eddy advection is negligible during this phase. At this longitude within the Gulf of Guinea TIWs have not been detected. Instead, Athie and Marin (2008) as well as the numerical simulation of Jouanno et al. (2013) suggest intraseasonal variability in the Gulf of Guinea dominated by wind-forced Yanai waves having long zonal wavelengths and a period between 10 and 20 days. From our analysis it seems that the effect of these waves with long zonal wavelengths on the eddy advection is significantly reduced compared to the effect of TIWs in the western equatorial ACT. Nevertheless, they might contribute to enhance vertical shear of horizontal velocity and thereby favor diapycnal mixing at this location as suggested by Jouanno et al. (2013).
During the mature phase of the ACT (August to the end of the year), a similar evolution towards the end of the year is observed as at 10°W on the equator. Net surface heat fluxes increase due to the increase in the incoming solar radiation. Meridional and eddy heat advection both contribute about 20 W m −2 during this phase but with opposite sign. The cooling by the diapycnal heat flux reduces, which leads together with the increased warming by the atmosphere, to the retraction of the ACT towards the end of the year.
Comparison of the sum of individual terms to the ML heat budget and the observed heat storage reveals a residual of up to 30 W m −2 when omitting the contribution of the diapycnal heat flux (Fig. 5f), which is within the uncertainties at this location. However, if the diapycnal heat flux is included, the residual reduces minimum by a factor of 2 (Fig. 5f). Even if reduced in magnitude compared to the western and central equatorial region, the diapycnal heat flux provides the largest subsurface cooling term within the development phase of the ACT and hence seems to supply the essential contribution to cool SSTs.
No detailed observational study of the ML heat budget as far east as 0°E in the Gulf of Guinea has been published so far. Wade et al. (2011) determined the individual terms of the ML heat budget for considerable larger regions (about 5° latitude and 9° longitude). Their box 5, representative for a region including 0°N, 0°E, shows a considerable larger residual term of maximum 80 W m −2 during ACT development, which they associated with the diapycnal heat flux. Note, though that their box 5 extends until 6°W, where the diapycnal heat flux may still be elevated compared to 0°E. In addition, their estimate of the net surface heat flux within this box ranges from 50 to 120 W m −2 and is above our estimates at 0°N, 0°E especially during ACT development (Fig. 2b). Thus, for their study a larger cooling by oceanic processes is required to match the observed heat storage.

10°S, 10°W
Investigation of the background setting at the different PIRATA sites already revealed fundamental differences between the equatorial and the southern ACT region (Fig. 2). The net surface heat flux actually cools the ML during ACT development in contrast to the equatorial ACT region, where net atmospheric warming is observed throughout the year and at all locations. In general, surface velocities are weaker at this location compared to the equatorial region (Fig. 2e, f) and are not associated with distinct current features, but with the Ekman flow. Meridional heat advection tends to warm the ML at this location, while entrainment, zonal advection and eddy advection mostly cool the ML.
During the absence of the ACT (January to April) the increase in ML heat storage is balanced by net atmospheric fluxes, warming the ML (Fig. 5g).
During the development of the ACT (May and August) net surface heat fluxes significantly cool the ML by up to 90 W m −2 (Figs. 2b, 5g) due to the increase of the latent heat flux associated with increased winds as well as a reduction in the incoming solar radiation (Figs. 2d, 5g). Surface currents associated with the Ekman flow are directed towards the southwest. Meridional advection contributes to warm the ML during this period due to the southward advection of warmer waters from the north (see Fig. 1). The other advection terms as well as entrainment are rather small during the development phase of the ACT. Zonal heat advection is small due to the lack of a significant zonal temperature gradient (not shown). The diapycnal heat flux inferred from microstructure observations during June 2006 provides a negligible contribution to the cooling (Fig. 5g). Thus, the cooling of the ML and hence SSTs during ACT development at this location is dominated by the net atmospheric forcing. Hence, the deepening of the ML from April to August can not be associated with mechanical mixing, but must be due to an increase in wind stress curl or surface buoyancy flux.
During the mature phase of the ACT (August to the end of the year) the net surface heat flux increases and gradually warms the ML again. This is due to the increase in incoming solar radiation and a reduction of the latent heat flux due to decreasing winds. The diapycnal heat flux has not been estimated at this location during this phase in the seasonal cycle. However, highest mixing activity throughout the equatorial ACT was observed in boreal summer (Hummels et al. 2013), when vertical shear of horizontal velocities is strong. As the vertical shear at this location is in general rather low due to the lack of strong current features, it seems unlikely that this term gives an important contribution to the ML heat budget during this phase.
The sum of contributing terms closely follows the evolution of the observed heat storage at this location (Fig. 5h). The residual is negligible throughout the year in agreement with the previous estimates of Foltz et al. (2003). This agrees with the fact that the diapycnal heat flux at this location was estimated rather low during the main mixing season of the equatorial ACT (Figs. 3, 5g). The decrease in SSTs at this location during ACT development is governed by the net atmospheric forcing.

Summary and discussion
A unique multi-cruise microstructure data set, observations from the PIRATA buoy network as well as climatological data sets are used to investigate the individual contributions to the seasonal ML heat budget at four different locations within the ACT region with a special focus on the role of the diapycnal heat flux. The ML budgets within different phases of the seasonal cycle associated with the absence, development and mature phase of the ACT are described. Microstructure observations allowed estimating the contribution of the diapycnal heat flux within the entire ACT region directly and not as a residual as has been done previously (Wang and McPhaden 1999;Foltz et al. 2003Foltz et al. , 2013Wade et al. 2011). Here, the heat budgets are presented for four PIRATA buoy sites at 0°N, 23°W; 0°N, 10°W; 0°N, 0°E and 10°S, 10°W thus extending the work of Hummels et al. (2013) focusing on the heat budget at 0°N, 10°W. The major result is that the diapycnal heat flux is a dominant cooling term for the ML heat budget in the entire equatorial Atlantic during ACT development. Towards the eastern equatorial ACT region, the discrepancy between the net surface heat flux and the observed heat storage diminishes. Likewise the magnitude of diapycnal heat flux is reduced from the western, equatorial ACT towards the east in observations (Hummels et al. 2013) and models (Jouanno et al. 2011b). However, even in the eastern equatorial ACT it was found to be the largest contribution to ML cooling.
The dominance of the diapycnal ML heat loss within the entire equatorial region has been suggested from modeling studies (Peter et al. 2006;Jouanno et al. 2011b), but up to now was only assessed from direct observations for 0°N, 10°W (Hummels et al. 2013). Note that the results for 0°N, 10°W presented here are essentially the same as in Hummels et al. (2013) although different data sets are used for the surface velocities and the net surface heat fluxes. This points towards the fact that the obtained results concerning the individual contributions at this location are rather robust. Outside of the equatorial region, the variability of ML heat content within the ACT is set by net atmospheric forcing and horizontal advection.
The ML heat budget during the development of the ACT from May to August in the equatorial ACT region consists of the warming by net atmospheric fluxes and eddy advection and cooling by subsurface processes and horizontal advection. The cooling is dominated by the diapycnal heat flux at the base of the ML within the entire equatorial ACT. Eddy advection, which is associated with the lateral advection of heat by intraseasonal waves, moderates the cooling within the western and central equatorial ACT. Within the southern ACT region the cooling of SSTs during ACT development is driven by net atmospheric forcing.
During the mature phase of the ACT (August to the end of the year) net atmospheric forcing increases throughout the ACT region. Additional warming by eddy advection is still present. In the equatorial belt the cooling by the diapycnal heat flux decreases, while meridional heat advection becomes more important in cooling SST towards the end of the year.
During the absence of the ACT (January to April), the increase in ML heat content at all locations within the ACT is dominated by net atmospheric forcing. The role of diapycnal heat fluxes during this period is hardly resolved from the multi-cruise data set. This will be further discussed below.
In general, the observed heat storage can be explained by the sum of terms, if the diapycnal heat flux is included. However, residuals between the sum of contributing terms remain especially within the western equatorial ACT at 0°N, 23°W. In addition, the uncertainties attributed to the sum of contributing terms are rather large. It should be noted though that here the 95 % confidence limits were presented, while in other studies only the standard error is considered. In the study of Foltz et al. (2003), the residual excluding the diapycnal ML heat loss at 23°W on the equator is of maximum 50 W m −2 and only present during the first half of the year. In a more recent study Foltz et al. (2013) revised this estimate using another surface velocity product. The differences between Foltz et al. (2003) and the results presented here can be attributed to extremely large zonal advection on the order of about 120 W m −2 during June-August and extremely low eddy advection from October to December, which reduce the residual in Foltz et al. (2003) during the second half of the year. Hence, the use of other surface velocity climatologies can lead to different results for the heat advection terms. This is also illustrated in Appendix 2 and demonstrates the necessity for improved surface velocity fields in order to get reliable estimates for the heat advection terms and their effect on the ML heat budget.
Unfortunately, and due to the scheduling of the microstructure measurement program the variability of the diapycnal heat flux during the absence of the ACT is not adequately resolved. However, Hummels et al. (2013) described a close agreement of the variability in the magnitude of turbulent parameters and the variability of large scale shear and stratification levels. Due to this good correspondence of vertical shear of horizontal velocities (S 2 ) and stratification (N 2 ) and the observed mixing intensities below the ML in the ACT region, Hummels (2012) proposed a simple relation of the form ε = a · Ri b in order to estimate the seasonal variability of dissipation rates from observations of shear and stratification only. Here Ri = N 2 /S 2 is the gradient Richardson number. The coefficients a and b were best fit to all available cruise data for the equatorial ACT region directly below the ML (MLD + 5 m-MLD + 20 m), which resulted in a = 4 × 10 −8 and b = −1.2. To obtain this fit, shear was calculated from horizontal velocities observed with the vessel mounted ADCPs, which had a bin size of 8 m. This simple relation seemed to provide reasonable results in terms of magnitude and seasonal variation of the dissipation rates and further inferred diapycnal ML heat losses when applied to independent shear and stratification observations. Independent observations of stratification were estimated from the subsurface measurements of temperature and salinity at fixed depths at the PIRATA buoys at 10°W and 23°W on the equator. Vertical shear of horizontal velocity was obtained from ADCPs moored in the proximity of the PIRATA buoy locations. Using the above mixing parameterization, the diapycnal ML heat loss on the equator at 23°W ranges on average between 22 W m −2 at the beginning of the year and 45 W m −2 during early summer (Fig. 6a). Here only data of 2002 and 2009-2012 were used, as within these years shear and stratification directly below the ML were adequately resolved. Maximum summer heat losses due to diapycnal mixing within the individual years from this method range between 33 and 67 W m −2 , which covers the range estimated directly from shipboard microstructure observations in summer (58 and 54 W m −2 ). Similar good correspondence is found at 10°W on the equator, where on average the parameterized diapycnal heat flux ranges between 40 and 65 W m −2 , whereas individual summer maxima range between 70 and 80 W m −2 compared to the maximum 90 W m −2 inferred from direct microstructure observations. Hence, this additional analysis tends to support the seasonal estimates obtained with the direct microstructure observations and could provide an estimate of the entire seasonal cycle of the diapycnal ML heat loss. In addition, the parameterized diapycnal ML heat losses suggest that even during the absence of the ACT the diapycnal ML heat loss is among the largest cooling terms in the western equatorial ACT. Incorporating these parameterized seasonal cycles of the diapycnal ML heat loss into the equatorial ML balances at 23°W and 10°W would lead to the closure of the ML heat budgets for the entire seasonal cycle within the uncertainties.
Despite the consistent results mentioned above further assessment of the magnitude and variability of the diapycnal heat flux via parameterization schemes is required. To fulfill this task, long term observations of turbulence and background flow and stratification parameters are needed. This could be accomplished either with moored microstructure observations (e.g. Chi-Pods, Moum et al. 2013) or by gliders mounted with microstructure probes profiling near subsurface ADCP moorings.
This study contributes to improve the understanding of the variability of SSTs in the eastern equatorial Atlantic.
The dominant contribution of the diapycnal heat flux to the ML heat budget within the equatorial ACT points out the importance of its correct representation in numerical simulations in order to improve biases in SSTs, which are of crucial relevance to reliable climate predictions.
Acknowledgements The authors thank Bernard Bourlés for the cooperation during the EGEE field campaign. We also want to thank the crews of the EGEE campaigns, R/V Meteor and Maria S. Merian as well as our technical group for their help with the field work. Financial support for this study was provided by the Deutsche Forschungsgemeinschaft as part of the Emmy Noether Program (RH and MD, grant DE1369/1-1), the PREFACE Project (grant 603521), the European Union 7th Framework Programme (FP7 2007-2013) GROOM Project (MS, grant 284321), the Sonderforschungsbereich 754 "Climate-Biogeochemistry Interactions in the Tropical Ocean" and through several research cruises with R/V Meteor and R/V Maria S. Merian as well as by the Bundesministerium für Bildung und Forschung as part of the Verbundvorhaben NORDATLANTIK (03F0443B) and RACE (03F0651B). The authors also acknowledge the PIRATA program for free data access.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.

Appendix 1: Comparison of TropFlux and PIRATA derived surface heat fluxes
To evaluate the TropFlux net surface heat fluxes, the seasonal cycles of the latent and sensible heat fluxes were estimated from the PIRATA observations using the COARE algorithm (Fairall et al. 2003). In addition, the contribution of the long-wave radiation as given by daSilva et al. (1994) was compared to TropFlux. In general, the different estimates for the individual terms agree rather well with maximum discrepancies of 20 W m −2 (Table 2; Fig. 7) for the combined net surface heat flux (shortwave and monthly averages of the individual years (dashed grey lines). Black circles mark the diapycnal ML heat loss estimated from microstructure observations together with their confidence limits as presented in Fig. 5a, c respectively Table 2 Average differences between the monthly estimates of individual contributions to the surface heat flux from TropFlux and PIRATA As a measure of the spread of these average differences, the standard deviation of the difference is given. As additional information the amount of daily PIRATA observations is given as percentage of the gap-free daily TropFlux time-series between 1997 and 2012 (5844 days) 0°N, 23°W 0°N, 10°W 0°N, 0°E 10°S, 10°W Average difference in surface heat flux components between TropFlux and PIRATA/daSilva (for outgoing longwave radiation) Incoming shortwave radiation (W m −2 ) 5.8 ± 4.5 14.4 ± 11 9.4 ± 9.1 13.2 ± 8.9 Outgoing longwave radiation (W m −2 ) 3.3 ± 1.3 4.6 ± 2.7 6.5 ± 4.0 7.2 ± 3.1 Latent heat flux (W m −2 ) 15.0 ± 7.1 6.8 ± 4.0 9.3 ± 4.5 7.6 ± 4.1 Sensible heat flux (W m −2 ) 1.1 ± 1.2 1.0 ± 0.5 1.7 ± 1.1 4.8 ± 3.0 Net heat flux (W m −2 ) 10.7 ± 9.8 18. radiation + outgoing long-wave radiation + sensible heat flux + latent heat flux). This does not come as a surprise as the PIRATA data are included in the TropFlux product and latent and sensible heat fluxes are also evaluated with the COARE algorithm. However, some discrepancies exist, which can be probably explained by data gaps in the PIRATA records (Table 2) and by quality control procedures used to derive the TropFlux product. As some of the PIRATA estimates are based on a considerably smaller data base (e.g. only 38 % of daily observations at 0°N, 10°W for latent and sensible heat flux compared to the gap-free Trop-Flux time series, Table 2) and TropFlux seems to provide reasonable estimates for all contributions to the net surface heat flux at these locations, TropFlux will be further used for the analysis of the ML heat budgets.

Appendix 2: Uncertainties in heat advection terms due to choice of surface velocity product
The variability of the estimates of the advection terms caused by the use of different surface velocity products is considered in the following. Two additional surface velocity products, namely the drifter climatology of Lumpkin and Garraffo (2005) and the OSCAR product (Bonjean and Lagerloef 2002) were used to calculate the advection terms. These two products are available on a 1° × 1° grid in longitude and latitude and climatological monthly values are used. To illustrate the spread in the resulting advection terms an error bar is constructed as the standard deviation of the monthly advection term values, when the advection terms are estimated with the different surface velocity products (Fig. 8, top panels). In general, the uncertainty of the advection terms is larger in the second half of the year at all locations. Dependent on the surface velocity product, not only the magnitude, but also the sign of the resulting advection term can differ during different periods within the year. This is due to opposing velocities among the products, e.g. in spring and autumn at 23°W for the zonal as well as the meridional velocity (Fig. 8, lower panels). Particularly meridional velocities at 23°W, which are negative (against the wind) in the drifter climatology of Lumpkin and Garraffo (2005) and hardly existent in the OSCAR product within the entire equatorial Atlantic, lead to the use of the combined (ARGO and drifter) surface velocity product for this study. In summary, it appears that additional work has to be done to improve surface velocities within the equatorial Atlantic as they alter the effect of the heat advection terms on the ML heat budget.

Appendix 3: Error estimates
Total errors of each monthly estimate of the individual terms contributing to the heat budget of the ML are comprised as err total = err 2 data + err 2 seasonal . err data is considered as an error inherent to the instrumentation or the algorithms used to derive the estimate and is usually given by the data provider. For the TropFlux product Foltz et al. (2013) find the monthly err data for the different heat flux components to be 8 W m −2 for latent heat flux and shortwave radiation, 2 W m −2 for outgoing long-wave radiation and 1 W m −2 for sensible heat flux, accordingly for the wind stress 0.004 N m −2 is calculated. The accuracy of PIRATA temperatures is given as 0.03 °C (Freitag 1994). For TMI SST an error of 0.1 °C is considered (http://www.ssmi.com/tmi/tmi_validation.html). err seasonal is the statistical uncertainty of each monthly estimate resulting from a finite length of the data record Fig. 7 a Comparison of the seasonal cycle at 0°N, 10°W of the different components of the net surface heat flux as obtained from Trop-Flux and PIRATA (incoming shortwave, latent and sensible heat flux)/daSilva (outgoing longwave), b daily estimates of the latent flux at 0°N, 10°W from TropFlux (red) and PIRATA (blue); the correlation coefficient for TropFlux and PIRATA is 0.57 (significant at the 95 % confidence level) during the periods when they coexist and is estimated as the standard error (SE) of each calendar month.
Concerning the surface velocities, no statement about the individual data error of drifter velocities could be found on the website (www.aoml.noaa.gov/phod/dac/). However, the natural oceanic variability in monthly surface currents, which is dominated by averaging over TIWs, is much larger than the data error for these individual observations. Hence for the surface velocities err total = err seas and is estimated as the SE for all independent velocity estimates (ARGO and drifter) falling in a grid box within each month of the year.
As the uncertainty for terms associated with turbulent mixing is described in terms of 95 % confidence limits (Hummels et al. 2013), all uncertainties (err total ) within this study are converted to 95 % confidence limits (CL95) for consistency. The conversion from SE to CL95 is via CL95 = x ± 1.96 · SE, whereas x is the monthly average of term x.
CL95 for turbulent dissipation rates are derived via bootstrapping (Hummels et al. 2013), as dissipation rates are not distributed log-normally and therefore the standard error should not be used.
To derive error estimates for terms comprised of different variables, which are attributed with their individual errors estimated as described above, standard error propagation is used. Fig. 8 Top The advection terms at the four locations within the ACT: zonal heat advection in grey, meridional heat advection in orange and eddy heat advection in cyan calculated for the surface velocities used in this study (combined ARGO and drifter velocities). The error bars are constructed as the standard deviation of monthly values derived by using the different surface climatology products (see text for details). Bottom Seasonal cycle of equatorial zonal (u) and meridional (v) surface velocities for the different velocity products