CR3.1 | Modelling ice sheets and glaciers
Modelling ice sheets and glaciers
Convener: Sainan Sun | Co-conveners: Fabien Gillet-Chaulet, Stephen Cornford, Mauro Werder, Rabea Sondershaus
Orals
| Tue, 25 Apr, 14:00–17:55 (CEST)
 
Room L2, Wed, 26 Apr, 08:30–10:10 (CEST)
 
Room 1.61/62
Posters on site
| Attendance Mon, 24 Apr, 14:00–15:45 (CEST)
 
Hall X5
Orals |
Tue, 14:00
Mon, 14:00
This session is intended to attract a broad range of ice-sheet and glacier modelling contributions, welcoming applied and theoretical contributions. Theoretical topics that are encouraged are higher-order mechanical models, data inversion and assimilation, representation of other earth sub-systems in ice-sheet models, and the incorporation of basal processes and novel constitutive relationships in these models.
Applications of newer modelling themes to ice-sheets and glaciers past and present are particularly encouraged, in particular those considering ice streams, rapid change, grounding line motion and ice-sheet model intercomparisons.

Orals: Tue, 25 Apr | Room L2

Chairperson: Sainan Sun
14:00–14:05
14:05–14:25
|
EGU23-2804
|
ECS
|
solicited
|
On-site presentation
Theo Clayton, Ravindra Duddu, Martin Siegert, and Emilio Martinez-Pañeda

Surface crevasses are predominately mode I fractures that penetrate tens of metres deep into grounded glaciers and floating ice shelves. However, elevated surrounding temperatures have resulted in the production of surface meltwater, which accumulates in neighbouring crevasses and applies additional tensile stresses to crack walls. This process is known as hydrofracture; and if sufficient, can promote full thickness crevasse propagation, and lead to iceberg calving events. Net ablation of ice sheets has become of great concern, as it has become the largest contributor to sea-level rise. To overcome the limitations of empirical and analytical approaches to crevasse predictions, we here propose a thermo-dynamically consistent phase field damage model to simulate damage growth in both ice sheets and floating ice shelves using the finite element method.

The model incorporates the three elements needed to mechanistically simulate hydrofracture of surface and basal crevasses: (i) a constitutive description of glacier flow, incorporating the non-linear viscous rheology of ice using Glen’s flow law, (ii) a phase field formulation capable of capturing cracking phenomena of arbitrary complexity, such as 3D crevasse interaction, and (iii) a poro-damage mechanics representation to account for the role of meltwater pressure on crevasse growth. To assess the suitability of the method, we simulated the propagation of surface and basal crevasses within grounded glaciers and floating ice shelves and compared the predicted crevasse depths with analytical methods such as linear elastic fracture mechanics and the Nye zero stress
method, with results showing good agreement for idealised conditions.

References
T. Clayton, R. Duddu, M. Siegert, E. Martínez-Pañeda, A stress-based poro-damage phase field
model for hydrofracturing of creeping glaciers and ice shelves, Engineering Fracture Mechanics
272 (2022) 108693

 

How to cite: Clayton, T., Duddu, R., Siegert, M., and Martinez-Pañeda, E.: Phase field modelling of glacial crevasses subject to meltwater-driven hydro-fracture, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2804, https://doi.org/10.5194/egusphere-egu23-2804, 2023.

Data driven modelling
14:25–14:35
|
EGU23-2607
|
ECS
|
On-site presentation
Samuel Cook, Guillaume Jouvet, Romain Millan, and Antoine Rabatel

Mountain glaciers are a major source of sea-level rise and also represent an important freshwater resource in many mountainous regions. Thus, accurate estimations of their thickness and, therefore, the total ice volume are important both in predicting and mitigating the global and local effects of climate change. However, to date, only 2% of the world’s glaciers outside the ice sheets have any thickness data, due to the logistical difficulties of obtaining such measurements, creating a large and policy-relevant scientific gap.

The recent development of a global-scale ice-velocity dataset, however, provides an ideal opportunity to fill this gap and determine ice thickness across the 98% of glaciers for which no thickness data is available. This can be done by inverting an ice-dynamics model to solve for the ice thickness. For accurate thickness results, this needs to be a full-Stokes model, but such a model is far too computationally cumbersome to apply on a global scale, and simpler, quicker methods usually based on the shallow ice approximation (SIA) are too inaccurate, particularly where sliding dominates glacier motion. The only attempt that has been made to leverage the global velocity dataset to retrieve ice thickness has, though, used the SIA, simply because higher-order approaches are not computationally realistic at this scale. Consequently, most of the widely-used global glacier models have made no concerted attempt to invert for global ice thickness, owing to these limitations.

As an additional related problem, failing to fully assimilate ice-velocity data into an ice-flow model necessarily introduces a shock when initialising prognostic glacier simulations, resulting in model glaciers and predictions that may diverge substantially from their real-world counterparts.

As a solution to these problems, we present results from a deep-learning-driven inversion model that emulates the performance of state-of-the-art full-Stokes models at a thousandth of the computational cost. This model, by solving an optimisation problem, can fully use and assimilate all available input datasets (surface velocity and topography, ice thickness, etc.) as components of its cost function to simultaneously invert for and optimise multiple control parameters (here, we focus on ice thickness). This approach also gives us the possibility of using the same ice-velocity field for inversion and forward modelling, reducing the magnitude of the shock inherent in traditional modelling approaches. With a view to a large-scale application to all the world’s 200,000 glaciers, we present initial thickness-inversion results for the relatively well-documented European Alps to help constrain model parameters and provide a test bed for extension to other glaciated regions, with initial extension to the Caucasus and the Southern Alps.

How to cite: Cook, S., Jouvet, G., Millan, R., and Rabatel, A.: Glacier ice thickness estimation using deep-learning-driven emulation of Stokes, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2607, https://doi.org/10.5194/egusphere-egu23-2607, 2023.

14:35–14:45
|
EGU23-15264
|
On-site presentation
Sebastian Rosier, Christopher Bull, Wai Woo, and Hilmar Gudmundsson

Through their role in buttressing upstream ice flow, Antarctic ice shelves play an important part in regulating future sea level change. Reduction in ice-shelf buttressing caused by increased ocean-induced melt along their undersides is now understood to be one of the key drivers of ice loss from the Antarctic Ice Sheet. However, despite the importance of this forcing mechanism, most ice-sheet simulations currently rely on simple melt-parametrisations of this ocean-driven process since a fully coupled ice-ocean modelling framework is prohibitively computationally expensive. Here, we provide an alternative approach that can capture the greatly improved physical description of this process provided by large-scale ocean-circulation models over currently employed melt-parameterisations, but with trivial computational expense.  This new method brings together deep learning and physical modelling to develop a deep neural network framework, MELTNET, that can emulate ocean model predictions of sub-ice shelf melt rates. We train MELTNET on synthetic geometries, using the NEMO ocean model as a ground-truth in lieu of observations to provide melt rates both for training and to evaluate the performance of the trained network. We show that MELTNET can accurately predict melt rates for a wide range of complex synthetic geometries, with a normalized root mean squared error of 0.11m/yr compared to the ocean model. MELTNET calculates melt rates several orders of magnitude faster than the ocean model and outperforms more traditional parameterisations for 96% of geometries tested. Furthermore, we find MELTNET's melt rate estimates show sensitivity to established physical relationships such as changes in thermal forcing and ice shelf slope. This study demonstrates the potential for a deep learning framework to calculate melt rates with almost no computational expense, that could in the future be used in conjunction with an ice sheet model to provide predictions for large-scale ice sheet models.

How to cite: Rosier, S., Bull, C., Woo, W., and Gudmundsson, H.: Predicting ocean-induced ice-shelf melt rates using deep learning, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-15264, https://doi.org/10.5194/egusphere-egu23-15264, 2023.

14:45–14:55
|
EGU23-5414
|
On-site presentation
James Maddison, Beatriz Recinos, and Daniel Goldberg

Ice sheet models, calibrated using observational data, provide a means of projecting our current best state of the knowledge of the system state into the future, so as to obtain information about possible future behaviour. However it is important to be able to estimate the uncertainty associated with these projections. The problem of quantifying ice sheet parametric uncertainty is considered, focusing specifically on the problem of quantifying the posterior uncertainty in inferred basal sliding and rheology coefficients. These measures of uncertainty are projected forwards in time to obtain measures of uncertainty in future quantities of interest. Automated code generation and automated differentiation tools are utilised, leading to an extensible approach. The role of the observational error model in defining parametric uncertainty is considered.

How to cite: Maddison, J., Recinos, B., and Goldberg, D.: A framework for quantifying parametric ice sheet model uncertainty, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-5414, https://doi.org/10.5194/egusphere-egu23-5414, 2023.

14:55–15:05
|
EGU23-6969
|
ECS
|
Virtual presentation
Anna Crawford, Greg Crocker, Derek Mueller, and Jesse Smith

The calving of ice tongues and ice shelves can generate large, tabular icebergs that have climatological implications given their role in dispersing freshwater from the Greenland and Antarctic ice sheets. These ‘ice islands’ also pose potential risk to marine industry. It is therefore critical that influential deterioration mechanisms be accurately represented in simulations of ice island drift and deterioration, both for risk mitigation in offshore industry and for climatological studies that are focused on the Polar Regions. The majority of ice island deterioration is the result of sidewall erosion, and specifically that which results from waterline wave-erosion leading to ram growth and buoyancy-forced fracture. This study therefore focuses on the inclusion of the buoyancy-driven “footloose” calving mechanism (Wagner et al., 2014) in simulations of ice island length and areal change. Using size and lineage information of ice islands tracked in the Canadian Ice Island Drift, Deterioration and Detection (CI2D3) Database, we quantitatively assess the performance of the footloose calving model by simulating the deterioration of 172 ice islands. The mean model error was +15 (+/- 400) m over 20 d and increased to +401 (+1400/-800) m for simulations that ran to 80 d. The performance of the footloose calving model is a substantial improvement when compared to simulations that did not include this calving mechanism. For example, a thermal-melt model had mean errors of -252 and -1403 m at 20 and 80 d of simulation, respectively, and the mean error of a zero-melt model was -281 and -1545 m over the same time periods. We also present a new approach to modelling ice island areal change resulting from footloose calving. This simple, two-parameter approach simulates discrete footloose calving events and adjusts the ice island surface area to maintain a constant aspect ratio. Mean model error remained under 1 km2 over 80 d of simulation, showing that the model performs well over numerous months. Using the CI2D3 Database, we were able to conduct the first large-scale assessment of the footloose model’s performance in simulating change to the ice island length dimension. The morphological data included in the database also provided the opportunity to develop an approach for modelling areal deterioration resulting from footloose calving events. The model assessments would benefit from more observations of long-duration ice island tracks, as there were a limited number of ice islands that were tracked in the CI2D3 Database for over 40 d. Future work can look to implement the presented approaches in operational and climatological modelling while the iceberg modelling community also develops an approach to simulate larger-scale ice island fracture.

How to cite: Crawford, A., Crocker, G., Mueller, D., and Smith, J.: Developments to 2D modelling of ice island deterioration using the CI2D3 Database, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6969, https://doi.org/10.5194/egusphere-egu23-6969, 2023.

15:05–15:15
|
EGU23-12705
|
ECS
|
On-site presentation
|
Reuben W. Nixon-Hill, Daniel Shapero, Colin J. Cotter, and David A. Ham
We present methods and tools which significantly improve the ability to estimate quantities and fields which are difficult to directly measure at large scales, such as the fluidity of ice, using point data sources from remote sensing. Our tools work with both sparse and dense point data with estimated quantities and fields becoming more accurate as the number of measurements are increased. These are often used as input variables to mathematical models that are used to make predictions so improving their accuracy is of vital importance.
 
The tool we introduce, Firedrake, generates highly optimised code for solving PDEs via the finite element method. It is easy to use, can be integrated with other python libraries, and is scalable for use on high performance computers. As proof, we highlight Icepack, a Python library for solving the equations of motion of glacier flow, which is written using Firedrake. Firedrake has an interface to the dolfin-adjoint/pyadjoint tool which allows these data assimilation problems to be solved with just a few lines of code by automatically generating the adjoint system of PDEs.
 
By carefully considering the nature of finite element method solutions, we show, using Firedrake and Icepack, how the choice of misfit functional (i.e. objective function) significantly impacts our inferred fields. This required the development of new infrastructure not previously available in these automated code generation tools.

How to cite: Nixon-Hill, R. W., Shapero, D., Cotter, C. J., and Ham, D. A.: Point Data Assimilation in Firedrake and Icepack, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-12705, https://doi.org/10.5194/egusphere-egu23-12705, 2023.

Greenland
15:15–15:25
|
EGU23-16383
|
ECS
|
On-site presentation
Ivan Utkin, Ludovic Räss, and Mauro Werder

Antarctic and Greenland ice sheets lose most of their mass by a few corridors of rapidly flowing ice. These ice conveyor belts constitute fast drainage routes whose flow velocities are undoubtedly sensitive to climate perturbations directly impacting sea-level. Observations suggest the ice is rather sliding than flowing, the key being where sliding is accommodated. Commonly, sliding occurs at the ice-bedrock interface, but recent studies favour englacial sliding to explain data from Western Margin of Greenland. 

Our aim is to demonstrate that the mechanisms controlling the spontaneous formation of englacial sliding explains the transition from slow flowing to fast sliding ice over Greenland. We employ a new thermo-mechanical ice flow model to predict thermally activated creep instability leading to the spontaneous rearrangement of ice motion in three dimensions. Accurately resolving these nonlinear interactions on regional to ice sheet scales requires high spatial and temporal resolution which can only be achieved using a supercomputer.

We present a new thermo-mechanical ice flow model, FastIce.jl, that is capable of predicting the evolution of ice sheet at unprecedented scale. The model uses the full-Stokes formulation for the ice flow and the enthalpy method for describing the polythermal ice behaviour. FastIce.jl uses GPU acceleration for solving the flow equations, resulting in close to ideal scaling in distributed computing benchmarks. We compare our simulation results with other full-Stokes models, and present the results of simulating the 120x120 km regions of Greenland ice sheet at 10m resolution. High resolution allowed us to capture the transition from slow to sliding flow regimes without any simplifying assumptions.

How to cite: Utkin, I., Räss, L., and Werder, M.: Large-scale thermo-mechanical modelling of Greenland ice sheet, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16383, https://doi.org/10.5194/egusphere-egu23-16383, 2023.

15:25–15:35
|
EGU23-8265
|
ECS
|
On-site presentation
Matt Trevers, Tony Payne, and Stephen Cornford

Jakobshavn Isbræ (JI), on the West coast of Greenland, is one of the fastest flowing outlet glaciers of the Greenland Ice Sheet, draining 7% of the ice sheet area. Since the late 1990s it has dramatically accelerated, thinned and retreated in a series of phases alternating with periods of quiescence. It exhibits strong seasonal variations in flow speed and calving front position. Between 2012 and 2015, JI attained its point of furthest retreat and flow speeds in excess of 17 km/yr. Since 2016 it has modestly thickened, concurrent with deceleration and readvance of the calving front.

The very fast flow and strong annual and interannual variability present significant challenges for ice sheet modellers. We model the evolution of JI between 2009 and 2018 using the BISICLES ice sheet model. The standard modelling technique of assimilating surface velocity observations to infer a power law basal friction coefficient for a snapshot in time fails to account for rapidly changing basal conditions, underestimating the annual variability. We implement a time-series inverse method in which regular velocity observations are assimilated throughout the study period to produce a time-evolving basal friction coefficient. This method is able to reproduce the large annual variations in flow speed much more accurately than the static method.

This reliance on regular observations to drive the model poses a problem for future projections. We compare a range of sliding laws applied with the normal snapshot inverse method. A modern regularized Coulomb friction sliding law is better able to reproduce JI’s annual variations in flow speed due to its ability to modulate the basal friction in response to movement of the grounding line. As a result, it may be a more appropriate choice of sliding law for modelling the future evolution of fast-flowing outlet glaciers.

How to cite: Trevers, M., Payne, T., and Cornford, S.: A comparison of inverse methods and basal sliding laws applied to a hindcast model of Jakobshavn Isbræ from 2009 to 2018., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8265, https://doi.org/10.5194/egusphere-egu23-8265, 2023.

Mountain and Valley Glacers
15:35–15:45
|
EGU23-9021
|
ECS
|
On-site presentation
Ludovic Räss, Christophe Ogier, Ivan Utkin, Mauro Werder, Andreas Bauder, and Daniel Farinotti

Rapid climate modifications perturb the long-term dynamic equilibrium of many natural systems. Polar and high-altitude regions such as alpine environments represent locations where perturbations such as glacier collapse features, become visible. Glacier collapse features are characterized by a circular depressions on the ice surface, are bounded by low-angle crevasses and are the surface expression of a cavity developing most often over a subglacial channel, commonly occurring at the glacier snout. Understanding the physical processes governing the collapse feature dynamics is essential to assess hazards and processes related to them, such as, rapid glacier length variations, snout collapses and sudden blockage of the subglacial drainage system.

Field observations from an on-going collapse feature developing at the snout of Rhonegletscher (Switzerland) in Summer 2022 suggest mechanical failure of ice lamellas from the underlying cavity roof to drive the collapse. In order to test this hypothesis, namely mechanical failure to drive glacier collapse features, we developed full-Stokes 2D and 3D mechanical models implementing a temperature and pressure dependent visco-elasto-plastic rheology. We use the extensive dataset from Rhonegletscher to constrain the numerical models to predict possible failure patterns as function of increasing cavity size. We use vertical displacement located in the centre of the collapse feature to validate our models. Preliminary results show the formation of tension failure patterns on the ice surface at locations coinciding with the low-angle circular crevasses. The model results will advance our understanding of the physics of collapse features and provide predictive tools to assess future occurrences and their related risks.

How to cite: Räss, L., Ogier, C., Utkin, I., Werder, M., Bauder, A., and Farinotti, D.: Mechanical failure to drive the glacier collapse feature at Rhonegletscher, Switzerland, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9021, https://doi.org/10.5194/egusphere-egu23-9021, 2023.

Coffee break
Chairperson: Fabien Gillet-Chaulet
16:15–16:25
|
EGU23-8982
|
ECS
|
On-site presentation
Erik Young, Gwenn Flowers, Hester Jiskoot, and Daniel Gibson

Glacier surges produce iconic valley-scale folds which encode a history of polyphase deformation resulting from switches between quiescent and surging flow. The folding is passive, resulting from disturbances to ice foliation during surging flow, and subsequently altered during quiescent flow. We investigate the kinematic evolution of these kilometre-scale folds, using Elmer/Ice, by modelling folds through multiple surge cycles using idealized synthetic glacier confluence configurations, and identifying how differences in glacier flow regimes imprint themselves onto three-dimensional fold geometry. The surge and quiescent phases are simulated by changing the basal conditions of one of the tributaries, and matching the scale of velocity variations observed in temperate glacier surges. We determine fold geometry using a particle tracking algorithm applied to the modelled velocity fields, where, mimicking a medial moraine, vertically-spaced particles are injected at the flow unit confluence and advected downglacier. Using structural analysis of the model outputs, we present an archetype of kinematic evolution that describes the transition from cylindrical, vertically plunging, gentle folds emplaced during the surge phase, to complex, non-cylindrical, depth-varying folds following multiple cycles of surging and quiescent flow. The initial fold geometry is controlled by longitudinal and lateral shear stress regimes during surging, while fold evolution is governed primarily by lateral shearing after emplacement. We examine the sensitivity of fold geometry to valley geometry, glacier dynamics, and mass balance. Finally, we illustrate the potential of our approach to reconstruct more complex fold geometries as observed in nature, by applying it to a large surge-type glacier in the St. Elias Mountains of Northern Canada.

How to cite: Young, E., Flowers, G., Jiskoot, H., and Gibson, D.: Modelling the kinematic evolution of valley-scale folding in surge-type glaciers using Elmer/Ice, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8982, https://doi.org/10.5194/egusphere-egu23-8982, 2023.

16:25–16:35
|
EGU23-10981
|
ECS
|
On-site presentation
Dongqi Lin, Marwan Katurji, and Heather Purdie

Observations worldwide have shown that glaciers are receding. Thinning snow cover can result in crevasses becoming exposed at the glacier surface for longer time periods. Crevasses can increase surface roughness, change the surface wind flow fields, and impact the air temperatures within and outside the crevasses (Purdie et al., 2022). Therefore, more can be learnt about regarding the impact of crevasses on energy exchange with glacier atmospheric boundary layer. In order to understand and investigate the atmosphere-crevasse energy exchange, we carried out numerical Large Eddy Simulation (LES) experiments using the PALM model system 6.0 with crevasse-resolving grid spacings less than 1 m. The PALM model system 6.0 has been used for atmosphere and marine boundary layer studies to understand complex processes of atmospheric dynamics and energy balance. Our preliminary results show that the air inside a crevasse does not cool as fast as the air outside a crevasse resulting in net warming from the crevasse relative to the glacier surface. These results agree with the field study conducted at Tasman Glacier (Purdie et al., 2022), and confirm that crevasses could lead to heat storage and increased melting. During the daytime, air temperature inside a crevasse could be 1 °C higher than the air above the glacier surface. After sunset, the presence of the crevasse entrains and traps the warm air such that the centre of the crevasse could still be warmer than the glacier surface during the first half of the evening.  During windy evenings, our results show that turbulent heat exchange associated with eddy entrainment leads to exchange of air mass between the crevasse and the glacier surface, which then causes glacier surface warming. Our preliminary simulations only included one crevasse, while future work will include a field of crevasses to investigate the impact of crevasses in a more realistic environment. This study highlights the importance of including crevasses energy balance in glacier modelling, neglecting which would lead to significant bias in snow melt and mass balance estimations.

 

Reference:

Purdie, H., et al. (2022). "Variability in the vertical temperature profile within crevasses at an alpine glacier." Journal of Glaciology: 1-5.

How to cite: Lin, D., Katurji, M., and Purdie, H.: Impact of crevasses on surface energy balance at an alpine glacier, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10981, https://doi.org/10.5194/egusphere-egu23-10981, 2023.

Processes
16:35–16:45
|
EGU23-3993
|
ECS
|
Virtual presentation
|
Johannes Feldmann, Ricarda Winkelmann, and Anders Levermann

Ice-shelf pinning points such as ice rises or ice rumples can have an important role in regulating the ice discharge of marine outlet glaciers. For instance, the observed gradual ungrounding of the ice shelf of West Antarctica's Thwaites Glacier from its last pinning points diminishes the buttressing effect of the ice shelf and thus contributes to the destabilization of the outlet. Here we use an idealized experimental setting to simulate the response of an Antarctic-type marine outlet glacier to a successive ungrounding of its ice shelf from a pinning point. This is realized by perturbing steady states by a step-wise lowering of the pinning point, which induces a buttressing reduction. After the complete detachment of the ice shelf from the pinning point the perturbation is reversed, i.e., the pinnning point is incrementally elevated toward its initial elevation. First results show that the glacier retreat down the landward down-sloping (retrograde) bed, induced by the loss in buttressing, can be reversible in case of a relatively flat retrograde bed slope. For steeper slopes, glacier retreat and re-advance show a hysteretic behavior. Thus, if the bed depression is sufficiently deep, the glacier does not recover from its fully retreated state even for pinning-point elevations that are higher than the initial elevation.

How to cite: Feldmann, J., Winkelmann, R., and Levermann, A.: Hysteresis of idealized marine outlet glaciers under variation of pinning-point buttressing, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3993, https://doi.org/10.5194/egusphere-egu23-3993, 2023.

16:45–16:55
|
EGU23-2908
|
On-site presentation
Olga Sergienko and Marianne Haseloff
The results of numerical studies suggest that the choice of sliding laws affect simulated dynamics of the grounding lines. Sliding laws that depend on the effective pressure imply that basal shear vanishes at the grounding line. Using analytical and numerical approaches we investigate what effects vanishing basal shear play in the dynamics and stability of the buttressed marine ice sheets. Our results show that the steady-state configurations of marine ice sheets with vanishing and non-vanishing basal shear do not differ greatly. In contrast, the time-variant behaviours of such marine ice sheets are drastically different. In response to a stochastic temporal variability in submarine melting, marine ice sheets with vanishing basal shear can exhibit unstoppable retreat while marine ice sheets with non-zero basal shear in the vicinity of the grounding line can exhibit intermittent advance and retreat. These results suggest that basal conditions in the vicinity of the grounding lines exhibit strong control of their dynamic behaviour. They also suggest the importance of the temporal variability of the basal shear in the dynamics of buttressed marine ice sheets.

How to cite: Sergienko, O. and Haseloff, M.: Basal controls of the grounding line dynamics of buttresses marine ice sheets, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2908, https://doi.org/10.5194/egusphere-egu23-2908, 2023.

16:55–17:05
|
EGU23-2844
|
On-site presentation
Daniel Martin, Samuel Kachuck, Joanna Millstein, and Brent Minchew

Ice is a non-Newtonian fluid whose rheology is typically described using Glen's flow law, a power-law relationship between stress and strain rate with a stress exponent, n, generally taken to be 3. Recent observation-based work suggests that a more accurate choice for the Glen’s law exponent in high-strain regions like ice shelves may be n=4, implying that ice viscosity is more sensitive to changes in stress than is generally assumed. The implications of a higher stress exponent for ice sheet models and their projections of ice sheet response to climate forcing are unclear and likely to be complex. Rheological parameters, such as ice viscosity, are fundamental to ice sheet dynamics and influence the evolution of marine ice sheets. 

 

Here, we present work that explores the rheological parameter space within the idealized MSIMIP+ marine ice sheet configuration using the BISICLES model. We explore the impacts of increasing  the stress exponent from 3 to 4, highlighting the considerable changes to the ice sheet system caused by increasing the stress exponent. Beyond dynamic changes in the ice behavior, changes become necessary to the other flow law parameters generally computed during initialization. For example, it may be that viscosity modifiers typically interpreted as “damage” may instead be indications of mismatches in rheology.  This study underscores the dynamic sensitivity of glacial ice to changes in the rheological parameters and calls attention to the key variables influencing ice sheet evolution. 

 

How to cite: Martin, D., Kachuck, S., Millstein, J., and Minchew, B.: Examining the sensitivity of ice sheet models to updates in rheology (n=4), EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2844, https://doi.org/10.5194/egusphere-egu23-2844, 2023.

17:05–17:15
|
EGU23-11460
|
ECS
|
On-site presentation
Hugo Rousseau, Johan Gaume, Lars Blatny, and Martin P. Lüthi

Understanding glaciers evolution is of major concern to evaluate their contribution to sea level rise in the context of global warming. Among the various processes involved in glacier dynamics, fractures like the calving of ice at the front of marine terminating glacier and crevasse formation affect the stress state, potentially modifying the glacier’s velocity. Crevasses also impact the melting rate. The fractures alter the roughness of the ground, increasing the amount of absorbed radiation, and open new networks in which the meltwater is likely to penetrate deeper toward the glacier bed.

In this work we propose to model fractures in glacier based on finite strain elastoplasticity, using the Material Point Method (Wolper et al. 2021): we solve the classical governing equations for ice deformation in an Eulerian-Lagrangian framework and we use a strain softening Drucker-Prager constitutive model to simulate plasticity. Thanks to the Lagrangian part of the model, the fractures appear explicitly where high levels of total plastic deformation are reached. The behaviour of a glacier flowing over a step in the bedrock is investigated. The simulations show that crevasse patterns appear with regular spacing between the fractures. We perform a parametric study to determine which parameters affect the length of these patterns and potential dimensionless numbers are outlined.

How to cite: Rousseau, H., Gaume, J., Blatny, L., and Lüthi, M. P.: Modelling discontinuities in ice flow using the Material Point Method and elastoplasticity, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11460, https://doi.org/10.5194/egusphere-egu23-11460, 2023.

17:15–17:25
|
EGU23-11076
|
ECS
|
On-site presentation
Gabriela Clara Racz, Kevin Mirng En Yeo, Aly Thobani, Sean Henry, Maryam Zarrinderakht, Camilo Rada, and Christian Schoof

Commonly, the parts of the glacier bed that are hydraulically connected to the surface experience significant diurnal variations in water pressure, in response to cycles of surface melting. Closely spaced points on the bed often exhibit nearly identical temporal variations in water pressure, suggesting that they are connected not only to the surface but to each other through conduits along the bed. This behaviour is typically observed directly through instrumented boreholes drilled to the glacier bed. A ‘switching event’ occurs when one of a pair of boreholes abruptly changes from being connected, in the sense of exhibiting the same diurnal oscillations as the other borehole, to being disconnected, or vice versa. A switching event is indicative of a connection through a subglacial conduit being closed, or opened, and therefore provides a limited but highly specific window into the evolution of subglacial conduits and permeability.

However, in most subglacial drainage models, conduits are not represented individually but averaged over a small area of the bed to produce a macroporous continuum representation as a ‘water sheet’, quantified by a mean conduit depth h. The most common assumption is that the water sheet consists of linked cavities and that these open due to basal sliding over bed roughness, and close due to viscous creep (e.g. Hewitt, 2011). Within that framework, the simplest mechanism for a switching event is that a connection is established or closed when the sheet thickness h passes through some percolation threshold hc (Rada and Schoof, 2018).

We want to test whether the observed switching events can be explained by that mechanism, which in turn implies that two conditions must be met: water sheet depth indeed evolves according to a competition between opening due to basal sliding and creep closure, and that a simple threshold in h suffices to capture the geometric complexity involved in creating or closing connections at the bed.

In a large dataset of borehole water-pressure time series, we identify borehole pairs that exhibit strong evidence of switching behaviour. We assume that switching events can be described by the evolution of a water sheet, with connections between boreholes being opened and closed as sheet thickness passes through a threshold value as described above. We use the switching event catalogue we have created to invert for parameters in the sheet evolution model using a binary indicator function for connectedness to compute the model data mismatch in the absence of any other direct measures of sheet thickness.

This procedure allows us to capture the majority of observed switching events with plausible parameter values. The exception is a set of short-lived periods of connectedness characterized by switching events that are clustered in space and time. In a complementary study (Racz et al, 2023 in prep.), we, therefore, investigate if this class of switching events can instead be explained by an alternative mechanism in which the sudden resumption of surface water supply, following a period of snow cover, drives the propagation of a hydrofracture (e.g. Tsai and Rice 2010, 2012).

How to cite: Racz, G. C., Yeo, K. M. E., Thobani, A., Henry, S., Zarrinderakht, M., Rada, C., and Schoof, C.: Physical processes driving 'switching events', EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11076, https://doi.org/10.5194/egusphere-egu23-11076, 2023.

Antarctica
17:25–17:35
|
EGU23-11337
|
On-site presentation
Benoit Urruty, Olivier Gagliardini, Fabien Gillet-Chaulet, Gael Durand, and Mondher Chekki

The stability of the grounding lines of Antarctica is a fundamental question in term of sea level rise The strong mass loss of the Antarctic Ice Sheet (AIS) in recent years has raised concerns about the possibility of the ice sheet reaching a tipping point, beyond which it would experience rapid and irreversible loss of mass. Such a tipping point could be triggered by a combination of external forcing factors, including continued warming of the ocean and atmosphere, as well as changes in ice sheet dynamics. As shown by Urruty et al. (in review), the current mass loss and retreat is mainly due to external forcing such as melt induced by the ocean and current grounding lines are not yet engaged in an unstable retreat. But if forcing remains similar or increases, some irreversible and fast mass loss may occur as a result of grounding lines crossing a tipping point.

As part of the TiPACCs project, we are conducting experiments to evaluate the stability of the grounding lines of the AIS in the future. Building on the stability experiment described in Urruty et al. (in review), we are using the same initial state created with Elmer/Ice to perform a new set of experiments. In these experiments, we are applying large-amplitude perturbations to the grounding line of a steady-state AIS by increasing ocean temperature by 1°C, 3°C, and 5°C for periods ranging from 20 to 100 years in order to push the grounding line far from its current position. After the perturbation is removed, we then apply 80 years of constant forcing to see if an unstable position is reached. These experiments will help us better understand the stability of the grounding line at different positions that could be reached in the near future if current observed forcing trends continue.

How to cite: Urruty, B., Gagliardini, O., Gillet-Chaulet, F., Durand, G., and Chekki, M.: Large-amplitude perturbation experiments to assess the unstable behaviour of AIS in the near future, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11337, https://doi.org/10.5194/egusphere-egu23-11337, 2023.

17:35–17:45
|
EGU23-13025
|
On-site presentation
G. Hilmar Gudmundsson, Jan De Rydt, Sebastian Rosier, Jowan Barnes, Daniel Goldberg, and Mathieu Morlighem

We use numerical modelling to address several questions related to the future evolution of Thwaites Glacier over the next 50 years. The importance of Thwaites Ice Shelf for upstream grounded flow is investigated by quantifying the buttressing stresses along the grounding line. Removing the ice shelf changes the stress regime along the grounding line by less than 20%. This change is small compared to many, if not most, grounding lines of the Antarctic Ice Sheet, and much smaller than corresponding changes for the neighboring Pine Island and Pope, Smith and Kohler glaciers.  Transient ice-flow modelling experiments show that mass loss from Thwaites Glacier over the next 50 years is insignificantly affected by removal of the ice shelf. We then explore the consequences of the proposed marine ice-cliff instability for Thwaites Glacier. For recently proposed calving laws, where the calving rate increases sharply with cliff height, we do not observe an onset of an unstable calving front retreat. Further numerical modelling experiments for future climatic forcing scenarios will be presented, including uncertainty quantification. Interactions between the ice and the ocean are studied using a coupled ice+ocean modelling framework. As shown before in several studies, we find when simulating its future evolution, that Thwaites Glacier can enter unstable periods of self-enhancing retreat. This appears to be a very robust result, and this behavior is found in all model runs, including coupled ice+ocean simulations.

How to cite: Gudmundsson, G. H., De Rydt, J., Rosier, S., Barnes, J., Goldberg, D., and Morlighem, M.: The Future of Thwaites Glacier, West Antarctica., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-13025, https://doi.org/10.5194/egusphere-egu23-13025, 2023.

17:45–17:55
|
EGU23-9198
|
ECS
|
On-site presentation
Cyrille Mosbeux, Nicolas Jourdain, Olivier Gagliardini, Peter Råback, and Adrien Gilbert

Ice mass loss from Antarctic Ice Sheet is increasing, accelerating its contribution to global sea level rise. In the Amundsen Sea sector, recent observations of rapid ice-shelf thinning and grounding line retreat have been attributed to increased basal melting driven by inflows of warm Circumpolar Deep Water. However, recent studies have shown that basal melting alone might not be sufficient to explain the recent acceleration, retreat and thinning of the outlet glaciers in the sector.

As part of the European Horizon 2020 research project PROTECT ­— that assesses and projects changes in the land-based cryosphere to produce robust projections of SLR we conduct numerical simulations to determine the role of damage on changes observed over the last two decades in the Amundsen Sea Sector. More particularly, we use a Stokes flow formulation combined with a Continuum Damage Mechanics model of the open-source ice flow model Elmer/Ice to simulate the ice flow evolution. We initialize our ice sheet model with data assimilation methods using 1996 observations of surface velocities as well as a corrected geometry based on the current ice-sheet geometry and the ice thickness rates of change observed over the past 20 years. From this initial state, we run forward simulations over 20 years with and without damage mechanics, and compare the model evolution to observed surface velocities and ice thickness rates of change, as well as observations of grounding line positions. Our results shed light on the importance of damage in the evolution of the region, in particular an acceleration of several hundred meters per year due to the decreasing buttressing effect of the ice shelves triggered by the increasing damage.

How to cite: Mosbeux, C., Jourdain, N., Gagliardini, O., Råback, P., and Gilbert, A.: On the effect of damage on the recent changes in the Amundsen Sea Sector, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9198, https://doi.org/10.5194/egusphere-egu23-9198, 2023.

Orals: Wed, 26 Apr | Room 1.61/62

Chairperson: Rabea Sondershaus
08:30–08:40
|
EGU23-11798
|
ECS
|
On-site presentation
Yanjun Li, Javier Blasco Navarro, Frank Pattyn, and Gang Qiao

Ice shelves around the coastal portion of the Antarctic ice sheet are sensitive indicators of climate change. The thinning of ice shelves diminishes buttressing, promotes longitudinal spreading, and increases ice flux across the grounding line, leading to accelerated glacier discharge into the ocean. Thwaites Glacier in the Amundsen Sea Embayment is one of the fastest-changing outlet glaciers in Antarctica. Damaged areas with densely distributed crevasses and open fractures on Thwaites Glacier are key to future ice shelf stability, grounding line retreat and sea level contribution. The damage feedback processes should be taken into consideration when simulating the evolution of Thwaites Glacier using ice sheet models.

Here, we add the continuum damage mechanics approach to the F.ETISh/Kori ice flow model, to simulate the present-day and near future behavior of the ice sheet and ice shelf system, including brittle ice physics. The damage field is described by equating it to the total crevasse depths used in Nick et al. (2010) and Sun et al. (2017). 100 years simulations under present-day climate conditions with and without damage in different scenarios have been conducted, and the change in ice velocity, ice thickness, the grounding line retreat and the sea level contribution of Thwaites Glacier have been analysed. Moreover, the change in ice velocity along four ice flow profiles in the first 20-year simulation has been analysed and the impact of damage on velocity has been assessed by comparing the simulated velocity fields with the observations (e.g., the MEaSUREs and ITS_LIVE ice velocity products).

Results indicate that damage drastically increases the ice velocity over the ice shelves and weakens them as such that grounding line retreat ~18 km after 20 years, accelerates (~1.5 times) compared to the observed increase in flow speed and contributes around an order of magnitude to the sea level rise in 50 years. Change in ice velocity profiles of Thwaites Glacier also show that local damage may overestimate ice velocity, especially in the grounded ice near the grounding line, while it underestimates observed ice flow when local damage is omitted. Through a series of further sensitivity experiments, an analysis on the timing and magnitude of damage has been carried out to gauge the current and near future state of the Thwaites glacier basin.

 

How to cite: Li, Y., Navarro, J. B., Pattyn, F., and Qiao, G.: Modelling the impacts of ice damage on the response of Thwaites Glacier, West Antarctica, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11798, https://doi.org/10.5194/egusphere-egu23-11798, 2023.

08:40–08:50
|
EGU23-11300
|
ECS
|
On-site presentation
Tim van den Akker, William H. Lipscomb, Gunter R. Leguy, Willem Jan van de Berg, and Roderik van de Wal

Projections of sea level rise are subject to large uncertainties in the contribution of the Antarctic Ice Sheet (AIS), as it is unclear how AIS dynamics will evolve over time. Ice sheet models use spin-up techniques to initialize the ice sheet to the present-day state. Previous attempts using the Community Ice Sheet Model (CISM) assumed that the ice sheet is in equilibrium at the end of the spin-up. This assumption limits the contribution of model drift, but does not match observations and might bias future projections.

 

For this reason, we have incorporated present-day thickness change rates from Smith et al. (2020) in our Antarctic spin-ups. As in previous spin-ups, we tune basal friction coefficients beneath grounded ice, and ocean temperatures beneath floating ice, to match observed present-day thickness. In the new procedure, CISM is also forced to match thickening and thinning rates, with the surface mass balance (SMB) adjusted to allow the AIS to maintain the observed thickness and grounding-line locations. This technique improves the modelled velocities in regions with substantial thinning. When the SMB adjustments are removed, the modelled ice sheet exhibits the observed thickness change rates. We use this initialised state to project AIS evolution without additional forcing (‘committed climate change’). For a range of parameter settings, this causes Thwaites glacier to collapse irreversibly, without further ocean warming. The time of initiation of collapse is sensitive to model parameters, but once initiated the collapse is largely complete within two to three centuries. The sensitivity tests are carried out for a range of parameters related to basal sliding and various ocean warming scenarios.

How to cite: van den Akker, T., Lipscomb, W. H., Leguy, G. R., van de Berg, W. J., and van de Wal, R.: New Antarctic spin-up method results in committed Thwaites glacier collapse, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11300, https://doi.org/10.5194/egusphere-egu23-11300, 2023.

08:50–09:00
|
EGU23-13533
|
ECS
|
On-site presentation
David Bett, Alexander Bradley, Rosie Williams, Paul Holland, Robert Arthern, and Daniel Goldberg

The Amundsen Sea Sector has some of the highest thinning rates of ice shelves in Antarctica, thought to be driven by high, but interannually variable, ocean driven melt rates. This thinning can lead to increased ice flow speeds, eventually leading to sea level rise. To fully represent these processes and other feedbacks, a fully coupled ice/ocean model must be used. Therefore, a fully synchronous mass conservative coupled ice-sheet/ocean model of the Amundsen Sea Sector has been developed. This new coupled model builds upon previous coupling developments and involves coupling of the WAVI ice-sheet model to the 3D ocean model MITgcm, via the Streamice ice-sheet model. Coupled model projections are presented, examining ice grounding line retreat rates and ice mass loss, along with ocean driven melt rate evolution. The sensitivity of these results to ocean forcings is shown, specifically the thickness of the relatively warm Circumpolar Deep Water layer and its variability. In addition, we discuss the impact of the present-day initialisation and tuning of the coupled model.

How to cite: Bett, D., Bradley, A., Williams, R., Holland, P., Arthern, R., and Goldberg, D.: Fully synchronous coupled ice/ocean modelling of future changes in the Amundsen Sea Sector, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-13533, https://doi.org/10.5194/egusphere-egu23-13533, 2023.

09:00–09:10
|
EGU23-3216
|
On-site presentation
Felicity McCormack, Bernd Kulessa, and Jason Roberts

Vanderford Glacier is one of the fastest retreating glaciers in East Antarctica, with approximately 18.6 km of grounding line retreat since 1996. Together with the Totten Glacier, the Vanderford Glacier is a key outlet glacier of the Aurora Subglacial Basin (ASB), which contains approximately 7 m of global sea level equivalent, of which ~3.5 m is vulnerable to ocean driven melting, and is rapidly losing mass. While the Totten Glacier currently discharges almost twice as much ice as the Vanderford Glacier, sediment records from the Sabrina and Knox Coast Sectors indicate that the Vanderford Glacier has had sedimentation rates over twice that at Totten in the past. Here, we examine the current flow configuration between Vanderford and Totten Glaciers and drivers of it, including interactions between the subglacial topography, hydraulic potential, climate, and ice sheet dynamics. We use the Ice-sheet and Sea-level System Model (ISSM) under experiments of heightened ocean warming concentrated at Vanderford Glacier, and heightened surface mass balance at Totten Glacier, to show that the present-day flow configuration between the Totten and Vanderford Glaciers is tenuous. Rerouting towards Vanderford Glacier could occur under even minor changes in surface elevation at both glaciers. Such rerouting potentially exposes large parts of the underbelly of the ASB to enhanced ocean-driven ice shelf melting in the event of rapid retreat of Vanderford Glacier, with implications for global sea level rise.

How to cite: McCormack, F., Kulessa, B., and Roberts, J.: Controls on the flow configuration of Vanderford Glacier, East Antarctica, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-3216, https://doi.org/10.5194/egusphere-egu23-3216, 2023.

09:10–09:20
|
EGU23-11803
|
ECS
|
On-site presentation
James ONeill, tamsin Edwards, and Lauren Gregoire

The warm Pliocene was a period of comparable atmospheric carbon dioxide concentrations to modern, but with sea levels up to ~20 m higher. High Pliocene sea level implies collapse of the West Antarctic ice sheet, and mass loss from East Antarctica. Modelling studies have sought to reproduce Pliocene deglaciation, and use sea level reconstructions as a constraint on future projections, despite their large uncertainties. We simulated the Pliocene Antarctic ice sheet under warm Pliocene climate with the BISICLES ice sheet model, capturing grounding line and ice stream dynamics down to 4 km resolution. Our perturbed parameter ensemble approach explores uncertainties in basal sliding, surface mass balance processes, bedrock-ice sheet interactions, ice shelf basal melt sensitivity to ocean forcing and choice of climate model. We simulated a mean Antarctic sea level contribution of 1.85 m and a range of -15.90 to 28.27 m, largely driven by uncertainty in the perturbed basal sliding parameter. We applied a joint calibration, combining a Pliocene Antarctic sea level contribution range and a comparison of regional grounding line with reconstructed Pliocene retreat. This reduced the mean to 1.46 m. The calibration reduced the simulated range by a factor of ~4, and was more effective in reducing uncertainty than comparing to sea level reconstructions alone. Further ensembles explored initial condition uncertainty, and the impact of perturbing the control. The Pliocene initial condition was tested for a subset of main ensemble simulations (mean = -2.35 m), increasing the mean contribution by 11.45 m with all simulations passing the joint calibration. We perturbed the control simulation for the same subset of ensemble members. This increased the mean Antarctic contribution by 7.78 m, and by 8.45 m in combination with the two Pliocene data constraints. We demonstrate a modelling framework that captures important interactions between the ice sheet and other components of the Earth system, whilst being efficient for ensemble studies. Moreover, we used two Pliocene data constraints to rule out ensemble members. This Pliocene-calibrated modelling framework can be run under future climate scenarios, to reduce uncertainty in projections of Antarctica’s long-term contribution to sea level under anthropogenic climate change.

How to cite: ONeill, J., Edwards, T., and Gregoire, L.: Pliocene Antarctic ice sheet model ensembles with joint constraints from reconstructed sea level and margin retreat, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-11803, https://doi.org/10.5194/egusphere-egu23-11803, 2023.

09:20–09:30
|
EGU23-12764
|
ECS
|
On-site presentation
Clara Henry, Reinhard Drews, Clemens Schannwell, Vjeran Višnjević, Inka Koch, Heiko Spiegel, Leah Muhle, Olaf Eisen, Daniela Jansen, Steven Franke, and Paul Bons

The geometry of englacial isochrones is a product of the past and present ice velocity field and is useful for our understanding of steady-state ice flow dynamics, flow regime re-organisation, and calibration of models. Ice rises contain various flow regimes (divide flow, flank flow, and grounding zones) on small spatial scales, meaning they are ideal locations to study ice-flow dynamics and stratigraphy to constrain model parameters. We run full Stokes, thermo-mechanically coupled simulations of Derwael Ice Rise in East Antarctica and simulate the three-dimensional stratigraphy of the ice rise and the surrounding ice shelf using the finite element model Elmer/Ice. Over the ice rise, we derive the accumulation rate from internal reflection horizons and use RACMO2.3 surface mass balance data over the surrounding ice shelf. Simulations are run for Glen's flow law exponents of n=3 and n=4 with appropriate values derived for the Arrhenius law.

To calibrate the model, comparisons are made with the BedMachine surface elevation and density-adjusted internal reflection horizons observed in many transects recorded by AWI’s ultra-wide band radar covering the divide, the flanks, and the grounding zones. To understand ice flow dynamics where the velocity field of the ice rise and the ice shelf converge in the compressive and shear zones, we analyse the modelled englacial stress and strain rate fields. Our results allow us to investigate isochronal structures where observed internal reflection horizons are too steep or obscured to be adequately picked up by radar. A comparison between the model and observed fracturing can be used to infer threshold stress and strain rates for fracture initiation. These simulations are a blueprint for the full Stokes, three-dimensional modelling of ice rises and have further relevance in the study of three-dimensional influences on Raymond arch evolution, the constrained coupling of the anisotropy equations, comparisons with ice core data and the automated inference of ice flow parameters from internal reflection horizons.

How to cite: Henry, C., Drews, R., Schannwell, C., Višnjević, V., Koch, I., Spiegel, H., Muhle, L., Eisen, O., Jansen, D., Franke, S., and Bons, P.: Modelling the three-dimensional stratigraphy of an ice rise, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-12764, https://doi.org/10.5194/egusphere-egu23-12764, 2023.

Methods and model comparisons
09:30–09:40
|
EGU23-7264
|
On-site presentation
James Jordan and Frank Pattyn

Ice shelf calving is responsible for roughly half the mass lost by Antarctic Ice Shelves and is of vital importance for determining ice-shelf stability. Despite this, it is currently poorly represented in numerical models and as such has seen limited inclusion in numerical simulations of the future Antarctic Ice Sheet. A first step towards improving this situation is assessing the capabilities and limitations of current numerical models regarding calving. The Calving Model Intercomparison Project (CaMIP) is, therefore, being undertaken to address this issue as part of the EU funded Horizon 2020 project PROTECT.

We present here an overview of the project, as well as preliminary results from the first round of experiments by a range of participating modelling groups from across the cryospheric community.

How to cite: Jordan, J. and Pattyn, F.: Calving Model Intercomparison Project (CaMIP): Overview and preliminary results, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-7264, https://doi.org/10.5194/egusphere-egu23-7264, 2023.

09:40–09:50
|
EGU23-6947
|
ECS
|
On-site presentation
Cruz Garcia-Molina, Fabien Gillet-Chaulet, Mondher Chekki, Gael Durand, Olivier Gagliardini, and Nicolas Jourdain

Ice-calving plays a major role in the mass balance of the water-ending glaciers. Thus, it is crucial to have a well-adapted calving law for simulations over long periods. Due to its dependence on several physical parameters, this phenomenon is usually poorly parametrized in long-term numerical simulations. A worldwide model intercomparison project on ice damage and calving, CalvingMIP (see https://github.com/JRowanJordan/CalvingMIP/wiki), is proposed as part of the European project, PROTECT. The CalvingMIP project aims to evaluate the uncertainties in modelling the ice and to provide recommendations to improve the calving laws in the ice-sheet models. This intercomparison project consists of five experiments using two topographic profiles: a hill and Thule bathymetry. For the first phase, a steady-state configuration is implemented for a fixed calving position. In the second phase: the front velocity is prescribed, forcing the front to advance and, then, to retreat. Finally, the last experiment aims to test a realistic calving law. We study the experiments of this configuration by using the community finite element code, Elmer/Ice (see http://elmerice.elmerfem.org/). We study the front evolution using a level-set function (φ), defined as a signed distance to the front. Here, we present the results obtained with our model for this intercomparison experiment, discuss the sensitivity to different physical and numerical parameters, and its application to a realistic configuration.

How to cite: Garcia-Molina, C., Gillet-Chaulet, F., Chekki, M., Durand, G., Gagliardini, O., and Jourdain, N.: Elmer/Ice results on the CalvingMIPintercomparison project using a level-set function, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6947, https://doi.org/10.5194/egusphere-egu23-6947, 2023.

09:50–10:00
|
EGU23-548
|
ECS
|
On-site presentation
Kevin Hank, Lev Tarasov, and Elisa Mantelli

Some ice sheets and glaciers experience long quiescent periods interspersed with short periods of rapid ice advance, such as the binge-purge-type cycling hypothesized to be associated with Heinrich Events. Modeling ice stream activation/de-activation, however, is numerically challenging given the relatively abrupt changes at surge onset and the high ice velocities. In spite of this, a number of high-profile modeling papers have explored Heinrich events and ice surges, but generally with very limited consideration of numerical aspects. Here we test the ability of the 3D Glacial Systems Model (GSM) and Parallel Ice Sheet Model (PISM) to simulate binge-purge-type surges and explore the stability of the simulations with respect to relevant numerical and discretization uncertainties. 

We find surge characteristics exhibit a resolution dependency but converge at higher horizontal grid resolutions (order 5 km). In accordance with theoretical and experimental work, our model results suggest that the thermal activation of basal sliding should start below the pressure melting point. A resolution-dependent basal temperature ramp for the thermal activation of basal sliding as well as a subglacial hydrology model can reduce the discrepancies between high and coarse horizontal grid resolutions. Furthermore, incorporating a bed thermal and at least a minimal complexity subglacial hydrology model significantly affects surge characteristics and is, therefore, essential for modeling large-scale ice stream cycling.

How to cite: Hank, K., Tarasov, L., and Mantelli, E.: Towards confidence in numerical modeling of ice stream cycling, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-548, https://doi.org/10.5194/egusphere-egu23-548, 2023.

10:00–10:10
|
EGU23-14849
|
ECS
|
On-site presentation
André Löfgren, Josefin Ahlkrona, Thomas Zwinger, Christian Helanow, and Denis Cohen

Ice flow models often suffer from numerical instabilities that restricts time-step sizes. For higher-order models this constitutes a severe bottleneck. We present a method for increasing the largest stable time step in full Stokes models, allowing for a significant speed-up of simulations.  This type of stabilisation was originally developed for mantle-convection simulations and is here extended to ice flow problems. The method is mimicking an implicit solver but the computational cost per time step is nearly as low as for an explicit solver. As it only consists of adding a stabilisation term to the gravitational force in the full Stokes equations, it is very easy to implement. We test the method using both Elmer/Ice and FEniCS on artificial glaciers with varying bedrock roughness, slip rate and surface inclination, as well as on a real world case.

How to cite: Löfgren, A., Ahlkrona, J., Zwinger, T., Helanow, C., and Cohen, D.: Increasing the largest stable time-step size in ice flow models, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14849, https://doi.org/10.5194/egusphere-egu23-14849, 2023.

Posters on site: Mon, 24 Apr, 14:00–15:45 | Hall X5

Chairperson: Mauro Werder
X5.213
|
EGU23-16754
Javier Blasco, Yanjun Li, and Frank Pattyn

As stated in the latest IPCC report, sea level will continue to rise at the end of this century and most likely well beyond, depending on future emission pathways. The Antarctic ice sheet plays an important role, as it is the largest ice sheet and thus the largest source of water storage on Earth. However, projections for Antarctica from ice-sheet models yield very mixed results due to ice-sheet-related processes that are difficult to assess. One of the main sources of uncertainty is the stability of floating ice shelves. Although ice shelves do not directly contribute to sea-level rise, they have been shown to play an important role, as they modulate the grounded ice flow via their buttressing effect. Therefore, it is necessary to assess the stability of ice shelves in a warmer climate to make more accurate predictions and define safe trajectory scenarios. Satellite images show the formation of crevasse in regions with a high deformation rate. These crevasses weaken the stability of the ice shelf, as damage enhances inland ice acceleration and promotes further shearing and retreat. However, most continental-scale ice-sheet models do not account for ice shelf damage and its consequent potential feedback mechanisms. Part of this statement is due to the fact that ice shelves at coarse resolutions show low stability to damage implementation even in simple domains. Here we force a three-dimensional ice-sheet-shelf model with various damage formulations from the literature. Given the high uncertainty in damage formation and propagation, several parameters affecting the stability of the ice shelf are evaluated. Experiments are performed in different domains to test their influence in simple and symmetric cases, such as MISMIP+, as well as in the Amundsen-Sea Embayment. Our results highlight the importance of further research on ice damage, as it has strong implications for projections but is poorly accounted for in ice-sheet models.

How to cite: Blasco, J., Li, Y., and Pattyn, F.: Assessing the role of ice-shelf damage on a three-dimensional ice-sheet model, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-16754, https://doi.org/10.5194/egusphere-egu23-16754, 2023.

X5.214
|
EGU23-15896
Jan De Rydt, Kaitlin Naughten, and Hilmar Gudmundsson

Ice-shelf cavities in the Amundsen Sea are expanding as the ice thins and grounding lines retreat. To sustain ice-shelf thinning, whilst accommodating a (up to) 2-fold increase in ice flux across the grounding line, basal melt and calving rates must respond. Changes in far-field ocean temperature are often evoked to explain the sustained thinning, but there is no indication of significant trends in ocean properties over the observational period. On the other hand, internal feedbacks between changes in ice-shelf geometry and basal melt could play a role in driving glacier retreat, but these processes remain poorly understood. Here we explore such melt-geometry feedbacks using a coupled ice+ocean model of the Amundsen Sea glaciers. Under present-day ocean conditions, all glaciers, including Pine Island and Thwaites glaciers, continue to retreat, and we find significant trends in ice-shelf melt rates despite the absence of trends in the far-field ocean forcing. For example, melt rates for Pine Island Glacier double over a 50-year period under constant ocean conditions. For all cavities, the trend in melt can be attributed to a reconfiguration of the ocean circulation beneath the ice shelf, in response to changes in cavity geometry. We argue that these melt-geometry feedbacks may play an important role in the evolution of the Amundsen Sea glaciers, and they should be adequately captured in numerical simulations – something simple basal-melt parameterizations are currently unable to do.

How to cite: De Rydt, J., Naughten, K., and Gudmundsson, H.: Feedbacks between basal melt and cavity geometry in coupled ice+ocean simulations of the Amundsen Sea glaciers, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-15896, https://doi.org/10.5194/egusphere-egu23-15896, 2023.

X5.215
|
EGU23-1401
|
ECS
Niall Gandy, Gemma Ives, Gwyneth Rivers, and Lachlan Astfalck

After running a large ensemble of palaeo ice sheet model runs, it is common to either rank the simulations or determine which simulations are an acceptable match to observations and which are ruled out. This task requires human judgement, which is usually left only to the authors of the research. Tools have been developed to compare ice sheet simulations to empirical reconstructions numerically, but even this approach requires some human input on values for match thresholds.

An alternative is to use expert elicitation to identify “good” ice sheet simulations. Expert elicitation normally captures expert’s beliefs in the form of a probability distribution; for something as complicated as ice sheet geometry this is much too arduous a task. Instead, we propose to elicit binary classifications of “good” and “bad” and find descriptions of plausible ice sheets through probabilistic inverse modelling. Experts can consider empirical ice sheet reconstructions, but also “soft-knowledge” about the sectors of the ice sheet it is most important to match, margin shapes considered to be glaciologically plausible, and an idea of the reasonable best-reconstruction a model will be able to provide. By seeking the input of many experts, it is possible to both lighten the task load of determining the quality of 100-1000s of simulations, and gain a wisdom of the crowd benefit to the results. Just like any other method of ranking ice sheet simulations, this method requires human judgement; in this case more explicitly than usual.

We are seeking expert input to rank an existing ensemble of North American Ice Sheet simulations. By asking experts at EGU 2023 to spend 3-5 minutes sorting simulations using an online interface we will build up an average community view on which LGM North American Ice Sheet simulations are “good”. This will provide a community resource to compare future ice sheet simulations against that is a justifiable representation of academic expert knowledge, adding to the current arsenal of model-data intercomparison tools.

How to cite: Gandy, N., Ives, G., Rivers, G., and Astfalck, L.: Ice Sheet Speed-dating: Using Expert Elicitation to identify “good” simulations of the LGM North American Ice Sheets, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-1401, https://doi.org/10.5194/egusphere-egu23-1401, 2023.

X5.216
|
EGU23-4200
|
ECS
|
Helen Werner, Dirk Scherler, Ricarda Winkelmann, and Guillaume Jouvet

For reconstructing paleoclimate or studying glacial isostatic effects on the Earth’s lithosphere, increasingly more studies focus on modeling the large-scale ice cover in mountainous regions over long time scales. However, balancing model complexity and the spatial extent with computational costs is challenging. Previous studies of large-scale ice cover simulation in mountain areas such as the European Alps, New Zealand, and the Tibetan Plateau, typically used 1-2 km spatial resolution. However, mountains are characterized by high peaks and steep slopes - topographic features that are crucial for glacier mass balance and dynamics, but poorly resolved in coarse resolution topography.

The Instructed Glacier Model (IGM) is a novel 3D ice model equipped with a Convolutional Neural Network which is trained from high order ice models to simulate ice flow. This results in a significant acceleration of run times, and thereby opening the possibility of running in higher spatial resolution. We use IGM to perform simulations of the entire European Alps (covering 480 240 km2), comparing models with 200 m and 2000 m resolution. We apply a linear cooling rate to today’s climate until 6 °C colder to mimic ice age conditions and model the expanding ice cover over a time period of 70,000 years.

Preliminary results indicate systematic, resolution-related differences: At the beginning of cooling, when ice accumulates at high elevations, the lower resolution yields slightly more ice volume. However, this trend reverses after ~ 41,000 years, right before the large valleys are filled with thick ice. When the Alps are fully ice covered, we find up to 15% more ice volume in the higher resolution model. The differences in ice volume are not uniformly distributed in space. The higher resolution model yields thicker and more extensive ice in some regions - mostly large valley systems - of up to 12,000 km2 , and thinner and less extensive ice in other, slightly smaller regions of the Alps. Currently, we analyze to what extent the glacier flow from steep slopes into larger, shallow valleys is represented at the different resolutions.

How to cite: Werner, H., Scherler, D., Winkelmann, R., and Jouvet, G.: Impact of spatial resolution on large-scale ice cover modeling of mountainous regions, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-4200, https://doi.org/10.5194/egusphere-egu23-4200, 2023.

X5.217
|
EGU23-6272
|
ECS
|
Lena Nicola, Julius Garbe, Ronja Reese, and Ricarda Winkelmann

The Antarctic Ice Sheet is currently losing mass through ocean induced melting at the underside of its large ice shelves. In the future, ice shelf cavities could switch from a ‘’cold” to a “warm” state, following a distinct increase in ocean temperatures e.g. by a redirection of coastal currents allowing warm circumpolar waters to access Antarctic grounding lines. With the ice-sheet model PISM, we delineate potential thresholds, at which the Antarctic Ice Sheet could experience ocean-induced non-linear ice loss. To this end, we apply circum-Antarctic ocean temperature perturbations of 1 to 5 K for different durations, ranging from tens to hundreds of years, and analyze the ice-sheet evolution after reversing the forcing over centennial to millennial timescales. Additionally, we perform ice-sheet simulations in which we slowly ramp up our forcing over similar timescales. Using these idealized overshoot scenarios, we analyze when and where critical thresholds that lead to large-scale, irreversible grounding line retreat are crossed.  We assess uncertainties of these thresholds by analysing the initial state uncertainty as well as parametric and structural uncertainties.

How to cite: Nicola, L., Garbe, J., Reese, R., and Winkelmann, R.: Identifying thresholds of ocean-induced Antarctic ice loss through idealized ice-sheet model simulations, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6272, https://doi.org/10.5194/egusphere-egu23-6272, 2023.

X5.218
|
EGU23-6509
|
ECS
|
Gustav Jungdal-Olesen, Vivi Kathrine Pedersen, Jane Lund Andersen, and Andreas Born

Ice sheets have shaped the Scandinavian landscape during numerous glacial periods throughout the Quaternary, but little is known about the effects of a changing landscape on the Scandinavian ice sheets in return. Here, we use a higher-order ice-sheet model (iSOSIA) to investigate how past morphological changes in the Scandinavian landscape may have affected ice-sheet extent and dynamics. Our preliminary results indicate that the Scandinavian ice sheet would have extended further south before the formation of the Norwegian Channel, which is believed to have been formed by glacial erosion during recent glacial periods (since ∼0.5 Ma). This suggests that landscape changes should be considered in addition to varying climate conditions, when exploring changes in ice-sheet dynamics and extent between glacial periods.

How to cite: Jungdal-Olesen, G., Pedersen, V. K., Andersen, J. L., and Born, A.: The influence of landscape evolution on Scandinavian Ice Sheet dynamics and extent, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-6509, https://doi.org/10.5194/egusphere-egu23-6509, 2023.

X5.219
|
EGU23-8932
|
ECS
Thomas Gregov, Frank Pattyn, and Maarten Arnst

Marine ice sheets are complex systems whose response to external forcing is the subject of much attention in the scientific community. In particular, the West Antarctic ice sheet, which could have a significant impact on future sea-level rise, is a major concern. One method of studying the response of marine ice sheets consists in investigating the relationship between the parameters and the equilibrium states of such systems. However, this is typically done by varying these parameters and letting the ice sheets evolve to new steady states, i.e., through transient simulations, which are computationally expensive.

An alternative is to consider continuation methods, where the equilibrium state of a system is studied directly as a function of the parameters. Such an approach has already been used in glaciology to study the mechanical behaviour of 1D marine ice sheets (Mulder et al., 2018), highlighting the hysteresis phenomena that had previously been obtained theoretically (Schoof, 2007). However, this study, because it only considers 1D geometries, does not allow to take into account the effect of lateral drag and of complex bedrock geometries, which are two factors that have the potential to stabilize the grounding line (Gudmundsson, 2013, Sergienko and Wingham, 2021).

Here we consider the continuation problem in the context of 2D marine ice sheets. This introduces several mathematical difficulties, notably related to the treatment of the distinction between grounded and floating parts. Mathematically, the problem takes the form of a contact problem between the bedrock and the lower part of the ice sheet. This leads to a system of equations that is not differentiable, which is challenging to solve numerically. We address these challenges in the context of the continuation problem, and propose several solutions, including a norm-based approach that is inspired from earlier studies (Mittelmann, 1987). Finally, we present some preliminary results which show that our numerical method is promising.

How to cite: Gregov, T., Pattyn, F., and Arnst, M.: Investigation of numerical continuation methods for marine ice-sheet systems formulated as contact problems, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-8932, https://doi.org/10.5194/egusphere-egu23-8932, 2023.

X5.220
|
EGU23-9343
|
ECS
Andrés Castillo-Llarena, Franco Retamal-Ramirez, Jorge Bernales, Martin Jacques-Coper, and Irina Rogozhina

During the Last Glacial Maximum (LGM, 23 to 19 thousand years ago), the Patagonian Ice Sheet (PIS) expanded along the Andes between ~ 38 °S to 55 °S. Existing paleoclimatic and paleoglacial evidence, especially that derived from glacial landforms, seems to indicate that the timing of maximum ice sheet expansions in the Southern and Northern Hemispheres was not synchronous. Moreover, significant uncertainties still exist in the onset of regional deglaciation and its major drivers. In this context, a combination of ice sheet modelling, glacial geochronology and paleoclimate reconstructions can provide important insights into the former PIS geometry and its contribution to the sea level low during the LGM. It can also help us infer likely paleoclimate scenarios and climate models that capture regional climate responses to global change in the most realistic manner.

Here we present an ensemble of numerical ice sheet simulations of the PIS at the LGM to constrain an envelope of probable atmospheric conditions derived from a range of model-based climate forcing products from the phases 3 and 4 of the Paleoclimate Modelling Intercomparison Project (PMIP). The resulting ensemble is then used as a guideline to identify sectors of the PIS where a significant disagreement between the field evidence and modelling results is obtained, highlighting a strong dependence of the PIS geometry on the uncertainties in near-surface air temperature forcing. We find that all ensemble members consistently fail to reproduce the ice sheet extent towards the northern part of Patagonia within the explored model parameter space. At the same time, the modelled PIS expands beyond its southeastern reconstructed boundary. Our analysis of the ice sheet’s mass budget seems to indicate that these discrepancies between the modelled and reconstructed PIS extents arise from poorly resolved topographic features within the global climate models and the general lack of observational data on ice thickness distribution during the LGM. We conclude that INM-CM4-8 and MPI-ESM1-2-LR produce the most realistic climate forcing across Patagonia at the LGM. It should be kept in mind that this analysis is based only on the evaluation of modelled ice sheet extents against geological evidence, as observational data on the former ice sheet thickness are still lacking. Nevertheless, our analysis suggests that the quality of the regional model-based climate reconstructions is directly linked to the horizontal resolution that must be capable of resolving topographic features of the Andes and ideally of the PIS itself.

How to cite: Castillo-Llarena, A., Retamal-Ramirez, F., Bernales, J., Jacques-Coper, M., and Rogozhina, I.: Patagonian Ice Sheet extent as an indicator of the regional climate regime at the Last Glacial Maximum, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9343, https://doi.org/10.5194/egusphere-egu23-9343, 2023.

X5.221
|
EGU23-9497
|
ECS
|
Koi McArthur, Felicity S. McCormack, and Christine Dow

The key process of basal sliding in Antarctic glaciers is often incorporated into ice dynamics models via the use of a friction law, which relates the basal shear stress to the effective pressure. With few ice dynamics models actively coupled to subglacial hydrology models, the effects of subglacial hydrology often manifest in the friction coefficient – an unknown parameter in the friction law. We investigate the impact of friction coefficients for Denman Glacier, East Antarctica, by comparing Ice-sheet and Sea-level System Model (ISSM) inversion simulations using the effective pressure produced from the Glacier Drainage System (GlaDS) model compared with  a typically prescribed effective pressure using a combination of ice overburden pressure and height above sea level (NO). We apply these comparative model runs for the Budd and Schoof friction laws. In regions of fast ice flow, we find a positive correlation between the GlaDS output effective pressure and the friction coefficient for the Schoof law. In addition, using the GlaDS output effective pressure compared to  NO leads to a smoother friction coefficient as well as smaller differences between the simulated and observed surface velocity. In general we find that spatial variations in the Schoof law match more closely with the known physics of subglacial hydrology than the Budd law and therefore suggest that using the GlaDS output effective pressure compared to NO produces more realistic results. This demonstrates the need to couple ice sheet and subglacial hydrological systems to accurately represent ice flow.

How to cite: McArthur, K., McCormack, F. S., and Dow, C.: Basal conditions of Denman Glacier from hydrology modeling and their application to various friction laws, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-9497, https://doi.org/10.5194/egusphere-egu23-9497, 2023.

X5.222
|
EGU23-10914
Su-jeong Lim and Byung-Dal So

Marine ice sheet dynamics plays a key role in the mechanical interaction between glaciers and the ocean. A steady-state marine ice sheet, due to the buttressing force of ice shelves, suppresses the glaciers from discharging into the ocean. Various atmospheric-oceanic-glacial mechanisms such as ice melting, buttressing, and snow accumulation lead to the instability of marine ice sheets. The unstable marine ice sheet cause thinning ice thickness and the resultant ice shelf collapse. Ice shelf collapse accelerates the rapid glaciers discharge and retreats, directly affecting sea level rise. The positive feedback loop of these processes significantly influences the ice sheets dynamics, especially the advance and retreat of the grounding line. The grounding line can be defined as a transition zone between the grounded ice sheet in bedrock and the floating ice shelf due to water buoyancy. Understanding grounding line dynamics may help us predict marine ice sheet instability and future sea level rise. It has been challenging to determine an accurate grounding line in the numerical simulation of the marine ice sheet due to mesh resolution. Here we performed two-dimensional ice flow modeling using open-source finite element software (i.e., Elmer/ice), to quantitatively evaluate the effect of diverse mesh sizes on the position of the grounding line. We tested a series of numerical models to precisely define varying the grounding line position influenced by bedrock topography, snow accumulation, and ice melting. The models with high mesh resolution required short time steps to obtain the accurate grounding line position. To consider both calculation efficiency and the accuracy of the position of the grounding line, we found adequate time steps corresponding to each mesh size ambient the grounding line. Our systematic results can provide inspiration for choosing a suitable mesh size and time step to determine a more accurate grounding line position.

How to cite: Lim, S. and So, B.-D.: Effect of mesh resolution on accurate grounding line definition using 2D finite element software, Elmer/ice, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10914, https://doi.org/10.5194/egusphere-egu23-10914, 2023.

X5.223
|
EGU23-10924
Yi-Jeong Baek and Byung-Dal So

Glacier mass loss related to ice flow is a quantitative factor that controls sea level rise and global warming. Ice flow is attributed to the complicated interaction between grounding line migration, calving basal friction, and topography. Previous studies using massive geophysical observation, including IPR (Ice Penetrating Radar), InSAR (Interferometric Synthetic Aperture Rada), and heat flow exploration, highlight the need to develop a sophisticated numerical model. Because it is difficult to numerically solve the full-stokes equation obtaining ice velocities, the simplified governing equation (i.e., HO, Higher Order model; SSA, Shallow Shelf Approximation; SIA, Shallow Ice Approximation) is widely used in the ice sheet dynamics community. The SSA approach, which assumes a depth-independent velocity model, has computational cost reduction based on simplified descriptions of the full-stokes equation. Here we developed the two-dimensional SSA numerical model to better understand ice dynamics using COMSOL Multiphysics® (hereafter COMSOL), which is a user-friendly finite element software providing convenient GUI, mesh generation, post-processing tools, various types of elements, and the order of shape function. To verify the application of COMSOL to ice sheet dynamics, we compared it with an open-source finite element package, Ice Sheet System Model (ISSM). The concise toy model successfully simulated the distribution of viscosity and velocity and the evolution of surface topography in both COMSOL and ISSM. Furthermore, we applied realistic bathymetry data (i.e., Bedmap2) of Pine Island Glacier, where the large ice mass loss occurs, to clarify more similarity of each model. The surface velocity patterns calculated by COMSOL and ISSM are significantly similar for various physical properties (e.g., ice viscosity, friction, and rate of accumulation and melt). We propose that COMSOL, which efficiently handles mesh generation and visualization, and various weak forms, can sufficiently be applied in ice dynamics.

How to cite: Baek, Y.-J. and So, B.-D.: Development of ice sheet model using COMSOL Multiphysics®: Comparison with Ice Sheet System Model (ISSM), EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-10924, https://doi.org/10.5194/egusphere-egu23-10924, 2023.

X5.224
|
EGU23-12477
Yoni Verhaegen and Philippe Huybrechts

We have modelled the influence of a supraglacial debris cover on the behavior of an idealized reference mountain glacier. A calibrated 3D coupled ice flow-mass balance-supraglacial debris cover model is used to assess the impact of the melt-altering effect of various supraglacial debris deposit rates on the overall steady state characteristics of the glacier. Additional experiments are also carried out to simulate the behavior of the debris-covered glacier in a warming future climate. The main results show that, when compared to its clean-ice version, the debris-covered version of the glacier exhibits longer but thinner ablation zones, accompanied by lower ice flow velocities, lower runoff production, as well as a dampening of the mass balance-elevation profile over the debris-covered ice. Experiments for warming climatic conditions primarily point out towards a significant delay of glacier retreat, as the dominant process for ice mass loss encompasses thinning out of the ablation zone rather than retreat. The above-mentioned effects are modelled to be increasingly pronounced with an increasing thickness and extent of the superimposed supraglacial debris cover.

How to cite: Verhaegen, Y. and Huybrechts, P.: Modelling the influence of a supraglacial debris cover on the mass balance and dynamics of mountain glaciers using a 3D higher order ice flow model, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-12477, https://doi.org/10.5194/egusphere-egu23-12477, 2023.

X5.225
|
EGU23-14854
|
ECS
Jowan Barnes and G. Hilmar Gudmundsson

Thwaites and Pine Island Glaciers in the Amundsen Sea Embayment, West Antarctica, are among the fastest evolving on the continent, and hold enough ice between them to raise sea levels by over a metre. In their present states, the two glaciers represent different configurations of floating ice, and therefore may not respond in the same way to changes in ocean forcing. Pine Island Ice Shelf is contained within a bay and provides a large amount of buttressing to its glacier. Thwaites Ice Shelf has two components; a heavily damaged ice tongue and a shelf which is only restrained by a single pinning point. Neither of these provide much buttressing. In our modelling experiments, we prescribe calving rates to the ice shelves of these two glaciers alongside thermal forcing from ISMIP6, to investigate the combined effects of warming oceans and continued calving on the future of the region. We demonstrate the potential impacts of adding the calving process into our model by using a range of constant calving rates. Examining the different responses can tell us how important the process is for each glacier, and how sensitive they are to changes in calving. This in turn can be used to determine whether or not significant effort should be invested in improving calving laws to more accurately predict the future shape of the ice front. 

How to cite: Barnes, J. and Gudmundsson, G. H.: How calving could affect the future of Thwaites and Pine Island Glaciers, EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-14854, https://doi.org/10.5194/egusphere-egu23-14854, 2023.

X5.226
|
EGU23-2365
Suzanne Bevan, Stephen Cornford, Inés Otosaka,, Trystan Surawy-Stepney, Lin Gilbert, and Daniel Martin

Mass loss from the Amundsen Sea Embayment of the West Antarctic Ice Sheet has been increasing over recent decades and is a major contributor to global sea level rise. Predictions of future sea level rise are increasingly modelled using ensembles of simulations within which model parameters and external forcings are varied widely then scored according to observations. Accurately reporting the uncertainty associated with these predictions is vital to enable effective planning for, and maybe construction of defences against, rising sea levels. Here we constrain, or calibrate, an ensemble of simulations of ice loss from the Amundsen Sea Embayment using the BISICLES ice sheet model with remotely sensed observations of surface elevation change and ice speed. The calibrations make it possible to reduce the 90% credibility bounds of predicted contributions to sea-level rise by 40%.

How to cite: Bevan, S., Cornford, S., Otosaka,, I., Surawy-Stepney, T., Gilbert, L., and Martin, D.: Bayesian calibration of an ice sheet model for the Amundsen Sea Embayment region., EGU General Assembly 2023, Vienna, Austria, 23–28 Apr 2023, EGU23-2365, https://doi.org/10.5194/egusphere-egu23-2365, 2023.