The Paleocene–Eocene Thermal Maximum (PETM), ca. 56 Ma, was a major global environmental perturbation attributed to a rapid rise in the concentration of greenhouse gases in the atmosphere. Geochemical records of tropical sea-surface temperatures (SSTs) from the PETM are rare and are typically affected by post-depositional diagenesis. To circumvent this issue, we have analyzed oxygen isotope ratios (δ18O) of single specimens of exceptionally well-preserved planktonic foraminifera from the PETM in Tanzania (∼19°S paleolatitude), which yield extremely low δ18O, down to <–5‰. After accounting for changes in seawater chemistry and pH, we estimate from the foraminifer δ18O that tropical SSTs rose by >3 °C during the PETM and may have exceeded 40 °C. Calcareous plankton are absent from a large part of the Tanzania PETM record; extreme environmental change may have temporarily caused foraminiferal exclusion.
During the Paleocene–Eocene Thermal Maximum (PETM), >2000 Gt of isotopically light carbon was released into the atmosphere in <60 k.y., possibly by the destabilization of deep-sea methane hydrates (Dickens, 2011) or soil organic carbon within permafrost (DeConto et al., 2012). The carbon release caused a substantial negative carbon isotope (δ13C) excursion (CIE), the magnitude of which remains uncertain (e.g., Zachos et al., 2007). The PETM was associated with globally averaged warming estimated to be >5 °C (Dunkley Jones et al., 2013). The absorption of such large quantities of carbon into the ocean resulted in a lowering of oceanic pH and a shoaling of the calcium carbonate compensation depth (Zachos et al., 2005). Records of tropical sea-surface temperatures (SSTs) from calcareous organisms are rare, as deep-sea sediments commonly have poor microfossil preservation (Zachos et al., 2007). Here we report new geochemical and faunal data from an expanded section from the continental margin of East Africa that provides information about the PETM in the tropics.
Tanzania Drilling Project Site 14 (TDP-14) (9°16′59.89″S, 39°30′45.04″E) (Nicholas et al., 2006) comprises two ∼30 m holes drilled 1 m apart in late Paleocene–early Eocene hemipelagic sediment. The site was at ∼19°S paleolatitude in an outer shelf to slope bathyal environment, estimated at water depths of >300 m (Nicholas et al., 2006). The principal lithologies are claystone and siltstone, with excellent microfossil preservation (e.g., Bown and Pearson, 2009). The site’s proximity to the continent (∼70 km to the paleoshoreline; Kent et al., 1971) and shallow burial history mean that TDP-14 contains abundant well-preserved organic biomarkers (van Dongen et al., 2006).
Planktonic foraminifer assemblage count samples were sieved at 63 μm, dried, and split. The first 300 (if present) specimens were counted from the 125–600 μm fraction. Calcareous nannofossils were viewed as smear slides (Bown and Young, 1998), using microscopy in cross-polarized and phase-contrast light on rock surfaces using scanning electron microscopy (Lees et al., 2004). Count data were from more than five fields of view until ∼400 specimens were counted.
Species of the genera Acarinina, Morozovella, and Subbotina from the 300–355 μm fraction were screened for preservation quality and infilling. Individuals with excellent preservation (see the GSA Data Repository1) were analyzed as single specimens for carbon and oxygen isotope composition (δ13Cforam and δ18Oforam), augmenting multiple-specimen Subbotina data published by Handley et al. (2008). Specimens were analyzed using an IsoPrime mass spectrometer, and data are reported to the Vienna Peedee belemnite scale.
For δ13C analyses of C25–C31n-alkanes (δ13Calk), 40–75 g of sediment was solvent extracted and separated into an alkane fraction by alumina flash column chromatography and urea adduction. Analyses were performed on a ThermoFinnigan DeltaplusXP coupled to a Trace 2000 GC via a GC-C III interface.
The abundances of organic carbon and calcium carbonate were measured from aliquots of sediment powders using a LECO CNS-2000 elemental analyzer. Decarbonated sample powders were also analyzed for their organic carbon isotope composition (δ13Corg) using a Thermo Flash HT elemental analyzer coupled to a ThermoFinnigan MAT 253 mass spectrometer (see the Data Repository for reproducibility and uncertainties of our methods).
Sedimentological and Biotic Records
Deposits of the PETM in TDP-14 are ∼11 m thick. The base of the PETM interval is ∼24 m below surface (mbs), as defined by the first occurrence of the excursion taxon Acarinina africana and a negative shift in δ13Calk, along with some lower δ13C and δ18O in single specimens of planktonic foraminifera (Fig. 1). The lower part of the PETM interval is sedimentologically complex, with fluctuating abundances of organic carbon and δ13Corg that likely reflect changes in the relative abundance of distinct sources of organic matter (Fig. 1). No mass transport features were observed in the core (Nicholas et al., 2006), but reworking, switching sources, or sediment mixing may account for this complexity. The pre-excursion sediments and reworked packets contain an abundant open-ocean planktonic foraminifer assemblage comprising more than 40 species (see the Data Repository). Stratigraphic horizons containing excursion taxa commonly have low abundance (commonly <1 foraminifer/g) and diversity of planktonic foraminifera. Between ∼24 and 18 mbs, planktonic foraminifer abundances fluctuate between pre-excursion levels and very low abundances, with a near complete absence of calcareous microfossils between ∼18 and 13 mbs. Where planktonic foraminifer abundances are lowest, the CaCO3 content drops to near zero (Fig. 1). Between 24 and 19 mbs (“mixed zone” in Fig. 1), levels with higher microfossil abundances and δ13Cforam are interpreted as dominated by transported pre-PETM sediments. The top of the PETM section is truncated by a hiatus (∼13 mbs), and a rich and diverse microfossil assemblage returns higher in the core (13–10 mbs).
Single-specimen pre-PETM δ13Cforam values for mixed layer– and thermocline-dwelling species are typically ∼4.9‰ (± 1.94‰, 2 standard deviations [SD]) and 1.5‰ (± 1.11‰, 2 SD), respectively (Figs. 1 and 2). Mean pre-PETM δ18Oforam values are –3.3‰ for mixed layer– and –2.7‰ for thermocline-dwelling species (Figs. 1 and 2). Two Morozovella specimens from within the CIE exhibit δ18Oforam lower than –5‰. However, not all specimens that exhibit the lowest δ13Cforam also record the lowest δ18Oforam (see the Data Repository). Due to the complex stratigraphy at TDP-14, PETM and pre-PETM specimens occur in the same stratigraphic intervals, which makes an estimation of the true magnitude of the CIE from foraminiferal calcite problematic.
To quantify paleo-SSTs, we explore variables that can impact estimates including (1) seawater δ18O (δ18Osw), (2) pH, and (3) the paleotemperature equation.
Large continental ice sheets are unlikely to have been present during the Paleocene-Eocene transition (e.g., Sluijs et al., 2008), so all SSTs are estimated using a global δ18Osw correction of –1.0‰ (Cramer et al., 2011). Zachos et al. (1994) produced a correction factor for modern latitudinal variations in δ18Osw as a function of the transport of water vapor from the low to high latitudes, giving SSTs that are up to ∼4 °C higher than uncorrected values (because the latitude of TDP-14 is one of net evaporation). Local hydrological cycle changes may have influenced the input of meteoric waters to the surface ocean near the continental margin during the PETM, affecting the δ18Osw. Deuterium enrichment of n-alkanes from TDP-14 implies that regional conditions during the PETM were hotter and more arid than today, punctuated by intense seasonal precipitation events (Handley et al., 2012). The effect of adding meteoric waters to the surface ocean directly or via rivers is small (–0.06‰ δ18O change for a 1‰ decrease in salinity) (Damassa et al., 2006), therefore a major drop in salinity sufficient to lower the local δ18Osw seems unlikely; if anything, more arid conditions during the PETM in this region would make the δ18Oforam SSTs slight underestimates. For these reasons we have not corrected for local δ18O variability.
The changing carbonate ion concentration and decrease in pH associated with elevated CO2 levels during the PETM may have resulted in higher δ13Cforam and δ18Oforam because of the “pH effect” (e.g., Spero et al., 1997). We have corrected for changes in the isotopic fractionation of oxygen due to potential pH shifts (Uchikawa and Zeebe, 2010) using the end-member values of pH decline for the PETM (–0.25 to –0.45) (cf. Hönisch et al., 2012) to inform our pH correction, resulting in SSTs that are higher by up to ∼1.5 °C than uncorrected values.
The paleotemperature equation of Kim and O’Neil (1997) is favored because it is the only one that is calibrated above 30 °C. For comparison, we have also applied the foraminifera equations of Erez and Luz (1983) and Bemis et al. (1998) (Table 1).
To explore the effects of the paleotemperature assumptions and equations outlined above, we calculated SSTs from single-specimen δ18Oforam of the mixed-layer planktonic foraminifer genus Morozovella using the mean pre-PETM, mean PETM, and lowest recorded PETM δ18Oforam (Table 1). Resulting PETM SSTs from TDP-14 range from ∼28 °C to 43 °C. Model simulations suggest that PETM tropical SSTs could have reached ∼35 °C (Huber and Caballero, 2011), but following corrections for both paleolatitude and ocean acidification, PETM SSTs at this site most likely reached between ∼36 °C and 43 °C (Table 1; Fig. 2; see the Data Repository). The calculations suggest that prior to the CIE, SSTs were warmer than modern SSTs at equivalent latitudes, and that during the PETM, SSTs were significantly higher. The change in temperature using the mean background and mean CIE δ18Oforam is ∼+3 °C (± 0.5 °C). The lowest δ18Oforam values likely reflect maxima in seasonal warmth on a background of already high temperatures; this may explain the absence of significant temperature change recorded by thermocline species that were living at greater depths.
Plankton Assemblage Changes
The absence of calcareous plankton from much of the main PETM interval at TDP-14 may be due to dissolution, exclusion due to environmental stress (e.g., Kelly et al., 1996), or sedimentary dilution. Dissolution effects are unlikely as there is no increase in planktonic foraminifer fragmentation (Hemleben et al., 1989) prior to their disappearance, and where present, diminutive and fragile nannofossils are exceptionally well-preserved (Bown and Pearson, 2009).
Shelf PETM sections show increased sedimentation rates related to hydrological cycle changes (e.g., Sluijs et al., 2008), and this likely also happened in Tanzania (Handley et al. 2012) leading to decreased calcareous plankton concentrations. The complex stratigraphy and short core length hinder development of an age model that could determine whether the decline in calcareous plankton abundance is the result of sediment dilution. However, the reduction in planktonic foraminifer specimens per gram (s/g) before (average = 86 s/g) and during the peak of the CIE (average = 1 s/g) (Fig. 1) would require sedimentation rates to increase by a factor of ∼80, whereas marine biomarker concentrations indicate a sedimentation increase by an order of magnitude less (∼6; Handley et al., 2012). Hence, we suggest that the declines in foraminifer abundances represent exclusion due to environmental pressure in combination with sedimentary dilution. As “tropical” foraminifer assemblages appear in higher latitudes during the PETM (Kelly, 2002) and culture experiments demonstrate upper temperature tolerances of ∼33 °C (Hemleben et al., 1989), extreme temperatures are likely to have been the principal environmental agent driving the ecological changes captured at TDP-14. The most extreme ecological conditions and the calcareous plankton exclusion occur higher in the core than the geochemical and biostratigraphic evidence for the onset of the CIE. This may reflect a slow response to the initial forcing, whereby plankton populations were able to tolerate environmental change for a significant period of time before tolerance thresholds were breached.
TDP-14 contains exceptionally well-preserved calcareous plankton, and although the stratigraphy is complex, a number of biotic and geochemical records document the environmental perturbation and biotic responses of the PETM at low latitudes. δ18O from some single-specimen planktonic foraminifera display low values (∼–5‰) and yield paleotemperatures in excess of 40 °C when corrected for paleolatitude and pH effects. A major decline in calcareous plankton abundance throughout the PETM interval likely represents exclusion due to environmental pressure in combination with sedimentary dilution. It is possible that SSTs near TDP-14 may have been even higher than have been inferred during the intervals of the PETM when foraminifera were absent.
Aze and Pearson were supported by the DEFRA Ocean Acidification Programme and the Natural Environment Research Council (NERC; grant IP-1275-1111). Dickson, Cohen, and Coe were supported by NERC grant NE/F021313/1. We thank A. Kuhl of the NERC Life Sciences Mass Spectrometry Facility for assistance with δ13Calk. Badger and Pancost were supported by NERC.