Articles | Volume 18, issue 9
https://doi.org/10.5194/tc-18-3991-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/tc-18-3991-2024
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Ice viscosity governs hydraulic fracture that causes rapid drainage of supraglacial lakes
Tim Hageman
Department of Engineering Science, University of Oxford, Oxford OX1 3PJ, UK
Jessica Mejía
Department of Geology, the University at Buffalo, Buffalo, NY 14260, USA
Ravindra Duddu
CORRESPONDING AUTHOR
Department of Civil and Environmental Engineering, Vanderbilt university, Nashville, TN37235 USA
Department of Earth and Environmental sciences, Vanderbilt university, Nashville, TN37235 USA
Emilio Martínez-Pañeda
CORRESPONDING AUTHOR
Department of Engineering Science, University of Oxford, Oxford OX1 3PJ, UK
Department of Civil and Environmental Engineering, Imperial College London, London SW7 2AZ, UK
Related authors
Theo Clayton, Ravindra Duddu, Tim Hageman, and Emilio Martínez-Pañeda
The Cryosphere, 18, 5573–5593, https://doi.org/10.5194/tc-18-5573-2024, https://doi.org/10.5194/tc-18-5573-2024, 2024
Short summary
Short summary
We develop and validate new analytical solutions that quantitatively consider how the properties of ice vary along the depth of ice shelves and that can be readily used in existing ice sheet models. Depth-varying firn properties are found to have a profound impact on ice sheet fracture and calving events. Our results show that grounded glaciers are less vulnerable than previously anticipated, while floating ice shelves are significantly more vulnerable to fracture and calving.
Theo Clayton, Ravindra Duddu, Tim Hageman, and Emilio Martínez-Pañeda
The Cryosphere, 18, 5573–5593, https://doi.org/10.5194/tc-18-5573-2024, https://doi.org/10.5194/tc-18-5573-2024, 2024
Short summary
Short summary
We develop and validate new analytical solutions that quantitatively consider how the properties of ice vary along the depth of ice shelves and that can be readily used in existing ice sheet models. Depth-varying firn properties are found to have a profound impact on ice sheet fracture and calving events. Our results show that grounded glaciers are less vulnerable than previously anticipated, while floating ice shelves are significantly more vulnerable to fracture and calving.
Meghana Ranganathan, Alexander A. Robel, Alexander Huth, and Ravindra Duddu
EGUsphere, https://doi.org/10.5194/egusphere-2024-1850, https://doi.org/10.5194/egusphere-2024-1850, 2024
Short summary
Short summary
The rate of ice loss from ice sheets is controlled by the flow of ice from the center of the ice sheet and internal fracturing of the ice. These are coupled – fractures can reduce the viscosity of ice and enable more rapid flow, and rapid flow can cause fracturing of ice. We present a simplified way of representing damage that is applicable to long-timescale flow estimates. Using this model, we find that including fracturing into an ice sheet simulation can increase the loss of ice by 13–29 %.
Celia Trunz, Kristin Poinar, Lauren C. Andrews, Matthew D. Covington, Jessica Mejia, Jason Gulley, and Victoria Siegel
The Cryosphere, 17, 5075–5094, https://doi.org/10.5194/tc-17-5075-2023, https://doi.org/10.5194/tc-17-5075-2023, 2023
Short summary
Short summary
Models simulating water pressure variations at the bottom of glaciers must use large storage parameters to produce realistic results. Whether that storage occurs englacially (in moulins) or subglacially is a matter of debate. Here, we directly simulate moulin volume to constrain the storage there. We find it is not enough. Instead, subglacial processes, including basal melt and input from upstream moulins, must be responsible for stabilizing these water pressure fluctuations.
Related subject area
Discipline: Ice sheets | Subject: Numerical Modelling
Two-way coupling between ice flow and channelized subglacial drainage enhances modeled marine-ice-sheet retreat
Sensitivity of the future evolution of the Wilkes Subglacial Basin ice sheet to grounding-line melt parameterizations
Antarctic sensitivity to oceanic melting parameterizations
Analytical solutions for the advective–diffusive ice column in the presence of strain heating
Biases in ice sheet models from missing noise-induced drift
Present-day mass loss rates are a precursor for West Antarctic Ice Sheet collapse
Modeling the timing of Patagonian Ice Sheet retreat in the Chilean Lake District from 22–10 ka
Using specularity content to evaluate eight geothermal heat flow maps of Totten Glacier
Surging of a Hudson Strait-scale ice stream: subglacial hydrology matters but the process details mostly do not
Regularization and L-curves in ice sheet inverse models: a case study in the Filchner–Ronne catchment
Quantifying the uncertainty in the Eurasian ice-sheet geometry at the Penultimate Glacial Maximum (Marine Isotope Stage 6)
The stability of present-day Antarctic grounding lines – Part 2: Onset of irreversible retreat of Amundsen Sea glaciers under current climate on centennial timescales cannot be excluded
The stability of present-day Antarctic grounding lines – Part 1: No indication of marine ice sheet instability in the current geometry
Geothermal heat flux is the dominant source of uncertainty in englacial-temperature-based dating of ice rise formation
Improving interpretation of sea-level projections through a machine-learning-based local explanation approach
Subglacial hydrology modulates basal sliding response of the Antarctic ice sheet to climate forcing
The predictive power of ice sheet models and the regional sensitivity of ice loss to basal sliding parameterisations: a case study of Pine Island and Thwaites glaciers, West Antarctica
Simulations of firn processes over the Greenland and Antarctic ice sheets: 1980–2021
Evaluation of six geothermal heat flux maps for the Antarctic Lambert–Amery glacial system
Impact of runoff temporal distribution on ice dynamics
Can changes in deformation regimes be inferred from crystallographic preferred orientations in polar ice?
Stabilizing effect of mélange buttressing on the marine ice-cliff instability of the West Antarctic Ice Sheet
Effective coefficient of diffusion and permeability of firn at Dome C and Lock In, Antarctica, and of various snow types – estimates over the 100–850 kg m−3 density range
The instantaneous impact of calving and thinning on the Larsen C Ice Shelf
Derivation of bedrock topography measurement requirements for the reduction of uncertainty in ice-sheet model projections of Thwaites Glacier
A comparison of the stability and performance of depth-integrated ice-dynamics solvers
A new vertically integrated MOno-Layer Higher-Order (MOLHO) ice flow model
On the contribution of grain boundary sliding type creep to firn densification – an assessment using an optimization approach
Marine ice sheet experiments with the Community Ice Sheet Model
The transferability of adjoint inversion products between different ice flow models
Inferring the basal sliding coefficient field for the Stokes ice sheet model under rheological uncertainty
The tipping points and early warning indicators for Pine Island Glacier, West Antarctica
Sensitivity of ice sheet surface velocity and elevation to variations in basal friction and topography in the full Stokes and shallow-shelf approximation frameworks using adjoint equations
Quantifying the effect of ocean bed properties on ice sheet geometry over 40 000 years with a full-Stokes model
Bayesian calibration of firn densification models
A kinematic formalism for tracking ice–ocean mass exchange on the Earth's surface and estimating sea-level change
Results of the third Marine Ice Sheet Model Intercomparison Project (MISMIP+)
Ocean-forced evolution of the Amundsen Sea catchment, West Antarctica, by 2100
Parameter sensitivity analysis of dynamic ice sheet models – numerical computations
Simulated retreat of Jakobshavn Isbræ during the 21st century
Development of physically based liquid water schemes for Greenland firn-densification models
Regional grid refinement in an Earth system model: impacts on the simulated Greenland surface mass balance
initMIP-Antarctica: an ice sheet model initialization experiment of ISMIP6
Modeling the response of northwest Greenland to enhanced ocean thermal forcing and subglacial discharge
Assessment of the Greenland ice sheet–atmosphere feedbacks for the next century with a regional atmospheric model coupled to an ice sheet model
Sensitivity of centennial mass loss projections of the Amundsen basin to the friction law
Retreat of Thwaites Glacier, West Antarctica, over the next 100 years using various ice flow models, ice shelf melt scenarios and basal friction laws
Comparison of four calving laws to model Greenland outlet glaciers
Neutral equilibrium and forcing feedbacks in marine ice sheet modelling
Representation of basal melting at the grounding line in ice flow models
George Lu and Jonathan Kingslake
The Cryosphere, 18, 5301–5321, https://doi.org/10.5194/tc-18-5301-2024, https://doi.org/10.5194/tc-18-5301-2024, 2024
Short summary
Short summary
Water below ice sheets affects ice-sheet motion, while the evolution of ice sheets likewise affects the water below. We create a model that allows for water and ice to affect each other and use it to see how this coupling or lack thereof may impact ice-sheet retreat. We find that coupling an evolving water system with the ice sheet results in more retreat than if we assume unchanging conditions under the ice, which indicates a need to better represent the effects of water in ice-sheet models.
Yu Wang, Chen Zhao, Rupert Gladstone, Thomas Zwinger, Benjamin K. Galton-Fenzi, and Poul Christoffersen
The Cryosphere, 18, 5117–5137, https://doi.org/10.5194/tc-18-5117-2024, https://doi.org/10.5194/tc-18-5117-2024, 2024
Short summary
Short summary
Our research delves into the future evolution of Antarctica's Wilkes Subglacial Basin (WSB) and its potential contribution to sea level rise, focusing on how basal melt is implemented at the grounding line in ice flow models. Our findings suggest that these implementation methods can significantly impact the magnitude of future ice loss projections. Under a high-emission scenario, the WSB ice sheet could undergo massive and rapid retreat between 2200 and 2300.
Antonio Juarez-Martinez, Javier Blasco, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4257–4283, https://doi.org/10.5194/tc-18-4257-2024, https://doi.org/10.5194/tc-18-4257-2024, 2024
Short summary
Short summary
We present sea level projections for Antarctica in the context of ISMIP6-2300 with several forcings but extend the simulations to 2500, showing that more than 3 m of sea level contribution could be reached. We also test the sensitivity on a basal melting parameter and determine the timing of the loss of ice in the west region. All the simulations were carried out with the ice sheet model Yelmo.
Daniel Moreno-Parada, Alexander Robinson, Marisa Montoya, and Jorge Alvarez-Solas
The Cryosphere, 18, 4215–4232, https://doi.org/10.5194/tc-18-4215-2024, https://doi.org/10.5194/tc-18-4215-2024, 2024
Short summary
Short summary
Our study tries to understand how the ice temperature evolves in a large mass as in the case of Antarctica. We found a relation that tells us the ice temperature at any point. These results are important because they also determine how the ice moves. In general, ice moves due to slow deformation (as if pouring honey from a jar). Nevertheless, in some regions the ice base warms enough and melts. The liquid water then serves as lubricant and the ice slides and its velocity increases rapidly.
Alexander A. Robel, Vincent Verjans, and Aminat A. Ambelorun
The Cryosphere, 18, 2613–2623, https://doi.org/10.5194/tc-18-2613-2024, https://doi.org/10.5194/tc-18-2613-2024, 2024
Short summary
Short summary
The average size of many glaciers and ice sheets changes when noise is added to the system. The reasons for this drift in glacier state is intrinsic to the dynamics of how ice flows and the bumpiness of the Earth's surface. We argue that not including noise in projections of ice sheet evolution over coming decades and centuries is a pervasive source of bias in these computer models, and so realistic variability in glacier and climate processes must be included in models.
Tim van den Akker, William H. Lipscomb, Gunter R. Leguy, Jorjo Bernales, Constantijn Berends, Willem Jan van de Berg, and Roderik S. W. van de Wal
EGUsphere, https://doi.org/10.5194/egusphere-2024-851, https://doi.org/10.5194/egusphere-2024-851, 2024
Short summary
Short summary
In this study, we present an improved way of representing ice thickness change rates into an ice sheet model. We apply this method using two ice sheet models on the Antarctic Ice Sheet. We found that the two largest outlet glaciers on the Antarctic Ice Sheet, the Thwaites Glacier and Pine Island Glacier, will collapse without further warming on a timescale of centuries. This would cause a sea level rise of about 1.2 meters globally.
Joshua Cuzzone, Matias Romero, and Shaun A. Marcott
The Cryosphere, 18, 1381–1398, https://doi.org/10.5194/tc-18-1381-2024, https://doi.org/10.5194/tc-18-1381-2024, 2024
Short summary
Short summary
We simulate the retreat history of the Patagonian Ice Sheet (PIS) across the Chilean Lake District from 22–10 ka. These results improve our understanding of the response of the PIS to deglacial warming and the patterns of deglacial ice margin retreat where gaps in the geologic record still exist, and they indicate that changes in large-scale precipitation during the last deglaciation played an important role in modulating the response of ice margin change across the PIS to deglacial warming.
Yan Huang, Liyun Zhao, Michael Wolovick, Yiliang Ma, and John C. Moore
The Cryosphere, 18, 103–119, https://doi.org/10.5194/tc-18-103-2024, https://doi.org/10.5194/tc-18-103-2024, 2024
Short summary
Short summary
Geothermal heat flux (GHF) is an important factor affecting the basal thermal environment of an ice sheet and crucial for its dynamics. But it is poorly defined for the Antarctic ice sheet. We simulate the basal temperature and basal melting rate with eight different GHF datasets. We use specularity content as a two-sided constraint to discriminate between local wet or dry basal conditions. Two medium-magnitude GHF distribution maps rank well, showing that most of the inland bed area is frozen.
Matthew Drew and Lev Tarasov
The Cryosphere, 17, 5391–5415, https://doi.org/10.5194/tc-17-5391-2023, https://doi.org/10.5194/tc-17-5391-2023, 2023
Short summary
Short summary
The interaction of fast-flowing regions of continental ice sheets with their beds governs how quickly they slide and therefore flow. The coupling of fast ice to its bed is controlled by the pressure of meltwater at its base. It is currently poorly understood how the physical details of these hydrologic systems affect ice speedup. Using numerical models we find, surprisingly, that they largely do not, except for the duration of the surge. This suggests that cheap models are sufficient.
Michael Wolovick, Angelika Humbert, Thomas Kleiner, and Martin Rückamp
The Cryosphere, 17, 5027–5060, https://doi.org/10.5194/tc-17-5027-2023, https://doi.org/10.5194/tc-17-5027-2023, 2023
Short summary
Short summary
The friction underneath ice sheets can be inferred from observed velocity at the top, but this inference requires smoothing. The selection of smoothing has been highly variable in the literature. Here we show how to rigorously select the best smoothing, and we show that the inferred friction converges towards the best knowable field as model resolution improves. We use this to learn about the best description of basal friction and to formulate recommended best practices for other modelers.
Oliver G. Pollard, Natasha L. M. Barlow, Lauren J. Gregoire, Natalya Gomez, Víctor Cartelle, Jeremy C. Ely, and Lachlan C. Astfalck
The Cryosphere, 17, 4751–4777, https://doi.org/10.5194/tc-17-4751-2023, https://doi.org/10.5194/tc-17-4751-2023, 2023
Short summary
Short summary
We use advanced statistical techniques and a simple ice-sheet model to produce an ensemble of plausible 3D shapes of the ice sheet that once stretched across northern Europe during the previous glacial maximum (140,000 years ago). This new reconstruction, equivalent in volume to 48 ± 8 m of global mean sea-level rise, will improve the interpretation of high sea levels recorded from the Last Interglacial period (120 000 years ago) that provide a useful perspective on the future.
Ronja Reese, Julius Garbe, Emily A. Hill, Benoît Urruty, Kaitlin A. Naughten, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, David Chandler, Petra M. Langebroek, and Ricarda Winkelmann
The Cryosphere, 17, 3761–3783, https://doi.org/10.5194/tc-17-3761-2023, https://doi.org/10.5194/tc-17-3761-2023, 2023
Short summary
Short summary
We use an ice sheet model to test where current climate conditions in Antarctica might lead. We find that present-day ocean and atmosphere conditions might commit an irreversible collapse of parts of West Antarctica which evolves over centuries to millennia. Importantly, this collapse is not irreversible yet.
Emily A. Hill, Benoît Urruty, Ronja Reese, Julius Garbe, Olivier Gagliardini, Gaël Durand, Fabien Gillet-Chaulet, G. Hilmar Gudmundsson, Ricarda Winkelmann, Mondher Chekki, David Chandler, and Petra M. Langebroek
The Cryosphere, 17, 3739–3759, https://doi.org/10.5194/tc-17-3739-2023, https://doi.org/10.5194/tc-17-3739-2023, 2023
Short summary
Short summary
The grounding lines of the Antarctic Ice Sheet could enter phases of irreversible retreat or advance. We use three ice sheet models to show that the present-day locations of Antarctic grounding lines are reversible with respect to a small perturbation away from their current position. This indicates that present-day retreat of the grounding lines is not yet irreversible or self-enhancing.
Aleksandr Montelli and Jonathan Kingslake
The Cryosphere, 17, 195–210, https://doi.org/10.5194/tc-17-195-2023, https://doi.org/10.5194/tc-17-195-2023, 2023
Short summary
Short summary
Thermal modelling and Bayesian inversion techniques are used to evaluate the uncertainties inherent in inferences of ice-sheet evolution from borehole temperature measurements. We show that the same temperature profiles may result from a range of parameters, of which geothermal heat flux through underlying bedrock plays a key role. Careful model parameterisation and evaluation of heat flux are essential for inferring past ice-sheet evolution from englacial borehole thermometry.
Jeremy Rohmer, Remi Thieblemont, Goneri Le Cozannet, Heiko Goelzer, and Gael Durand
The Cryosphere, 16, 4637–4657, https://doi.org/10.5194/tc-16-4637-2022, https://doi.org/10.5194/tc-16-4637-2022, 2022
Short summary
Short summary
To improve the interpretability of process-based projections of the sea-level contribution from land ice components, we apply the machine-learning-based
SHapley Additive exPlanationsapproach to a subset of a multi-model ensemble study for the Greenland ice sheet. This allows us to quantify the influence of particular modelling decisions (related to numerical implementation, initial conditions, or parametrisation of ice-sheet processes) directly in terms of sea-level change contribution.
Elise Kazmierczak, Sainan Sun, Violaine Coulon, and Frank Pattyn
The Cryosphere, 16, 4537–4552, https://doi.org/10.5194/tc-16-4537-2022, https://doi.org/10.5194/tc-16-4537-2022, 2022
Short summary
Short summary
The water at the interface between ice sheets and underlying bedrock leads to lubrication between the ice and the bed. Due to a lack of direct observations, subglacial conditions beneath the Antarctic ice sheet are poorly understood. Here, we compare different approaches in which the subglacial water could influence sliding on the underlying bedrock and suggest that it modulates the Antarctic ice sheet response and increases uncertainties, especially in the context of global warming.
Jowan M. Barnes and G. Hilmar Gudmundsson
The Cryosphere, 16, 4291–4304, https://doi.org/10.5194/tc-16-4291-2022, https://doi.org/10.5194/tc-16-4291-2022, 2022
Short summary
Short summary
Models must represent how glaciers slide along the bed, but there are many ways to do so. In this paper, several sliding laws are tested and found to affect different regions of the Antarctic Ice Sheet in different ways and at different speeds. However, the variability in ice volume loss due to sliding-law choices is low compared to other factors, so limited empirical knowledge of sliding does not prevent us from making predictions of how an ice sheet will evolve.
Brooke Medley, Thomas A. Neumann, H. Jay Zwally, Benjamin E. Smith, and C. Max Stevens
The Cryosphere, 16, 3971–4011, https://doi.org/10.5194/tc-16-3971-2022, https://doi.org/10.5194/tc-16-3971-2022, 2022
Short summary
Short summary
Satellite altimeters measure the height or volume change over Earth's ice sheets, but in order to understand how that change translates into ice mass, we must account for various processes at the surface. Specifically, snowfall events generate large, transient increases in surface height, yet snow fall has a relatively low density, which means much of that height change is composed of air. This air signal must be removed from the observed height changes before we can assess ice mass change.
Haoran Kang, Liyun Zhao, Michael Wolovick, and John C. Moore
The Cryosphere, 16, 3619–3633, https://doi.org/10.5194/tc-16-3619-2022, https://doi.org/10.5194/tc-16-3619-2022, 2022
Short summary
Short summary
Basal thermal conditions are important to ice dynamics and sensitive to geothermal heat flux (GHF). We estimate basal thermal conditions of the Lambert–Amery Glacier system with six GHF maps. Recent GHFs inverted from aerial geomagnetic observations produce a larger warm-based area and match the observed subglacial lakes better than the other GHFs. The modelled basal melt rate is 10 to hundreds of millimetres per year in fast-flowing glaciers feeding the Amery Ice Shelf and smaller inland.
Basile de Fleurian, Richard Davy, and Petra M. Langebroek
The Cryosphere, 16, 2265–2283, https://doi.org/10.5194/tc-16-2265-2022, https://doi.org/10.5194/tc-16-2265-2022, 2022
Short summary
Short summary
As temperature increases, more snow and ice melt at the surface of ice sheets. Here we use an ice dynamics and subglacial hydrology model with simplified geometry and climate forcing to study the impact of variations in meltwater on ice dynamics. We focus on the variations in length and intensity of the melt season. Our results show that a longer melt season leads to faster glaciers, but a more intense melt season reduces glaciers' seasonal velocities, albeit leading to higher peak velocities.
Maria-Gema Llorens, Albert Griera, Paul D. Bons, Ilka Weikusat, David J. Prior, Enrique Gomez-Rivas, Tamara de Riese, Ivone Jimenez-Munt, Daniel García-Castellanos, and Ricardo A. Lebensohn
The Cryosphere, 16, 2009–2024, https://doi.org/10.5194/tc-16-2009-2022, https://doi.org/10.5194/tc-16-2009-2022, 2022
Short summary
Short summary
Polar ice is formed by ice crystals, which form fabrics that are utilised to interpret how ice sheets flow. It is unclear whether fabrics result from the current flow regime or if they are inherited. To understand the extent to which ice crystals can be reoriented when ice flow conditions change, we simulate and evaluate multi-stage ice flow scenarios according to natural cases. We find that second deformation regimes normally overprint inherited fabrics, with a range of transitional fabrics.
Tanja Schlemm, Johannes Feldmann, Ricarda Winkelmann, and Anders Levermann
The Cryosphere, 16, 1979–1996, https://doi.org/10.5194/tc-16-1979-2022, https://doi.org/10.5194/tc-16-1979-2022, 2022
Short summary
Short summary
Marine cliff instability, if it exists, could dominate Antarctica's contribution to future sea-level rise. It is likely to speed up with ice thickness and thus would accelerate in most parts of Antarctica. Here, we investigate a possible mechanism that might stop cliff instability through cloaking by ice mélange. It is only a first step, but it shows that embayment geometry is, in principle, able to stop marine cliff instability in most parts of West Antarctica.
Neige Calonne, Alexis Burr, Armelle Philip, Frédéric Flin, and Christian Geindreau
The Cryosphere, 16, 967–980, https://doi.org/10.5194/tc-16-967-2022, https://doi.org/10.5194/tc-16-967-2022, 2022
Short summary
Short summary
Modeling gas transport in ice sheets from surface to close-off is key to interpreting climate archives. Estimates of the diffusion coefficient and permeability of snow and firn are required but remain a large source of uncertainty. We present a new dataset of diffusion coefficients and permeability from 20 to 120 m depth at two Antarctic sites. We suggest predictive formulas to estimate both properties over the entire 100–850 kg m3 density range, i.e., anywhere within the ice sheet column.
Tom Mitcham, G. Hilmar Gudmundsson, and Jonathan L. Bamber
The Cryosphere, 16, 883–901, https://doi.org/10.5194/tc-16-883-2022, https://doi.org/10.5194/tc-16-883-2022, 2022
Short summary
Short summary
We modelled the response of the Larsen C Ice Shelf (LCIS) and its tributary glaciers to the calving of the A68 iceberg and validated our results with observations. We found that the impact was limited, confirming that mostly passive ice was calved. Through further calving experiments we quantified the total buttressing provided by the LCIS and found that over 80 % of the buttressing capacity is generated in the first 5 km of the ice shelf downstream of the grounding line.
Blake A. Castleman, Nicole-Jeanne Schlegel, Lambert Caron, Eric Larour, and Ala Khazendar
The Cryosphere, 16, 761–778, https://doi.org/10.5194/tc-16-761-2022, https://doi.org/10.5194/tc-16-761-2022, 2022
Short summary
Short summary
In the described study, we derive an uncertainty range for global mean sea level rise (SLR) contribution from Thwaites Glacier in a 200-year period under an extreme ocean warming scenario. We derive the spatial and vertical resolutions needed for bedrock data acquisition missions in order to limit global mean SLR contribution from Thwaites Glacier to ±2 cm in a 200-year period. We conduct sensitivity experiments in order to present the locations of critical regions in need of accurate mapping.
Alexander Robinson, Daniel Goldberg, and William H. Lipscomb
The Cryosphere, 16, 689–709, https://doi.org/10.5194/tc-16-689-2022, https://doi.org/10.5194/tc-16-689-2022, 2022
Short summary
Short summary
Here we investigate the numerical stability of several commonly used methods in order to determine which of them are capable of resolving the complex physics of the ice flow and are also computationally efficient. We find that the so-called DIVA solver outperforms the others. Its representation of the physics is consistent with more complex methods, while it remains computationally efficient at high resolution.
Thiago Dias dos Santos, Mathieu Morlighem, and Douglas Brinkerhoff
The Cryosphere, 16, 179–195, https://doi.org/10.5194/tc-16-179-2022, https://doi.org/10.5194/tc-16-179-2022, 2022
Short summary
Short summary
Projecting the future evolution of Greenland and Antarctica and their potential contribution to sea level rise often relies on computer simulations carried out by numerical ice sheet models. Here we present a new vertically integrated ice sheet model and assess its performance using different benchmarks. The new model shows results comparable to a three-dimensional model at relatively lower computational cost, suggesting that it is an excellent alternative for long-term simulations.
Timm Schultz, Ralf Müller, Dietmar Gross, and Angelika Humbert
The Cryosphere, 16, 143–158, https://doi.org/10.5194/tc-16-143-2022, https://doi.org/10.5194/tc-16-143-2022, 2022
Short summary
Short summary
Firn is the interstage product between snow and ice. Simulations describing the process of firn densification are used in the context of estimating mass changes of the ice sheets and past climate reconstructions. The first stage of firn densification takes place in the upper few meters of the firn column. We investigate how well a material law describing the process of grain boundary sliding works for the numerical simulation of firn densification in this stage.
Gunter R. Leguy, William H. Lipscomb, and Xylar S. Asay-Davis
The Cryosphere, 15, 3229–3253, https://doi.org/10.5194/tc-15-3229-2021, https://doi.org/10.5194/tc-15-3229-2021, 2021
Short summary
Short summary
We present numerical features of the Community Ice Sheet Model in representing ocean termini glaciers. Using idealized test cases, we show that applying melt in a partly grounded cell is beneficial, in contrast to recent studies. We confirm that parameterizing partly grounded cells yields accurate ice sheet representation at a grid resolution of ~2 km (arguably 4 km), allowing ice sheet simulations at a continental scale. The choice of basal friction law also influences the ice flow.
Jowan M. Barnes, Thiago Dias dos Santos, Daniel Goldberg, G. Hilmar Gudmundsson, Mathieu Morlighem, and Jan De Rydt
The Cryosphere, 15, 1975–2000, https://doi.org/10.5194/tc-15-1975-2021, https://doi.org/10.5194/tc-15-1975-2021, 2021
Short summary
Short summary
Some properties of ice flow models must be initialised using observed data before they can be used to produce reliable predictions of the future. Different models have different ways of doing this, and the process is generally seen as being specific to an individual model. We compare the methods used by three different models and show that they produce similar outputs. We also demonstrate that the outputs from one model can be used in other models without introducing large uncertainties.
Olalekan Babaniyi, Ruanui Nicholson, Umberto Villa, and Noémi Petra
The Cryosphere, 15, 1731–1750, https://doi.org/10.5194/tc-15-1731-2021, https://doi.org/10.5194/tc-15-1731-2021, 2021
Short summary
Short summary
We consider the problem of inferring unknown parameter fields under additional uncertainty for an ice sheet model from synthetic surface ice flow velocity measurements. Our results indicate that accounting for model uncertainty stemming from the presence of nuisance parameters is crucial. Namely our findings suggest that using nominal values for these parameters, as is often done in practice, without taking into account the resulting modeling error can lead to overconfident and biased results.
Sebastian H. R. Rosier, Ronja Reese, Jonathan F. Donges, Jan De Rydt, G. Hilmar Gudmundsson, and Ricarda Winkelmann
The Cryosphere, 15, 1501–1516, https://doi.org/10.5194/tc-15-1501-2021, https://doi.org/10.5194/tc-15-1501-2021, 2021
Short summary
Short summary
Pine Island Glacier has contributed more to sea-level rise over the past decades than any other glacier in Antarctica. Ice-flow modelling studies have shown that it can undergo periods of rapid mass loss, but no study has shown that these future changes could cross a tipping point and therefore be effectively irreversible. Here, we assess the stability of Pine Island Glacier, quantifying the changes in ocean temperatures required to cross future tipping points using statistical methods.
Gong Cheng, Nina Kirchner, and Per Lötstedt
The Cryosphere, 15, 715–742, https://doi.org/10.5194/tc-15-715-2021, https://doi.org/10.5194/tc-15-715-2021, 2021
Short summary
Short summary
We present an inverse modeling approach to improve the understanding of spatiotemporally variable processes at the inaccessible base of an ice sheet by determining the sensitivity of direct surface observations to perturbations of basal conditions. Time dependency is proved to be important in these types of problems. The effect of perturbations is analyzed based on analytical and numerical solutions.
Clemens Schannwell, Reinhard Drews, Todd A. Ehlers, Olaf Eisen, Christoph Mayer, Mika Malinen, Emma C. Smith, and Hannes Eisermann
The Cryosphere, 14, 3917–3934, https://doi.org/10.5194/tc-14-3917-2020, https://doi.org/10.5194/tc-14-3917-2020, 2020
Short summary
Short summary
To reduce uncertainties associated with sea level rise projections, an accurate representation of ice flow is paramount. Most ice sheet models rely on simplified versions of the underlying ice flow equations. Due to the high computational costs, ice sheet models based on the complete ice flow equations have been restricted to < 1000 years. Here, we present a new model setup that extends the applicability of such models by an order of magnitude, permitting simulations of 40 000 years.
Vincent Verjans, Amber A. Leeson, Christopher Nemeth, C. Max Stevens, Peter Kuipers Munneke, Brice Noël, and Jan Melchior van Wessem
The Cryosphere, 14, 3017–3032, https://doi.org/10.5194/tc-14-3017-2020, https://doi.org/10.5194/tc-14-3017-2020, 2020
Short summary
Short summary
Ice sheets are covered by a firn layer, which is the transition stage between fresh snow and ice. Accurate modelling of firn density properties is important in many glaciological aspects. Current models show disagreements, are mostly calibrated to match specific observations of firn density and lack thorough uncertainty analysis. We use a novel calibration method for firn models based on a Bayesian statistical framework, which results in improved model accuracy and in uncertainty evaluation.
Surendra Adhikari, Erik R. Ivins, Eric Larour, Lambert Caron, and Helene Seroussi
The Cryosphere, 14, 2819–2833, https://doi.org/10.5194/tc-14-2819-2020, https://doi.org/10.5194/tc-14-2819-2020, 2020
Short summary
Short summary
The mathematical formalism presented in this paper aims at simplifying computational strategies for tracking ice–ocean mass exchange in the Earth system. To this end, we define a set of generic, and quite simple, descriptions of evolving land, ocean and ice interfaces and present a unified method to compute the sea-level contribution of evolving ice sheets. The formalism can be applied to arbitrary geometries and at all timescales.
Stephen L. Cornford, Helene Seroussi, Xylar S. Asay-Davis, G. Hilmar Gudmundsson, Rob Arthern, Chris Borstad, Julia Christmann, Thiago Dias dos Santos, Johannes Feldmann, Daniel Goldberg, Matthew J. Hoffman, Angelika Humbert, Thomas Kleiner, Gunter Leguy, William H. Lipscomb, Nacho Merino, Gaël Durand, Mathieu Morlighem, David Pollard, Martin Rückamp, C. Rosie Williams, and Hongju Yu
The Cryosphere, 14, 2283–2301, https://doi.org/10.5194/tc-14-2283-2020, https://doi.org/10.5194/tc-14-2283-2020, 2020
Short summary
Short summary
We present the results of the third Marine Ice Sheet Intercomparison Project (MISMIP+). MISMIP+ is one in a series of exercises that test numerical models of ice sheet flow in simple situations. This particular exercise concentrates on the response of ice sheet models to the thinning of their floating ice shelves, which is of interest because numerical models are currently used to model the response to contemporary and near-future thinning in Antarctic ice shelves.
Alanna V. Alevropoulos-Borrill, Isabel J. Nias, Antony J. Payne, Nicholas R. Golledge, and Rory J. Bingham
The Cryosphere, 14, 1245–1258, https://doi.org/10.5194/tc-14-1245-2020, https://doi.org/10.5194/tc-14-1245-2020, 2020
Gong Cheng and Per Lötstedt
The Cryosphere, 14, 673–691, https://doi.org/10.5194/tc-14-673-2020, https://doi.org/10.5194/tc-14-673-2020, 2020
Short summary
Short summary
We present a time-dependent inverse method for ice sheet modeling. By investigating the sensitivity of the observations of the velocity and the height at the surface to the basal conditions of the ice, we show that if the basal parameters are time dependent, then time cannot be ignored in the inversion. By looking at the numerical features, we conclude that adding the height information of an ice sheet in the velocity inversion procedure could improve the robustness of the inference.
Xiaoran Guo, Liyun Zhao, Rupert M. Gladstone, Sainan Sun, and John C. Moore
The Cryosphere, 13, 3139–3153, https://doi.org/10.5194/tc-13-3139-2019, https://doi.org/10.5194/tc-13-3139-2019, 2019
Vincent Verjans, Amber A. Leeson, C. Max Stevens, Michael MacFerrin, Brice Noël, and Michiel R. van den Broeke
The Cryosphere, 13, 1819–1842, https://doi.org/10.5194/tc-13-1819-2019, https://doi.org/10.5194/tc-13-1819-2019, 2019
Short summary
Short summary
Firn models rely on empirical approaches for representing the percolation and refreezing of meltwater through the firn column. We develop liquid water schemes of different levels of complexity for firn models and compare their performances with respect to observations of density profiles from Greenland. Our results demonstrate that physically advanced water schemes do not lead to better agreement with density observations. Uncertainties in other processes contribute more to model discrepancy.
Leonardus van Kampenhout, Alan M. Rhoades, Adam R. Herrington, Colin M. Zarzycki, Jan T. M. Lenaerts, William J. Sacks, and Michiel R. van den Broeke
The Cryosphere, 13, 1547–1564, https://doi.org/10.5194/tc-13-1547-2019, https://doi.org/10.5194/tc-13-1547-2019, 2019
Short summary
Short summary
A new tool is evaluated in which the climate and surface mass balance (SMB) of the Greenland ice sheet are resolved at 55 and 28 km resolution, while the rest of the globe is modelled at ~110 km. The local refinement of resolution leads to improved accumulation (SMB > 0) compared to observations; however ablation (SMB < 0) is deteriorated in some regions. This is attributed to changes in cloud cover and a reduced effectiveness of a model-specific vertical downscaling technique.
Hélène Seroussi, Sophie Nowicki, Erika Simon, Ayako Abe-Ouchi, Torsten Albrecht, Julien Brondex, Stephen Cornford, Christophe Dumas, Fabien Gillet-Chaulet, Heiko Goelzer, Nicholas R. Golledge, Jonathan M. Gregory, Ralf Greve, Matthew J. Hoffman, Angelika Humbert, Philippe Huybrechts, Thomas Kleiner, Eric Larour, Gunter Leguy, William H. Lipscomb, Daniel Lowry, Matthias Mengel, Mathieu Morlighem, Frank Pattyn, Anthony J. Payne, David Pollard, Stephen F. Price, Aurélien Quiquet, Thomas J. Reerink, Ronja Reese, Christian B. Rodehacke, Nicole-Jeanne Schlegel, Andrew Shepherd, Sainan Sun, Johannes Sutter, Jonas Van Breedam, Roderik S. W. van de Wal, Ricarda Winkelmann, and Tong Zhang
The Cryosphere, 13, 1441–1471, https://doi.org/10.5194/tc-13-1441-2019, https://doi.org/10.5194/tc-13-1441-2019, 2019
Short summary
Short summary
We compare a wide range of Antarctic ice sheet simulations with varying initialization techniques and model parameters to understand the role they play on the projected evolution of this ice sheet under simple scenarios. Results are improved compared to previous assessments and show that continued improvements in the representation of the floating ice around Antarctica are critical to reduce the uncertainty in the future ice sheet contribution to sea level rise.
Mathieu Morlighem, Michael Wood, Hélène Seroussi, Youngmin Choi, and Eric Rignot
The Cryosphere, 13, 723–734, https://doi.org/10.5194/tc-13-723-2019, https://doi.org/10.5194/tc-13-723-2019, 2019
Short summary
Short summary
Many glaciers along the coast of Greenland have been retreating. It has been suggested that this retreat is triggered by the presence of warm water in the fjords, and surface melt at the top of the ice sheet is exacerbating this problem. Here, we quantify the vulnerability of northwestern Greenland to further warming using a numerical model. We find that in current conditions, this sector alone will contribute more than 1 cm to sea rise level by 2100, and up to 3 cm in the most extreme scenario.
Sébastien Le clec'h, Sylvie Charbit, Aurélien Quiquet, Xavier Fettweis, Christophe Dumas, Masa Kageyama, Coraline Wyard, and Catherine Ritz
The Cryosphere, 13, 373–395, https://doi.org/10.5194/tc-13-373-2019, https://doi.org/10.5194/tc-13-373-2019, 2019
Short summary
Short summary
Quantifying the future contribution of the Greenland ice sheet (GrIS) to sea-level rise in response to atmospheric changes is important but remains challenging. For the first time a full representation of the feedbacks between a GrIS model and a regional atmospheric model was implemented. The authors highlight the fundamental need for representing the GrIS topography change feedbacks with respect to the atmospheric component face to the strong impact on the projected sea-level rise.
Julien Brondex, Fabien Gillet-Chaulet, and Olivier Gagliardini
The Cryosphere, 13, 177–195, https://doi.org/10.5194/tc-13-177-2019, https://doi.org/10.5194/tc-13-177-2019, 2019
Short summary
Short summary
Here, we apply a synthetic perturbation to the most active drainage basin of Antarctica and show that centennial mass loss projections obtained through ice flow models depend strongly on the implemented friction law, i.e. the mathematical relationship between basal drag and sliding velocities. In particular, the commonly used Weertman law considerably underestimates the sea-level contribution of this basin in comparison to two water pressure-dependent laws which rely on stronger physical bases.
Hongju Yu, Eric Rignot, Helene Seroussi, and Mathieu Morlighem
The Cryosphere, 12, 3861–3876, https://doi.org/10.5194/tc-12-3861-2018, https://doi.org/10.5194/tc-12-3861-2018, 2018
Short summary
Short summary
Thwaites Glacier, West Antarctica, has experienced rapid grounding line retreat and mass loss in the past decades. In this study, we simulate the evolution of Thwaites Glacier over the next century using different model configurations. Overall, we estimate a 5 mm contribution to global sea level rise from Thwaites Glacier in the next 30 years. However, a 300 % uncertainty is found over the next 100 years, ranging from 14 to 42 mm, depending on the model setup.
Youngmin Choi, Mathieu Morlighem, Michael Wood, and Johannes H. Bondzio
The Cryosphere, 12, 3735–3746, https://doi.org/10.5194/tc-12-3735-2018, https://doi.org/10.5194/tc-12-3735-2018, 2018
Short summary
Short summary
Calving is an important mechanism that controls the dynamics of Greenland outlet glaciers. We test and compare four calving laws and assess which calving law has better predictive abilities. Overall, the calving law based on von Mises stress is more satisfactory than other laws, but new parameterizations should be derived to better capture the detailed processes involved in calving.
Rupert M. Gladstone, Yuwei Xia, and John Moore
The Cryosphere, 12, 3605–3615, https://doi.org/10.5194/tc-12-3605-2018, https://doi.org/10.5194/tc-12-3605-2018, 2018
Short summary
Short summary
Computer models for the simulation of marine ice sheets (ice sheets resting on bedrock below sea level) historically show poor numerical convergence for grounding line (the boundary between grounded and floating parts of the ice sheet) movement. We have further characterised the nature of the numerical problems leading to poor convergence and highlighted implications for the design of computer experiments that test grounding line movement.
Hélène Seroussi and Mathieu Morlighem
The Cryosphere, 12, 3085–3096, https://doi.org/10.5194/tc-12-3085-2018, https://doi.org/10.5194/tc-12-3085-2018, 2018
Cited articles
Andrews, L. C., Hoffman, M. J., Neumann, T. A., and Catania, G. A.: Seasonal Evolution of the Subglacial Hydrologic System Modified by Supraglacial Lake Drainage in Western Greenland, J. Geophys. Res.-Earth, 123, 1479–1496, https://doi.org/10.1029/2017JF004585, 2018. a, b
Andrews, L. C., Poinar, K., and Trunz, C.: Controls on Greenland moulin geometry and evolution from the Moulin Shape model, The Cryosphere, 16, 2421–2448, https://doi.org/10.5194/tc-16-2421-2022, 2022. a, b, c, d
Aschwanden, A., Bartholomaus, T. C., Brinkerhoff, D. J., and Truffer, M.: Brief communication: A roadmap towards credible projections of ice sheet contribution to sea level, The Cryosphere, 15, 5705–5715, https://doi.org/10.5194/tc-15-5705-2021, 2021. a
Bamber, J. L., Oppenheimer, M., Kopp, R. E., Aspinall, W. P., and Cooke, R. M.: Ice sheet contributions to future sea-level rise from structured expert judgment, P. Natl. Acad. Sci. USA, 166, 11195–11200, https://doi.org/10.1073/pnas.1817205116, 2019. a
Benn, D. I., Hulton, N. R., and Mottram, R. H.: 'Calving laws', 'sliding laws' and the stability of tidewater glaciers, Ann. Glaciol., 46, 123–130, https://doi.org/10.3189/172756407782871161, 2007. a
Bevis, M., Harig, C., Khan, S. A., Brown, A., Simons, F. J., Willis, M. J., Fettweis, X., van den Broeke, M. R., Madsen, F. B., Kendrick, E., Caccamise, D. J., van Dam, T., Knudsen, P., and Nylen, T.: Accelerating changes in ice mass within Greenland, and the ice sheet's sensitivity to atmospheric forcing, P. Natl. Acad. Sci. USA, 116, 1934–1939, https://doi.org/10.1073/pnas.1806562116, 2019. a
Boone, T. J. and Ingraffea, A. R.: A numerical procedure for simulation of hydraulically-driven fracture propagation in poroelastic media, Int. J. Numer. Anal. Met., 14, 27–47, https://doi.org/10.1002/nag.1610140103, 1990. a
Bueler, E. and van Pelt, W.: Mass-conserving subglacial hydrology in the Parallel Ice Sheet Model version 0.6, Geosci. Model Dev., 8, 1613–1635, https://doi.org/10.5194/gmd-8-1613-2015, 2015. a
Carrier, B. and Granet, S.: Numerical modeling of hydraulic fracture problem in permeable medium using cohesive zone model, Eng. Fract. Mech., 79, 312–328, https://doi.org/10.1016/j.engfracmech.2011.11.012, 2012. a
Christoffersen, P., Bougamont, M., Hubbard, A., Doyle, S. H., Grigsby, S., and Pettersson, R.: Cascading lake drainage on the Greenland Ice Sheet triggered by tensile shock and fracture, Nat. Commun., 9, 1–12, https://doi.org/10.1038/s41467-018-03420-8, 2018. a
Chudley, T. R., Christoffersen, P., Doyle, S. H., Bougamont, M., Schoonman, C. M., Hubbard, B., and James, M. R.: Supraglacial lake drainage at a fast-flowing Greenlandic outlet glacier, P. Natl. Acad. Sci. USA, 116, 25468–25477, https://doi.org/10.1073/PNAS.1913685116, 2019. a, b, c
Clayton, T., Duddu, R., Siegert, M., and Martínez-Pañeda, E.: A stress-based poro-damage phase field model for hydrofracturing of creeping glaciers and ice shelves, Eng. Fract. Mech., 272, 108693, https://doi.org/10.1016/J.ENGFRACMECH.2022.108693, 2022. a, b
Crawford, A. J., Benn, D. I., Todd, J., Åström, J. A., Bassis, J. N., and Zwinger, T.: Marine Ice-Cliff Instability Modeling Shows Mixed-Mode Ice-Cliff Failure and Yields Calving Rate Parameterization, Nat. Commun., 12, 1–9, https://doi.org/10.1038/s41467-021-23070-7, 2021. a
de Borst, R.: Fluid flow in fractured and fracturing porous media: A unified view, Mech. Res. Commun., 80, 47–57, https://doi.org/10.1016/j.mechrescom.2016.05.004, 2017. a
de Fleurian, B., Gagliardini, O., Zwinger, T., Durand, G., Le Meur, E., Mair, D., and Råback, P.: A double continuum hydrological model for glacier applications, The Cryosphere, 8, 137–153, https://doi.org/10.5194/tc-8-137-2014, 2014. a
Desroches, J., Detournay, E., Lenoach, B., Papanastasiou, P., Pearson, J. R. A., Thiercelin, M., and Cheng, A.: The crack tip region in hydraulic fracturing, P. Roy. Soc. Lond. A-Mat., 447, 39–48, 1994. a
Doyle, S. H., Hubbard, A. L., Dow, C. F., Jones, G. A., Fitzpatrick, A., Gusmeroli, A., Kulessa, B., Lindback, K., Pettersson, R., and Box, J. E.: Ice tectonic deformation during the rapid in situ drainage of a supraglacial lake on the Greenland Ice Sheet, The Cryosphere, 7, 129–140, https://doi.org/10.5194/tc-7-129-2013, 2013. a, b, c
Duddu, R. and Waisman, H.: A temperature dependent creep damage model for polycrystalline ice, Mech. Mater., 46, 23–41, https://doi.org/10.1016/J.MECHMAT.2011.11.007, 2012. a, b
Duddu, R., Jiménez, S., and Bassis, J.: A non-local continuum poro-damage mechanics model for hydrofracturing of surface crevasses in grounded glaciers, J. Glaciol., 66, 415–429, 2020. a
Gauckler, P.: Etudes Théoriques et Pratiques sur l'Ecoulement et le Mouvement des Eaux, C. R. Acad. Sci. Paris, 64, 818–822, 1867. a
Ghosh, G., Duddu, R., and Annavarapu, C.: A stabilized finite element method for enforcing stiff anisotropic cohesive laws using interface elements, Computer Method. Appl. M., 348, 1013–1038, 2019. a
Glen, J. W.: The creep of polycrystalline ice, P. Roy. Soc. Lond. A-Mat., 228, 519–538, https://doi.org/10.1098/rspa.1955.0066, 1955. a, b
Goelzer, H., Robinson, A., Seroussi, H., and van de Wal, R. S.: Recent Progress in Greenland Ice Sheet Modelling, Current Climate Change Reports, 3, 291–302, https://doi.org/10.1007/s40641-017-0073-y, 2017. a
Goelzer, H., Nowicki, S., Payne, A., Larour, E., Seroussi, H., Lipscomb, W. H., Gregory, J., Abe-Ouchi, A., Shepherd, A., Simon, E., Agosta, C., Alexander, P., Aschwanden, A., Barthel, A., Calov, R., Chambers, C., Choi, Y., Cuzzone, J., Dumas, C., Edwards, T., Felikson, D., Fettweis, X., Golledge, N. R., Greve, R., Humbert, A., Huybrechts, P., Le clec'h, S., Lee, V., Leguy, G., Little, C., Lowry, D. P., Morlighem, M., Nias, I., Quiquet, A., Rückamp, M., Schlegel, N.-J., Slater, D. A., Smith, R. S., Straneo, F., Tarasov, L., van de Wal, R., and van den Broeke, M.: The future sea-level contribution of the Greenland ice sheet: a multi-model ensemble study of ISMIP6, The Cryosphere, 14, 3071–3096, https://doi.org/10.5194/tc-14-3071-2020, 2020. a
Greve, R., Zwinger, T., and Gong, Y.: On the pressure dependence of the rate factor in Glen's flow law, J. Glaciol., 60, 397–399, https://doi.org/10.3189/2014JOG14J019, 2014. a
Hageman, T.: T-Hageman/MATLAB_IceHydroFrac: Final version (final), Zenodo [code], https://doi.org/10.5281/zenodo.13629447, 2024. a
Hageman, T. and de Borst, R.: Flow of non-Newtonian fluids in fractured porous media: Isogeometric vs standard finite element discretisation, Int. J. Numer. Anal. Met., 43, 2020–2037, https://doi.org/10.1002/nag.2948, 2019. a
Hageman, T. and de Borst, R.: A refined two-scale model for Newtonian and non-Newtonian fluids in fractured poroelastic media, J. Comput. Phys., 441, 110424, https://doi.org/10.1016/j.jcp.2021.110424, 2021. a, b
Hageman, T. and de Borst, R.: Direct simulation vs subgrid scale modelling of fluid flow in fractured or fracturing porous media, Computat. Geosci., 26, 503–515, https://doi.org/10.1007/s10596-022-10138-6, 2022. a
Hageman, T., Pervaiz Fathima, K. M., and de Borst, R.: Isogeometric analysis of fracture propagation in saturated porous media due to a pressurised non-Newtonian fluid, Comput. Geotech., 112, 272–283, https://doi.org/10.1016/j.compgeo.2019.04.030, 2019. a
Hageman, T., Mejía, J., Duddu, R., and Martínez-Pañeda, E.: Ice viscosity governs hydraulic fracture causing rapid drainage of supraglacial lakes, TIB AV-Portal [video], https://doi.org/10.5446/68350, 2024. a
Hewitt, I. J., Schoof, C., and Werder, M. A.: Flotation and free surface flow in a model for subglacial drainage. Part 2. Channel flow, J. Fluid Mech., 702, 157–187, https://doi.org/10.1017/JFM.2012.166, 2012. a
Hofer, S., Lang, C., Amory, C., Kittel, C., Delhasse, A., Tedstone, A., and Fettweis, X.: Greater Greenland Ice Sheet contribution to global sea level rise in CMIP6, Nat. Commun., 11, 6289, https://doi.org/10.1038/s41467-020-20011-8, 2020. a
Hoffman, M. J., Catania, G. A., Neumann, T. A., Andrews, L. C., and Rumrill, J. A.: Links between acceleration, melting, and supraglacial lake drainage of the western Greenland Ice Sheet, J. Geophys. Res.-Earth, 116, 1–16, https://doi.org/10.1029/2010JF001934, 2011. a, b
Hoffman, M. J., Perego, M., Andrews, L. C., Catania, G. A., Price, S. F., Lüthi, M. P., Neumann, T. A., and Johnson, J. V.: Widespread Moulin Formation During Supraglacial Lake Drainages in Greenland, Geophys. Res. Lett., 45, 778–788, https://doi.org/10.1002/2017GL075659, 2018. a
Howat, I. M., de la Peña, S., van Angelen, J. H., Lenaerts, J. T. M., and van den Broeke, M. R.: Brief Communication “Expansion of meltwater lakes on the Greenland Ice Sheet”, The Cryosphere, 7, 201–204, https://doi.org/10.5194/tc-7-201-2013, 2013. a
Imperial College London research computing service: https://doi.org/10.14469/hpc/2232, 2024. a
Jellinek, H. H. G. and Brill, R.: Viscoelastic Properties of Ice, J. Appl. Phys., 27, 1198–1209, https://doi.org/10.1063/1.1722231, 1956. a
Jezek, K. C.: A modified theory of bottom crevasses used as a means for measuring the buttressing effect of ice shelves on inland ice sheets, J. Geophys. Res.-Sol. Ea., 89, 1925–1931, https://doi.org/10.1029/JB089IB03P01925, 1984. a
Jimenez, S. and Duddu, R.: On the evaluation of the stress intensity factor in calving models using linear elastic fracture mechanics, J. Glaciol., 64, 759–770, 2018. a
Krawczynski, M. J., Behn, M. D., Das, S. B., and Joughin, I.: Constraints on the lake volume required for hydro-fracture through ice sheets, Geophys. Res. Lett., 36, 10501, https://doi.org/10.1029/2008GL036765, 2009. a
Lai, C. Y., Stevens, L. A., Chase, D. L., Creyts, T. T., Behn, M. D., Das, S. B., and Stone, H. A.: Hydraulic Transmissivity Inferred from Ice-Sheet Relaxation Following Greenland Supraglacial Lake Drainages, Nat. Commun., 12, 1–10, https://doi.org/10.1038/s41467-021-24186-6, 2021. a, b, c
Larour, E., Seroussi, H., Morlighem, M., and Rignot, E.: Continental scale, high order, high spatial resolution, ice sheet modeling using the Ice Sheet System Model (ISSM), J. Geophys. Res.-Earth, 117, 1022, https://doi.org/10.1029/2011JF002140, 2012. a
Liang, Y. L., Colgan, W., Lv, Q., Steffen, K., Abdalati, W., Stroeve, J., Gallaher, D., and Bayou, N.: A decadal investigation of supraglacial lakes in West Greenland using a fully automatic detection and tracking algorithm, Remote Sens. Environ., 123, 127–138, https://doi.org/10.1016/J.RSE.2012.03.020, 2012. a
Lipscomb, W. H., Price, S. F., Hoffman, M. J., Leguy, G. R., Bennett, A. R., Bradley, S. L., Evans, K. J., Fyke, J. G., Kennedy, J. H., Perego, M., Ranken, D. M., Sacks, W. J., Salinger, A. G., Vargo, L. J., and Worley, P. H.: Description and evaluation of the Community Ice Sheet Model (CISM) v2.1, Geosci. Model Dev., 12, 387–424, https://doi.org/10.5194/gmd-12-387-2019, 2019. a, b
Litwin, K. L., Zygielbaum, B. R., Polito, P. J., Sklar, L. S., and Collins, G. C.: Influence of temperature, composition, and grain size on the tensile failure of water ice: Implications for erosion on Titan, J. Geophys. Res.-Planets, 117, 8013, https://doi.org/10.1029/2012JE004101, 2012. a
Mejía, J. Z., Gulley, J. D., Trunz, C., Covington, M. D., Bartholomaus, T. C., Xie, S., and Dixon, T.: Isolated cavities dominate Greenland Ice Sheet dynamic response to lake drainage, Geophys. Res. Lett., 48, 1–11, https://doi.org/10.1029/2021gl094762, 2021. a, b, c
Nick, F. M., Van Der Veen, C. J., Vieli, A., and Benn, D. I.: A physically based calving model applied to marine outlet glaciers and implications for the glacier dynamics, J. Glaciol., 56, 781–794, https://doi.org/10.3189/002214310794457344, 2010. a
Pattyn, F.: The Paradigm Shift in Antarctic Ice Sheet Modelling, Nat. Commun., 9, 1–3, https://doi.org/10.1038/s41467-018-05003-z, 2018. a
Pimentel, S. and Flowers, G. E.: A numerical study of hydrologically driven glacier dynamics and subglacial flooding, P. Roy. Soc. A-Math. Phy., 467, 537–558, https://doi.org/10.1098/RSPA.2010.0211, 2011. a
Poinar, K., Joughin, I., Lilien, D., Brucker, L., Kehrl, L., and Nowicki, S.: Drainage of southeast Greenland firn aquifer water through crevasses to the bed, Front. Earth Sci., 5, https://doi.org/10.3389/feart.2017.00005, 2017. a
Réthoré, J., de Borst, R., and Abellan, M.-A.: A two-scale approach for fluid flow in fractured porous media, Int. J. Numer. Meth. Eng., 71, 780–800, https://doi.org/10.1002/nme.1962, 2006. a
Rice, J. R., Tsai, V. C., Fernandes, M. C., and Platt, J. D.: Time scale for rapid draining of a surficial lake into the Greenland ice sheet, J. Appl. Mech.-T. ASME, 82, 071001, https://doi.org/10.1115/1.4030325, 2015. a
Ryser, C., Lüthi, M. P., Andrews, L. C., Catania, G. A., Funk, M., Hawley, R., Hoffman, M., and Neumann, T. A.: Caterpillar-like ice motion in the ablation zone of the Greenland ice sheet, J. Geophys. Res.-Earth, 119, 2258–2271, https://doi.org/10.1002/2013JF003067, 2014a. a
Ryser, C., Luthi, M. P., Andrews, L. C., Hoffman, M. J., Catania, G. A., Hawley, R. L., Neumann, T. A., and Kristensen, S. S.: Sustained high basal motion of the Greenland ice sheet revealed by borehole deformation, J. Glaciol., 60, 647–660, https://doi.org/10.3189/2014JOG13J196, 2014b. a, b
Selmes, N., Murray, T., and James, T. D.: Fast draining lakes on the Greenland Ice Sheet, Geophys. Res. Lett., 38, 15501, https://doi.org/10.1029/2011GL047872, 2011. a, b
Seroussi, H., Morlighem, M., Rignot, E., Larour, E., Aubry, D., Ben Dhia, H., and Kristensen, S. S.: Ice Flux Divergence Anomalies on 79north Glacier, Greenland, Geophys. Res. Lett., 38, https://doi.org/10.1029/2011GL047338, 2011. a
Shannon, S. R., Payne, A. J., Bartholomew, I. D., Broeke, M. R. V. D., Edwards, T. L., Fettweis, X., Gagliardini, O., Gillet-Chaulet, F., Goelzer, H., Hoffman, M. J., Huybrechts, P., Mair, D. W., Nienow, P. W., Perego, M., Price, S. F., Smeets, C. J. P., Sole, A. J., Wal, R. S. V. D., and Zwinger, T.: Enhanced basal lubrication and the contribution of the Greenland ice sheet to future sea-level rise, P. Natl. Acad. Sci. USA, 110, 14156–14161, https://doi.org/10.1073/pnas.1212647110, 2013. a
Smith, L. C., Chu, V. W., Yang, K., Gleason, C. J., Pitcher, L. H., Rennermalm, A. K., Legleiter, C. J., Behar, A. E., Overstreet, B. T., Moustafa, S. E., Tedesco, M., Forster, R. R., LeWinter, A. L., Finnegan, D. C., Sheng, Y., and Balog, J.: Efficient meltwater drainage through supraglacial streams and rivers on the southwest Greenland ice sheet, P. Natl. Acad. Sci. USA, 112, 1001–1006, https://doi.org/10.1073/pnas.1413024112, 2015. a
Smith, R. A.: The Application of Fracture Mechanics to the Problem of Crevasse Penetration, J. Glaciol., 17, 223–228, https://doi.org/10.3189/S0022143000013563, 1976. a
Stevens, L. A., Behn, M. D., McGuire, J. J., Das, S. B., Joughin, I., Herring, T., Shean, D. E., and King, M. A.: Greenland supraglacial lake drainages triggered by hydrologically induced basal slip, Nature, 522, 73–76, https://doi.org/10.1038/nature14480, 2015. a, b, c, d
Stevens, L. A., Nettles, M., Davis, J. L., Creyts, T. T., Kingslake, J., Hewitt, I. J., and Stubblefield, A.: Tidewater-glacier response to supraglacial lake drainage, Nat. Commun., 13, 6065, https://doi.org/10.1038/s41467-022-33763-2, 2022. a
Sun, X., Duddu, R., and Hirshikesh: A poro-damage phase field model for hydrofracturing of glacier crevasses, Extreme Mechanics Letters, 45, 101277, https://doi.org/10.1016/j.eml.2021.101277, 2021. a
Trunz, C., Covington, M. D., Poinar, K., Andrews, L. C., Mejia, J., and Gulley, J.: Modeling the Influence of Moulin Shape on Subglacial Hydrology, J. Geophys. Res.-Earth, 127, e2022JF006674, https://doi.org/10.1029/2022JF006674, 2022. a
Tsai, V. C. and Rice, J. R.: A model for turbulent hydraulic fracture and application to crack propagation at glacier beds, J. Geophys. Res.-Earth, 115, 3007, https://doi.org/10.1029/2009JF001474, 2010. a, b, c
Tsai, V. C. and Rice, J. R.: Modeling Turbulent Hydraulic Fracture Near a Free Surface, J. Appl. Mech., 79, 031003, https://doi.org/10.1115/1.4005879, 2012. a, b, c
Van Der Veen, C. J.: Fracture mechanics approach to penetration of surface crevasses on glaciers, Cold Regions Science and Technology, 27, 31–47, https://doi.org/10.1016/S0165-232X(97)00022-0, 1998. a
van der Veen, C. J.: Fracture propagation as means of rapidly transferring surface meltwater to the base of glaciers, Geophys. Res. Lett., 34, L01501, https://doi.org/10.1029/2006GL028385, 2007. a, b
Weertman, J.: Theory of water-filled crevasses in glaciers applied to vertical magma transport beneath oceanic ridges, J. Geophys. Res., 76, 1171–1183, https://doi.org/10.1029/JB076I005P01171, 1971. a, b
Weertman, J.: Can a water filled crevasse reach the bottom surface of a glacier?, IASH Publ., 95, 139–145, https://doi.org/10.1017/CBO9781107415324.004, 1973. a, b
Weertman, J.: Creep deformation of ice, Annu. Rev. Earth Pl. Sc., 11, 215–240, https://doi.org/10.1146/annurev.ea.11.050183.001243, 1983. a
Witherspoon, P. A., Wang, J. S. Y., Iwai, K., and Gale, J. E.: Validity of Cubic Law for fluid flow in a deformable rock fracture, Water Resour. Res., 16, 1016–1024, https://doi.org/10.1029/WR016i006p01016, 1980. a
Zarrinderakht, M., Schoof, C., and Peirce, A.: The effect of hydrology and crevasse wall contact on calving, The Cryosphere, 16, 4491–4512, https://doi.org/10.5194/tc-16-4491-2022, 2022. a, b
Co-editor-in-chief
The study is one of the first to model fractures in ice sheets - a fascinating and visually stunning aspect of ice sheets. The model shows that crevasses may transport large volumes of water to the bed of a glacier very quickly and captures the opening of the crevasses due to the water inflow. The impact of surface lakes on the Greenland ice sheet dynamics and mass loss is now better described.
The study is one of the first to model fractures in ice sheets - a fascinating and visually...
Short summary
Due to surface melting, meltwater lakes seasonally form on the surface of glaciers. These lakes drive hydrofractures that rapidly transfer water to the base of ice sheets. This paper presents a computational method to capture the complicated hydrofracturing process. Our work reveals that viscous ice rheology has a great influence on the short-term propagation of fractures, enabling fast lake drainage, whereas thermal effects (frictional heating, conduction, and freezing) have little influence.
Due to surface melting, meltwater lakes seasonally form on the surface of glaciers. These lakes...