|Oscillators and relaxation phenomena in Pleistocene climate theory.|
|Jump to Full Text|
|PMID: 22291227 Owner: NLM Status: PubMed-not-MEDLINE|
|Ice sheets appeared in the northern hemisphere around 3 Ma (million years) ago and glacial-interglacial cycles have paced Earth's climate since then. Superimposed on these long glacial cycles comes an intricate pattern of millennial and sub-millennial variability, including Dansgaard-Oeschger and Heinrich events. There are numerous theories about these oscillations. Here, we review a number of them in order to draw a parallel between climatic concepts and dynamical system concepts, including, in particular, the relaxation oscillator, excitability, slow-fast dynamics and homoclinic orbits. Namely, almost all theories of ice ages reviewed here feature a phenomenon of synchronization between internal climate dynamics and astronomical forcing. However, these theories differ in their bifurcation structure and this has an effect on the way the ice age phenomenon could grow 3 Ma ago. All theories on rapid events reviewed here rely on the concept of a limit cycle excited by changes in the surface freshwater balance of the ocean. The article also reviews basic effects of stochastic fluctuations on these models, including the phenomenon of phase dispersion, shortening of the limit cycle and stochastic resonance. It concludes with a more personal statement about the potential for inference with simple stochastic dynamical systems in palaeoclimate science.|
Related Documents :
|12006007 - Probabilistic formalism and hierarchy of models for polydispersed turbulent two-phase f...
20650717 - In situ visualization for large-scale combustion simulations.
20866867 - On the validity of the local diffusive paradigm in turbulent plasma transport.
20922007 - Beam scintillations for ground-to-space propagation. part i: path integrals and analyti...
24104277 - Full-duplex lightwave transport systems based on long-haul smf and optical free-space t...
22418337 - 175 fs tm:lu2o3 laser at 2.07 µm mode-locked using single-walled carbon nanotubes.
|Type: Journal Article|
|Title: Philosophical transactions. Series A, Mathematical, physical, and engineering sciences Volume: 370 ISSN: 1364-503X ISO Abbreviation: Philos Trans A Math Phys Eng Sci Publication Date: 2012 Mar|
|Created Date: 2012-01-31 Completed Date: 2012-03-27 Revised Date: 2013-04-24|
Medline Journal Info:
|Nlm Unique ID: 101133385 Medline TA: Philos Trans A Math Phys Eng Sci Country: England|
|Languages: eng Pagination: 1140-65 Citation Subset: -|
|Georges Lemaître Centre for Earth and Climate Research, Earth and Life Institute, Université Catholique de Louvain, Louvain-la-Neuve, Belgium. firstname.lastname@example.org|
|APA/MLA Format Download EndNote Download BibTex|
Journal ID (nlm-ta): Philos Transact A Math Phys Eng Sci
Journal ID (publisher-id): RSTA
Journal ID (hwp): roypta
Publisher: The Royal Society Publishing
This journal is © 2012 The Royal Society
Print publication date: Day: 13 Month: 3 Year: 2012
pmc-release publication date: Day: 13 Month: 3 Year: 2012
Volume: 370 Issue: 1962
First Page: 1140 Last Page: 1165
PubMed Id: 22291227
Publisher Id: rsta20110315
|Oscillators and relaxation phenomena in Pleistocene climate theory Alternate Title:Review. Oscillators in Pleistocene|
|Georges Lemaître Centre for Earth and Climate Research, Earth and Life Institute, Université catholique de Louvain, Louvain-la-Neuve, Belgium
[other] One contribution of 13 to a Theme Issue ‘Climate predictions: the influence of nonlinearity and randomness’.
The Pliocene and the Pleistocene cover approximately the past 5 Myr. The climatic fluctuations that characterized this period may be reconstructed from numerous natural archives, including marine, continental and ice core records. These archives show a complex climate history. Ice sheets appeared in the northern hemisphere around 3–3.5 Ma (million years) ago [1,2]. The volume of these ice sheets fluctuated with the variations of the seasonal and spatial distributions of incoming solar radiation (insolation), which are induced by changes in the geometry of the Earth's orbit and the angle (obliquity) between Earth's equator and the ecliptic [3–5]. This is called the astronomical forcing.1 Glacial cycles had an average duration of about 40 000 years  until about 800 000 years ago. The dominant period of glacial cycles increased around 800 000 years ago and this is referred to as the Middle Pleistocene Transition. Data and models about the Middle Pleistocene Transition are reviewed by Clark et al. . Time-series analyses based on band-pass filtering provide further evidence of the nonlinear nature of the climate response to the astronomical forcing, from about 1.4 Ma ago . The latest four glacial cycles, in particular, are distinguished by a pronounced saw-tooth time structure: ice accumulates over the continents during about 80 000 years and then melts in about 10 000 years (figure 1).
Superimposed on these long glacial cycles comes a complex pattern of millennial and sub-millennial variability . For example, the Greenland record features at least 20 events of abrupt rise and slower decline in oxygen isotopic ratio (a proxy for temperature) [12,13] and methane  during the latest glacial epoch. These events are known as Dansgaard–Oeschger events. They were found to occur from at least the last glacial inception  and the Antarctic ice core records provide evidence that they are characteristic of Pleistocene glacial climates . Some of these events follow pulses of iceberg discharges into the North Atlantic Ocean, called Heinrich events [17–19]. Heinrich events and Dansgaard–Oeschger events have left climatic footprints all over the globe , including in Antarctica . The current interglacial period is referred to as the Holocene. It is also characterized by millennial and centennial variability, mainly observed in the North Atlantic [21–23], but of a much weaker amplitude than during the preceding glacial period.
This paper reviews attempts to explain these fluctuations with concepts that originate in dynamical system theory. These are the concepts of limit cycle, synchronization and excitability. The central message of the paper is that current theories of ice ages and rapid events may often be interpreted in terms of generic deterministic models, which are also used in other areas of science such as biology and ecology. However, stochastic parametrizations are an essential part of any complex system model, and their effects on climatic oscillations have to be taken into account.
Dynamical system theory entered palaeoclimate science with idealized models representing the response of ice sheets to the astronomical forcing. These models were directly derived from the physics of the ice sheet–atmosphere system [24–27]. Ghil & Childress , in particular, insisted on the interest of analysing such models in terms of bifurcation theory. For modelling the complex carbon cycle response, authors sometimes adopted a more heuristic approach by considering simple models and confronting the results of palaeoclimate evidence .
Nowadays, climate research is largely oriented towards large climate simulators (typically general circulation models), which are developed to include as many climate processes as possible. However, thinking in terms of dynamical system theory remains insightful. Indeed, the behaviour of a complex system at a certain spatio-temporal scale is in practice often dominated by a few leading modes, of which the dynamics may be captured fairly convincingly with a low-order dynamical system. Climate scientists are increasingly using this property. For example, they formulate simpler models to explain the seemingly complex behaviours observed in ocean–atmosphere simulators. Examples have been provided in recent years by focusing on interannual , centennial  and millennial [32,33] variability. In parallel, so-called hysteresis experiments, which aim at identifying the number of stable states in individual components of the climate system such as the ocean circulation  or ice sheets , contribute to a dynamical system-founded understanding of the climate system. This approach may also help us to predict and communicate about the proximity of bifurcations, which may result in catastrophic climatic changes. Timmermann & Jin  termed our ability to anticipate bifurcation phenomena as predictability of the third kind, by reference to the predictabilities of the first and second kind originally introduced by Lorenz .
The article is structured as follows. Section 2 reviews some of the basic concepts of oscillator theory. This is no substitute for proper textbook reading, but the reader will find essential notions and definitions needed to understand the remainder. Section 3 reviews how these concepts enter theories of ice ages and rapid events. Section 4 discusses effects of stochastic fluctuations and, finally, §5 is a more personal statement about the potential for inference with simple stochastic dynamical systems in palaeoclimate science.
The reader will find an accessible introduction to dynamical system theory and concepts in Strogatz . More formal background on oscillator theory, albeit a bit dated, is available in Guckenheimer & Holmes . Bifurcation and oscillator theory is explicitly connected to climate theory in Ghil & Childress , in particular, ch. 12 and Saltzman , ch. 7. Background on synchronization and an introduction to the phenomenon of excitability is available in Pikovski et al. . Finally, the Scholarpaedia peer-reviewed website is an increasingly rich and authoritative source of information on dynamical systems. Only the notions essential for the present article are summarized here.
Oscillator. The oscillator is a dynamical system that has a globally attracting limit cycle. In more simple terms, it oscillates even in the absence of an external drive. Here, we are interested in oscillators to describe climate phenomena, which involve dissipation of energy. The minimal model for a dissipative oscillator includes two ordinary differential equations, of which at least one is nonlinear.
Relaxation oscillator. The relaxation oscillator is a particular kind of oscillator featuring an interplay between relaxation dynamics (generally fast) and a destabilization process (generally slow). The relaxation is the process by which the system is attracted to a region of the phase space. This evokes the relaxation of a spring. In a relaxation oscillator, the system continues to evolve slowly after the relaxation phase. During this slow evolution phase, the system's stability diminishes gradually until the system is ejected out of its relaxation state, either towards another relaxation state or to the same relaxation state via a dissipative loop. In this review, we will encounter three kinds of relaxation oscillators (figure 2): relaxation founded on slow–fast dynamics (involving a slow manifold); relaxations structured by a homoclinic orbit (involving only one relaxation state); and relaxations structured around a focus. More details are provided in the caption of figure 2.
Excitability. An excitable system has a globally attracting fixed point (it does not oscillate spontaneously). However, an external perturbation may have the effect of exciting it. During this excitation, the system is being ejected far from its fixed point and then returns to it.
Link between relaxation dynamics and excitability. In practice, it is often found that a relaxation oscillator may be transformed into an excitable system by a mere change in parameter, and vice versa. The reason is the following. A relaxation oscillation is often structured globally in the phase space; for example, by a slow manifold (figure 2a) or by one or several saddle points (figure 2b). Suppose now that the oscillation displayed by such a system ceases because a parameter has been changed. The system is then no longer an oscillator, but the ‘backbone’ of the oscillation dynamics are still latent in the phase space because the elements that structured the limit cycle (the slow manifold or the saddle points) have not disappeared. Consequently, the system may be run on a trajectory close to the defunct limit cycle if it is being pushed by some external force (the excitation) into the region of the phase space previously occupied by this limit cycle. This point is illustrated on the basis of slow–fast dynamics in figure 2d, but similar excitation dynamics generally occur near any kind of ‘explosive bifurcation’, that is, bifurcations that give rise rapidly to a fully developed limit cycle. This includes homoclinic, heteroclinic and certain Hopf bifurcations (two examples follow and are illustrated in figure 6).
Saltzman established a theory in which ice ages are interpreted as a limit cycle synchronized on the astronomical forcing. Saltzman et al.  wrote a series of articles on the subject, starting with the introduction of the limit cycle idea  and synchronization hypothesis, the interpretation of the Middle Pleistocene Transition as a bifurcation , and the more complete models in the mid-1990s . The full theory is developed in a book . Here, we concentrate on two intermediate models [46,47]. They are called SM90 and SM91, by reference to the authors (Saltzman & Maasch) and the year of publication. The variables I, μ and θ are the continental ice mass, CO2 concentration and deep-ocean temperature, respectively. The reader is referred to the original publications for the meaning and value of the different parameters. They are not crucial here; it suffices to know that they are all positive.
In both models, the first equation describes the ice mass response to changes in CO2 (μ) and the astronomical forcing (FI(t)); Saltzman adopts the so-called Milankovitch view2 that an increase in insolation causes a decrease in ice mass. Increases in CO2 or in ocean temperature have the same effects.
The other two equations describe the dynamics of CO2 and the response of deep-ocean temperature to changes in ice volume. It is further assumed that the mean state of the climate varied slowly throughout the Pliocene–Pleistocene; in particular, in response to a ‘tectonically driven’ decline in the average concentration in CO2, consistently with an earlier proposal . This tectonically driven decline is modelled here as a slow decrease in the forcing term Fμ(t) throughout the Pleistocene.
Consider the bifurcation diagrams of the SM90 and SM91 models with respect to Fμ, assuming no astronomical forcing (FI=0) (figure 3). The systems are then said to be free or autonomous. Depending on Fμ, both models show regimes with a stable fixed point, and regimes for which the fixed point is unstable, so that the system orbits along a limit cycle.
These considerations led Saltzman to interpret the Middle Pleistocene Transition as a bifurcation between a ‘quasi-linear’ response regime to the astronomical forcing (in the fixed-point regime) to a regime of nonlinear synchronization (resonance) on the astronomical forcing. He concluded that ice ages would occur today even in the absence of astronomical forcing. The main effect of the astronomical forcing is to control the timing of glaciations.
Saltzman's theory is seductive because it explains in a consistent framework several aspects of the Pleistocene climate history, including the change from linear to nonlinear regime , the presence of 100 000 year periodicity in climate records , the lack of a 400 000 year spectral peak in the ice volume record (such a peak appears in the simple piece-wise linear model devised by Imbrie & Imbrie , owing to rectification of the precession signal), the synchronization of deglaciations on the astronomical forcing [53,54], and the occurrence of large climatic transitions even when eccentricity, which modulates the effect of precession on insolation, is at its lowest.
The difficulty for accepting Saltzman's models as a definitive theory lies in the physical interpretation of the CO2 equation. This equation encapsulates all the interesting dynamics of the system and it is thus crucial to the theory. Some semi-empirical justification for the CO2 equation is given in Saltzman & Sutera , but the form of this equation has undergone somewhat ad hoc adjustments in SM90. The form present in SM91 is again different, with important effects on the bifurcation structure, while the authors did not justify this latter change based on physical or biogeochemical considerations.
To better appreciate the structural differences between the two models, let us return to the bifurcation diagrams. Consider SM90. As the forcing is decreased, the fixed point gives rise to a locally unstable limit cycle. The system must, therefore, find a stable limit cycle further away from the fixed point, but in this case not much further. This stable limit cycle is under the influence of the unstable fixed point and, in particular, the system slows down when it passes near it (figure 4). This scenario is depicted in figure 2c. The limit cycle evolves as the tectonic forcing is further decreased, until it shrinks smoothly around a perpetually glaciated state.
In SM91, the bifurcation induced by the decrease in tectonic forcing is much more explosive. The system lands on a stable limit cycle that turns out to be little affected by the position of the unstable point. The cycle dynamics do not show clear phases of acceleration and the system cannot be regarded as a relaxation system. The limit cycle disappears abruptly as the tectonic forcing is further decreased, through a phenomenon called a saddle bifurcation of cycles. The consequences of the difference between the bifurcation structures of SM90 and SM91 may be further appreciated in the transient experiments shown in figure 5.
Paillard & Labeyrie  have been advocating the concept of relaxation for understanding palaeoclimate dynamics, both ice ages and the more abrupt events, since the publication of a seminal paper in 1994. We return to this article later on, and concentrate on another article published in 1998, in which Paillard  introduces a conceptual model of ice ages. Ice volume dynamics respond to an ordinary differential equation,57] because the relaxation equation involves a discrete state variable, here denoted by y. Its state may be ‘deep glacial’ (G), ‘mild glacial’ (g) or ‘interglacial’ (i). The numerical values of xR and τR depend on this climate state. Climate states y follow a sequence according to a set of conditions formulated on the level of glaciation x and insolation. Namely, the transition is triggered when the forcing F(t) exceeds a certain threshold. Occurrence of G drives climate quickly into an interglacial state i because xR(G) and τR(G) are specified in the model to be low.
Paillard is not very specific about the physical meaning of the discrete variable, but it accommodates the paradigm that the Atlantic Ocean circulation has gone through three different states during the latest glacial period: intermediate circulation, shut down of the circulation and modern, deep-sinking circulation. The system (P98) features the concept of slow–fast relaxation dynamics. However, this is not an oscillator because the shift from g to G is determined by the course of the external forcing. The Middle Pleistocene Transition is induced in (P98) in a fashion similar to that in Saltzman, and on the basis of similar physical assumptions (tectonically driven decline in CO2). The drift in climatic conditions induced by tectonics is accounted for by a term added to the astronomical forcing. In a later review, Paillard  further emphasizes empirical evidence for the relevance of the relaxation concept in the phenomenon of deglaciation.
Gildor & Tziperman  take a moderate step towards higher model complexity by considering a slightly more explicit representation of atmosphere, ocean, sea ice and land ice dynamics. Namely, the ocean is divided into eight boxes, and the atmosphere into four boxes. The sea ice fraction responds to standard energy balance equations. More crucially, land ice growth is influenced by a somewhat controversial feedback between sea ice and precipitation. The feedback is controversial because it is assumed that cold climate results in a reduction in ice volume: sea ice growth causes a reduction in precipitation in ice-covered areas and, by this mechanism, almost suppresses accumulation of snow on ice sheets. The latter then no longer compensates for ice ablation and ice volume shrinks.
A free oscillation arises from the fact that the ice volume thresholds for switching sea ice cover ‘on’ and ‘off’ differ. In other words, sea ice displays a hysteresis response to variations in ice volume. This is exactly the principle of the slow–fast relaxation oscillator depicted in figure 2a: the curve of equilibrium of sea ice with respect to ice volume is the slow manifold, and ice volume integrates the state of sea ice in time. In turn, this oscillation can be synchronized on the astronomical forcing.
The Gildor–Tziperman model is coupled to a biogeochemical cycle in a companion paper , but the essential dynamics of the glacial oscillation are unchanged. Tziperman et al.  further comment on the model and its property of synchronization on the astronomical forcing, and find that its behaviour is essentially reducible to a hybrid dynamical system.
Paillard & Parrenin  proposed yet another relaxation model in 2004 (PP04). The prognostic variables are ice volume I, the area of the Antarctic continental ice sheet A and the atmospheric concentration in CO2 (μ) (a…j are parameters):
As in the other ice age models, ice volume is a slow variable driven by the astronomical forcing. It is here coupled to a variable with a similar time scale (τA∼τI) and a faster one (τμ=τI/3). The term H(−D), where H is the Heaviside function, represents the ventilation of the Southern Ocean. CO2 is released into the atmosphere when the Southern Ocean is ventilated (D<0), which drives deglaciation. Ice then grows slowly, until a Southern Ocean ventilation flush sends the system back to interglacial conditions. Ocean ventilation is thus the fast process in this model and it is the only nonlinear process accounted for. However, contrary to the Gildor–Tziperman model, it does not present a hysteresis behaviour. Consequently, the glacial cycles featured by this model cannot be interpreted in terms of shifts between the branches of a slow manifold.
To better understand the dynamics of glacial cycles in this model, we consider the bifurcation diagram along typical solutions in the phase space for the free (i.e. unforced) system (figure 6). The parameter g is taken in this example as the control parameter, in order to preserve Saltzman's idea that ice age cycles appear as the consequence of a slow perturbation of the carbon cycle. As in SM90, PP04 exhibits a limit cycle arising from a subcritical Hopf bifurcation. The dynamics along the limit cycle close to the bifurcation point are strongly influenced by the presence of the unstable focus. This is the configuration shown in figure 2c. Depending on g, the focus is either on the low-ice-volume side of the limit cycle (i.e. the system spends most of its time with high CO2) or on the high-volume side of the limit cycle (i.e. the system spends most of its time in low CO2). Parrenin and Paillard estimate that we are currently in the second configuration.
It has been claimed [61,63] that any model that has some form of 100 000 year internal periodicity could be used to reproduce the course of ice volume over the last 800 000 years. Taking the argument at face value, Crucifix  used one of the simplest possible slow–fast oscillators—the van der Pol (VDP) oscillator—with minimal modifications to account for the astronomical forcing and the asymmetry between the phase of ice build-up and melt during the Late Pleistocene (α, β, γ and τ are parameters; F(t) is the astronomical forcing),
The system dynamics are determined by the structure of the slow manifold x=y3/3−y. The parameter β controls the position of the fixed point on the slow manifold and, consequently, the ratio of times spent by the system in the two branches (‘glacial’ and ‘interglacial’) of the slow manifold. The ice age curve can be captured with some tuning (figure 7), although it is fair to add that a small change in parameters may shift the timing of one or several ice age cycles. This minimal model was used to challenge intuitive arguments about the predictability of ice ages .
Welander  introduced the concept of relaxation oscillations in the context of ocean dynamics. He described a heat–salt oscillator involving exchanges of heat and salt within a single oceanic column, coupled to a phenomenon of surface temperature relaxation. The destabilization process needed for the relaxation oscillation to appear is here related to diffusion between the deep ocean and the mixed layer. The system dynamics are further controlled by the mean freshwater flux at the top of the ocean column. It determines the transitions from a regime characterized by perpetual convection in the oceanic column, to a regime with intermittent convection (oscillation) and finally to a regime with no convection . The bifurcations between the different regimes bear the character of global bifurcations, with the oscillation period approaching infinity near the bifurcation points (in particular, the second bifurcation bears the character of a homoclinic bifurcation). The heat–salt oscillator belongs thus to the class shown in figure 2b.
Welander  and Winton & Sarachik  later introduce the concept of another kind of relaxation oscillator in the ocean. It involves the meridional structure of the ocean thermohaline circulation, and the key nonlinear process is the meridional advection of heat and salt. The oscillations featured by this model are termed ‘deep-decoupling’ oscillations . Given that the slow process now relates to heat accumulating in the global ocean, the characteristic time of deep-decoupling oscillations is of the order of 1000 years. The net flux of freshwater delivered to the North Atlantic acts as a bifurcation parameter controlling the transition between non-oscillating and oscillating regimes in the Winton–Sarachik model .
Millennial oscillations have since been observed across a hierarchy of ocean models, including three-dimensional ocean models with prescribed freshwater flux and restoring conditions to surface temperature , and three-dimensional models coupled to a simple atmosphere [71,72]. Sakai & Peltier  proposed that millennial deep-decoupling oscillations could explain Dansgaard–Oeschger events. Colin de Verdière et al. [33,74,75] complement this early proposal with a fairly complete theory based on ocean circulation model experiments. The oscillations described by Colin de Verdière et al. involve the processes of turbulent vertical mixing (neglected in Winton & Sarachik ), advection and convection, which unify the salt oscillator with the deep-decoupling oscillation model. Incidentally, Colin de Verdière et al.  dismiss the nonlinearity of the equation of state as the cause of the oscillations.
There is, across the model hierarchy, consistency about the fact that the transition between the oscillating circulation regime and the so-called diffusive, haline regime (without deep convection) is associated with a homoclinic bifurcation [33,76]. The nature of the bifurcation between the convective regime and the oscillation is more model-dependent. Timmermann et al. , based on experiments with the eight-ocean box Gildor–Tziperman model, find a Hopf bifurcation; salt-conserving experiments with a two-dimensional ocean model show a transition towards a finite-period cycle, but of increasing period as the bifurcation is approached; experiments with a more idealized model, formulated as a two-equation dynamical system, reveal the signature of an infinite-period bifurcation .3 The latter implies that Dansgaard–Oeschger events, at the time when they appear soon after the glacial inception process, should be very long but of a similar amplitude to the Dansgaard–Oeschger events coming later in the glacial cycle. This feature is consistent with the Greenland ice core record (figure 1). More specifically, the first Dansgaard–Oeschger cycles that appeared at the beginning of the glacial era were characterized by a long ‘plateau’ phase (also called interstadial) during which the thermohaline circulation was certainly very active . In the Colin de Verdière et al.  theory, the plateau phase is the phase of the trajectory influenced by the ‘ghost of the saddle point’.
Given the explosive nature of the bifurcations involved in ocean dynamics it is no surprise to find excitability properties in ocean models. Weaver & Hughes  discuss this effect in salt-conserving experiments with an idealized-geometry, ocean model. The ocean–atmosphere model of intermediate complexity, CLIMBER (CLIMate BiosphERe model), was shown to exhibit excitability properties when boundary conditions are set to be typical of the latest glacial era . The ocean circulation has then one stable state, with a moderate Atlantic overturning and a ‘quasi-stable state’ with more intense overturning. The conceptual sketch of the excitation cycles shown by Ganopolski & Rahmstorf , fig. 1 can be interpreted in terms of slow–fast dynamics, in which the different states of the ocean circulation constitute the different branches of a slow manifold. The intense overturning state, which is the ‘plateau’ phase of the Dansgaard–Oeschger event, may thus be viewed as the repelling branch of the slow manifold in the excitable regime (figure 2d). The excitable Dansgaard–Oeschger hypothesis was used as a possible basis to explain how a weak forcing, exogenous to the system, could explain the observed 1500 year periodicity of Dansgaard–Oeschger cycles (on this periodicity refer to earlier studies [79,80], but see the other view in the study of Ditlevsen & Ditlevsen ). Two such theories were developed on the basis of experiments with CLIMBER. One suggests that Dansgaard–Oeschger events are excited by stochastic fluctuations, modulated by a weak, hypothetical solar periodic forcing  (more on the effects of stochastic fluctuations in §4). The alternative theory suggests that the excitation is induced by the interference between two solar forcings with periods close to 1470/7(=210) and 1470/17(≈87) years , possibly combined with noise .
MacAyeal  proposed an ice binge–purge theory to explain Heinrich events. The theory rests on experiments with a one-spatial direction model of ice flow dynamics. Suppose, as a starting point, that ice volume grows in response to net accumulation of snow. The growth continues until the accumulated effect of geothermal heat flux causes basal sliding. A volume of ice is then released into the ocean (this is the ‘purge’), causing the release of icebergs characteristic of Heinrich events. Ice volume thus decreases, until ice accumulation wins over so that ice volume can grow again. The ice binge–purge model is thus a relaxation oscillator combining a slow integrating process (ice mass accumulation) with a fast lateral discharge process.
To what extent may Heinrich events interfere with Dansgaard–Oeschger dynamics? Paillard [55,86] investigated this question by coupling the MacAyeal ice model—but reduced to ordinary differential equations by Galerkin truncation—with a three-ocean box model. The coupling simply assumes that ice released into the ocean causes a net freshening of the surface of the North Atlantic that alters the deep-ocean circulation. Paillard realized that this coupling could lead to fairly non-intuitive and complex effects, such as the succession of Dansgaard–Oeschger events of decreasing amplitude between Heinrich events. This succession is known in the literature on palaeoclimate records as Bond cycles . Paillard also found that the oscillations are aperiodic in this model under certain parameter configurations.
The issue is further explored in Schulz et al. , based on the Winton–Sarachik ocean model, and in Timmermann et al. , based on the slightly more sophisticated Gildor–Tziperman  ocean model. The objective was to study the response of deep-decoupling ocean oscillations to prescribed Heinrich cycles. Schulz et al.  noted that deep-decoupling oscillations could be synchronized on the Heinrich cycles. Timmermann et al.  then proposed, on the basis of numerical experiments with a fairly idealized model, that Heinrich events excite Dansgaard–Oeschger cycles because the variation in ice volume caused by a Heinrich event modifies slowly the amount of net freshwater released in the ocean. In turn, they suggested, Dansgaard–Oeschger may have a control on ice volume growth. This yields a two-way coupling between Dansgaard–Oeschger and Heinrich events.
Experiments with more comprehensive models of the ocean–ice sheet–atmosphere system [87,88] generally support the idea that the different water and heat fluxes involved in the different phases of ice build-up and iceberg release are quantitatively sufficient to support a coupling between ice sheets and ocean circulation during the latest glacial era. However, it was also noted that ‘three-dimensional thermomechanical ice-sheet models are unable to satisfactorily reproduce the binge–purge mechanism without an ad hoc basal parameterisation’ .
To address this difficulty, a theory in which Dansgaard–Oeschger events trigger Heinrich events was recently proposed . The ice shelf plays a key role, in blocking the ice stream flow from the ice sheet to the oceans. Heinrich events occur when this ice shelf is broken; for example, under the influence of ocean sub-surface warming associated with a Dansgaard–Oeschger event. The resulting model is a system displaying a slow ice build-up—the Heinrich release cycle excited by fluctuations in ocean sub-surface temperature.
The much smaller ocean oscillations that characterized the Holocene period may also be a relaxation phenomenon. Schulz et al.  observed oscillations in the atmosphere–ocean model of intermediate complexity ECBILT-CLIO. These oscillations are related to the convective activity in the Labrador Sea.
Schulz et al.  considered the existence of such an oscillator in an earlier reference and speculated on the possible interactions between the centennial oscillations, millennial oscillations and Heinrich cycles. They considered a model in which each of these three kinds of oscillations is modelled as a Morris–Lecar relaxation oscillator. Their working hypothesis is that glacial conditions induce a coupling between these oscillators. They then observed that a very stable 1500 year oscillation appears, which they interpreted as a model equivalent of Dansgaard–Oeschger events.
The myriad of chaotic motions that characterize the dynamics of the ocean and the atmosphere may be taken into account in the form of parametrizations involving stochastic time processes. The method was introduced in climatology in the 1970s  and the theoretical justifications, which allow one to model chaotic motions as a (linear) stochastic process, are reviewed by Penland [92,93]. In a statistical inferential framework, the stochastic parametrizations may also be viewed as a way to account for the distance necessarily existing between the concepts and dynamics featured by the model, and the complex system being observed.
The effects of stochastic processes on relaxation oscillators and excitable systems are generally well documented in the literature because this is a topic of general interest . Here, we review some of them in the specific context of palaeoclimate dynamics.
One of the basic effects of noise on oscillators is the phenomenon of phase dispersion: a weak stochastic forcing on an oscillator causes a fading out of the memory of the exact initial conditions, even though the gross structure of the oscillation visualized in the phase space is conserved. The phenomenon is well known and it is an immediate consequence of the neutral stability of the phase of a free oscillator with respect to fluctuations. It was previously suggested that this phenomenon of phase dispersion may concern ice ages , but it is more commonly believed that ice ages are phase locked on the astronomical forcing. This phase locking should act against dispersion and permit a very long predictability horizon of ice ages. However, a phenomenon of phase dispersion may happen in oscillators that are locked on a periodic forcing. A stochastic fluctuation may momentarily cause a burst of desynchronization, called phase slip, during which the system is unhooked from its corresponding deterministic trajectory and attracted to another trajectory, which leads or lags the original one by one forcing period , §3.1.3. The difference between phase diffusion in a free system and in a periodic forcing-driven oscillator is that the diffusion effect has, in the latter, a quantum nature. More formally, it is said that the stochastic forcing disperses the system states around the different attractors that are compatible with the forcing. In a work in preparation, we suggest that the astronomically forced climate system may satisfy the conditions for a similar phenomenon of phase dispersion to occur (B. De Saedeleer, M. Crucifix & S. Wieczorek 2011, unpublished data). Given that the astronomical forcing is aperiodic the description of the phenomenon requires a suitable theoretical framework, which relies on the notion of a ‘local pullback attractor’. The equivalent of a phase slip is, in the aperiodic forcing context, a stochastic shift from one of the deterministic pullback attractors to another one. The phenomenon is illustrated based on experiments with the VDP model in figure 8.
Additive fluctuations generally reduce the period of relaxation oscillators. In an oscillator presenting a homoclinic orbit such as the Duffing oscillator, additive fluctuations reduce the time spent near the unstable focus . This implies that, even at the corresponding bifurcation point in the deterministic system, the return time of oscillations in the stochastically perturbed system remains finite. In a slow–fast oscillator such as the VDP oscillator, additive fluctuations generally result in early escapes of the branch of the slow manifold on which the system lies (figure 9). The period of the oscillator is thus affected by a correction that increases approximately linearly with the noise variance in the slow–fast VDP oscillator .4 This property was used at least once in Pleistocene theory, in the silicate weathering hypothesis advanced by Toggweiler . Additive fluctuations reduce the limit cycle period from 800 000 to about 100 000 years. The reasons for the period reduction being so dramatic are left for another study.
Noise may naturally act as an excitation agent in an excitable system (figure 9). The topic is extensively reviewed by Lindner et al. . Excitation loops are sporadic if the noise amplitude is weak, in which case the recurrence time of the events is set by the noise amplitude, whereas the amplitude of the events is set by the structure of the deterministic vector field. The frequency of excitation loops increases with the noise amplitude, until the system behaviour is qualitatively similar to a limit cycle regime. This is the coherent resonance regime. For yet higher noise amplitudes, the limit cycle structure is destroyed.
The concept of stochastic excitation has been considered several times in Dansgaard–Oeschger theories. The idea is introduced based on experiments with an ocean general circulation model with idealized geometry and forcing . The effect of stochastic fluctuations is not only to excite oscillations in a system normally at rest, but also to reduce the period of these oscillations when the system is in oscillatory regime.
A phenomenon of non-autonomous stochastic resonance may occur if the noise is superimposed on a weak external drive. For this to happen, the autonomous system needs to be stable but excitable. The external forcing must be too weak to cause excitation by itself. The role of the noise is to provide the additional power to induce excitation. The timing of the excitation is then related to the phase of the external forcing. The mechanism was proposed several times [76,82] to explain the 1500 year recurrence time of Dansgaard–Oeschger events. The idea remains questioned, either on the grounds that the 1500 year recurrence time observed in palaeoclimate records is coincidental  or on the grounds that the 1500 year external forcing is unidentified [33,83]. A more subtle case of stochastic resonance involves the combination of noise with two solar cycles of 210 and 87 years, which yields the concept of ‘ghost resonance’  for which some support, albeit not conclusive, is found in the observations .
Coupled oscillators may exhibit, collectively, a resonance period that is more robust to external fluctuations than the uncoupled oscillators. Schulz et al.  used this property to explain the stability of the Dansgaard–Oeschger recurrence period of 1500 years in the presence of random fluctuations, without having to invoke an external forcing.
Finally, a behaviour reminiscent of oscillations may occur in a system that is neither oscillating nor excitable, but which presents several stable states. Noise then simply induces jumps between these different states. The simplest mathematical model is the Langevin equation and this is on this basis that Schulz et al.  interpret the Holocene oscillations observed in the ECBILT-CLIO climate model.
The review has shown that relaxation oscillations are a popular and powerful model to explain oscillations observed in the Pleistocene record. The concept of relaxation implies some form of slow–fast separation, in the sense that at least one component of the system spends most of its time in ‘quasi-equilibrium’ states (this may be a ‘slow manifold branch’ or a region influenced by a saddle point, depending on the system structure), with acceleration phases.
Some of these models were constructed following a fairly careful procedure of truncation of a system of partial differential equations, which describes some of the fluid dynamics of the climate system. Others were proposed on a more conceptual basis, the idea being precisely to test a hypothesis based on palaeoclimate observations. The latter approach is sometimes criticized, on the grounds that box models, for example, cannot reasonably be taken as an adequate representation of the complex dynamics of the oceans .
This leads us to the last question of this review: can dynamical systems be used for inference on palaeoclimates? Inference implies that something is being learned by confronting a model with observations. This inference process may take the form of a calibration procedure (update our knowledge on parameters on the basis of observations) or a model selection procedure (which model, among different alternatives, explains the observations best).
The position taken here is that there is not such a thing as an ‘attractor’ of the climate system that is to be ‘discovered’. The hope is that some of its modes of behaviour are sufficiently decoupled from the rest of the variability to justify the fact that simple dynamical systems may capture the fundamental dynamical properties of these modes, and we want to learn about these modes from palaeoclimate observations.
The programme is challenging. Indeed, it was emphasized that different physical assumptions may lead to dynamical systems with dynamical properties that are similar enough to produce a convincing visual fit on palaeoclimate data . The message is largely echoed in this review. The modeller's challenge is, therefore, to operate a model selection on more stringent criteria than just fitting some standard time series. For example, palaeoclimate observations may yield constraints on the bifurcation structure of the system. The Middle Pleistocene Transition is an attractive test case in this respect.
In a statistical inference process, the observations should be a plausible outcome or realization of the model. This makes sense only if the model has a stochastic component, which describes its uncertainties, limitations, and the noise that emerges from the chaotic motions of the atmosphere and oceans.
Stochastic dynamical systems are beginning to be used for inference on palaeoclimate time series. In a method called ‘potential analysis’, the climate system is modelled as a Langevin equation, that is, the combination of a down-gradient drift with a Wiener additive process and inference is made on the number of wells of the potential function . The method was applied on Pleistocene climate records, yielding the conclusion that the number of wells increased from two to three over the course of the Pleistocene .
However, our position so far has been to favour a Bayesian methodology, because it allows one to encode physical constraints in the form of prior distributions on model parameters. The Bayesian formalism is also naturally designed for model calibration, selection and probabilistic predictions.
The fact is that Bayesian methods for selection and calibration of dynamical systems on noisy observations are only emerging. In a recent attempt, we considered a particle filter for parameter and state estimation . To be honest, there is ample room for progress. Whether the process of inference with simple dynamical systems on palaeoclimate data will lead new insight in this context still needs to be demonstrated.
1The astronomical forcing will generally be taken into account here in the form of a normalized measure of insolation during the month or on the day of summer solstice at a northerly latitude, typically 60 or 65° N. This is a fairly complex, aperiodic signal, with dominant harmonics corresponding to the phenomena of precession (23 716, 22 428 and 18 976 years), and obliquity (41 000 years) .
3This particular case is not illustrated in figure 2. There is no saddle point along or near the orbit, but there is a combination of parameters for which a fixed point appears on the limit cycle. Beyond this particular parameter value, this fixed point splits into a saddle and a node. This particular parameter value corresponds to an ‘infinite-period’ bifurcation. In practice, as long as the limit cycle exists, the trajectory slows down near the point where the saddle–node will appear. Some authors then refer to the influence of the ‘ghost’ of the saddle point .
4This result is established assuming fluctuations added to the slow variable. A much more general theory, suitable for different slow manifolds and additive fluctuations to slow and fast variables, is now available .
This review benefited from stimulating discussions during the Isaac Newton Programme Mathematical and Statistical Approaches to Climate Modelling and Prediction held in Cambridge during the summer–autumn 2010. Guillaume Lenoir, Bernard De Saedeleer and Jonathan Rougier provided helpful comments on a first version of the article. Thanks are also due to Didier Paillard, Olivier Arzel, Alain Colin de Verdière, André Paul and Sebastian Wieczorek for e-mail correspondence during this review. The author is Research Associate with the Belgian National Fund of Scientific Research. This research is supported by the FP7-ERC starting grant ITOP ERG-StG-2009-239604. The editor (Jan Sieber) and two reviewers are acknowledged.
|1.||Shackleton N. J.,et al. Year: 1984Oxygen isotope calibration of the onset of ice-rafting and history of glaciation in the north Atlantic regionNature30762062310.1038/307620a0 (doi:10.1038/307620a0)|
|2.||Meyers S. R.,Hinnov L. A.. Year: 2010Northern Hemisphere glaciation and the evolution of Plio-Pleistocene climate noisePaleoceanography2511110.1029/2009PA001834 (doi:10.1029/2009PA001834)|
|3.||Croll J.. Year: 1875Climate and time in their geological relations: a theory of secular changes of the Earth's climateNew York, NYAppleton|
|4.||Milankovitch M.. Year: 1998Canon of insolation and the ice-age problemBeograd, SerbiaNarodna Biblioteka Srbije .(English translation of the original 1941 publication.).|
|5.||Berger A. L.. Year: 1978Long-term variations of daily insolation and quaternary climatic changesJ. Atmos. Sci.352362236710.1175/1520-0469(1978)035%3C2362:LTVODI%3E2.0.CO;2 (doi:10.1175/1520-0469(1978)035<2362:LTVODI>2.0.CO;2)|
|6.||Ruddiman W. F.,Raymo M.,McIntyre A.. Year: 1986Matuyama 41 000-year cycles: North Atlantic Ocean and northern hemisphere ice sheetsEarth Planet. Sci. Lett.8011712910.1016/0012-821X(86)90024-5 (doi:10.1016/0012-821X(86)90024-5)|
|7.||Clark P. U.,Archer D.,Pollard D.,Blum J. D.,Rial J. A.,Brovkin V.,Mix A. C.,Pisias N. G.,Roy M.. Year: 2006The Middle Pleistocene transition: characteristics. mechanisms, and implications for long-term changes in atmospheric pCO2Quat. Sci. Rev.253150318410.1016/j.quascirev.2006.07.008 (doi:10.1016/j.quascirev.2006.07.008)|
|8.||Lisiecki L. E.,Raymo M. E.. Year: 2007Plio-Pleistocene climate evolution: trends and transitions in glacial cycles dynamicsQuat. Sci. Rev.26566910.1016/j.quascirev.2006.09.005 (doi:10.1016/j.quascirev.2006.09.005)|
|9.||Lisiecki L. E.,Raymo M. E.. Year: 2005A Pliocene–Pleistocene stack of 57 globally distributed benthic δ18O recordsPaleoceanography201710.1029/2004PA001071 (doi:10.1029/2004PA001071)|
|10.||North Greenland Ice Core Project MembersYear: 2004High-resolution record of northern hemisphere climate extending into the last interglacial periodNature43114715110.1038/nature02805 (doi:10.1038/nature02805)15356621|
|11.||McManus J. F.,Oppo D. W.,Cullen J. L.. Year: 1999A 0.5-million-year record of millennial-scale climate variability in the North AtlanticScience28397197510.1126/science.283.5404.971 (doi:10.1126/science.283.5404.971)9974387|
|12.||Johnsen S. J.,et al. Year: 1992Irregular glacial interstadials recorded in a new Greenland ice coreNature35931131310.1038/359311a0 (doi:10.1038/359311a0)|
|13.||Dansgaard W.,et al. Year: 1993Evidence for general instability of past climate from a 250-kyr ice-core recordNature36421822010.1038/364218a0 (doi:10.1038/364218a0)|
|14.||Chappellaz J.,Bluniert T.,Raynaud D.,Barnola J. M.,Schwander J.,Stauffert B.. Year: 1993Synchronous changes in atmospheric CH4 and Greenland climate between 40 and 8 kyr BPNature36644344510.1038/366443a0 (doi:10.1038/366443a0)|
|15.||Capron E.,et al. Year: 2010Millennial and sub-millennial scale climatic variations recorded in polar ice cores over the last glacial periodClim. Past634536510.5194/cp-6-345-2010 (doi:10.5194/cp-6-345-2010)|
|16.||Loulergue L.,et al. Year: 2008Orbital and millennial-scale features of atmospheric CH4 over the past 800,000 yearsNature45338338610.1038/nature06950 (doi:10.1038/nature06950)18480822|
|17.||Heinrich H.. Year: 1988Origin and consequences of cyclic ice rafting in the Northeast Atlantic Ocean during the past 130 000 yearsQuat. Res.2914215210.1016/0033-5894(88)90057-9 (doi:10.1016/0033-5894(88)90057-9)|
|18.||Bond G.,et al. Year: 1992Evidence for massive discharges of icebergs into the north atlantic ocean during the last glacial periodNature36024524910.1038/360245a0 (doi:10.1038/360245a0)|
|19.||Grousset F. E.,Labeyrie L.,Sinko J. A.,Cremer M.,Bond G.,Duprat J.,Cortijo E.,Huon S.. Year: 1993Patterns of ice-rafted detritus in the glacial North Atlantic (40-55N)Paleoceanography817519210.1029/92PA02923 (doi:10.1029/92PA02923)|
|20.||Voelker A. H. L.. Year: 2002Global distribution of centennial-scale records for marine isotope stage (MIS) 3: a databaseQuat. Sci. Rev.211185121210.1016/S0277-3791(01)00139-1 (doi:10.1016/S0277-3791(01)00139-1)|
|21.||Bond G.,et al. Year: 1997A pervasive millennial-scale cycle in North Atlantic Holocene and glacial climatesScience2781257126610.1126/science.278.5341.1257 (doi:10.1126/science.278.5341.1257)|
|22.||Bianchi G. G.,McCave I. N.. Year: 1999Holocene periodicity in north Atlantic climate and deep-ocean flow south of IcelandNature39751551710.1038/17362 (doi:10.1038/17362)|
|23.||Paul A.,Schulz M.. Wefer G.,Berger W. H.,Behre K. E.,Jansen E.Year: 2002Holocene climate variability on centennial-to-millennial time scales. I. Climate records from the North Atlantic realmClimate development and history of the North Atlantic realm4154Berlin, GermanySpringer|
|24.||Oerlemans J.. Year: 1980Model experiments on the 100,000-yr glacial cycleNature28743043210.1038/287430a0 (doi:10.1038/287430a0)|
|25.||Oerlemans J.. Year: 1982Glacial cycles and ice sheet modellingClim. Change4353374|
|26.||Ghil M.,Le Treut H.. Year: 1981A climate model with cryodynamics and geodynamicsJ. Geophys. Res.865262527010.1029/JC086iC06p05262 (doi:10.1029/JC086iC06p05262)|
|27.||Le Treut H.,Ghil M.. Year: 1983Orbital forcing, climatic interactions and glaciation cyclesJ. Geophys. Res.885167519010.1029/JC088iC09p05167 (doi:10.1029/JC088iC09p05167)|
|28.||Ghil M.,Childress S.. Year: 1987Topics in geophysical fluid dynamics: atmospheric dynamics, dynamo theory and climate dynamicsBerlin, GermanySpringer Applied Mathematical Sciences vol. 60.|
|29.||Saltzman B.,Maasch K. A.. Year: 1988Carbon cycle instability as a cause of the late Pleistocene ice age oscillations: modeling the asymmetric responseGlobal Biogeochem. Cycles211718510.1029/GB002i002p00177 (doi:10.1029/GB002i002p00177)|
|30.||Timmermann A.. Year: 2003Decadal ENSO amplitude modulations: a nonlinear paradigmGlobal Planet. Change3713515610.1016/S0921-8181(02)00194-7 (doi:10.1016/S0921-8181(02)00194-7)|
|31.||Schulz M.,Prange M.,Klocker A.. Year: 2007Low-frequency oscillations of the Atlantic Ocean meridional overturning circulation in a coupled climate modelClim. Past39710710.5194/cp-3-97-2007 (doi:10.5194/cp-3-97-2007)|
|32.||Sakai K.,Peltier W. R.. Year: 1999A dynamical systems model of the Dansgaard–Oeschger oscillation and the origin of the Bond cycleJ. Clim.122238225510.1175/1520-0442(1999)012%3C2238:ADSMOT%3E2.0.CO;2 (doi:10.1175/1520-0442(1999)012<2238:ADSMOT>2.0.CO;2)|
|33.||Colin de Verdière A.. Year: 2007A simple model of millennial oscillations of the thermohaline circulationJ. Phys. Oceanogr.371142115510.1175/JPO3056.1 (doi:10.1175/JPO3056.1)|
|34.||Rahmstorf S.,et al. Year: 2005Thermohaline circulation hysteresis: a model intercomparisonGeophys. Res. Lett.32L2360510.1029/2005GL023655 (doi:10.1029/2005GL023655)|
|35.||Calov R.,Ganopolski A.. Year: 2005Multistability and hysteresis in the climate-cryosphere system under orbital forcingGeophys. Res. Lett.32L2171710.1029/2005GL024518 (doi:10.1029/2005GL024518)|
|36.||Timmermann A.,Jin F. F.. Palmer T.,Hagedorn R.Year: 2006Predictability of coupled processesPredictability of weather and climate251274Cambridge, UKCambridge University Press|
|37.||Lorenz E. N.. Year: 1975Climate predictabilityThe physical basis of climate and climate modelling132136Geneva, SwitzerlandWorld Meteorological Organisation WMO GARP Publication Series, No. 16, Appendix 2.1.|
|38.||Strogatz S. H.. Year: 1994Nonlinear dynamics and chaos, with applications to physics, biology, chemistry, and engineering (studies in nonlinearity)Reading, MAAddison-Wesley|
|39.||Guckenheimer J.,Holmes P.. Year: 1983Nonlinear oscillations, dynamical systems, and bifurcations of vector fieldsNew York, NYSpringer Applied Mathematical Sciences vol. 42.|
|40.||Saltzman B.. Year: 2001Dynamical paleoclimatology: generalized theory of global climate change: international geophysicsNew York, NYAcademic Press International Geophysics Series vol. 80.|
|41.||Pikovski A.,Rosenblum M.,Kurths J.. Year: 2001Synchronization: a universal concept in nonlinear sciencesCambridge, UKCambridge University Press Cambridge Nonlinear Science Series vol. 12.|
|42.||Saltzman B.,Hansen A. R.,Maasch K. A.. Year: 1984The late Quaternary glaciations as the response of a 3-component feedback-system to Earth-orbital forcingJ. Atmos. Sci.413380338910.1175/1520-0469(1984)041%3C3380:TLQGAT%3E2.0.CO;2 (doi:10.1175/1520-0469(1984)041<3380:TLQGAT>2.0.CO;2)|
|43.||Saltzman B.,Sutera A.. Year: 1984A model of the internal feedback system involved in late Quaternary climatic variationsJ. Atmos. Sci.4173674510.1175/1520-0469(1984)041%3C0736:AMOTIF%3E2.0.CO;2 (doi:10.1175/1520-0469(1984)041<0736:AMOTIF>2.0.CO;2)|
|44.||Saltzman B.,Sutera A.. Year: 1987The mid-Quaternary climatic transition as the free response of a three-variable dynamical modelJ. Atmos. Sci.4423624110.1175/1520-0469(1987)044%3C0236:TMQCTA%3E2.0.CO;2 (doi:10.1175/1520-0469(1987)044<0236:TMQCTA>2.0.CO;2)|
|45.||Saltzman B.,Verbitsky M. Y.. Year: 1993Multiple instabilities and modes of glacial rhythmicity in the Plio-Pleistocene: a general theory of late Cenozoic climatic changeClim. Dyn.911510.1007/BF00208010 (doi:10.1007/BF00208010)|
|46.||Saltzman B.,Maasch K. A.. Year: 1990A first-order global model of late Cenozoic climateTrans. R. Soc. Edinburgh Earth Sci.8131532510.1017/S0263593300020824 (doi:10.1017/S0263593300020824)|
|47.||Saltzman B.,Maasch K. A.. Year: 1991A first-order global model of late Cenozoic climate II. Further analysis based on a simplification of the CO2 dynamicsClim. Dyn.520121010.1007/BF00210005 (doi:10.1007/BF00210005)|
|48.||Murphy J. J.. Year: 1876The glacial climate and the polar ice-capQ. J. Geol. Soc. Lond.3240040610.1144/GSL.JGS.1876.032.01-04.45 (doi:10.1144/GSL.JGS.1876.032.01-04.45)|
|49.||Raymo M. E.,Ruddiman W. F.,Froelich P. N.. Year: 1988Influence of late Cenozoic mountain building on ocean geochemical cyclesGeology1664965310.1130/0091-7613(1988)016%3C0649:IOLCMB%3E2.3.CO;2 (doi:10.1130/0091-7613(1988)016<0649:IOLCMB>2.3.CO;2)|
|50.||Doedel E. J.,Oldeman B. E.. Year: 2009AUTO-07P: continuation and bifurcation software for ordinary differential equations Montreal, Canada:Concordia University, Technical report, See http://indy.cs.concordia.ca/auto.|
|51.||Hays J. D.,Imbrie J.,Shackleton N. J.. Year: 1976Variations in the Earth's orbit: pacemaker of ice agesScience1941121113210.1126/science.194.4270.1121 (doi:10.1126/science.194.4270.1121)17790893|
|52.||Imbrie J.,Imbrie J. Z.. Year: 1980Modelling the climatic response to orbital variationsScience20794395310.1126/science.207.4434.943 (doi:10.1126/science.207.4434.943)17830447|
|53.||Raymo M.. Year: 1997The timing of major climate terminationsPaleoceanography1257758510.1029/97PA01169 (doi:10.1029/97PA01169)|
|54.||Huybers P.. Year: 2007Glacial variability over the last two million years: an extended depth-derived age model, continuous obliquity pacing, and the Pleistocene progressionQuat. Sci. Rev.26375510.1016/j.quascirev.2006.07.013 (doi:10.1016/j.quascirev.2006.07.013)|
|55.||Paillard D.,Labeyrie L.. Year: 1994Role of the thermohaline circulation in the abrupt warming after Heinrich eventsNature37216216410.1038/372162a0 (doi:10.1038/372162a0)|
|56.||Paillard D.. Year: 1998The timing of Pleistocene glaciations from a simple multiple-state climate modelNature39137838110.1038/34891 (doi:10.1038/34891)|
|57.||Guckenheimer J.,Hoffman K.,Weckesser W.. Year: 2003The forced van der Pol equation. I. The slow flow and its bifurcationsSIAM J. Appl. Dyn. Syst.213510.1137/S1111111102404738 (doi:10.1137/S1111111102404738)|
|58.||Paillard D.. Year: 2001Glacial cycles: toward a new paradigmRev. Geophys.3932534610.1029/2000RG000091 (doi:10.1029/2000RG000091)|
|59.||Gildor H.,Tziperman E.. Year: 2000Sea ice as the glacial cycles climate switch: role of seasonal and orbital forcingPaleoceanography1560561510.1029/1999PA000461 (doi:10.1029/1999PA000461)|
|60.||Gildor H.,Tziperman E.. Year: 2001Physical mechanisms behind biogeochemical glacial-interglacial CO2 variationsGeophys. Res. Lett.282421242410.1029/2000GL012571 (doi:10.1029/2000GL012571)|
|61.||Tziperman E.,Raymo M. E.,Huybers P.,Wunsch C.. Year: 2006Consequences of pacing the Pleistocene 100 kyr ice ages by nonlinear phase locking to Milankovitch forcingPaleoceanography211110.1029/2005PA001241 (doi:10.1029/2005PA001241)|
|62.||Paillard D.,Parrenin F.. Year: 2004The Antarctic ice sheet and the triggering of deglaciationsEarth Planet. Sci. Lett.22726327110.1016/j.epsl.2004.08.023 (doi:10.1016/j.epsl.2004.08.023)|
|63.||Cane M. A.,et al. Year: 2006Progress in paleoclimate modelingJ. Clim.195031505710.1175/JCLI3899.1 (doi:10.1175/JCLI3899.1)|
|64.||Crucifix M.. Year: 2011How can a glacial inception be predicted?Holocene2183184210.1177/0959683610394883 (doi:10.1177/0959683610394883)|
|65.||Welander P.. Year: 1982A simple heat–salt oscillatorDyn. Atmos. Oceans623324210.1016/0377-0265(82)90030-6 (doi:10.1016/0377-0265(82)90030-6)|
|66.||Cessi P.. Year: 1996Convective adjustment and thermohaline excitabilityJ. Phys. Oceanogr.2648149110.1175/1520-0485(1996)026%3C0481:CAATE%3E2.0.CO;2 (doi:10.1175/1520-0485(1996)026<0481:CAATE>2.0.CO;2)|
|67.||Welander P.. Willebrand J.,Anderson D. L. T.Year: 1986Thermohaline effects in the ocean circulation and related simple modelsLarge-scale transport processes in oceans and atmosphere163200Dordrecht, GermanyReidel|
|68.||Winton M.,Sarachik E. S.. Year: 1993Thermohaline oscillations induced by strong steady salinity forcing of ocean general circulation modelsJ. Phys. Oceanogr.231389141010.1175/1520-0485(1993)023%3C1389:TOIBSS%3E2.0.CO;2 (doi:10.1175/1520-0485(1993)023<1389:TOIBSS>2.0.CO;2)|
|69.||Schulz M.,Paul A.,Timmermann A.. Year: 2002Relaxation oscillators in concert: a framework for climate change at millennial timescales during the late PleistoceneGeophys. Res. Lett.29219310.1029/2002GL016144 (doi:10.1029/2002GL016144)|
|70.||Weaver A. J.,Hughes T. M. C.. Year: 1994Rapid interglacial climate fluctuations driven by north Atlantic ocean circulationNature36744745010.1038/367447a0 (doi:10.1038/367447a0)|
|71.||Haarsma R. J.,Opsteegh J. D.,Selten F. M.,Wang X.. Year: 2001Rapid transitions and ultra-low frequency behaviour in a 40 kyr integration with a coupled climate model of intermediate complexityClim. Dyn.1755957010.1007/s003820000129 (doi:10.1007/s003820000129)|
|72.||Meissner K.,Eby M.,Weaver A.,Saenko O.. Year: 2008CO2 threshold for millennial-scale oscillations in the climate system: implications for global warming scenariosClim. Dyn.3016117410.1007/s00382-007-0279-0 (doi:10.1007/s00382-007-0279-0)|
|73.||Sakai K.,Peltier W. R.. Year: 1997Dansgaard–Oeschger oscillations in a coupled atmosphere–ocean climate modelJ. Clim.1094997010.1175/1520-0442(1997)010%3C0949:DOOIAC%3E2.0.CO;2 (doi:10.1175/1520-0442(1997)010<0949:DOOIAC>2.0.CO;2)|
|74.||Colin de Verdière A.,Ben Jelloul M.,Sévellec F.. Year: 2006Bifurcation structure of thermohaline millennial oscillationsJ. Clim.195777579510.1175/JCLI3950.1 (doi:10.1175/JCLI3950.1)|
|75.||Colin de Verdière A.,Te Raa L.. Year: 2010Weak oceanic heat transport as a cause of the instability of glacial climatesClim. Dyn.351237125610.1007/s00382-009-0675-8 (doi:10.1007/s00382-009-0675-8)|
|76.||Timmermann A.,Gildor H.,Schulz M.,Tziperman E.. Year: 2003Coherent resonant millennial-scale climate oscillations triggered by massive meltwater pulsesJ. Clim.162569258510.1175/1520-0442(2003)016%3C2569:CRMCOT%3E2.0.CO;2 (doi:10.1175/1520-0442(2003)016<2569:CRMCOT>2.0.CO;2)|
|77.||Ganopolski A.,Rahmstorf S.. Year: 2002Abrupt glacial climate changes due to stochastic resonancePhys. Rev. Lett.8803850110.1103/PhysRevLett.88.038501 (doi:10.1103/PhysRevLett.88.038501)11801092|
|78.||Ganopolski A.,Rahmstorf S.. Year: 2001Rapid changes of glacial climate simulated in a coupled climate modelNature40915315810.1038/35051500 (doi:10.1038/35051500)11196631|
|79.||Schulz M.. Year: 2002On the 1470-year pacing of Dansgaard–Oeshger warm eventsPaleoceanography17101410.1029/2000PA000571 (doi:10.1029/2000PA000571)|
|80.||Braun H.,Kurths J.. Year: 2010Were Dansgaard–Oeschger events forced by the sun?Eur. Phys. J. Spec. Top.19111712910.1140/epjst/e2010-01345-5 (doi:10.1140/epjst/e2010-01345-5)|
|81.||Ditlevsen P. D.,Ditlevsen O. D.. Year: 2009On the stochastic nature of the rapid climate shifts during the last ice ageJ. Clim.2244645710.1175/2008JCLI2430.1 (doi:10.1175/2008JCLI2430.1)|
|82.||Ganopolski A.,Rahmstorf S.. Year: 2002Abrupt glacial climate changes due to stochastic resonancePhys. Rev. Lett.8803850110.1103/PhysRevLett.88.038501 (doi:10.1103/PhysRevLett.88.038501)11801092|
|83.||Braun H.,Christl M.,Rahmstorf S.,Ganopolski A.,Mangini A.,Kubatzki C.,Roth K.,Kromer B.. Year: 2005Possible solar origin of the 1,470-year glacial climate cycle demonstrated in a coupled modelNature43820821110.1038/nature04121 (doi:10.1038/nature04121)16281042|
|84.||Braun H.,Ditlevsen P.,Chialvo D. R.. Year: 2008Solar forced Dansgaard–Oeschger events and their phase relation with solar proxiesGeophys. Res. Lett.35510.1029/2008GL033414 (doi:10.1029/2008GL033414)|
|85.||MacAyeal D.. Year: 1993Binge/purge oscillations of the Laurentide ice sheet as a cause of the North Atlantic's Heinrich eventsPaleoceanography877578410.1029/93PA02200 (doi:10.1029/93PA02200)|
|86.||Paillard D.. Year: 1995The hierarchical structure of glacial climatic oscillations: interactions between ice-sheet dynamics and climateClim. Dyn.1116217710.1007/BF00223499 (doi:10.1007/BF00223499)|
|87.||Schmittner A.,Yoshimori M.,Weaver A. J.. Year: 2002Instability of glacial climate in a model of the ocean–atmosphere–cryosphere systemScience2951489149310.1126/science.1066174 (doi:10.1126/science.1066174)11823604|
|88.||Ganopolski A.,Calov R.,Claussen M.. Year: 2010Simulation of the last glacial cycle with a coupled climate ice-sheet model of intermediate complexityClim. Past622924410.5194/cp-6-229-2010 (doi:10.5194/cp-6-229-2010)|
|89.||Alvarez-Solas J.,Charbit S.,Ritz C.,Paillard D.,Ramstein G.. Year: 2010Links between ocean temperature and iceberg discharge during Heinrich eventsNat. Geosci.312212610.1038/NGEO752 (doi:10.1038/NGEO752)|
|90.||Schulz M.,Paul A.,Timmermann A.. Year: 2004Glacial-interglacial contrast in climate variability at centennial-to-millennial timescales: observations and conceptual modelQuat. Sci. Rev.232219223010.1016/j.quascirev.2004.08.014 (doi:10.1016/j.quascirev.2004.08.014)|
|91.||Hasselmann K.. Year: 1976Stochastic climate models. I. TheoryTellus2847348510.1111/j.2153-3490.1976.tb00696.x (doi:10.1111/j.2153-3490.1976.tb00696.x)|
|92.||Penland C.. Year: 2003Noise out of chaos and why it won't go awayBull. Am. Meteorol. Soc.8492192510.1175/BAMS-84-7-921 (doi:10.1175/BAMS-84-7-921)|
|93.||Penland C.. Tsonis A. A.,Elsner J. B.Year: 2007Stochastic linear models of nonlinear geosystemsNonlinear dynamics in geosciences485515New York, NYSpringer10.1007/978-0-387-34918-3_27 (doi:10.1007/978-0-387-34918-3_27)|
|94.||Lindner B.,Garcìa-Ojalvo J.,Neiman A.,Schimansky-Geier L.. Year: 2004Effects of noise in excitable systemsPhys. Rep.39232142410.1016/j.physrep.2003.10.015 (doi:10.1016/j.physrep.2003.10.015)|
|95.||Nicolis C.. Nicolis C.,Nicolis G.Year: 1987Climate predictability and dynamical systemsIrreversible phenomena and dynamical system analysis in the geosciencesDordrecht, The NetherlandsKluwer321354 NATO ASI Series C: Mathematical and Physical Sciences, vol. 192.|
|96.||Stone E.,Holmes P.. Year: 1990Random perturbations of heteroclinic attractorsSIAM J. Appl. Math.5072674310.1137/0150043 (doi:10.1137/0150043)|
|97.||Grasman J.,Roerdink J. B. T. M.. Year: 1989Stochastic and chaotic relaxation oscillationsJ. Stat. Phys.5494997010.1007/BF01019783 (doi:10.1007/BF01019783)|
|98.||Berglund N.,Gentz B.. Year: 2002The effect of additive noise on dynamical hysteresisNonlinearity1560510.1088/0951-7715/15/3/305 (doi:10.1088/0951-7715/15/3/305)|
|99.||Toggweiler J. R.. Year: 2008Origin of the 100,000-year timescale in Antarctic temperatures and atmospheric CO2Paleoceanography231710.1029/2006PA001405 (doi:10.1029/2006PA001405)|
|100.||Braun H.,Ganopolski A.,Christl M.,Chialvo D. R.. Year: 2007A simple conceptual model of abrupt glacial climate eventsNonlinear Process. Geophys.1470972110.5194/npg-14-709-2007 (doi:10.5194/npg-14-709-2007)|
|101.||Wunsch C.. Year: 2010Towards understanding the paleoceanQuat. Sci. Rev.291960196710.1016/j.quascirev.2010.05.020 (doi:10.1016/j.quascirev.2010.05.020)|
|102.||Livina V. N.,Kwasniok F.,Lenton T. M.. Year: 2010Potential analysis reveals changing number of climate states during the last 60 kyrClim. Past6778210.5194/cp-6-77-2010 (doi:10.5194/cp-6-77-2010)|
|103.||Livina V. N.,Kwasniok F.,Lohmann G.,Kantelhardt J. W.,Lenton T. M.. Year: 2011Changing climate states and stability: from Pliocene to presentClim. Dyn372437245310.1007/s00382-010-0980-2 (doi:10.1007/s00382-010-0980-2)|
|104.||Crucifix M.,Rougier J.. Year: 2009On the use of simple dynamical systems for climate predictions: a Bayesian prediction of the next glacial inceptionEur. Phys. J. Spec. Top.174113110.1140/epjst/e2009-01087-5 (doi:10.1140/epjst/e2009-01087-5)|
[Figure ID: RSTA20110315F1]
Climatic fluctuations over the Late Pleistocene. Ice ages are reconstructed using the oxygen isotopic ratio of the calcite shells of benthic foraminifera . Within the last ice ages, large temperature fluctuations were recorded in Greenland, here depicted by changes in the oxygen isotopic ratio of ice .
[Figure ID: RSTA20110315F2]
Schematic of several forms of relaxation oscillations. (a) The vector space is structured by a slow manifold with several stable branches. All points of the state space are attracted towards the stable branches of the slow manifold (solid lines) along the fast direction, which is here the horizontal. The slow evolution consists of an upward or downward course along the slow manifold depending on whether the system lives on the right- or left-hand side of the null-cline of the slow variable. The relaxation oscillation consists of alternate jumps between the two branches of the slow manifold. (b) Trajectories are rapidly attracted towards a region of the phase space influenced by a saddle point. In the scenario displayed here, there exists a combination of parameters for which the limit cycle crosses the saddle point. In that case, the period of the orbit, which includes the saddle–node, is infinitely long. It is a ‘homoclinic orbit’, hence this particular bifurcation is named a homoclinic bifurcation. The homoclinic orbit only exists at the bifurcation point, but it influences the orbit when the parameter is close to the bifurcation. This is the reason why one refers to the ‘ghost’ of the homoclinic orbit. There is another scenario for which different saddle–nodes are connected to each other. The orbit and the associated bifurcation are then said to be ‘heteroclinic’. (c) The relaxation oscillation is organized around a fixed point, with complex eigenvalues with a positive real part. The bifurcation giving rise to this orbit is a Hopf bifurcation. (d) One example of excitability, here depicted for a slow–fast system. The system resides in a stable space, but a fluctuation may cause an ejection out of the unstable (dashed lines) branch of the slow manifold. The system then loops all the way through the slow manifold before coming back to rest. (Online version in colour.)
[Figure ID: RSTA20110315F3]
Bifurcation diagrams of the Saltzman and Maasch models (a) SM90 and (b) SM91 as a function of Fμ, here treated as a constant control parameter. Note the difference in scales on both axes. Black lines are fixed points. Continental ice mass (I) is shown as a function of tectonic forcing Fμ. The lines denoted ‘max’ and ‘min’ indicate the boundaries of the limit cycles. Unstable fixed points or limit cycles are denoted by dashed lines. Calculations and figures were made using the pseudo-arc-length continuation software package AUTO . (Online version in colour.)
[Figure ID: RSTA20110315F4]
(a,b) Periodic and fixed-point solutions for SM90 (with Fμ=1.2 ppm ka1) and (c,d) SM91 (with Fμ=0.5 ppm ka1), near the sub-critical Hopf bifurcation points. The dynamics of SM90 slow down near the unstable fixed point, whereas the limit cycle of SM91 is much more decoupled from the position of the fixed point. (Online version in colour.)
[Figure ID: RSTA20110315F5]
Two histories of ice volume generated with the same models (a) SM90 and (b) SM91, using the astronomical forcing and a decline scenario for the tectonic forcing of CO2. The scenarios were chosen here to evidence the rise and decline of 100 000 year ice ages and are not the same as those used by Saltzman. Observe the explosive character of the appearance and decline of 100 000 year cycles in SM91, with early and late excitations of the limit cycle by the astronomical forcing, and the smoother character of the evolution on SM90. Time is running from right to left.
[Figure ID: RSTA20110315F6]
Bifurcation diagrams and phase-space trajectories of the free Paillard–Parrenin model (a,b) PP04 and (c,d) the van der Pol (VDP) model. Both systems display explosive bifurcation scenarios, but the details are different. PP04 exhibits a subcritical Hopf bifurcation, whereas the limit cycle of VDP is explicitly framed by a slow manifold. Phase-space trajectories are drawn near the bifurcation points, that is: g=0.4 in PP04 and β=0.9 in VDP. The Heaviside function in PP04 is approximated as H(x)=a tan(500x)/π+0.5 for analysis by AUTO. (Online version in colour.)
[Figure ID: RSTA20110315F7]
Astronomical forcing, x and y trajectories obtained using system (van der Pol, VDP) with α=30, β=0.75, γ=0.4 and τ=36 ka (1 ka=1000 years). Black dots are an authoritative natural archive thought to mainly represent fluctuations in ice volume and deep-ocean temperature, and compiled in Lisiecki & Raymo .
[Figure ID: RSTA20110315F8]
The phenomenon of trajectory shift illustrated with the model (van der Pol, VDP) with astronomical forcing, with parameters τ=36, β=0.75, γ= 0.4, α=30 ka. The black line represents the deterministic system. The generated history reproduces reasonably well the fluctuations of ice volume of the last 800 000 years. The red (grey) line represents the same system but with an additive Wiener process added to the fast variable y, with variance b=0.2 (ka)−1/2. One possible realization of the stochastic system is shown. Observe the solution shift around 450 000 years ago. Depending on the realization, the shift may occur at different places. (Online version in colour.)
[Figure ID: RSTA20110315F9]
Two possible effects of additive noisy fluctuations in a slow–fast system. (a) If the slow–fast system is in the excitable regime, then noisy fluctuations can cause excitations. Excitations will be sporadic if the amplitude of the oscillations is weak (as shown here), or regular, as in a limit cycle, in the coherent resonance regime. (b) If the slow–fast system is oscillating, then additive fluctuations may cause early or delayed escapes from the slow branches, compared with the deterministic system. The net effect is a shortening of the cycle duration. (Online version in colour.)
Keywords: palaeoclimates, dynamical systems, limit cycle, ice ages, Dansgaard–Oeschger events.
Previous Document: ELPIS-JP: a dataset of local-scale daily climate change scenarios for Japan.
Next Document: Tipping points in open systems: bifurcation, noise-induced and rate-dependent examples in the climat...