Global cooling during the Eocene Oligocene Transition

About 34 million years ago, Earth's climate shifted from a relatively ice-free world to one with glacial conditions on Antarctica characterized by substantial ice sheets. How Earth's temperature changed during this climate transition remains poorly understood, and evidence for Northern Hemisphere polar ice is controversial. Here, we report proxy records of sea surface temperatures from multiple ocean localities and show that the high-latitude temperature decrease was substantial and heterogeneous. High-latitude (45 degrees to 70 degrees in both hemispheres) temperatures before the climate transition were ~20°C and cooled an average of ~5°C. Our results, combined with ocean and ice-sheet model simulations and benthic oxygen isotope records, indicate that Northern Hemisphere glaciation was not required to accommodate the magnitude of continental ice growth during this time.

Chert intervals in DSDP and ODP sites (Table 1)

Radiolarian cherts in the Tethyan realm of Jurassic age were recently interpreted as resulting from high biosiliceous productivity along upwelling zones in subequatorial paleolatitudes the locations of which were confirmed by revised paleomagnetic estimates. However, the widespread occurrence of cherts in the Eocene suggests that cherts may not always be reliable proxies of latitude and upwelling zones. In a new survey of the global spatio-temporal distribution of Cenozoic cherts in Deep Sea Drilling Project (DSDP) and Ocean Drilling Program (ODP) sediment cores, we found that cherts occur most frequently in the Paleocene and early Eocene, with a peak in occurrences at ~50 Ma that is coincident with the time of highest bottom water temperatures of the early Eocene climatic optimum (EECO) when the global ocean was presumably characterized by reduced upwelling efficiency and biosiliceous productivity. Cherts occur less commonly during the subsequent Eocene global cooling trend. Primary paleoclimatic factors rather than secondary diagenetic processes seem therefore to control chert formation. This timing of peak Eocene chert occurrence, which is supported by detailed stratigraphic correlations, contradicts currently accepted models that involve an initial loading of large amounts of dissolved silica from enhanced weathering and/or volcanism in a supposedly sluggish ocean of the EECO, followed during the subsequent middle Eocene global cooling by more vigorous oceanic circulation and consequent upwelling that made this silica reservoir available for enhanced biosilicification, with the formation of chert as a result of biosilica transformation during diagenesis. Instead, we suggest that basin-basin fractionation by deep-sea circulation could have raised the concentration of EECO dissolved silica especially in the North Atlantic, where an alternative mode of silica burial involving widespread direct precipitation and/or absorption of silica by clay minerals could have been operative in order to maintain balance between silica input and output during the upwelling-deficient conditions of the EECO. Cherts may therefore not always be proxies of biosiliceous productivity associated with latitudinally focused upwelling zones.

Concentrations and isotopic ratios of thorium and uranium in pelagic sediments

The predictable in situ production of 230Th from the decay of uranium in seawater, and its subsequent removal by scavenging onto falling particles, provides a valuable tool for normalizing fluxes to the seafloor. We describe a new application, determination of the 232Th that dissolves in the water column and is removed to the seafloor. 232Th is supplied to the ocean in continental minerals, dissolution of which leads to a measurable standing stock in the water column. Sedimentary adsorbed 232Th/230Th ratios have the potential to provide a proxy for estimating the amount of dissolved material that enters the ocean, both today and in the past. Ten core top samples were treated with up to eight different leaching techniques in order to determine the best method for the separating adsorbed from lattice bound thorium. In addition, separate components of the sediments were analyzed to test whether clay dissolution was an important contribution to the final measurement. There was no systematic correlation between the strength of acid used in the leach and the measured 232Th/230Th ratios. In all cases clean foraminifera produced the same ratio as leaches on bulk sediment. In three out of five samples leaches performed on non-carbonate detritus in the <63 µm size fraction were also identical. Without additional water column data it is not yet clear whether there is a simple one to one correlation between the expected deep-water 232Th/230Th and that produced by leaching, especially in carbonate-rich sediments. However, higher ratios, and associated high 232Th adsorbed fluxes, were observed in areas with high expected detrital inputs. The adsorbed fraction was ~35-50% of the total 232Th in seven out of ten samples. Our 230Th normalized 232Th fluxes are reasonable by comparison to global estimates of detrital inputs to the ocean. In nine cases out of ten, the total 230Th-normalized 232Th flux is greater than predicted from the annual dust fall at each specific location, but lower than the average global detrital input from all sources.

Strontium isotope ratios and concentrations in fossil fish teeth

We analyzed 87Sr/86Sr ratios in foraminifera, pore fluids, and fish teeth for samples ranging in age from Eocene to Pleistocene from four Ocean Drilling Program sites distributed around the globe: Site 1090 in the Cape Basin of the Southern Ocean, Site 757 on the Ninetyeast Ridge in the Indian Ocean, Site 807 on the Ontong-Java Plateau in the western equatorial Pacific, and Site 689 on the Maud Rise in the Southern Ocean. Sr isotopic ratios for dated foraminifera consistently plot on the global seawater Sr isotope curve. For Sites 1090, 757, and 807 Sr isotopic values of the pore fluids are generally less radiogenic than contemporaneous seawater values, as are values for fossil fish teeth. In contrast, pore fluid 87Sr/86Sr values at Site 689 are more radiogenic than contemporaneous seawater, and the corresponding fish teeth also record more radiogenic values. Thus, Sr isotopic values preserved in fossil fish teeth are consistently altered in the direction of the pore fluid values; furthermore, there is a correlation between the magnitude of the offset between the pore fluids and the seawater curve, and the associated offset between the fish teeth and the seawater curve. These data suggest that the hydroxyfluorapatite of the fossil fish teeth continues to recrystallize and exchange Sr with its surroundings during burial and diagenesis. Therefore, Sr chemostratigraphy can be used to determine rough ages for fossil fish teeth in these cores, but cannot be used to fine-tune age models. In contrast to the Sr isotopic system, our Nd concentration data, combined with published isotopic and rare earth element data, suggest that fish teeth acquire Nd during early diagenesis while they are still in direct contact with seawater. The concentrations of Nd acquired at this stage are extremely high relative to the concentrations in surrounding pore fluids. As a result, Nd isotopes are not altered during burial and later diagenesis. Therefore, fossil fish teeth from a variety of marine environments preserve a reliable and robust record of deep seawater Nd isotopic compositions from the time of deposition.

Distribution of planktonic foraminifera and paleotemperature reconstruction for the Subantarctic Zone of the South Atlantic

Pleistocene summer sea-surface temperatures (SSST) have been reconstructed on a composite core section recovered in the Subantarctic Zone of the Southern Ocean from planktonic foraminifers applying the Modern Analog Technique. The composite consists of Core PS2489-2 and the sections recovered at ODP Site 1090, and documents the last 1.83 Ma. Three distinct climatic periods can be identified that mirror the Pleistocene development of the Southern Ocean hydrography. Cold climatic conditions prevailed at 43°S during glacial as well as during interglacial periods during the early Pleistocene (1.83-0.87 Ma), indicating a northward shift of isotherms that characterize the present-day Polar Front Zone by about 7° of latitude. Evidence shows a strong linkage between Southern Ocean and low latitude climate during that interval time. Between the Mid-Pleistocene Revolution (ca. 0.9 Ma) and the Mid-Brunhes Event (ca. 0.4 Ma), we observe higher amplitude fluctuations in the SSST between glacial and interglacial periods, corresponding to the temperature range between the present Polar Front and Subantarctic Front. These climatic variations have been related to changes in the northern hemisphere ice sheets. The past 0.4 Ma are characterized by strong SSST variations, of up to 8°C, between glacials and interglacials. Only during the climatic optima (stages 11.3, 9.3, 7.5, 7.1, 5.5, and the early Holocene), SSST exceeded present SSST at the core locality (10.2°C). Although the carbonate dissolution record exhibits high variability during the Pleistocene, it can be shown that SSST estimates were not significantly biased. The Mid-Brunhes dissolution cycle as well as the Mid-Pleistocene enhanced carbonate preservation appear to belong to a global long-term variability in carbonate preservation.

Nd isotopes an geochemistry of bulk deep sea sediments of the Atlantic Ocean

Nd isotopes preserved in fossil fish teeth and ferromanganese crusts have become a common tool for tracking variations in water mass composition and circulation through time. Studies of Nd isotopes extracted from Pleistocene to Holocene bulk sediments using hydroxylamine hydrochloride (HH) solution yield high resolution records of Nd isotopes that can be interpreted in terms of deep water circulation, but concerns about diagenesis and potential contamination of the seawater signal limit application of this technique to geologically young samples. In this study we demonstrate that Nd extracted from the > 63 µm, decarbonated fraction of older Ocean Drilling Program (ODP) sediments using a 0.02 M HH solution produces Nd isotopic ratios that are within error of values from cleaned fossil fish teeth collected from the same samples, indicating that the HH-extractions are robust recorders of deep sea Nd isotopes. This excellent correlation was achieved for 94 paired fish teeth and HH-extraction samples ranging in age from the Miocene to Cretaceous, distributed throughout the north, tropical and south Atlantic, and composed of a range of lithologies including carbonate-rich oozes/chalks and black shales. The strong Nd signal recovered from Cretaceous anoxic black shale sequences is unlikely to be associated with ferromanganese oxide coatings, but may be derived from abundant phosphatic fish teeth and debris or organic matter in these samples. In contrast to the deep water Nd isotopic signal, Sr isotopes from HH-extractions are often offset from seawater values, suggesting that evaluation of Sr isotopes is a conservative test for the integrity of Nd isotopes in the HH fraction. However, rare earth elements (REE) from the HH-extractions and fish teeth produce distinctive middle REE bulge patterns that may prove useful for evaluating whether the Nd isotopic signal represents uncontaminated seawater. Alternatively, a few paired HH-extraction and cleaned fish teeth samples from each site of interest can be used to verify the seawater composition of the HH-extractions. The similarity between isotopic values for the HH-extraction and fish teeth illustrates that the extensive cleaning protocol applied to fish teeth samples is not necessary in typical, carbonate-rich, deep sea sediments.

Relative abundance and ranges of selected diatoms from Pliocene-Pleistocene sections of ODP Leg 177, Atlantic sector of the Southern Ocean

During Ocean Drilling Program (ODP) Leg 177, seven sites were drilled aligned on a transect across the Antarctic Circumpolar Current in the Atlantic sector of the Southern Ocean. The primary scientific objective of Leg 177 was the study of the Cenozoic paleoceanographic and paleoclimatic history of the southern high latitudes and its relationship with the Antarctic cryosphere development. Of special emphasis was the recovery of Pliocene-Pleistocene sections, allowing paleoceanographic studies at millennial or higher time resolution, and the establishment of refined biostratigraphic zonations tied to the geomagnetic polarity record and stable isotope records. At most sites, multiple holes were drilled to ensure complete recovery of the section. A description of the recovered sections and the construction of a multihole splice for the establishment of a continuous composite is presented in the Leg 177 Initial Reports volume for each of the sites (Gersonde, Hodell, Blum, et al., 1999). Here we present the relative abundance pattern and the stratigraphic ranges of diatom taxa encountered from shore-based light microscope studies completed on the Pliocene-Pleistocene sequences from six of the drilled sites (Sites 1089-1094). No shore-based diatom studies have been conducted on the Pliocene-Pleistocene sediments obtained at Site 1088, located on the northern crest of the Agulhas Ridge, because of the scattered occurrence and poor preservation of diatoms in these sections (Shipboard Scientific Party, 1999b). The data included in our report present the baseline of a diatom biostratigraphic study of Zielinski and Gersonde (2002), which (1) includes a refinement of the southern high-latitude Pliocene-Pleistocene diatom zonation, in particular for the middle and late Pleistocene, and (2) presents a biostratigraphic framework for the establishment of age models of the recovered sediment sections. Zielinski and Gersonde (2002) correlated the diatom ranges with the geomagnetic polarity record established shipboard (Sites 1090 and 1092) (Shipboard Scientific Party, 1999c, 1999d) and on shore (Sites 1089, 1091, 1093, and 1094) by Channell and Stoner (2002). The Pliocene-Pleistocene diatom zonation proposed by Zielinski and Gersonde (2002) relies on a diatom zonation from Gersonde and Bárcena (1998) for the northern belt of the Southern Ocean. Because of latitudinal differentiation of sea-surface temperature, nutrients, and salinity between Antarctic and Subantarctic/subtropical water masses, the Pliocene-Pleistocene stratigraphic marker diatoms are not uniformly distributed in the Southern Ocean (Fenner, 1991; Gersonde and Bárcena, 1998). As a consequence, Zielinski and Gersonde (2002) propose two diatom zonations for application in the Antarctic Zone south of the Polar Front (Southern Zonation, Sites 1094 and 1093) and the area encompassing the Polar Front Zone (PFZ) and the Subantarctic Zone (Northern Zonation, Sites 1089-1092). This accounts especially for the Pleistocene zonation where Hemidiscus karstenii, whose first abundant occurrence datum and last occurrence datum defines the subzonation of the northern Thalassiosira lentiginosa Zone, occurs only sporadically in the cold-water realm south of the PFZ and thus is not applicable in sections from this area. However, newly established marker species assigned to the genus Rouxia (Rouxia leventerae and Rouxia constricta) are more related to cold-water environments and allow a refinement of the Pleistocene stratigraphic zonation for the southern cold areas. A study relying on quantitative counts of both Rouxia species confirms the utility of these stratigraphic markers for the identification of sequences attributed to marine isotope Stages 6 and 8 in the southern Southern Ocean (Zielinski et al., 2002).

Miocene to Pliocene calcareous nannofossil biostratigraphy at ODP Leg 117 Sites 1088 and 1090, Southern Ocean

The calcareous nannofossil biostratigraphy of ODP Leg 177 Sites 1088 and 1090 (Subantarctic sector from the Atlantic Ocean) is discussed. Most nannofossil zonal boundaries of Martini (1971) and Okada and Bukry (1980) were recognized for the studied mid-high-latitude sediments. Conventional low-latitude marker species such as Amaurolithus spp., Discoaster spp., Triquetrorhabdulus spp., Ceratolithus spp. were recorded as rare and scattered, which impeded the development of a detailed nannofossil biostratigraphic zonation of some Miocene and Pliocene intervals. Because of the absence of some primary biostratigraphic marker species, additional second-order bioevents, such as the first occurrence of Calcidiscus macintyrei and the last occurrence of Coccolithus miopelagicus, have been used to approximate the base of zones NN7 and NN8, respectively. Several disconformities disturbing the Pliocene and Miocene intervals of Site 1090 could be determined based on nannofossil distribution although the occurrence of intervals with dissolved nannofloras and low species diversity prevented a reliable age assignment. An acme of small Gephyrocapsa was recognized near the lower/middle Pliocene boundary, close to the NN15-NN16 zonal boundary, presenting an event for further improvement of the calcareous nannofossil biostratigraphy of this interval time. The first occurrence of Pseudoemiliania lacunosa (>4 µm) occurs close to this interval, representing a fairly reliable event to approximate the base of NN15 zone when other biozonal events are absent. A paracme of R. pseudoumbilicus (>7 µm) was detected in the lower Pliocene NN12 and in the upper Miocene NN11. These temporary absences of the species seem to be isochronous between high-latitude and low-middle-latitude areas.

Calcareous nannofossils of ODP Hole 177-1090D

Calcareous nannofossil assemblages have been investigated at Ocean Drilling Program (ODP) Site 1090 located in the modern Subantarctic Zone, through the Pleistocene Marine Isotope Stages (MIS) 34-29, between 1150 and 1000 ka. A previously developed age model and new biostratigraphic constraints provide a reliable chronological framework for the studied section and allow correlation with other records. Two relevant biostratigraphic events have been identified: the First Common Occurrence of Reticulofenestra asanoi, distinctly correlated to MIS 31-32; the re-entry of medium Gephyrocapsa at MIS 29, unexpectedly similar to what was observed at low latitude sites. The composition of the calcareous nannofossil assemblage permits identification of three intervals (I-III). Intervals I and III, correlated to MIS 34-32 and MIS 30-29 respectively, are identified as characteristic of water masses located south of the Subtropical Front and reflecting the southern border of Subantarctic Zone, at the transition with the Polar Front Zone. This evidence is consistent with the hypothesis of a northward shift of the frontal system in the early Pleistocene with respect to the present position and therefore a northernmost location of the Subantarctic Front. During interval II, which is correlated to MIS 31, calcareous nannofossil assemblages display the most significant change, characterized by a distinct increase of Syracosphaera spp. and Helicosphaera carteri, lasting about 20 ky. An integrated analysis of calcareous nannofossil abundances and few mineralogical proxies suggests that during interval II, Site 1090 experienced the influence of subtropical waters, possibly related to a southward migration of the Subtropical Front, coupled with an expansion of the warmer Agulhas Current at the core location. This pronounced warming event is associated to a minimum in the austral summer insolation. The present results provide a broader framework on the Mid-Pleistocene dynamic of the ocean frontal system in the Atlantic sector of the Southern Ocean, as well as additional evidence on the variability of the Indian-Atlantic ocean exchange.

Ice-rafted debris in ODP Holes 177-1090D and 188-1165B

Constraining the nature of Antarctic Ice Sheet (AIS) response to major past climate changes may provide a window onto future ice response and rates of sea level rise. One approach to tracking AIS dynamics, and differentiating whole system versus potentially heterogeneous ice sheet sector changes, is to integrate multiple climate proxies for a specific time slice across widely distributed locations. This study presents new iceberg-rafted debris (IRD) data across the interval that includes Marine Isotope Stage 31 (MIS 31: 1.081-1.062 Ma, a span of ~19 kyr; Lisiecki and Raymo, 2005), which lies on the cusp of the mid-Brunhes climate transition (as glacial cycles shifted from ~41,000 yr to ~100,000 yr duration). Two sites are studied - distal Ocean Drilling Program (ODP) Leg 177 Site 1090 (Site 1090) in the eastern subantarctic sector of the South Atlantic Ocean, and proximal ODP Leg 188 Site 1165 (Site 1165), near Prydz Bay, in the Indian Ocean sector of the Antarctic margin. At each of these sites, MIS 31 is marked by the presence of the Jaramillo Subchron (0.988-1.072 Ma; Lourens et al., 2004) which provides a time-marker to correlate these two sites with relative precision. At both sites, records of multiple climate proxies are available to aid in interpretation. The presence of IRD in sediments from our study areas, which include garnets indicating a likely East Antarctic Ice Sheet (EAIS) origin, supports the conclusion that although the EAIS apparently withdrew significantly over MIS 31 in the Prydz Bay region and other sectors, some sectors of the EAIS must still have maintained marine margins capable of launching icebergs even through the warmest intervals. Thus, the EAIS did not respond in complete synchrony even to major climate changes such as MIS 31. Further, the record at Site 1090 (supported by records from other subantarctic locations) indicates that the glacial MIS 32 should be reduced to no more than a stadial, and the warm interval of Antarctic ice retreat that includes MIS 31 should be expanded to MIS 33-31. This revised warm interval lasted about 52 kyr, in line with several other interglacials in the benthic d18O records stack of Lisiecki and Raymo (2005), including the super-interglacials MIS 11 (duration of 50 kyr) and MIS 5 (duration of 59 kyr). The record from Antarctica-proximal Site 1165, when interpreted in accord with the record from ANDRILL-1B, indicates that in these southern high latitude sectors, ice sheet retreat and the effects of warming lasted longer than at Site 1090, perhaps until MIS 27. In the current interpretations of the age models of the proximal sites, ice sheet retreat began relatively slowly, and was not really evident until the start of MIS 31. In another somewhat more speculative interpretation, ice sheet retreat began noticeably with MIS 33, and accelerated during MIS 31. Ice sheet inertia (the lag-times in the large-scale responses of major ice sheets to a forcing) likely plays an important part in the timing and scale of these events in vulnerable sectors of the AIS.