Recent ice ages on Mars

A key pacemaker of ice ages on the Earth is climatic forcing due to variations in planetary orbital parameters. Recent Mars exploration has revealed dusty, water-ice-rich mantling deposits that are layered, metres thick and latitude dependent, occurring in both hemispheres from mid-latitudes to the poles. Here we show evidence that these deposits formed during a geologically recent ice age that occurred from about 2.1 to 0.4 Myr ago. The deposits were emplaced symmetrically down to latitudes of approximately 30 degrees--equivalent to Saudi Arabia and the southern United States on the Earth--in response to the changing stability of water ice and dust during variations in obliquity (the angle between Mars' pole of rotation and the ecliptic plane) reaching 30-35 degrees. Mars is at present in an 'interglacial' period, and the ice-rich deposits are undergoing reworking, degradation and retreat in response to the current instability of near-surface ice. Unlike the Earth, martian ice ages are characterized by warmer polar climates and enhanced equatorward transport of atmospheric water and dust to produce widespread smooth deposits down to mid-latitudes.     

Evidence for recent climate change on Mars from the identification of youthful near-surface ground ice

Ground ice in the crust and soil may be one of the largest reservoirs of water on Mars. Near-surface ground ice is predicted to be stable at latitudes higher than 40 degrees (ref. 4), where a number of geomorphologic features indicative of viscous creep and hence ground ice have been observed. Mid-latitude soils have also been implicated as a water-ice reservoir, the capacity of which is predicted to vary on a 100,000-year timescale owing to orbitally driven variations in climate. It is uncertain, however, whether near-surface ground ice currently exists at these latitudes, and how it is changing with time. Here we report observational evidence for a mid-latitude reservoir of near-surface water ice occupying the pore space of soils. The thickness of the ice-occupied soil reservoir (1-10 m) and its distribution in the 30 degrees to 60 degrees latitude bands indicate a reservoir of (1.5-6.0) x 104 km3, equivalent to a global layer of water 10-40 cm thick. We infer that the reservoir was created during the last phase of high orbital obliquity less than 100,000 years ago, and is now being diminished.     

Tropical to mid-latitude snow and ice accumulation, flow and glaciation on Mars

Images from the Mars Express HRSC (High-Resolution Stereo Camera) of debris aprons at the base of massifs in eastern Hellas reveal numerous concentrically ridged lobate and pitted features and related evidence of extremely ice-rich glacier-like viscous flow and sublimation. Together with new evidence for recent ice-rich rock glaciers at the base of the Olympus Mons scarp superposed on larger Late Amazonian debris-covered piedmont glaciers, we interpret these deposits as evidence for geologically recent and recurring glacial activity in tropical and mid-latitude regions of Mars during periods of increased spin-axis obliquity when polar ice was mobilized and redeposited in microenvironments at lower latitudes. The data indicate that abundant residual ice probably remains in these deposits and that these records of geologically recent climate changes are accessible to future automated and human surface exploration.     

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.     

North-south geological differences between the residual polar caps on Mars

Polar processes can be sensitive indicators of global climate, and the geological features associated with polar ice caps can therefore indicate evolution of climate with time. The polar regions on Mars have distinctive morphologic and climatologic features: thick layered deposits, seasonal CO2 frost caps extending to mid latitudes, and near-polar residual frost deposits that survive the summer. The relationship of the seasonal and residual frost caps to the layered deposits has been poorly constrained, mainly by the limited spatial resolution of the available data. In particular, it has not been known if the residual caps represent simple thin frost cover or substantial geologic features. Here we show that the residual cap on the south pole is a distinct geologic unit with striking collapse and erosional topography; this is very different from the residual cap on the north pole, which grades into the underlying layered materials. These findings indicate that the differences between the caps are substantial (rather than reflecting short-lived differences in frost cover), and so support the idea of long-term asymmetry in the polar climates of Mars.     

Geology and insolation-driven climatic history of Amazonian north polar materials on Mars

Mariner 9 and Viking spacecraft images revealed that the polar regions of Mars, like those of Earth, record the planet's climate history. However, fundamental uncertainties regarding the materials, features, ages and processes constituting the geologic record remained. Recently acquired Mars Orbiter Laser Altimeter data and Mars Orbiter Camera high-resolution images from the Mars Global Surveyor spacecraft and moderately high-resolution Thermal Emission Imaging System visible images from the Mars Odyssey spacecraft permit more comprehensive geologic and climatic analyses. Here I map and show the history of geologic materials and features in the north polar region that span the Amazonian period (approximately 3.0 Gyr ago to present). Erosion and redeposition of putative circumpolar mud volcano deposits (formed by eruption of liquefied, fine-grained material) led to the formation of an Early Amazonian polar plateau consisting of dark layered materials. Crater ejecta superposed on pedestals indicate that a thin mantle was present during most of the Amazonian, suggesting generally higher obliquity and insolation conditions at the poles than at present. Brighter polar layered deposits rest unconformably on the dark layers and formed mainly during lower obliquity over the past 4-5 Myr (ref. 20). Finally, the uppermost layers post-date the latest downtrend in obliquity <20,000 years ago.     

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.     

Perennial water ice identified in the south polar cap of Mars

The inventory of water and carbon dioxide reservoirs on Mars are important clues for understanding the geological, climatic and potentially exobiological evolution of the planet. From the early mapping observation of the permanent ice caps on the martian poles, the northern cap was believed to be mainly composed of water ice, whereas the southern cap was thought to be constituted of carbon dioxide ice. However, recent missions (NASA missions Mars Global Surveyor and Odyssey) have revealed surface structures, altimetry profiles, underlying buried hydrogen, and temperatures of the south polar regions that are thermodynamically consistent with a mixture of surface water ice and carbon dioxide. Here we present the first direct identification and mapping of both carbon dioxide and water ice in the martian high southern latitudes, at a resolution of 2 km, during the local summer, when the extent of the polar ice is at its minimum. We observe that this south polar cap contains perennial water ice in extended areas: as a small admixture to carbon dioxide in the bright regions; associated with dust, without carbon dioxide, at the edges of this bright cap; and, unexpectedly, in large areas tens of kilometres away from the bright cap.     

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

A topographically forced asymmetry in the martian circulation and climate

Large seasonal and hemispheric asymmetries in the martian climate system are generally ascribed to variations in solar heating associated with orbital eccentricity. As the orbital elements slowly change (over a period of >104 years), characteristics of the climate such as dustiness and the vigour of atmospheric circulation are thought to vary, as should asymmetries in the climate (for example, the deposition of water ice at the northern versus the southern pole). Such orbitally driven climate change might be responsible for the observed layering in Mars' polar deposits by modulating deposition of dust and water ice. Most current theories assume that climate asymmetries completely reverse as the angular distance between equinox and perihelion changes by 180 degrees. Here we describe a major climate mechanism that will not precess in this way. We show that Mars' global north-south elevation difference forces a dominant southern summer Hadley circulation that is independent of perihelion timing. The Hadley circulation, a tropical overturning cell responsible for trade winds, largely controls interhemispheric transport of water and the bulk dustiness of the atmosphere. The topography therefore imprints a strong handedness on climate, with water ice and the active formation of polar layered deposits more likely in the north.     

Nonlinear dynamics of ice-wedge networks and resulting sensitivity to severe cooling events

Patterns of subsurface wedges of ice that form along cooling-induced tension fractures, expressed at the ground surface by ridges or troughs spaced 10 30 m apart, are ubiquitous in polar lowlands. Fossilized ice wedges, which are widespread at lower latitudes, have been used to infer the duration and mean temperature of cold periods within Proterozoic and Quaternary climates, and recent climate trends have been inferred from fracture frequency in active ice wedges. Here we present simulations from a numerical model for the evolution of ice-wedge networks over a range of climate scenarios, based on the interactions between thermal tensile stress, fracture and ice wedges. We find that short-lived periods of severe cooling permanently alter the spacing between ice wedges as well as their fracture frequency. This affects the rate at which the widths of ice wedges increase as well as the network's response to subsequent climate change. We conclude that wedge spacing and width in ice-wedge networks mainly reflect infrequent episodes of rapidly falling ground temperatures rather than mean conditions.     

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

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.     

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

No evidence for thick deposits of ice at the lunar south pole

Shackleton crater at the Moon's south pole has been suggested as a possible site of concentrated deposits of water ice, on the basis of modelling of bi-static radar polarization properties and interpretations of earlier Earth-based radar images. This suggestion, and parallel assumptions about other topographic cold traps, is a significant element in planning for future lunar landings. Hydrogen enhancements have been identified in the polar regions, but these data do not identify the host species or its local distribution. The earlier Earth-based radar data lack the resolution and coverage for detailed studies of the relationship between radar scattering properties, cold traps in permanently shadowed areas, and local terrain features such as the walls and ejecta of small craters. Here we present new 20-m resolution, 13-cm-wavelength radar images that show no evidence for concentrated deposits of water ice in Shackleton crater or elsewhere at the south pole. The polarization properties normally associated with reflections from icy surfaces in the Solar System were found at all the observed latitudes and are strongly correlated with the rock-strewn walls and ejecta of young craters, including the inner wall of Shackleton. There is no correlation between the polarization properties and the degree of solar illumination. If the hydrogen enhancement observed by the Lunar Prospector orbiter indicates the presence of water ice, then our data are consistent with the ice being present only as disseminated grains in the lunar regolith.     

Possible methane-induced polar warming in the early Eocene

Reconstructions of early Eocene climate depict a world in which the polar environments support mammals and reptiles, deciduous forests, warm oceans and rare frost conditions. At the same time, tropical sea surface temperatures are interpreted to have been the same as or slightly cooler than present values. The question of how to warm polar regions of Earth without noticeably warming the tropics remains unresolved; increased amounts of greenhouse gases would be expected to warm all latitudes equally. Oceanic heat transport has been postulated as a mechanism for heating high latitudes, but it is difficult to explain the dynamics that would achieve this. Here we consider estimates of Eocene wetland areas and suggest that the flux of methane, an important greenhouse gas, may have been substantially greater during the Eocene than at present. Elevated methane concentrations would have enhanced early Eocene global warming, and also might specifically have prevented severe winter cooling of polar regions because of the potential of atmospheric methane to promote the formation of optically thick, polar stratospheric ice clouds.     

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.     

Dust-climate couplings over the past 800,000 years from the EPICA Dome C ice core

Dust can affect the radiative balance of the atmosphere by absorbing or reflecting incoming solar radiation; it can also be a source of micronutrients, such as iron, to the ocean. It has been suggested that production, transport and deposition of dust is influenced by climatic changes on glacial-interglacial timescales. Here we present a high-resolution record of aeolian dust from the EPICA Dome C ice core in East Antarctica, which provides an undisturbed climate sequence over the past eight climatic cycles. We find that there is a significant correlation between dust flux and temperature records during glacial periods that is absent during interglacial periods. Our data suggest that dust flux is increasingly correlated with Antarctic temperature as the climate becomes colder. We interpret this as progressive coupling of the climates of Antarctic and lower latitudes. Limited changes in glacial-interglacial atmospheric transport time suggest that the sources and lifetime of dust are the main factors controlling the high glacial dust input. We propose that the observed approximately 25-fold increase in glacial dust flux over all eight glacial periods can be attributed to a strengthening of South American dust sources, together with a longer lifetime for atmospheric dust particles in the upper troposphere resulting from a reduced hydrological cycle during the ice ages.     

Orbital forcing of the martian polar layered deposits

Since the first images of polar regions on Mars revealed alternating bright and dark layers, there has been speculation that their formation might be tied to the planet's orbital climate forcing. But uncertainties in the deposition timescale exceed two orders of magnitude: estimates based on assumptions of dust deposition, ice formation and sublimation, and their variations with orbital forcing suggest a deposition rate of 10(-3) to 10(-2) cm yr(-1) (refs 5, 6), whereas estimates based on cratering rate result in values as high as 0.1 to 0.2 cm yr(-1) (ref. 7). Here we use a combination of high-resolution images of the polar layered terrains, high-resolution topography and revised calculations of the orbital and rotational parameters of Mars to show that a correlation exists between ice-layer radiance as a function of depth (obtained from photometric data of the images of the layered terrains) and the insolation variations in summer at the martian north pole, similar to what has been shown for palaeoclimate studies of the Earth. For the best fit between the radiance profile and the simulated insolation parameters, we obtain an average deposition rate of 0.05 cm yr(-1) for the top 250 m of deposits on the ice cap of the north pole of Mars.     

A 780-year record of explosive volcanism from DT263 ice core in east Antarctica

Ice cores recovered from polar ice sheet received and preserved sulfuric acid fallout from explosive volcanic eruptions. DT263 ice core was retrieved from an east Antarctic location. The ice core is dated using a combination of annual layer counting and volcanic time stratigraphic horizon as 780 years (1215-1996 A.D.). The ice core record demonstrates that during the period of approximately 1460-1800 A.D., the accumulation is sharply lower than the levels prior to and after this period. This period coincides with the most recent neoglacial climatic episode, the "Little Ice Age (LIA)", that has been found in numerous Northern Hemisphere proxy and historic records. The non-sea-salt SO42- concentrations indicate seventeen volcanic events in DT263 ice core. Compared with those from previous Antarctic ice cores, significant discrepancies are found between these records in relative volcanic flux of several well-known events. The discrepancies among these records may be explained by the differences in surface topography, accumulation rate, snow drift and distribution which highlight the potential impact of local glaci-ology on ice core volcanic records, analytical techniques used for sulfate measurement, etc. Volcanic eruptions in middle and high southern latitudes affect volcanic records in Antarctic snow more intensively than those in the low latitudes.