A 15-million-year surface- and subsurface-integrated TEX86 temperature record from the eastern equatorial Atlantic
- 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 (email@example.com)
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 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 (U and ) 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 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 ∼5 ∘C 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.
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 (foraminifer calcite), U (alkenones) and TEX86 (glycerol dialkyl glycerol tetraethers; GDGTs). Each proxy has specific confounding factors and calibration issues. For example, paleothermometry requires corrections for the ratio of seawater and diagenesis (Dekens et al., 2002; Evans et al., 2016), and the U proxy is insensitive to temperatures above ∼28 ∘C (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 ratio; Hernández-Sánchez et al., 2014; Pearson et al., 2016; Taylor et al., 2013; Villanueva et al., 2015). GDGT 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 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).
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 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 evolution on million-year and glacial–interglacial timescales during the late Pliocene. We complement these records with U 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 -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.
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. U 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, ) 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, ), hexane : DCM (1:1, ) and DCM : methanol (1:1, ) 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 (). 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. U 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).
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 (2×) and 40 % cold HF (2×) 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 ratio (Versteegh, 1994) as a paleoproductivity index. The 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.1 TEX86, U and GDGT evolution at Site 959
We explore our TEX86 results using several calibration approaches in Fig. 2, including an SST calibration using the exponential TEX (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 TEX calibration for SST (Figs. 3–7), while subsurface calibrations are discussed in Figs. 7 and 8.
Between 15 and 11 Ma, TEX-based SST values fluctuate around an average value of ∼29 ∘C (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 ∼2 ∘C each, around 4.9 and 3.4 Ma, in which TEX 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 TEX records a ∼4 ∘C warming between 0.5 Ma and present, with a core top TEX of 24.0 ∘C. U is at saturation (∼28 ∘C) 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 ratio is 7.0 in the core top. In the oldest part of the studied interval, GDGT values vary between 4 and 5.5. Following an abrupt increase at ∼13.6 Ma, GDGT varies largely between 5 and 8. After 1.7 Ma, GDGT 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 TEX-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 paleoproductivity index is 0.48 on average, with minimum and maximum values of 0.23 and 0.70, respectively (Fig. 4). values are primarily driven by the abundance of Brigantedinium spp. The ratio shows no significant relationship with SST as derived from TEX (R2=0.05, p=0.14 in linear regression analysis).
4.1 Identifying the source depth of the TEX86 signal
The 15-million-year TEX record at Site 959 shows a cooling from >30 ∘C between 11 and 15 Ma to <20 ∘C between 0 and 3.3 Ma (Fig. 3). The U record, on the other hand, is near or at proxy saturation (>27 ∘C) for most of the interval, with some lower values after 1.8 Ma. Assuming that both proxies reflect SST, the relatively small offset between U and TEX-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 TEX and U 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 U 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 TEX (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.
4.1.1 GDGT values
A substantial contribution of GDGTs from below the oceanic surface layer is supported by the relatively high GDGT 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 TEX-SST and GDGT (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 increase (Fig. 3), signaling systematic deepening of the source of the sedimentary GDGT signal.
The ratio of GDGT to GDGT 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 increases relative to that of GDGT with export depth, then samples with a higher GDGT ratio should be associated with lower TEX86 values. The positive relationship between TEX86 and GDGT (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 ratios in sediments, recently supported by SPM analyses (Hurley et al., 2018). They also conclude that the positive correlation between SST and GDGT, 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 shift at 13.6 Ma, TEX does not systematically change until ∼11 Ma (Fig. 3). Moreover, late Pliocene dinocyst assemblages (Fig. 4) are characterized by a highly variable ratio, indicating that upwelling and nutrient supply were highly variable on sub-Milankovitch timescales. In contrast to TEX (Fig. 5), the ratio shows no glacial–interglacial variability. It is therefore unlikely that TEX86 variability was primarily driven by non-thermal factors.
4.1.2 Surface, thermocline and sub-thermocline temperature variability
The late Pliocene TEX record follows a different evolution than SST records of U and the recent of the surface-dwelling Globigerinoides ruber ((G. ruber)) and thermocline-dwelling Neogloboquadrina dutertrei ((N. dutertrei)) generated on the same samples (Fig. 5). The records (van der Weijst et al., 2022) are based on very well-preserved foraminifera and apply a constant Pliocene correction for the of seawater, a species-specific core top calibration with a constant partitioning coefficient and embedded dissolution correction (Dekens et al., 2002). Whereas TEX decreases by ∼2 ∘C, U and -based SSTs show no significant net change in this interval. Moreover, although and U 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, TEX cooling leads to increasing abundances of cool dinocysts, as well as (G. ruber) and U 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, TEX 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 TEX 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 TEX record shows weak coherence with the Site 959 (N. dutertrei) thermocline temperature record (Fig. 5; van der Weijst et al., 2022). Whereas TEX86 registers late Pliocene cooling, (N. dutertrei) registers warming in relation to thermocline deepening (van der Weijst et al., 2022). Moreover, while the range of (N. dutertrei) temperatures is similar to TEX on glacial–interglacial timescales, (N. dutertrei) does not follow the glacial–interglacial cyclicity as recorded in the TEX86.
In contrast, the variability in the TEX-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 TEX leads δ18O by ∼5 kyr (Fig. 6), which indicates that, despite the similarities between TEX-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 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.
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 TEX-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 TEX-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 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).
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 (TEX-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 U 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 ∼5 ∘C 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 TEX-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.
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 TEX-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.
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 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 (U and (G. ruber) SST and (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 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.
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: https://doi.org/10.5194/cp-18-1947-2022-supplement.
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.
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.
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.
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, https://doi.org/10.1007/s00382-005-0086-4, 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, https://doi.org/10.5194/cp-15-1603-2019, 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, https://doi.org/10.1016/j.orggeochem.2019.06.008, 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, https://doi.org/10.1002/jqs.654, 2001.
Dekens, P. S., Lea, D. W., Pak, D. K., and Spero, H. J.: Core top calibration of in tropical foraminifera: Refining paleotemperature estimation, Geochem. Geophy. Geosy., 3, 1–29, https://doi.org/10.1029/2001GC000200, 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, https://doi.org/10.1038/s41598-020-67154-8, 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, https://doi.org/10.1029/2008PA001725, 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, https://doi.org/10.1371/journal.pone.0081508, 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, https://doi.org/10.1016/j.epsl.2018.07.035, 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, https://doi.org/10.1175/JPO-D-16-0126.1, 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, https://doi.org/10.1016/j.gloplacha.2015.02.001, 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 , temperature and sea level change, Earth Planet. Sc. Lett., 438, 139–148, https://doi.org/10.1016/j.epsl.2016.01.013, 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, https://doi.org/10.1098/rsta.2012.0294, 2013.
Haug, G. H., Tiedemann, R., Zahn, R., and Ravelo, A. C.: Role of Panama uplift on oceanic freshwater balance, Geology, 29, 207–210, https://doi.org/10.1130/0091-7613(2001)029<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, https://doi.org/10.5194/cp-16-2095-2020, 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, https://doi.org/10.1017/S0022336000027153, 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, https://doi.org/10.1126/science.1185435, 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, https://doi.org/10.1016/j.gca.2014.02.009, 2014.
Ho, S. L. and Laepple, T.: Glacial cooling as inferred from marine temperature proxies TEX and U, Earth Planet. Sci. Lett., 409, 15–22, https://doi.org/10.1016/j.epsl.2014.10.033, 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, https://doi.org/10.1038/ngeo2763, 2016.
Ho, S. L. and Laepple, T.: Reply to “Eocene temperature gradients,” Nat. Geosci., 10, 539–540, https://doi.org/10.1038/ngeo2998, 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, https://doi.org/10.1016/j.epsl.2004.05.012, 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, https://doi.org/10.1016/j.orggeochem.2015.12.006, 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 U), Paleoceanography, 21, PA3003, https://doi.org/10.1029/2005PA001215, 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, https://doi.org/10.1073/pnas.1518534113, 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, https://doi.org/10.1016/j.orggeochem.2017.11.010, 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, https://doi.org/10.1073/pnas.0510157103, 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, https://doi.org/10.1016/j.marmicro.2020.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, https://doi.org/10.1016/j.gca.2007.12.010, 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, https://doi.org/10.1016/j.gca.2010.05.027, 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, https://doi.org/10.1016/j.epsl.2012.05.018, 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 TEX paleothermometer in the Mediterranean Sea, Geochim. Cosmochim. Ac., 150, 125–141, https://doi.org/10.1016/j.gca.2014.11.017, 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, https://doi.org/10.1016/j.gca.2015.09.010, 2016.
Laskar, J.: The chaotic motion of the solar system: A numerical estimate of the size of the chaotic zones, Icarus, 88, 266–291, https://doi.org/10.1016/0019-1035(90)90084-M, 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, https://doi.org/10.1016/j.gca.2014.07.013, 2014.
Lisiecki, L. E. and Raymo, M. E.: A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18O records, Paleoceanography, 20, PA1003, https://doi.org/10.1029/2004PA001071, 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., https://doi.org/10.7289/V55X26VD, 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, https://doi.org/10.1016/j.epsl.2010.10.030, 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, https://doi.org/10.1016/0198-0149(87)90086-0, 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, https://doi.org/10.1038/nature10310, 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, https://doi.org/10.2973/odp.proc.ir.159.1996, 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, https://doi.org/10.1029/2011PA002226, 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, https://doi.org/10.1073/pnas.1112248109, 2012.
Middelburg, J. J.: Marine Carbon Biogeochemistry, Briefs in Earth System Sciences, Springer, Cham, 1–8, https://doi.org/10.1007/978-3-030-10822-9_1, 2019.
Müller, P. J., Kirst, G., Ruhland, G., Von Storch, I., and Rosell-Melé, A.: Calibration of the alkenone paleotemperature index U based on core-tops from the eastern South Atlantic and the global ocean (60∘ N–60∘ S), Geochim. Cosmochim. Ac., 62, 1757–1772, https://doi.org/10.1016/S0016-7037(98)00097-0, 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, https://doi.org/10.2973/odp.proc.sr.159.021.1998, 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, https://doi.org/10.2973/odp.proc.sr.159.036.1998, 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, https://doi.org/10.1016/j.orggeochem.2018.07.015, 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, https://doi.org/10.1038/ngeo2273, 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, https://doi.org/10.1016/j.gca.2016.02.034, 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, https://doi.org/10.1038/330367a0, 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, https://doi.org/10.1029/2010PA002074, 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, https://doi.org/10.1016/S0012-821X(02)00979-2, 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, https://doi.org/10.1016/j.gca.2012.05.002, 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, https://doi.org/10.2973/odp.proc.sr.159.022.1998, 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, https://doi.org/10.1016/j.gca.2008.06.021, 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, https://doi.org/10.3389/fmicb.2016.00077, 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, https://doi.org/10.1016/S0967-0645(97)00010-6, 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, https://doi.org/10.1016/j.gloplacha.2013.06.011, 2013.
Tierney, J. E. and Tingley, M. P.: A Bayesian, spatially-varying calibration model for the TEX86 proxy, Geochim. Cosmochim. Ac., 127, 83–106, https://doi.org/10.1016/j.gca.2013.11.026, 2014.
Tierney, J. E. and Tingley, M. P.: A TEX86 surface sediment database and extended Bayesian calibration, Sci. Data, 2, 150029, https://doi.org/10.1038/sdata.2015.29, 2015.
Tierney, J. E. and Tingley, M. P.: BAYSPLINE: A new calibration for the alkenone paleothermometer, Palaeogeogr. Palaeocl., 33, 281–301, https://doi.org/10.1002/2017PA003201, 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, https://doi.org/10.1038/ngeo2997, 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, https://doi.org/10.1126/science.aay3701, 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, https://doi.org/10.1007/s00531-016-1350-z, 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, https://doi.org/10.1029/2019PA003804, 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, https://doi.org/10.5194/cp-18-961-2022, 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, https://doi.org/10.1016/0377-8398(94)90005-1, 1994.
Verstraete, J. M.: The seasonal upwellings in the Gulf of Guinea, Prog. Oceanogr., 29, 1–60, https://doi.org/10.1016/0079-6611(92)90002-H, 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, https://doi.org/10.1111/1462-2920.12508, 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, https://doi.org/10.2973/odp.proc.sr.159.018.1998, 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, https://doi.org/10.1016/j.orggeochem.2006.07.018, 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, https://doi.org/10.1029/2019pa003812, 2020.
Wiafe, G. and Nyadjro, E. S.: Satellite observations of upwelling in the gulf of Guinea, IEEE Geosci. Remote S., 12, 1066–1070, https://doi.org/10.1109/LGRS.2014.2379474, 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, https://doi.org/10.1029/2004PA001110, 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, https://doi.org/10.1007/s00376-016-5215-y, 2016.
Zhang, Y. G. and Liu, X.: Export Depth of the TEX 86 Signal, Paleoceanogr. Paleocl., 666–671, https://doi.org/10.1029/2018PA003337, 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, https://doi.org/10.1016/j.epsl.2011.05.031, 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, https://doi.org/10.1002/2015PA002848, 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, https://doi.org/10.1016/j.revpalbo.2012.08.003, 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, https://doi.org/10.7289/V5251G4D, 2013.