Analysis of an ancient global warming event, the Paleocene-Eocene thermal maximum, under examination from ODP Holes 113-689B and 113-690B

Spatiotemporal patterns of carbonate dissolution provide a critical constraint on carbon input during an ancient (~55.5 Ma) global warming event known as the Paleocene-Eocene thermal maximum (PETM), yet the magnitude of lysocline shoaling in the Southern Ocean is poorly constrained due to limited spatial coverage in the circum-Antarctic region. This shortcoming is partially addressed by comparing patterns of carbonate sedimentation at the Site 690 PETM reference section to those herein reconstructed for nearby Site 689. Biochemostratigraphic correlation of the two records reveals that the first ~36 ka of the carbon isotope excursion (CIE) signaling PETM conditions is captured by the Site 689 section, while the remainder of the CIE interval and nearly all of the CIE recovery are missing due to a coring gap. A relatively expanded stratigraphy and higher carbonate content at mid-bathyal Site 689 indicate that dissolution was less severe than at Site 690. Thus, the bathymetric transect delimited by these two PETM records indicates that the lysocline shoaled above Site 689 (~1,100 m) while the calcite compensation depth remained below Site 690 (~1,900 m) in the Weddell Sea region. The ensuing recovery of carbonate sedimentation conforms to a bathymetric trend best explained by gradual lysocline deepening as negative feedback mechanisms neutralized ocean acidification. Further, biochemostratigraphic evidence indicates the tail end of the CIE recovery interval at both sites has been truncated by a hiatus most likely related to vigorous production and advection of intermediate waters.

Tomographic measurements of O. umbonatus and N. truempyi during the PETM and ETM-2 at ODP Sites 208-1262 and 208-1263

Predicting the impact of ongoing anthropogenic CO2 emissions on calcifying marine organisms is complex, owing to the synergy between direct changes (acidification) and indirect changes through climate change (e.g., warming, changes in ocean circulation, and deoxygenation). Laboratory experiments, particularly on longer-lived organisms, tend to be too short to reveal the potential of organisms to acclimatize, adapt, or evolve and usually do not incorporate multiple stressors. We studied two examples of rapid carbon release in the geological record, Eocene Thermal Maximum 2 (~53.2 Ma) and the Paleocene Eocene Thermal Maximum (PETM, ~55.5 Ma), the best analogs over the last 65 Ma for future ocean acidification related to high atmospheric CO2 levels. We use benthic foraminifers, which suffered severe extinction during the PETM, as a model group. Using synchrotron radiation X-ray tomographic microscopy, we reconstruct the calcification response of survivor species and find, contrary to expectations, that calcification significantly increased during the PETM. In contrast, there was no significant response to the smaller Eocene Thermal Maximum 2, which was associated with a minor change in diversity only. These observations suggest that there is a response threshold for extinction and calcification response, while highlighting the utility of the geological record in helping constrain the sensitivity of biotic response to environmental change.

Sr, Nd and Pb isotopic composition of volcanic and volcaniclastic samples of ODP Leg 134

A suite of volcanic and volcaniclastic rocks selected from Ocean Drilling Program Leg 134 Sites 832 and 833 in the North Aoba Basin (Central New Hebrides Island Arc) has been analyzed for Sr, Nd, and Pb isotopes to investigate the temporal evolution of the arc magmatism. This arc shows two unusual features with respect to other western Pacific arcs: 1) subduction is eastdirected; and 2) a major submarine ridge, the d'Entrecasteaux Zone, has been colliding almost perpendicularly with the central part of the arc since about 3 Ma. Volcanic rocks from the upper parts of both holes, generated during the last 2 m.y., show higher 87Sr/86Sr and significantly lower 206Pb/204Pb and 143Nd/144 Nd values compared to those volcanics erupted before the collision of this ridge, as represented by samples from the lower section of both holes, or remote from the collisional region, in the southern part of the arc. These isotopic differences in the respective mantle sources cannot be interpreted in terms of geochemical input into the mantle wedge induced by the collision itself. Rather, they require long term (>500 m.y.) enrichment processes. The enriched mantle source could be, on a regional scale, a DUPAL-type reservoir with strong similarities to the source of Indian Ocean basalts. Isotopic analyses of drilled rocks from the DEZ show that the anomalous, enriched mantle component is not derived from this feature. We currently cannot identify a source for this enriched component, but note that it also exists in Lau Basin backarc volcanics, lavas from the West Philippine Sea, and also some lavas from the Mariana-Izu-Bonin arc.

Geochemistry of basalts from ODP Hole 113-690C

Basalts from Maud Rise, Weddell Sea, are vesicular and olivine-phyric. Major, trace, and rare earth element concentrations are similar to those of alkali basalts from ocean islands and seamounts. The rocks are low in MgO, Cr, Ni, and Sc, and high in TiO2, K2O, P2O5, Zr, and LREE contents. The abundance of "primary" biotite and apatite in the matrix indicates the melting of a hydrous mantle. Prevalence of olivine and absence of plagioclase in the rocks suggests that the volatile in the melt was an H2O-CO2 mixture, where H2O was <0.5. Mantle derived xenocrysts in the basalt include corroded orthopyroxene, chromite, apatite, and olivine. Olivine (Fo90) is too magnesian to be in equilibrium with the basalts, as they contain only 5-6 wt% MgO. Based on the presence of mantle xenocrysts, the high concentration of incompatible elements, the spatial and chemical affinity with other ocean island basalts from the area, and the relative age of the basalt (overlain by late Campanian sediments), it is suggested that Maud Rise was probably generated by hot-spot activity, possible during a ridge crest jump prior to 84 Ma (anomaly 34 time). Iddingsite, a complex intergrowth of montmorillonite and goethite, is the major alteration product of second generation olivine. It is suggested that iddingsite crystallized at low temperatures (<200°C) from an oxidized fluid during deuteric alteration. Vesicles are commonly filled by zeolites which have been replaced by K-feldspars.

Helium isotope ratios and sedimentation rate models of ODP holes

A continuous age model for the brief climate excursion at the Paleocene-Eocene boundary has been constructed by assuming a constant flux of extraterrestrial 3He (3He[ET]) to the seafloor. 3He[ET] measurements from ODP Site 690 provide quantitative evidence for the rapid onset (age models indicating extremely rapid release of isotopically light carbon, possibly from seafloor methane hydrate, as the proximal cause of the event. However, the 3He[ET] technique indicates a previously unrecognized and extreme increase in sedimentation rate coincident with the return of climate proxies to pre-event values. The 3He[ET]-based age model thus suggests a far more rapid recovery from the climatic perturbation than previously proposed or predicted on the basis of the modern carbon cycle, and so may indicate additional or accelerated mechanisms of carbon removal from the ocean-atmosphere system during this period. 3He[ET] was also measured at ODP Site 1051 to test the validity of the Site 690 chronology. Comparison of these data sets seems to require removal of several tens of kyr of sediment within the climatic excursion at Site 1051, an observation consistent with sediment structures and previous age modeling efforts. The Site 1051 age model shows a ~30 kyr period in which climate proxies return toward pre-event values, after which they remain invariant for ~80 kyr. If this rise represents the recovery interval identified at Site 690, then the 3HeET-based age models of the two sites are in good agreement. However, alternative interpretations are possible, and work on less disrupted sites is required to evaluate the reliability of the proposed new chronology of the climate excursion. Regardless of these details, this work shows that the 3HeET technique can provide useful independent evidence for the development and testing of astronomically calibrated age models.

Faunal turnovers in the central Pacific benthic foraminifera during the Paleocene-Eocene thermal maximum

Although it is well known that the Paleocene/Eocene thermal maximum (PETM) coincided with a major benthic foraminiferal extinction event, the detailed pattern of the faunal turnover has not yet been clarified. Our high-resolution benthic foraminiferal and carbon isotope analyses at the low latitude Pacific Ocean Shatsky Rise have revealed the following record of major faunal transitions: (1) An initial turnover which involved the benthic foraminiferal extinction event (BFE). The BFE, marked by a sharp transition from Pre-extinction fauna to Disaster fauna represented by small-sized Bolivina gracilis, expresses the onset of the PETM and the abrupt extinction of about 30% of taxa. This faunal transition lasted about 45-74 kyr after the initiation of the PETM and was followed by: (2) the appearance of Opportunistic fauna represented by Quadrimorphina profunda, which existed for about 74-91 kyr after the initiation of the PETM. These two faunas, which appeared after the extinction event, are characterized by low diversity and dwarfism, possibly due to lowered oxygen condition and decreased surface productivity. The second pronounced turnover involved the gradual recovery from Opportunistic Fauna to the establishment of Recovery fauna, which coincided with the recovery about 83-91 kyr after its initiation.

Nannofossil and stable isotope record across the Paleocene/Eocene Thermal Maximum in ODP Hole 183-1135A

The Paleocene/Eocene Thermal Maximum (PETM, ca. 55 Ma) is an abrupt, profound perturbation of climate and the carbon cycle associated with a massive injection of isotopically light carbon into the ocean-atmosphere system. As such, it provides an analogue for understanding the interplay between phytoplankton and climate under modern anthropogenic global-warming conditions. However, the accompanying enhanced dissolution poses uncertainty on the reconstruction of the affected ecology and productivity. We present a high-resolution record of bulk isotopes and nannofossil absolute abundance from Ocean Drilling Program (ODP) Site 1135 on the Kerguelen Plateau, Southern Indian Ocean to quantitatively constrain for the first time the influence of dissolution on paleoecological reconstruction. Our bulk-carbonate isotope record closely resembles that of the classic PETM site at ODP Site 690 on the opposite side of the Antarctic continent, and its correlation with those from ODP Sites 690, 1262 and 1263 records allows recognition of 14 precessional cycles upsection from the onset of the carbon isotopic excursion (CIE). This, together with a full range of common Discoasteraraneus and an abundance crossover between Fasciculithus and Zygrhablithusbijugatus, indicates the presence of the PETM at Site 1135, a poorly known record with calcareous fossils throughout the interval. The strong correlation between the absolute abundances of Chiasmolithus and coccolith assemblages reveals a dominant paleoecological signal in the poorly preserved fossil assemblages, while the influence of dissolution is only strong during the CIE. This suggests that r-selected taxa can preserve faithful ecological information even in the severely-altered assemblages studied here, and therefore provide a strong case for the application of nannofossils to paleoecological studies in better-preserved PETM sections. The inferred nannoplankton productivity drops abruptly at the CIE onset, but rapidly increases after the CIE peak, both of which may be driven by nutrient availability related to ocean stratification and vertical mixing due to changed sea-surface temperatures.

Geochemical traces across the Paleocene-Eocene thermal maximum on Walvis Ridge

An understanding of sediment redox conditions across the Paleocene-Eocene thermal maximum (PETM) (?55 Ma) is essential for evaluating changes in processes that control deep-sea oxygenation, as well as identifying the mechanisms responsible for driving the benthic foraminifera extinction. Sites cored on the flanks of Walvis Ridge (Ocean Drilling Program Leg 208, Sites 1262, 1266, and 1263) allow us to examine changes in bottom and pore water redox conditions across a ~2 km depth transect of deep-sea sediments of PETM age recovered from the South Atlantic. Here we present measurements of the concentrations of redox-sensitive trace metals manganese (Mn) and uranium (U) in bulk sediment as proxies for redox chemistry at the sediment-water interface and below. All three Walvis Ridge sites exhibit bulk Mn enrichment factors (EF) ranging between 4 and 12 prior to the warming, values at crustal averages (Mn EF = 1) during the warming interval, and a return to pre-event values during the recovery period. U enrichment factors across the PETM remains at crustal averages (U EF = 1) at Site 1262 (deep) and Site 1266 (intermediate depth). U enrichment factors at Site 1263 (shallow) peaked at 5 immediately prior to the PETM and dropped to values near crustal averages during and after the event. All sites were lower in dissolved oxygen content during the PETM. Before and after the PETM, the deep and intermediate sites were oxygenated, while the shallow site was suboxic. Our geochemical results indicate that oxygen concentrations did indeed drop during the PETM but not sufficiently to cause massive extinction of benthic foraminifera.

Core alignment and composite depth scale of DSDP Site 86-577

Deep Sea Drilling Project Site 577 on Shatsky Rise (North Pacific Ocean) recovered a series of cores at three holes that contain calcareous nannofossil ooze of latest Cretaceous (late Maastrichtian) through early Eocene age. Several important records have been generated using samples from these cores, but the stratigraphy has remained outdated and confusing. Here we revise the stratigraphy at Site 577. This includes refining several age datums, realigning cores in the depth domain, and placing all stratigraphic markers on a current time scale. The work provides a template for appropriately bringing latest Cretaceous and Paleogene data sets at old drill sites into current paleoceanographic literature for this time interval. While the Paleocene Eocene Thermal Maximum (PETM) lies within core gaps at Holes 577* and 577A, the sedimentary record at the site holds other important events and remains crucially relevant to understanding changes in oceanographic conditions from the latest Cretaceous through early Paleogene.

Chemical composition of volcanic glass of some tephra layers in Northeastern Pacific Ocean

Five widespread upper Cenozoic tephra layers that are found within continental sediments of the western United States have been correlated with tephra layers in marine sediments in the Humboldt and Ventura basins of coastal California by similarities in major-and trace-element abundances; four of these layers have also been identified in deep-ocean sediments at DSDP sites 34, 36, 173, and 470 in the northeastern Pacific Ocean. These layers, erupted from vents in the Yellowstone National Park area of Wyoming and Idaho (Y), the Cascade Range of the Pacific Northwest (C), and the Long Valley area, California (L), are the Huckleberry Ridge ash bed (2.0 Ma, Y), Rio Dell ash bed (ca. 1.5 Ma, C), Bishop ash bed (0.74 Ma, L), Lava Creek B ash bed (0.62 Ma, Y), and Loleta ash bed (ca. 0.4 Ma, C). The isochronous nature of these beds allows direct comparison of chronologic and climatic data in a variety of depositional environments. For example, the widespread Bishop ash bed is correlated from proximal localities near Bishop in east-central California, where it is interbedded with volcanic and glacial deposits, to lacustrine beds near Tecopa, southeastern California, to deformed on-shore marine strata near Ventura, southwestern California, to deep-ocean sediments at site 470 in the eastern Pacific Ocean west of northern Mexico. The correlations allow us to compare isotopic ages determined for the tephra layers with ages of continental and marine biostratigraphic zones determined by magnetostratigraphy and other numerical age control and also provide iterative checks for available age control. Relative age variations of as much as 0.5 m.y. exist between marine biostratigraphic datums [for example, highest occurrence level of Discoaster brouweri and Calcidiscus tropicus (= C. macintyrei)], as determined from sedimentation rate curves derived from other age control available at each of several sites. These discrepancies may be due to several factors, among which are (1) diachronism of the lowest and highest occurrence levels of marine faunal and floral species with latitude because of ecologic thresholds, (2) upward reworking of older forms in hemipelagic sections adjacent to the tectonically active coast of the western United States and other similar analytical problems in identification of biostratigraphic and magnetostratigraphic datums, (3) dissolution of microfossils or selective diagenesis of some taxa, (4) lack of precision in isotopic age calibration of these datums, (5) errors in isotopic ages of tephra beds, and (6) large variations in sedimentation rates or hiatuses in stratigraphic sections that result in age errors of interpolated datums. Correlation of tephra layers between on-land marine and deep-ocean deposits indicates that some biostratigraphic datums (diatom and calcareous nannofossil) may be truly time transgressive because at some sites, they are found above and, at other sites, below the same tephra layers.

Atmospheric carbon dioxide, methane, deuterium, and calculated Antarctic temperature of EPICA Dome C ice core

Understanding the role of atmospheric CO2 during past climate changes requires clear knowledge of how it varies in time relative to temperature. Antarctic ice cores preserve highly resolved records of atmospheric CO2 and Antarctic temperature for the past 800,000 years. Here we propose a revised relative age scale for the concentration of atmospheric CO2 and Antarctic temperature for the last deglacial warming, using data from five Antarctic ice cores. We infer the phasing between CO2 concentration and Antarctic temperature at four times when their trends change abruptly. We find no significant asynchrony between them, indicating that Antarctic temperature did not begin to rise hundreds of years before the concentration of atmospheric CO2, as has been suggested by earlier studies.

Total number of benthic foraminifera in sediments below the K/T boundary in the Atlantic and Pacific Ocean

Most species of Late Cretaceous deep-sea benthic foraminifera are believed to be cosmopolitan and therefore to exhibit only minor biogeographical differences. In this preliminary report, six Deep Sea Drilling Project (DSDP) sites from different oceans, paleolatitudes, and paleodepths were analyzed for terminal Cretaceous abyssal-bathyal benthic foraminifera in order to investigate their assumed cosmopolitan distribution and the question of whether different faunal compositions are related to time, different paleolatitudes, and/or different paleodepths. The material studied was obtained from the low-latitude Site 465 (Pacific Ocean), and the intermediate-latitude Sites 384 (North Atlantic) and 356, 516, 525, and 527 (South Atlantic). The material analyzed represents a time slice encompassing the last 20-50 k.y. of the Cretaceous. The faunas contain numerous "Velasco-type" species, such as Gavelinella beccariiformis (White), Cibicidoides velascoensis (Cushman), Nuttallides truempyi (Nuttall), Gaudryina pyramidata Cushman, and various gyroidinoids and buliminids. The results contradict the general assumption of the cosmopolitan nature of Late Cretaceous deep-sea benthic foraminifera advocated in the literature. Only about 9% of the taxa identified were found to be truly "cosmopolitan" through their occurrence at all the sites analyzed. On the basis of correspondence analysis and relative abundance data, three assemblages and three subassemblages were recognized: (1) a bathyal-abyssal assemblage [Nuttallinella sp. A, Cibicidoides hyphalus (Fisher), Valvalabamina sp. evolute form, and Gyroidinoides spp.] at the South Atlantic Sites 356, 516, 525, and 527, divided into three subassemblages, namely (a) a middle bathyal subassemblage [Eouvigerina subsculptura McNeil and Caldwell, Truaxia aspera (Cushman), and G. pyramidata] at Sites 516 and 525, (b) a lower bathyal subassemblage [Osangularia? sp., Pyramidina rudita (Cushman and Parker), and Quadrimorphina camerata (Brotzen)] at Site 356, and (c) an abyssal subassemblage [Gyroidinoides sp. C, Hyperammina-Bathysiphon, Gyroidinoides beisseli (White), and Globorotalites sp. B] at Site 527; (2) an abyssal assemblage [Buliminella cf. plana (Cushman and Parker) and Bulimina incisa Cushman] at the North Atlantic Site 384; and (3) a middle bathyal assemblage [Vulvulina sp. A, Osangularia navarroana (Cushman), Alabamina? sp., Bulimina velascoensis (Cushman), Spiroplectammina spp. calcareous forms, and Bulimina trinitatensis Cushman and Jarvis] at the Pacific Site 465.

Barium intensities and tie points between holes of Leg 208

The Paleocene - Eocene thermal maximum (PETM) is one of the best known examples of a transient climate perturbation, associated with a brief, but intense, interval of global warming and a massive perturbation of the global carbon cycle from injection of isotopically light carbon into the ocean-atmosphere system. One key to quantifying the mass of carbon released, identifying the source(s), and understanding the ultimate fate of this carbon is to develop high-resolution age models. Two independent strategies have been employed, cycle stratigraphy and analysis of extraterrestrial Helium (HeET), both of which were first tested on Ocean Drilling Program (ODP) Site 690. Both methods are in agreement for the onset of the PETM and initial recovery, or the clay layer ("main body"), but seem to differ in the final recovery phase of the event above the clay layer, where the carbonate contents rise and carbon isotope values return toward background values. Here we present a state-of-the-art age model for the PETM derived from a new orbital chronology developed with cycle stratigraphic records from sites drilled during ODP Leg 208 (Walvis Ridge, Southeastern Atlantic) integrated with published records from Site 690 (Weddell Sea, Southern Ocean, ODP Leg 113). During Leg 208, five Paleocene - Eocene (P-E) boundary sections (Sites 1262 to 1267) were recovered in multiple holes over a depth transect of more than 2200 m at the Walvis Ridge yielding the first stratigraphically complete P-E deep-sea sequence with moderate to relatively high sedimentation rates (1 to 3 cm/kyr). A detailed chronology was developed with non-destructive X-ray fluorescence (XRF) core scanning records on the scale of precession cycles, with a total duration of the PETM now estimated to be ~ 170 kyr. The revised cycle stratigraphic record confirms original estimates for the duration of the onset and initial recovery, but suggests a new duration for the final recovery that is intermediate to the previous estimates by cycle stratigraphy and HeET.

Natural remanent magnetization (NRM) from ODP Site 1263 covering magnetochron C20r and C21n and high-resolution bulk carbon isotope (d13C) records from Sites 702 and 1263

To explore cause and consequences of past climate change, very accurate age models such as those provided by the astronomical timescale (ATS) are needed. Beyond 40 million years the accuracy of the ATS critically depends on the correctness of orbital models and radioisotopic dating techniques. Discrepancies in the age dating of sedimentary successions and the lack of suitable records spanning the middle Eocene have prevented development of a continuous astronomically calibrated geological timescale for the entire Cenozoic Era. We now solve this problem by constructing an independent astrochronological stratigraphy based on Earth's stable 405 kyr eccentricity cycle between 41 and 48 million years ago (Ma) with new data from deep-sea sedimentary sequences in the South Atlantic Ocean. This new link completes the Paleogene astronomical timescale and confirms the intercalibration of radioisotopic and astronomical dating methods back through the Paleocene-Eocene Thermal Maximum (PETM, 55.930 Ma) and the Cretaceous-Paleogene boundary (66.022 Ma). Coupling of the Paleogene 405 kyr cyclostratigraphic frameworks across the middle Eocene further paves the way for extending the ATS into the Mesozoic.

Geochemical data derived from an absolutely dated shell-based chronology for the Irish Sea

The identification in various proxy records of periods of rapid (decadal scale) climate change over recent millennia, together with the possibility that feedback mechanisms may amplify climate system responses to increasing atmospheric CO2, highlights the importance of a detailed understanding, at high spatial and temporal resolutions, of forcings and feedbacks within the system. Such an understanding has hitherto been limited because the temperate marine environment has lacked an absolute timescale of the kind provided by tree-rings for the terrestrial environment and by corals for the tropical marine environment. Here we present the first annually resolved, multi-centennial (489-year), absolutely dated, shell-based marine master chronology. The chronology has been constructed by detrending and averaging annual growth increment widths in the shells of multiple specimens of the very long-lived bivalve mollusc Arctica islandica, collected from sites to the south and west of the Isle of Man in the Irish Sea. The strength of the common environmental signal expressed in the chronology is fully comparable with equivalent statistics for tree-ring chronologies. Analysis of the 14C signal in the shells shows no trend in the marine radiocarbon reservoir correction (DR), although it may be more variable before ~1750. The d13C signal shows a very significant (R**2 = 0.456, p < 0.0001) trend due to the 13C Suess effect.