Knowledge of frontal ablation from marine-terminating glaciers (i.e., mass lost at the calving face) is critical for constraining glacier mass balance, improving projections of mass change, and identifying the processes that govern frontal mass loss. Here, we discuss the challenges involved in computing frontal ablation and the unique issues pertaining to both glaciers and ice sheets. Frontal ablation estimates require numerous datasets, including glacier terminus area change, thickness, surface velocity, density, and climatic mass balance. Observations and models of these variables have improved over the past decade, but significant gaps and regional discrepancies remain, and better quantification of temporal variability in frontal ablation is needed. Despite major advances in satellite-derived large-scale datasets, large uncertainties remain with respect to ice thickness, depth-averaged velocities, and the bulk density of glacier ice close to calving termini or grounding lines. We suggest ways in which we can move toward globally complete frontal ablation estimates, highlighting areas where we need improved datasets and increased collaboration.

Water pressure beneath glaciers influences ice velocity. Subglacial hydrology models are helpful for gaining insight into basal conditions, but models depend on unconstrained parameters, and a current challenge is reproducing elevated water pressures in winter. We eliminate terms related to englacial storage, opening by sliding, and melt due to changes in the pressure-melting-point temperature, to create a minimalist version of the Subglacial Hydrology And Kinetic, Transient Interactions (SHAKTI) model, and apply this model to Helheim Glacier in east Greenland to explore the winter base state of the subglacial drainage system. Our results suggest that meltwater produced at the bed alone supports active winter drainage with large areas of elevated water pressure and preferential drainage pathways, using a continuum approach that allows for transitions between flow regimes. Transmissivity varies spatially over several orders of magnitude from 10−4 to 103 m2s−1, with regions of weak transmissivity representing poorly connected regions of the system. Bed topography controls the location of primary drainage pathways, and high basal melt rates occur along the steep valley walls. Frictional heat from sliding is a dominant source of basal melt; different approaches for calculating basal shear stress produce significantly different basal melt rates and subglacial discharge.

Frontal ablation, the combination of submarine melting and iceberg calving, changes the geometry of a glacier's terminus, influencing glacier dynamics, the fate of upwelling plumes and the distribution of submarine meltwater input into the ocean. Directly observing frontal ablation and terminus morphology below the waterline is difficult, however, limiting our understanding of these coupled ice–ocean processes. To investigate the evolution of a tidewater glacier's submarine terminus, we combine 3-D multibeam point clouds of the subsurface ice face at LeConte Glacier, Alaska, with concurrent observations of environmental conditions during three field campaigns between 2016 and 2018. We observe terminus morphology that was predominately overcut (52% in August 2016, 63% in May 2017 and 74% in September 2018), accompanied by high multibeam sonar-derived melt rates (4.84 m d−1 in 2016, 1.13 m d−1 in 2017 and 1.85 m d−1 in 2018). We find that periods of high subglacial discharge lead to localized undercut discharge outlets, but adjacent to these outlets the terminus maintains significantly overcut geometry, with an ice ramp that protrudes 75 m into the fjord in 2017 and 125 m in 2018. Our data challenge the assumption that tidewater glacier termini are largely undercut during periods of high submarine melting.

Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica.

We perform Bayesian inference of the parameters of a time-dependent model of ice flow and calving at Helheim Glacier, East Greenland. We find that, while a time-independent calving parameterization can recover the mean observed terminus position, such a model is unable to recover sub-annual variability, even when forced with seasonally varying climate. To address this, we develop a simple stochastic model relating surface runoff rates and calving threshold. Again inferring model parameters from observations, we find that this parameterization is able to reproduce observations with respect to both mean position and characteristic temporal variability. This result demonstrates the importance of considering potential sub-annual controls on calving rates in numerical models, which may include variable undercutting rates or surface runoff-dependent surface crevasse propagation.

The Drygalski Ice Tongue in East Antarctica stretches 90 km into the Ross Sea and influences the local ocean circulation, and persistence of the Terra Nova Bay Polynya. We examine the controls on the size of this floating ice body by comparing the propagation of six large fractures on the ice tongue's northern side using 21 years of Landsat imagery with hydrostatic ice thickness maps and strain rate calculations. We also apply a subglacial hydrology model to estimate the location and discharge from subglacial channels over the grounding line and compare these with basal channels identified along the ice tongue using remote sensing and airborne radar data. Our results suggest that large fractures are inhibited from full-width propagation by thicker ice between basal channels. We hypothesize that only once the ice tongue thins towards the terminus, can fractures propagate and cause large calving events. This suggests an important relationship between the melting of floating ice from subglacial and ocean sources and the expansion of fractures that lead to ice tongue calving.

Antarctica's ice shelves modulate the grounded ice flow, and weakening of ice shelves due to climate forcing will decrease their ‘buttressing’ effect, causing a response in the grounded ice. While the processes governing ice-shelf weakening are complex, uncertainties in the response of the grounded ice sheet are also difficult to assess. The Antarctic BUttressing Model Intercomparison Project (ABUMIP) compares ice-sheet model responses to decrease in buttressing by investigating the ‘end-member’ scenario of total and sustained loss of ice shelves. Although unrealistic, this scenario enables gauging the sensitivity of an ensemble of 15 ice-sheet models to a total loss of buttressing, hence exhibiting the full potential of marine ice-sheet instability. All models predict that this scenario leads to multi-metre (1–12 m) sea-level rise over 500 years from present day. West Antarctic ice sheet collapse alone leads to a 1.91–5.08 m sea-level rise due to the marine ice-sheet instability. Mass loss rates are a strong function of the sliding/friction law, with plastic laws cause a further destabilization of the Aurora and Wilkes Subglacial Basins, East Antarctica. Improvements to marine ice-sheet models have greatly reduced variability between modelled ice-sheet responses to extreme ice-shelf loss, e.g. compared to the SeaRISE assessments.

Ice flow numerical models are essential for predicting the evolution of ice sheets in a warming climate. Recent research emphasizes the need for higher-order and even full-Stokes flow models, instead of the traditional shallow-ice approximation, whose assumptions are not valid in certain critical areas. These higher-order models are, however, computationally intensive and difficult to use at the continental scale. Here we present a new technique, the Tiling method, to couple ice flow models of varying orders of complexity. The goal of the method is to limit the spatial extent of where higherorder models are applied to reduce the computational cost, while maintaining the model precision. We apply this method on synthetic geometries to demonstrate its practical use. We first use a geometry for which all models yield the same results to check the consistency of the method. Then we apply our method to a geometry for which a full-Stokes model is required in the vicinity of the ice front. Our results show that the hybrid models present significant improvements over mono-model approaches and reduce computational times.

Ten ice-sheet models are used to study sensitivity of the Greenland and Antarctic ice sheets to prescribed changes of surface mass balance, sub-ice-shelf melting and basal sliding. Results exhibit a large range in projected contributions to sea-level change. In most cases, the ice volume above flotation lost is linearly dependent on the strength of the forcing. Combinations of forcings can be closely approximated by linearly summing the contributions from single forcing experiments, suggesting that nonlinear feedbacks are modest. Our models indicate that Greenland is more sensitive than Antarctica to likely atmospheric changes in temperature and precipitation, while Antarctica is more sensitive to increased ice-shelf basal melting. An experiment approximating the Intergovernmental Panel on Climate Change’s RCP8.5 scenario produces additional first-century contributions to sea level of 22.3 and 8.1 cm from Greenland and Antarctica, respectively, with a range among models of 62 and 14 cm, respectively. By 200 years, projections increase to 53.2 and 26.7 cm, respectively, with ranges of 79 and 43 cm. Linear interpolation of the sensitivity results closely approximates these projections, revealing the relative contributions of the individual forcings on the combined volume change and suggesting that total ice-sheet response to complicated forcings over 200 years can be linearized.

Tidewater glaciers in Greenland experienced widespread retreat during the last century. Information on their behaviour prior to this is often poorly constrained due to lack of observations, while determining the drivers prior to instrumental records is also problematic. Here we present a record of the dynamics of Kangiata Nunaata Sermia (KNS), southwest Greenland, from its Little Ice Age maximum (LIAmax) to 1859 – the period before continuous air temperature observations began at Nuuk in 1866. Using glacial geomorphology, historical accounts, photographs and GIS analyses, we provide evidence KNS was at its LIAmax by 1761, had retreated by ~5 km by 1808 and a further 7 km by 1859. This predates retreat at Jakobshavn Isbræ by 43–113 years, demonstrating the asynchroneity of tidewater glacier terminus response following the LIA. We use a one-dimensional flowband model to determine the relative sensitivity of KNS to atmospheric and oceanic climate forcing. Results demonstrate that terminus forcing rather than surface mass balance drove the retreat. Modelled glacier sensitivity to submarine melt rates is also insufficient to explain the retreat observed. However, moderate increases in crevasse water depth, driving an increase in calving, are capable of causing terminus retreat of the observed magnitude and timing.

Predictions of marine ice-sheet behaviour require models able to simulate grounding-line migration. We present results of an intercomparison experiment for plan-view marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no buttressing effects from lateral drag). Perturbation experiments specifying spatial variation in basal sliding parameters permitted the evolution of curved grounding lines, generating buttressing effects. The experiments showed regions of compression and extensional flow across the grounding line, thereby invalidating the boundary layer theory. Steady-state grounding-line positions were found to be dependent on the level of physical model approximation. Resolving grounding lines requires inclusion of membrane stresses, a sufficiently small grid size (<500 m), or subgrid interpolation of the grounding line. The latter still requires nominal grid sizes of <5 km. For larger grid spacings, appropriate parameterizations for ice flux may be imposed at the grounding line, but the short-time transient behaviour is then incorrect and different from models that do not incorporate grounding-line parameterizations. The numerical error associated with predicting grounding-line motion can be reduced significantly below the errors associated with parameter ignorance and uncertainties in future scenarios.

The mechanical properties of glacier beds play a fundamental role in regulating the sensitivity of glaciers to environmental forcing across a wide range of timescales. Glaciers are commonly underlain by deformable till whose mechanical properties and influence on ice flow are not well understood but are critical for reliable projections of future glacier states. Using synoptic-scale observations of glacier motion in different seasons to constrain numerical ice flow models, we study the mechanics of the bed beneath Hofsjökull, a land-terminating ice cap in central Iceland. Our results indicate that the bed deforms plastically and weakens following incipient summertime surface melt. Combining the inferred basal shear traction fields with a Coulomb-plastic bed model, we estimate the spatially distributed effective basal water pressure and show that changes in basal water pressure and glacier accelerations are non-local and non-linear. These results motivate an idealized physical model relating mean basal water pressure and basal slip rate wherein the sensitivity of glacier flow to changes in basal water pressure is inversely related to the ice surface slope.