(Supplement 3) Radar sounding observations including bedrock and surface elevation and ice thickness in the ASB region, Antarctica

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 under¬standing 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 ice- flow regimes of different scales and directions, each of which would have over-deepened existing topographic depressions, reversing valley floor slopes.

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

Marine ice-rafted debris records and age determinations along the northern Eurasian margin

Ice-rafted debris (IRD) (>2 mm), input in eight sediment cores along the Eurasian continental margin (Arctic Ocean), have been studied over the last two glacial/interglacial cycles. Together with the revised chronologies and new micropaleontological data of two cores from the northern Barents Sea (PS2138) and northeastern Kara Sea (PS2741) spanning Marine Isotope Stages (MIS) 6 to 1, the IRD data give new insights into the glacial history of northern Eurasian ice-sheets over the last 150 ka. The chronologies of the cores are based on stable isotope records, AMS 14C datings, paleomagnetic and biostratigraphic data. Extensive episodes of northern Barents Sea ice-sheet growth, probably to the shelf edge, occurred during the late Weichselian (MIS 2) and the Saalian (MIS 6). Major IRD discharge at the MIS 4/3-transition hints to another severe glaciation, probably onto the outer shelf, during MIS 4. IRD-based instabilities of the marine-based ice margin along the northern Barents Sea between MIS 4 and 2 are similar in timing with North Atlantic Heinrich events and Nordic Seas IRD events, suggesting similar atmospheric cooling over a broad region or linkage of ice-sheet fluctuations through small sea-level events. In the relatively low-precipitation areas of eastern Eurasia, IRD peak values during Termination II and MIS 4/3-transition suggest a Kara Sea ice-sheet advance onto the outer shelf, probably to the shelf edge, during glacial MIS 6 and 4. This suggests that during the initial cooling following the interglacials MIS 5, and possibly MIS 7, the combined effect of sustained inflow of Atlantic water into the Arctic Ocean and penetration of moisture-bearing cyclones into easterly direction supported major ice build-up during Saalian (MIS 6) and Mid-Weichselian (MIS 4) glaciation. IRD peak values in MIS 5 indicate at least two advances of the Severnaya Semlya ice-sheet to the coast line during the Early Weichselian. In contrast, a distinct Kara Sea ice advance during the Late Weichselian (MIS 2) is not documented by the IRD records along the northeastern Kara Sea margin.

(Table 1) Grain size, sedimentation rates and accumulation rates of ODP Hole 152-919A coarse-sand ice-rafted debris

At mid- to high-latitude marine sites, ice-rafted debris (IRD) is commonly recognized as anomalously coarse-grained terrigenous material contained within a fine-grained hemipelagic or pelagic matrix (e.g., Conolly and Ewing, 1970; Ruddiman, 1977, doi:10.1130/0016-7606(1977)88<1813:LQDOIS>2.0.CO;2; Krissek, 1989, doi:10.2973/odp.proc.sr.104.114.1989; Jansen et al., 1990; Bond et al., doi:10.1038/360245a0, 1992; Krissek, 1995, doi:10.2973/odp.proc.sr.145.118.1995). The presence of such ice-rafted material is a valuable indicator of the presence of glacial ice at sea level on an adjacent continent, whereas the composition of the IRD can often be used to identify the location of the source area (e.g., Goldschmidt, 1995, doi:10.1016/0025-3227(95)00098-J). Because the amount of core recovered during Leg 163 was very limited, this shore-based, postcruise study focuses on materials recovered at a nearby site during Leg 152. In particular, this study examines sediments recovered at Site 919; these sediments were described as containing a significant ice-rafted component in the Leg 152 Initial Reports volume (Larsen, Saunders, Clift, et al., 1994, doi:10.2973/odp.proc.ir.152.1994). In this study, the sedimentary section from Site 919 has been examined with the goal of providing a detailed history of glaciations on Greenland and other landmasses adjacent to the Norwegian-Greenland Sea; this history ultimately will be calibrated using an oxygen isotope stratigraphy (Flower, 1998, doi:10.2973/odp.proc.sr.152.219.1998), although that calibration has not been completed at this time. Because ice-core studies of the Greenland Ice Sheet (GIS) have shown that the GIS changed dramatically, and in some cases extremely rapidly, during at least the last interglacial stage (GRIP Members, 1993, doi:10.1038/364203a0), a detailed IRD record from the Southeast Greenland margin should provide insight into the longer term behavior of this sensitive component of the Northern Hemisphere climate system.

(Table S1) Brominated flame retardant concentrations in an ice core from Holtedahlfonna, Svalbard

Brominated flame retardants (BFRs) have been found in Arctic wildlife, lake sediment, and air. To identify the atmospheric BFR deposition history on Svalbard, Norway, we analyzed 19 BFRs, including hexabromocyclododecane (HBCD), 1,2-bis(2,4,6-tribromophenoxy)ethane (BTBPE), decabromodiphenyl ethane (DBDPE), pentabromoethylbenzene (PBEB),and 15 polybrominated diphenyl ether congeners (PBDE) in the upper 34 m of an ice core (representing 1953-2005) from Holtedahlfonna, the western-most ice sheet on Svalbard. All of the non-PBDE compounds were detected in nearly continuous profiles in the core. Seven PBDEs were not observed above background (28,47,66,100,99,154,153), while 4 were found in 1 or 2 of 6 segments (17,85,138,183). BDEs-49,71,190,209 had nearly continuous profiles but only BDE-209 in large amounts. The greatest inputs were HBCD and BDE-209, 910, and 320 pg/cm**2/yr from 1995-2005. DBDPE, BTBPE, and PBEB show nearly continuous input growth in recent core segments, but all were <6 pg/cm**2/yr. Long-range atmospheric processes may have moved these particle-bound BFRs to the site, probably during the Arctic haze season. Average air mass trajectories over 10 years show >75% of atmospheric flow to Holtedahlfonna coming from Eurasia during haze periods (March and April).

Sedimentology on cores off Filchner Shelf Ice, Weddell Sea, Antarctica

The sediments of 14 box cores and 7 gravity cores, mainly taken directly in front of the Filchner(-Ronne) ice shelf northwest of Berkner Island (Weddell Sea), allowed to distinguish six sediment types. On the one hand,the retreat of the at first grounded and then floated ice from the last glacial maximum is documented. On the other hand,the sediments give an insight into extensive Holocene sediment deposition and remobilization northwest of Berkner Island. The ortho till was deposited directly by the grounded ice sheet and is lacking any marine influence. After floating of the ice shelf, partly very weIl stratified, partly unstratified, non-bioturbated paratill is deposited beneath the ice shelf. Lack of IRD-content in the paratill immediately above the orthotill indicates freezing at the bottom of the ice, at least for a short period after the ice became afloat. The orthotill and paratill contain small amounts of fragmented Tertiary diatoms, which allow the conclusion, that glacial-marine sediments in the accumulation area of the Ronne ice shelf will be eroded and later deposited by ice in the investigation area. Starting of bioturbation and therefore change in sedimentation from paratill to bioturbated paratill,is caused by the retreat of the ice shelf to its actual position. Isostatic uplift of the sea-bed after the Ice Age causes minor water depths with higher current velocities. The fine-fraction is eroding and mean particle-size will increase. Maybe, also isostatic uplift is responsible for repeated great advances of the floated ice shelf as shown in an erosional horizon in some cores containing bioturbated paratill. Postglacial sediment-thicknesses exceed 3 m. Assuming floating of the ice 15.000 YBP, accumulation rates reach nearly 20cm/lOOO years. Following the theories about sediment input in front of wide ice shelves, this was not expected. In the shallower water depths of Berkner Bank, the oscillations of the ice shelf are recorded in the sediments. Sorting and redistribution by high current velocities from beneath the ice up to the calving line, lead to the deposition of the weIl to very weIl sorted sandy till. In front of the calving line the finer fraction will settle down. Remobilization is possible by bioturbation and increasing current-velocity. According to the intensity of mixing of the sandy till with the fine fraction, modified till or muddy till results.

Dinocyst-derived sea-surface characteristics, palynomorph abundances, ice-rafted debris, and planktonic foraminifera abundances and d18O of the 35-41 ka BP interval of core MD99-2285

The last glacial millennial climatic events (i.e. Dansgaard-Oeschger and Heinrich events) constitute outstanding case studies of coupled atmosphere-ocean-cryosphere interactions. Here, we investigate the evolution of sea-surface and subsurface conditions, in terms of temperature, salinity and sea ice cover, at very high-resolution (mean resolution between 55 and 155 years depending on proxies) during the 35-41 ka cal BP interval covering three Dansgaard-Oeschger cycles and including Heinrich event 4, in a new unpublished marine record, i.e. the MD99-2285 core (62.69°N; -3.57s°E). We use a large panel of complementary tools, which notably includes dinocyst-derived sea-ice cover duration quantifications. The high temporal resolution and multiproxy approach of this work allows us to identify the sequence of processes and to assess ocean-cryosphere interactions occurring during these periodic ice-sheet collapse events. Our results evidence a paradoxical hydrological scheme where (i) Greenland interstadials are marked by a homogeneous and cold upper water column, with intensive winter sea ice formation and summer sea ice melting, and (ii) Greenland and Heinrich stadials are characterized by a very warm and low saline surface layer with iceberg calving and reduced sea ice formation, separated by a strong halocline from a less warm and saltier subsurface layer. Our work also suggests that this stadial surface/subsurface warming started before massive iceberg release, in relation with warm Atlantic water advection. These findings thus support the theory that upper ocean warming might have triggered European ice-sheet destabilization. Besides, previous paleoceanographic studies conducted along the Atlantic inflow pathways close to the edge of European ice-sheets suggest that such a feature might have occurred in this whole area. Nonetheless, additional high resolution paleoreconstructions are required to confirm such a regional scheme.

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

Isotopic ratios and ion composition of massive ground ice on Herschel Island

Herschel Island in the southern Beaufort Sea is a push moraine at the northwestern-most limit of the Laurentide Ice Sheet. Stable water isotope (d18O, dD) and hydrochemical studies were applied to two tabular massive ground ice bodies to unravel their genetic origin. Buried glacier ice or basal regelation ice was encountered beneath an ice-rich diamicton with strong glaciotectonic deformation structures. The massive ice isotopic composition was highly depleted in heavy isotopes (mean d18O: -33 per mil; mean dD: -258 per mil), suggesting full-glacial conditions during ice formation. Other massive ice of unknown origin with a very large d18O range (from -39 to -21 per mil) was found adjacent to large, striated boulders. A clear freezing slope was present with progressive depletion in heavy isotopes towards the centre of the ice body. Fractionation must have taken place during closed-system freezing, possibly of a glacial meltwater pond. Both massive ground ice bodies exhibited a mixed ion composition suggestive of terrestrial waters with a marine influence. Hydrochemical signatures resemble the Herschel Island sediments that are derived from nearshore marine deposits upthrust by the Laurentide ice. A prolonged contact between water feeding the ice bodies and the surrounding sediment is therefore inferred.

Table 3.5 Content and composition of air in fresh ice

Natural ice is formed by freezing of water or by sintering of dry or wet snow. Each of these processes causes atmospheric air to be enclosed in ice as bubbles. The air amount and composition as well as the bubble sizes and density depend not only on the kind of process but also on several environmental conditions. The ice in the deepest layers of the Greenland and thc Antarctic ice sheet was formed more than 100 000 years ago. In the bubbles of this ice, samples of atmospheric air from that time are preserved. The enclosure of air is discussed for each of the three processes. Of special interest are the parameters which control the amount and composition of the enclosed air. If the ice is formed by sintering of very cold dry snow, the air composition in the bubbles corresponds with good accuracy to the composition of atmospheric air.