Internal Climate Dynamics Research Articles

Orbital- and millennial-scale Asian winter monsoon variability across the Pliocene–Pleistocene glacial intensification

Intensification of northern hemisphere glaciation (iNHG), ~2.7 million years ago (Ma), led to establishment of the Pleistocene to present-day bipolar icehouse state. Here we document evolution of orbital- and millennial-scale Asian winter monsoon (AWM) variability across the iNHG using a palaeomagnetically dated centennial-resolution grain size record between 3.6 and 1.9 Ma from a previously undescribed loess-palaeosol/red clay section on the central Chinese Loess Plateau. We find that the late Pliocene–early Pleistocene AWM was characterized by combined 41-kyr and ~100-kyr cycles, in response to ice volume and atmospheric CO2 forcing. Northern hemisphere ice sheet expansion, which was accompanied by an atmospheric CO2 concentration decline, substantially increased glacial AWM intensity and its orbitally oscillating amplitudes across the iNHG. Superposed on orbital variability, we find that millennial AWM intensity fluctuations persisted during both the warmer (higher-CO2) late Pliocene and colder (lower-CO2) early Pleistocene, in response to both external astronomical forcing and internal climate dynamics.

A Systematic Role for Extreme Ocean‐Atmosphere Oscillations in the Development of Glacial Conditions Since the Mid Pleistocene Transition

AbstractWe introduce a new hypothesis concerning the role of internal climate dynamics in the non‐linear transitions from interglacial to glacial (IG‐G) state since the Mid Pleistocene Transition (MPT). These transitions encompass large and abrupt changes in atmospheric CO2, ice volume, and temperature that we suggest involve critical interactions between insolation and high amplitude oscillations in ocean/atmosphere circulation patterns. Specifically, we highlight the large amplitude of millennial‐scale climate oscillations across the transition from Marine Isotope Stage (MIS) 5 to 4, which we argue led to amplified cooling of the deep ocean and we demonstrate that analogous episodes of extreme cooling systematically preceded glacial periods of the last 800 kyr. We suggest that such cooling necessitates a reconfiguration of the deep ocean to avoid a density paradox between northern and southern‐sourced deep waters (SSW), which could be accomplished by increasing the relative volume and or salinity of SSW, thus providing the necessary storage capacity for the subsequent (delayed) and relatively abrupt drawdown of CO2. We therefore explain the transient decoupling of Antarctic temperature from CO2 across MIS 5/4 as a direct consequence of millennial activity at that time. We further show that similar climatic decoupling typically occurred during times of low obliquity and was a ubiquitous feature of IG‐G transitions over the past 800 kyr, producing the appearance of bimodality in records of CO2, benthic δ18O and others. Finally we argue that the apparent lack of bimodality in the pre‐MPT record of benthic δ18O implies that the dynamics associated with IG‐G transitions changed across the MPT.

A 1.5-million-year record of orbital and millennial climate variability in the North Atlantic

Abstract. Climate during the last glacial period was marked by abrupt instability on millennial timescales that included large swings of temperature in and around Greenland (Daansgard–Oeschger events) and smaller, more gradual changes in Antarctica (AIM events). Less is known about the existence and nature of similar variability during older glacial periods, especially during the early Pleistocene when glacial cycles were dominantly occurring at 41 kyr intervals compared to the much longer and deeper glaciations of the more recent period. Here, we report a continuous millennially resolved record of stable isotopes of planktic and benthic foraminifera at IODP Site U1385 (the “Shackleton Site”) from the southwestern Iberian margin for the last 1.5 million years, which includes the Middle Pleistocene Transition (MPT). Our results demonstrate that millennial climate variability (MCV) was a persistent feature of glacial climate, both before and after the MPT. Prior to 1.2 Ma in the early Pleistocene, the amplitude of MCV was modulated by the 41 kyr obliquity cycle and increased when axial tilt dropped below 23.5∘ and benthic δ18O exceeded ∼3.8 ‰ (corrected to Uvigerina), indicating a threshold response to orbital forcing. Afterwards, MCV became focused mainly on the transitions into and out of glacial states (i.e. inceptions and terminations) and during times of intermediate ice volume. After 1.2 Ma, obliquity continued to play a role in modulating the amplitude of MCV, especially during times of glacial inceptions, which are always associated with declining obliquity. A non-linear role for obliquity is also indicated by the appearance of multiples (82, 123 kyr) and combination tones (28 kyr) of the 41 kyr cycle. Near the end of the MPT (∼0.65 Ma), obliquity modulation of MCV amplitude wanes as quasi-periodic 100 kyr and precession power increase, coinciding with the growth of oversized ice sheets on North America and the appearance of Heinrich layers in North Atlantic sediments. Whereas the planktic δ18O of Site U1385 shows a strong resemblance to Greenland temperature and atmospheric methane (i.e. Northern Hemisphere climate), millennial changes in benthic δ18O closely follow the temperature history of Antarctica for the past 800 kyr. The phasing of millennial planktic and benthic δ18O variation is similar to that observed for MIS 3 throughout much of the record, which has been suggested to mimic the signature of the bipolar seesaw – i.e. an interhemispheric asymmetry between the timing of cooling in Antarctica and warming in Greenland. The Iberian margin isotopic record suggests that bipolar asymmetry was a robust feature of interhemispheric glacial climate variations for at least the past 1.5 Ma despite changing glacial boundary conditions. A strong correlation exists between millennial increases in planktic δ18O (cooling) and decreases in benthic δ13C, indicating that millennial variations in North Atlantic surface temperature are mirrored by changes in deep-water circulation and remineralization of carbon in the abyssal ocean. We find strong evidence that climate variability on millennial and orbital scales is coupled across different timescales and interacts in both directions, which may be important for linking internal climate dynamics and external astronomical forcing.

Forced signal and predictability in a prototype climate model: Implications for fingerprinting based detection in the presence of multidecadal natural variability.

Advanced numerical models used for climate prediction are known to exhibit biases in their simulated climate response to variable concentrations of the atmospheric greenhouse gases and aerosols that force a non-uniform, in space and time, secular global warming. We argue here that these biases can be particularly pronounced due to misrepresentation, in these models, of the multidecadal internal climate variability characterized by large-scale, hemispheric-to-global patterns. This point is illustrated through the development and analysis of a prototype climate model comprised of two damped linear oscillators, which mimic interannual and multidecadal internal climate dynamics and are set into motion via a combination of stochastic driving, representing weather noise, and deterministic external forcing inducing a secular climate change. The model time series are paired with pre-specified patterns in the physical space and form, conceptually, a spatially extended time series of the zonal-mean near-surface temperature, which is further contaminated by a spatiotemporal noise simulating the rest of climate variability. The choices of patterns and model parameters were informed by observations and climate-model simulations of the 20th century near-surface air temperature. Our main finding is that the intensity and spatial patterns of the internal multidecadal variability associated with the slow-oscillator model component greatly affect (i) the ability of modern pattern-recognition/fingerprinting methods to isolate the forced response of the climate system in the 20th century ensemble simulations and (ii) climate-system predictability, especially decadal predictability, as well as the estimates of this predictability using climate models in which the internal multidecadal variability is underestimated or otherwise misrepresented.

Variability of Multisite Decadal Running Means Arising from Year-To-Year Fluctuations

AbstractDecadal mean variables are frequently used to characterize decadal climate variabilities. Decadal means are often calculated using yearly data, which can represent variability at time scales from annual to centennial. Residuals from interannual fluctuations may contribute to the variability in decadal time series. Such variability is more difficult to be predicted at the long range. Removing it from the decadal variability means that the remaining variability is more likely to arise from slowly varying multidecadal or longer time scale external forcing and internal climate dynamics, which are more likely to be predicted. Here, a new approach is proposed to understand the uncertainty, potential predictability, and drivers of decadal mean variables. The covariance matrix of multivariate decadal running means is decomposed into unpredictable fast decadal variability and the potentially predictable slow decadal variability. EOF analysis is then applied to the decomposed matrices to find the dominant modes, which may be related to the drivers of the two types of variabilities in the multivariate decadal means. The methodology has been applied to 140-yr datasets of North Pacific sea surface temperature and the Northern Hemisphere 1000-hPa geopotential height. For sea surface temperature, the Pacific decadal oscillation is the major driver of the fast decadal variability, while the radiative forcing and the Atlantic multidecadal oscillation are major drivers of the slow decadal variability. For the 1000-hPa geopotential height, fast decadal variability is associated with the northern annular mode, the east Atlantic mode, and the Pacific decadal oscillation. Slow decadal variability is associated with the northern annular mode and the Atlantic multidecadal oscillation.

2,200-Year tree-ring and lake-sediment based snowpack reconstruction for the northern Rocky Mountains highlights the historic magnitude of recent snow drought

In recent decades, Rocky Mountain accumulated snowpack levels have experienced rapid declines, yet long-term records of snowpack prior to the installation of snowpack observation stations in the early and mid 20th century are limited. To date, a small number of tree-ring based reconstructions of April 1 Snow Water Equivalent (SWE) in the northern Rocky Mountains have extended modern records of snowpack variability to ∼1200 C.E. Carbonate isotope lake sediment records, provide an opportunity to further extend tree-ring based reconstructions through the Holocene, providing a millennial-scale temporal record that allows for an evaluation of multi-scale drivers of snowpack variability, from internal climate dynamics to orbital-scale forcings. Here we present a ∼2200 year preliminary reconstruction of northern Rockies snowpack based on δ18O measurements of sediment carbonates collected from Foy Lake, Montana. We explore the statistical calibration of lake sediment δ18O to an annually resolved snowpack reconstruction from tree rings, and develop an approach to assess and quantify potential sources of error in this reconstruction approach. The sediment-based snowpack reconstruction shows strong low-frequency variability in snowpack over the last two millennia with few snow droughts approaching the magnitude of recent snowpack declines. Given the growing availability of high-resolution, carbonate-rich lake sediment records, such reconstructions could help improve our understanding of how snowpack conditions varied under previous climatic events (mid-Holocene climate optimum ca. 9−6 ka), providing critical insights for anticipating future snowpack conditions.

Drought Cycles Over the Last 8,200 Years Recorded in Maar Lake Twintaung, Myanmar

AbstractThe tropical Asian regions are particularly prone to catastrophic droughts due to high temperature and evaporation and their sensitivity to variable oceanic‐atmospheric circulation. The extent to which future global warming may intensify droughts in this populous region is a matter of utmost concern. Palaeoclimatic data can help inform scientists learn more about the temporal patterns and drivers of monsoon change over geological timescales. In this paper, we present a new 8,200 year biomarker record preserved in the annually laminated sediments of Maar Lake Twintaung, Myanmar. The δ13C27‐35 sedimentary record is mainly derived from leaf wax lipids in the lake catchment, which is presently dominated by C3 vegetation in the modern lake catchment. The δ13C27‐35 is mainly regulated by plant physiological and biochemical responses to drought stress and, as such, can be used as a proxy of past changes in moisture. The δ13C27–35 record reveals distinct decadal‐to‐centennial‐scale droughts superimposed on a trend of gradually decreasing summer monsoon intensity from the mid‐Holocene to late Holocene. Within the limits of the dating uncertainties, these decadal‐to‐centennial‐scale droughts are found to be well correlated with the southward shift of the Intertropical Convergence Zone. It highlights the importance of the Intertropical Convergence Zone shift (tropical monsoon trough) in regulating monsoon rainfall in the region. In addition, most of the droughts are likely to be linked with active volcanism and solar minima and suggest a coupled process between external drivers and internal climate dynamics.

The Origin and Propagation of the Antarctic Centennial Oscillation

The Antarctic Centennial Oscillation (ACO) is a paleoclimate temperature cycle that originates in the Southern Hemisphere, is the presumptive evolutionary precursor of the contemporary Antarctic Oscillation (AAO), and teleconnects to the Northern Hemisphere to influence global temperature. In this study we investigate the internal climate dynamics of the ACO over the last 21 millennia using stable water isotopes frozen in ice cores from 11 Antarctic drill sites as temperature proxies. Spectral and time series analyses reveal that ACOs occurred at all 11 sites over all time periods evaluated, suggesting that the ACO encompasses all of Antarctica. From the Last Glacial Maximum through the Last Glacial Termination (LGT), ACO cycles propagated on a multicentennial time scale from the East Antarctic coastline clockwise around Antarctica in the streamline of the Antarctic Circumpolar Current (ACC). The velocity of teleconnection (VT) is correlated with the geophysical characteristics of drill sites, including distance from the ocean and temperature. During the LGT, the VT to coastal sites doubled while the VT to inland sites decreased fourfold, correlated with increasing solar insolation at 65°N. These results implicate two interdependent mechanisms of teleconnection, oceanic and atmospheric, and suggest possible physical mechanisms for each. During the warmer Holocene, ACOs arrived synchronously at all drill sites examined, suggesting that the VT increased with temperature. Backward extrapolation of ACO propagation direction and velocity places its estimated geographic origin in the Southern Ocean east of Antarctica, in the region of the strongest sustained surface wind stress over any body of ocean water on Earth. ACO period is correlated with all major cycle parameters except cycle symmetry, consistent with a forced, undamped oscillation in which the driving energy affects all major cycle metrics. Cycle period and symmetry are not discernibly different for the ACO and AAO over the same time periods, suggesting that they are the same climate cycle. We postulate that the ACO/AAO is generated by relaxation oscillation of Westerly Wind velocity forced by the equator-to-pole temperature gradient and propagated regionally by identified air-sea-ice interactions.

Early Interglacial Legacy of Deglacial Climate Instability

AbstractThroughout the last glacial cycle millennial timescale variations in atmospheric CO2 occurred in parallel with perturbations in deep ocean circulation, which were themselves reflected by observable changes in surface conditions across the North Atlantic region. Here we use continuous proxy records to argue that an equivalent relationship has held throughout the last 800 kyr, that is, since before the first occurrence of Heinrich events (strictly defined). Our results highlight the importance of internal climate dynamics in amplifying external (insolation) forcing on the climate system to produce the large amplitude of glacial terminations (deglaciations) during the middle to late Pleistocene. We show that terminations are characterized by an interval of intense ice rafting followed by a subsequent and abrupt shift to anomalously warm surface conditions (with respect to the more gradually evolving background state), which we interpret to reflect an abrupt recovery of deep ocean circulation in the Atlantic. According to our synthesis, this is followed by a period of enhanced (or at least anomalous) overturning lasting thousands of years until equilibrium interglacial conditions are attained and during which atmospheric CO2 is likely to decrease. Our results therefore suggest that deglacial oscillations in ocean circulation can have a lasting influence on early interglacial climate and highlight the transient nature of atmospheric CO2 overshoots associated with the onset of some previous interglacials. Accordingly, we suggest that these intervals should be considered as a part of the deglacial process. This has implications for studies concerned with the evolution of atmospheric CO2 during interglacial periods including the Holocene.

Bayesian Data Analysis for Revealing Causes of the Middle Pleistocene Transition

Currently, causes of the middle Pleistocene transition (MPT) – the onset of large-amplitude glacial variability with 100 kyr time scale instead of regular 41 kyr cycles before – are a challenging puzzle in Paleoclimatology. Here we show how a Bayesian data analysis based on machine learning approaches can help to reveal the main mechanisms underlying the Pleistocene variability, which most likely explain proxy records and can be used for testing existing theories. We construct a Bayesian data-driven model from benthic δ18O records (LR04 stack) accounting for the main factors which may potentially impact climate of the Pleistocene: internal climate dynamics, gradual trends, variations of insolation, and millennial variability. In contrast to some theories, we uncover that under long-term trends in climate, the strong glacial cycles have appeared due to internal nonlinear oscillations induced by millennial noise. We find that while the orbital Milankovitch forcing does not matter for the MPT onset, the obliquity oscillation phase-locks the climate cycles through the meridional gradient of insolation.

ЭМПИРИЧЕСКОЕ ИССЛЕДОВАНИЕ КРИТИЧЕСКОГО ПЕРЕХОДА В КЛИМАТЕ ПЛЕЙСТОЦЕНА НА ОСНОВЕ НЕЛИНЕЙНОДИНАМИЧЕСКОЙ РЕКОНСТРУКЦИИ

Until now, cause of the Mid-Pleistocene Transition (MPT), when the dominant periodicity of climate cycles changed from 41,000 to 100,000 years in the absence of significant change in orbital forcing, are still an open question in Paleoclimatology. Here we show how a Bayesian data analysis and nonlinear dynamical reconstruction methods can help to reveal the main mechanisms underlying the Pleistocene variability. Our Bayesian data-driven model from benthic d18O records (LR04 stack) accounts the main factors which may potentially impact climate of the Pleistocene: internal climate dynamics, gradual trends, variations of insolation, and millennial variability. In contrast to some theories, we uncover that under long-term trends in climate, the strong glacial cycles have appeared due to internal nonlinear oscillations induced by millennial noise. We find that while the orbital Milankovitch forcing does not matter for the MPT onset, the obliquity oscillation phase-locks the climate cycles through the meridional gradient of insolation. The research was supported by the RAS Presidium Program «Nonlinear dynamics: fundamental problems and applications».

The Mid-Pleistocene Transition induced by delayed feedback and bistability

The Mid-Pleistocene Transition, the shift from 41 kyr to 100 kyr glacial-interglacial cycles that occurred roughly 1 Myr ago, is often considered as a change in internal climate dynamics. Here we revisit the model of Quaternary climate dynamics that was proposed by Saltzman and Maasch (1988). We show that it is quantitatively similar to a scalar equation for the ice dynamics only when combining the remaining components into a single delayed feedback term. The delay is the sum of the internal times scales of ocean transport and ice sheet dynamics, which is on the order of 10 kyr. We find that, in the absence of astronomical forcing, the delayed feedback leads to bistable behaviour, where stable large-amplitude oscillations of ice volume and an equilibrium coexist over a large range of values for the delay. We then apply astronomical forcing. We perform a systematic study to show how the system response depends on the forcing amplitude. We find that over a wide range of forcing amplitudes the forcing leads to a switch from small-scale oscillations of 41 kyr to large-amplitude oscillations of roughly 100 kyr without any change of other parameters. The transition in the forced model consistently occurs near the time of the Mid-Pleistocene Transition as observed in data records. This provides evidence that the MPT could have been primarily a forcing-induced switch between attractors of the internal dynamics. Small additional random disturbances make the forcing-induced transition near 800 kyr BP even more robust. We also find that the forced system forgets its initial history during the small-scale oscillations, in particular, nearby initial conditions converge prior to transitioning. In contrast to this, in the regime of large-amplitude oscillations, the oscillation phase is very sensitive to random perturbations, which has a strong effect on the timing of the deglaciation events.

Modulation of hydroxyl variability by ENSO in the absence of external forcing

The hydroxyl radical (OH) is the primary oxidant in the troposphere, and the impact of its fluctuations on the methane budget has been disputed in recent years, however measurements of OH are insufficient to characterize global interannual fluctuations relevant for methane. Here, we use a 6,000-y control simulation of preindustrial conditions with a chemistry-climate model to quantify the natural variability in OH and internal feedbacks governing that variability. We find that, even in the absence of external forcing, maximum OH changes are 3.8 ± 0.8% over a decade, which is large in the context of the recent methane growth from 2007-2017. We show that the OH variability is not a white-noise process. A wavelet analysis indicates that OH variability exhibits significant feedbacks with the same periodicity as the El Niño-Southern Oscillation (ENSO). We find intrinsically generated modulation of the OH variability, suggesting that OH may show periods of rapid or no change in future decades that are solely due to the internal climate dynamics (as opposed to external forcings). An empirical orthogonal function analysis further indicates that ENSO is the dominant mode of OH variability, with the modulation of OH occurring primarily through lightning [Formula: see text] La Niña is associated with an increase in convection in the Tropical Pacific, which increases the simulated occurrence of lightning and allows for more OH production. Understanding this link between OH and ENSO may improve the predictability of the oxidative capacity of the troposphere and assist in elucidating the causes of current and historical trends in methane.

Modelling Climate and Societal Resilience in the Eastern Mediterranean in the Last Millennium

This article analyses high-quality hydroclimate proxy records and spatial reconstructions from the Central and Eastern Mediterranean and compares them with two Earth System Model simulations (CCSM4, MPI-ESM-P) for the Crusader period in the Levant (1095–1290 CE), the Mamluk regime in Transjordan (1260–1516 CE) and the Ottoman crisis and Celâlî Rebellion (1580–1610 CE). During the three time intervals, environmental and climatic stress tested the resilience of complex societies. We find that the multidecadal precipitation and drought variations in the Central and Eastern Mediterranean cannot be explained by external forcings (solar variations, tropical volcanism); rather they were driven by internal climate dynamics. Our research emphasises the challenges, opportunities and limitations of linking proxy records, palaeoreconstructions and model simulations to better understand how climate can affect human history.

Assessing reconstruction techniques of the Atlantic Ocean circulation variability during the last millennium

We assess the use of the meridional thermal-wind transport estimated from zonal density gradients to reconstruct the oceanic circulation variability during the last millennium in a forced simulation with the ECHO-G coupled climate model. Following a perfect-model approach, model-based pseudo-reconstructions of the Atlantic meridional overturning circulation (AMOC) and the Florida Current volume transport (FCT) are evaluated against their true simulated variability. The pseudo-FCT is additionally verified as proxy for AMOC strength and compared with the available proxy-based reconstruction. The thermal-wind component reproduces most of the simulated AMOC variability, which is mostly driven by internal climate dynamics during the preindustrial period and by increasing greenhouse gases afterwards. The pseudo-reconstructed FCT reproduces well the simulated FCT and reasonably well the variability of the AMOC strength, including the response to external forcing. The pseudo-reconstructed FCT, however, underestimates/overestimates the simulated variability at deep/shallow levels. Density changes responsible for the pseudo-reconstructed FCT are mainly driven by zonal temperature differences; salinity differences oppose but play a minor role. These results thus support the use of the thermal-wind relationship to reconstruct the oceanic circulation past variability, in particular at multidecadal timescales. Yet model-data comparison highlights important differences between the simulated and the proxy-based FCT variability. ECHO-G simulates a prominent weakening in the North Atlantic circulation that contrasts with the reconstructed enhancement. Our model results thus do not support the reconstructed FC minimum during the Little Ice Age. This points to a failure in the reconstruction, misrepresented processes in the model, or an important role of internal ocean dynamics.

Testing competing forms of the Milankovitch hypothesis: A multivariate approach

AbstractWe test competing forms of the Milankovitch hypothesis by estimating the coefficients and diagnostic statistics for a cointegrated vector autoregressive model that includes 10 climate variables and four exogenous variables for solar insolation. The estimates are consistent with the physical mechanisms postulated to drive glacial cycles. They show that the climate variables are driven partly by solar insolation, determining the timing and magnitude of glaciations and terminations, and partly by internal feedback dynamics, pushing the climate variables away from equilibrium. We argue that the latter is consistent with a weak form of the Milankovitch hypothesis and that it should be restated as follows: internal climate dynamics impose perturbations on glacial cycles that are driven by solar insolation. Our results show that these perturbations are likely caused by slow adjustment between land ice volume and solar insolation. The estimated adjustment dynamics show that solar insolation affects an array of climate variables other than ice volume, each at a unique rate. This implies that previous efforts to test the strong form of the Milankovitch hypothesis by examining the relationship between solar insolation and a single climate variable are likely to suffer from omitted variable bias.

A Pacific Centennial Oscillation Predicted by Coupled GCMs*

Abstract Internal climate variability at the centennial time scale is investigated using long control integrations from three state-of-the-art global coupled general circulation models. In the absence of external forcing, all three models produce centennial variability in the mean zonal sea surface temperature (SST) and sea level pressure (SLP) gradients in the equatorial Pacific with counterparts in the extratropics. The centennial pattern in the tropical Pacific is dissimilar to that of the interannual El Niño–Southern Oscillation (ENSO), in that the most prominent expression in temperature is found beneath the surface of the western Pacific warm pool. Some global repercussions nevertheless are analogous, such as a hemispherically symmetric atmospheric wave pattern of alternating highs and lows. Centennial variability in western equatorial Pacific SST is a result of the strong asymmetry of interannual ocean heat content anomalies, while the eastern equatorial Pacific exhibits a lagged, Bjerknes-like response to temperature and convection in the west. The extratropical counterpart is shown to be a flux-driven response to the hemispherically symmetric circulation anomalies emanating from the tropical Pacific. Significant centennial-length trends in the zonal SST and SLP gradients rivaling those estimated from observations and model simulations forced with increasing CO2 appear to be inherent features of the internal climate dynamics simulated by all three models. Unforced variability and trends on the centennial time scale therefore need to be addressed in estimated uncertainties, beyond more traditional signal-to-noise estimates that do not account for natural variability on the centennial time scale.

Oscillators and relaxation phenomena in Pleistocene climate theory

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.

Surface changes in the North Atlantic meridional overturning circulation during the last millennium

Despite numerous investigations, the dynamical origins of the Medieval Climate Anomaly and the Little Ice Age remain uncertain. A major unresolved issue relating to internal climate dynamics is the mode and tempo of Atlantic meridional overturning circulation variability, and the significance of decadal-to-centennial scale changes in Atlantic meridional overturning circulation strength in regulating the climate of the last millennium. Here we use the time-constrained high-resolution local radiocarbon reservoir age offset derived from an absolutely dated annually resolved shell chronology spanning the past 1,350 years, to reconstruct changes in surface ocean circulation and climate. The water mass tracer data presented here from the North Icelandic shelf, combined with previously published data from the Arctic and subtropical Atlantic, show that surface Atlantic meridional overturning circulation dynamics likely amplified the relatively warm conditions during the Medieval Climate Anomaly and the relatively cool conditions during the Little Ice Age within the North Atlantic sector.

El Niño–like pattern in the Pacific during marine isotope stages (MIS) 13 and 11?

Paleoceanographic evidence from midlatitudes bearing on the long‐term development of Pacific Ocean surface circulation through the middle and late Pleistocene is scant because of a lack of marine paleorecords exceeding the last four or five glacial cycles. Here we compare benthic and planktonic foraminiferal stable carbon isotope data from a 1‐Myr marine record in the southeastern subtropical Pacific with records from the equatorial Pacific and Southern Ocean in order to reconstruct sea surface circulation in this part of the world ocean and its possible link to the El Niño/Southern Oscillation (ENSO). We especially address marine isotope stages (MIS) 13 and 11, when internal climate dynamics are suggested to produce somewhat unusual interglacial periods. Our results show that the Southern Ocean controlled the circulation in the subtropical Pacific throughout the middle and late Pleistocene and that the hitherto existing hypotheses cannot explain the peak δ13C values in both MIS 13 and 11. We argue that the surface circulation pattern in the southeast Pacific during MIS 13 and 11 should have differed from any other interglacial, and we suggest that the subsequent changes in the latitudinal and meridional heat and moisture transfer contributed to the unusual paleoceanographic settings during these stages. We further hypothesize unusually strong El Niño–like conditions during MIS 13 and 11. This hypothesis, if true, challenges a direct forcing of ENSO by orbital parameters.