Sulfur Isotopic Measurements from a West Antarctic Ice Core: Implications for Sulfate Source and Transport

Measurements of delta(34)S covering the years 1935-76 and including the 1963 Agung (Indonesia) eruption were made on a West Antarctic firn core, RIDSA (78.73 degrees S, 116.33 degrees W; 1740m a.s.l.), and results are used to unravel potential source functions in the sulfur cycle over West Antarctica. The delta(34)S values Of SO42- range from 3.1 parts per thousand to 9.9 parts per thousand. These values are lower than those reported for central Antarctica, from near South Pole station, of 9.3-18.1 parts per thousand (Patris and others, 2000). While the Agung period is isotopically distinct at South Pole, it is not in the RIDSA dataset, suggesting differences in the source associations for the sulfur cycle between these two regions. Given the relatively large input of marine aerosols at RIDSA (determined from Na+ data and the seasonal SO42- cycle), there is likely a large marine biogenic SO42- influence. The delta(34)S values indicate, however, that this marine biogenic SO42-, with a well-established delta(34)S of 18 parts per thousand, is mixing with SO42- that has extremely negative delta(34)S values to produce the measured isotope values in the RIDSA core. We suggest that the transport and deposition of stratospheric SO42- in West Antarctica, combined with local volcanic input, accounts for the observed variance in delta(34)S values.

A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum

A robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.

(Table 1) Basin area, average late-summer speedup, number of sinks and runoff of West Greenland Ice Sheet catchments

We use interferometric synthetic aperture radar observations recorded in a land-terminating sector of western Greenland to characterise the ice sheet surface hydrology and to quantify spatial variations in the seasonality of ice sheet flow. Our data reveal a non-uniform pattern of late-summer ice speedup that, in places, extends over 100 km inland. We show that the degree of late-summer speedup is positively correlated with modelled runoff within the 10 glacier catchments of our survey, and that the pattern of late-summer speedup follows that of water routed at the ice sheet surface. In late-summer, ice within the largest catchment flows on average 48% faster than during winter, whereas changes in smaller catchments are less pronounced. Our observations show that the routing of seasonal runoff at the ice sheet surface plays an important role in shaping the magnitude and extent of seasonal ice sheet speedup.

Simulated temperature changes as a result of Antarctic ice sheet growth and CO2 reduction

During the Middle Miocene climate transition about 14 million years ago, the Antarctic ice sheet expanded to near-modern volume. Surprisingly, this ice sheet growth was accompanied by a warming in the surface waters of the Southern Ocean, whereas a slight deep-water temperature increase was delayed by more than 200 thousand years. Here we use a coupled atmosphere-ocean model to assess the relative effects of changes in atmospheric CO2 concentration and ice sheet growth on regional and global temperatures. In the simulations, changes in the wind field associated with the growth of the ice sheet induce changes in ocean circulation, deep-water formation and sea-ice cover that result in sea surface warming and deep-water cooling in large swaths of the Atlantic and Indian ocean sectors of the Southern Ocean. We interpret these changes as the dominant ocean surface response to a 100-thousand-year phase of massive ice growth in Antarctica. A rise in global annual mean temperatures is also seen in response to increased Antarctic ice surface elevation. In contrast, the longer-term surface and deep-water temperature trends are dominated by changes in atmospheric CO2 concentration. We therefore conclude that the climatic and oceanographic impacts of the Miocene expansion of the Antarctic ice sheet are governed by a complex interplay between wind field, ocean circulation and the sea-ice system.

Planimetric ice-flow convergence and flow-orthonormal strain rate across the Antarctic Ice Sheet, with links to model result files

Fast-flowing ice streams discharge most of the ice from the interior of the Antarctic Ice Sheet coastward. Understanding how their tributary organisation is governed and evolves is essential for developing reliable models of the ice sheet's response to climate change. Despite much research on ice-stream mechanics, this problem is unsolved, because the complexity of flow within and across the tributary networks has hardly been interrogated. Here I present the first map of planimetric flow convergence across the ice sheet, calculated from satellite measurements of ice surface velocity, and use it to explore this complexity. The convergence map of Antarctica elucidates how ice-stream tributaries draw ice from the interior. It also reveals curvilinear zones of convergence along lateral shear margins of streaming, and abundant convergence ripples associated with nonlinear ice rheology and changes in bed topography and friction. Flow convergence on ice-stream tributaries and their feeding zones is markedly uneven, and interspersed with divergence at distances of the order of kilometres. For individual drainage basins as well as the ice sheet as a whole, the range of convergence and divergence decreases systematically with flow speed, implying that fast flow cannot converge or diverge as much as slow flow. I therefore deduce that flow in ice-stream networks is subject to mechanical regulation that limits flow-orthonormal strain rates. These properties and the gridded data of convergence and flow-orthonormal strain rate in this archive provide targets for ice- sheet simulations and motivate more research into the origin and dynamics of tributarization.

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.

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.

(Table 1) Accumulation and ice discharge of West Antarctic glaciers draining into the Amundsen Sea between 1974 and 2007

The Advanced Land Observation System (ALOS) Phased-Array Synthetic-Aperture Radar (PALSAR) is an L-band frequency (1.27 GHz) radar capable of continental-scale interferometric observations of ice sheet motion. Here, we show that PALSAR data yield excellent measurements of ice motion compared to C-band (5.6 GHz) radar data because of greater temporal coherence over snow and firn. We compare PALSAR velocities from year 2006 in Pine Island Bay, West Antarctica with those spanning years 1974 to 2007. Between 1996 and 2007, Pine Island Glacier sped up 42% and ungrounded over most of its ice plain. Smith Glacier accelerated 83% and ungrounded as well. Their largest speed up are recorded in 2007. Thwaites Glacier is not accelerating but widening with time and its eastern ice shelf doubled its speed. Total ice discharge from these glaciers increased 30% in 12 yr and the net mass loss increased 170% from 39 ± 15 Gt/yr to 105 ± 27 Gt/yr. Longer-term velocity changes suggest only a moderate loss in the 1970s. As the glaciers unground into the deeper, smoother beds inland, the mass loss from this region will grow considerably larger in years to come.

A 12,000 Year Record of Explosive Volcanism in the Siple Dome Ice Core, West Antarctica

Air mass trajectories in the Southern Hemisphere provide a mechanism for transport to and deposition of volcanic products on the Antarctic ice sheet from local volcanoes and from tropical and subtropical volcanic centers. This study extends the detailed record of Antarctic, South American, and equatorial volcanism over the last 12,000 years using continuous glaciochemical series developed from the Siple Dome A (SDMA) ice core, West Antarctica. The largest volcanic sulfate spike ( 280 mu g/L) occurs at 5881 B. C. E. Other large signals with unknown sources are observed around 325 B. C. E. ( 270 mu g/L) and 2818 B. C. E. ( 191 mu g/L). Ages of several large equatorial or Southern Hemisphere volcanic eruptions are synchronous with many sulfate peaks detected in the SDMA volcanic ice chemistry record. The microprobe "fingerprinting'' of glass shards in the SDMA core points to the following Antarctic volcanic centers as sources of tephra found in the SDMA core: Balenny Island, Pleiades, Mount Berlin, Mount Takahe, and Mount Melbourne as well as Mount Hudson and possibly Mount Burney volcanoes of South America. Identified volcanic sources provide an insight into the poorly resolved transport history of volcanic products from source volcanoes to the West Antarctic ice sheet.

High-Resolution Ice Cores from US ITASE (West Antarctica): Development and Validation of Chronologies and Determination of Precision and Accuracy

Shallow ice cores were obtained from widely distributed sites across the West Antarctic ice sheet, as part of the United States portion of the International Trans-Antarctic Scientific Expedition (US ITASE) program. The US ITASE cores have been dated by annual-layer counting, primarily through the identification of summer peaks in non-sea-salt sulfate (nssSO(4)(2-)) concentration. Absolute dating accuracy of better than 2 years and relative dating accuracy better than 1 year is demonstrated by the identification of multiple volcanic marker horizons in each of the cores, Tambora, Indonesia (1815), being the most prominent. Independent validation is provided by the tracing of isochronal layers from site to site using high-frequency ice-penetrating radar observations, and by the timing of mid-winter warming events in stable-isotope ratios, which demonstrate significantly better than 1 year accuracy in the last 20 years. Dating precision to 1 month is demonstrated by the occurrence of summer nitrate peaks and stable-isotope ratios in phase with nssSO(4)(2-), and winter-time sea-salt peaks out of phase, with phase variation of < 1 month. Dating precision and accuracy are uniform with depth, for at least the last 100 years.

Corrected ICESat altimetry data, surface mass balance, and firn elevation change on Antarctic ice shelves

Accurate prediction of global sea-level rise requires that we understand the cause of recent, widespread and intensifying glacier acceleration along Antarctic ice-sheet coastal margins. Floating ice shelves buttress the flow of grounded tributary glaciers and their thickness and extent are particularly susceptible to changes in both climate and ocean forcing. Recent ice-shelf collapse led to retreat and acceleration of several glaciers on the Antarctic Peninsula. However, the extent and magnitude of ice-shelf thickness change, its causes and its link to glacier flow rate are so poorly understood that its influence on the future of the ice sheets cannot yet be predicted. Here we use satellite laser altimetry and modelling of the surface firn layer to reveal for the first time the circum-Antarctic pattern of ice-shelf thinning through increased basal melt. We deduce that this increased melt is the primary driver of Antarctic ice-sheet loss, through a reduction in buttressing of the adjacent ice sheet that has led to accelerated glacier flow. The highest thinning rates (~7 m/a) occur where warm water at depth can access thick ice shelves via submarine troughs crossing the continental shelf. Wind forcing could explain the dominant patterns of both basal melting and the surface melting and collapse of Antarctic ice shelves, through ocean upwelling in the Amundsen and Bellingshausen Seas and atmospheric warming on the Antarctic Peninsula. This implies that climate forcing through changing winds influences Antarctic Ice Sheet mass balance, and hence global sea-level, on annual to decadal timescales.

A 200 Year Sulfate Record from Sixteen Antarctic Ice Cores and Associations With Southern Ocean Sea-Ice Extent

Chemistry data from 16, 50-115 m deep, sub-annually dated ice cores are used to investigate spatial and temporal concentration variability of sea-salt (ss) SO42- and excess (xs) SO42- over West Antarctica and the South Pole for the last 200 years. Low-elevation ice-core sites in western West Antarctica contain higher concentrations Of SO42- as a result of cyclogenesis over the Ross Ice Shelf and proximity to the Ross Sea Polynya. Linear correlation analysis of 15 West Antarctic ice-core SO42- time series demonstrates that at several sites concentrations Of ssSO(4)(2-) are higher when sea-ice (SIE) extent is greater, and the inverse for XSS04. Concentrations Of XSS04 from the South Pole site (East Antarctica) are associated with SIE from the Weddell region, and West Antarctic XSSO42- concentrations are associated with SIE from the Bellingshausen-Amundsen-Ross region. The only notable rise of the last 200 years in xsSO(4)(2-), around 1940, is not related to SIE fluctuations and is most likely a result of increased xsSO(4)(2-) production in the mid-low latitudes and/or an increase in transport efficiency from the mid-low latitudes to central West Antarctica. These high-resolution records show that the source types and source areas Of ssSO(4)(2-) and xsSO(4)(2-) delivered to eastern and western West Antarctica and the South Pole differ from site to site but can best be resolved using records from spatial ice-core arrays such as the International Trans-Antarctic Scientific Expedition (ITASE).

ICESat Measurements Reveal Complex Pattern of Elevation Changes on Siple Coast Ice Streams, Antarctica

We compare ICESat data (2003-2004) to airborne laser altimetry data (1997-98 and 1999-2000) to monitor surface changes over portions of Van der Veen (VdVIS), Whillans (WIS) and Kamb ice streams (KIS) in the Ross Embayment of the West Antarctic Ice Sheet. The spatial pattern of detected surface changes is generally consistent with earlier observations. However, important changes have occurred during the past decade. For example, areas on the VdVIS and WIS, where large thinning was detected by the airborne surveys, are now closer to being in balance. The upper trunk of KIS continues to build up with thickening rates reaching 0.4 m/year. Our results provide new evidence that the overall mass balance of the region is becoming more positive, but a significant spatial variability exists. They also demonstrate the potential of ICESat data for detecting spatial patterns of surface elevation change in Antarctica.

On the Disintegration of Ice Shelves: The Role of Fracture

Crevasses can be ignored in studying the dynamics of most glaciers because they are only about 20 m deep, a small fraction of ice thickness. In ice shelves, however, s urface crevasses 20 m deep often reach sealevel and bottom crevasses can move upward to sea-level (Clough, 1974; Weertman, 1980). The ice shelf is fractured completely through if surface and basal crevasses meet (Barrett, 1975; Hughes, 1979). This is especially likely if surface melt water fills surface crevasses (Weertman, 1973; Pfeffer, 1982; Fastook and Schmidt, 1982). Fracture may therefore play an important role i n the disintegration of ice shelves. Two fracture criteria which can be evaluated experimentally and applied to ice shelves, are presented. Fracture is then examined for the general strain field of an ice shelf and for local strain fields caused by shear rupture alongside ice streams entering the ice shelf, fatigue rupture along ice shelf grounding lines, and buckling up-stream from ice rises. The effect of these fracture patterns on the stability of Antarctic ice shelves and the West Antarctic ice sheet is then discussed.