High-latitude biomes and rock weathering mediate climate-carbon cycle feedbacks on eccentricity timescales.

The International Ocean Discovery Programme (IODP) and its predecessors generated a treasure trove of Cenozoic climate and carbon cycle dynamics. Yet, it remains unclear how climate and carbon cycle interacted under changing geologic boundary conditions. Here, we present the carbon isotope (δ13C) megasplice, documenting deep-ocean δ13C evolution since 35 million years ago (Ma). We juxtapose the δ13C megasplice with its δ18O counterpart and determine their phase-difference on ~100-kyr eccentricity timescales. This analysis reveals that 2.4-Myr eccentricity cycles modulate the δ13C-δ18O phase relationship throughout the Oligo-Miocene (34-6 Ma), potentially through changes in continental weathering. At 6 Ma, a striking switch from in-phase to anti-phase behaviour occurs, signalling a reorganization of the climate-carbon cycle system. We hypothesize that this transition is consistent with Arctic cooling: Prior to 6 Ma, low-latitude continental carbon reservoirs expanded during astronomically-forced cool spells. After 6 Ma, however, continental carbon reservoirs contract rather than expand during cold periods due to competing effects between Arctic biomes (ice, tundra, taiga). We conclude that, on geologic timescales, System Earth experienced state-dependent modes of climate-carbon cycle interaction.

An astronomically dated record of Earth's climate and its predictability over the last 66 million years.

Much of our understanding of Earth's past climate comes from the measurement of oxygen and carbon isotope variations in deep-sea benthic foraminifera. Yet, long intervals in existing records lack the temporal resolution and age control needed to thoroughly categorize climate states of the Cenozoic era and to study their dynamics. Here, we present a new, highly resolved, astronomically dated, continuous composite of benthic foraminifer isotope records developed in our laboratories. Four climate states-Hothouse, Warmhouse, Coolhouse, Icehouse-are identified on the basis of their distinctive response to astronomical forcing depending on greenhouse gas concentrations and polar ice sheet volume. Statistical analysis of the nonlinear behavior encoded in our record reveals the key role that polar ice volume plays in the predictability of Cenozoic climate dynamics.

Temperate rainforests near the South Pole during peak Cretaceous warmth.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years1-5, driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume6. In the near absence of proximal geological records from south of the Antarctic Circle, it is disputed whether polar ice could exist under such environmental conditions. Here we use a sedimentary sequence recovered from the West Antarctic shelf-the southernmost Cretaceous record reported so far-and show that a temperate lowland rainforest environment existed at a palaeolatitude of about 82° S during the Turonian-Santonian age (92 to 83 million years ago). This record contains an intact 3-metre-long network of in situ fossil roots embedded in a mudstone matrix containing diverse pollen and spores. A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric carbon dioxide concentrations of 1,120-1,680 parts per million by volume and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo under high levels of atmospheric carbon dioxide.

Causes of ice age intensification across the Mid-Pleistocene Transition.

During the Mid-Pleistocene Transition (MPT; 1,200-800 kya), Earth's orbitally paced ice age cycles intensified, lengthened from ∼40,000 (∼40 ky) to ∼100 ky, and became distinctly asymmetrical. Testing hypotheses that implicate changing atmospheric CO2 levels as a driver of the MPT has proven difficult with available observations. Here, we use orbitally resolved, boron isotope CO2 data to show that the glacial to interglacial CO2 difference increased from ∼43 to ∼75 µatm across the MPT, mainly because of lower glacial CO2 levels. Through carbon cycle modeling, we attribute this decline primarily to the initiation of substantive dust-borne iron fertilization of the Southern Ocean during peak glacial stages. We also observe a twofold steepening of the relationship between sea level and CO2-related climate forcing that is suggestive of a change in the dynamics that govern ice sheet stability, such as that expected from the removal of subglacial regolith or interhemispheric ice sheet phase-locking. We argue that neither ice sheet dynamics nor CO2 change in isolation can explain the MPT. Instead, we infer that the MPT was initiated by a change in ice sheet dynamics and that longer and deeper post-MPT ice ages were sustained by carbon cycle feedbacks related to dust fertilization of the Southern Ocean as a consequence of larger ice sheets.

Very large release of mostly volcanic carbon during the Palaeocene-Eocene Thermal Maximum.

The Palaeocene-Eocene Thermal Maximum (PETM) was a global warming event that occurred about 56 million years ago, and is commonly thought to have been driven primarily by the destabilization of carbon from surface sedimentary reservoirs such as methane hydrates. However, it remains controversial whether such reservoirs were indeed the source of the carbon that drove the warming. Resolving this issue is key to understanding the proximal cause of the warming, and to quantifying the roles of triggers versus feedbacks. Here we present boron isotope data-a proxy for seawater pH-that show that the ocean surface pH was persistently low during the PETM. We combine our pH data with a paired carbon isotope record in an Earth system model in order to reconstruct the unfolding carbon-cycle dynamics during the event. We find strong evidence for a much larger (more than 10,000 petagrams)-and, on average, isotopically heavier-carbon source than considered previously. This leads us to identify volcanism associated with the North Atlantic Igneous Province, rather than carbon from a surface reservoir, as the main driver of the PETM. This finding implies that climate-driven amplification of organic carbon feedbacks probably played only a minor part in driving the event. However, we find that enhanced burial of organic matter seems to have been important in eventually sequestering the released carbon and accelerating the recovery of the Earth system.

Evolution of the early Antarctic ice ages.

Understanding the stability of the early Antarctic ice cap in the geological past is of societal interest because present-day atmospheric CO2 concentrations have reached values comparable to those estimated for the Oligocene and the Early Miocene epochs. Here we analyze a new high-resolution deep-sea oxygen isotope (δ18O) record from the South Atlantic Ocean spanning an interval between 30.1 My and 17.1 My ago. The record displays major oscillations in deep-sea temperature and Antarctic ice volume in response to the ∼110-ky eccentricity modulation of precession. Conservative minimum ice volume estimates show that waxing and waning of at least ∼85 to 110% of the volume of the present East Antarctic Ice Sheet is required to explain many of the ∼110-ky cycles. Antarctic ice sheets were typically largest during repeated glacial cycles of the mid-Oligocene (∼28.0 My to ∼26.3 My ago) and across the Oligocene-Miocene Transition (∼23.0 My ago). However, the high-amplitude glacial-interglacial cycles of the mid-Oligocene are highly symmetrical, indicating a more direct response to eccentricity modulation of precession than their Early Miocene counterparts, which are distinctly asymmetrical-indicative of prolonged ice buildup and delayed, but rapid, glacial terminations. We hypothesize that the long-term transition to a warmer climate state with sawtooth-shaped glacial cycles in the Early Miocene was brought about by subsidence and glacial erosion in West Antarctica during the Late Oligocene and/or a change in the variability of atmospheric CO2 levels on astronomical time scales that is not yet captured in existing proxy reconstructions.

Lessons on Climate Sensitivity From Past Climate Changes.

Over the last decade, our understanding of climate sensitivity has improved considerably. The climate system shows variability on many timescales, is subject to non-stationary forcing and it is most likely out of equilibrium with the changes in the radiative forcing. Slow and fast feedbacks complicate the interpretation of geological records as feedback strengths vary over time. In the geological past, the forcing timescales were different than at present, suggesting that the response may have behaved differently. Do these insights constrain the climate sensitivity relevant for the present day? In this paper, we review the progress made in theoretical understanding of climate sensitivity and on the estimation of climate sensitivity from proxy records. Particular focus lies on the background state dependence of feedback processes and on the impact of tipping points on the climate system. We suggest how to further use palaeo data to advance our understanding of the currently ongoing climate change.

A Cenozoic record of the equatorial Pacific carbonate compensation depth.

Atmospheric carbon dioxide concentrations and climate are regulated on geological timescales by the balance between carbon input from volcanic and metamorphic outgassing and its removal by weathering feedbacks; these feedbacks involve the erosion of silicate rocks and organic-carbon-bearing rocks. The integrated effect of these processes is reflected in the calcium carbonate compensation depth, which is the oceanic depth at which calcium carbonate is dissolved. Here we present a carbonate accumulation record that covers the past 53 million years from a depth transect in the equatorial Pacific Ocean. The carbonate compensation depth tracks long-term ocean cooling, deepening from 3.0-3.5 kilometres during the early Cenozoic (approximately 55 million years ago) to 4.6 kilometres at present, consistent with an overall Cenozoic increase in weathering. We find large superimposed fluctuations in carbonate compensation depth during the middle and late Eocene. Using Earth system models, we identify changes in weathering and the mode of organic-carbon delivery as two key processes to explain these large-scale Eocene fluctuations of the carbonate compensation depth.

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.

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.

The heartbeat of the Oligocene climate system.

A 13-million-year continuous record of Oligocene climate from the equatorial Pacific reveals a pronounced "heartbeat" in the global carbon cycle and periodicity of glaciations. This heartbeat consists of 405,000-, 127,000-, and 96,000-year eccentricity cycles and 1.2-million-year obliquity cycles in periodically recurring glacial and carbon cycle events. That climate system response to intricate orbital variations suggests a fundamental interaction of the carbon cycle, solar forcing, and glacial events. Box modeling shows that the interaction of the carbon cycle and solar forcing modulates deep ocean acidity as well as the production and burial of global biomass. The pronounced 405,000-year eccentricity cycle is amplified by the long residence time of carbon in the oceans.

The Cenozoic palaeoenvironment of the Arctic Ocean.

The history of the Arctic Ocean during the Cenozoic era (0-65 million years ago) is largely unknown from direct evidence. Here we present a Cenozoic palaeoceanographic record constructed from >400 m of sediment core from a recent drilling expedition to the Lomonosov ridge in the Arctic Ocean. Our record shows a palaeoenvironmental transition from a warm 'greenhouse' world, during the late Palaeocene and early Eocene epochs, to a colder 'icehouse' world influenced by sea ice and icebergs from the middle Eocene epoch to the present. For the most recent approximately 14 Myr, we find sedimentation rates of 1-2 cm per thousand years, in stark contrast to the substantially lower rates proposed in earlier studies; this record of the Neogene reveals cooling of the Arctic that was synchronous with the expansion of Greenland ice (approximately 3.2 Myr ago) and East Antarctic ice (approximately 14 Myr ago). We find evidence for the first occurrence of ice-rafted debris in the middle Eocene epoch (approximately 45 Myr ago), some 35 Myr earlier than previously thought; fresh surface waters were present at approximately 49 Myr ago, before the onset of ice-rafted debris. Also, the temperatures of surface waters during the Palaeocene/Eocene thermal maximum (approximately 55 Myr ago) appear to have been substantially warmer than previously estimated. The revised timing of the earliest Arctic cooling events coincides with those from Antarctica, supporting arguments for bipolar symmetry in climate change.

Centennial-scale climate cooling with a sudden cold event around 8,200 years ago.

The extent of climate variability during the current interglacial period, the Holocene, is still debated. Temperature records derived from central Greenland ice cores show one significant temperature anomaly between 8,200 and 8,100 years ago, which is often attributed to a meltwater outflow into the North Atlantic Ocean and a slowdown of North Atlantic Deep Water formation--this anomaly provides an opportunity to study such processes with relevance to present-day freshening of the North Atlantic. Anomalies in climate proxy records from locations around the globe are often correlated with this sharp event in Greenland. But the anomalies in many of these records span 400 to 600 years, start from about 8,600 years ago and form part of a repeating pattern within the Holocene. More sudden climate changes around 8,200 years ago appear superimposed on this longer-term cooling. The compounded nature of the signals implies that far-field climate anomalies around 8,200 years ago cannot be used in a straightforward manner to assess the impact of a slowdown of North Atlantic Deep Water formation, and the geographical extent of the rapid cooling event 8,200 years ago remains to be determined.

Rapid stepwise onset of Antarctic glaciation and deeper calcite compensation in the Pacific Ocean.

The ocean depth at which the rate of calcium carbonate input from surface waters equals the rate of dissolution is termed the calcite compensation depth. At present, this depth is approximately 4,500 m, with some variation between and within ocean basins. The calcite compensation depth is linked to ocean acidity, which is in turn linked to atmospheric carbon dioxide concentrations and hence global climate. Geological records of changes in the calcite compensation depth show a prominent deepening of more than 1 km near the Eocene/Oligocene boundary (approximately 34 million years ago) when significant permanent ice sheets first appeared on Antarctica, but the relationship between these two events is poorly understood. Here we present ocean sediment records of calcium carbonate content as well as carbon and oxygen isotopic compositions from the tropical Pacific Ocean that cover the Eocene/Oligocene boundary. We find that the deepening of the calcite compensation depth was more rapid than previously documented and occurred in two jumps of about 40,000 years each, synchronous with the stepwise onset of Antarctic ice-sheet growth. The glaciation was initiated, after climatic preconditioning, by an interval when the Earth's orbit of the Sun favoured cool summers. The changes in oxygen-isotope composition across the Eocene/Oligocene boundary are too large to be explained by Antarctic ice-sheet growth alone and must therefore also indicate contemporaneous global cooling and/or Northern Hemisphere glaciation.