Search

The Amundsen Sea sector has some of the fastest-thinning ice shelves in Antarctica, caused by high, ocean-driven basal melt rates, which can lead to increased ice stream flow, causing increased sea level rise (SLR) contributions. In this study, we present the results of a new synchronously coupled ice-sheet/ocean model of the Amundsen Sea sector. We use the WAVI ice sheet model to solve for ice velocities and the MITgcm to solve for ice thickness and three-dimensional ocean properties, allowing for full mass conservation in the coupled ice/ocean system. The coupled model is initialised in the present day and run forward under idealised warm and cold ocean conditions. We find that Thwaites Glacier dominates the future SLR from the Amundsen Sea sector, with a SLR that is approximately quadratic in time. The future evolution of Thwaites Glacier depends on the life-span of small pinning points that form during the retreat. The rate of melting around these pinning points provides the link between future ocean conditions and the SLR from this sector and will be difficult to capture without a coupled ice/ocean model. Grounding-line retreat leads to a progressively larger Thwaites ice-shelf cavity, leading to a positive trend in total melting, resulting from the increased ice basal surface area. Despite these important sensitivities, Thwaites Glacier retreats even in a scenario with zero ocean-driven melting. This demonstrates that a tipping point may have been passed and some SLR from this sector is now committed.

The retreating Pine Island Glacier (PIG), West Antarctica, presently contributes ~5–10% of global sea-level rise. PIG’s retreat rate has increased in recent decades with associated thinning migrating upstream into tributaries feeding the main glacier trunk. To project future change requires modelling that includes robust parameterisation of basal traction, the resistance to ice flow at the bed. However, most ice-sheet models estimate basal traction from satellite-derived surface velocity, without a priori knowledge of the key processes from which it is derived, namely friction at the ice-bed interface and form drag, and the resistance to ice flow that arises as ice deforms to negotiate bed topography. Here, we present high-resolution maps, acquired using ice-penetrating radar, of the bed topography across parts of PIG. Contrary to lower-resolution data currently used for ice-sheet models, these data show a contrasting topography across the ice-bed interface. We show that these diverse subglacial landscapes have an impact on ice flow, and present a challenge for modelling ice-sheet evolution and projecting global sea-level rise from ice-sheet loss., Projecting the future retreat and thus global sea level contributions of Antarctica’s Pine Island Glacier is hampered by a poor grasp of what controls flow at the ice base. Here, via high-resolution ice-radar imaging, the authors show diverse landscapes beneath the glacier fundamentally influence ice flow.

West Antarctic Ice Sheet loss is a significant contributor to sea level rise. While the ice loss is thought to be triggered by fluctuations in oceanic heat at the ice shelf bases, ice sheet response to ocean variability remains poorly understood. Using a synchronously coupled ice-ocean model permitting grounding line migration, this study evaluates the response of an ice sheet to periodic variations in ocean forcing. Resulting oscillations in grounded ice volume amplitude is shown to grow as a nonlinear function of ocean forcing period. This implies that slower oscillations in climatic forcing are disproportionately important to ice sheets. The ice shelf residence time offers a critical time scale, above which the ice response amplitude is a linear function of ocean forcing period and below which it is quadratic. These results highlight the sensitivity of West Antarctic ice streams to perturbations in heat fluxes occurring at decadal time scales.

Accurate dynamical models of the Antarctic ice sheet with carefully specified initial conditions and well-calibrated rheological parameters are needed to forecast global sea level. By adapting an inverse method previously used in electric impedance tomography, we infer present-day flow speeds within the ice sheet. This inversion uses satellite observations of surface velocity, snow accumulation rate, and rate of change of surface elevation to estimate the basal drag coefficient and an ice stiffness parameter that influences viscosity. We represent interior ice motion using a vertically integrated approximation to incompressible Stokes flow. This model represents vertical shearing within the ice and membrane stresses caused by horizontal stretching and shearing. Combining observations and model, we recover marked geographical variations in the basal drag coefficient. Relative changes in basal shear stress are smaller. No simple sliding law adequately represents basal shear stress as a function of sliding speed. Low basal shear stress predominates in central East Antarctica, where thick insulating ice allows liquid water at the base to lubricate sliding. Higher shear stress occurs in coastal East Antarctica, where a frozen bed is more likely. Examining Thwaites glacier in more detail shows that the slowest sliding often coincides with elevated basal topography. Differences between our results and a similar adjoint-based inversion suggest that inversion or regularization methods can influence recovered parameters for slow sliding and finer scales; on broader scales we recover a similar pattern of low basal drag underneath major ice streams and extensive regions in East Antarctica that move by basal sliding.

[1] The air content of glacial firn determines the effect and attribution of observed changes in ice surface elevation, but is currently measurable only using labor‐intensive ground‐based techniques. Here a novel method is presented for using radar sounding measurements to decompose the total thickness of floating ice shelves into thicknesses of solid ice and firn air (or firn water). The method is applied to a 1997/98 airborne survey of Larsen Ice Shelf, revealing large spatial gradients in air content that are consistent with existing measurements and local meteorology. The gradients appear to be governed by meltwater‐induced firn densification. We find sufficient air in Larsen C Ice Shelf for increased densification to account for its previously observed surface lowering, and the rate of lowering superficially agrees with published trends in melting. This does not preclude a contribution to the lowering from oceanic basal melting, but a modern repeat of the survey could conclusively distinguish atmosphere‐led from ocean‐led change. The technique also holds promise for the calibration of firn‐density models, derivation of ice thickness from surface elevation measurements, and calculation of the sea‐level contribution of changes in grounded‐ice discharge.

We date a recent ice-flow reorganization of an ice divide in the Weddell Sea Sector, West Antarctica, using a novel combination of inverse methods and ice-penetrating radars. We invert for two-dimensional ice flow within an ice divide from data collected with a phase-sensitive ice-penetrating radar while accounting for the effect of firn on radar propagation and ice flow. By comparing isochronal layers simulated using radar-derived flow velocities with internal layers observed with an impulse radar, we show that the divide's internal structure is not in a steady state but underwent a disturbance, potentially implying a regional ice-flow reorganization, 2.5 (1.8–2.9) kyr B.P. Our data are consistent with slow ice flow in this location before the reorganization and the ice divide subsequently remaining stationary. These findings increase our knowledge of the glacial history of a region that lacks dated constraints on late-Holocene ice-sheet retreat and provides a key target for models that reconstruct and predict ice-sheet behavior.