Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Age determination of sediment from Lake Komi

Beach and shoreface sediments deposited in the more than 800-km long ice-dammed Lake Komi in northern European Russia have been investigated and dated. The lake flooded the lowland areas between the Barents-Kara Ice Sheet in the north and the continental drainage divide in the south. Shoreline facies have been dated by 18 optical stimulated luminescence (OSL) dates, most of which are closely grouped in the range 80-100 ka, with a mean of 88 +/- 3 ka. This implies that that the Barents-Kara Ice Sheet had its Late Pleistocene maximum extension during the Early Weichselian, probably in the cold interval (Rederstall) between the Brørup and Odderade interstadials of western Europe, correlated with marine isotope stage 5b. This is in strong contrast to the Scandinavian and North American ice sheets, which had their maxima in isotope stage 2, about 20 ka. Field and air photo interpretations suggest that Lake Komi was dammed by the ice advance, which formed the Harbei-Harmon-Sopkay Moraines. These has earlier been correlated with the Markhida moraine across the Pechora River Valley and its western extension. However, OSL dates on fluvial sediments below the Markhida moraine have yielded ages as young as 60 ka. This suggests that the Russian mainland was inundated by two major ice sheet advances from the Barents-Kara seas after the last interglacial: one during the Early Weichselian (about 90 ka) that dammed Lake Komi and one during the Middle Weichselian (about 60 ka). Normal fluvial drainage prevailed during the Late Weichselian, when the ice front was located offshore.

Ice-rafted debris record in the Prydz Bay

Ice-rafted debris mass accumulation rates (IRD MAR) at a drill site on the Antarctic continental margin are investigated to evaluate the linkages between East Antarctic Ice Sheet extent and Southern Ocean temperatures in the early to mid-Pliocene. ODP Site 1165 is within 400 km of the Antarctic coastline and in the direct pathway of icebergs released by the Amery Ice Shelf. The Amery Ice Shelf is the largest ice shelf in East Antarctica and it buttresses the Lambert Glacier drainage system, which accounts for 14% of the outflow from the East Antarctic Ice Sheet. IRD MAR were low during peak Southern Ocean warming in the early Pliocene. After a brief precursor, a tenfold increase in IRD MAR at 3.3 Ma marks the termination of the early Pliocene ice sheet minimum, coincident with the M2 glacial. For the mid-Pliocene, a strong correlation exists between the high-amplitude signal in the LR04 benthic stack and IRD MAR, suggesting linkages between East Antarctic ice extent, global ice volume and deep-water temperatures. The IRD record at Site 1165 provides evidence of greater sensitivity of the Lambert Glacier-Amery Ice Shelf system to Southern Ocean warming than is currently predicted by ice sheet models, which may relate to uncertainties in the understanding of ocean heat uptake, poleward heat transport and ice sheet-ocean interactions.

Neogene ice-rafted debris record of ODP Hole 120-751A

One of the primary objectives of Leg 120 was to obtain a high-resolution Neogene stratigraphic section from the Kerguelen Plateau. Site 751, located in the central part of the Raggatt Basin on the Southern Kerguelen Plateau in 1633.8 m of water (57°43.56'S; 79°48.89'E), was selected as the dedicated Neogene site for this objective. High-resolution sampling at Site 751 was used to delineate in detail the Neogene ice-rafted debris (IRD) occurrences on the Kerguelen Plateau. The oldest IRD found at Site 751 was approximately 9.9 Ma, and it was not until approximately 8.5 Ma that significant concentrations of IRD were detected. The first major IRD event at this site occurred in the uppermost Miocene between 6.0 and 5.5 Ma. During this time period, a general climatic cooling and glacial expansion occurred on Antarctica. The late Miocene IRD event was followed by a continuous episode of elevated IRD deposition in the lowermost Pliocene between 4.5 and 4.1 Ma. The 0.4-m.y. duration and the timing of the early Pliocene IRD event on the Kerguelen Plateau corresponds with IRD fluxes observed on the Falkland Plateau and in the Weddell Abyssal Plain. This correspondence of data indicates that a major global climatic event occurred during the early Pliocene. The East Antarctic Ice Sheet may have experienced deglaciation between 4.5 and 4.1 Ma and, as a result, released large volumes of sediment-laden ice into the Southern Ocean.

EPICA Dronning Maud Land EDML ice core drilling protocol

We report on the EPICA Dronning Maud Land (East Antarctica) deep drilling operation. Starting with the scientific questions that led to the outline of the EPICA project, we introduce the setting of sister drillings at NorthGRIP and EPICA Dome C within the European ice-coring community. The progress of the drilling operation is described within the context of three parallel, deep-drilling operations, the problems that occurred and the solutions we developed. Modified procedures are described, such as the monitoring of penetration rate via cable weight rather than motor torque, and modifications to the system (e.g. closing the openings at the lower end of the outer barrel to reduce the risk of immersing the drill in highly concentrated chip suspension). Parameters of the drilling (e.g. core-break force, cutter pitch, chips balance, liquid level, core production rate and piece number) are discussed. We also review the operational mode, particularly in the context of achieved core length and piece length, which have to be optimized for drilling efficiency and core quality respectively. We conclude with recommendations addressing the design of the chip-collection openings and strictly limiting the cable-load drop with respect to the load at the start of the run.

Ice rafted debris in the Nordic Seas

A continuous 3.5 Myr IRD record was produced from Ocean Drilling Program (ODP) Site 907. A timescale based on magnetic polarity chrons, oxygen isotope stratigraphy (for the last 1Myr) and orbital tuning was developed. The record documents a stepwise inception of large-scale glacial cycles in the Nordic Seas region, the first being a marked expansion of the Greenland ice sheet at 3.3 Ma. A second step occurred at 2.74 Ma by an expansion of large scale ice sheets in the Northern Hemisphere. Ice sheet variability around the Nordic Seas was tightly coupled to global ice volume over the past 3.3 Myr. Between 3 and 1 Ma, most of the variance of the IRD signal is in the 41 kyr band, whereas the last 1 Myr is characterized by stronger 100 kyr variance. The Gamma Ray Porosity Evaluator (GRAPE) density record is closely linked with IRD variations and documents sub orbital variability resembling the late Quaternary Heinrich/Bond cycles.

(Table 1) Two-way travel times and depths of picked internal reflection horizons (IRHs) of the NGRIP ice core

During the summer of 2003, a ground-penetrating radar survey around the North Greenland Icecore Project (NorthGRIP) deep ice-core drilling site (75°06' N, 42°20' W; 2957 m a.s.l.) was carried out using a shielded 250 MHz radar system. The drill site is located on an ice divide, roughly 300 km north-northwest of the summit of the Greenland ice sheet. More than 430 km of profiles were measured, covering a 10 km by 10 km area, with a grid centered on the drilling location, and eight profiles extending beyond this grid. Seven internal horizons within the upper 120 m of the ice sheet were continuously tracked, containing the last 400 years of accumulation history. Based on the age-depth and density-depth distribution of the deep core, the internal layers have been dated and the regional and temporal distribution of accumulation rate in the vicinity of NorthGRIP has been derived. The distribution of accumulation shows a relatively smoothly increasing trend from east to west from 145 kg/m**2/a to 200 kg/m**2/a over a distance of 50 km across the ice divide. The general trend is overlain by small-scale variations on the order of 2.5 kg/m**2/a/km, i.e. around 1.5% of the accumulation mean. The temporal variations of the seven periods defined by the seven tracked isochrones are on the order of +-4% of the mean of the last 400 years, i.e. at NorthGRIP ±7 kg/m**2/a. If the regional accumulation pattern has been stable for the last several thousand years during the Holocene, and ice flow has been comparable to today, advective effects along the particle trajectory upstream of NorthGRIP do not have a significant effect on the interpretation of climatically induced changes in accumulation rates derived from the deep ice core over the last 10 kyr.

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.

Synthesized grounding line and ice shelf mask for Antarctica

Iceberg calving has been assumed to be the dominant cause of mass loss for the Antarctic ice sheet, with previous estimates of the calving flux exceeding 2,000 gigatonnes per year. More recently, the importance of melting by the ocean has been demonstrated close to the grounding line and near the calving front. So far, however, no study has reliably quantified the calving flux and the basal mass balance (the balance between accretion and ablation at the ice-sheet base) for the whole of Antarctica. The distribution of fresh water in the Southern Ocean and its partitioning between the liquid and solid phases is therefore poorly constrained. Here we estimate the mass balance components for all ice shelves in Antarctica, using satellite measurements of calving flux and grounding-line flux, modelled ice-shelf snow accumulation rates and a regional scaling that accounts for unsurveyed areas. We obtain a total calving flux of 1,321 ± 144 gigatonnes per year and a total basal mass balance of -1,454 ± 174 gigatonnes per year. This means that about half of the ice-sheet surface mass gain is lost through oceanic erosion before reaching the ice front, and the calving flux is about 34 per cent less than previous estimates derived from iceberg tracking. In addition, the fraction of mass loss due to basal processes varies from about 10 to 90 per cent between ice shelves. We find a significant positive correlation between basal mass loss and surface elevation change for ice shelves experiencing surface lowering and enhanced discharge. We suggest that basal mass loss is a valuable metric for predicting future ice-shelf vulnerability to oceanic forcing.

Strontium and Neodymium isotope composition of sediments from the North Atlantic

Sr and Nd isotopic compositions have been measured on the lithic fraction of last climatic cycle sediments from the North Atlantic (~40°N/~60°N), in order to identify the origins of the particles. From the reconstruction of their transport pathways, we deduce the mechanisms that explain their distributions. The main source regions are the Canadian shield (mostly the area of Baffin Bay and western Greenland), the Scandinavian shield, the European region (British Isles and Bay of Biscay), and Iceland. We observe a significant glacial/interglacial contrast, characterized by a dominant Icelandic input via near-bottom transport by North Atlantic Deep Water (NADW) during the interglacials and a largely continent-derived contribution of surface-transported, ice-rafted detritus (IRD) during the glacial period. During the last glacial period, the Heinrich events (abrupt, massive discharges of IRD) originated not only from the Laurentide ice sheet as heretofore envisioned but also from other sources. Three other major North Atlantic ice sheets (Fennoscandian, British Isles, and Icelandic) probably surged simultaneously, discharging ice and IRD into the North Atlantic. As opposed to theories implying a unique, Laurentide origin [Gwiazda et al., 1996 doi:10.1029/95PA03135] driven by an internal mechanism [MacAyeal, 1993 doi:10.1029/93PA02200], we confirm that the Icelandic and the Fennoscandian ice sheets also surged as recently proposed by other authors, and we here also distinguish a possible detrital contribution from the British Isles ice sheet. This pan-North Atlantic phenomenon thus requires a common regional, external forcing.

Ages and geochemical and geophysical characteristics of sediment core MD99-2294, Lofoten Contourite Drift

A sediment core from the Lofoten Contourite Drift on the continental slope off Northern Norway, proximal to the former Vestfjorden-Trsnadjupet Ice Stream, details the development, variability and decline of marine margins of the northwestern Fennoscandian Ice Sheet during the time interval 25.3-14 cal ka BP, including the Last Glacial Maximum and onset of the deglaciation based on high-resolution IRD records. From the core interval between 25.3 and 17.7 cal ka BP we report data points with a mean time step of 10 years, between 17.7 cal ka BP and the Holocene time steps are typically 50 years. The core is divided into 7 informal ice-rafted debris (IRD) zones based on the variations in IRD including 7 major IRD maxima (A-G), inferred to represent periods of high iceberg production. Petrological identification reveals dominance of crystalline IRD (monocrystalline, plutonic and metamorphic rock fragments) accounting for 75-80% of total IRD assemblages, while sedimentary fragments generally account for 15-20%. The crystalline fragments (including eclogite and mangerite from a nearby terrestrial source) increase across the IRD peaks while the sedimentary fragments remain constant. This points to the importance of erosional products from icebergs originating from fast-flowing paleo-ice streams including the Vestfjorden-Trsnadjupet Ice Stream draining from the Fennoscandian mainland during the IRD maxima periods. Increased temperature of the adjacent surface water masses was probably an important external forcing factor on the Fennoscandian Ice Sheet behavior because some IRD maxima and plumite deposition from meltwater plumes post-date periods of increased sea surface temperatures. The peak IRD depositions occur in centennial and millennial time cycles (~200, 1030 and 3900 year) indicating some external forcing by solar variation. Both mechanisms could explain the observed synchronous instability of the northwestern Fennoscandian Ice Sheet to other European Ice Sheets.