Atmospheric carbon dioxide concentrations over the past 60 million years

Knowledge of the evolution of atmospheric carbon dioxide concentrations throughout the Earth's history is important for a reconstruction of the links between climate and radiative forcing of the Earth's surface temperatures. Although atmospheric carbon dioxide concentrations in the early Cenozoic era (about 60 Myr ago) are widely believed to have been higher than at present, there is disagreement regarding the exact carbon dioxide levels, the timing of the decline and the mechanisms that are most important for the control of CO2 concentrations over geological timescales. Here we use the boron-isotope ratios of ancient planktonic foraminifer shells to estimate the pH of surface-layer sea water throughout the past 60 million years, which can be used to reconstruct atmospheric CO2 concentrations. We estimate CO2 concentrations of more than 2,000 p.p.m. for the late Palaeocene and earliest Eocene periods (from about 60 to 52 Myr ago), and find an erratic decline between 55 and 40 Myr ago that may have been caused by reduced CO2 outgassing from ocean ridges, volcanoes and metamorphic belts and increased carbon burial. Since the early Miocene (about 24 Myr ago), atmospheric CO2 concentrations appear to have remained below 500 p.p.m. and were more stable than before, although transient intervals of CO2 reduction may have occurred during periods of rapid cooling approximately 15 and 3 Myr ago.     

Warm tropical ocean surface and global anoxia during the mid-Cretaceous period

The middle of the Cretaceous period (about 120 to 80 Myr ago) was a time of unusually warm polar temperatures, repeated reef-drowning in the tropics and a series of oceanic anoxic events (OAEs) that promoted both the widespread deposition of organic-carbon-rich marine sediments and high biological turnover. The cause of the warm temperatures is unproven but widely attributed to high levels of atmospheric greenhouse gases such as carbon dioxide. In contrast, there is no consensus on the climatic causes and effects of the OAEs, with both high biological productivity and ocean 'stagnation' being invoked as the cause of ocean anoxia. Here we show, using stable isotope records from multiple species of well-preserved foraminifera, that the thermal structure of surface waters in the western tropical Atlantic Ocean underwent pronounced variability about 100 Myr ago, with maximum sea surface temperatures 3-5 degrees C warmer than today. This variability culminated in a collapse of upper-ocean stratification during OAE-1d (the 'Breistroffer' event), a globally significant period of organic-carbon burial that we show to have fundamental, stratigraphically valuable, geochemical similarities to the main OAEs of the Mesozoic era. Our records are consistent with greenhouse forcing being responsible for the warm temperatures, but are inconsistent both with explanations for OAEs based on ocean stagnation, and with the traditional view (reviewed in ref. 12) that past warm periods were more stable than today's climate.     

Late Miocene decoupling of oceanic warmth and atmospheric carbon dioxide forcing

Deep-time palaeoclimate studies are vitally important for developing a complete understanding of climate responses to changes in the atmospheric carbon dioxide concentration (that is, the atmospheric partial pressure of CO(2), p(co(2))). Although past studies have explored these responses during portions of the Cenozoic era (the most recent 65.5 million years (Myr) of Earth history), comparatively little is known about the climate of the late Miocene (～12-5 Myr ago), an interval with p(co(2)) values of only 200-350?parts per million by volume but nearly ice-free conditions in the Northern Hemisphere and warmer-than-modern temperatures on the continents. Here we present quantitative geochemical sea surface temperature estimates from the Miocene mid-latitude North Pacific Ocean, and show that oceanic warmth persisted throughout the interval of low p(co(2)) ～12-5 Myr ago. We also present new stable isotope measurements from the western equatorial Pacific that, in conjunction with previously published data, reveal a long-term trend of thermocline shoaling in the equatorial Pacific since ～13?Myr ago. We propose that a relatively deep global thermocline, reductions in low-latitude gradients in sea surface temperature, and cloud and water vapour feedbacks may help to explain the warmth of the late Miocene. Additional shoaling of the thermocline after 5?Myr ago probably explains the stronger coupling between p(co(2)), sea surface temperatures and climate that is characteristic of the more recent Pliocene and Pleistocene epochs.     

Evidence for decoupling of atmospheric CO2 and global climate during the Phanerozoic eon

Atmospheric carbon dioxide concentrations are believed to drive climate changes from glacial to interglacial modes, although geological and astronomical mechanisms have been invoked as ultimate causes. Additionally, it is unclear whether the changes between cold and warm modes should be regarded as a global phenomenon, affecting tropical and high-latitude temperatures alike, or if they are better described as an expansion and contraction of the latitudinal climate zones, keeping equatorial temperatures approximately constant. Here we present a reconstruction of tropical sea surface temperatures throughout the Phanerozoic eon (the past approximately 550 Myr) from our database of oxygen isotopes in calcite and aragonite shells. The data indicate large oscillations of tropical sea surface temperatures in phase with the cold-warm cycles, thus favouring the idea of climate variability as a global phenomenon. But our data conflict with a temperature reconstruction using an energy balance model that is forced by reconstructed atmospheric carbon dioxide concentrations. The results can be reconciled if atmospheric carbon dioxide concentrations were not the principal driver of climate variability on geological timescales for at least one-third of the Phanerozoic eon, or if the reconstructed carbon dioxide concentrations are not reliable.     

The influence of rivers on marine boron isotopes and implications for reconstructing past ocean pH

Ocean pH is particularly sensitive to atmospheric carbon dioxide content. Records of ocean pH can therefore be used to estimate past atmospheric carbon dioxide concentrations. The isotopic composition of boron (delta11B) contained in the carbonate shells of marine organisms varies according to pH, from which ocean pH can be reconstructed. This requires independent estimates of the delta11B of dissolved boron in sea water through time. The marine delta11B budget, however, is still largely unconstrained. Here we show that, by incorporating the global flux of riverine boron (as estimated from delta11B measurements in 22 of the world's main rivers), the marine boron isotope budget can be balanced. We also derive ocean delta11B budgets for the past 120 Myr. Estimated isotope compositions of boron in sea water show a remarkable consistency with records of delta11B in foraminiferal carbonates, suggesting that foraminifera delta11B records may in part reflect changes in the marine boron isotope budget rather than changes in ocean pH over the Cenozoic era.     

Termination of global warmth at the Palaeocene/Eocene boundary through productivity feedback

The onset of the Palaeocene/Eocene thermal maximum (about 55 Myr ago) was marked by global surface temperatures warming by 5-7 degrees C over approximately 30,000 yr (ref. 1), probably because of enhanced mantle outgassing and the pulsed release of approximately 1,500 gigatonnes of methane carbon from decomposing gas-hydrate reservoirs. The aftermath of this rapid, intense and global warming event may be the best example in the geological record of the response of the Earth to high atmospheric carbon dioxide concentrations and high temperatures. This response has been suggested to include an intensified flux of organic carbon from the ocean surface to the deep ocean and its subsequent burial through biogeochemical feedback mechanisms. Here we present firm evidence for this view from two ocean drilling cores, which record the largest accumulation rates of biogenic barium--indicative of export palaeoproductivity--at times of maximum global temperatures and peak excursion values of delta13C. The unusually rapid return of delta13C to values similar to those before the methane release and the apparent coupling of the accumulation rates of biogenic barium to temperature, suggests that the enhanced deposition of organic matter to the deep sea may have efficiently cooled this greenhouse climate by the rapid removal of excess carbon dioxide from the atmosphere.     

Stable sea surface temperatures in the western Pacific warm pool over the past 1.75 million years

About 850,000 years ago, the period of the glacial cycles changed from 41,000 to 100,000 years. This mid-Pleistocene climate transition has been attributed to global cooling, possibly caused by a decrease in atmospheric carbon dioxide concentrations. However, evidence for such cooling is currently restricted to the cool upwelling regions in the eastern equatorial oceans, although the tropical warm pools on the western side of the ocean basins are particularly sensitive to changes in radiative forcing. Here we present high-resolution records of sea surface temperatures spanning the past 1.75 million years, obtained from oxygen isotopes and Mg/Ca ratios in planktonic foraminifera from the western Pacific warm pool. In contrast with the eastern equatorial regions, sea surface temperatures in the western Pacific warm pool are relatively stable throughout the Pleistocene epoch, implying little long-term change in the tropical net radiation budget. Our results challenge the hypothesis of a gradual decrease in atmospheric carbon dioxide concentrations as a dominant trigger of the longer glacial cycles since 850,000 years ago. Instead, we infer that the temperature contrast across the equatorial Pacific Ocean increased, which might have had a significant influence on the mid-Pleistocene climate transition.     

Persistence of full glacial conditions in the central Pacific until 15,000 years ago

The magnitude of atmospheric cooling during the Last Glacial Maximum and the timing of the transition into the current interglacial period remain poorly constrained in tropical regions, partly because of a lack of suitable climate records. Glacial moraines provide a method of reconstructing past temperatures, but they are relatively rare in the tropics. Here we present a reconstruction of atmospheric temperatures in the central Pacific during the last deglaciation on the basis of cosmogenic 3He ages of moraines and numerical modelling of the ice cap on Mauna Kea volcano, Hawaii--the only highland in the central Pacific on which moraines that formed during the last glacial period are preserved. Our reconstruction indicates that the Last Glacial Maximum occurred between 19,000 and 16,000 years ago in this region and that temperatures at high elevations were about 7 degrees C lower than today during this interval. Glacial retreat began about 16,000 years ago, but temperatures were still about 6.5 degrees C lower than today until 15,000 years ago. When combined with estimates of sea surface temperatures in the central Pacific Ocean, our reconstruction indicates that the lapse rate during the Last Glacial Maximum was higher than at present, which is consistent with the proposal that the atmosphere was drier at that time. Furthermore, the persistence of full glacial conditions until 15,000 years ago is consistent with the relatively late and abrupt transition to warmer temperatures in Greenland, indicating that there may have been an atmospheric teleconnection between the central Pacific and North Atlantic regions during the last deglaciation.     

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.     

A carbon isotope challenge to the snowball Earth

The snowball Earth hypothesis postulates that the planet was entirely covered by ice for millions of years in the Neoproterozoic era, in a self-enhanced glaciation caused by the high albedo of the ice-covered planet. In a hard-snowball picture, the subsequent rapid unfreezing resulted from an ultra-greenhouse event attributed to the buildup of volcanic carbon dioxide (CO(2)) during glaciation. High partial pressures of atmospheric CO(2) (pCO2; from 20,000 to 90,000?p.p.m.v.) in the aftermath of the Marinoan glaciation (～635?Myr ago) have been inferred from both boron and triple oxygen isotopes. These pCO2 values are 50 to 225 times higher than present-day levels. Here, we re-evaluate these estimates using paired carbon isotopic data for carbonate layers that cap Neoproterozoic glacial deposits and are considered to record post-glacial sea level rise. The new data reported here for Brazilian cap carbonates, together with previous ones for time-equivalent units, provide estimates lower than 3,200?p.p.m.v.--and possibly as low as the current value of ～400?p.p.m.v. Our new constraint, and our re-interpretation of the boron and triple oxygen isotope data, provide a completely different picture of the late Neoproterozoic environment, with low atmospheric concentrations of carbon dioxide and oxygen that are inconsistent with a hard-snowball Earth.     

A 'snowball Earth' climate triggered by continental break-up through changes in runoff

Geological and palaeomagnetic studies indicate that ice sheets may have reached the Equator at the end of the Proterozoic eon, 800 to 550 million years ago, leading to the suggestion of a fully ice-covered 'snowball Earth'. Climate model simulations indicate that such a snowball state for the Earth depends on anomalously low atmospheric carbon dioxide concentrations, in addition to the Sun being 6 per cent fainter than it is today. However, the mechanisms producing such low carbon dioxide concentrations remain controversial. Here we assess the effect of the palaeogeographic changes preceding the Sturtian glacial period, 750 million years ago, on the long-term evolution of atmospheric carbon dioxide levels using the coupled climate-geochemical model GEOCLIM. In our simulation, the continental break-up of Rodinia leads to an increase in runoff and hence consumption of carbon dioxide through continental weathering that decreases atmospheric carbon dioxide concentrations by 1,320 p.p.m. This indicates that tectonic changes could have triggered a progressive transition from a 'greenhouse' to an 'icehouse' climate during the Neoproterozoic era. When we combine these results with the concomitant weathering effect of the voluminous basaltic traps erupted throughout the break-up of Rodinia, our simulation results in a snowball glaciation.     

Evolution of leaf-form in land plants linked to atmospheric CO2 decline in the Late Palaeozoic era

The widespread appearance of megaphyll leaves, with their branched veins and planate form, did not occur until the close of the Devonian period at about 360 Myr ago. This happened about 40 Myr after simple leafless vascular plants first colonized the land in the Late Silurian/Early Devonian, but the reason for the slow emergence of this common feature of present-day plants is presently unresolved. Here we show, in a series of quantitative analyses using fossil leaf characters and biophysical principles, that the delay was causally linked with a 90% drop in atmospheric pCO2 during the Late Palaeozoic era. In contrast to simulations for a typical Early Devonian land plant, possessing few stomata on leafless stems, those for a planate leaf with the same stomatal characteristics indicate that it would have suffered lethal overheating, because of greater interception of solar energy and low transpiration. When planate leaves first appeared in the Late Devonian and subsequently diversified in the Carboniferous period, they possessed substantially higher stomatal densities. This observation is consistent with the effects of the pCO2 on stomatal development and suggests that the evolution of planate leaves could only have occurred after an increase in stomatal density, allowing higher transpiration rates that were sufficient to maintain cool and viable leaf temperatures.     

High CO2 levels in the Proterozoic atmosphere estimated from analyses of individual microfossils

Solar luminosity on the early Earth was significantly lower than today. Therefore, solar luminosity models suggest that, in the atmosphere of the early Earth, the concentration of greenhouse gases such as carbon dioxide and methane must have been much higher. However, empirical estimates of Proterozoic levels of atmospheric carbon dioxide concentrations have not hitherto been available. Here we present ion microprobe analyses of the carbon isotopes in individual organic-walled microfossils extracted from a Proterozoic ( approximately 1.4-gigayear-old) shale in North China. Calculated magnitudes of the carbon isotope fractionation in these large, morphologically complex microfossils suggest elevated levels of carbon dioxide in the ancient atmosphere--between 10 and 200 times the present atmospheric level. Our results indicate that carbon dioxide was an important greenhouse gas during periods of lower solar luminosity, probably dominating over methane after the atmosphere and hydrosphere became pervasively oxygenated between 2 and 2.2 gigayears ago.     

Climate sensitivity constrained by CO2 concentrations over the past 420 million years

A firm understanding of the relationship between atmospheric carbon dioxide concentration and temperature is critical for interpreting past climate change and for predicting future climate change. A recent synthesis suggests that the increase in global-mean surface temperature in response to a doubling of the atmospheric carbon dioxide concentration, termed 'climate sensitivity', is between 1.5 and 6.2 degrees C (5-95 per cent likelihood range), but some evidence is inconsistent with this range. Moreover, most estimates of climate sensitivity are based on records of climate change over the past few decades to thousands of years, when carbon dioxide concentrations and global temperatures were similar to or lower than today, so such calculations tend to underestimate the magnitude of large climate-change events and may not be applicable to climate change under warmer conditions in the future. Here we estimate long-term equilibrium climate sensitivity by modelling carbon dioxide concentrations over the past 420 million years and comparing our calculations with a proxy record. Our estimates are broadly consistent with estimates based on short-term climate records, and indicate that a weak radiative forcing by carbon dioxide is highly unlikely on multi-million-year timescales. We conclude that a climate sensitivity greater than 1.5 degrees C has probably been a robust feature of the Earth's climate system over the past 420 million years, regardless of temporal scaling.     

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).     

Interdecadal variation in the extent of South Pacific tropical waters during the Younger Dryas event

During the Younger Dryas event, about 12,000 years ago, the Northern Hemisphere cooled by between 2 and 10 degrees C (refs 1, 2) whereas East Antarctica experienced warming. But the spatial signature of the event in the southern mid-latitudes and tropics is less well known, as records are sparse and inconclusive. Here we present high-resolution analyses of skeletal Sr/Ca and 18O/16O ratios for a giant fossil Diploastrea heliopora coral that was preserved in growth position on the raised reef terraces of Espiritu Santo Island, Vanuatu, in the southwestern tropical Pacific Ocean. Our data indicate that sea surface temperatures in Vanuatu were on average 4.5 +/- 1.3 degrees C cooler during the Younger Dryas event than today, with a significant interdecadal modulation. The amplified annual cycle of sea surface temperatures, relative to today, indicates that cooling was caused by the compression of tropical waters towards the Equator. The positive correlation in our record between the oxygen isotope ratios of sea water and sea surface temperatures suggests that the South Pacific convergence zone, which brings 18O-depleted precipitation to the area today, was not active during the Younger Dryas period.     

Oxidation of the Ediacaran ocean

Oxygenation of the Earth's surface is increasingly thought to have occurred in two steps. The first step, which occurred approximately 2,300 million years (Myr) ago, involved a significant increase in atmospheric oxygen concentrations and oxygenation of the surface ocean. A further increase in atmospheric oxygen appears to have taken place during the late Neoproterozoic period ( approximately 800-542 Myr ago). This increase may have stimulated the evolution of macroscopic multicellular animals and the subsequent radiation of calcified invertebrates, and may have led to oxygenation of the deep ocean. However, the nature and timing of Neoproterozoic oxidation remain uncertain. Here we present high-resolution carbon isotope and sulphur isotope records from the Huqf Supergroup, Sultanate of Oman, that cover most of the Ediacaran period (approximately 635 to approximately 548 Myr ago). These records indicate that the ocean became increasingly oxygenated after the end of the Marinoan glaciation, and they allow us to identify three distinct stages of oxidation. When considered in the context of other records from this period, our data indicate that certain groups of eukaryotic organisms appeared and diversified during the second and third stages of oxygenation. The second stage corresponds with the Shuram excursion in the carbon isotope record and seems to have involved the oxidation of a large reservoir of organic carbon suspended in the deep ocean, indicating that this event may have had a key role in the evolution of eukaryotic organisms. Our data thus provide new insights into the oxygenation of the Ediacaran ocean and the stepwise restructuring of the carbon and sulphur cycles that occurred during this significant period of Earth's history.     

Oceanic Cd/P ratio and nutrient utilization in the glacial Southern Ocean

During glacial periods, low atmospheric carbon dioxide concentration has been associated with increased oceanic carbon uptake, particularly in the southern oceans. The mechanism involved remains unclear. Because ocean productivity is strongly influenced by nutrient levels, palaeo-oceanographic proxies have been applied to investigate nutrient utilization in surface water across glacial transitions. Here we show that present-day cadmium and phosphorus concentrations in the global oceans can be explained by a chemical fractionation during particle formation, whereby uptake of cadmium occurs in preference to uptake of phosphorus. This allows the reconstruction of past surface water phosphate concentrations from the cadmium/calcium ratio of planktonic foraminifera. Results from the Last Glacial Maximum show similar phosphate utilization in the subantarctic to that of today, but much smaller utilization in the polar Southern Ocean, in a model that is consistent with the expansion of glacial sea ice and which can reconcile all proxy records of polar nutrient utilization. By restricting communication between the ocean and atmosphere, sea ice expansion also provides a mechanism for reduced CO2 release by the Southern Ocean and lower glacial atmospheric CO2.     

Permanent El Niño during the Pliocene warm period not supported by coral evidence

The El Ni?o/Southern Oscillation (ENSO) system during the Pliocene warm period (PWP; 3-5 million years ago) may have existed in a permanent El Ni?o state with a sharply reduced zonal sea surface temperature (SST) gradient in the equatorial Pacific Ocean. This suggests that during the PWP, when global mean temperatures and atmospheric carbon dioxide concentrations were similar to those projected for near-term climate change, ENSO variability--and related global climate teleconnections-could have been radically different from that today. Yet, owing to a lack of observational evidence on seasonal and interannual SST variability from crucial low-latitude sites, this fundamental climate characteristic of the PWP remains controversial. Here we show that permanent El Ni?o conditions did not exist during the PWP. Our spectral analysis of the δ(18)O SST and salinity proxy, extracted from two 35-year, monthly resolved PWP Porites corals in the Philippines, reveals variability that is similar to present ENSO variation. Although our fossil corals cannot be directly compared with modern ENSO records, two lines of evidence suggest that Philippine corals are appropriate ENSO proxies. First, δ(18)O anomalies from a nearby live Porites coral are correlated with modern records of ENSO variability. Second, negative-δ(18)O events in the fossil corals closely resemble the decreases in δ(18)O seen in the live coral during El Ni?o events. Prior research advocating a permanent El Ni?o state may have been limited by the coarse resolution of many SST proxies, whereas our coral-based analysis identifies climate variability at the temporal scale required to resolve ENSO structure firmly.     

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.