1.2 Ma stacked benthic foraminiferal carbon isotope record

Pleistocene stable carbon isotope (d13C) records from surface and deep dwelling foraminifera in all major ocean basins show two distinct long-term carbon isotope fluctuations since 1.00 Ma. The first started around 1.00 Ma and was characterised by a 0.35 per mil decrease in d13C values until 0.90 Ma, followed by an increase of 0.60 per mil lasting until 0.50 Ma. The subsequent fluctuation started with a 0.40 per mil decrease between 0.50 and 0.25 Ma, followed by an increase of 0.30 per mil between 0.25 and 0.10 Ma. Here, we evaluate existing evidence and various hypotheses for these global Pleistocene d13C fluctuations and present an interpretation, where the fluctuations most likely resulted from concomitant changes in the burial fluxes of organic and inorganic carbon due to ventilation changes and/or changes in the production and export ratio. Our model indicates that to satisfy the long-term 'stability' of the Pleistocene lysocline, the ratio between the amounts of change in the organic and inorganic carbon burial fluxes would have to be close to a 1:1 ratio, as deviations from this ratio would lead to sizable variations in the depth of the lysocline. It is then apparent that the mid-Pleistocene climate transition, which, apart from the glacial cycles, represents the most fundamental change in the Pleistocene climate, was likely not associated with a fundamental change in atmospheric pCO2. While recognising that high frequency glacial/interglacial cycles are associated with relatively large (100 ppmv) changes in pCO2, our model scenario (with burial changes close to a 1:1 ratio) produces a maximum long-term variability of only 20 ppmv over the fluctuation between 1.00 and 0.50 Ma.

Stable isotope record and age determination of ODP Site 121-758 in the northeastern Indian Ocean

Numerous studies have shown that delta18O records from benthic and planktonic foraminifera, primarily a proxy of global ice volume variations, reflect Milankovitch periodicities. To study climatic response to orbital forcing at Ocean Drilling Program site 758, we have generated continuous delta18O and delta13C records from a single benthic foraminiferal species Cibicides wuellerstorfi for the last 3.6 m.y. and extended the planktonic foraminiferal isotope records of Farrell and Janecek (1991, doi:10.2973/odp.proc.sr.121.124.1991) (0-2.5 Ma, based on Globigerinoides sacculifer) to 3.6 Ma (Chen, 1994). We then constructed an age model by matching, correlating and tuning the benthic delta18O record to a model simulation of ice volume (Imbrie and Imbrie, 1980, doi:10.1126/science.207.4434.943). The filtered 41- and 23-kyr signals based on the resultant astronomically tuned age model are highly correlated to obliquity (r=0.83) and precession (r=0.75), respectively. Although derived with methodology different from Shackleton et al. (1990) and Hilgen (1991, doi:10.1016/0012-821X(91)90206-W, 1991, doi:10.1016/0012-821X(91)90082-S), our results generally agree with their published astronomical timescales for the time interval from 0 to 3.0 Ma, providing additional support for the newly emerging chronology based on orbital tuning. Slight discrepancies exist in the time interval from 3.0 to 3.6 Ma, suggesting several possibilities, including differences in the approaches of orbital tuning and the relatively low amplitude of delta18O variations in our record. However, even if the discrepancies are due to the relatively low amplitude of the isotope signals in our record at 3.0-3.6 Ma, our resultant timescale as a whole does not adversely affect our evaluation of the paleoclimatology and paleoceanography of the Indian Ocean, such as the evolution of the 100-, 41- and 23-kyr cycles, and variation of global ice volume and deepwater temperature during the past 3.6 m.y.