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

Comparative responses of two dominant Antarctic phytoplankton taxa to interactions between ocean acidification, warming, irradiance, and iron availability

We investigated the responses of the ecologically dominant Antarctic phytoplankton species Phaeocystis antarctica (a prymnesiophyte) and Fragilariopsis cylindrus (a diatom) to a clustered matrix of three global change variables (CO2, mixed-layer depth, and temperature) under both iron (Fe)-replete and Fe-limited conditions based roughly on the Intergovernmental Panel on Climate Change (IPCC) A2 scenario: (1) Current conditions, 39 Pa (380 ppmv) CO2, 50 µmol photons/m**2/s light, and 2°C; (2) Year 2060, 61 Pa (600 ppmv) CO2, 100 µmol photons/m**2/s light, and 4°C; (3) Year 2100, 81 Pa (800 ppmv) CO2, 150 µmol photons/m**2/s light, and 6°C. The combined interactive effects of these global change variables and changing Fe availability on growth, primary production, and cell morphology are species specific. A competition experiment suggested that future conditions could lead to a shift away from P. antarctica and toward diatoms such as F. cylindrus. Along with decreases in diatom cell size and shifts from prymnesiophyte colonies to single cells under the future scenario, this could potentially lead to decreased carbon export to the deep ocean. Fe : C uptake ratios of both species increased under future conditions, suggesting phytoplankton of the Southern Ocean will increase their Fe requirements relative to carbon fixation. The interactive effects of Fe, light, CO2, and temperature on Antarctic phytoplankton need to be considered when predicting the future responses of biology and biogeochemistry in this region.

Spatial distribution patterns of ascidians (Ascidiacea: Tunicata) in combination with information on bathymetry, oceanography, chlorophyll-a, and sea-ice on the continental shelves off the northern Antarctic Peninsula during POLARSTERN cruise ANT-XXIX/3

Ascidians (Ascidiacea: Tunicata) are sessile suspension feeders that represent dominant epifaunal components of the Southern Ocean shelf benthos and play a significant role in the pelagic-benthic coupling. Here, we report the results of a first study on the relationship between the distribution patterns of eight common and/or abundant (putative) ascidian species, and environmental drivers in the waters off the northern Antarctic Peninsula. During RV Polarstern cruise XXIX/3 (PS81) in January-March 2013, we used seabed imaging surveys along 28 photographic transects of 2 km length each at water depths from 70 to 770 m in three regions (northwestern Weddell Sea, southern Bransfield Strait and southern Drake Passage), differing in their general environmental setting, primarily oceanographic characteristics and sea-ice dynamics, to comparatively analyze the spatial patterns in the abundance of the selected ascidians, reliably to be identified in the photographs, at three nested spatial scales. At a regional (100-km) scale, the ascidian assemblages of the Weddell Sea differed significantly from those of the other two regions, whereas at an intermediate 10-km scale no such differences were detected among habitat types (bank, upper slope, slope, deep/canyon) on the shelf and at the shelf break within each region. These spatial patterns were superimposed by a marked small-scale (10-m) patchiness of ascidian distribution within the 2-km-long transects. Among the environmental variables considered in our study, a combination of water-mass characteristics, sea-ice dynamics (approximated by 5-year averages in sea-ice cover in the region of or surrounding the photographic stations), as well as the seabed ruggedness, was identified as explaining best the distribution patterns of the ascidians.

Tab. 1: Zircon isotopic data and rock ages from the Antarctic Peninsula

Morphological and U-Pb isotope studies on sedimentary zircons reflect the orogenic evolution of their former host rocks. The orogenic history of detrital zircons from the Trinity Peninsula Formation (TPF) defines the former geological surrounding of the sedimentation basin of the TPF. Same few weil rounded, polycyclic zircons of Precambrian age and Cambrian overprint give hints for an old cratonic source rock. Because of their very low frequency compared with euhedral types, the contribution of an cratonic shield area to the bulk of the sedimentary debris is neglectable low. Euhedral zircons of granitoid origin and Carboniferous age indicate a derivation from an area of widespread Carboniferous intrusions. Except for southern South America and unsurveyed regions in the Antarctic Peninsula itself, no region could deliver zircons with a Carboniferous age record. The only acceptable explanation for the origin of these zircons is a position of the Antarctic Peninsula during the sedimentation of the TPF approximately southwest of southern South America.

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.

Table 2.- Bryozoa collected during the LAMPOS expeditions between Punta Arenas and the northernmost tip of the Antarctic Peninsula

The 78 bryozoan species collected by the German R/V "Polarstern" during the LAMPOS cruise in April 2002, encompassing the Scotia Arc archipelagos between Tierra del Fuego and the Antarctic Peninsula, were studied to discern the biogeographical patterns of the Magellan region of South America, the Scotia Arc archipelagos and the Antarctic. The resulting dendrogram shows three clusters: an isolated one with the three easternmost archipelagos and the other two linking some of the northern and southern Scotia Arc archipelagos with Tierra del Fuego. A more comprehensive analysis using all the species previously recorded from the Scotia Arc archipelagos and adjacent areas (214 spp.) produced a clearer zoogeographical pattern without isolated clusters of localities. The Antarctic Peninsula plus the Scotia Arc archipelagos form a large cluster distinct from the Magellan-Falkland Subantarctic area. A third analysis making use of 78 genera present in the study area plus Australia and New Zealand reinforces this pattern, showing two clusters: one uniting South America and the Australian-New Zealand realm and the other linking the Scotia Arc archipelagos with the Antarctic Peninsula. These results indicate that the Scotia Arc archipelagos represent merely a very narrow bridge connecting two different bryozoan faunas with only a few bryozoan species in common between the study areas.

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.

Thermal measurements on sub-Antarctic megaherbs and environmental characteristics on Campbell Island in 2010

We evaluated whether heating occurs in sub-Antarctic megaherbs, and the relation of heating to relevant environmental variables. We measured leaf and inflorescence temperature in six sub-Antarctic megaherb species on Campbell Island, latitude 52.3°S, New Zealand Biological Region. Using thermal imaging camera (Fluke TI20, http://www.fluke.com/fluke/caen/support/software/ti-update) and thermal probe (Fluke 51 II digital thermal probe), in combination with measurement of solar radiation, ambient air temperature, wind speed, wind chill and humidity.

Seawater carbonate chemistry and hatch rates of Antarctic krill

Marine ecosystems of the Southern Ocean are particularly vulnerable to ocean acidification. Antarctic krill (Euphausia superba; hereafter krill) is the key pelagic species of the region and its largest fishery resource. There is therefore concern about the combined effects of climate change, ocean acidification and an expanding fishery on krill and ultimately, their dependent predators-whales, seals and penguins. However, little is known about the sensitivity of krill to ocean acidification. Juvenile and adult krill are already exposed to variable seawater carbonate chemistry because they occupy a range of habitats and migrate both vertically and horizontally on a daily and seasonal basis. Moreover, krill eggs sink from the surface to hatch at 700-1,000 m, where the carbon dioxide partial pressure (pCO2) in sea water is already greater than it is in the atmosphere. Krill eggs sink passively and so cannot avoid these conditions. Here we describe the sensitivity of krill egg hatch rates to increased CO2, and present a circumpolar risk map of krill hatching success under projected pCO2 levels. We find that important krill habitats of the Weddell Sea and the Haakon VII Sea to the east are likely to become high-risk areas for krill recruitment within a century. Furthermore, unless CO2 emissions are mitigated, the Southern Ocean krill population could collapse by 2300 with dire consequences for the entire ecosystem.

Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo 1.0.5-beta)

Sub-ice shelf circulation and freezing/melting rates in ocean general circulation models depend critically on an accurate and consistent representation of cavity geometry. Existing global or pan-Antarctic data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional fields into a global data set. We use the S-2004 global 1-minute bathymetry as the backbone and add an improved version of the BEDMAP topography for an area that roughly coincides with the Antarctic continental shelf. Locations of the merging line have been carefully adjusted in order to get the best out of each data set. High-resolution gridded data for upper and lower ice surface topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier have been carefully merged into the ambient ice and ocean topographies. Multibeam survey data for bathymetry in the former Larsen B cavity and the southeastern Bellingshausen Sea have been obtained from the data centers of Alfred Wegener Institute (AWI), British Antarctic Survey (BAS) and Lamont-Doherty Earth Observatory (LDEO), gridded, and again carefully merged into the existing bathymetry map. The global 1-minute dataset (RTopo-1 Version 1.0.5) has been split into two netCDF files. The first contains digital maps for global bedrock topography, ice bottom topography, and surface elevation. The second contains the auxiliary maps for data sources and the surface type mask. A regional subset that covers all variables for the region south of 50 deg S is also available in netCDF format. Datasets for the locations of grounding and coast lines are provided in ASCII format.

Opal chemistry of costal Antarctic sediments

During the last 50 years, the Antarctic Peninsula has experienced rapid warming with associated retreat of 87% of marine and tidewater glacier fronts. Accelerated glacial retreat and iceberg calving may have a significant impact on the freshwater and nutrient supply to the phytoplankton communities of the highly productive coastal regions. However, commonly used biogenic carbonate proxies for nutrient and salinity conditions are not preserved in sediments from coastal Antarctica. Here we describe a method for the measurement of zinc to silicon ratios in diatom opal, (Zn/Si)opal, which is a potential archive in Antarctic marine sediments. A core top calibration from the West Antarctic Peninsula shows (Zn/Si)opal is a proxy for mixed layer salinity. We present down-core (Zn/Si)opal paleosalinity records from two rapidly accumulating sites taken from nearshore environments off the West Antarctic Peninsula which show an increase in meltwater input in recent decades. Our records show that the recent melting in this region is unprecedented for over 120 years.