Middle to late Miocene oxygen isotope stratigraphy of ODP site 1085 (SE Atlantic): new constrains on Miocene climate variability and sea-level fluctuations

The middle Miocene delta18O increase represents a fundamental change in earth's climate system due to a major expansion and permanent establishment of the East Antarctic Ice Sheet accompanied by some effect of deepwater cooling. The long-term cooling trend in the middle to late Miocene was superimposed by several punctuated periods of glaciations (Mi-Events) characterized by oxygen isotopic shifts that have been related to the waxing and waning of the Antarctic ice-sheet and bottom water cooling. Here, we present a high-resolution benthic stable oxygen isotope record from ODP Site 1085 located at the southwestern African continental margin that provides a detailed chronology for the middle to late Miocene (13.9-7.3 Ma) climate transition in the eastern South Atlantic. A composite Fe intensity record obtained by XRF core scanning ODP Sites 1085 and 1087 was used to construct an astronomically calibrated chronology based on orbital tuning. The oxygen isotope data exhibit four distinct delta18O excursions, which have astronomical ages of 13.8, 13.2, 11.7, and 10.4 Ma and correspond to the Mi3, Mi4, Mi5, and Mi6 events. A global climate record was extracted from the oxygen isotopic composition. Both long- and short-term variabilities in the climate record are discussed in terms of sea-level and deep-water temperature changes. The oxygen isotope data support a causal link between sequence boundaries traced from the shelf and glacioeustatic changes due to ice-sheet growth. Spectral analysis of the benthic delta18O record shows strong power in the 400-kyr and 100-kyr bands documenting a paleoceanographic response to eccentricity-modulated variations in precession. A spectral peak around 180-kyr might be related to the asymmetry of the obliquity cycle indicating that the response of the dominantly unipolar Antarctic ice-sheet to obliquityinduced variations probably controlled the middle to late Miocene climate system. Maxima in the delta18O record, interpreted as glacial periods, correspond to minima in 100-kyr eccentricity cycle and minima in the 174-kyr obliquity modulation. Strong middle to late Miocene glacial events are associated with 400-kyr eccentricity minima and obliquity modulation minima. Thus, fluctuations in the amplitude of obliquity and eccentricity seem to be the driving force for the middle to late Miocene climate variability.

Benthic foraminifera and sea level analyses from IODP Hole 310-M0005D

The course of sea-level fluctuations during Termination II (TII; the penultimate deglaciation), which is critical for understanding ice-sheet dynamics and suborbital climate variability, has yet to be established. This is partly because most shallow-water sequences encompassing TII were eroded during sea-level lowstands of the last glacial period or were deposited below the present sea level. Here we report a new sequence recording sea-level changes during TII in the Pleistocene sequence at Hole M0005D (water depth: 59.63 m below sea level [mbsl]) off Tahiti, French Polynesia, which was drilled during Integrated Ocean Drilling Program Expedition 310. Lithofacies variations and stratigraphic changes in the taxonomic composition, preservation states, and intraspecific test morphology of large benthic foraminifers indicate a deepening-upward sequence in the interval from Core 310-M0005D-26R (core depth: 134 mbsl) through -16R (core depth: 106 mbsl). Reconstruction of relative sea levels, based on paleodepth estimations using large benthic foraminifers, indicated a rise in sea level of about 90 m during this interval, suggesting its correlation with one of the terminations. Assuming that this rise in sea level corresponds to that during TII, after correcting for subsidence since the time of deposition, a highstand sea-level position would be 2 ± 15 m above present sea level (masl), which is generally consistent with highstand sea-level positions in MIS 5e (4 ± 2 masl). If this rise in sea level corresponds to that during older terminations, the subsidence-corrected highstand sea-level positions (30 ± 15 masl for Termination III and 54 ± 15 masl for Termination IV) are not consistent with reported ranges of interglacial sea-level highstands (-18 to 15 masl). Therefore, the studied interval likely records the rise in sea level and associated environmental changes during TII. In particular, the intervening cored materials between the two episodes of sea-level rise found in the studied interval might record the sea-level reversal event during TII. This conclusion is consistent with U/Th ages of around 133 ka, which were obtained from slightly diagenetically altered (i.e., < 1% calcite) in situ corals in the studied interval (Core 310-M0005D-20R [core depth: 118 mbsl]). This study also suggests that our inverse approach to correlate a stratigraphic interval with an approximate time frame could be useful as an independent check on the accuracy of uranium-series dating, which has been applied extensively to fossil corals in late Quaternary sea-level studies.

Orbital-resolution Pliocene-Pleistocene simulations of CO2, sea level, global temperature and benthic d18O, and d11B-based proxy-CO2 data

During the past five million yrs, benthic d18O records indicate a large range of climates, from warmer than today during the Pliocene Warm Period to considerably colder during glacials. Antarctic ice cores have revealed Pleistocene glacial-interglacial CO2 variability of 60-100 ppm, while sea level fluctuations of typically 125 m are documented by proxy data. However, in the pre-ice core period, CO2 and sea level proxy data are scarce and there is disagreement between different proxies and different records of the same proxy. This hampers comprehensive understanding of the long-term relations between CO2, sea level and climate. Here, we drive a coupled climate-ice sheet model over the past five million years, inversely forced by a stacked benthic d18O record. We obtain continuous simulations of benthic d18O, sea level and CO2 that are mutually consistent. Our model shows CO2 concentrations of 300 to 470 ppm during the Early Pliocene. Furthermore, we simulate strong CO2 variability during the Pliocene and Early Pleistocene. These features are broadly supported by existing and new d11B-based proxy CO2 data, but less by alkenone-based records. The simulated concentrations and variations therein are larger than expected from global mean temperature changes. Our findings thus suggest a smaller Earth System Sensitivity than previously thought. This is explained by a more restricted role of land ice variability in the Pliocene. The largest uncertainty in our simulation arises from the mass balance formulation of East Antarctica, which governs the variability in sea level, but only modestly affects the modeled CO2 concentrations.