Filtered state estimates with the SM91 model constrained by the SPECMAP data (squres), and the Antarctic ice core data (pluses). The state estimates are represented by shades, dark and light gray representing the 25 th ; 75 th and 5 th ; 95 th quantiles of the particle weighted distributions, respectively. The lower graph represents, for each data, the model predictive probability that the data would have been lower than it actually was, given the previous parameter and state estimates. The repetition of probabilities below 0.05 or above 0.95 tend to invalidate the model.... Leads and lags between CO 2 and ice volume is another difficult problem where risks posed by hidden dating assumptions and circular reasoning lie at every corner. Here is one typical illustration: Saltzman and Verbitsky showed at several occasions (e.g.: ) a phase diagram showing the SPECMAP δ 18 O stack versus the first full ice-core records of CO 2 available at that time . It is reproduced here (Figure fig:phase, left). The phase diagram clearly suggests that CO 2 leads ice volume at the 100-kyr time scale. However, a detailed inspection of the original publications reveals that the SPECMAP record was astronomically tuned, and that the Vostok time-scale uses a conventional date of isotopic stage 5.4 of 110 kyr BP …by reference to SPECMAP ! The hysteresis is therefore partly conditioned by arbitrary choices. On Figure fig:phase we further illustrate the fact that the shape of the hysteresis depends on the stack record itself. The situation today is that there is no clear consensus about the phase relationship between ice volume and CO 2 at the glacial interglacial time scale (compare ). According to the quite careful analysis of , CO 2 leads ice volume at the precession (20 kyr) period, but CO 2 and ice volume are roughly synchronous at the obliquity (40 kyr) period. Current evidence about the latest termination is that decrease in ice volume and the rise in CO 2 were grossly simultaneous and began around 19,000 years ago... Visual inspection of the LR04 stack (Figure fig:LR04) nicely evidences the gradual transition from the Pliocene — warm and fairly stable — to the spectacular oscillations of the late Pleistocene. The globally averaged temperature at the early Pliocene was about 5 ∘ C higher than today ( and references therein); that one at the last glacial maximum (20,000 years ago) was roughly 5 ∘ C lower. The central research question we are busy with is to characterise these oscillations, understand their origin and qualify their predictability.... In response to this problem Ghil and colleagues promoted single-spectrum analysis, in which a time series is linearly decomposed into a number of prominent modes (which need not be harmonic), plus a number of small-amplitude modes. Assuming that the two groups are indeed separated by an amplitude gap, the first group provides the low-order backbone of the signal dynamics while the second group is interpreted as stochastic noise. Single spectrum analysis was applied with a certain success to various sediment and ice-core records of the few last-glacial interglacial cycles and have in general confirmed that the backbone of climate oscillations may be captured as a linear combination of a small number of amplitude and / or frequency-modulated oscillations. Single-spectrum analysis of the last million years of LR04 (Figure fig:LR04ssa) confirms this statement.... Saltzman and Maasch explored the different solution regimes of this system and they observed that climate trajectories converged to a limit cycle characterised by saw-tooth-shaped oscillations for a realistic range of parameters. When the model is further forced by the astronomical forcing, the uncertainty left on the empirical parameters of equations ( dIdt – dthetadt) provides the freedom to obtain very convincing solutions for the variations in ice volume and CO 2 during the lateQuaternary. Figure fig:Saltzman15 reproduces the original solution , using the parameters published at the time. As in the original publication, the solution is compared with Imbrie’s δ 18 O-stack interpreted as a proxy for ice volume, and CO 2 record extracted from the Vostok and EPICA(Antarctica) ice cores .... Response of the palaeoclimate model of Saltzman and Maasch (1991) . Shown are the insolation forcing, taken as the summer solstice incoming solar radiation at 65 ∘ N after ; the ice volume anomaly (full), overlain with the SPECMAP planctonic δ 18 O c stack (dashed), the CO 2 atmospheric concentration, overlain with the Antarctic ice core data from Vostok and EPICA , and finally deep-ocean temperature. Note that I ' and θ ' are anomalies to the tectonic average. A similar figure was shown in the original article by Saltzman and Maasch... The concentration in CO 2 measured in the Vostok ice core record over the last glacial-interglacial cycle is plotted versus two proxies of continental ice volume: (left) : The planctonic δ 18 O stack by Imbrie et al. (1984) and (right) : The benthic δ 18 O stack by Lisiecki and Raymo (2004). Numbers are dates, expressed in kyr BP (before present). While the Imbrie stack suggests an hysteresis behaviour with CO 2 leading ice-volume variations, the picture based on LR04 is not so obvious.