Articles | Volume 18, issue 8
Clim. Past, 18, 1947–1962, 2022
Clim. Past, 18, 1947–1962, 2022
Research article
26 Aug 2022
Research article | 26 Aug 2022

A 15-million-year surface- and subsurface-integrated TEX86 temperature record from the eastern equatorial Atlantic

A 15-million-year surface- and subsurface-integrated TEX86 temperature record from the eastern equatorial Atlantic
Carolien M. H. van der Weijst1, Koen J. van der Laan1, Francien Peterse1, Gert-Jan Reichart1,2, Francesca Sangiorgi1, Stefan Schouten1,2, Tjerk J. T. Veenstra1, and Appy Sluijs1 Carolien M. H. van der Weijst et al.
  • 1Department of Earth Sciences, Utrecht University, 3584 CB Utrecht, the Netherlands
  • 2NIOZ Royal Netherlands Institute for Sea Research, 1797 SZ `t Horntje, the Netherlands

Correspondence: Carolien M. H. van der Weijst (


TEX86 is a paleothermometer based on Thaumarcheotal glycerol dialkyl glycerol tetraether (GDGT) lipids and is one of the most frequently used proxies for sea-surface temperature (SST) in warmer-than-present climates. However, GDGTs are not exclusively produced in and exported from the mixed layer, so sedimentary GDGTs may contain a depth-integrated signal that is also sensitive to local subsurface temperature variability. In addition, the correlation between TEX86 and SST is not significantly stronger than that to depth-integrated mixed-layer to subsurface temperatures. The calibration of TEX86 to SST is therefore controversial. Here we assess the influence of subsurface temperature variability on TEX86 using a downcore approach. We present a 15 Myr TEX86 record from Ocean Drilling Program Site 959 in the Gulf of Guinea and use additional proxies to elucidate the source of the recorded TEX86 variability. Relatively high GDGT[2/3] ratio values from 13.6 Ma indicate that sedimentary GDGTs were partly sourced from deeper (>200 m) waters. Moreover, late Pliocene TEX86 variability is highly sensitive to glacial–interglacial cyclicity, as is also recorded by benthic δ18O, while the variability within dinoflagellate assemblages and surface/thermocline temperature records (U37k and Mg/Ca) is not primarily explained by glacial–interglacial cyclicity. Combined, these observations are best explained by TEX86 sensitivity to sub-thermocline temperature variability. We conclude that TEX86 represents a depth-integrated signal that incorporates a SST and a deeper component, which is compatible with the present-day depth distribution of Thaumarchaeota and with the GDGT[2/3] distribution in core tops. The depth-integrated TEX86 record can potentially be used to infer SST variability, because subsurface temperature variability is generally tightly linked to SST variability. Using a subsurface calibration with peak calibration weight between 100 and 350 m, we estimate that east equatorial Atlantic SST cooled by ∼5C between the Late Miocene and Pleistocene. On shorter timescales, we use the TEX86 record as a proxy for South Atlantic Central Water (SACW), which originates from surface waters in the South Atlantic Gyre and mixes at depth with Antarctic Intermediate Water (AAIW). Leads and lags around the Pliocene M2 glacial (∼3.3 Ma) in our record, combined with published information, suggest that the M2 glacial was marked by SACW cooling during an austral summer insolation minimum and that decreasing CO2 levels were a feedback, not the initiator, of glacial expansion.

1 Introduction

Accurate tropical sea-surface temperature (SST) reconstructions are needed to assess oceanic heat distribution and ocean–atmosphere circulation in warmer-than-present climates, such as during the Pliocene. Tropical warming, even of a small magnitude, can lead to major changes in atmospheric circulation, with effects including intensification of monsoonal precipitation (Haywood et al., 2020; Zhang et al., 2016) and warming of the extratropics via teleconnections (Barreiro et al., 2006). At present, the most frequently used proxies for past SSTs are Mg/Ca (foraminifer calcite), U37k (alkenones) and TEX86 (glycerol dialkyl glycerol tetraethers; GDGTs). Each proxy has specific confounding factors and calibration issues. For example, Mg/Ca paleothermometry requires corrections for the Mg/Ca ratio of seawater and diagenesis (Dekens et al., 2002; Evans et al., 2016), and the U37k proxy is insensitive to temperatures above ∼28C (Müller et al., 1998). TEX86 is potentially one of the best proxies to reconstruct SSTs above 28 C, but one critical concern is that GDGTs are produced by marine Thaumarchaeota throughout the water column and often dominantly below the mixed layer (Ingalls et al., 2006; Kim et al., 2015; Lengger et al., 2014; Shah et al., 2008).

Although GDGTs may be dominantly produced at depths below the mixed layer, the efficiency of organic carbon export towards the sediment decreases exponentially with depth below the photic zone as shown with sediment trap data and simulations (e.g., Martin et al., 1987; Middelburg, 2019), which for GDGTs is supported by in situ measurements (e.g., Wüchter et al., 2005). However, as the production of GDGTs in the mixed layer is generally low relative to production within and below the thermocline, GDGTs produced below the mixed layer might still comprise a significant component of the sedimentary assemblage (e.g., Taylor et al., 2013; Ho and Laepple, 2017). The export depth of downcore GDGTs can potentially be traced with the fractional abundance of GDGT-2 to GDGT-3 (GDGT[2/3] ratio; Hernández-Sánchez et al., 2014; Pearson et al., 2016; Taylor et al., 2013; Villanueva et al., 2015). GDGT[2/3] values are typically low (<5) in suspended particulate matter (SPM) sampled from shallow (<100 m) water and sharply increase to >25 in deeper water (Hurley et al., 2018). This division is currently best explained by depth differentiation of Thaumarcheotal ecotypes with different cyclization patterns (Kim et al., 2016; Villanueva et al., 2015). Core tops do not have GDGT[2/3] values indicative of a purely “deep” Thaumarcheotal GDGT source (Fig. 1). Instead, they show a gradual transition from low to moderately high values, which can be attributed to the integration of GDGTs from a range of shallow to intermediate depths. Likewise, TEX86 (Eq. 1) may reflect a depth-integrated surface to subsurface signal (SubST) that is sensitive to temperature variability at and below the surface ocean (Kim et al., 2008, 2012; Schouten et al., 2002; Tierney and Tingley, 2015).

(1) TEX 86 = [ GDGT-2 + GDGT-3 + cren ] [ GDGT-1 + GDGT-2 + GDGT-3 + cren ]

Indeed, several multi-proxy paleostudies suggest that dissimilarities between TEX86-based SST reconstructions and other temperature records may be explained by a SubST signal recorded in TEX86 variability (e.g., Lopes dos Santos et al., 2010; McClymont et al., 2012; Rommerskirchen et al., 2011; White and Ravelo, 2020). Together with the GDGT[2/3] distribution in core tops (Fig. 1), these downcore records support the use of a SubST calibration. However, the targeted depth interval and weight distribution of the temperature integration remain the subject of discussion (Ho and Laepple, 2016, 2017; Tierney et al., 2014). The correlation of core top TEX86 to SST is not significantly stronger than to SubST, as obtained from a wide range of depth intervals (Ho and Laepple, 2016; Schouten et al., 2002). Therefore, the sensitivity of sedimentary TEX86 to temperature cannot be inferred from variability in the spatial dimension, but variability in the temporal dimension may provide solutions. Here, we present a new 0–15 Ma TEX86 record from Ocean Drilling Program (ODP) Site 959 in the eastern equatorial Atlantic. We evaluate the TEX86 and GDGT[2/3] evolution on million-year and glacial–interglacial timescales during the late Pliocene. We complement these records with U37k SST estimates and dinoflagellate cyst (dinocyst) assemblages and compare them to the Site 959 benthic δ18O record (Norris, 1998a; van der Weijst et al., 2020) and Mg/Ca-based SST and thermocline temperature records (van der Weijst et al., 2022) to elucidate the local source of downcore TEX86 index variability and, consequently, paleoclimatic change.

Figure 1(a) Depth distribution of GDGT[2/3] ratio in suspended particulate matter (SPM; solid lines) at several stations in the Atlantic Ocean (Hurley et al., 2018) and in core tops (circles) between 65 N and 65 S (Tierney and Tingley, 2015). (b) GDGT[2/3] vs. TEX86, with logarithmic fits for GDGT[2/3]>5 (blue line) and GDGT[2/3]<5 (orange line). (c) Depth-integrated (linear) GDGT[2/3] from SPM data (station 23; Hurley et al., 2018). Integration start depths at 0, 50, 100 and 200 m and continuous end depth on the x axis. Range of typical Site 959 sedimentary GDGT[2/3] values between orange dashed lines and corresponding range of potential integration end depths in the vertical. Note that the assumed linear depth weighing may not be realistic but that an integration start depth at 200 m is incompatible with the sedimentary GDGT[2/3] values at Site 959, regardless of depth weighing.


2 Material and methods

2.1 Site, age model and sampling

ODP Site 959 was drilled during Leg 159 in the Gulf of Guinea, ∼160 km offshore Ghana and Côte D'Ivoire (3.62 N, 2.73 W; 2090 m depth; Mascle et al., 1996). It is presently located below the eastward-flowing Guinea Current, which originates from the North Equatorial Counter Current and the Canary Current (Norris, 1998a). SST currently varies between 25.3 and 28.6 C on a seasonal timescale (Locarnini et al., 2013), and surface salinity ranges between 34.6 and 35.0 (Zweng et al., 2013). The water column is characterized by a shallow thermocline, with the 20 C isotherm depth annually varying between  40 and 60 m in response to coastal upwelling (Locarnini et al., 2013; van der Weijst et al., 2022). A minor upwelling event occurs in boreal winter, and a longer and stronger event occurs in boreal summer (Djakouré et al., 2017; Verstraete, 1992; Wiafe and Nyadjro, 2015).

We studied the 0–160 m interval of Hole C and 162–194 m interval of Hole A. The lithology gradually changes from nannofossil/foraminifer chalk at the bottom to nannofossil/foraminifer ooze at the top (Mascle et al., 1996). The age model is based on nannofossil biostratigraphy (Shin et al., 1998) and benthic δ18O stratigraphy (Wagner, 1998) between 0 and 23 meters below seafloor (m b.s.f.), benthic δ18O stratigraphy between 33 and 46 mbsf (van der Weijst et al., 2020), astronomical tuning of high-resolution X-ray fluorescence data (Vallé et al., 2016) between 51 and 100 mbsf, and a spline function through planktic foraminifer and nannofossil biostratigraphy between 162 and 194 mbsf (Norris, 1998b; Shafik et al., 1998). We plot the data between 21.35 and 97.43 mbsf on the revised meters composite depth (rMCD) scale of Vallé et al. (2016). Because outside this interval a splice is unavailable, we constructed a revised meters below sea floor (rMBSF) scale in coherence with Vallé et al. (2016), assuming ∼5 % loss of sediment between consecutive cores (Supplement).

The core was sampled at a relatively low resolution between 0 and 15 Ma for biomarker analysis and at a higher resolution in the late Pliocene (2.7–3.5 Ma) interval for multi-proxy reconstructions. A total of 219 samples were extracted for TEX86, of which 105 are closely spaced in the 2.7–3.5 Ma interval. U37k data were generated for 79 of these samples, of which 52 are from the 2.7–3.5 Ma interval. Additionally, we performed palynological assessments on 48 late Pliocene splits of samples that were also used for biomarker analysis.

2.2 Biomarker proxies

At Utrecht University, sediment samples were freeze-dried, and the outer surfaces were removed using a clean knife to prevent contamination. The sediment was finely ground using a mortar and pestle. In the process, a split of coarsely ground material was taken for palynology. Biomarkers were extracted at 100 C and 7.6×106 Pa using a dichloromethane (DCM) : methanol (9:1, v/v) solvent mixture in an accelerated solvent extractor (ASE 350, Dionex). The total lipid extract was separated over an activated Al2O3 column into apolar, ketone and polar fractions using hexane : DCM (9:1, v/v), hexane : DCM (1:1, v/v) and DCM : methanol (1:1, v/v) as eluents, respectively. For GDGT analysis, a C46 GDGT (99 ng) standard was added to the polar fraction before it was filtered over a 0.45 µm PTFE filter and dissolved in hexane : isopropanol 99:1 (v/v). Measurements were performed on an ultrahigh-performance liquid chromatograph–mass spectrometer (UHPLC-MS) using the method of Hopmans et al. (2016). Samples with a branched to isoprenoid tetraether (BIT) index of >0.3 could be influenced by terrestrially produced GDGTs (Hopmans et al., 2004; Weijers et al., 2006) and were excluded from further analysis (n=3). All samples were within the reliable range of the methane index (MI; Zhang et al., 2011) and ring index (RI; Zhang et al., 2016), suggesting no appreciable inputs from methanogenic archaea or other non-thermal factors. We will present our results using various calibrations, including the conventional approach of calibrating TEX86 to SST, but also using various models using the relation between TEX86 and various depth integrations.

The ketone fraction was dissolved in hexane and analyzed on a gas chromatograph (GC) coupled to a flame ionization detector (GC-FID, Hewlett Packard 6890 series) equipped with a CP-Sil 5 fused silica capillary column (25 m×0.32 mm; film thickness 0.12 µm) and a 0.53 mm precolumn. Samples were injected on-column at 70 C with helium as a carrier gas and a flow rate of 2 mL min−1. The oven program was as follows: 70 C for 1 min, then ramped to 130 C at 20 C min−1, then to 320 C at 4 C min−1, and then held isothermal for 10 min. U37k values were calculated from the fractional abundances of the C37:2 and C37:3 alkenones following Prahl and Wakeham (1987) and calibrated to SST using the calibration of Müller et al. (1998). We chose not to adopt the Bayspline approach (Tierney and Tingley, 2018) or other tropical calibrations (e.g., Sonzogni et al., 1997, DSR II) as they may amplify noise into signal because of difficulty in accurately determining the very small C37:3 peak areas at the high range of SST (Herbert et al., 2020).

2.3 Palynology

Coarsely crushed freeze-dried samples were spiked with a known amount of Lycopodium clavatum spores to quantify absolute palynomorph (Stockmarr, 1972) and treated with 30 % HCl () and 40 % cold HF () to remove carbonates and silicates. The residue was sieved over nylon mesh, and from the 10–250 µm fraction, the lighter organic fraction was separated from the heavy mineral fraction by suspension. The remaining palynological residue was mounted on glass microscopic slides with glycerine jelly. Samples were analyzed using an optical microscope under 400× magnification and approximately 300 dinoflagellate cysts (dinocysts) were counted per sample. Taxonomy follows that of Williams et al. (2017). Capisocysta was only identified at the genus level, because known Capisocysta species are best identified by the number of antapical plates (Head, 1998), which is complicated by the common dissociation of the hypocystal plates. We calculate the relative abundances of a taxon based on the total dinocyst sum. Dinocysts with affinities for cold-temperate waters in the modern ocean (Boessenkool et al., 2001; Zonneveld et al., 2013) were grouped as “cool species” and used for estimating surface water temperature. This species are Ataxodinium choane, Corrudinium devernaliae, Corrudinium harlandii, Impagidinium pallidum, Nematosphaeropsis labyrinthus and Pyxidinopsis reticulata. We also use the P/(P+G) ratio (Versteegh, 1994) as a paleoproductivity index. The P/(P+G) ratio quantifies the number of cysts produced by heterotrophic dinoflagellates, which are characterized by peridinioid-type tabulation (P) over the total number of dinocysts, including those produced by autotrophic or mixotrophic species with a gonyaulacoid-type tabulation (G).

3 Results

3.1 TEX86, U37k and GDGT[2/3] evolution at Site 959

We explore our TEX86 results using several calibration approaches in Fig. 2, including an SST calibration using the exponential TEX86H (Kim et al., 2010). It should be noted that this calibration yields significant residuals in tropical regions and suffers from regression dilution (Tierney and Tingley, 2014). We also use a linear calibration model using a Bayesian method, a model specific for the tropical Atlantic and models assuming contributions from GDGTs produced below the mixed layer. Below, we first explore our TEX86 record in the context of other GDGT-based indicators and other data using the TEX86H calibration for SST (Figs. 3–7), while subsurface calibrations are discussed in Figs. 7 and 8.

Figure 2Site 959 TEX86 record, calibrated using TEX86H and alternative SST calibrations and with the HL16-SubST calibration (Ho and Laepple, 2015), compared to the U37k record. Late Pliocene glacial–interglacial variability in panel (a) and 15 Myr evolution in panel (b).


Between 15 and 11 Ma, TEX86H-based SST values fluctuate around an average value of ∼29C (Fig. 3), similar to the present-day local annual maximum surface temperature of 28.6 C. A long-term gradual cooling starts around 11 Ma and is punctuated by two major cooling steps of ∼2C each, around 4.9 and 3.4 Ma, in which TEX86H drops well below the present-day local annual minimum surface temperature of 25.3 C. The cooling trend seems to end in the late Pleistocene, where TEX86H records a ∼4C warming between 0.5 Ma and present, with a core top TEX86H of 24.0 C. U37k is at saturation (∼28C) for most of the interval but registers some colder temperatures in the late Pliocene (27–28 C) and drops below 27 C after 1.8 Ma. The core top sample registers a temperature of 26.6 C, which is within calibration error (1.5 C; Müller et al., 1998) of the 27.5 C modern mean annual SST (Locarnini et al., 2013). The GDGT[2/3] ratio is 7.0 in the core top. In the oldest part of the studied interval, GDGT[2/3] values vary between 4 and 5.5. Following an abrupt increase at ∼13.6 Ma, GDGT[2/3] varies largely between 5 and 8. After 1.7 Ma, GDGT[2/3] persistently increases to >7.

3.2 Pliocene dinocyst assemblages

Palynological assemblages are dominated by dinocysts with only minor amounts of terrestrially derived pollen and spores. Dinocyst preservation is generally good. The assemblages are dominated by Brigantedinium spp. (19 %–60 %), closely followed by Spiniferites spp. (8 %–36 %) (Fig. 4; see the Supplement for a complete overview of the assemblages). The interval between 3.33 and 3.16 Ma was sampled at a higher resolution (∼5 kyr on average) and displays strong variability of the negatively correlated Brigantedinium spp. and Spiniferites spp. on timescales <10 kyr. Other relatively abundant and consistently present taxa are Operculodinium spp., Pentapharsodinium dalei and Impagidinium spp. Peak abundances of Lingulodinium spp. were recorded at ∼3.26 Ma (Lingulodinium machaerophorum) and at 2.77 Ma (Lingulodinium hemicystum). Capisocysta sp. is consistently present in the lower part of the interval, with peak abundances (>20 %) at 3.32 and 3.20 Ma, but has its last occurrence at 3.17 Ma. There is general coherence between TEX86H-SST variability and the relative abundance of the “cool species” group, which is dominated by N. labyrinthus. On average, the dinocysts with affinities for colder-temperate waters make up 6 % of the total assemblage, but peak abundances of 10 %–13 % occur during the M2 glacial. The P/(P+G) paleoproductivity index is 0.48 on average, with minimum and maximum values of 0.23 and 0.70, respectively (Fig. 4). P/(P+G) values are primarily driven by the abundance of Brigantedinium spp. The P/(P+G) ratio shows no significant relationship with SST as derived from TEX86H (R2=0.05, p=0.14 in linear regression analysis).

Figure 3New TEX86H-SST (Kim et al., 2010 calibration), U37k-SST (Müller et al., 1998 calibration) and GDGT[2/3] records at Site 959, compared to benthic δ18O records from Norris (1998a) and van der Weijst et al. (2020). Shaded interval and arrow indicate annual SST range and mean annual SST in the Gulf of Guinea, respectively (Locarnini et al., 2013).


Figure 4Relative abundances of the major groups of dinoflagellate cysts, the sum of dinocysts with cold affinities (“Cool species”) and P/(P+G) paleoproductivity index compared to TEX86H-SST. Shaded bands indicate TEX86H-SST minima during the late Pliocene M2, KM2 and G20 glacials.


4 Discussion

4.1 Identifying the source depth of the TEX86 signal

The 15-million-year TEX86H record at Site 959 shows a cooling from >30C between 11 and 15 Ma to <20C between 0 and 3.3 Ma (Fig. 3). The U37k record, on the other hand, is near or at proxy saturation (>27C) for most of the interval, with some lower values after 1.8 Ma. Assuming that both proxies reflect SST, the relatively small offset between U37k and TEX86H-SST offset in the core top sample (2.7 C) could be explained by a combination of calibration errors and seasonality. However, the average downcore offset between TEX86H and U37k is 5.4 C and peaks at 7–8 C during the late Pliocene interglacials (Figs. 5 and 6), which demonstrates that TEX86 underestimates past SST at this site, assuming that U37k is a more accurate representation of SST. This could potentially be resolved by using an alternative regional SST calibration. The BAYSPAR-SST calibration (Tierney and Tingley, 2014) accounts for spatial variability of the TEX86-SST relationship and is the only available calibration that raises absolute SST estimates. However, similar to TEX86H (Kim et al., 2010) and low-latitude Atlantic (Zhang et al., 2018) calibrations, it likely overestimates SST variability (Fig. 2), indicating that the calibration slope is too steep. This is potentially a problem inherent to the calibration of TEX86 to SST and could reflect downcore TEX86 sensitivity to SubST instead of SST variability (Ho and Laepple, 2016). In the following sections, we explore multi-proxy data from Site 959 in search of indications for the water depth the downcore TEX86 signal is reflecting.

Figure 5Compilation of late Pliocene multi-proxy temperature data and benthic δ18O at Site 959. Mg/Ca(G. ruber) (pink) and Mg/Ca(N. dutertrei) (yellow) from van der Weijst et al. (2022), benthic δ18O from van der Weijst et al. (2020), and LR04 δ18O stack (Lisiecki and Raymo, 2005) for reference. Dashed lines are LOESS-smoothed trends, and vertical bands indicate the M2, KM2 and G20 glacial stages.


4.1.1 GDGT[2/3] values

A substantial contribution of GDGTs from below the oceanic surface layer is supported by the relatively high GDGT[2/3] ratio, which fluctuates between 5 and 8 in the majority of the record (Fig. 3). Such values are rarely observed in modern-day shallow water (<100 m) SPM and surface sediments (e.g., Besseling et al., 2019; Hernández-Sánchez et al., 2014; Hurley et al., 2018). These values are best explained as a mixed signal of GDGTs from the upper 200 m and intermediate waters (Fig. 1c), i.e., likely dominantly from below the mixed layer (∼50 m) at Site 959. A cross-plot of TEX86H-SST and GDGT[2/3] (Fig. 7) shows that most data are clustered around a positive regression slope, but the oldest (>13.6 Ma) and youngest (<1.7 Ma) data plot outside this main cluster, potentially signaling systematic changes in GDGT production and export depth. Both 13.6 and 1.7 Ma mark the onset of a strong GDGT[2/3] increase (Fig. 3), signaling systematic deepening of the source of the sedimentary GDGT signal.

The ratio of GDGT[2] to GDGT[3] in samples does not directly affect the TEX86 index because they are both included in the denominator and divisor (Eq. 1). If the abundance of GDGT[2] increases relative to that of GDGT[3] with export depth, then samples with a higher GDGT[2/3] ratio should be associated with lower TEX86 values. The positive relationship between TEX86 and GDGT[2/3] (Fig. 7) is therefore an interesting feature. Taylor et al. (2013) found a similar relation in the calibration dataset and in the paleo-domain. Based on their analyses, they conclude that water depth is the dominant control of GDGT[2/3] ratios in sediments, recently supported by SPM analyses (Hurley et al., 2018). They also conclude that the positive correlation between SST and GDGT[2/3], notably the opposite relation to that expected based on homeoviscous adaptation, is an oceanographical artifact. Rather, the positive relation suggests that the long-term TEX86 trends are not driven by GDGT export depth but reflect overall temperature changes instead. Alternatively, non-thermal factors such as water column oxygenation and nutrient supply may have influenced GDGT cyclization: high nutrient availability and ammonia oxidation rates have been linked to low TEX86 values (Hurley et al., 2016; Park et al., 2018). However, following the GDGT[2/3] shift at 13.6 Ma, TEX86H does not systematically change until ∼11 Ma (Fig. 3). Moreover, late Pliocene dinocyst assemblages (Fig. 4) are characterized by a highly variable P/(P+G) ratio, indicating that upwelling and nutrient supply were highly variable on sub-Milankovitch timescales. In contrast to TEX86H (Fig. 5), the P/(P+G) ratio shows no glacial–interglacial variability. It is therefore unlikely that TEX86 variability was primarily driven by non-thermal factors.

Figure 6Close-up of G20, KM2 and M2 glacial stages in 100 kyr windows. Multi-proxy temperature data and benthic δ18O as in Fig. 5. The resolution of the dinocyst record was too low for meaningful comparisons in KM2 and M2 glacials.


4.1.2 Surface, thermocline and sub-thermocline temperature variability

The late Pliocene TEX86H record follows a different evolution than SST records of U37k and the recent Mg/Ca of the surface-dwelling Globigerinoides ruber (Mg/Ca(G. ruber)) and thermocline-dwelling Neogloboquadrina dutertrei (Mg/Ca(N. dutertrei)) generated on the same samples (Fig. 5). The Mg/Ca records (van der Weijst et al., 2022) are based on very well-preserved foraminifera and apply a constant Pliocene correction for the Mg/Ca of seawater, a species-specific core top calibration with a constant partitioning coefficient and embedded dissolution correction (Dekens et al., 2002). Whereas TEX86H decreases by ∼2C, U37k and Mg/Ca-based SSTs show no significant net change in this interval. Moreover, although Mg/Ca and U37k register some cooling during the late Pliocene glacials, their variability seems not primarily driven by glacial–interglacial cyclicity, as is best observed in the high-resolution interval from M2 to KM2 (Fig. 5). This is also true for the relative abundance of dinocysts with affinities for colder temperatures (Figs. 4 and 5). Furthermore, around the KM2 and M2 glacials, TEX86H cooling leads to increasing abundances of cool dinocysts, as well as Mg/Ca(G. ruber) and U37k cooling by  5–10 ka (Fig. 6), which underscores the independent evolution of these records.

At sites where TEX86 is suspected to be affected by GDGTs produced below the surface ocean, TEX86H has occasionally been interpreted to reflect thermocline temperature variability (e.g., Lopes dos Santos et al., 2010). Thermocline temperatures are lower and can be highly variable through time, which potentially explains both absolute SST underestimation and amplitude overestimation in specific TEX86H records. However, both could also be explained as a general artifact of an overestimated calibration slope that results from calibrating a depth-integrated GDGT signal to SST (Ho and Laepple, 2015, 2016). Our TEX86H record shows weak coherence with the Site 959 Mg/Ca(N. dutertrei) thermocline temperature record (Fig. 5; van der Weijst et al., 2022). Whereas TEX86 registers late Pliocene cooling, Mg/Ca(N. dutertrei) registers warming in relation to thermocline deepening (van der Weijst et al., 2022). Moreover, while the range of Mg/Ca(N. dutertrei) temperatures is similar to TEX86H on glacial–interglacial timescales, Mg/Ca(N. dutertrei) does not follow the glacial–interglacial cyclicity as recorded in the TEX86.

In contrast, the variability in the TEX86H-SST and benthic δ18O records at Site 959 is very similar (Figs. 3 and 5). Benthic δ18O is a faithful recorder of glacial–interglacial cyclicity because it registers a combined signal of deep ocean temperature and polar ice sheet volume. Could TEX86 at Site 959 also be sensitive to glacial–interglacial cyclicity in deeper waters? At the onset of the M2 and KM2 glacials, the glacial expression of TEX86H leads δ18O by ∼5 kyr (Fig. 6), which indicates that, despite the similarities between TEX86H-SST and benthic δ18O, these proxies record variability in different water masses. Bottom waters at Site 959 are predominantly ventilated by North Atlantic Deep Water (NADW), whereas depths below the thermocline are occupied by South Atlantic Central Water (SACW; Fig. 8). SACW originates from surface waters in the South Atlantic gyre, which at depth mixes with Antarctic Intermediate Water (AAIW). It is therefore sensitive to climate variability in the mid-latitudes in the Southern Hemisphere and, like benthic δ18O, sensitive to higher-latitude climate change. Although our record samples SACW, AAIW dynamics have previously been suggested to be the main driver of TEX86 variability in the Arabian Sea (Huguet et al., 2006) and in the southeast Atlantic (Rommerskirchen et al.,  2011). If TEX86 at Site 959 records a mixed signal from surface and subsurface waters, as is supported by the relatively high GDGT[2/3] ratio (Figs. 1 and 3), it is likely sensitive to SACW dynamics, which explains the strong coherence with the benthic δ18O record. Collectively, these observations suggest that the downcore TEX86 record at Site 959 is substantially affected by temperature variability below the surface ocean and that calibration to SubST is more appropriate than to SST.

Figure 7Cross-plot of TEX86H-SST and GDGT[2/3]. The oldest (>13.6 Ma; dark blue) and youngest (<1.7 Ma; dark red) samples largely stand out from the majority of the data. Linear regression lines determined from full dataset (stippled line) and samples dated between 13.6 and 1.7 Ma (solid line; 88 % samples, R2=0.37, p<0.001).


4.2 Calibration of the depth-integrated TEX86 record

Several TEX86 calibrations for subsurface temperature (TEX86-SubST) are available, with each assuming, among other things, a different depth integration of the water column (Fig. 8). The BAYSPAR-SubST (Tierney and Tingley, 2014) and TEX86H-SubST (Kim et al., 2012) calibrations both target the upper 200 m of the water column, but the weight of the BAYSPAR-SubST calibration is centered at relatively shallow depths compared to the linearly weighted TEX86H-SubST calibration (Fig. 8). The HL16-SubST (Ho and Laepple, 2016) calibration targets the upper 950 m. The validity of this depth interval was questioned (Tierney et al., 2017), but with peak calibration weight at 100–350 m (Fig. 8), it is compatible with the documented vertical distribution of Thaumarcheotal cell counts and GDGT concentrations (e.g., Hernández-Sánchez et al., 2014; Schouten et al., 2012; Wuchter et al., 2005; Sintes et al., 2016), leaving export efficiency as the major uncertainty. Moreover, based on the large proportion of core tops with higher (>5) GDGT[2/3] values in the global calibration set (Fig. 1), a slightly deeper calibration target (including depths below 200 m) is arguably preferable. This deeper calibration target lowers the slope of the calibration, which dampens reconstructed Miocene–Pliocene temperature variability at Site 959 (Fig. 2).

Figure 8(a) Temperature and salinity distribution along a depth transect in the modern Gulf of Guinea (numbers: meters below sea surface), data from Locarnini et al. (2013) and Zweng et al. (2013). (b) Vertical weight distribution of the discussed SubST calibrations (silhouettes) and the water column composition in the modern Gulf of Guinea. (c) Site 959 TEX86H-SubST evolution according to different SubST calibrations, compared to global mean temperature evolution (Tierney et al., 2020). Δ Temperature estimated from linear trends (dashed lines) between the Middle Miocene and Pleistocene.


The ratio between temperature change in the surface and subsurface ocean is 1:1 when averaged across many sites and on longer timescales (Ho and Laepple, 2016). Under the assumption that the (integrated) depth of GDGT sourcing and the structure of the water column was stable, a depth-integrated TEX86 record might still be of value to assess SST variability. Depending on the calibration, the magnitude of Late Miocene to Pleistocene cooling at Site 959 is 10 C (BAYSPAR-SubST), 7 C (TEX86H-SubST) or 5 C (HL16-SubST; Fig. 8). Assuming a 1:1 ratio, the long-term SST trend should be of a similar magnitude. Currently available tropical SST estimates on the studied timescale are either compromised by saturation of the U37k proxy or based on TEX86 and are therefore not suitable for independent comparison. However, global mean surface temperature estimates based on benthic δ18O data suggest a ∼5C cooling across the same interval (Hansen et al., 2013; Tierney et al., 2020). In the absence of major oceanographic changes, it is unlikely that the local cooling trend at Site 959 was larger, because temperature variability is generally lower in the tropics compared to high latitudes where deep-ocean waters derive from (i.e., polar amplification). This suggests that the BAYSPAR-SubST and TEX86H-SubST calibrations may overestimate the magnitude of long-term cooling (Fig. 8). From the currently available calibrations, HL16-SubST best approximates both expected glacial–interglacial variability and multi-million year trends (Figs. 2 and 8) and could potentially be used in combination with other SST records and/or modern water column data to correct for the local SST-SubST offset (Fig. 2). It is possible that the 1:1 ratio between SST and SubST did not persist across the entire record, e.g., due to changes in SACW production. To further improve tropical SST estimates in warmer climates, it should be further explored under which conditions a TEX86 record that is driven by subsurface variability can be used to reconstruct SST variability and how to correct for the surface–subsurface temperature offset.

4.3 Late Neogene TEX86 as an intermediate ocean signal at Site 959: new insights into M2 glacial inception

It is currently unclear what caused M2 glacial inception, but hypothesized mechanisms include declining atmospheric CO2 levels (Berends et al., 2019; Dolan et al., 2015) and reduced latitudinal heat transport by the North Atlantic Current (De Schepper et al., 2009, 2013) in the Northern Hemisphere (NH) and by Indonesian throughflow (De Vleeschouwer et al., 2018) in the Southern Hemisphere (SH). Multi-proxy data from Site 999 in the Caribbean Sea show that declining CO2 levels lagged the benthic δ18O glacial expression at the onset of M2 (de la Vega et al., 2020). This indicates that declining CO2 levels did not initiate the M2 glacial, although they might have affected its intensity and duration.

Figure 9Exploring the chronology around the M2 glacial with records aligned on peak glacial benthic δ18O values (original age models in Fig. S2). (a) Daily insolation at 60 N/S during NH/SH summer solstice (Laskar, 1990). (b) East equatorial Atlantic SST at Site 959 (this study) and Site 662 (Herbert et al., 2010) and Caribbean Sea SST at Site 999 (de la Vega et al., 2020). (c) Site 959 TEX86 HL16-SubST, indicative of intermediate ocean temperatures. (d) Atmospheric CO2 reconstructions from Site 999 (de la Vega et al., 2020). (e)  Benthic δ18O records of Sites 959 (van der Weijst et al., 2020), 662 (Lisiecki and Raymo, 2005) and 999 (Haug et al., 2001; de la Vega et al., 2020), aligned to LR04 (Lisiecki and Raymo, 2005) on peak M2 values.


Because the TEX86 record at Site 959 is affected by SACW, it can potentially be used to study the connection between our low-latitude site and higher latitudes. The TEX86 and benthic δ18O records at Site 959 correspond closely but are out of sync during the M2 glacial (Fig. 6). Both records, generated on the same samples, display a double peak around M2, but TEX86 leads δ18O by ∼5 kyr. Whereas the depth-integrated TEX86 signal at Site 959 is sensitive to SH conditions through SACW, bottom waters at Site 959 consisted predominantly of NADW sourced from the NH during the late Pliocene (van der Weijst et al., 2020). The TEX86 lead therefore indicates a lead of SH over NH cooling/glaciation at the onset of M2, in agreement with transient ice sheet simulation (Berends et al., 2019).

We further explore the relative chronology around M2 in Fig. 9. The age models were aligned with peak glacial δ18O values, and the same records are plotted in the original age models in Fig. S1 in the Supplement. Although it is difficult to discriminate between signal and noise at this resolution, some patterns seem to emerge from this compilation. Directly before M2, benthic δ18O, SubST and SST records from the eastern equatorial Atlantic (Sites 959 and 662) and Caribbean Sea (Site 999) indicate gradual, widespread cooling, while CO2 levels were rising. A sharp drop in TEX86H-SubST indicates intermediate ocean cooling at the onset of M2, leading pronounced glacial cooling (benthic δ18O) of the deep ocean and tropical surface ocean (best observed at Site 662). According to the chronology in Fig. 9, this intermediate ocean cooling occurred during a SH summer insolation minimum, suggesting that astronomical forcing may have amplified Southern Ocean cooling and tipped the system to glacial conditions. The sensitivity of the Antarctic ice sheet to astronomical forcing is also reflected by the 3.3 Ma onset of ∼20 kyr cyclicity in ice-rafted debris (IRD) at East Antarctic IODP Site U1361, with peak IRD mass accumulation rates at SH summer insolation minima (Patterson et al., 2014). Based on the chronology of Fig. 9, SH cooling led atmospheric CO2 drawdown by  15–20 kyr. Lags of a similar magnitude were not unusual during Pleistocene glacial inceptions (Petit et al., 1999). Marginal Antarctic diatom assemblages and lithology suggest that considerable sea ice expansion occurred during M2, which could have inhibited CO2 outgassing from the deep ocean (Ishino and Suto, 2020; McKay et al., 2012). Moreover, a slight increase in dust-mediated iron fertilization of the Southern Ocean may have facilitated export productivity (Martínez-Garcia et al., 2011). In combination with reduced mixing between Northern Component Water and Southern Component Water (van der Weijst et al., 2020), these processes would have led to increased carbon storage in the deep ocean in response to M2 cooling. Higher-resolution records and correlations would be required to fully test this hypothesis.

5 Conclusions

Several lines of evidence show that the TEX86 record at Site 959 in the eastern equatorial Atlantic is best explained as a depth-integrated signal that is substantially affected by temperature variability below the thermocline. Relatively high GDGT[2/3] ratios since ∼13.6 Ma indicate that the sedimentary GDGTs were partly sourced from deeper (>200 m) waters. Moreover, late Pliocene multi-proxy data show that TEX86 variability is predominantly driven by glacial–interglacial variability, in contrast to other SST and thermocline records (U37k and Mg/Ca(G. ruber) SST and Mg/Ca(N. duterrei)) and dinocyst abundances. The TEX86 record strongly resembles the benthic δ18O record, indicating that it dominantly records temperature changes in a sub-thermocline water mass. At Site 959, intermediate depths are occupied by SACW, which derives from the mid-latitudes of the South Atlantic. According to the present-day water column composition, a substantial contribution of GDGTs from deeper (>200) waters is needed to obtain TEX86 sensitivity to SACW. This favors a subsurface calibration that integrates across a wider range of water depths, such as HL16-SubST (Ho and Laepple, 2016), in which calibration weight peaks 100–350 m. This calibration target interval is compatible with pelagic Thaumarcheotal cell counts and GDGT concentrations and with core top GDGT[2/3] values in the global calibration set.

Even if TEX86 mainly reflects SubST variability, it may be used to reconstruct past SSTs if the temporal relationship between SST–SubST in a certain region is well-understood. Assuming a 1:1 relationship between long-term SubST and SST trends, our TEX86 record suggests 5 C tropical SST cooling between the Late Miocene and Pleistocene when calibrated with the Ho and Laepple (2016) subsurface calibration. Additionally, TEX86 is also highly informative as a SubST proxy, because the sensitivity of TEX86 to SACW at Site 959 offers a chance to explore connections between high- and low-latitude climate variability. The TEX86 and benthic δ18O records at Site 959 are highly similar, but TEX86 is sensitive to high-latitude SH conditions through AAIW, whereas the benthic δ18O record is primarily forced by high-latitude NH conditions through NADW. A ∼5 kyr lead of TEX86 relative to δ18O at the onset of the M2 glacial stage indicates a lead of SH over NH cooling. Albeit limited by the resolution of the datasets, multi-proxy data from Site 959, Site 999 (Caribbean Sea) and Site 662 (eastern equatorial Atlantic) aligned based on LR04 peak glacial δ18O values indicate that SH cooling also led tropical ocean cooling and CO2 levels by 15–20 kyr, suggesting that CO2 drawdown was a consequence of glacial conditions. Instead, glacial expansion and SACW cooling at SH insolation minima suggest that orbital forcing played a pivotal role in M2 glacial inception.

Data availability

New Site 959 data are available as the Supplement to this paper and have been submitted to the PANGAEA online data repository.


The supplement related to this article is available online at:

Author contributions

CMHvdW, FP, FS and AS designed the study. CMHvdW, KJvdL and TJTV generated the data. All authors contributed to interpretations. CMHvdW wrote the paper with input from FP, FS and AS and incorporated feedback from all authors.

Competing interests

At least one of the (co-)authors is a member of the editorial board of Climate of the Past. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


We thank the International Ocean Discovery Program and the predecessor drilling programs for samples and data. We thank Arnold van Dijk, Giovanni Dammers, Natasja Welters, Klaas van Nierop and Dominika Kasjaniuk (Utrecht University), and Wim Boer (NIOZ) for analytical support. We thank the two reviewers for their thorough reviews and valuable insights.

Financial support

This work was carried out under the program of the Netherlands Earth System Science Centre (NESSC), financially supported by the Ministry of Education, Culture and Science (OCW) (grant no. 024.002.001). Appy Sluijs was funded by the European Research Council under consolidator grant no. 771497.

Review statement

This paper was edited by Arne Winguth and reviewed by two anonymous referees.


Barreiro, M., Philander, S. G., Pacanowski, R., and Fedorov, A. V.: Simulations of warm tropical conditions with application to middle Pliocene atmospheres, Clim. Dynam., 26, 349–365,, 2006. 

Berends, C. J., de Boer, B., Dolan, A. M., Hill, D. J., and van de Wal, R. S. W.: Modelling ice sheet evolution and atmospheric CO2 during the Late Pliocene, Clim. Past, 15, 1603–1619,, 2019. 

Besseling, M. A., Hopmans, E. C., Koenen, M., van der Meer, M. T. J., Vreugdenhil, S., Schouten, S., Sinninghe Damsté, J. S., and Villanueva, L.: Depth-related differences in archaeal populations impact the isoprenoid tetraether lipid composition of the Mediterranean Sea water column, Org. Geochem., 135, 16–31,, 2019. 

Boessenkool, K. P., Van Gelder, M. J., Brinkhuis, H., and Troelstra, S. R.: Distribution of organic-walled dinoflagellate cysts in surface sediments from transects across the Polar Front offshore Southeast Greenland, J. Quaternary Sci., 16, 661–666,, 2001. 

Dekens, P. S., Lea, D. W., Pak, D. K., and Spero, H. J.: Core top calibration of Mg/Ca in tropical foraminifera: Refining paleotemperature estimation, Geochem. Geophy. Geosy., 3, 1–29,, 2002. 

de la Vega, E., Chalk, T. B., Wilson, P. A., Bysani, R. P., and Foster, G. L.: Atmospheric CO2 during the Mid-Piacenzian Warm Period and the M2 glaciation, Sci. Rep., 10, 11002,, 2020. 

De Schepper, S., Head, M. J., and Groeneveld, J.: North Atlantic Current variability through marine isotope stage M2 (circa 3.3 Ma) during the mid-Pliocene, Paleoceanography, 24, PA4206,, 2009. 

De Schepper, S., Groeneveld, J., Naafs, B. D. A., Van Renterghem, C., Hennissen, J., Head, M. J., Louwye, S., and Fabian, K.: Northern Hemisphere Glaciation during the Globally Warm Early Late Pliocene, edited by: Smith, V. C., PLoS One, 8, e81508,, 2013. 

De Vleeschouwer, D., Auer, G., Smith, R., Bogus, K., Christensen, B., Groeneveld, J., Petrick, B., Henderiks, J., Castañeda, I. S., O'Brien, E., Ellinghausen, M., Gallagher, S. J., Fulthorpe, C. S. and Pälike, H.: The amplifying effect of Indonesian Throughflow heat transport on Late Pliocene Southern Hemisphere climate cooling, Earth Planet. Sc. Lett., 500, 15–27,, 2018. 

Djakouré, S., Penven, P., Bourlès, B., Koné, V., and Veitch, J.: Respective roles of the Guinea current and local winds on the coastal upwelling in the northern Gulf of Guinea, J. Phys. Oceanogr., 47, 1367–1387,, 2017. 

Dolan, A. M., Haywood, A. M., Hunter, S. J., Tindall, J. C., Dowsett, H. J., Hill, D. J., and Pickering, S. J.: Modelling the enigmatic Late Pliocene Glacial Event – Marine Isotope Stage M2, Global Planet. Change, 128, 47–60,, 2015. 

Evans, D., Brierley, C. M., Raymo, M. E., Erez, J., and Müller, W.: Planktic foraminifera shell chemistry response to seawater chemistry: Pliocene–Pleistocene seawater Mg/Ca, temperature and sea level change, Earth Planet. Sc. Lett., 438, 139–148,, 2016. 

Hansen, J., Sato, M., Russell, G., and Kharecha, P.: Climate sensitivity, sea level and atmospheric carbon dioxide, Philos. T. Roy. Soc. A, 371, 20120294,, 2013. 

Haug, G. H., Tiedemann, R., Zahn, R., and Ravelo, A. C.: Role of Panama uplift on oceanic freshwater balance, Geology, 29, 207–210,<0207:ROPUOO>2.0.CO;2, 2001. 

Haywood, A. M., Tindall, J. C., Dowsett, H. J., Dolan, A. M., Foley, K. M., Hunter, S. J., Hill, D. J., Chan, W.-L., Abe-Ouchi, A., Stepanek, C., Lohmann, G., Chandan, D., Peltier, W. R., Tan, N., Contoux, C., Ramstein, G., Li, X., Zhang, Z., Guo, C., Nisancioglu, K. H., Zhang, Q., Li, Q., Kamae, Y., Chandler, M. A., Sohl, L. E., Otto-Bliesner, B. L., Feng, R., Brady, E. C., von der Heydt, A. S., Baatsen, M. L. J., and Lunt, D. J.: The Pliocene Model Intercomparison Project Phase 2: large-scale climate features and climate sensitivity, Clim. Past, 16, 2095–2123,, 2020. 

Head, M. J.: New goniodomacean dinoflagellates with a compound hypotractal archeopyle from the late cenozoic: Capisocysta Warny and Wrenn, emend., J. Paleontol., 72, 797–809,, 1998. 

Herbert, T. D., Peterson, L. C., Lawrence, K. T., and Liu, Z.: Tropical ocean temperatures over the past 3.5 million years, Science, 328, 1530–1534,, 2010. 

Hernández-Sánchez, M. T., Woodward, E. M. S., Taylor, K. W. R., Henderson, G. M., and Pancost, R. D.: Variations in GDGT distributions through the water column in the South East Atlantic Ocean, Geochim. Cosmochim. Ac., 132, 337–348,, 2014. 

Ho, S. L. and Laepple, T.: Glacial cooling as inferred from marine temperature proxies TEX86H and U37K, Earth Planet. Sci. Lett., 409, 15–22,, 2015. 

Ho, S. L. and Laepple, T.: Flat meridional temperature gradient in the early Eocene in the subsurface rather than surface ocean, Nat. Geosci., 9, 606–610,, 2016. 

Ho, S. L. and Laepple, T.: Reply to “Eocene temperature gradients,” Nat. Geosci., 10, 539–540,, 2017. 

Hopmans, E. C., Weijers, J. W. H., Schefuß, E., Herfort, L., Sinninghe Damsté, J. S., and Schouten, S.: A novel proxy for terrestrial organic matter in sediments based on branched and isoprenoid tetraether lipids, Earth Planet. Sc. Lett., 224, 107–116,, 2004. 

Hopmans, E. C., Schouten, S., and Sinninghe Damsté, J. S.: The effect of improved chromatography on GDGT-based palaeoproxies, Org. Geochem., 93, 1–6,, 2016. 

Huguet, C., Kim, J. H., Damsté, J. S. S., and Schouten, S.: Reconstruction of sea surface temperature variations in the Arabian Sea over the last 23 kyr using organic proxies (TEX86 and U37K), Paleoceanography, 21, PA3003,, 2006. 

Hurley, S. J., Elling, F. J., Könneke, M., Buchwald, C., Wankel, S. D., and Santoro, A. E.: Influence of ammonia oxidation rate on thaumarchaeal lipid composition and the TEX86 temperature proxy, P. Natl. Acad. Sci. USA, 113, 7762–7767,, 2016. 

Hurley, S. J., Lipp, J. S., Close, H. G., Hinrichs, K. U., and Pearson, A.: Distribution and export of isoprenoid tetraether lipids in suspended particulate matter from the water column of the Western Atlantic Ocean, Org. Geochem., 116, 90–102,, 2018. 

Ingalls, A. E., Shah, S. R., Hansman, R. L., Aluwihare, L. I., Santos, G. M., Druffel, E. R. M., and Pearson, A.: Quantifying archaeal community autotrophy in the mesopelagic ocean usinq natural radiocarbon, P. Natl. Acad. Sci. USA, 103, 6442–6447,, 2006. 

Ishino, S. and Suto, I.: Late Pliocene sea-ice expansion and its influence on diatom species turnover in the Southern Ocean, Mar. Micropaleontol., 160, 101895,, 2020. 

Kim, J.-H., Schouten, S., Hopmans, E. C., Donner, B., and Sinninghe Damsté, J. S.: Global sediment core-top calibration of the TEX86 paleothermometer in the ocean, Geochim. Cosmochim. Ac., 72, 1154–1173,, 2008. 

Kim, J.-H., van der Meer, J., Schouten, S., Helmke, P., Willmott, V., Sangiorgi, F., Koç, N., Hopmans, E. C., and Damsté, J. S. S.: New indices and calibrations derived from the distribution of crenarchaeal isoprenoid tetraether lipids: Implications for past sea surface temperature reconstructions, Geochim. Cosmochim. Ac., 74, 4639–4654,, 2010. 

Kim, J.-H., Romero, O. E., Lohmann, G., Donner, B., Laepple, T., Haam, E., and Sinninghe Damsté, J. S.: Pronounced subsurface cooling of North Atlantic waters off Northwest Africa during Dansgaard–Oeschger interstadials, Earth Planet. Sc. Lett., 339–340, 95–102,, 2012. 

Kim, J.-H., Schouten, S., Rodrigo-Gámiz, M., Rampen, S. W., Marino, G., Huguet, C., Helmke, P., Buscail, R., Hopmans, E. C., Pross, J., Sangiorgi, F., Middelburg, J. B. M., and Sinninghe Damsé, J. S.: Influence of deep-water derived isoprenoid tetraether lipids on the TEX86H paleothermometer in the Mediterranean Sea, Geochim. Cosmochim. Ac., 150, 125–141,, 2015. 

Kim, J.-H., Villanueva, L., Zell, C., and Sinninghe Damsté, J. S.: Biological source and provenance of deep-water derived isoprenoid tetraether lipids along the Portuguese continental margin, Geochim. Cosmochim. Ac., 172, 177–204,, 2016. 

Laskar, J.: The chaotic motion of the solar system: A numerical estimate of the size of the chaotic zones, Icarus, 88, 266–291,, 1990. 

Lengger, S. K., Hopmans, E. C., Sinninghe Damsté, J. S., and Schouten, S.: Impact of sedimentary degradation and deep water column production on GDGT abundance and distribution in surface sediments in the Arabian Sea: Implications for the TEX86 paleothermometer, Geochim. Cosmochim. Ac., 142, 386–399,, 2014. 

Lisiecki, L. E. and Raymo, M. E.: A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18O records, Paleoceanography, 20, PA1003,, 2005. 

Locarnini, R. A., Mishonov, A. V., Antonov, J. I., Boyer, T. P., Garcia, H. E., Baranova, O. K., Zweng, M. M., Paver, C. R., Reagan, J. R., Johnson, D. R., Hamilton, M., and Seidov, D.: World Ocean Atlas 2013. Vol. 1: Temperature, edited by: Levitus, S. and Mishonov, A., NOAA Atlas NESDIS, 73, 40 pp.,, 2013. 

Lopes dos Santos, R. A., Prange, M., Castañeda, I. S., Schefuß, E., Mulitza, S., Schulz, M., Niedermeyer, E. M., Sinninghe Damsté, J. S., and Schouten, S.: Glacial-interglacial variability in Atlantic meridional overturning circulation and thermocline adjustments in the tropical North Atlantic, Earth Planet. Sc. Lett., 300, 407–414,, 2010. 

Martin, J. H., Knauer, G. A., Karl, D. M., and Broenkow, W. W.: VERTEX: carbon cycling in the northeast Pacific, Deep-Sea Res., 34, 267–285,, 1987. 

Martínez-Garcia, A., Rosell-Melé, A., Jaccard, S. L., Geibert, W., Sigman, D. M., and Haug, G. H.: Southern Ocean dust–climate coupling over the past four million years, Nature, 476, 312–315,, 2011. 

Mascle, J., Lohmann, G., Clift, P. D., and and Shipboard Scientific Party: Proceedings of the Ocean Drilling Program, Initial Reports, Vol. 159, College Station, TX, Ocean Drilling Program,, 1996. 

McClymont, E. L., Ganeshram, R. S., Pichevin, L. E., Talbot, H. M., Van Dongen, B. E., Thunell, R. C., Haywood, A. M., Singarayer, J. S., and Valdes, P. J.: Sea-surface temperature records of Termination 1 in the Gulf of California: Challenges for seasonal and interannual analogues of tropical Pacific climate change, Paleoceanography, 27, PA2202,, 2012. 

McKay, R., Naish, T., Carter, L., Riesselman, C., Dunbar, R., Sjunneskog, C., Winter, D., Sangiorgi, F., Warren, C., Pagani, M., Schouten, S., Willmott, V., Levy, R., DeConto, R., and Powell, R. D.: Antarctic and Southern Ocean influences on Late Pliocene global cooling, P. Natl. Acad. Sci. USA, 109, 6423–6428,, 2012. 

Middelburg, J. J.: Marine Carbon Biogeochemistry, Briefs in Earth System Sciences, Springer, Cham, 1–8,, 2019. 

Müller, P. J., Kirst, G., Ruhland, G., Von Storch, I., and Rosell-Melé, A.: Calibration of the alkenone paleotemperature index U37K based on core-tops from the eastern South Atlantic and the global ocean (60 N–60 S), Geochim. Cosmochim. Ac., 62, 1757–1772,, 1998. 

Norris, R. D.: 40. Miocene-Pliocene surface-water hydrography of the eastern equatorial Atlantic, in: Proceedings of the Ocean Drilling Program, Scientific Results, edited by: Mascle, J., Lohmann, G. P., and Moullade, M., College Station, TX, Ocean Drilling Program, 159, 539–555,, 1998a. 

Norris, R. D.: 39. Planktonic foraminifer biostratigraphy: Eastern equatorial Atlantic in: Proceedings of the Ocean Drilling Program, Scientific Results, edited by: Mascle, J., Lohmann, G. P., and Moullade, M., College Station, TX, Ocean Drilling Program, 159, 445-479,, 1998b. 

Park, E., Hefter, J., Fischer, G., and Mollenhauer, G.: TEX86 in sinking particles in three eastern Atlantic upwelling regimes, Org. Geochem., 124, 151–163,, 2018. 

Patterson, M. O., McKay, R., Naish, T. R., Escutia, C., Jimenez-Espejo, F. J., Raymo, M. E., Meyers, S. R., Tauxe, L., Brinkhuis, H., and IODP Expedition 318 Scientists: Orbital forcing of the East Antarctic ice sheet during the Pliocene and Early Pleistocene, Nat. Geosci, 7, 841–847,, 2014. 

Pearson, A., Hurley, S. J., Walter, S. R. S., Kusch, S., Lichtin, S., and Zhang, Y. G.: Stable carbon isotope ratios of intact GDGTs indicate heterogeneous sources to marine sediments, Geochim. Cosmochim. Ac., 181, 18–35,, 2016. 

Petit, J. R., Jouzel, J., Raynaud, D., Barkov, N. I., Barnola, J.-M., Basile, I., Bender, M., Chappellaz, J., Davisk, M., Delaygue, G., Delmotte, M., Kotlyakov, V. M., Legrand, M., Lipenkov, V. Y., Lorius, C., Pé, L., Ritz, C., Saltzmank, E., and Stievenard, M.: Climate and atmospheric history of the past 420,000 years from the Vostok ice core, Antarctica The recent completion of drilling at Vostok station in East, Nature, 399, 429–436, 1999. 

Prahl, F. G. and Wakeham, S. G.: Calibration of unsaturation patterns in long-chain ketone compositions for palaeotemperature assessment, Nature, 330, 367–369,, 1987. 

Rommerskirchen, F., Condon, T., Mollenhauer, G., Dupont, L., and Schefuss, E.: Miocene to Pliocene development of surface and subsurface temperatures in the Benguela Current system, Paleoceanography, 26, PA3216,, 2011. 

Schouten, S., Hopmans, E. C., Schefuß, E., and Sinninghe Damsté, J. S.: Distributional variations in marine crenarchaeotal membrane lipids: a new tool for reconstructing ancient sea water temperatures?, Earth Planet. Sc. Lett., 204, 265–274,, 2002. 

Schouten, S., Pitcher, A., Hopmans, E. C., Villanueva, L., van Bleijswijk, J., and Sinninghe Damsté, J. S.: Intact polar and core glycerol dibiphytanyl glycerol tetraether lipids in the Arabian Sea oxygen minimum zone: I. Selective preservation and degradation in the water column and consequences for the TEX86, Geochim. Cosmochim. Ac., 98, 228–243,, 2012. 

Shafik, S., Watkins, D. K. and Shin, I. C.: 37. Upper Cenozoic calcareous nannofossil biostratigraphy, Cote d'Ivoire-Ghana Margin, eastern equatorial Atlantic, Proceedings of the Ocean Drilling Program, Scientific Results, 159, 509–523,, 1998. 

Shah, S. R., Mollenhauer, G., Ohkouchi, N., Eglinton, T. I., and Pearson, A.: Origins of archaeal tetraether lipids in sediments: Insights from radiocarbon analysis, Geochim. Cosmochim. Ac., 72, 4577–4594,, 2008. 

Shin, I. C., Shafik, S., and Watkins, D. K.: High-resolution Pliocene-Pleistocene biostratigraphy of Site 959, eastern equatorial Atlantic Ocean, Proceedings of the Ocean Drilling Program, Scientific Results, 159, 533–538, 1998. 

Sintes, E., De Corte, D., Haberleitner, E., and Herndl, G. J.: Geographic Distribution of Archaeal Ammonia Oxidizing Ecotypes in the Atlantic Ocean, Front. Microbiol., 7, 77,, 2016. 

Sonzogni, C., Bard, E., Rostek, F., Lafont, R., Rosell-Mele, A., and Eglinton, G.: Core-top calibration of the alkenone index vs sea surface temperature in the Indian Ocean, Deep-Sea Res. Pt. II, 44, 1445–1460,, 1997. 

Stockmarr, J.: Tablets with spores used in absolute pollen analysis, Pollen et Spores, 13, 615–621, 1972. 

Taylor, K. W. R., Huber, M., Hollis, C. J., Hernandez-Sanchez, M. T., and Pancost, R. D.: Re-evaluating modern and Palaeogene GDGT distributions: Implications for SST reconstructions, Global Planet. Change, 108, 158–174,, 2013. 

Tierney, J. E. and Tingley, M. P.: A Bayesian, spatially-varying calibration model for the TEX86 proxy, Geochim. Cosmochim. Ac., 127, 83–106,, 2014. 

Tierney, J. E. and Tingley, M. P.: A TEX86 surface sediment database and extended Bayesian calibration, Sci. Data, 2, 150029,, 2015. 

Tierney, J. E. and Tingley, M. P.: BAYSPLINE: A new calibration for the alkenone paleothermometer, Palaeogeogr. Palaeocl., 33, 281–301,, 2018. Tierney, J. E., Sinninghe Damsté, J. S., Pancost, R. D., Sluijs, A., and Zachos, J. C.: Eocene temperature gradients, Nat. Geosci., 10, 538–539,, 2017. 

Tierney, J. E., Poulsen, C. J., Montañez, I. P., Bhattacharya, T., Feng, R., Ford, H. L., Hönisch, B., Inglis, G. N., Petersen, S. V., Sagoo, N., Tabor, C. R., Thirumalai, K., Zhu, J., Burls, N. J., Foster, G. L., Goddéris, Y., Huber, B. T., Ivany, L. C., Turner, S. K., Lunt, D. J., McElwain, J. C., Mills, B. J. W., Otto-Bliesner, B. L., Ridgwell, A., and Zhang, Y. G.: Past climates inform our future, Science, 370, eaay3701,, 2020. 

Vallé, F., Westerhold, T., and Dupont, L.: Orbital-driven environmental changes recorded at ODP Site 959 (eastern equatorial Atlantic) from the Late Miocene to the Early Pleistocene, Int. J. Earth Sci., 106, 1161–1174,, 2016. 

van der Weijst, C. M. H., Winkelhorst, J., Lourens, L., Raymo, M. E., Sangiorgi, F., and Sluijs, A.: A Ternary Mixing Model Approach Using Benthic Foraminifer δ13C-δ18O Data to Reconstruct Late Pliocene Deep Atlantic Water Mass Mixing, Paleoceanogr. Paleocl., 35, e2019PA003804,, 2020. 

van der Weijst, C. M. H., Winkelhorst, J., de Nooijer, W., von der Heydt, A., Reichart, G.-J., Sangiorgi, F., and Sluijs, A.: Pliocene evolution of the tropical Atlantic thermocline depth, Clim. Past, 18, 961–973,, 2022. 

Versteegh, G. J. M.: Recognition of cyclic and non-cyclic environmental changes in the Mediterranean Pliocene: A palynological approach, Mar. Micropaleontol., 23, 147–183,, 1994. 

Verstraete, J. M.: The seasonal upwellings in the Gulf of Guinea, Prog. Oceanogr., 29, 1–60,, 1992. 

Villanueva, L., Schouten, S., and Sinninghe Damsté, J. S.: Depth-related distribution of a key gene of the tetraether lipid biosynthetic pathway in marine Thaumarchaeota, Environ. Microbiol., 17, 3527–3539,, 2015. 

Wagner, T.: 41. Pliocene-Pleistocene deposition of carbonate and organic carbon at Site 959: Paleoenvironmental implications for the eastern equatorial Atlantic off the Ivory Coast/Ghana, in: Proceedings of the Ocean Drilling Program, Scientific Results, edited by: Mascle, J., Lohmann, G. P., and Moullade, M., College Station, TX (Ocean Drilling Program), 159, 557–574,, 1998. 

Weijers, J. W. H., Schouten, S., Spaargaren, O. C., and Sinninghe Damsté, J. S.: Occurrence and distribution of tetraether membrane lipids in soils: Implications for the use of the TEX86 proxy and the BIT index, Org. Geochem., 37, 1680–1693,, 2006. 

White, S. M. and Ravelo, A. C.: The benthic B/Ca record at Site 806: new constraints on the temperature of the West Pacific Warm Pool and the “El Padre” state in the Pliocene, Paleoceanogr. Paleocl., 35, e2019PA003812,, 2020. 

Wiafe, G. and Nyadjro, E. S.: Satellite observations of upwelling in the gulf of Guinea, IEEE Geosci. Remote S., 12, 1066–1070,, 2015. 

Williams, G. L., Fensome, R. A., and MacRae, R. A.: The Lentin and Williams Index of Fossil Dinoflagellates, 2017 Edition, AASP Contributions Series Number 48, American Association of Stratigraphic Palynologists Foundation, 2017.  

Wuchter, C., Schouten, S., Wakeham, S. G., and Sinninghe Damsté, J. S.: Temporal and spatial variation in tetraether membrane lipids of marine Crenarchaeota in particulate organic matter: Implications for TEX86 paleothermometry, Paleoceanography, 20, PA3013,, 2005. 

Zhang, R., Zhang, Z., Jiang, D., Yan, Q., Zhou, X., and Cheng, Z.: Strengthened African summer monsoon in the mid-Piacenzian, Adv. Atmos. Sci., 33, 1061–1070,, 2016. 

Zhang, Y. G. and Liu, X.: Export Depth of the TEX 86 Signal, Paleoceanogr. Paleocl., 666–671,, 2018. 

Zhang, Y. G., Zhang, C. L., Liu, X. L., Li, L., Hinrichs, K. U., and Noakes, J. E.: Methane Index: A tetraether archaeal lipid biomarker indicator for detecting the instability of marine gas hydrates, Earth Planet. Sc. Lett., 307, 525–534,, 2011. 

Zhang, Y. G., Pagani, M., and Wang, Z.: Ring Index: A new strategy to evaluate the integrity of TEX86 paleothermometry, Paleoceanography, 31, 220–232,, 2016. 

Zonneveld, K. A. F., Marret, F., Versteegh, G., Bogus, K., Bonnet, S., Bouimetarhan, I., Crouch, E., de Vernal, A., Elshanawany, R., Edwards, L., Esper, O., Forke, S., Grøsfjeld, K., Henry, M., Holzwarth, U., Kielt, J. F., Kim, S. Y., Ladouceur, S., Ledu, D., Chen, L., Limoges, A., Londeix, L., Lu, S. H., Mahmoud, M. S., Marino, G., Matsouka, K., Matthiessen, J., Mildenhal, D. C., Mudie, P., Neil, H. L., Pospelova, V., Qi, Y., Radi, T., Richerol, T., Rochon, A., Sangiorgi, F., Solignac, S., Turon, J. L., Verleye, T., Wang, Y., Wang, Z., and Young, M.: Atlas of modern dinoflagellate cyst distribution based on 2405 data points, Rev. Palaeobot. Palyno., 191, 1–197,, 2013. 

Zweng, M. M., Reagan, J. R., Antonov, J. I., Locarnini, R. A., Mishonov, A. V., Boyer, T. P., Garcia, H. E., Baranova, O. K., Johnson, D. R., Seidov, D., and Biddle, M. M.: World Ocean Atlas 2013, Volume 2: Salinity, NOAA Atlas NESDIS, 74, 1–39,, 2013. 

Short summary
The TEX86 proxy is often used by paleoceanographers to reconstruct past sea-surface temperatures. However, the origin of the TEX86 signal in marine sediments has been debated since the proxy was first proposed. In our paper, we show that TEX86 carries a mixed sea-surface and subsurface temperature signal and should be calibrated accordingly. Using our 15-million-year record, we subsequently show how a TEX86 subsurface temperature record can be used to inform us on past sea-surface temperatures.