scholarly journals Adaptive mesh refinement versus subgrid friction interpolation in simulations of Antarctic ice dynamics

2016 ◽  
Vol 57 (73) ◽  
pp. 1-9 ◽  
Author(s):  
S. L. Cornford ◽  
D. F. Martin ◽  
V. Lee ◽  
A. J. Payne ◽  
E. G. Ng
Keyword(s):  
Adaptive Mesh Refinement ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
Ice Streams ◽  
Order Error ◽  
Ice Dynamics ◽  
West Antarctic ◽  
Grounding Line ◽  
Improve Accuracy ◽  
Mesh Spacing

ABSTRACTAt least in conventional hydrostatic ice-sheet models, the numerical error associated with grounding line dynamics can be reduced by modifications to the discretization scheme. These involve altering the integration formulae for the basal traction and/or driving stress close to the grounding line and exhibit lower – if still first-order – error in the MISMIP3d experiments. MISMIP3d may not represent the variety of real ice streams, in that it lacks strong lateral stresses, and imposes a large basal traction at the grounding line. We study resolution sensitivity in the context of extreme forcing simulations of the entire Antarctic ice sheet, using the BISICLES adaptive mesh ice-sheet model with two schemes: the original treatment, and a scheme, which modifies the discretization of the basal traction. The second scheme does indeed improve accuracy – by around a factor of two – for a given mesh spacing, but $\lesssim 1$ km resolution is still necessary. For example, in coarser resolution simulations Thwaites Glacier retreats so slowly that other ice streams divert its trunk. In contrast, with $\lesssim 1$ km meshes, the same glacier retreats far more quickly and triggers the final phase of West Antarctic collapse a century before any such diversion can take place.

scholarly journals Contrasting the modelled sensitivity of the Amundsen Sea Embayment ice streams

2016 ◽  
Vol 62 (233) ◽  
pp. 552-562 ◽  
Author(s):  
ISABEL J. NIAS ◽  
STEPHEN L. CORNFORD ◽  
ANTONY J. PAYNE
Keyword(s):  
Sea Level Rise ◽  
Sea Level ◽  
Adaptive Mesh Refinement ◽  
Adaptive Mesh ◽  
Mass Conservation ◽  
Physical Parameters ◽  
Ice Streams ◽  
Amundsen Sea ◽  
West Antarctic ◽  
Grounding Line

AbstractPresent-day mass loss from the West Antarctic ice sheet is centred on the Amundsen Sea Embayment (ASE), primarily through ice streams, including Pine Island, Thwaites and Smith glaciers. To understand the differences in response of these ice streams, we ran a perturbed parameter ensemble, using a vertically-integrated ice flow model with adaptive mesh refinement. We generated 71 sets of three physical parameters (basal traction coefficient, ice viscosity stiffening factor and sub-shelf melt rate), which we used to simulate the ASE for 50 years. We also explored the effects of different bed geometries and basal sliding laws. The mean rate of sea-level rise across the ensemble of simulations is comparable with current observed rates for the ASE. We found evidence that grounding line dynamics are sensitive to features in the bed geometry: simulations using BedMap2 geometry resulted in a higher rate of sea-level rise than simulations using a rougher geometry, created using mass conservation. Modelled grounding-line retreat of all the three ice streams was sensitive to viscosity and basal traction, while the melt rate was more important in Pine Island and Smith glaciers, which flow through more confined ice shelves than Thwaites, which has a relatively unconfined shelf.

scholarly journals Ocean-forced evolution of the Amundsen Sea catchment, West Antarctica, by 2100

2020 ◽  
Vol 14 (4) ◽  
pp. 1245-1258
Author(s):  
Alanna V. Alevropoulos-Borrill ◽  
Isabel J. Nias ◽  
Antony J. Payne ◽  
Nicholas R. Golledge ◽  
Rory J. Bingham
Keyword(s):  
Sea Level ◽  
21St Century ◽  
Adaptive Mesh Refinement ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
Ocean Temperature ◽  
Ice Shelf ◽  
Ice Streams ◽  
Amundsen Sea ◽  
Grounding Line

Abstract. The response of ice streams in the Amundsen Sea Embayment (ASE) to future climate forcing is highly uncertain. Here we present projections of 21st century response of ASE ice streams to modelled local ocean temperature change using a subset of Coupled Model Intercomparison Project (CMIP5) simulations. We use the BISICLES adaptive mesh refinement (AMR) ice sheet model, with high-resolution grounding line resolving capabilities, to explore grounding line migration in response to projected sub-ice-shelf basal melting. We find a contribution to sea level rise of between 2.0 and 4.5 cm by 2100 under RCP8.5 conditions from the CMIP5 subset, where the mass loss response is linearly related to the mean ocean temperature anomaly. To account for uncertainty associated with model initialization, we perform three further sets of CMIP5-forced experiments using different parameterizations that explore perturbations to the prescription of initial basal melt, the basal traction coefficient and the ice stiffening factor. We find that the response of the ASE to ocean temperature forcing is highly dependent on the parameter fields obtained in the initialization procedure, where the sensitivity of the ASE ice streams to the sub-ice-shelf melt forcing is dependent on the choice of parameter set. Accounting for ice sheet model parameter uncertainty results in a projected range in sea level equivalent contribution from the ASE of between −0.02 and 12.1 cm by the end of the 21st century.

scholarly journals Ocean forced evolution of the Amundsen Sea catchment, West Antarctica, by 2100

2019 ◽  
Author(s):  
Alanna V. Alevropoulos-Borrill ◽  
Isabel J. Nias ◽  
Antony J. Payne ◽  
Nicholas R. Golledge ◽  
Rory J. Bingham
Keyword(s):  
Sea Level ◽  
21St Century ◽  
Adaptive Mesh Refinement ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
Ocean Temperature ◽  
Ice Shelf ◽  
Ice Streams ◽  
Amundsen Sea ◽  
Grounding Line

Abstract. The response of ice streams in the Amundsen Sea Embayment (ASE) to future climate forcing is highly uncertain. Here we present projections of 21st century response of ASE ice streams to modelled local ocean temperature change using a subset of Coupled Model Intercomparison Project (CMIP5) simulations. We use the BISICLES adaptive mesh refinement (AMR) ice sheet model, with high resolution grounding line resolving capabilities, to explore grounding line migration in response to projected sub-ice shelf basal melting. We find a contribution to sea level rise of between 2.0 cm and 4.5 cm by 2100 under RCP8.5 conditions from the CMIP5 subset, where the mass loss response is linearly related to the mean ocean temperature anomaly. To account for uncertainty associated with model initialisation, we perform three further sets of CMIP5 forced experiments using different parameterisations that explore perturbations to the prescription of initial basal melt, the basal traction coefficient, and the ice stiffening factor. We find that the response of the ASE to ocean temperature forcing is highly dependent on the parameter fields obtained in the initialisation procedure, where the sensitivity of the ASE ice streams to the sub-ice shelf melt forcing is dependent on the choice of parameter set. Accounting for ice sheet model parameter uncertainty results in a projected range in sea level equivalent contribution from the ASE of between −0.02 cm and 12.1 cm by the end of the 21st century.

scholarly journals Results of the Marine Ice Sheet Model Intercomparison Project, MISMIP

2012 ◽  
Vol 6 (3) ◽  
pp. 573-588 ◽  
Author(s):  
F. Pattyn ◽  
C. Schoof ◽  
L. Perichon ◽  
R. C. A. Hindmarsh ◽  
E. Bueler ◽  
...  
Keyword(s):  
Steady State ◽  
Adaptive Mesh Refinement ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
Stable Steady State ◽  
Fixed Grid ◽  
Approximate Analytical Solutions ◽  
Grounding Line ◽  
Intercomparison Project ◽  
Line Positions

Abstract. Predictions of marine ice-sheet behaviour require models that are able to robustly simulate grounding line migration. We present results of an intercomparison exercise for marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no effects of lateral buttressing). Unique steady state grounding line positions exist for ice sheets on a downward sloping bed, while hysteresis occurs across an overdeepened bed, and stable steady state grounding line positions only occur on the downward-sloping sections. Models based on the shallow ice approximation, which does not resolve extensional stresses, do not reproduce the approximate analytical results unless appropriate parameterizations for ice flux are imposed at the grounding line. For extensional-stress resolving "shelfy stream" models, differences between model results were mainly due to the choice of spatial discretization. Moving grid methods were found to be the most accurate at capturing grounding line evolution, since they track the grounding line explicitly. Adaptive mesh refinement can further improve accuracy, including fixed grid models that generally perform poorly at coarse resolution. Fixed grid models, with nested grid representations of the grounding line, are able to generate accurate steady state positions, but can be inaccurate over transients. Only one full-Stokes model was included in the intercomparison, and consequently the accuracy of shelfy stream models as approximations of full-Stokes models remains to be determined in detail, especially during transients.

scholarly journals Dynamical processes involved in the retreat of marine ice sheets

2001 ◽  
Vol 47 (157) ◽  
pp. 271-282 ◽  
Author(s):  
Richard C.A. Hindmarsh ◽  
E. Le Meur
Keyword(s):  
Ice Sheets ◽  
Ice Sheet ◽  
Ice Shelf ◽  
Ice Streams ◽  
Antarctic Ice Sheet ◽  
West Antarctic Ice Sheet ◽  
West Antarctic ◽  
Grounding Line ◽  
Neutral Equilibrium ◽  
Generally True

AbstractMarine ice sheets with mechanics described by the shallow-ice approximation by definition do not couple mechanically with the shelf. Such ice sheets are known to have neutral equilibria. We consider the implications of this for their dynamics and in particular for mechanisms which promote marine ice-sheet retreat. The removal of ice-shelf buttressing leading to enhanced flow in grounded ice is discounted as a significant influence on mechanical grounds. Sea-level rise leading to reduced effective pressures under ice streams is shown to be a feasible mechanism for producing postglacial West Antarctic ice-sheet retreat but is inconsistent with borehole evidence. Warming thins the ice sheet by reducing the average viscosity but does not lead to grounding-line retreat. Internal oscillations either specified or generated via a MacAyeal–Payne thermal mechanism promote migration. This is a noise-induced drift phenomenon stemming from the neutral equilibrium property of marine ice sheets. This migration occurs at quite slow rates, but these are sufficiently large to have possibly played a role in the dynamics of the West Antarctic ice sheet after the glacial maximum. Numerical experiments suggest that it is generally true that while significant changes in thickness can be caused by spatially uniform changes, spatial variability coupled with dynamical variability is needed to cause margin movement.

scholarly journals Implementation and performance of adaptive mesh refinement in the Ice Sheet System Model (ISSM v4.14)

2019 ◽  
Vol 12 (1) ◽  
pp. 215-232 ◽  
Author(s):  
Thiago Dias dos Santos ◽  
Mathieu Morlighem ◽  
Hélène Seroussi ◽  
Philippe Remy Bernard Devloo ◽  
Jefferson Cardia Simões
Keyword(s):  
Adaptive Mesh Refinement ◽  
Mesh Refinement ◽  
Computational Cost ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
System Model ◽  
Computational Time ◽  
Error Estimator ◽  
Fine Mesh ◽  
Grounding Line

Abstract. Accurate projections of the evolution of ice sheets in a changing climate require a fine mesh/grid resolution in ice sheet models to correctly capture fundamental physical processes, such as the evolution of the grounding line, the region where grounded ice starts to float. The evolution of the grounding line indeed plays a major role in ice sheet dynamics, as it is a fundamental control on marine ice sheet stability. Numerical modeling of a grounding line requires significant computational resources since the accuracy of its position depends on grid or mesh resolution. A technique that improves accuracy with reduced computational cost is the adaptive mesh refinement (AMR) approach. We present here the implementation of the AMR technique in the finite element Ice Sheet System Model (ISSM) to simulate grounding line dynamics under two different benchmarks: MISMIP3d and MISMIP+. We test different refinement criteria: (a) distance around the grounding line, (b) a posteriori error estimator, the Zienkiewicz–Zhu (ZZ) error estimator, and (c) different combinations of (a) and (b). In both benchmarks, the ZZ error estimator presents high values around the grounding line. In the MISMIP+ setup, this estimator also presents high values in the grounded part of the ice sheet, following the complex shape of the bedrock geometry. The ZZ estimator helps guide the refinement procedure such that AMR performance is improved. Our results show that computational time with AMR depends on the required accuracy, but in all cases, it is significantly shorter than for uniformly refined meshes. We conclude that AMR without an associated error estimator should be avoided, especially for real glaciers that have a complex bed geometry.

scholarly journals Isostatic Equilibrium Grounding Line Between The West Antarctic Ice Sheet And The Ross Ice Shelf

1979 ◽  
Vol 24 (90) ◽  
pp. 500 ◽  
Author(s):  
C. R. Bentley ◽  
L. Greischar
Keyword(s):  
Ice Sheet ◽  
Shelf Area ◽  
Ice Shelf ◽  
Ice Streams ◽  
West Antarctic Ice Sheet ◽  
Ross Ice Shelf ◽  
Isostatic Rebound ◽  
West Antarctic ◽  
Grounding Line ◽  
Eastern Half

Abstract Taking various retreat-rates for the presumed grounded ice sheet in the Ross embayment during Wisconsin time, as calculated by Thomas (Thomas and Bentley, 1978), and assuming a time constant of 4400 years for isostatic rebound, a sea-floor uplift of 100±50 m still to be expected in the grid western part of the Ross Ice Shelf can be calculated. The expected uplift diminishes from grid west to grid east, and is probably negligible in the eastern half of the shelf area. There are extensive areas near the present grounding line where the water depth beneath the shelf is less than 100 m, so that uplift would lead to grounding. As grounding occurred, the neighboring ice shelf would thicken, causing grounding to advance farther. This process would probably extend the grounding line to a position running grid north-eastward across the shelf from the seaward end of Roosevelt Island, deeply indented by the extensions of the present ice streams. Floating ice would remain in the grid south-eastern half of the shelf.

scholarly journals Implementation and Performance of Adaptive Mesh Refinement in the Ice Sheet System Model (ISSM v4.14)

2018 ◽  
Author(s):  
Thiago Dias dos Santos ◽  
Mathieu Morlighem ◽  
Hélène Seroussi ◽  
Philippe Remy Bernard Devloo ◽  
Jefferson Cardia Simões
Keyword(s):  
Adaptive Mesh Refinement ◽  
Mesh Refinement ◽  
Computational Cost ◽  
Ice Sheet ◽  
Adaptive Mesh ◽  
System Model ◽  
Computational Time ◽  
Error Estimator ◽  
Grounding Line ◽  
Refined Meshes

Abstract. Accurate projections of the evolution of ice sheets in a changing climate require a fine mesh/grid resolution to correctly capture fundamental physical processes, such as the evolution of the grounding line, the region where grounded ice starts to float. The evolution of the grounding line indeed plays a major role in ice sheet dynamics, as it is a fundamental control on marine ice sheet stability. Numerical modeling of grounding line requires significant computational resources since the accuracy of its position depends on grid or mesh resolution. A technique that improves accuracy with reduced computational cost is the adaptive mesh refinement approach, AMR. We present here the implementation of the AMR technique in the finite element Ice Sheet System Model (ISSM) to simulate grounding line dynamics under two different benchmarks, MISMIP3d and MISMIP+. We test different refinement criteria: (a) distance around grounding line, (b) a posteriori error estimator, the Zienkiewicz-Zhu (ZZ) error estimator, and (c) different combinations of (a) and (b). We find that for MISMIP3d setup, refining 5 km around the grounding line, both on grounded and floating ice, is sufficient to produce AMR results similar to the ones obtained with uniformly refined meshes. However, for the MISMIP+ setup, we note that there is a minimum distance of 30 km around the grounding line required to produce accurate results. We find this AMR mesh-dependency is linked to the complex bedrock topography of MISMIP+. In both benchmarks, the ZZ error estimator presents high values around the grounding line. Particularly for MISMIP+ setup, the estimator also presents high values in the grounded part of the ice sheet, following the complex shape of the bedrock geometry. This estimator helps guide the refinement procedure such that AMR performance is improved. Our results show that computational time with AMR depends on the required accuracy, but in all cases, it is significantly shorter than for uniformly refined meshes. We conclude that AMR without an associated error estimator should be avoided, especially for real glaciers that have a complex bed geometry.

scholarly journals The Dynamics of Marine Ice Sheets

1979 ◽  
Vol 24 (90) ◽  
pp. 167-177 ◽  
Author(s):  
Robert H. Thomas
Keyword(s):  
Sea Level ◽  
Ice Sheets ◽  
Ice Sheet ◽  
Ice Shelf ◽  
Ice Streams ◽  
Ice Shelves ◽  
The North ◽  
West Antarctic ◽  
Grounding Line ◽  
Ice Sheet Dynamics

AbstractMarine ice sheets rest on land that, for the most part, is below sea-level. Ice that flows across the grounding line, where the ice sheet becomes afloat, either calves into icebergs or forms a floating ice shelf joined to the ice sheet. At the grounding line there is a transition from ice-sheet dynamics to ice-shelf dynamics, and the creep-thinning rate in this region is very sensitive to sea depth; rising sea-level causes increased thinning-rates and grounding-line retreat, falling sea-level has the reverse effect. If the bedrock slopes down towards the centre of the ice sheet there may be only two stable modes: a freely-floating ice shelf or a marine ice sheet that extends to the edge of the continental shelf. Once started, collapse of such an ice sheet to form an ice shelf may take place extremely rapidly. Ice shelves which form in embayments of a marine ice sheet, or which are partially grounded, have a stabilizing influence since ice flowing across the grounding line has to push the ice shelf past its sides. Retreat of the grounding line tends to enlarge the ice shelf, which ultimately may become large enough to prevent excessive outflow from the ice sheet so that a new equilibrium grounding line is established; removal of the ice shelf would allow retreat to continue. During the late-Wisconsin glacial maximum there may have been marine ice sheets in the northern hemisphere but the only current example is the West Antarctic ice sheet. This is buttressed by the Ross and Ronne Ice Shelves, and if climatic warming were to prohibit the existence of these ice shelves then the ice sheet would collapse. Field observations suggest that, at present, the ice sheet may be advancing into parts of the Ross Ice Shelf. Such advance, however, would not ensure the security of the ice sheet since ice streams that drain to the north appear to flow directly into the sea with little or no ice shelf to buttress them. If these ice streams do not flow over a sufficiently high bedrock sill then they provide the most likely avenues for ice-sheet retreat.

scholarly journals Century-scale simulations of the response of the West Antarctic Ice Sheet to a warming climate

2015 ◽  
Vol 9 (2) ◽  
pp. 1887-1942 ◽  
Author(s):  
S. L. Cornford ◽  
D. F. Martin ◽  
A. J. Payne ◽  
E. G. Ng ◽  
A. M. Le Brocq ◽  
...  
Keyword(s):  
Ice Sheet ◽  
Ice Shelf ◽  
Ice Streams ◽  
The West ◽  
Amundsen Sea ◽  
West Antarctic Ice Sheet ◽  
West Antarctic ◽  
Ice Accumulation ◽  
Grounding Line ◽  
Mesh Resolution

Abstract. We use the BISICLES adaptive mesh ice sheet model to carry out one, two, and three century simulations of the fast-flowing ice streams of the West Antarctic Ice Sheet. Each of the simulations begins with a geometry and velocity close to present day observations, and evolves according to variation in meteoric ice accumulation, ice shelf melting, and mesh resolution. Future changes in accumulation and melt rates range from no change, through anomalies computed by atmosphere and ocean models driven by the E1 and A1B emissions scenarios, to spatially uniform melt rates anomalies that remove most of the ice shelves over a few centuries. We find that variation in the resulting ice dynamics is dominated by the choice of initial conditions, ice shelf melt rate and mesh resolution, although ice accumulation affects the net change in volume above flotation to a similar degree. Given sufficient melt rates, we compute grounding line retreat over hundreds of kilometers in every major ice stream, but the ocean models do not predict such melt rates outside of the Amundsen Sea Embayment until after 2100. Sensitivity to mesh resolution is spurious, and we find that sub-kilometer resolution is needed along most regions of the grounding line to avoid systematic under-estimates of the retreat rate, although resolution requirements are more stringent in some regions – for example the Amundsen Sea Embayment – than others – such as the Möller and Institute ice streams.

Sign in / Sign up

Export Citation Format

Share Document