Benchmark experiments for higher-order and full Stokes ice sheet models (ISMIP-HOM)

Pattyn, Frank; Perichon, Laura; Aschwanden, A.; Breuer, B.; Smedt, Brigitte De; Gagliardini, Olivier; Gudmundsson, G. Hilmar; Hindmarsh, Richard C. A.; Hubbard, Alun; Johnson, Jesse V.; Kleiner, Thomas; Konovalov, Y. V.; Martin, Charlotte; Payne, A. J.; Pollard, David; Price, Stephen; Rückamp, Martin; Saito, Fuyuki; Souček, Ondřej; Sugiyama, Shin; Zwinger, Thomas

doi:10.5194/tcd-2-111-2008

Cited by 62 publications

(101 citation statements)

References 32 publications

(33 reference statements)

Supporting

Mentioning

100

Contrasting

Order By: Relevance

“…as ice gradually filled the depression. Concave upward reflectors would be the result of large-scale flow, as is found with terrestrial ice sheets [Robin et al, 1969;Pattyn et al, 2008]. However, the reflector geometry in Korolev neither mimics hypothesized basal topography nor continues unperturbed laterally across the crater.…”

Section: 1002/2015gl066440mentioning

confidence: 92%

Three‐dimensional structure and origin of a 1.8 km thick ice dome within Korolev Crater, Mars

Holt

2016

Geophysical Research Letters

View full text Add to dashboard Cite

Korolev is an 80 km diameter impact crater located at 72.7°N, 164.5°E containing a large domed deposit. Perennial water ice is thought to be currently unstable there; however, using a 3‐D analysis of the dome's internal radar stratigraphy, we estimate that Korolev's central mound contains between 1400 and 3500 km3 of water ice that is up to 1.8 km thick. Furthermore, the stratigraphic structure of this ice dome is strikingly similar to the north polar layered deposits (NPLD) on Planum Boreum, approximately 600 km to the north. Additionally, our stratigraphic analysis suggests that Korolev ice was not previously part of a once larger polar ice sheet but rather it was deposited locally. We conclude that Korolev Crater's ice likely deposited during the same climate regime as Planum Boreum's NPLD, but independently. This implies that the incorporation of circumpolar ice deposits such as Korolev may prove useful in reconstructing a unique Mars polar climate history.

show abstract

Section: 1002/2015gl066440mentioning

confidence: 92%

Three‐dimensional structure and origin of a 1.8 km thick ice dome within Korolev Crater, Mars

Holt

2016

Geophysical Research Letters

View full text Add to dashboard Cite

show abstract

“…[11] Many higher-order approximations exist [e.g., Hindmarsh, 2004;Pattyn et al, 2008] that provide increased accuracy in situations where is larger. In most cases the higher-order approximations include the longitudinal stress components (membrane stresses [Hindmarsh, 2006]) and their spatial gradients, which seems sufficient in many "ice sheet scenarios" where the ice surface is only gently sloping [e.g., Blatter, 1995;Pattyn, 2002;Schoof and Hindmarsh, 2010;Goldberg, 2011].…”

Section: The Isosiamentioning

confidence: 99%

“…These new components call for careful validation. In this section, we thus make use of existing benchmark experiments [Pattyn et al, 2008] for comparing iSOSIA simulation results to similar predictions based on a wide range of ice sheet models, including full-Stokes formulations as well as multilayer and one-layer higher-order approximations.…”

Section: Validationmentioning

confidence: 99%

Modeling the flow of glaciers in steep terrains: The integrated second‐order shallow ice approximation (iSOSIA)

et al. 2011

View full text Add to dashboard Cite

[1] The rugged topography of mountain ranges represents a special challenge to computational ice sheet models simulating past or present glaciations. In mountainous regions, the topography steers glaciers through relatively narrow and steep valleys, and as a consequence hereof, the flow rate of alpine-style glaciers varies significantly at length scales comparable to that of the topography. Localized flow rate undulations generate longitudinal and transverse stress gradients within the ice, which, in turn, are of known importance to the flow itself, whether by internal deformation of the ice or by basal sliding, and to the interaction with topography through glacial erosion. Standard models capable of simulating mountain range-scale glaciations build on the so-called shallow ice approximation, which, among other parameters, neglects the longitudinal and transverse stress gradients, and therefore fails to capture the full effects of the rugged topography and related feedbacks between erosion by glacial sliding and the extent and style of glaciation. Here we present and test a new depth-integrated model framework which, on the one hand, takes into account the "higher-order" effects related to steep and rugged bed topography while it still, on the other hand, provides the computational efficiency needed for three-dimensional simulations of glaciation and landscape evolution in response to, for example, long-term climate variations. The model framework (integrated second-order shallow ice approximation (iSOSIA)) is based on depth integration of the second-order shallow ice approximation. Although depth integration is not guaranteed to maintain so-called second-order accuracy, the results of computational benchmark experiments show that iSOSIA, in spite of its efficient depth-integrated structure, performs well in situations with steep and rapidly undulating bed topography.

show abstract

“…Accurate representation of long-distance stress transmission requires the incorporation of membrane stresses to model the effects of ice-shelf or frontal changes on inland ice flow. Such stresses are now incorporated in some higher order large-scale ice-sheet models (Blatter 1995; Pattyn 2003; Pattyn et al 2008; Bueler & Brown 2009; Price et al 2011). One scheme for including these stresses in ice-sheet models, the vertically integrated membrane stress approximation (MSA), has proved useful (Kamb & Echelmeyer 1986; Muszynski & Birchfield 1987; MacAyeal 1989; Hindmarsh 2006 a ).…”

Section: Introductionmentioning

confidence: 99%

Frequency response of ice streams

2012

View full text Add to dashboard Cite

Changes at the grounding line of ice streams have consequences for inland ice dynamics and hence sea level. Despite substantial evidence documenting upstream propagation of frontal change, the mechanisms by which these changes are transmitted inland are not well understood. In this vein, the frequency response of an idealized ice stream to periodic forcing in the downstream strain rate is examined for basally and laterally resisted ice streams using a one-dimensional, linearized membrane stress approximation. This reveals two distinct behavioural branches, which we find to correspond to different mechanisms of upstream velocity and thickness propagation, depending on the forcing frequency. At low frequencies (centennial to millennial periods), slope and thickness covary hundreds of kilometres inland, and the shallow-ice approximation is sufficient to explain upstream propagation, which occurs through changes in grounding-line flow and geometry. At high frequencies (decadal to sub-decadal periods), penetration distances are tens of kilometres; while velocity adjusts rapidly to such forcing, thickness varies little and upstream propagation occurs through the direct transmission of membrane stresses. Propagation properties vary significantly between 29 Antarctic ice streams considered. A square-wave function in frontal stress is explored by summing frequency solutions, simulating some aspects of the dynamical response to sudden ice-shelf change.

show abstract

Benchmark experiments for higher-order and full Stokes ice sheet models (ISMIP-HOM)

Cited by 62 publications

References 32 publications

Three‐dimensional structure and origin of a 1.8 km thick ice dome within Korolev Crater, Mars

Three‐dimensional structure and origin of a 1.8 km thick ice dome within Korolev Crater, Mars

Modeling the flow of glaciers in steep terrains: The integrated second‐order shallow ice approximation (iSOSIA)

Frequency response of ice streams

Contact Info

Product

Resources

About