Eocene bipolar glaciation associated with global carbon cycle changes

The transition from the extreme global warmth of the early Eocene 'greenhouse' climate approximately 55 million years ago to the present glaciated state is one of the most prominent changes in Earth's climatic evolution. It is widely accepted that large ice sheets first appeared on Antarctica approximately 34 million years ago, coincident with decreasing atmospheric carbon dioxide concentrations and a deepening of the calcite compensation depth in the world's oceans, and that glaciation in the Northern Hemisphere began much later, between 10 and 6 million years ago. Here we present records of sediment and foraminiferal geochemistry covering the greenhouse-icehouse climate transition. We report evidence for synchronous deepening and subsequent oscillations in the calcite compensation depth in the tropical Pacific and South Atlantic oceans from approximately 42 million years ago, with a permanent deepening 34 million years ago. The most prominent variations in the calcite compensation depth coincide with changes in seawater oxygen isotope ratios of up to 1.5 per mil, suggesting a lowering of global sea level through significant storage of ice in both hemispheres by at least 100 to 125 metres. Variations in benthic carbon isotope ratios of up to approximately 1.4 per mil occurred at the same time, indicating large changes in carbon cycling. We suggest that the greenhouse-icehouse transition was closely coupled to the evolution of atmospheric carbon dioxide, and that negative carbon cycle feedbacks may have prevented the permanent establishment of large ice sheets earlier than 34 million years ago.     

No extreme bipolar glaciation during the main Eocene calcite compensation shift

Major ice sheets were permanently established on Antarctica approximately 34 million years ago, close to the Eocene/Oligocene boundary, at the same time as a permanent deepening of the calcite compensation depth in the world's oceans. Until recently, it was thought that Northern Hemisphere glaciation began much later, between 11 and 5 million years ago. This view has been challenged, however, by records of ice rafting at high northern latitudes during the Eocene epoch and by estimates of global ice volume that exceed the storage capacity of Antarctica at the same time as a temporary deepening of the calcite compensation depth approximately 41.6 million years ago. Here we test the hypothesis that large ice sheets were present in both hemispheres approximately 41.6 million years ago using marine sediment records of oxygen and carbon isotope values and of calcium carbonate content from the equatorial Atlantic Ocean. These records allow, at most, an ice budget that can easily be accommodated on Antarctica, indicating that large ice sheets were not present in the Northern Hemisphere. The records also reveal a brief interval shortly before the temporary deepening of the calcite compensation depth during which the calcite compensation depth shoaled, ocean temperatures increased and carbon isotope values decreased in the equatorial Atlantic. The nature of these changes around 41.6 million years ago implies common links, in terms of carbon cycling, with events at the Eocene/Oligocene boundary and with the 'hyperthermals' of the Early Eocene climate optimum. Our findings help to resolve the apparent discrepancy between the geological records of Northern Hemisphere glaciation and model results that indicate that the threshold for continental glaciation was crossed earlier in the Southern Hemisphere than in the Northern Hemisphere.     

Thinning of the ice sheet in northwest Greenland over the past forty years.

Thermal expansion of the oceans, as well as melting of glaciers, ice sheets and ice caps have been the main contributors to global sea level rise over the past century. The greatest uncertainty in predicting future sea level changes lies with our estimates of the mass balance of the ice sheets in Greenland and Antarctica. Satellite measurements have been used to determine changes in these ice sheets on short timescales, demonstrating that surface-elevation changes on timescales of decades or less result mainly from variations in snow accumulation. Here we present direct measurements of the changes in surface elevation between 1954 and 1995 on a traverse across the north Greenland ice sheet. Measurements over a time interval of this length should reflect changes in ice flow-the important quantity for predicting changes in sea level-relatively unperturbed by short-term fluctuations in snow accumulation. We find only small changes in the eastern part of the transect, except for some thickening of the north ice stream. On the west side, however, the thinning rates of the ice sheet are significantly higher and thinning extends to higher elevations than had been anticipated from previous studies.     

Similar meltwater contributions to glacial sea level changes from Antarctic and northern ice sheets

The period between 75,000 and 20,000 years ago was characterized by high variability in climate and sea level. Southern Ocean records of ice-rafted debris suggest a significant contribution to the sea level changes from melt water of Antarctic origin, in addition to likely contributions from northern ice sheets, but the relative volumes of melt water from northern and southern sources have yet to be established. Here we simulate the first-order impact of a range of relative meltwater releases from the two polar regions on the distribution of marine oxygen isotopes, using an intermediate complexity model. By comparing our simulations with oxygen isotope data from sediment cores, we infer that the contributions from Antarctica and the northern ice sheets to the documented sea level rises between 65,000 and 35,000 years ago were approximately equal, each accounting for a rise of about 15 m. The reductions in Antarctic ice volume implied by our analysis are comparable to that inferred previously for the Antarctic contribution to meltwater pulse 1A (refs 16, 17), which occurred about 14,200 years ago, during the last deglaciation.     

Extreme winds and waves in the aftermath of a Neoproterozoic glaciation

The most severe excursions in the Earth's climatic history are thought to be associated with Proterozoic glaciations. According to the 'Snowball Earth' hypothesis, the Marinoan glaciation, which ended about 635 million years ago, involved global or nearly global ice cover. At the termination of this glacial period, rapid melting of continental ice sheets must have caused a large rise in sea level. Here we show that sediments deposited during this sea level rise contain remarkable structures that we interpret as giant wave ripples. These structures occur at homologous stratigraphic levels in Australia, Brazil, Canada, Namibia and Svalbard. Our hydrodynamic analysis of these structures suggests maximum wave periods of 21 to 30 seconds, significantly longer than those typical for today's oceans. The reconstructed wave conditions could only have been generated under sustained high wind velocities exceeding 20 metres per second in fetch-unlimited ocean basins. We propose that these extraordinary wind and wave conditions were characteristic of the climatic transit, and provide observational targets for atmospheric circulation models.     

Modelled atmospheric temperatures and global sea levels over the past million years

Marine records of sediment oxygen isotope compositions show that the Earth's climate has gone through a succession of glacial and interglacial periods during the past million years. But the interpretation of the oxygen isotope records is complicated because both isotope storage in ice sheets and deep-water temperature affect the recorded isotopic composition. Separating these two effects would require long records of either sea level or deep-ocean temperature, which are currently not available. Here we use a coupled model of the Northern Hemisphere ice sheets and ocean temperatures, forced to match an oxygen isotope record for the past million years compiled from 57 globally distributed sediment cores, to quantify both contributions simultaneously. We find that the ice-sheet contribution to the variability in oxygen isotope composition varied from ten per cent in the beginning of glacial periods to sixty per cent at glacial maxima, suggesting that strong ocean cooling preceded slow ice-sheet build-up. The model yields mutually consistent time series of continental mean surface temperatures between 40 and 80 degrees N, ice volume and global sea level. We find that during extreme glacial stages, air temperatures were 17 +/- 1.8 degrees C lower than present, with a 120 +/- 10 m sea level equivalent of continental ice present.     

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.     

Continental ice in Greenland during the Eocene and Oligocene

The Eocene and Oligocene epochs (approximately 55 to 23 million years ago) comprise a critical phase in Earth history. An array of geological records supported by climate modelling indicates a profound shift in global climate during this interval, from a state that was largely free of polar ice caps to one in which ice sheets on Antarctica approached their modern size. However, the early glaciation history of the Northern Hemisphere is a subject of controversy. Here we report stratigraphically extensive ice-rafted debris, including macroscopic dropstones, in late Eocene to early Oligocene sediments from the Norwegian-Greenland Sea that were deposited between about 38 and 30 million years ago. Our data indicate sediment rafting by glacial ice, rather than sea ice, and point to East Greenland as the likely source. Records of this type from one site alone cannot be used to determine the extent of ice involved. However, our data suggest the existence of (at least) isolated glaciers on Greenland about 20 million years earlier than previously documented, at a time when temperatures and atmospheric carbon dioxide concentrations were substantially higher.     

A dynamic early East Antarctic Ice Sheet suggested by ice-covered fjord landscapes

The first Cenozoic ice sheets initiated in Antarctica from the Gamburtsev Subglacial Mountains and other highlands as a result of rapid global cooling ～34 million years ago. In the subsequent 20 million years, at a time of declining atmospheric carbon dioxide concentrations and an evolving Antarctic circumpolar current, sedimentary sequence interpretation and numerical modelling suggest that cyclical periods of ice-sheet expansion to the continental margin, followed by retreat to the subglacial highlands, occurred up to thirty times. These fluctuations were paced by orbital changes and were a major influence on global sea levels. Ice-sheet models show that the nature of such oscillations is critically dependent on the pattern and extent of Antarctic topographic lowlands. Here we show that the basal topography of the Aurora Subglacial Basin of East Antarctica, at present overlain by 2-4.5?km of ice, is characterized by a series of well-defined topographic channels within a mountain block landscape. The identification of this fjord landscape, based on new data from ice-penetrating radar, provides an improved understanding of the topography of the Aurora Subglacial Basin and its surroundings, and reveals a complex surface sculpted by a succession of ice-sheet configurations substantially different from today's. At different stages during its fluctuations, the edge of the East Antarctic Ice Sheet lay pinned along the margins of the Aurora Subglacial Basin, the upland boundaries of which are currently above sea level and the deepest parts of which are more than 1?km below sea level. Although the timing of the channel incision remains uncertain, our results suggest that the fjord landscape was carved by at least two iceflow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.     

Regional climate shifts caused by gradual global cooling in the Pliocene epoch

The Earth's climate has undergone a global transition over the past four million years, from warm conditions with global surface temperatures about 3 degrees C warmer than today, smaller ice sheets and higher sea levels to the current cooler conditions. Tectonic changes and their influence on ocean heat transport have been suggested as forcing factors for that transition, including the onset of significant Northern Hemisphere glaciation approximately 2.75 million years ago, but the ultimate causes for the climatic changes are still under debate. Here we compare climate records from high latitudes, subtropical regions and the tropics, indicating that the onset of large glacial/interglacial cycles did not coincide with a specific climate reorganization event at lower latitudes. The regional differences in the timing of cooling imply that global cooling was a gradual process, rather than the response to a single threshold or episodic event as previously suggested. We also find that high-latitude climate sensitivity to variations in solar heating increased gradually, culminating after cool tropical and subtropical upwelling conditions were established two million years ago. Our results suggest that mean low-latitude climate conditions can significantly influence global climate feedbacks.     

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

Recent mass balance of polar ice sheets inferred from patterns of global sea-level change

Global sea level is an indicator of climate change, as it is sensitive to both thermal expansion of the oceans and a reduction of land-based glaciers. Global sea-level rise has been estimated by correcting observations from tide gauges for glacial isostatic adjustment--the continuing sea-level response due to melting of Late Pleistocene ice--and by computing the global mean of these residual trends. In such analyses, spatial patterns of sea-level rise are assumed to be signals that will average out over geographically distributed tide-gauge data. But a long history of modelling studies has demonstrated that non-uniform--that is, non-eustatic--sea-level redistributions can be produced by variations in the volume of the polar ice sheets. Here we present numerical predictions of gravitationally consistent patterns of sea-level change following variations in either the Antarctic or Greenland ice sheets or the melting of a suite of small mountain glaciers. These predictions are characterized by geometrically distinct patterns that reconcile spatial variations in previously published sea-level records. Under the--albeit coarse--assumption of a globally uniform thermal expansion of the oceans, our approach suggests melting of the Greenland ice complex over the last century equivalent to -0.6 mm yr(-1) of sea-level rise.     

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.     

Dynamics of ice ages on Mars

Unlike Earth, where astronomical climate forcing is comparatively small, Mars experiences dramatic changes in incident sunlight that are capable of redistributing ice on a global scale. The geographic extent of the subsurface ice found poleward of approximately +/-60 degrees latitude on both hemispheres of Mars coincides with the areas where ice is stable. However, the tilt of Mars' rotation axis (obliquity) changed considerably in the past several million years. Earlier work has shown that regions of ice stability, which are defined by temperature and atmospheric humidity, differed in the recent past from today's, and subsurface ice is expected to retreat quickly when unstable. Here I explain how the subsurface ice sheets could have evolved to the state in which we see them today. Simulations of the retreat and growth of ground ice as a result of sublimation loss and recharge reveal forty major ice ages over the past five million years. Today, this gives rise to pore ice at mid-latitudes and a three-layered depth distribution in the high latitudes of, from top to bottom, a dry layer, pore ice, and a massive ice sheet. Combined, these layers provide enough ice to be compatible with existing neutron and gamma-ray measurements.     

Deglacial rapid sea level rises caused by ice-sheet saddle collapses

The last deglaciation (21 to 7 thousand years ago) was punctuated by several abrupt meltwater pulses, which sometimes caused noticeable climate change. Around 14 thousand years ago, meltwater pulse 1A (MWP-1A), the largest of these events, produced a sea level rise of 14-18?metres over 350?years. Although this enormous surge of water certainly originated from retreating ice sheets, there is no consensus on the geographical source or underlying physical mechanisms governing the rapid sea level rise. Here we present an ice-sheet modelling simulation in which the separation of the Laurentide and Cordilleran ice sheets in North America produces a meltwater pulse corresponding to MWP-1A. Another meltwater pulse is produced when the Labrador and Baffin ice domes around Hudson Bay separate, which could be associated with the '8,200-year' event, the most pronounced abrupt climate event of the past nine thousand years. For both modelled pulses, the saddle between the two ice domes becomes subject to surface melting because of a general surface lowering caused by climate warming. The melting then rapidly accelerates as the saddle between the two domes gets lower, producing nine metres of sea level rise over 500 years. This mechanism of an ice 'saddle collapse' probably explains MWP-1A and the 8,200-year event and sheds light on the consequences of these events on climate.     

North Pacific seasonality and the glaciation of North America 2.7 million years ago

In the context of gradual Cenozoic cooling, the timing of the onset of significant Northern Hemisphere glaciation 2.7 million years ago is consistent with Milankovitch's orbital theory, which posited that ice sheets grow when polar summertime insolation and temperature are low. However, the role of moisture supply in the initiation of large Northern Hemisphere ice sheets has remained unclear. The subarctic Pacific Ocean represents a significant source of water vapour to boreal North America, but it has been largely overlooked in efforts to explain Northern Hemisphere glaciation. Here we present alkenone unsaturation ratios and diatom oxygen isotope ratios from a sediment core in the western subarctic Pacific Ocean, indicating that 2.7 million years ago late-summer sea surface temperatures in this ocean region rose in response to an increase in stratification. At the same time, winter sea surface temperatures cooled, winter floating ice became more abundant and global climate descended into glacial conditions. We suggest that the observed summer warming extended into the autumn, providing water vapour to northern North America, where it precipitated and accumulated as snow, and thus allowed the initiation of Northern Hemisphere glaciation.     

Eocene/Oligocene ocean de-acidification linked to Antarctic glaciation by sea-level fall

One of the most dramatic perturbations to the Earth system during the past 100 million years was the rapid onset of Antarctic glaciation near the Eocene/Oligocene epoch boundary (approximately 34 million years ago). This climate transition was accompanied by a deepening of the calcite compensation depth--the ocean depth at which the rate of calcium carbonate input from surface waters equals the rate of dissolution. Changes in the global carbon cycle, rather than changes in continental configuration, have recently been proposed as the most likely root cause of Antarctic glaciation, but the mechanism linking glaciation to the deepening of calcite compensation depth remains unclear. Here we use a global biogeochemical box model to test competing hypotheses put forward to explain the Eocene/Oligocene transition. We find that, of the candidate hypotheses, only shelf to deep sea carbonate partitioning is capable of explaining the observed changes in both carbon isotope composition and calcium carbonate accumulation at the sea floor. In our simulations, glacioeustatic sea-level fall associated with the growth of Antarctic ice sheets permanently reduces global calcium carbonate accumulation on the continental shelves, leading to an increase in pelagic burial via permanent deepening of the calcite compensation depth. At the same time, fresh limestones are exposed to erosion, thus temporarily increasing global river inputs of dissolved carbonate and increasing seawater delta13C. Our work sheds new light on the mechanisms linking glaciation and ocean acidity change across arguably the most important climate transition of the Cenozoic era.     

Dual modes of the carbon cycle since the Last Glacial Maximum

The most conspicuous feature of the record of past climate contained in polar ice is the rapid warming which occurs after long intervals of gradual cooling. During the last four transitions from glacial to interglacial conditions, over which such abrupt warmings occur, ice records indicate that the CO2 concentration of the atmosphere increased by roughly 80 to 100 parts per million by volume. But the causes of the atmospheric CO2 concentration increases are unclear. Here we present the stable-carbon-isotope composition (delta 13 CO2) of CO2 extracted from air trapped in ice at Taylor Dome, Antarctica, from the Last Glacial Maximum to the onset of Holocene times. The global carbon cycle is shown to have operated in two distinct primary modes on the timescale of thousands of years, one when climate was changing relatively slowly and another when warming was rapid, each with a characteristic average stable-carbon-isotope composition of the net CO2 exchanged by the atmosphere with the land and oceans. delta 13 CO2 increased between 16.5 and 9 thousand years ago by slightly more than would be estimated to be caused by the physical effects of a 5 degrees C rise in global average sea surface temperature driving a CO2 efflux from the ocean, but our data do not allow specific causes to be constrained.     

Subtropical Arctic Ocean temperatures during the Palaeocene/Eocene thermal maximum

The Palaeocene/Eocene thermal maximum, approximately 55 million years ago, was a brief period of widespread, extreme climatic warming, that was associated with massive atmospheric greenhouse gas input. Although aspects of the resulting environmental changes are well documented at low latitudes, no data were available to quantify simultaneous changes in the Arctic region. Here we identify the Palaeocene/Eocene thermal maximum in a marine sedimentary sequence obtained during the Arctic Coring Expedition. We show that sea surface temperatures near the North Pole increased from 18 degrees C to over 23 degrees C during this event. Such warm values imply the absence of ice and thus exclude the influence of ice-albedo feedbacks on this Arctic warming. At the same time, sea level rose while anoxic and euxinic conditions developed in the ocean's bottom waters and photic zone, respectively. Increasing temperature and sea level match expectations based on palaeoclimate model simulations, but the absolute polar temperatures that we derive before, during and after the event are more than 10 degrees C warmer than those model-predicted. This suggests that higher-than-modern greenhouse gas concentrations must have operated in conjunction with other feedback mechanisms--perhaps polar stratospheric clouds or hurricane-induced ocean mixing--to amplify early Palaeogene polar temperatures.     

High temperatures in the Late Cretaceous Arctic Ocean

To understand the climate dynamics of the warm, equable greenhouse world of the Late Cretaceous period, it is important to determine polar palaeotemperatures. The early palaeoceanographic history of the Arctic Ocean has, however, remained largely unknown, because the sea floor and underlying deposits are usually inaccessible beneath a cover of floating ice. A shallow piston core taken from a drifting ice island in 1970 fortuitously retrieved unconsolidated Upper Cretaceous organic-rich sediment from Alpha ridge, a submarine elevated feature of probable oceanic origin. A lack of carbonate in the sediments from this core has prevented the use of traditional oxygen-isotope palaeothermometry. Here we determine Arctic palaeotemperatures from these Upper Cretaceous deposits using TEX86, a new palaeothermometer that is based on the composition of membrane lipids derived from a ubiquitous component of marine plankton, Crenarchaeota. From these analyses we infer an average sea surface temperature of approximately 15 degrees C for the Arctic Ocean about 70 million years ago. This calibration point implies an Equator-to-pole gradient in sea surface temperatures of approximately 15 degrees C during this interval and, by extrapolation, we suggest that polar waters were generally warmer than 20 degrees C during the middle Cretaceous (approximately 90 million years ago).