Early onset and tropical forcing of 100,000-year Pleistocene glacial cycles

Between 1.5 and 0.6 Myr ago, the period of the Earth's glacial cycles changed from 41 kyr, the period of the Earth's obliquity cycles, to 100 kyr, the period of the Earth's orbital eccentricity, which has a much smaller effect on global insolation. The timing of this transition and its causes pose one of the most perplexing problems in palaeoclimate research. Here we use complex demodulation to examine the phase evolution of precession and semiprecession cycles--the latter of which are phase-coupled to both precession and eccentricity--in the tropical and extratropical Atlantic Ocean. We find that about 1.5 Myr ago, tropical semiprecession cycles (with periods of about 11.5 kyr) started to propagate to higher latitudes, coincident with a growing amplitude envelope of the 100-kyr cycles. Evidence from numerical models suggests that cycles of about 10 kyr in length may be required to explain the high amplitude of the 100-kyr cycles. Combining our results with consideration of a modern analogue, we conclude that increased heat flow across the equator or from the tropics to higher latitudes around 1.5 Myr ago strengthened the semiprecession cycle in the Northern Hemisphere, and triggered the transition to sustained 100-kyr glacial cycles.     

Atlantic overturning responses to Late Pleistocene climate forcings

The factors driving glacial changes in ocean overturning circulation are not well understood. On the basis of a comparison of 20 climate variables over the past four glacial cycles, the SPECMAP project proposed that summer insolation at high northern latitudes (that is, Milankovitch forcing) drives the same sequence of ocean circulation and other climate responses over 100-kyr eccentricity cycles, 41-kyr obliquity cycles and 23-kyr precession cycles. SPECMAP analysed the circulation response at only a few sites in the Atlantic Ocean, however, and the phase of circulation response has been shown to vary by site and orbital band. Here we test the SPECMAP hypothesis by measuring the phase of orbital responses in benthic delta(13)C (a proxy indicator of ocean nutrient content) at 24 sites throughout the Atlantic over the past 425 kyr. On the basis of delta(13)C responses at 3,000-4,010 m water depth, we find that maxima in Milankovitch forcing are associated with greater mid-depth overturning in the obliquity band but less overturning in the precession band. This suggests that Atlantic overturning is strongly sensitive to factors beyond ice volume and summer insolation at high northern latitudes. A better understanding of these processes could lead to improvements in model estimates of overturning rates, which range from a 40 per cent increase to a 40 per cent decrease at the Last Glacial Maximum and a 10-50 per cent decrease over the next 140 yr in response to projected increases in atmospheric CO(2) (ref. 4).     

High-latitude influence on the eastern equatorial Pacific climate in the early Pleistocene epoch

Many records of tropical sea surface temperature and marine productivity exhibit cycles of 23 kyr (orbital precession) and 100 kyr during the past 0.5 Myr (refs 1-5), whereas high-latitude sea surface temperature records display much more pronounced obliquity cycles at a period of about 41 kyr (ref. 6). Little is known, however, about tropical climate variability before the mid-Pleistocene transition about 900 kyr ago, which marks the change from a climate dominated by 41-kyr cycles (when ice-age cycles and high-latitude sea surface temperature variations were dictated by changes in the Earth's obliquity) to the more recent 100-kyr cycles of ice ages. Here we analyse alkenones from marine sediments in the eastern equatorial Pacific Ocean to reconstruct sea surface temperatures and marine productivity over the past 1.8 Myr. We find that both records are dominated by the 41-kyr obliquity cycles between 1.8 and 1.2 Myr ago, with a relatively small contribution from orbital precession, and that early Pleistocene sea surface temperatures varied in the opposite sense to local annual insolation in the eastern equatorial Pacific Ocean. We conclude that during the early Pleistocene epoch, climate variability at our study site must have been determined by high-latitude processes that were driven by orbital obliquity forcing.     

Combined obliquity and precession pacing of late Pleistocene deglaciations

Milankovitch proposed that Earth resides in an interglacial state when its spin axis both tilts to a high obliquity and precesses to align the Northern Hemisphere summer with Earth's nearest approach to the Sun. This general concept has been elaborated into hypotheses that precession, obliquity or combinations of both could pace deglaciations during the late Pleistocene. Earlier tests have shown that obliquity paces the late Pleistocene glacial cycles but have been inconclusive with regard to precession, whose shorter period of about 20,000 years makes phasing more sensitive to timing errors. No quantitative test has provided firm evidence for a dual effect. Here I show that both obliquity and precession pace late Pleistocene glacial cycles. Deficiencies in time control that have long stymied efforts to establish orbital effects on deglaciation are overcome using a new statistical test that focuses on maxima in orbital forcing. The results are fully consistent with Milankovitch's proposal but also admit the possibility that long Southern Hemisphere summers contribute to deglaciation.     

Proterozoic low orbital obliquity and axial-dipolar geomagnetic field from evaporite palaeolatitudes

Palaeomagnetism of climatically sensitive sedimentary rock types, such as glacial deposits and evaporites, can test the uniformitarianism of ancient geomagnetic fields and palaeoclimate zones. Proterozoic glacial deposits laid down in near-equatorial palaeomagnetic latitudes can be explained by 'snowball Earth' episodes, high orbital obliquity or markedly non-uniformitarian geomagnetic fields. Here I present a global palaeomagnetic compilation of the Earth's entire basin-scale evaporite record. Magnetic inclinations are consistent with low orbital obliquity and a geocentric-axial-dipole magnetic field for most of the past two billion years, and the snowball Earth hypothesis accordingly remains the most viable model for low-latitude Proterozoic ice ages. Efforts to reconstruct Proterozoic supercontinents are strengthened by this demonstration of a consistently axial and dipolar geomagnetic reference frame, which itself implies stability of geodynamo processes on billion-year timescales.     

Orbital and millennial-scale features of atmospheric CH4 over the past 800,000 years

Atmospheric methane is an important greenhouse gas and a sensitive indicator of climate change and millennial-scale temperature variability. Its concentrations over the past 650,000 years have varied between approximately 350 and approximately 800 parts per 10(9) by volume (p.p.b.v.) during glacial and interglacial periods, respectively. In comparison, present-day methane levels of approximately 1,770 p.p.b.v. have been reported. Insights into the external forcing factors and internal feedbacks controlling atmospheric methane are essential for predicting the methane budget in a warmer world. Here we present a detailed atmospheric methane record from the EPICA Dome C ice core that extends the history of this greenhouse gas to 800,000 yr before present. The average time resolution of the new data is approximately 380 yr and permits the identification of orbital and millennial-scale features. Spectral analyses indicate that the long-term variability in atmospheric methane levels is dominated by approximately 100,000 yr glacial-interglacial cycles up to approximately 400,000 yr ago with an increasing contribution of the precessional component during the four more recent climatic cycles. We suggest that changes in the strength of tropical methane sources and sinks (wetlands, atmospheric oxidation), possibly influenced by changes in monsoon systems and the position of the intertropical convergence zone, controlled the atmospheric methane budget, with an additional source input during major terminations as the retreat of the northern ice sheet allowed higher methane emissions from extending periglacial wetlands. Millennial-scale changes in methane levels identified in our record as being associated with Antarctic isotope maxima events are indicative of ubiquitous millennial-scale temperature variability during the past eight glacial cycles.     

Orbitally induced oscillations in the East Antarctic ice sheet at the Oligocene/Miocene boundary

Between 34 and 15 million years (Myr) ago, when planetary temperatures were 3-4 degrees C warmer than at present and atmospheric CO2 concentrations were twice as high as today, the Antarctic ice sheets may have been unstable. Oxygen isotope records from deep-sea sediment cores suggest that during this time fluctuations in global temperatures and high-latitude continental ice volumes were influenced by orbital cycles. But it has hitherto not been possible to calibrate the inferred changes in ice volume with direct evidence for oscillations of the Antarctic ice sheets. Here we present sediment data from shallow marine cores in the western Ross Sea that exhibit well dated cyclic variations, and which link the extent of the East Antarctic ice sheet directly to orbital cycles during the Oligocene/Miocene transition (24.1-23.7 Myr ago). Three rapidly deposited glacimarine sequences are constrained to a period of less than 450 kyr by our age model, suggesting that orbital influences at the frequencies of obliquity (40 kyr) and eccentricity (125 kyr) controlled the oscillations of the ice margin at that time. An erosional hiatus covering 250 kyr provides direct evidence for a major episode of global cooling and ice-sheet expansion about 23.7 Myr ago, which had previously been inferred from oxygen isotope data (Mi1 event).     

Long-period astronomical forcing of mammal turnover

Mammals are among the fastest-radiating groups, being characterized by a mean species lifespan of the order of 2.5 million years (Myr). The basis for this characteristic timescale of origination, extinction and turnover is not well understood. Various studies have invoked climate change to explain mammalian species turnover, but other studies have either challenged or only partly confirmed the climate-turnover hypothesis. Here we use an exceptionally long (24.5-2.5 Myr ago), dense, and well-dated terrestrial record of rodent lineages from central Spain, and show the existence of turnover cycles with periods of 2.4-2.5 and 1.0 Myr. We link these cycles to low-frequency modulations of Milankovitch oscillations, and show that pulses of turnover occur at minima of the 2.37-Myr eccentricity cycle and nodes of the 1.2-Myr obliquity cycle. Because obliquity nodes and eccentricity minima are associated with ice sheet expansion and cooling and affect regional precipitation, we infer that long-period astronomical climate forcing is a major determinant of species turnover in small mammals and probably other groups as well.     

Geological constraints on tidal dissipation and dynamical ellipticity of the Earth over the past three million years

The evolution of the Solar System has been shown to be chaotic, which limits our ability to retrace the orbital and precessional motion of the Earth over more than 35-50 Myr (ref. 2). Moreover, the precession, obliquity and insolation parameters can also be influenced by secular variations in the tidal dissipation and dynamical ellipticity of the Earth induced by glacial cyclicity and mantle convection. Here we determine the average values of these dissipative effects over the past three million years. We have computed the optimal fit between an exceptional palaeoclimate record from the eastern Mediterranean Sea and a model of the astronomical and insolation history by testing a number of values for the tidal dissipation and dynamical ellipticity parameters. We find that the combined effects of dynamical ellipticity and tidal dissipation were, on average, significantly lower over the past three million years, compared to their present-day values (determined from artificial satellite data and lunar ranging). This secular variation associated with the Plio-Pleistocene ice load history has caused an average acceleration in the Earth's rotation over the past 3 Myr, which needs to be considered in the construction of astronomical timescales and in research into the stationarity of phase relations in the ocean-climate system through time.     

Thresholds for Cenozoic bipolar glaciation

The long-standing view of Earth's Cenozoic glacial history calls for the first continental-scale glaciation of Antarctica in the earliest Oligocene epoch ( approximately 33.6 million years ago), followed by the onset of northern-hemispheric glacial cycles in the late Pliocene epoch, about 31 million years later. The pivotal early Oligocene event is characterized by a rapid shift of 1.5 parts per thousand in deep-sea benthic oxygen-isotope values (Oi-1) within a few hundred thousand years, reflecting a combination of terrestrial ice growth and deep-sea cooling. The apparent absence of contemporaneous cooling in deep-sea Mg/Ca records, however, has been argued to reflect the growth of more ice than can be accommodated on Antarctica; this, combined with new evidence of continental cooling and ice-rafted debris in the Northern Hemisphere during this period, raises the possibility that Oi-1 represents a precursory bipolar glaciation. Here we test this hypothesis using an isotope-capable global climate/ice-sheet model that accommodates both the long-term decline of Cenozoic atmospheric CO(2) levels and the effects of orbital forcing. We show that the CO(2) threshold below which glaciation occurs in the Northern Hemisphere ( approximately 280 p.p.m.v.) is much lower than that for Antarctica ( approximately 750 p.p.m.v.). Therefore, the growth of ice sheets in the Northern Hemisphere immediately following Antarctic glaciation would have required rapid CO(2) drawdown within the Oi-1 timeframe, to levels lower than those estimated by geochemical proxies and carbon-cycle models. Instead of bipolar glaciation, we find that Oi-1 is best explained by Antarctic glaciation alone, combined with deep-sea cooling of up to 4 degrees C and Antarctic ice that is less isotopically depleted (-30 to -35 per thousand) than previously suggested. Proxy CO(2) estimates remain above our model's northern-hemispheric glaciation threshold of approximately 280 p.p.m.v. until approximately 25 Myr ago, but have been near or below that level ever since. This implies that episodic northern-hemispheric ice sheets have been possible some 20 million years earlier than currently assumed (although still much later than Oi-1) and could explain some of the variability in Miocene sea-level records.     

Insolation-driven changes in atmospheric circulation over the past 116,000 years in subtropical Brazil

During the last glacial period, large millennial-scale temperature oscillations--the 'Dansgaard/Oeschger' cycles--were the primary climate signal in Northern Hemisphere climate archives from the high latitudes to the tropics. But whether the influence of these abrupt climate changes extended to the tropical and subtropical Southern Hemisphere, where changes in insolation are thought to be the main direct forcing of climate, has remained unclear. Here we present a high-resolution oxygen isotope record of a U/Th-dated stalagmite from subtropical southern Brazil, covering the past 116,200 years. The oxygen isotope signature varies with shifts in the source region and amount of rainfall in the area, and hence records changes in atmospheric circulation and convective intensity over South America. We find that these variations in rainfall source and amount are primarily driven by summer solar radiation, which is controlled by the Earth's precessional cycle. The Dansgaard/Oeschger cycles can be detected in our record and therefore we confirm that they also affect the tropical hydrological cycle, but that in southern subtropical Brazil, millennial-scale climate changes are not as dominant as they are in the Northern Hemisphere.     

Impacts of orbital forcing and atmospheric carbon dioxide on Miocene ice-sheet expansion

The processes causing the middle Miocene global cooling, which marked the Earth's final transition into an 'icehouse' climate about 13.9 million years ago (Myr ago), remain enigmatic. Tectonically driven circulation changes and variations in atmospheric carbon dioxide levels have been suggested as driving mechanisms, but the lack of adequately preserved sedimentary successions has made rigorous testing of these hypotheses difficult. Here we present high-resolution climate proxy records, covering the period from 14.7 to 12.7 million years ago, from two complete sediment cores from the northwest and southeast subtropical Pacific Ocean. Using new chronologies through the correlation to the latest orbital model, we find relatively constant, low summer insolation over Antarctica coincident with declining atmospheric carbon dioxide levels at the time of Antarctic ice-sheet expansion and global cooling, suggesting a causal link. We surmise that the thermal isolation of Antarctica played a role in providing sustained long-term climatic boundary conditions propitious for ice-sheet formation. Our data document that Antarctic glaciation was rapid, taking place within two obliquity cycles, and coincided with a striking transition from obliquity to eccentricity as the drivers of climatic change.     

Tropical climate changes at millennial and orbital timescales on the Bolivian Altiplano

Tropical South America is one of the three main centres of the global, zonal overturning circulation of the equatorial atmosphere (generally termed the 'Walker' circulation). Although this area plays a key role in global climate cycles, little is known about South American climate history. Here we describe sediment cores and down-hole logging results of deep drilling in the Salar de Uyuni, on the Bolivian Altiplano, located in the tropical Andes. We demonstrate that during the past 50,000 years the Altiplano underwent important changes in effective moisture at both orbital (20,000-year) and millennial timescales. Long-duration wet periods, such as the Last Glacial Maximum--marked in the drill core by continuous deposition of lacustrine sediments--appear to have occurred in phase with summer insolation maxima produced by the Earth's precessional cycle. Short-duration, millennial events correlate well with North Atlantic cold events, including Heinrich events 1 and 2, as well as the Younger Dryas episode. At both millennial and orbital timescales, cold sea surface temperatures in the high-latitude North Atlantic were coeval with wet conditions in tropical South America, suggesting a common forcing.     

Astronomical pacing of late Palaeocene to early Eocene global warming events

At the boundary between the Palaeocene and Eocene epochs, about 55 million years ago, the Earth experienced a strong global warming event, the Palaeocene-Eocene thermal maximum. The leading hypothesis to explain the extreme greenhouse conditions prevalent during this period is the dissociation of 1,400 to 2,800 gigatonnes of methane from ocean clathrates, resulting in a large negative carbon isotope excursion and severe carbonate dissolution in marine sediments. Possible triggering mechanisms for this event include crossing a threshold temperature as the Earth warmed gradually, comet impact, explosive volcanism or ocean current reorganization and erosion at continental slopes, whereas orbital forcing has been excluded. Here we report a distinct carbonate-poor red clay layer in deep-sea cores from Walvis ridge, which we term the Elmo horizon. Using orbital tuning, we estimate deposition of the Elmo horizon at about 2 million years after the Palaeocene-Eocene thermal maximum. The Elmo horizon has similar geochemical and biotic characteristics as the Palaeocene-Eocene thermal maximum, but of smaller magnitude. It is coincident with carbon isotope depletion events in other ocean basins, suggesting that it represents a second global thermal maximum. We show that both events correspond to maxima in the approximately 405-kyr and approximately 100-kyr eccentricity cycles that post-date prolonged minima in the 2.25-Myr eccentricity cycle, implying that they are indeed astronomically paced.     

Past extreme warming events linked to massive carbon release from thawing permafrost

Between about 55.5 and 52 million years ago, Earth experienced a series of sudden and extreme global warming events (hyperthermals) superimposed on a long-term warming trend. The first and largest of these events, the Palaeocene-Eocene Thermal Maximum (PETM), is characterized by a massive input of carbon, ocean acidification and an increase in global temperature of about 5 °C within a few thousand years. Although various explanations for the PETM have been proposed, a satisfactory model that accounts for the source, magnitude and timing of carbon release at the PETM and successive hyperthermals remains elusive. Here we use a new astronomically calibrated cyclostratigraphic record from central Italy to show that the Early Eocene hyperthermals occurred during orbits with a combination of high eccentricity and high obliquity. Corresponding climate-ecosystem-soil simulations accounting for rising concentrations of background greenhouse gases and orbital forcing show that the magnitude and timing of the PETM and subsequent hyperthermals can be explained by the orbitally triggered decomposition of soil organic carbon in circum-Arctic and Antarctic terrestrial permafrost. This massive carbon reservoir had the potential to repeatedly release thousands of petagrams (10(15) grams) of carbon to the atmosphere-ocean system, once a long-term warming threshold had been reached just before the PETM. Replenishment of permafrost soil carbon stocks following peak warming probably contributed to the rapid recovery from each event, while providing a sensitive carbon reservoir for the next hyperthermal. As background temperatures continued to rise following the PETM, the areal extent of permafrost steadily declined, resulting in an incrementally smaller available carbon pool and smaller hyperthermals at each successive orbital forcing maximum. A mechanism linking Earth's orbital properties with release of soil carbon from permafrost provides a unifying model accounting for the salient features of the hyperthermals.     

Muted climate variations in continental Siberia during the mid-Pleistocene epoch

The large difference in carbon and oxygen isotope data from the marine record between marine oxygen isotope stage 12 (MIS 12) and MIS 11, spanning the interval between about 480 and 380 kyr ago, has been interpreted as a transition between an extremely cold glacial period and an unusually warm interglacial period, with consequences for global ice volume, sea level and the global carbon cycle. The extent of the change is intriguing, because orbital forcing is predicted to have been relatively weak at that time. Here we analyse a continuous sediment record from Lake Baikal, Siberia, which reveals a virtually continuous interglacial diatom assemblage, a stable littoral benthic diatom assemblage and lithogenic sediments with 'interglacial' characteristics for the period from MIS 15a to MIS 11 (from about 580 to 380 kyr ago). From these data, we infer significantly weaker climate contrasts between MIS 12 and 11 than during more recent glacial-interglacial transitions in the late Pleistocene epoch (about 130 to 10 kyr ago). For the period from MIS 15a to MIS 11, we also infer an apparent lack of extensive mountain glaciation.     

North American ice-sheet dynamics and the onset of 100,000-year glacial cycles

The onset of major glaciations in the Northern Hemisphere about 2.7 million years ago was most probably induced by climate cooling during the late Pliocene epoch. These glaciations, during which the Northern Hemisphere ice sheets successively expanded and retreated, are superimposed on this long-term climate trend, and have been linked to variations in the Earth's orbital parameters. One intriguing problem associated with orbitally driven glacial cycles is the transition from 41,000-year to 100,000-year climatic cycles that occurred without an apparent change in insolation forcing. Several hypotheses have been proposed to explain the transition, both including and excluding ice-sheet dynamics. Difficulties in finding a conclusive answer to this palaeoclimatic problem are related to the lack of sufficiently long records of ice-sheet volume or sea level. Here we use a comprehensive ice-sheet model and a simple ocean-temperature model to extract three-million-year mutually consistent records of surface air temperature, ice volume and sea level from marine benthic oxygen isotopes. Although these records and their relative phasings are subject to considerable uncertainty owing to limited availability of palaeoclimate constraints, the results suggest that the gradual emergence of the 100,000-year cycles can be attributed to the increased ability of the merged North American ice sheets to survive insolation maxima and reach continental-scale size. The oversized, wet-based ice sheet probably responded to the subsequent insolation maximum by rapid thinning through increased basal-sliding, thereby initiating a glacial termination. Based on our assessment of the temporal changes in air temperature and ice volume during individual glacials, we demonstrate the importance of ice dynamics and ice-climate interactions in establishing the 100,000-year glacial cycles, with enhanced North American ice-sheet growth and the subsequent merging of the ice sheets being key elements.     

Possible solar origin of the 1,470-year glacial climate cycle demonstrated in a coupled model

Many palaeoclimate records from the North Atlantic region show a pattern of rapid climate oscillations, the so-called Dansgaard-Oeschger events, with a quasi-periodicity of approximately 1,470 years for the late glacial period. Various hypotheses have been suggested to explain these rapid temperature shifts, including internal oscillations in the climate system and external forcing, possibly from the Sun. But whereas pronounced solar cycles of approximately 87 and approximately 210 years are well known, a approximately 1,470-year solar cycle has not been detected. Here we show that an intermediate-complexity climate model with glacial climate conditions simulates rapid climate shifts similar to the Dansgaard-Oeschger events with a spacing of 1,470 years when forced by periodic freshwater input into the North Atlantic Ocean in cycles of approximately 87 and approximately 210 years. We attribute the robust 1,470-year response time to the superposition of the two shorter cycles, together with strongly nonlinear dynamics and the long characteristic timescale of the thermohaline circulation. For Holocene conditions, similar events do not occur. We conclude that the glacial 1,470-year climate cycles could have been triggered by solar forcing despite the absence of a 1,470-year solar cycle.     

Eocene global warming events driven by ventilation of oceanic dissolved organic carbon

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (～65-34 million years (Myr) ago). The most extreme hyperthermal was the ～170 thousand year (kyr) interval of 5-7 °C global warming during the Palaeocene-Eocene Thermal Maximum (PETM, 56?Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs, and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon. Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (～40?kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was re-sequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM. Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources, but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.     

Northern Hemisphere forcing of climatic cycles in Antarctica over the past 360,000 years

The Milankovitch theory of climate change proposes that glacial-interglacial cycles are driven by changes in summer insolation at high northern latitudes. The timing of climate change in the Southern Hemisphere at glacial-interglacial transitions (which are known as terminations) relative to variations in summer insolation in the Northern Hemisphere is an important test of this hypothesis. So far, it has only been possible to apply this test to the most recent termination, because the dating uncertainty associated with older terminations is too large to allow phase relationships to be determined. Here we present a new chronology of Antarctic climate change over the past 360,000 years that is based on the ratio of oxygen to nitrogen molecules in air trapped in the Dome Fuji and Vostok ice cores. This ratio is a proxy for local summer insolation, and thus allows the chronology to be constructed by orbital tuning without the need to assume a lag between a climate record and an orbital parameter. The accuracy of the chronology allows us to examine the phase relationships between climate records from the ice cores and changes in insolation. Our results indicate that orbital-scale Antarctic climate change lags Northern Hemisphere insolation by a few millennia, and that the increases in Antarctic temperature and atmospheric carbon dioxide concentration during the last four terminations occurred within the rising phase of Northern Hemisphere summer insolation. These results support the Milankovitch theory that Northern Hemisphere summer insolation triggered the last four deglaciations.