Is the astronomical forcing a reliable and unique pacemaker for climate? A conceptual model study

There is evidence that ice age cycles are paced by astronomical forcing, suggesting some kind of synchronisation phenomenon. Here, we identify the type of such synchronisation and explore systematically its uniqueness and robustness using a simple paleoclimate model akin to the van der Pol relaxation oscillator and dynamical system theory. As the insolation is quite a complex quasiperiodic signal involving different frequencies, the traditional concepts used to define synchronisation to periodic forcing are no longer applicable. Instead, we explore a different concept of generalised synchronisation in terms of (coexisting) synchronised solutions for the forced system, their basins of attraction and instabilities. We propose a clustering technique to compute the number of synchronised solutions, each of which corresponds to a different paleoclimate history. In this way, we uncover multistable synchronisation (reminiscent of phase- or frequency-locking to individual periodic components of astronomical forcing) at low forcing strength, and monostable or unique synchronisation at stronger forcing. In the multistable regime, different initial conditions may lead to different paleoclimate histories. To study their robustness, we analyse Lyapunov exponents that quantify the rate of convergence towards each synchronised solution (local stability), and basins of attraction that indicate critical levels of external perturbations (global stability). We find that even though synchronised solutions are stable on a long term, there exist short episodes of desynchronisation where nearby climate trajectories diverge temporarily (for about 50 kyr). As the attracting trajectory can sometimes lie close to the boundary of its basin of attraction, a small perturbation could quite easily make climate to jump between different histories, reducing the predictability. Our study brings new insight into paleoclimate dynamics and reveals a possibility for the climate system to wander throughout different climatic histories related to preferential synchronisation regimes on obliquity, precession or combinations of both, all over the history of the Pleistocene.