We describe the evolution of climate system dynamics by examining the climate response to changes in obliquity and precession over the last 5.3 Myr. In particular, we examine changes in the shape of glacial cycles and the power of obliquity and precession response in benthic δ18O. When the exponential trend in δ18O variance is removed, its spectral power exhibits strong, proportional responses to amplitude modulations in orbital forcing over most of the Plio-Pleistocene. Precession responses correlate with modulations in forcing for the last 5 Myr, but 41-kyr response is sensitive to obliquity modulation only before 1.4 Myr. Where responses are sensitive to modulations in forcing, we demonstrate that glacial cycles are orbitally forced rather than being self-sustained or paced by orbital changes. The shapes of glacial cycles have several nonlinear properties, which may be indicative of glacial-interglacial differences in climate sensitivity or response time. The "saw-tooth" asymmetry of glacial cycles first appears shortly after the onset of major northern hemisphere glaciation, and the relative duration of interglacial stages decreases at 1.4 Myr. Collectively, trends in the shape of glacial cycles and the sensitivity of δ18O to obliquity and precession are suggestive of major transitions in climate dynamics at approximately 2.5 and 1.4 Myr but show no significant change associated with the appearance of strong 100-kyr cycles during the mid-Pleistocene transition. © 2006 Elsevier Ltd. All rights reserved.

Benthic δ 18O is often used as a stratigraphic tool to place marine records on a common age model and as a proxy for the timing of ice volume/sea level change. However, Skinner and Shackleton (2005) found that the timing of benthic δ 18O change at the last termination differed by 3900 years between one Atlantic site and one Pacific site. These results suggest that benthic δ 18O change may not always accurately record the timing of deglaciation. We compare benthic δ 18O records from 20 Atlantic sites and 14 Pacific sites to evaluate systematic differences in the timing of terminations in benthic δ 18O. Analysis of sedimentation rates derived from the alignment of benthic δ 18O suggests a statistically significant Atlantic lead over Pacific benthic δ 18O change during the last six terminations. We estimate an average Pacific benthic δ 18O lag of 1600 years for Terminations 1-5, slightly larger than the delay expected from ocean mixing rates given that most glacial meltwater probably enters the North Atlantic. We additionally find evidence of ∼4000-year Pacific δ 18O lags at approximately 128 ka and 330 ka, suggesting that stratigraphic correlation of δ 18O has the potential to generate age model errors of several thousand years during terminations. A simple model demonstrates that these lags can be generated by diachronous temperature changes and do not require slower circulation rates. Most importantly, diachronous benthic δ 18O responses must be taken into account when comparing Atlantic and Pacific benthic δ 18O records or when using benthic δ 18O records as a proxy for the timing of ice volume change. Copyright 2009 by the American Geophysical Union.

We present a 5.3-Myr stack (the "LR04" stack) of benthic δ18O records from 57 globally distributed sites aligned by an automated graphic correlation algorithm. This is the first benthic δ18O 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 δ18O 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 δ18O 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 δ18O from 2.7-1.6 Ma is primarily a deep-water temperature signal and that the phase of δ 18O precession response changed suddenly at 1.6 Ma. Copyright 2005 by the American Geophysical Union.

The sequence of feedbacks that characterized 100-kyr glacial cycles of the past million years remains uncertain, hampering an understanding of the interconnections between insolation, ice sheets, greenhouse gas forcing, and climate. Critical to addressing this issue is an accurate interpretation of the marine δ18O record, the main template for the Ice Ages. This study uses a global compilation of 49 paired sea surface temperature-planktonic δ18O records to extract the mean δ18O of surface ocean seawater over the past 800 kyr, which we interpret to dominantly reflect global ice volume. The results indicate that global surface temperature, inferred deep ocean temperature, and atmospheric CO2 decrease early during each glacial cycle in close association with one another, whereas major ice sheet growth occurs later in glacial cycles. These relationships suggest that ice volume may have exhibited a threshold response to global cooling, and that global deglaciations do not occur until after the growth of large ice sheets. This phase sequence also suggests that the ice sheets had relatively little feedback on global cooling. Simple modeling shows that the rate of ice volume change through time is largely determined by the combined influence of insolation, temperature, and ice sheet size, with possible implications for the evolution of glacial cycles over the past three million years.

Ice loss in the Southern Hemisphere has been greatest over the past 30 years in West Antarctica. The high sensitivity of this region to climate change has motivated geologists to examine marine sedimentary records for evidence of past episodes of West Antarctic Ice Sheet (WAIS) instability. Sediments accumulating in the Scotia Sea are useful to examine for this purpose because they receive iceberg-rafted debris (IBRD) sourced from the Pacific- and Atlantic-facing sectors of West Antarctica. Here we report on the sedimentology and provenance of the oldest of three cm-scale coarse-grained layers recovered from this sea at International Ocean Discovery Program Site U1538. These layers are preserved in opal-rich sediments deposited ∼1.2 Ma during a relatively warm regional climate. Our microCT-based analysis of the layer's in-situ fabric confirms its ice-rafted origin. We further infer that it is the product of an intense but short-lived episode of IBRD deposition. Based on the petrography of its sand fraction and the Phanerozoic ⁴⁰Ar/³⁹Ar ages of hornblende and mica it contains, we conclude that the IBRD it contains was likely sourced from the Weddell Sea and/or Amundsen Sea embayment(s) of West Antarctica. We attribute the high concentrations of IBRD in these layers to "dirty" icebergs calved from the WAIS following its retreat inland from its modern grounding line. These layers also sit at the top of a ∼366-m thick Pliocene and early Pleistocene sequence that is much more dropstone-rich than its overlying sediments. We speculate this fact may reflect that WAIS mass-balance was highly dynamic during the ∼41-kyr (inter)glacial world.