Pliocene-Pleistocene stack of globally distributed benthic stable oxygen isotope records

We present a 5.3-Myr stack (the ''LR04'' stack) of benthic d18O records from 57 globally distributed sites aligned by an automated graphic correlation algorithm. This is the first benthic delta18O stack composed of more than three records to extend beyond 850 ka, and we use its improved signal quality to identify 24 new marine isotope stages in the early Pliocene. We also present a new LR04 age model for the Pliocene-Pleistocene derived from tuning the delta18O stack to a simple ice model based on 21 June insolation at 65 N. Stacked sedimentation rates provide additional age model constraints to prevent overtuning. Despite a conservative tuning strategy, the LR04 benthic stack exhibits significant coherency with insolation in the obliquity band throughout the entire 5.3 Myr and in the precession band for more than half of the record. The LR04 stack contains significantly more variance in benthic delta18O than previously published stacks of the late Pleistocene as the result of higher resolution records, a better alignment technique, and a greater percentage of records from the Atlantic. Finally, the relative phases of the stack's 41- and 23-kyr components suggest that the precession component of delta18O from 2.7-1.6 Ma is primarily a deep-water temperature signal and that the phase of d18O precession response changed suddenly at 1.6 Ma.

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.

Pliocene to Pleistocene calcium carbonate and organic carbon record of ODP Hole 117-722B

Site 722 provides high resolution records of percent CaCO3, magnetic susceptibility, d18O, organic carbon, and coarse fraction for the past 3.4 m.y. from the crest of the Owen Ridge, northwestern Arabian Sea. Within this time interval, most of the carbonate percent variations can be attributed to terrigenous dilution and do not reflect changes in the carbonate system. From the late Pliocene to Present, the average rate of calcium carbonate accumulation increases from 1 to 3 g/cm**2/k.y. and the average accumulation of organic carbon decreases from 75 to 30 mg/cm**2/k.y. The carbonate component is more dissolved in the older interval. The long-term variations in carbonate accumulation may reflect a greater input of organic matter in the late Pliocene, which decomposes to produce CO2 and dissolve carbonate. Magnetic susceptibility and % noncarbonate (100 - CaCO3%) reflect changes in the amount of the lithogenic component in the sediments. The period of variation of lithogenic material is the same period as the original forcing of the regional summer monsoon, however, the timing matches global aridity patterns and global ice volume (sea level) changes. This preliminary analysis suggests that the high frequency variation of lithogenic material persists for at least the last 3.4 m.y. Within the last million years, calcium carbonate accumulation has a large amplitude signal that covaries with major changes in ice volume. Both calcium carbonate and noncarbonate (mostly terrigenous) accumulation are greatest during glacial stages. Interglacial intervals are characterized by low mass accumulation rates, increased foraminifer fragmentation, and increased opal concentration. The accumulation of organic carbon matches the high frequency changes in sedimentation rates. We attribute this high correlation to enhanced preservation of organic carbon by increased sedimentation rate. Of the three major biological components studied, only opal exhibits the variations expected for a biological productivity system forced by monsoonal upwelling driven by changes in northern hemisphere summer radiation.