Search

The Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2) is a natural progression of previous and ongoing model intercomparison exercises that have focused on the simulation of ice-sheet and ocean processes in Antarctica. The previous exercises motivate the move towards realistic configurations, as well as more diverse model parameters and resolutions. The main objective of MISOMIP2 is to investigate the performance of existing ocean and coupled ice-sheet–ocean models in a range of Antarctic environments through comparisons to observational data. We will assess the status of ice-sheet–ocean modelling as a community and identify common characteristics of models that are best able to capture observed features. As models are highly tuned based on present-day data, we will also compare their sensitivity to prescribed abrupt atmospheric perturbations leading to either very warm or slightly warmer ocean conditions compared to the present day. The approach of MISOMIP2 is to welcome contributions of models as they are, including global and regional configurations, but we request standardized variables and common grids for the outputs. We target the analysis at two specific regions, the Amundsen Sea and the Weddell Sea, since they describe two different ocean environments and have been relatively well observed compared to other areas of Antarctica. An observational “MIPkit” synthesizing existing ocean and ice-sheet observations for a common period is provided to evaluate ocean and ice-sheet models in these two regions.

The Pine Island and Thwaites glaciers are the two largest contributors to sea level rise from Antarctica. Here we examine the influence of basal friction and ice shelf basal melt in determining projected losses. We examine both Weertman and Coulomb friction laws with explicit weakening as the ice thins to flotation, which many friction laws include implicitly via the effective pressure. We find relatively small differences with the choice of friction law (Weertman or Coulomb) but find losses to be highly sensitive to the rate at which the basal traction is reduced as the area upstream of the grounding line thins. Consistent with earlier work on Pine Island Glacier, we find sea level contributions from both glaciers to vary linearly with the melt volume averaged over time and space, with little influence from the spatial or temporal distribution of melt. Based on recent estimates of melt from other studies, our simulations suggest that the combined melt-driven and sea level rise contribution from both glaciers may not exceed 10 cm by 2200, although the uncertainty in model parameters allows for larger increases. We do not include other factors, such as ice shelf breakup, that might increase loss, or factors such as increased accumulation and isostatic uplift that may mitigate loss.

Glaciers in the Amundsen Sea Embayment of West Antarctica are rapidly retreating and contributing to sea level rise. Ice loss is occurring primarily via exposure to warm ocean water, which varies in response to local wind variability. There is evidence that retreat was initiated in the mid-20th century, but the perturbation that may have triggered retreat remains unknown. A leading hypothesis is that large pressure and wind anomalies in the 1940s drove exceptionally strong oceanic ice-shelf melting. However, the characteristics, drivers, and rarity of the atmospheric event remain poorly constrained. We investigate the 1940s atmospheric event using paleoclimate reconstructions and climate model simulations. The reconstructions show that large westerly wind anomalies occurred from ∼1938–1942, a combined response to the very large El Niño event from 1940–1942 and other variability beginning years earlier. Climate model simulations provide evidence that events of similar magnitude and duration may occur tens to hundreds of times per 10 kyr of internal climate variability (∼0.2 to 2.5 occurrences per century). Our results suggest that the 1940s westerly event is unlikely to have been exceptional enough to be the sole explanation for the initiation of Amundsen Sea glacier retreat. Additional factors are likely needed to explain the onset of retreat in West Antarctica, such as naturally arising variability in ocean conditions prior to the 1940s or anthropogenically driven trends since the 1940s.

Ocean-driven ice loss from the West Antarctic Ice Sheet is a significant contributor to sea-level rise. Recent ocean variability in the Amundsen Sea is controlled by near-surface winds. We combine palaeoclimate reconstructions and climate model simulations to understand past and future influences on Amundsen Sea winds from anthropogenic forcing and internal climate variability. The reconstructions show strong historical wind trends. External forcing from greenhouse gases and stratospheric ozone depletion drove zonally uniform westerly wind trends centred over the deep Southern Ocean. Internally generated trends resemble a South Pacific Rossby wave train and were highly influential over the Amundsen Sea continental shelf. There was strong interannual and interdecadal variability over the Amundsen Sea, with periods of anticyclonic wind anomalies in the 1940s and 1990s, when rapid ice-sheet loss was initiated. Similar anticyclonic anomalies probably occurred prior to the 20th century but without causing the present ice loss. This suggests that ice loss may have been triggered naturally in the 1940s but failed to recover subsequently due to the increasing importance of anthropogenic forcing from greenhouse gases (since the 1960s) and ozone depletion (since the 1980s). Future projections also feature strong wind trends. Emissions mitigation influences wind trends over the deep Southern Ocean but has less influence on winds over the Amundsen Sea shelf, where internal variability creates a large and irreducible uncertainty. This suggests that strong emissions mitigation is needed to minimise ice loss this century but that the uncontrollable future influence of internal climate variability could be equally important.

Abstract At the Trans‐Atlantic Geotraverse hydrothermal field, metalliferous sediments cover extinct hydrothermal mounds and the surrounding seafloor. Here, we report the morphological, mineralogical and geochemical processes that deposit these sediments, remobilize their metals, and affect their preservation. We found that the initial sediment metal tenor is controlled by physical transport of hydrothermal material from its source, followed by diagenetic redistribution and potentially diffuse fluid flow after high‐temperature hydrothermal activity has ceased. We distinguished three different environments: (a) proximal metalliferous sediments on top of extinct mounds are mainly derived from oxidative weathering of primary sulfide structures and are predominantly composed of Fe oxyhydroxides with low contents of Cu, Co, and Zn; metal enrichments in specific layers are likely related to upward flow of low‐temperature hydrothermal fluids; (b) medial distant metalliferous sediments found at the base of the mounds, deposited by mass transport, contain cm‐thick layers of unsorted sulfide sands with high base metal contents (e.g., up to 28% Cu); these buried sulfides continue to undergo dissolution, resulting in metal release into porewaters; (c) distal metalliferous sediments, found in depositional basins a few hundreds of meters from the extinct mounds, include fining‐upwards sequences of thin sulfide sand layers with Fe oxyhydroxides and were deposited by recurrent turbiditic flows. Dissolved metals (e.g., Cu2+ and Mn2+) diffuse upwards under reducing conditions and precipitate within the sediment. Hence, when using hydrothermal sediments to construct reliable geochronological records of hydrothermal activity, distance from source, local seafloor morphology, mass‐transport and depositional, and diagenetic modification should all be considered.

Ocean-induced basal melting is responsible for much of the Amundsen Sea Embayment ice loss in recent decades, but the total magnitude and spatiotemporal evolution of this melt is poorly constrained. To address this problem, we generated a record of high-resolution digital elevation models (DEMs) for Pine Island Glacier (PIG) using commercial sub-meter satellite stereo imagery and integrated additional 2002–2015 DEM and altimetry data. We implemented a Lagrangian elevation change (Dh∕Dt) framework to estimate ice shelf basal melt rates at 32–256 m resolution. We describe this methodology and consider basal melt rates and elevation change over the PIG ice shelf and lower catchment from 2008 to 2015. We document the evolution of Eulerian elevation change (dh∕dt) and upstream propagation of thinning signals following the end of rapid grounding line retreat around 2010. Mean full-shelf basal melt rates for the 2008–2015 period were ∼82–93 Gt yr−1, with ∼200–250 m yr−1 basal melt rates within large channels near the grounding line, ∼10–30 m yr−1 over the main shelf, and ∼0–10 m yr−1 over the North shelf and South shelf, with the notable exception of a small area with rates of ∼50–100 m yr−1 near the grounding line of a fast-flowing tributary on the South shelf. The observed basal melt rates show excellent agreement with, and provide context for, in situ basal melt-rate observations. We also document the relative melt rates for kilometer-scale basal channels and keels at different locations on the ice shelf and consider implications for ocean circulation and heat content. These methods and results offer new indirect observations of ice–ocean interaction and constraints on the processes driving sub-shelf melting beneath vulnerable ice shelves in West Antarctica.

By thinning and accelerating, West Antarctic ice streams are contributing about 10% of the observed global sea level rise. Much of this ice loss is from Pine Island Glacier, which has thinned since at least 1992, driven by changes in ocean heat transport beneath its ice shelf and retreat of the grounding line. Details of the processes driving this change, however, remain largely elusive, hampering our ability to predict the future behaviour of this and similar systems. Here, a Lagrangian methodology is developed to measure oceanic melting of such rapidly advecting ice. High-resolution satellite and airborne observations of ice surface velocity and elevation are used to quantify patterns of basal melt under the Pine Island Glacier ice shelf and the associated adjustments to ice flow. At the broad scale, melt rates of up to 100 m yr−1 occur near the grounding line, reducing to 30 m yr−1 just 20 km downstream. Between 2008 and 2011, basal melting was largely compensated by ice advection, allowing us to estimate an average loss of ice to the ocean of 87 km3 yr−1, in close agreement with 2009 oceanographically constrained estimates. At smaller scales, a network of basal channels typically 500 m to 3 km wide is sculpted by concentrated melt, with kilometre-scale anomalies reaching 50% of the broad-scale basal melt. Basal melting enlarges the channels close to the grounding line, but farther downstream melting tends to diminish them. Kilometre-scale variations in melt are a key component of the complex ice–ocean interaction beneath the ice shelf, implying that greater understanding of their effect, or very high resolution models, are required to predict the sea-level contribution of the region.

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 topography data sets have turned out to contain various inconsistencies and inaccuracies. The goal of this work is to compile independent regional surveys and maps into a global data set. We use the S-2004 global 1-min bathymetry as the backbone and add an improved version of the BEDMAP topography (ALBMAP bedrock topography) for an area that roughly coincides with the Antarctic continental shelf. The position of the merging line is individually chosen in different sectors in order to capture the best of both data sets. High-resolution gridded data for ice shelf topography and cavity geometry of the Amery, Fimbul, Filchner-Ronne, Larsen C and George VI Ice Shelves, and for Pine Island Glacier are 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 blended into the existing bathymetry map. The resulting global 1-min Refined Topography data set (RTopo-1) contains self-consistent maps for upper and lower ice surface heights, bedrock topography, and surface type (open ocean, grounded ice, floating ice, bare land surface). The data set is available in NetCDF format from the PANGAEA database at doi:10.1594/pangaea.741917.