ISMIP-HOM benchmark experiments using Underworld

[dansand]

-> submitter ORCID (or name)

0000-0002-2207-6837

-> slug

sachau-2022-icesheet

-> license

CC-BY-4.0

-> alternative license URL

No response

-> model category

community benchmark, forward model

-> model status

completed

-> associated publication DOI

https://doi.org/10.5194/gmd-15-8749-2022

-> model creators

0000-0002-3790-1385
0000-0002-8628-3704
Justin Lang
0000-0002-6469-3526
https://orcid.org/0000-0003-3685-174X

-> title

No response

-> description

Knowledge of the internal structures of the major continental ice sheets is improving, thanks to new investigative techniques. These structures are an essential indication of the flow behavior and dynamics of ice transport, which in turn is important for understanding the actual impact of the vast amounts of water trapped in continental ice sheets on global sea-level rise. The software studied here is specifically designed to simulate such structures and their evolution.

-> abstract

No response

-> scientific keywords

ice-sheet, benchmark, Stokes, anisotropy, mechanical

-> funder

https://ror.org/04s1m4564
https://ror.org/03a1kwz48

-> model embargo?

No response

-> include model code ?

• yes

-> model code/inputs DOI

https://doi.org/10.5281/zenodo.7384424

-> model code/inputs notes

No response

-> include model output data?

• yes

-> data creators

No response

-> model output data DOI

https://doi.org/10.5281/zenodo.7384424

-> model output data notes

No response

-> model output data size

2.9 Mb

-> software framework DOI/URI

https://doi.org/10.5281/zenodo.5935717

-> software framework source repository

No response

-> name of primary software framework (e.g. Underworld, ASPECT, Badlands, OpenFOAM)

No response

-> software framework authors

No response

-> software & algorithm keywords

Python, Finite-Element, Particle-in-cell

-> computer URI/DOI

No response

-> add landing page image and caption

Marker lines prior to (a) and after 750 years of flow of (b) isotropic and (c) anisotropic ice. The axial plane of the resulting shear fold in isotropic ice mimics the bedrock topography, while it is controlled by shearing along a horizontal shear zone in the case of anisotropic ice. Green: bedrock, flow to the right.

-> add an animation (if relevant)

No response

-> add a graphic abstract figure (if relevant)

Velocity field and strain rate field in isotropic ice (1) and inisotropic ice (2). Large strain rates and velocities occur in the vicinity of the bottleneck formed by the crest of the hill. Green: bedrock. For velocity, red is 70 m a−1 and blue is 0 m a−1. For strain rate, red is 0.032 $a^{−1}$, and blue is 0 $a^{−1}$.

-> add a model setup figure (if relevant)

(a) 2D geometry of Experiment B. This is identical to a section parallel X located at yˆ = 0.25 in Experiment A (right). Sloping angle α is given in degrees. Also depicted is the velocity field of the flowing ice, resulting for a model width L of 5000 m from the simulations described below. Color and arrow length visualize the amount of velocity. (b) Bedrock topography for Experiment A and general naming scheme for the axes of 3D experiments.

-> add a description of your model setup

No response

Please provide any feedback on the model submission process?

No response

[m-te-bot]

Model Report

Thank you for submitting.

Using Github actions, we have regenerated a report summarising information about your model

• Please check the report below, including the Errors and Warnings section

• You can update any information, by editing the markdown file at the top of the issue

• these edits will trigger the report will be regenerated

• once you are satisfied with the results, please add a review requested Review requested label

Parsed data

Section 1: Summary of your model

Model Submitter:

Dan Sandiford (0000-0002-2207-6837)

Model Creator(s):

• Till Sachau (0000-0002-3790-1385)
• Haibin Yang (0000-0002-8628-3704)
• Paul D. Bons (0000-0002-6469-3526)
• Louis-Noel Moresi (0000-0003-3685-174X)

Model slug:

sachau-2022-icesheet

(this will be the name of the model repository when created)

Model name:

ISMIP-HOM benchmark experiments using Underworld

License:

Creative Commons Attribution 4.0 International

Model Category:

• community benchmark
• forward model

Model Status:

• completed

Associated Publication title:

ISMIP-HOM benchmark experiments using Underworld

Short description:

Knowledge of the internal structures of the major continental ice sheets is improving, thanks to new investigative techniques. These structures are an essential indication of the flow behavior and dynamics of ice transport, which in turn is important for understanding the actual impact of the vast amounts of water trapped in continental ice sheets on global sea-level rise. The software studied here is specifically designed to simulate such structures and their evolution.

Abstract:

Abstract. Numerical models have become an indispensable tool for
understanding and predicting the flow of ice sheets and glaciers. Here we
present the full-Stokes software package Underworld to the glaciological
community. The code is already well established in simulating complex
geodynamic systems. Advantages for glaciology are that it provides a
full-Stokes solution for elastic–viscous–plastic materials and includes
mechanical anisotropy. Underworld uses a material point method to track the
full history information of Lagrangian material points, of stratigraphic
layers and of free surfaces. We show that Underworld successfully reproduces
the results of other full-Stokes models for the benchmark experiments of the Ice Sheet Model Intercomparison Project for Higher-Order Models
(ISMIP-HOM). Furthermore, we test finite-element meshes with different
geometries and highlight the need to be able to adapt the finite-element
grid to discontinuous interfaces between materials with strongly different
properties, such as the ice–bedrock boundary.

Scientific Keywords:

• ice-sheet
• benchmark
• Stokes
• anisotropy
• mechanical

Funder(s):

• AuScope (https://ror.org/04s1m4564)
• University of Tübingen (https://ror.org/03a1kwz48)

Section 2: your model code, output data

No embargo on model contents requested

Include model code:

True

Model code existing URL/DOI:

https://doi.org/10.5281/zenodo.7384424

Include model output data:

True

Model output data, existing URL/DOI:

https://doi.org/10.5281/zenodo.7384424

Section 3: software framework and compute details

Software Framework DOI/URL:

Found software: Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Name of primary software framework:

Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Software & algorithm keywords:

• Python
• Finite-Element
• Particle-in-cell

Section 4: web material (for mate.science)

Landing page image:

Filename: gmd-15-8749-2022-f09.png
Caption: Marker lines prior to (a) and after 750 years of flow of (b) isotropic and (c) anisotropic ice. The axial plane of the resulting shear fold in isotropic ice mimics the bedrock topography, while it is controlled by shearing along a horizontal shear zone in the case of anisotropic ice. Green: bedrock, flow to the right.

Animation:

Filename:

Graphic abstract:

Filename: graphic_abstract.png
Caption: Velocity field and strain rate field in isotropic ice (1) and inisotropic ice (2). Large strain rates and velocities occur in the vicinity of the bottleneck formed by the crest of the hill. Green: bedrock. For velocity, red is 70 m a−1 and blue is 0 m a−1. For strain rate, red is 0.032 $a^{−1}$, and blue is 0 $a^{−1}$.

Model setup figure:

Filename: gmd-15-8749-2022-f01-web.png
Caption: (a) 2D geometry of Experiment B. This is identical to a section parallel X located at yˆ = 0.25 in Experiment A (right). Sloping angle α is given in degrees. Also depicted is the velocity field of the flowing ice, resulting for a model width L of 5000 m from the simulations described below. Color and arrow length visualize the amount of velocity. (b) Bedrock topography for Experiment A and general naming scheme for the axes of 3D experiments.

Errors and Warnings

Associated Publication
Error fetching metadata with application/ld+json from https://api.crossref.org/works/https://doi.org/10.5194/gmd-15-8749-2022: 406 Client Error: Not Acceptable for url: https://api.crossref.org/works/https://doi.org/10.5194/gmd-15-8749-2022
Software Framework DOI/URI
doi.org metadata record succesfully extracted in json-ld format
Software Repository
Warning: no repository URL provided.
Submitter
ORCID metadata record succesfully extracted in json-ld format

Model creators
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format

• Error: name Justin Lang in unexpected format. Expected last name(s), first name(s) or ORCID.
  ORCID metadata record succesfully extracted in json-ld format
  ORCID metadata record succesfully extracted in json-ld format
  Could not parse Embargo date. Check format is
  Model code/inputs notes
  Warning: No notes provided.
  Model creators
  Error: no data creators found
  Model data notes
  Warning: No notes provided.
  Computer URI/DOI
  Warning: No URI/DOI provided.
  Animation
  Warning: No animation uploaded.

Model setup description
Warning: No description given

Next steps

• once the model_reviewers team has approved the model, we will create a repository for your model

[m-te-bot]

Review Requested

A review of this submission has been requested from @ModelAtlasofTheEarth/model_reviewers

[m-te-bot]

Model repository created at https://github.com/ModelAtlasofTheEarth/sachau-2022-icesheet