Clay mineralogy of ODP Leg 178 sites

The clay mineral compositions of upper Miocene to Quaternary sediments recovered at Ocean Drilling Program (ODP) Leg 178, Sites 1095 and 1096, from the continental rise west of the Antarctic Peninsula were analyzed in order to reconstruct the Neogene and Quaternary Antarctic paleoclimate and ice dynamics. The clay mineral assemblages are dominated by smectite, illite, and chlorite. Kaolinite occurs only in trace amounts. Analysis of a surface-sample data set facilitates the assignment of these clay minerals to particular source areas on the Antarctic Peninsula and, thus, the reconstruction of transport pathways. In the ODP cores, clay mineral composition cyclically alternates between two end-member assemblages. One assemblage is characterized by <20% smectite and >40% chlorite. The other assemblage has >20% smectite and <40% chlorite. Illite fluctuates between 30% and 50% without a significant affinity to one end-member assemblage. By comparison with a Quaternary sediment sequence from gravity core PS1565, the clay mineral fluctuations can be ascribed to glacial and interglacial periods, respectively. The cyclic changes in the clay mineral composition suggest that glacial-interglacial cycles, repeated ice advances and retreats, and changes in the Antarctic ice volume were already a main control of the global climate in late Miocene time. Throughout the late Neogene and Quaternary, the clay mineral records in the drift sediments exhibit only slight long-term changes predominantly attributed to local changes in glacial erosion and supply of source rocks. The absence of clear long-term trends associated with major climatic or glaciological changes points to an onset of vast glaciation in the Antarctic Peninsula region before ~9 Ma and to relative stability of the Antarctic ice sheet since then.

Nd-Sr isotopes, diatom and biogenic opal content of Pliocene sediments from IODP Site 318-U1361

Warm intervals within the Pliocene epoch (5.33-2.58 million years ago) were characterized by global temperatures comparable to those predicted for the end of this century (Haywood and Valdes, doi:10.1016/S0012-821X(03)00685-X) and atmospheric CO2 concentrations similar to today (Seki et al., 2010, doi:10.1016/j.epsl.2010.01.037; Bartoli et al., 2011, doi:10.1029/2010PA002055; Pagani et al., 2010, doi:10.1038/ngeo724). Estimates for global sea level highstands during these times (Miller et al., 2012, doi:10.1130/G32869.1) imply possible retreat of the East Antarctic ice sheet, but ice-proximal evidence from the Antarctic margin is scarce. Here we present new data from Pliocene marine sediments recovered offshore of Adélie Land, East Antarctica, that reveal dynamic behaviour of the East Antarctic ice sheet in the vicinity of the low-lying Wilkes Subglacial Basin during times of past climatic warmth. Sedimentary sequences deposited between 5.3 and 3.3 million years ago indicate increases in Southern Ocean surface water productivity, associated with elevated circum-Antarctic temperatures. The geochemical provenance of detrital material deposited during these warm intervals suggests active erosion of continental bedrock from within the Wilkes Subglacial Basin, an area today buried beneath the East Antarctic ice sheet. We interpret this erosion to be associated with retreat of the ice sheet margin several hundreds of kilometres inland and conclude that the East Antarctic ice sheet was sensitive to climatic warmth during the Pliocene.

(Table 1) Accumulation and ice discharge of West Antarctic glaciers draining into the Amundsen Sea between 1974 and 2007

The Advanced Land Observation System (ALOS) Phased-Array Synthetic-Aperture Radar (PALSAR) is an L-band frequency (1.27 GHz) radar capable of continental-scale interferometric observations of ice sheet motion. Here, we show that PALSAR data yield excellent measurements of ice motion compared to C-band (5.6 GHz) radar data because of greater temporal coherence over snow and firn. We compare PALSAR velocities from year 2006 in Pine Island Bay, West Antarctica with those spanning years 1974 to 2007. Between 1996 and 2007, Pine Island Glacier sped up 42% and ungrounded over most of its ice plain. Smith Glacier accelerated 83% and ungrounded as well. Their largest speed up are recorded in 2007. Thwaites Glacier is not accelerating but widening with time and its eastern ice shelf doubled its speed. Total ice discharge from these glaciers increased 30% in 12 yr and the net mass loss increased 170% from 39 ± 15 Gt/yr to 105 ± 27 Gt/yr. Longer-term velocity changes suggest only a moderate loss in the 1970s. As the glaciers unground into the deeper, smoother beds inland, the mass loss from this region will grow considerably larger in years to come.

(Data S1) Palynological analysis of IODP and ODP sediment cores around the Antarctic margin

The circum-Antarctic Southern Ocean is an important region for global marine food webs and carbon cycling because of sea-ice formation and its unique plankton ecosystem. However, the mechanisms underlying the installation of this distinct ecosystem and the geological timing of its development remain unknown. Here, we show, on the basis of fossil marine dinoflagellate cyst records, that a major restructuring of the Southern Ocean plankton ecosystem occurred abruptly and concomitant with the first major Antarctic glaciation in the earliest Oligocene (~33.6 million years ago). This turnover marks a regime shift in zooplankton-phytoplankton interactions and community structure, which indicates the appearance of eutrophic and seasonally productive environments on the Antarctic margin. We conclude that earliest Oligocene cooling, ice-sheet expansion, and subsequent sea-ice formation were important drivers of biotic evolution in the Southern Ocean.

Age tie points and geochemistry of sediment cores from the Antarctic continental margin

Few high-latitude terrestrial records document the timing and nature of the Cenozoic "Greenhouse" to "Icehouse" transition. Here we exploit the bulk geochemistry of marine siliciclastic sediments from drill cores on Antarctica's continental margin to extract a unique semiquantitative temperature and precipitation record for Eocene to mid-Miocene (~54-13 Ma). Alkaline elements are strongly enriched in the detrital mineral fraction in fine-grained siliciclastic marine sediments and only occur as trace metals in the biogenic fraction. Hence, terrestrial climofunctions similar to the chemical index of alteration (CIA) can be applied to the alkaline major element geochemistry of marine sediments on continental margins in order to reconstruct changes in precipitation and temperature. We validate this approach by comparison with published paleotemperature and precipitation records derived from fossil wood, leaves, and pollen and find remarkable agreement, despite uncertainties in the calibrations of the different proxies. A long-term cooling on the order of >=8°C is observed between the Early Eocene Climatic Optimum (~54-52 Ma) and the middle Miocene (~15-13 Ma) with the onset of transient cooling episodes in the middle Eocene at ~46-45 Ma. High-latitude stratigraphic records currently exhibit insufficient temporal resolution to reconstruct continental aridity and inferred ice-sheet development during the middle to late Eocene (~45-37 Ma). However, we find an abrupt aridification of East Antarctica near the Eocene-Oligocene transition (~34 Ma), which suggests that ice coverage influenced high-latitude atmospheric circulation patterns through albedo effects from the earliest Oligocene onward.

EPICA Dome C Ice Core 800KYr deuterium data and temperature estimates

A high-resolution deuterium profile is now available along the entire European Project for Ice Coring in Antarctica Dome C ice core, extending this climate record back to marine isotope stage 20.2, ~800,000 years ago. Experiments performed with an atmospheric general circulation model including water isotopes support its temperature interpretation. We assessed the general correspondence between Dansgaard-Oeschger events and their smoothed Antarctic counterparts for this Dome C record, which reveals the presence of such features with similar amplitudes during previous glacial periods. We suggest that the interplay between obliquity and precession accounts for the variable intensity of interglacial periods in ice core records. Temperature was estimated after correction for sea-water isotopic composition (Bintanja et al, 2005) and for ice sheet elevation (Parrenin et al, 2007) on EDC3 age scale (Parrenin et al, 2007).

Atmospheric carbon dioxide, methane, deuterium, and calculated Antarctic temperature of EPICA Dome C ice core

Understanding the role of atmospheric CO2 during past climate changes requires clear knowledge of how it varies in time relative to temperature. Antarctic ice cores preserve highly resolved records of atmospheric CO2 and Antarctic temperature for the past 800,000 years. Here we propose a revised relative age scale for the concentration of atmospheric CO2 and Antarctic temperature for the last deglacial warming, using data from five Antarctic ice cores. We infer the phasing between CO2 concentration and Antarctic temperature at four times when their trends change abruptly. We find no significant asynchrony between them, indicating that Antarctic temperature did not begin to rise hundreds of years before the concentration of atmospheric CO2, as has been suggested by earlier studies.

Characteristics of ice-rafted pebbles from ODP Leg 178 holes

Pebbles (>10 mm) sampled from three drill sites on the continental rise west of the Antarctic Peninsula during Ocean Drilling Program Leg 178 were classified by shape and roundness. In addition, pebble lithology and surface texture were visually identified. To increase the pebble sample number to 331, three sites that were drilled 94 to 213 km from the continental shelf edge were integrated into the data set using magnetostratigraphy for core correlation. Pebbles were compared in three groups defined by the same stratigraphic intervals at each site: 3.1-2.2 Ma (late Pliocene), 2.2-0.76 Ma (late Pliocene-late Pleistocene), and 0.76 Ma to the Holocene. Pebble lithologies originate from sources on the Antarctic Peninsula margin. Most pebbles are metamorphic and sedimentary pebbles are rare (<6%), whereas mafic volcanic and intrusive igneous lithologies increase in abundance upsection. Pebbles from 3.1 to 0.76 Ma, plotted on sphericity-roundness diagrams, indicate original transport as basal and supraglacial/englacial debris. Pebbles are abundant and of diverse lithology. From 0.76 Ma to the present, the number of pebbles is low and their shape characteristics indicate they originated as basal debris. Observed changes in ice-rafted pebbles can be explained by growth of an ice sheet and inundation of the Antarctic Peninsula topography by ice ~0.76 Ma. Prior to this, outlet and valley glaciers transported debris at high levels within and at the base of the ice. The mass accumulation rate of sand fluctuates and includes rounded quartz grains. Ice-sheet growth may have been accompanied by overall cooling from subpolar to polar glacial regimes, which halted meltwater production and enhanced the growth of ice shelves, which consequently reduced sediment supply to icebergs.

(Table 1) Ice thickness of new subglacial lakes identified in the ICECAP RES data and of previously known subglacial lakes

Subglacial hydrology in East Antarctica is poorly understood, yet may be critical to the manner in which ice flows. Data from a new regional airborne geophysical survey (ICECAP) have transformed our understanding of the topography and glaciology associated with the 287,000 km**2 Aurora Subglacial Basin in East Antarctica. Using these data, in conjunction with numerical ice sheet modeling, we present a suite of analyses that demonstrate the potential of the 1000 km-long basin as a route for subglacial water drainage from the ice sheet interior to the ice sheet margin. We present results from our analysis of basal topography, bed roughness and radar power reflectance and from our modeling of ice sheet flow and basal ice temperatures. Although no clear-cut subglacial lakes are found within the Aurora Basin itself, dozens of lake-like reflectors are observed that, in conjunction with other results reported here, support the hypothesis that the basin acts as a pathway allowing discharge from subglacial lakes near the Dome C ice divide to reach the coast via the Totten Glacier.

Sedimentological and paleomagnetic investigations on two gravity cores and ODP sites 113-690 and 693 off Kapp Norvegia, Antarctic continental margin

Oligocene to Quaternary sediments were recovered from the Antarctic continental margin in the eastern Weddell Sea during ODP Leg 113 and Polarstern expedition ANT-VI. Clay mineral composition and grain size distribution patterns are useful for distinguishing sediments that have been transported by ocean currents from those that were ice-rafted. This, in turn, has assisted in providing insights about the changing late Paleogene to Neogene sedimentary environment as the cryosphere developed in Antarctica. During the middle Oligocene, increasing glacial conditions on the continent are indicated by the presence of glauconite sands, that are interpreted to have formed on the shelf and then transported down the continental slope by advancing glaciers or as a result of sea-level lowering. The dominance of illite and a relatively high content of chlorite suggest predominantly physical weathering conditions on the continent. The high content of biogenic opal from the late Miocene to the late Pliocene resulted from increased upwelling processes at the continental margin due to increased wind strength related to global cooling. Partial melting of the ice-sheet occurred during an early Pliocene climate optimum as is shown by an increasing supply of predominantly current-derived sediment with a low mean grain size and peak values of smectite. Primary productivity decreased at ~ 3 Ma due to the development of a permanent sea-ice cover close to the continent. Late Pleistocene sediments are characterized by planktonic foraminifers and biogenic opal, concentrated in distinct horizons reflecting climatic cycles. Isotopic analysis of AT. pachyderma produced a stratigraphy which resulted in a calculated sedimentation rate of 1 cm/k.y. during the Pleistocene. Primary productivity was highest during the last three interglacial maxima and decreased during glacial episodes as a result of increasing sea-ice coverage.

Relative sea-level history of Marguerite Bay, Antarctic Peninsula derived from optically stimulated luminescence-dated beach cobbles

Calmette Bay within Marguerite Bay along the western side of the Antarctic Peninsula contains one of the most continuous flights of raised beaches described to date in Antarctica. Raised beaches extend to 40.8 m above sea level (masl) and are thought to reflect glacial isostatic adjustment due to the retreat of the Antarctic Peninsula Ice Sheet. Using optically stimulated luminescence (OSL), we dated quartz extracts from cobble surfaces buried in raised beaches at Calmette Bay. The beaches are separated into upper and lower beaches based on OSL ages, geomorphology, and sedimentary fabric. The two sets of beaches are separated by a prominent scarp. One of our OSL ages from the upper beaches dates to 9.3 thousand years ago (ka; as of 1950) consistent with previous extrapolation of sea-level data and the time of ice retreat from inner Marguerite Bay. However, four of the seven ages from the upper beaches date to the timing of glaciation. We interpret these ages to represent reworking of beaches deposited prior to the Last Glacial Maximum (LGM) by advancing and retreating LGM ice. Ages from the lower beaches record relative sea-level fall due to Holocene glacial-isostatic adjustment. We suggest a Holocene marine limit of 21.7 masl with an age of 5.5-7.3 ka based on OSL ages from Calmette Bay and other sea-level constraints in the area. A marine limit at 21.7 masl implies half as much relative sea-level change in Marguerite Bay during the Holocene as suggested by previous sea-level reconstructions. No evidence for a relative sea-level signature of neoglacial events, such as a decrease followed by an increase in RSL fall due to ice advance and retreat associated with the Little Ice Age, is found within Marguerite Bay indicating either: (1) no significant neoglacial advances occurred within Marguerite Bay; (2) rheological heterogeneity allows part of the Antarctic Peninsula (i.e. the South Shetland Islands) to respond to rapid ice mass changes while other regions are incapable of responding to short-lived ice advances; or (3) the magnitude of neoglacial events within Marguerite Bay is too small to resolve through relative sea-level reconstructions. Although the application of reconstructing sea-level histories using OSL-dated raised beach deposits provides a better understanding of the timing and nature of relative sea-level change in Marguerite Bay, we highlight possible problems associated with using raised beaches as sea-level indices due to post-depositional reworking by storm waves.

(Table S1) Diatom abundance according to ice cycles in ANDRILL AND1-1B drill core

The influence of Antarctica and the Southern Ocean on Late Pliocene global climate reconstructions has remained ambiguous due to a lack of well-dated Antarctic-proximal, paleoenvironmental records. Here we present ice sheet, sea-surface temperature, and sea ice reconstructions from the ANDRILL AND-1B sediment core recovered from beneath the Ross Ice Shelf. We provide evidence for a major expansion of an ice sheet in the Ross Sea that began at ~3.3 Ma, followed by a coastal sea surface temperature cooling of ~2.5°C, a stepwise expansion of sea ice, and polynya-style deep mixing in the Ross Sea between 3.3 and 2.5 Ma. The intensification of Antarctic cooling resulted in strengthened westerly winds and invigorated ocean circulation. The associated northward migration of Southern Ocean fronts has been linked with reduced Atlantic Meridional Overturning Circulation by restricting surface water connectivity between the ocean basins, with implications for heat transport to the high latitudes of the North Atlantic. While our results do not exclude low-latitude mechanisms as drivers for Pliocene cooling, they indicate an additional role played by southern high-latitude cooling during development of the bipolar world.

Cenozoic global ice volume and temperature simulations over the past 40 million years

Variations in global ice volume and temperature over the Cenozoic era have been investigated with a set of one-dimensional (1-D) ice-sheet models. Simulations include three ice sheets representing glaciation in the Northern Hemisphere, i.e. in Eurasia, North America and Greenland, and two separate ice sheets for Antarctic glaciation. The continental mean Northern Hemisphere surface-air temperature has been derived through an inverse procedure from observed benthic d18O records. These data have yielded a mutually consistent and continuous record of temperature, global ice volume and benthic d18O over the past 35 Ma. The simple 1-D model shows good agreement with a comprehensive 3-D ice-sheet model for the past 3 Ma. On average, differences are only 1.0°C for temperature and 6.2 m for sea level. Most notably, over the 35 Ma period, the reconstructed ice volume-temperature sensitivity shows a transition from a climate controlled by Southern Hemisphere ice sheets to one controlled by Northern Hemisphere ice sheets. Although the transient behaviour is important, equilibrium experiments show that the relationship between temperature and sea level is linear and symmetric, providing limited evidence for hysteresis. Furthermore, the results show a good comparison with other simulations of Antarctic ice volume and observed sea level.