Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2013S1, deployed during Antarctic Fast Ice Network 2013 (AFIN 2013)

Snow height was measured by the Snow Depth Buoy 2013S1, an autonomous platform, installed close to Neumayer III Base, Antarctic during Antarctic Fast Ice Network 2013 (AFIN 2013). The resulting time series describes the evolution of snow height as a function of place and time between 2013-02-11 and 2013-04-29 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on the ice shelf. In addition to snow height, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow height occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow height may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow height, because no initial readings of absolute snow height are available.

Snow height on sea ice and sea ice drift from autonomous measurements from buoy 2014S24, deployed during Antarctic Fast Ice Network 2014 (AFIN 2014)

Snow height was measured by the Snow Depth Buoy 2014S24, an autonomous platform, installed close to Neumayer III Base, Antarctic during Antarctic Fast Ice Network 2014 (AFIN 2014). The resulting time series describes the evolution of snow depth as a function of place and time between 2014-03-07 and 2014-05-16 in sample intervals of 1 hour. The Snow Depth Buoy consists of four independent sonar measurements representing the area (approx. 10 m**2) around the buoy. The buoy was installed on the ice shelf. In addition to snow depth, geographic position (GPS), barometric pressure, air temperature, and ice surface temperature were measured. Negative values of snow depth occur if surface ablation continues into the sea ice. Thus, these measurements describe the position of the sea ice surface relative to the original snow-ice interface. Differences between single sensors indicate small-scale variability of the snow pack around the buoy. The data set has been processed, including the removal of obvious inconsistencies (missing values). Records without any snow depth may still be used for sea ice drift analyses. Note: This data set contains only relative changes in snow depth, because no initial readings of absolute snow depth are available.

Surface velocity fields, digital elevation models and ice front positions derived from multi-mission SAR remote sensing data at Sjögren Inlet glaciers, Antarctic Peninsula

Substantial retreat or disintegration of numerous ice shelves have been observed on the Antarctic Peninsula. The ice shelf in the Prince Gustav Channel retreated gradually since the late 1980's and broke-up in 1995. Tributary glaciers reacted with speed-up, surface lowering and increased ice discharge, consequently contributing to sea level rise. We present a detailed long-term study (1993-2014) on the dynamic response of Sjögren Inlet glaciers to the disintegration of Prince Gustav Ice Shelf. We analyzed various remote sensing datasets to observe the reactions of the glaciers to the loss of the buttressing ice shelf. A strong increase in ice surface velocities was observed with maximum flow speeds reaching 2.82±0.48 m/d in 2007 and 1.50±0.32 m/d in 2004 at Sjögren and Boydell glaciers respectively. Subsequently, the flow velocities decelerated, however in late 2014, we still measured about two times the values of our first measurements in 1996. The tributary glaciers retreated 61.7±3.1 km² behind the former grounding line of the ice shelf. In regions below 1000 m a.s.l., a mean surface lowering of -68±10 m (-3.1 m/a) was observed in the period 1993-2014. The lowering rate decreased to -2.2 m/a in recent years. Based on the surface lowering rates, geodetic mass balances of the glaciers were derived for different time steps. High mass loss rate of -1.21±0.36 Gt/a was found in the earliest period (1993-2001). Due to the dynamic adjustments of the glaciers to the new boundary conditions the ice mass loss reduced to -0.59±0.11 Gt/a in the period 2012-2014, resulting in an average mass loss rate of -0.89±0.16 Gt/a (1993-2014). Including the retreat of the ice front and grounding line, a total mass change of -38.5±7.7 Gt and a contribution to sea level rise of 0.061±0.013 mm were computed. Analysis of the ice flux revealed that available bedrock elevation estimates at Sjögren Inlet are too shallow and are the major uncertainty in ice flux computations. This temporally dense time series analysis of Sjögren Inlet glaciers shows that the adjustments of tributary glaciers to ice shelf disintegration are still going on and provides detailed information of the changes in glacier dynamics.

Sedimentological, mineralogical, biogeochemical and micropaleontological parameters on sediment core PS69/849-2, Burton Basin, MacRobertson Shelf, East Antarctica

A multi-proxy study including sedimentological, mineralogical, biogeochemical and micropaleontological methods was conducted on sediment core PS69/849-2 retrieved from Burton Basin, MacRobertson Shelf, East Antarctica. The goal of this study was to depict the deglacial and Holocene environmental history of the MacRobertson Land-Prydz Bay region. A special focus was put on the timing of ice-sheet retreat and the variability of bottom-water formation due to sea ice formation through the Holocene. Results from site PS69/849-2 provide the first paleo-environmental record of Holocene variations in bottom-water production probably associated to the Cape Darnley polynya, which is the second largest polynya in the Antarctic. Methods included end-member modeling of laser-derived high-resolution grain size data to reconstruct the depositional regimes and bottom-water activity. The provenance of current-derived and ice-transported material was reconstructed using clay-mineral and heavy-mineral analysis. Conclusions on biogenic production were drawn by determination of biogenic opal and total organic carbon. It was found that the ice shelf front started to retreat from the site around 12.8 ka BP. This coincides with results from other records in Prydz Bay and suggests warming during the early Holocene optimum next to global sea level rise as the main trigger. Ice-rafted debris was then supplied to the site until 5.5 cal. ka BP, when Holocene global sea level rise stabilized and glacial isostatic rebound on MacRobertson Land commenced. Throughout the Holocene, three episodes of enhanced bottom-water activity probably due to elevated brine rejection in Cape Darnley polynya occured between 11.5 and 9 cal. ka BP, 5.6 and 4.5 cal. ka BP and since 1.5 cal. ka BP. These periods are related to shifts from warmer to cooler conditions at the end of Holocene warm periods, in particular the early Holocene optimum, the mid-Holocene warm period and at the beginning of the neoglacial. In contrast, between 7.7 and 6.7 cal. ka BP, brine rejection shut down, maybe owed to warm conditions and pronounced open-water intervals.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Sedimentology of cores from the shelf of Lazarev Sea, East Antarctica

Distinct facies types, classified in radiocarbon-dated sediments from the shelf of the Lazarev Sea, East Antarctica, reveal a detailed history of processes that have controlled sedimentation during the deglaciation over the last 10,000 yr. The ice retreat on this part of the Antarctic shelf started 9500 yr BP, marked by the deposition of laminated sediments, deposited from a floating ice shelf. These laminites, which occur on top of diamictons laid down from a grounded ice sheet, are the basal sediments of the postglacial sequence. The intensity of the Antarctic Coastal Current (ACC), directed by shelf morphology, controlled sedimentation of the postglacial facies. A residual glaciomarine sediment with the fine fraction winnowed by strong currents developed from 9000-8000 yr BP in the western part of the investigation area and from 9000-5000 yr BP in the eastern part, closer to the prominent 'Fenno Deep' trough. Current velocities apparently decreased between 8000 and 2000 yr BP due to a deflection of the ACC by advancing ice tongues to the east of the investigation area during the 'Hypsithermal'. This led to a deposition of fine-grained sediments, and clay mineralogy suggests a continental source, possibly near the grounding line of the Nivl Ice Shelf, rather than a winnowing of sediments near the shelf break or advection from deeper water. Current velocities intensified after 2000 yr BP, removed fine material from these sediments and led to a relict sediment, consisting of coarse bryozoan and molluscan debris.

Surface velocity fields, digital elevation models, ice front positions and grounding line derived from remote sensing data at Dinsmoor-Bombardier-Edgeworth glacier system, Antarctic Peninsula (1992-2014)

The northern Antarctic Peninsula is one of the fastest changing regions on Earth. The disintegration of the Larsen-A Ice Shelf in 1995 caused tributary glaciers to adjust by speeding up, surface lowering, and overall increased ice-mass discharge. In this study, we investigate the temporal variation of these changes at the Dinsmoor-Bombardier-Edgeworth glacier system by analyzing dense time series from various spaceborne and airborne Earth observation missions. Precollapse ice shelf conditions and subsequent adjustments through 2014 were covered. Our results show a response of the glacier system some months after the breakup, reaching maximum surface velocities at the glacier front of up to 8.8 m/d in 1999 and a subsequent decrease to ~1.5 m/d in 2014. Using a dense time series of interferometrically derived TanDEM-X digital elevation models and photogrammetric data, an exponential function was fitted for the decrease in surface elevation. Elevation changes in areas below 1000 m a.s.l. amounted to at least 130±15 m130±15 m between 1995 and 2014, with change rates of ~3.15 m/a between 2003 and 2008. Current change rates (2010-2014) are in the range of 1.7 m/a. Mass imbalances were computed with different scenarios of boundary conditions. The most plausible results amount to -40.7±3.9 Gt-40.7±3.9 Gt. The contribution to sea level rise was estimated to be 18.8±1.8 Gt18.8±1.8 Gt, corresponding to a 0.052±0.005 mm0.052±0.005 mm sea level equivalent, for the period 1995-2014. Our analysis and scenario considerations revealed that major uncertainties still exist due to insufficiently accurate ice-thickness information. The second largest uncertainty in the computations was the glacier surface mass balance, which is still poorly known. Our time series analysis facilitates an improved comparison with GRACE data and as input to modeling of glacio-isostatic uplift in this region. The study contributed to a better understanding of how glacier systems adjust to ice shelf disintegration.

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