Search

Over the past two decades, ice loss from the Greenland ice sheet (GrIS) has increased owing to enhanced surface melting and ice discharge to the ocean1-5. Whether continuing increased ice loss will accelerate further, and by how much, remains contentious6-9. A main contributor to future ice loss is the Northeast Greenland Ice Stream (NEGIS), Greenland's largest basin and a prominent feature of fast-flowing ice that reaches the interior of the GrIS10-12. Owing to its topographic setting, this sector is vulnerable to rapid retreat, leading to unstable conditions similar to those in the marine-based setting of ice streams in Antarctica13-20. Here we show that extensive speed-up and thinning triggered by frontal changes in 2012 have already propagated more than 200 km inland. We use unique global navigation satellite system (GNSS) observations, combined with surface elevation changes and surface speeds obtained from satellite data, to select the correct basal conditions to be used in ice flow numerical models, which we then use for future simulations. Our model results indicate that this marine-based sector alone will contribute 13.5-15.5 mm sea-level rise by 2100 (equivalent to the contribution of the entire ice sheet over the past 50 years) and will cause precipitous changes in the coming century. This study shows that measurements of subtle changes in the ice speed and elevation inland help to constrain numerical models of the future mass balance and higher-end projections show better agreement with observations.

The Southern Ocean surrounding Antarctica is a region that is key to a range of climatic and oceanographic processes with worldwide effects, and is characterised by high biological productivity and biodiversity. Since 2013, the International Bathymetric Chart of the Southern Ocean (IBCSO) has represented the most comprehensive compilation of bathymetry for the Southern Ocean south of 60°S. Recently, the IBCSO Project has combined its efforts with the Nippon Foundation - GEBCO Seabed 2030 Project supporting the goal of mapping the world's oceans by 2030. New datasets initiated a second version of IBCSO (IBCSO v2). This version extends to 50°S (covering approximately 2.4 times the area of seafloor of the previous version) including the gateways of the Antarctic Circumpolar Current and the Antarctic circumpolar frontal systems. Due to increased (multibeam) data coverage, IBCSO v2 significantly improves the overall representation of the Southern Ocean seafloor and resolves many submarine landforms in more detail. This makes IBCSO v2 the most authoritative seafloor map of the area south of 50°S.

Floating ice shelves buttress inland ice and curtail grounded-ice discharge. Climate warming causes melting and ultimately breakup of ice shelves, which could escalate ocean-bound ice discharge and thereby sea-level rise. Should ice shelves collapse, it is unclear whether they could recover, even if we meet the goals of the Paris Agreement. Here, we use a numerical ice-sheet model to determine if Petermann Ice Shelf in northwest Greenland can recover from a future breakup. Our experiments suggest that post-breakup recovery of confined ice shelves like Petermann's is unlikely, unless iceberg calving is greatly reduced. Ice discharge from Petermann Glacier also remains up to 40% higher than today, even if the ocean cools below present-day temperatures. If this behaviour is not unique for Petermann, continued near-future ocean warming may push the ice shelves protecting Earth's polar ice sheets into a new retreated high-discharge state which may be exceedingly difficult to recover from.

The land ice contribution to global mean sea level rise has not yet been predicted1 using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models2-8, but primarily used previous-generation scenarios9 and climate models10, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios11,12 using statistical emulation of the ice sheet and glacier models. We find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.

Humboldt Gletscher is a 100-km wide, slow-moving glacier in north Greenland which holds a 19-cm global sea level equivalent. Humboldt has been the fourth largest contributor to sea level rise since 1972 but the cause of its mass loss has not been elucidated. Multi-beam echo sounding data collected in 2019 indicate a seabed 200 m deeper than previously known. Conductivity temperature depth data reveal the presence of warm water of Atlantic origin at 0°C at the glacier front and a warming of the ocean waters by 0.9 ± 0.1°C since 1962. Using an ocean model, we reconstruct grounded ice undercutting by the ocean, combine it with calculated retreat caused by ice thinning to floatation, and are able to fully explain the observed retreat. Two thirds of the retreat are caused by undercutting of grounded ice, which is a physical process not included in most ice sheet models.

Greenland's outlet glaciers have been a leading source of mass loss and accompanying sea-level rise from the Greenland Ice Sheet (GrIS) over the last 25 years. The dynamic component of outlet glacier mass loss depends on both the ice flux through the terminus and the inland extent of glacier thinning, initiated at the ice-ocean interface. Here, we find limits to the inland spread of thinning that initiates at glacier termini for 141 ocean-terminating outlet glaciers around the GrIS. Inland diffusion of thinning is limited by steep reaches of bed topography that we call "knickpoints." We show that knickpoints exist beneath the majority of outlet glaciers but they are less steep in regions of gentle bed topography, giving glaciers in gentle bed topography the potential to contribute to ongoing and future mass loss from the GrIS by allowing the diffusion of thinning far into the ice sheet interior.

The retreat and acceleration of Greenland glaciers since the mid-1990s have been attributed to the enhanced intrusion of warm Atlantic Waters (AW) into fjords, but this assertion has not been quantitatively tested on a Greenland-wide basis or included in models. Here, we investigate how AW influenced retreat at 226 marine-terminating glaciers using ocean modeling, remote sensing, and in situ observations. We identify 74 glaciers in deep fjords with AW controlling 49% of the mass loss that retreated when warming increased undercutting by 48%. Conversely, 27 glaciers calving on shallow ridges and 24 in cold, shallow waters retreated little, contributing 15% of the loss, while 10 glaciers retreated substantially following the collapse of several ice shelves. The retreat mechanisms remain undiagnosed at 87 glaciers without ocean and bathymetry data, which controlled 19% of the loss. Ice sheet projections that exclude ocean-induced undercutting may underestimate mass loss by at least a factor of 2.

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.

The Greenland Ice Sheet is the largest land ice contributor to sea level rise. This will continue in the future but at an uncertain rate and observational estimates are limited to the last few decades. Understanding the long-term glacier response to external forcing is key to improving projections. Here we use historical photographs to calculate ice loss from 1880-2012 for Jakobshavn, Helheim, and Kangerlussuaq glacier. We estimate ice loss corresponding to a sea level rise of 8.1 ± 1.1 millimetres from these three glaciers. Projections of mass loss for these glaciers, using the worst-case scenario, Representative Concentration Pathways 8.5, suggest a sea level contribution of 9.1-14.9 mm by 2100. RCP8.5 implies an additional global temperature increase of 3.7 °C by 2100, approximately four times larger than that which has taken place since 1880. We infer that projections forced by RCP8.5 underestimate glacier mass loss which could exceed this worst-case scenario.

On behalf of the journal, AGU, and the scientific community, the editors would like to sincerely thank those who reviewed the manuscripts for Geophysical Research Letters in 2019. The hours reading and commenting on manuscripts not only improve the manuscripts but also increase the scientific rigor of future research in the field. We particularly appreciate the timely reviews in light of the demands imposed by the rapid review process at Geophysical Research Letters. With the revival of the “major revisions” decisions, we appreciate the reviewers' efforts on multiple versions of some manuscripts. With the advent of AGU's data policy, many reviewers have helped immensely to evaluate the accessibility and availability of data associated with the papers they have reviewed, and many have provided insightful comments that helped to improve the data presentation and quality. We greatly appreciate the assistance of the reviewers in advancing open science, which is a key objective of AGU's data policy. Many of those listed below went beyond and reviewed three or more manuscripts for our journal, and those are indicated in italics.

Bathymetry (seafloor depth), is a critical parameter providing the geospatial context for a multitude of marine scientific studies. Since 1997, the International Bathymetric Chart of the Arctic Ocean (IBCAO) has been the authoritative source of bathymetry for the Arctic Ocean. IBCAO has merged its efforts with the Nippon Foundation-GEBCO-Seabed 2030 Project, with the goal of mapping all of the oceans by 2030. Here we present the latest version (IBCAO Ver. 4.0), with more than twice the resolution (200 × 200 m versus 500 × 500 m) and with individual depth soundings constraining three times more area of the Arctic Ocean (∼19.8% versus 6.7%), than the previous IBCAO Ver. 3.0 released in 2012. Modern multibeam bathymetry comprises ∼14.3% in Ver. 4.0 compared to ∼5.4% in Ver. 3.0. Thus, the new IBCAO Ver. 4.0 has substantially more seafloor morphological information that offers new insights into a range of submarine features and processes; for example, the improved portrayal of Greenland fjords better serves predictive modelling of the fate of the Greenland Ice Sheet.

The sea level contribution of the Antarctic ice sheet constitutes a large uncertainty in future sea level projections. Here we apply a linear response theory approach to 16 state-of-the-art ice sheet models to estimate the Antarctic ice sheet contribution from basal ice shelf melting within the 21st century. The purpose of this computation is to estimate the uncertainty of Antarctica's future contribution to global sea level rise that arises from large uncertainty in the oceanic forcing and the associated ice shelf melting. Ice shelf melting is considered to be a major if not the largest perturbation of the ice sheet's flow into the ocean. However, by computing only the sea level contribution in response to ice shelf melting, our study is neglecting a number of processes such as surface-mass-balance-related contributions. In assuming linear response theory, we are able to capture complex temporal responses of the ice sheets, but we neglect any self-dampening or self-amplifying processes. This is particularly relevant in situations in which an instability is dominating the ice loss. The results obtained here are thus relevant, in particular wherever the ice loss is dominated by the forcing as opposed to an internal instability, for example in strong ocean warming scenarios. In order to allow for comparison the methodology was chosen to be exactly the same as in an earlier study (Levermann et al., 2014) but with 16 instead of 5 ice sheet models. We include uncertainty in the atmospheric warming response to carbon emissions (full range of CMIP5 climate model sensitivities), uncertainty in the oceanic transport to the Southern Ocean (obtained from the time-delayed and scaled oceanic subsurface warming in CMIP5 models in relation to the global mean surface warming), and the observed range of responses of basal ice shelf melting to oceanic warming outside the ice shelf cavity. This uncertainty in basal ice shelf melting is then convoluted with the linear response functions of each of the 16 ice sheet models to obtain the ice flow response to the individual global warming path. The model median for the observational period from 1992 to 2017 of the ice loss due to basal ice shelf melting is 10.2 mm, with a likely range between 5.2 and 21.3 mm. For the same period the Antarctic ice sheet lost mass equivalent to 7.4mm of global sea level rise, with a standard deviation of 3.7mm (Shepherd et al., 2018) including all processes, especially surface-mass-balance changes. For the unabated warming path, Representative Concentration Pathway 8.5 (RCP8.5), we obtain a median contribution of the Antarctic ice sheet to global mean sea level rise from basal ice shelf melting within the 21st century of 17 cm, with a likely range (66th percentile around the mean) between 9 and 36 cm and a very likely range (90th percentile around the mean) between 6 and 58 cm. For the RCP2.6 warming path, which will keep the global mean temperature below 2 °C of global warming and is thus consistent with the Paris Climate Agreement, the procedure yields a median of 13 cm of global mean sea level contribution. The likely range for the RCP2.6 scenario is between 7 and 24 cm, and the very likely range is between 4 and 37 cm. The structural uncertainties in the method do not allow for an interpretation of any higher uncertainty percentiles.We provide projections for the five Antarctic regions and for each model and each scenario separately. The rate of sea level contribution is highest under the RCP8.5 scenario. The maximum within the 21st century of the median value is 4 cm per decade, with a likely range between 2 and 9 cm per decade and a very likely range between 1 and 14 cm per decade.

On behalf of the journal, AGU, and the scientific community, the Editors would like to sincerely thank those who reviewed manuscripts for Geophysical Research Letters in 2018. The hours reading and commenting on manuscripts not only improves the manuscripts but also increases the scientific rigor of future research in the field. We particularly appreciate the timely reviews, in light of the demands imposed by the rapid review process at Geophysical Research Letters. With the revival of the “major revisions” decisions, we appreciate the reviewers' efforts on multiple versions of some manuscripts. Many of those listed below went beyond and reviewed three or more manuscripts for our journal, and those are indicated in italics. In total, 4,484 referees contributed to 7,557 individual reviews in journal. Thank you again. We look forward to the coming year of exciting advances in the field and communicating those advances to our community and to the broader public.

We reconstruct the mass balance of the Greenland Ice Sheet using a comprehensive survey of thickness, surface elevation, velocity, and surface mass balance (SMB) of 260 glaciers from 1972 to 2018. We calculate mass discharge, D, into the ocean directly for 107 glaciers (85% of D) and indirectly for 110 glaciers (15%) using velocity-scaled reference fluxes. The decadal mass balance switched from a mass gain of +47 ± 21 Gt/y in 1972-1980 to a loss of 51 ± 17 Gt/y in 1980-1990. The mass loss increased from 41 ± 17 Gt/y in 1990-2000, to 187 ± 17 Gt/y in 2000-2010, to 286 ± 20 Gt/y in 2010-2018, or sixfold since the 1980s, or 80 ± 6 Gt/y per decade, on average. The acceleration in mass loss switched from positive in 2000-2010 to negative in 2010-2018 due to a series of cold summers, which illustrates the difficulty of extrapolating short records into longer-term trends. Cumulated since 1972, the largest contributions to global sea level rise are from northwest (4.4 ± 0.2 mm), southeast (3.0 ± 0.3 mm), and central west (2.0 ± 0.2 mm) Greenland, with a total 13.7 ± 1.1 mm for the ice sheet. The mass loss is controlled at 66 ± 8% by glacier dynamics (9.1 mm) and 34 ± 8% by SMB (4.6 mm). Even in years of high SMB, enhanced glacier discharge has remained sufficiently high above equilibrium to maintain an annual mass loss every year since 1998.

Earlier large-scale Greenland ice sheet sea-level projections (e.g., those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of the initMIP-Greenland intercomparison exercise is to compare, evaluate and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project - phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of 1) the initial present-day state of the ice sheet and 2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly), and should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap, but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

We use updated drainage inventory, ice thickness, and ice velocity data to calculate the grounding line ice discharge of 176 basins draining the Antarctic Ice Sheet from 1979 to 2017. We compare the results with a surface mass balance model to deduce the ice sheet mass balance. The total mass loss increased from 40 ± 9 Gt/y in 1979-1990 to 50 ± 14 Gt/y in 1989-2000, 166 ± 18 Gt/y in 1999-2009, and 252 ± 26 Gt/y in 2009-2017. In 2009-2017, the mass loss was dominated by the Amundsen/Bellingshausen Sea sectors, in West Antarctica (159 ± 8 Gt/y), Wilkes Land, in East Antarctica (51 ± 13 Gt/y), and West and Northeast Peninsula (42 ± 5 Gt/y). The contribution to sea-level rise from Antarctica averaged 3.6 ± 0.5 mm per decade with a cumulative 14.0 ± 2.0 mm since 1979, including 6.9 ± 0.6 mm from West Antarctica, 4.4 ± 0.9 mm from East Antarctica, and 2.5 ± 0.4 mm from the Peninsula (i.e., East Antarctica is a major participant in the mass loss). During the entire period, the mass loss concentrated in areas closest to warm, salty, subsurface, circumpolar deep water (CDW), that is, consistent with enhanced polar westerlies pushing CDW toward Antarctica to melt its floating ice shelves, destabilize the glaciers, and raise sea level.

Large uncertainties in model parameterizations and input data sets make projections of future sea level rise contributions of outlet glaciers challenging. Here we introduce a novel technique for weighing large ensemble model simulations that uses information of key observables. The approach is robust to input errors and yields calibrated means and error estimates of a glacier's mass balance. We apply the technique to Jakobshavn Isbræ, using a model that includes a dynamic calving law, and closely reproduce the observed behavior from 1985 to 2018 by forcing the model with ocean temperatures only. Our calibrated projection suggests that the glacier will continue to retreat and contribute about 5.1 mm to eustatic sea level rise by 2100 under present-day climatic forcing. Our analysis shows that the glacier's future evolution will strongly depend on the ambient oceanic setting.

We report the discovery of a large impact crater beneath Hiawatha Glacier in northwest Greenland. From airborne radar surveys, we identify a 31-kilometer-wide, circular bedrock depression beneath up to a kilometer of ice. This depression has an elevated rim that cross-cuts tributary subglacial channels and a subdued central uplift that appears to be actively eroding. From ground investigations of the deglaciated foreland, we identify overprinted structures within Precambrian bedrock along the ice margin that strike tangent to the subglacial rim. Glaciofluvial sediment from the largest river draining the crater contains shocked quartz and other impact-related grains. Geochemical analysis of this sediment indicates that the impactor was a fractionated iron asteroid, which must have been more than a kilometer wide to produce the identified crater. Radiostratigraphy of the ice in the crater shows that the Holocene ice is continuous and conformable, but all deeper and older ice appears to be debris rich or heavily disturbed. The age of this impact crater is presently unknown, but from our geological and geophysical evidence, we conclude that it is unlikely to predate the Pleistocene inception of the Greenland Ice Sheet.

Satellite imagery reveals flowstripes on Foundation Ice Stream parallel to ice flow, and meandering features on the ice-shelf that cross-cut ice flow and are thought to be formed by water exiting a well-organised subglacial system. Here, ice-penetrating radar data show flow-parallel hard-bed landforms beneath the grounded ice, and channels incised upwards into the ice shelf beneath meandering surface channels. As the ice transitions to flotation, the ice shelf incorporates a corrugation resulting from the landforms. Radar reveals the presence of subglacial water alongside the landforms, indicating a well-organised drainage system in which water exits the ice sheet as a point source, mixes with cavity water and incises upwards into a corrugation peak, accentuating the corrugation downstream. Hard-bedded landforms influence both subglacial hydrology and ice-shelf structure and, as they are known to be widespread on formerly glaciated terrain, their influence on the ice-sheet-shelf transition could be more widespread than thought previously.

Understanding grounding line dynamics is critical for projecting glacier evolution and sea level rise. Observations from satellite radar interferometry reveal rapid grounding line migration forced by oceanic tides that are several kilometers larger than predicted by hydrostatic equilibrium, indicating the transition from grounded to floating ice is more complex than previously thought. Recent studies suggest seawater intrusion beneath grounded ice may play a role in driving rapid ice loss. Here, we investigate its impact on the evolution of Petermann Glacier, Greenland, using an ice sheet model. We compare model results with observed changes in grounding line position, velocity, and ice elevation between 2010 and 2022. We match the observed retreat, speed up, and thinning using 3‐km‐long seawater intrusion that drive peak ice melt rates of 50 m/yr; but we cannot obtain the same agreement without seawater intrusion. Including seawater intrusion in glacier modeling will increase the sensitivity to ocean warming. Plain Language Summary: Relatively warm seawater melts marine‐terminating glaciers from below. Recent observations suggest that seawater flows below grounded ice at high tide. The presence of seawater at this boundary, referred to as seawater intrusion, has the potential to increase glacier mass loss. We test this hypothesis on Petermann Glacier, Greenland, using an ice sheet flow model. We run the model to reconstruct the glacier's behavior from 2010 to 2022 with and without seawater intrusion. We compare the results with satellite observations of velocity, grounding line position, and ice thinning. When we use enhanced ice melt rates from kilometer‐scale seawater intrusion, we match the observed retreat, speed up, and thinning. When we do not, the model fails to replicate the observations. Seawater intrusion may play a critical role in glacier evolution. Adding this process to ice flow models will increase their sensitivity to ocean warming and projections of ice mass loss and sea level rise. Key Points: Ice melt caused by seawater intrusion in the grounding zone explains the observed grounding line retreat of Petermann GlacierWithout seawater intrusion in the grounding zone, we do not replicate the full extent of observed retreatIncluding seawater intrusion in the grounding zone increases glacier mass loss [ABSTRACT FROM AUTHOR]

On behalf of the journal, AGU, and the scientific community, the editors of Geophysical Research Letters would like to sincerely thank those who reviewed manuscripts for us in 2023. The hours reading and commenting on manuscripts not only improve the manuscripts, but also increase the scientific rigor of future research in the field. With the advent of AGU's data policy, many reviewers have also helped immensely to evaluate the accessibility and availability of data, and many have provided insightful comments that helped to improve the data presentation and quality. We greatly appreciate the assistance of the reviewers in advancing open science, which is a key objective of AGU's data policy. We particularly appreciate the timely reviews in light of the demands imposed by the rapid review process at Geophysical Research Letters. We received 4,512 submissions in 2023 and 5,112 reviewers contributed to their evaluation by providing 8,587 reviews in total. We deeply appreciate their contributions. Plain Language Summary: Individuals in italics provided three or more reviews for GRL in 2023. Key Points: The editors thank the 2023 peer‐reviewers [ABSTRACT FROM AUTHOR]

The Wordie Ice Shelf-Fleming Glacier system in the southern Antarctic Peninsula has experienced a long-term retreat and disintegration of its ice shelf in the past 50 years. Increases in the glacier velocity and dynamic thinning have been observed over the past two decades, especially after 2008 when only a small ice shelf remained at the Fleming Glacier front. It is important to know whether the substantial further speed-up and greater surface draw-down of the glacier since 2008 is a direct response to ocean forcing, or driven by feedbacks within the grounded marine-based glacier system, or both. Recent observational studies have suggested the 2008-2015 velocity change was due to the ungrounding of the Fleming Glacier front. To explore the mechanisms underlying the recent changes, we use a full-Stokes ice sheet model to simulate the basal shear stress distribution of the Fleming system in 2008 and 2015. This study is part of the first high resolution modelling campaign of this system. Comparison of inversions for basal shear stresses for 2008 and 2015 suggests the migration of the grounding line ĝ1/4 9 km upstream by 2015 from the 2008 ice front/grounding line positions, which virtually coincided with the 1996 grounding line position. This migration is consistent with the change in floating area deduced from the calculated height above buoyancy in 2015. The retrograde submarine bed underneath the lowest part of the Fleming Glacier may have promoted retreat of the grounding line. Grounding line retreat may also be enhanced by a feedback mechanism upstream of the grounding line by which increased basal lubrication due to increasing frictional heating enhances sliding and thinning. Improved knowledge of bed topography near the grounding line and further transient simulations with oceanic forcing are required to accurately predict the future movement of the Fleming Glacier system grounding line and better understand its ice dynamics and future contribution to sea level..

Many glaciers in the Antarctic Peninsula are now rapidly losing mass. Understanding of the dynamics of these fast-flowing glaciers, and their potential future behaviour, can be improved through ice sheet modelling studies. Inverse methods are commonly used in ice sheet models to infer the spatial distribution of a basal friction coefficient, which has a large effect on the basal velocity and ice deformation. Here we use the full-Stokes Elmer/Ice model to simulate the Wordie Ice Shelf-Fleming Glacier system in the southern Antarctic Peninsula. With an inverse method, we infer the pattern of the basal friction coefficient from surface velocities observed in 2008. We propose a multi-cycle spin-up scheme to reduce the influence of the assumed initial englacial temperature field on the final inversion. This is particularly important for glaciers like the Fleming Glacier, which have areas of strongly temperature-dependent deformational flow in the fast-flowing regions. Sensitivity tests using various bed elevation datasets, ice front positions and boundary conditions demonstrate the importance of high-accuracy ice thickness/bed geometry data and precise location of the ice front boundary..

While a lot of attention has been given to the numerical implementation of grounding lines and basal friction in the grounding zone, little has been done about the impact of the numerical treatment of ocean-induced basal melting in this region. Several strategies are currently being employed in the ice sheet modeling community, and the resulting grounding line dynamics may differ strongly, which ultimately adds significant uncertainty to the projected contribution of marine ice sheets to sea level rise. We investigate here several implementations of basal melt parameterization on partially floating elements in a finite-element framework, based on the Marine Ice Sheet-Ocean Model Intercomparison Project (MISOMIP) setup: (1) melt applied only to entirely floating elements, (2) melt applied over all elements that are crossed by the grounding line, and (3) melt integrated partially over the floating portion of a finite element using two different sub-element integration methods. All methods converge towards the same state when the mesh resolution is fine enough. However, (2) and (3) will systematically overestimate the rate of grounding line retreat in coarser resolutions, while (1) converges faster to the solution in most cases. The differences between sub-element parameterizations are exacerbated for experiments with high melting rates in the vicinity of the grounding line and for a Weertman sliding law. As most real-world simulations use horizontal mesh resolutions of several hundreds of meters at best, and high melt rates are generally present close to the grounding lines, we recommend not using (3) to avoid overestimating the rate of grounding line retreat and to carefully assess the impact of mesh resolution and sub-element melt parameterizations on all simulation results.

Antarctica and Greenland hold enough ice to raise sea level by more than 65 m if both ice sheets were to melt completely. Predicting future ice sheet mass balance depends on our ability to model these ice sheets, which is limited by our current understanding of several key physical processes, such as iceberg calving. Large-scale ice flow models either ignore this process or represent it crudely. To model fractured zones, an important component of many calving models, continuum damage mechanics as well as linear fracture mechanics are commonly used. However, these methods have a large number of uncertainties when applied across the entire Antarctic continent because the models were typically tuned to match processes seen on particular ice shelves. Here we present an alternative, statistics-based method to model the most probable zones of the location of fractures and demonstrate our approach on all main ice shelf regions in Antarctica, including the Antarctic Peninsula. We can predict the location of observed fractures with an average success rate of 84 % for grounded ice and 61 % for floating ice and a mean overestimation error rate of 26 % and 20 %, respectively. We found that Antarctic ice shelves can be classified into groups based on the factors that control fracture location.

Estimating the future evolution of the Antarctic Ice Sheet (AIS) is critical for improving future sea level rise (SLR) projections. Numerical ice sheet models are invaluable tools for bounding Antarctic vulnerability; yet, few continental-scale projections of century-scale AIS SLR contribution exist, and those that do vary by up to an order of magnitude. This is partly because model projections of future sea level are inherently uncertain and depend largely on the model's boundary conditions and climate forcing, which themselves are unknown due to the uncertainty in the projections of future anthropogenic emissions and subsequent climate response. Here, we aim to improve the understanding of how uncertainties in model forcing and boundary conditions affect ice sheet model simulations. With use of sampling techniques embedded within the Ice Sheet System Model (ISSM) framework, we assess how uncertainties in snow accumulation, ocean-induced melting, ice viscosity, basal friction, bedrock elevation, and the presence of ice shelves impact continental-scale 100-year model simulations of AIS future sea level contribution. Overall, we find that AIS sea level contribution is strongly affected by grounding line retreat, which is driven by the magnitude of ice shelf basal melt rates and by variations in bedrock topography. In addition, we find that over 1.2 m of AIS global mean sea level contribution over the next century is achievable, but not likely, as it is tenable only in response to unrealistically large melt rates and continental ice shelf collapse. Regionally, we find that under our most extreme 100-year warming experiment generalized for the entire ice sheet, the Amundsen Sea sector is the most significant source of model uncertainty (1032 mm 6σ spread) and the region with the largest potential for future sea level contribution (297 mm). In contrast, under a more plausible forcing informed regionally by literature and model sensitivity studies, the Ronne basin has a greater potential for local increases in ice shelf basal melt rates. As a result, under this more likely realization, where warm waters reach the continental shelf under the Ronne ice shelf, it is the Ronne basin, particularly the Evans and Rutford ice streams, that are the greatest contributors to potential SLR (161 mm) and to simulation uncertainty (420 mm 6spread).

Earth system models (ESMs) are becoming increasingly complex, requiring extensive knowledge and experience to deploy and use in an efficient manner. They run on high-performance architectures that are significantly different from the everyday environments that scientists use to pre- and post-process results (i.e., MATLAB, Python). This results in models that are hard to use for non-specialists and are increasingly specific in their application. It also makes them relatively inaccessible to the wider science community, not to mention to the general public. Here, we present a new software/model paradigm that attempts to bridge the gap between the science community and the complexity of ESMs by developing a new JavaScript application program interface (API) for the Ice Sheet System Model (ISSM). The aforementioned API allows cryosphere scientists to run ISSM on the client side of a web page within the JavaScript environment. When combined with a web server running ISSM (using a Python API), it enables the serving of ISSM computations in an easy and straightforward way. The deep integration and similarities between all the APIs in ISSM (MATLAB, Python, and now JavaScript) significantly shortens and simplifies the turnaround of state-of-the-art science runs and their use by the larger community. We demonstrate our approach via a new Virtual Earth System Laboratory (VESL) website (http://vesl.jpl.nasa.gov , VESL(2017)).

Thwaites Glacier (TG), West Antarctica, has been losing mass and retreating rapidly in the past few decades. Here, we present a study of its calving dynamics combining a two-dimensional flow-band full-Stokes (FS) model of its viscous flow with linear elastic fracture mechanics (LEFM) theory to model crevasse propagation and ice fracturing. We compare the results with those obtained with the higher-order (HO) and the shallow-shelf approximation (SSA) models coupled with LEFM. We find that FS/LEFM produces surface and bottom crevasses that are consistent with the distribution of depth and width of surface and bottom crevasses observed by NASA's Operation IceBridge radar depth sounder and laser altimeter, whereas HO/LEFM and SSA/LEFM do not generate crevasses that are consistent with observations. We attribute the difference to the nonhydrostatic condition of ice near the grounding line, which facilitates crevasse formation and is accounted for by the FS model but not by the HO or SSA models. We find that calving is enhanced when pre-existing surface crevasses are present, when the ice shelf is shortened or when the ice shelf front is undercut. The role of undercutting depends on the timescale of calving events. It is more prominent for glaciers with rapid calving rates than for glaciers with slow calving rates. Glaciers extending into a shorter ice shelf are more vulnerable to calving than glaciers developing a long ice shelf, especially as the ice front retreats close to the grounding line region, which leads to a positive feedback to calving events. We conclude that the FS/LEFM combination yields substantial improvements in capturing the stress field near the grounding line of a glacier for constraining crevasse formation and iceberg calving.

We employ airborne gravity data from NASA's Operation IceBridge collected in 2009–2014 to infer the bathymetry of sub–ice shelf cavities in front of Pine Island, Thwaites, Smith, and Kohler glaciers, West Antarctica. We use a three-dimensional inversion constrained by multibeam echo sounding data offshore and bed topography from a mass conservation reconstruction on land. The seamless bed elevation data refine details of the Pine Island sub–ice shelf cavity, a slightly thinner cavity beneath Thwaites, and previously unknown deep (>1200 m) channels beneath the Crosson and Dotson ice shelves that shallow (500 m and 750 m, respectively) near the ice shelf fronts. These sub–ice shelf channels define the natural pathways for warm, circumpolar deep water to reach the glacier grounding lines, melt the ice shelves from below, and constrain the pattern of past and future glacial retreat.

Identifying fast and robust numerical solvers is a critical issue that needs to be addressed in order to improve projections of polar ice sheets evolving in a changing climate. This work evaluates the impact of using advanced numerical solvers for transient ice-flow simulations conducted with the JPL-UCI Ice Sheet System Model (ISSM). We identify optimal numerical solvers by testing a broad suite of readily available solvers, ranging from direct sparse solvers to preconditioned iterative methods, on the commonly used Ice Sheet Model Intercomparison Project for Higher-Order ice sheet Models benchmark tests. Three types of analyses are considered: mass transport, horizontal stress balance, and incompressibility. The results of the fastest solvers for each analysis type are ranked based on their scalability across mesh size and basal boundary conditions. We find that the fastest iterative solvers are ĝ1/4ĝ€1.5-100 times faster than the default direct solver used in ISSM, with speed-ups improving rapidly with increased mesh resolution. We provide a set of recommendations for users in search of efficient solvers to use for transient ice-flow simulations, enabling higher-resolution meshes and faster turnaround time. The end result will be improved transient simulations for short-term, highly resolved forward projections (10-100 year time scale) and also improved long-term paleo-reconstructions using higher-order representations of stresses in the ice. This analysis will also enable a new generation of comprehensive uncertainty quantification assessments of forward sea-level rise projections, which rely heavily on ensemble or sampling approaches that are inherently expensive.

Knowledge of the ice thickness distribution of glaciers and ice caps is an important prerequisite for many glaciological and hydrological investigations. A wealth of approaches has recently been presented for inferring ice thickness from characteristics of the surface. With the Ice Thickness Models Intercomparison eXperiment (ITMIX) we performed the first coordinated assessment quantifying individual model performance. A set of 17 different models showed that individual ice thickness estimates can differ considerably - locally by a spread comparable to the observed thickness. Averaging the results of multiple models, however, significantly improved the results: on average over the 21 considered test cases, comparison against direct ice thickness measurements revealed deviations on the order of 10 ± 24% of the mean ice thickness (1σ estimate). Models relying on multiple data sets - such as surface ice velocity fields, surface mass balance, or rates of ice thickness change - showed high sensitivity to input data quality. Together with the requirement of being able to handle large regions in an automated fashion, the capacity of better accounting for uncertainties in the input data will be a key for an improved next generation of ice thickness estimation approaches.

Understanding of long-term dynamics of glaciers and ice caps is vital to assess their recent and future changes, yet few long-term reconstructions using ice flow models exist. Here we present simulations of the maritime Hardangerjøkulen ice cap in Norway from the mid-Holocene through the Little Ice Age (LIA) to the present day, using a numerical ice flow model combined with glacier and climate reconstructions. In our simulation, under a linear climate forcing, we find that Hardangerjøkulen grows from ice-free conditions in the mid-Holocene to its maximum extent during the LIA in a nonlinear, spatially asynchronous fashion. During its fastest stage of growth (2300-1300 BP), the ice cap triples its volume in less than 1000 years. The modeled ice cap extent and outlet glacier length changes from the LIA until today agree well with available observations. Volume and area for Hardangerjøkulen and several of its outlet glaciers vary out-of-phase for several centuries during the Holocene. This volume-area disequilibrium varies in time and from one outlet glacier to the next, illustrating that linear relations between ice extent, volume and glacier proxy records, as generally used in paleoclimatic reconstructions, have only limited validity. We also show that the present-day ice cap is highly sensitive to surface mass balance changes and that the effect of the ice cap hypsometry on the mass balance- altitude feedback is essential to this sensitivity. A mass balance shift by +0.5m w.e. relative to the mass balance from the last decades almost doubles ice volume, while a decrease of 0.2 mw.e. or more induces a strong mass balance-altitude feedback and makes Hardangerjøkulen disappear entirely. Furthermore, once disappeared, an additional +0.1m w.e. relative to the present mass balance is needed to regrow the ice cap to its present-day extent. We expect that other ice caps with comparable geometry in, for example, Norway, Iceland, Patagonia and peripheral Greenland may behave similarly, making them particularly vulnerable to climate change.

Melting of the Greenland Ice Sheet has the potential to raise sea level by 7.36 m and is already contributing to global sea level rise at a rate higher than 1 mm yr–1. Computer models are our best tools to make projections of the mass balance of Greenland over the next centuries, but these models rely on bed topography data that remain poorly constrained near glacier termini. Accurate bed topography in the vicinity of calving fronts is critical for numerical models, as the shapes of the glacier bed and of the nearby bathymetry control both the ocean circulation in the fjord and the stability and response of the ice sheet to climate warming. NASA’s Oceans Melting Greenland (OMG) mission is collecting bathymetry data along Greenland fjords at several glacier termini. Here, we show that these measurements are transforming our knowledge of fjord and glacier depths. Using a mass conservation approach, we combine OMG bathymetry with observations of ice velocity and thickness to produce estimates of bed depth and ice thickness across the ice-ocean boundary with unprecedented accuracy and reliability. Our results along the northwest coast of Greenland reveal complex structural features in bed elevation, such as valleys, ridges, bumps, and hollows. These features have important implications for both channeling ice flow toward the continental margin, and for controlling the amount of warm, salty Atlantic Water that reaches the glaciers.

Abstract. Thwaites Glacier (TG), West Antarctica, has been losing mass and retreating rapidly in the past few decades. Here, we present a study of its calving dynamics combining a two-dimensional flowband Full Stokes (FS) model of its viscous flow with linear elastic fracture mechanics (LEFM) theory to model crevasse propagation and ice fracturing. We compare the results with those obtained with the higher-order (HO) and the shallow-shelf approximation (SSA) models coupled with LEFM. We find that FS/LEFM produces surface and bottom crevasses that match the distribution of crevasse depth and width observed from NASA's Operation IceBridge radar depth sounders, whereas HO/LEFM and SSA/LEFM do not generate crevasses that match observations. We attribute the difference to the non-hydrostatic condition of ice near the grounding line, which facilitates crevasse formation, and is accounted for by the FS model but not by the HO or SSA model. We also find that calving is enhanced when pre-existing surface crevasses are present, when the ice shelf is shortened or when the ice shelf front is undercut. The role of undercutting depends on the time scale of calving events. It is more prominent for glaciers with rapid calving rates than glaciers with slow calving rates. Glaciers extending into a shorter ice shelf are more vulnerable to calving than glaciers developing a long ice shelf, especially as the ice front retreats close to the grounding line region, which leads to a positive feedback. We conclude that the FS/LEFM combination yields substantial improvements in capturing the stress field near the grounding line for constraining crevasse formation and iceberg calving.

Basal sliding is a main control on glacier flow primarily driven by water pressure at the glacier base. The ongoing increase in surface melting of the Greenland Ice Sheet warrants an examination of its impact on basal water pressure and in turn on basal sliding. Here we examine the case of Russell Glacier, in West Greenland, where an extensive set of observations has been collected. These observations suggest that the recent increase in melt has had an equivocal impact on the annual velocity, with stable flow on the lower part of the drainage basin but accelerated flow above the Equilibrium Line Altitude (ELA). These distinct behaviors have been attributed to different evolutions of the subglacial draining system during and after the melt season. Here we use a high-resolution subglacial hydrological model forced by reconstructed surface runoff for the period 2008 to 2012 to investigate the cause of these distinct behaviors. We find that the increase in meltwater production at low elevation yields a more efficient drainage system compatible with the observed stagnation of the mean annual flow below the ELA. At higher elevation, the model indicates that the drainage system is mostly inefficient and is therefore strongly sensitive to an increase in meltwater availability, which is consistent with the observed increase in ice velocity.

We reconstruct the history of the Greenland Ice Sheet margin on the Nuussuaq Peninsula in central West Greenland through the Holocene using lake sediment analysis and cosmogenic 10Be exposure dating of the prominent Drygalski Moraines. Erratics perched on bedrock outboard of the Drygalski Moraines constrain local deglaciation to ∼9.9 ± 0.6 ka (n = 2). Three Drygalski Moraine crests yield mean 10Be ages of 8.6 ± 0.4 ka (n = 2), 8.5 ± 0.2 ka (n = 3), and 7.6 ± 0.1 ka (n = 2) from outer to inner. Perched erratics between the inner two moraines average 7.8 ± 0.1 ka (n = 2) and are consistent with the moraine ages. Sediments from a proglacial lake with a catchment area extending an estimated 2 km beneath (inland of) the present ice sheet terminus constrain an ice sheet minimum extent from 5.4 ka to 0.6 ka. The moraine chronology paired with the lake sediment stratigraphy reveals that the ice margin likely remained within ∼2 km of its present position from ∼9.9 to 5.4 ka. This unexpected early Holocene stability, preceded by rapid ice retreat and followed by minimum ice extent between ∼5.4 and 0.6 ka, contrasts with many records of early Holocene warmth and the Northern Hemisphere summer insolation maximum. We suggest ice margin stability may instead be tied to adjacent ocean temperatures, which reached an optimum in the middle Holocene.