(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 of various cores from the West Antarctic continental margin

Three megascopic and disseminated tephra layers (which we refer to as layers A, B, and C) occur in late Quaternary glaciomarine sediments deposited on the West Antarctic continental margin. The stratigraphical positions of the distal tephra layers in 28 of the 32 studied sediment cores suggest their deposition during latest Marine Isotopic Stage (MIS) 6 and MIS 5. One prominent tephra layer (layer B), which was deposited subsequent to the penultimate deglaciation (Termination II), is present in almost all of the cores. Geochemical analyses carried out on the glass shards of the layers reveal a uniform trachytic composition and indicate Marie Byrd Land (MBL), West Antarctica, as the common volcanic source. The geochemical composition of the marine tephra is compared to that of ash layers of similar age described from Mount Moulton and Mount Takahe in MBL and from ice cores drilled at Dome Fuji, Vostok and EPICA Dome C in East Antarctica. The three tephra layers in the marine sediments are chemically indistinguishable. Also five englacial ash layers from Mt. Moulton, which originated from highly explosive Plinian eruptions of the Mt. Berlin volcano in MBL between 142 ka and 92 ka ago, are chemically very similar, as are two tephra layers erupted from Mt. Takahe at ca. 102 ka and ca. 93 ka. Statistical analysis of the chemical composition of the glass shards indicates that the youngest tephra (layer A) in the marine cores matches the ash layer erupted from Mt. Berlin at 92 ka, which was previously correlated with tephra layers in the EPICA Dome C and the Dome Fuji ice cores. A tephra erupted from Mt. Berlin at 136 ka seems to correspond to a tephra layer deposited at 1733 m in the EPICA Dome C ice core. Additionally, the oldest tephra (layer C) in the marine sediments resembles an ash layer deposited at Vostok around 142 ka, but statistical evidence for the validity of this correlation is inconclusive. Although our results underscore the potential of tephrostratigraphy for correlating terrestrial and marine palaeoclimate archives, our study also reveals limitations of this technique, which may result in the miscorrelation of tephra. Such pitfalls comprise failure to recognise the occurrence of various tephra layers in marine sediment cores, 'swamping' of records with chemically indistinguishable tephra from a single volcanic source, and exclusive use of 'geochemical fingerprinting' for correlating ash layers.

Water temperature reconstructions based on paired Uk'37 and Tex86 records, and paired surface and subsurface Mg/Ca records

This dataset contains the collection of available published paired Uk'37 and Tex86 records spanning multi-millennial to multi-million year time scales, as well as a collection of Mg/Ca-derived temperatures measured in parallel on surface and subsurface dwelling foraminifera, both used in the analyses of Ho and Laepple, Nature Geoscience 2016. As the signal-to-noise ratios of proxy-derived Holocene temperatures are relatively low, we selected records that contain at least the last deglaciation (oldest sample >18kyr BP).

(Table 2) Radiocarbon dates from lakes on Fildes Peninsula, King George Island

Precise relative sea level (RSL) data are important for inferring regional ice sheet histories, as well as helping to validate numerical models of ice sheet evolution and glacial isostatic adjustment. Here we develop a new RSL curve for Fildes Peninsula, South Shetland Islands (SSIs), a sub-Antarctic archipelago peripheral to the northern Antarctic Peninsula ice sheet, by integrating sedimentary evidence from isolation basins with geomorphological evidence from raised beaches. This combined approach yields not only a Holocene RSL curve, but also the spatial pattern of how RSL change varied across the archipelago. The curve shows a mid-Holocene RSL highstand on Fildes Peninsula at 15.5 m above mean sea level between 8000 and 7000 cal a BP. Subsequently RSL gradually fell as a consequence of isostatic uplift in response to regional deglaciation. We propose that isostatic uplift occurred at a non-steady rate, with a temporary pause in ice retreat ca. 7200 cal a BP, leading to a short-lived RSL rise of ~1 m and forming a second peak to the mid-Holocene highstand. Two independent approaches were taken to constrain the long-term tectonic uplift rate of the SSIs at 0.22-0.48 m/ka, placing the tectonic contribution to the reconstructed RSL highstand between 1.4 and 2.9 m. Finally, we make comparisons to predictions from three global sea level models.

Follow-up of the sponge and asteroid populations at the station Larsen A south (PS69/724-1 and PS77/253-1), Antarctic Peninsula, between 2007 and 2011

Over 30% of the Antarctic continental shelf is permanently covered by floating ice shelves, providing aphotic conditions for a depauperate fauna sustained by laterally advected food. In much of the remaining Antarctic shallows (<300 m depth), seasonal sea-ice melting allows a patchy primary production supporting rich megabenthic communities dominated by glass sponges (Porifera, Hexactinellida). The catastrophic collapse of ice shelves due to rapid regional warming along the Antarctic Peninsula in recent decades has exposed over 23,000 km**2 of seafloor to local primary production. The response of the benthos to this unprecedented flux of food is, however, still unknown. In 2007, 12 years after disintegration of the Larsen A ice shelf, a first biological survey interpreted the presence of hexactinellids as remnants of a former under-ice fauna with deep-sea characteristics. Four years later, we revisited the original transect, finding 2- and 3-fold increases in glass sponge biomass and abundance, respectively, after only two favorable growth periods. Our findings, along with other long-term studies, suggest that Antarctic hexactinellids, locked in arrested growth for decades, may undergo boom-and-bust cycles, allowing them to quickly colonize new habitats. The cues triggering growth and reproduction in Antarctic glass sponges remain enigmatic.

Investigations on sediment core Co1010 from Rauer Group, Antarctica

The Rauer Group is an archipelago in Prydz Bay, East Antarctica. The ice-free islands and the surrounding shallow marine areas provide valuable archives for the reconstruction of the late Pleistocene and Holocene environmental and climatic history of the region. Two sediment records from two marine inlets of Rauer Group have been studied for their sedimentological, geochemical, and biological characteristics. Radiocarbon ages from one of the inlets indicate ice-free conditions within the last glacial cycle, probably during the second half of Marine Isotope Stage 3. Subsequent ice sheet coverage of Rauer Group during the Last Glacial Maxiumum (LGM) can be inferred from a till layer recovered in one of the basins. The inlets became ice-free prior to 11,200 cal yr BP, when biogenic sedimentation started. Deglacial processes in the catchments, however, influenced the inlets until ~9200 cal. yr BP as evidenced by the input of minerogenic material. Marine productivity under relatively open water conditions indicates an early Holocene climate optimum until 8200 cal. yr BP, which is followed by a cooler period with increased sea ice. Warmer conditions are inferred for the mid Holocene, when both basins experienced an input of freshwater between ~5700-3500 cal. yr BP, probably due to ice-sheet melting and increased precipitation on the islands. Neoglacial cooling in the late Holocene since c. 3500 cal yr BP is reflected by an increase in sea ice in both inlets.

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.

Strontium and neodymium concentrations and isotopic compositions of Arctic marine sediments

Sr and Nd isotopic compositions of Arctic marine sediments characterize changes of sediment source regions and trace shelf-ocean particle pathways during glacial-interglacial transitions in the eastern Arctic Ocean. In the 140-ka sedimentary record of a marine core from Yermak Plateau, north of Svalbard, 87Sr/86Sr ratios and epsion-Nd values vary between 0.717 and 0.740 and 39.3 and 314.9, respectively. Sr and Nd isotopic composition both change characteristically during glacial-interglacial cycles and are correlated with the extension of the Svalbard/Barents Sea ice sheet (SBIS). The downcore variation in Sr and Nd isotopic composition indicates climatically induced changes in sediment provenance from two isotopically distinct end-members: (1) Eurasian shelf sediments as a distal source; and (2) Svalbard bedrock as a proximal source that coincide with a change in transport mechanism from sea ice to glacial ice. During glacier advance from Svalbard and intensified glacial bedrock erosion, epsion-Nd values decrease gradually to a minimum value of 314.9 due to increased input of crystalline Svalbard bedrock material. During glacial maxima, the SBIS covered the entire Barents Sea shelf and supplied increasing amounts of Eurasian shelf material to the Arctic Ocean as ice rafted detritus (IRD). Epsion-Nd values in glacial sediments reach maximum values that are comparable to the average value of modern Eurasian shelf and sea ice sediments (epsion-Nd = 310.3). This confirms ice rafting as a major sediment transport mechanism for Eurasian shelf sediments into the Arctic Ocean and trace a sediment origin from the Kara Sea/Laptev Sea shelf area. After the decay of the shelf-based SBIS, the glacial shelf sediment spikes during glacial terminations I (epsion-Nd = 310.6) and II (epsion-Nd = 310.1) epsion-Nd values rapidly decrease to values of 312.5 typical for interglacial averages. The downcore Sr isotopic composition is anticorrelated to the Nd isotopic composition, but may be also influenced by grain-size effects. In contrast, the Nd isotopic composition in clay- to silt-size fractions of one bulk sediment sample is similar to within 0.3-0.8 epsion-Nd units and seems to be a grain-size independent provenance tracer.

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.

Sea ice and wind variability during the Holocene in East Antarctica

A high-resolution continuous record of diatom census counts and diatom specific biomarkers in sediment core NBP0101-JPC24 allows assessment of oceanographic and environmental conditions in eastern Prydz Bay during the deglaciation (11 100-9000 cal yr BP) at decadal timescale. Our study improves previous snapshots investigations based on resin-embedded thin sections and presents a new proxy that compliments the diatom census counts. Our results suggest that the ice sheet retreat over the core site is dated at ~11 100 cal yr BP, setting the onset of local deglaciation and subsequent open marine conditions. The glacial retreat in Prydz Bay is due to global warming initiated at 18 cal ka BP and the regional development of the Prydz Bay cyclonic gyre. Our results further demonstrate that the deglaciation in eastern Prydz Bay can be separated in four phases: the first between 11 100 and 10 900 cal yr BP when the ice shelf was proximal and sea ice was almost perennial; the second and the third phases between 10 900-10 400 cal yr BP and 10 400-9900 cal yr BP, respectively, when the ice shelf retreated and seasonal sea ice cycle consequently developed promoting warmer water to pump into the bay within the gyre, which in turn forced the ice shelf recession and the yearly sea ice cycle establishment; and the fourth between 9900 and 9000 cal yr BP when Holocene condition were set with a recurrent seasonal sea ice cycle and a well established Prydz Bay gyre.

Pre-glacial to glacial sediment thickness grids for the Southern Pacific Margin of West Antarctica, NetCDF files

Circum-Antarctic sediment thickness grids provide constraints for basin evolution and paleotopographic reconstructions, which are important for paleo-ice sheet formation histories. By compiling old and new seismic data, we identify sequences representing pre-glacial, transitional and full glacial deposition processes along the Pacific margin of West Antarctica. The pre-glacial sediment grid depicts 1.3 to 4.0 km thick depocenters, relatively evenly distributed along the margin. The depocenters change markedly in the transitional phase at, or after, the Eocene/Oligocene boundary, when the first major ice sheets reached the shelf. Full glacial sequences, starting in the middle Miocene, indicate new depocenter formation North of the Amundsen Sea Embayment and localized eastward shifts in the Bellingshausen Sea and Antarctic Peninsula basins. Using present-day drainage paths and source areas on the continent, our calculations indicate an estimated observed total sedimentary volume of ~10 x 10**6 km**3 was eroded from West Antarctica since the separation of New Zealand in the Late Cretaceous. Of this 4.9 x 10**6 km**3 predates the onset of glaciation and need to be considered for a paleotopography reconstruction of 34 Ma. Whereas 5.1 x 10**6 km**3 postdate the onset of glaciation, of which 2.5 x 10**6 km**3 were deposited in post mid-Miocene full glacial conditions.

Mass balances of different sectors in Antarctica in 1996, 2000 and 2006

Large uncertainties remain in the current and future contribution to sea level rise from Antarctica. Climate warming may increase snowfall in the continent's interior, but enhance glacier discharge at the coast where warmer air and ocean temperatures erode the buttressing ice shelves. Here, we use satellite interferometric synthetic-aperture radar observations from 1992 to 2006 covering 85% of Antarctica's coastline to estimate the total mass flux into the ocean. We compare the mass fluxes from large drainage basin units with interior snow accumulation calculated from a regional atmospheric climate model for 1980 to 2004. In East Antarctica, small glacier losses in Wilkes Land and glacier gains at the mouths of the Filchner and Ross ice shelves combine to a near-zero loss of 4 ± 61 Gt/yr. In West Antarctica, widespread losses along the Bellingshausen and Amundsen seas increased the ice sheet loss by 59% in 10 years to reach 132 ± 60 Gt/yr in 2006. In the Peninsula, losses increased by 140% to reach 60 ± 46 Gt/yr in 2006. Losses are concentrated along narrow channels occupied by outlet glaciers and are caused by ongoing and past glacier acceleration. Changes in glacier flow therefore have a significant, if not dominant impact on ice sheet mass balance.