Numerical Modelling of Ice Sheets, Streams, and Shelves

Bueler, Ed

doi:10.1007/978-3-030-42584-5_8

Cited by 3 publications

(4 citation statements)

References 46 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…From this point of view, a first important task that a DT has to achieve is to remove any possibility of misunderstanding for users though a rigorous definition of a reference glossary. A good example, explicitly aimed at making unskilled users properly conscious about glacier modelling, is the one from [75].…”

Section: Dt Theoretical Frameworkmentioning

confidence: 99%

See 1 more Smart Citation

Towards a Digital Twin Prototype of Alpine Glaciers: Proposal for a Possible Theoretical Framework

Fissore¹,

Bovio

Perotti

et al. 2023

Remote Sensing

View full text Add to dashboard Cite

The Destination Earth (DestinE) European initiative has recently brought into the scientific community the concept of the Digital Twin (DT) applied to Earth Sciences. Within 2030, a very high precision digital model of the Earth, continuously fed and powered by Earth Observation (EO) data, will provide as many digital replicas (DTs) as the different domains of the earth sciences are. Considering that a DT is driven by use cases, depending on the selected application, the provided information has to change. It follows that, to achieve a reliable representation of the selected use case, a reasonable and complete a priori definition of the needed elements that DT must contain is mandatory. In this work, we define a possible theoretical framework for a future DT of the Italian Alpine glaciers, trying to define and describe all those information (both EO and in situ data) and relationships that necessarily have to enter the process as building blocks of the DT itself. Two main aspects of glaciers were considered and investigated: (i) the “metric quantification” of their spatial dynamics (achieved through measures) and (ii) the “qualitative (semantic) description” of their health status as definable through observations from domain experts. After the first identification of the building blocks, the work proceeds focusing on existing EO data sources providing their essential elements, with specific focus on open access high-resolution (HR) and very-high-resolution (VHR) images. This last issue considered two scales of analysis: local (single glacier) and regional (Italian Alps). Some considerations were furtherly reported about the expected glaciers-related applications enabled by the availability of a DT at regional level. Applications involving both metric and semantic information were considered and grouped in three main clusters: Glaciers Evolution Modelling (GEM), 4D Multi Reality, and Virtual Reality. Limitations were additionally explored, mainly related to the technical characteristics of available EO VHR open data and some conclusions provided.

show abstract

Section: Dt Theoretical Frameworkmentioning

confidence: 99%

“…Area and changes of the Italian glaciers (time span 1957-2015) sorted according to the administrative Italian region where they are located. Data source[71][72][73][74][75][76][77][78][79][80].…”

mentioning

confidence: 99%

Towards a Digital Twin Prototype of Alpine Glaciers: Proposal for a Possible Theoretical Framework

Fissore¹,

Bovio

Perotti

et al. 2023

Remote Sensing

View full text Add to dashboard Cite

show abstract

“…An ad hoc approach to gaining the benefits of computational speed and validity in multiple flow regimes is to combine the SIA and SSA by summing their contributions to obtain the "hybrid" horizontal velocity field (hereafter called the Hybrid solver). This approach was proposed by Bueler (2009), who used a weighting function to transition between the two regimes. Later, Winkelmann et al (2011) recognized that, in light of the approximations involved and the uncertainty in basal friction, it is more straightforward simply to sum the contributions.…”

Section: Introductionmentioning

confidence: 99%

A comparison of the stability and performance of depth-integrated ice-dynamics solvers

2022

View full text Add to dashboard Cite

Abstract. In the last decade, the number of ice-sheet models has increased substantially, in line with the growth of the glaciological community. These models use solvers based on different approximations of ice dynamics. In particular, several depth-integrated dynamics solvers have emerged as fast solvers capable of resolving the relevant physics of ice sheets at the continental scale. However, the numerical stability of these schemes has not been studied systematically to evaluate their effectiveness in practice. Here we focus on three such solvers, the so-called Hybrid, L1L2-SIA and DIVA solvers, as well as the well-known SIA and SSA solvers as boundary cases. We investigate the numerical stability of these solvers as a function of grid resolution and the state of the ice sheet for an explicit time discretization scheme of the mass conservation step. Under simplified conditions with constant viscosity, the maximum stable time step of the Hybrid solver, like the SIA solver, has a quadratic dependence on grid resolution. In contrast, the DIVA solver has a maximum time step that is independent of resolution as the grid becomes increasingly refined, like the SSA solver. A simple 1D implementation of the L1L2-SIA solver indicates that it should behave similarly, but in practice, the complexity of its implementation appears to restrict its stability. In realistic simulations of the Greenland Ice Sheet with a nonlinear rheology, the DIVA and SSA solvers maintain superior numerical stability, while the SIA, Hybrid and L1L2-SIA solvers show markedly poorer performance. At a grid resolution of Δx=4 km, the DIVA solver runs approximately 20 times faster than the Hybrid and L1L2-SIA solvers as a result of a larger stable time step. Our analysis shows that as resolution increases, the ice-dynamics solver can act as a bottleneck to model performance. The DIVA solver emerges as a clear outlier in terms of both model performance and its representation of the ice-flow physics itself.

show abstract

“…An ad-hoc approach to gaining the benefits of computational speed and validity in multiple flow regimes is to combine the SIA and SSA by summing their contributions to obtain the "hybrid" horizontal velocity field (hereafter called the Hybrid solver). This approach was proposed by Bueler (2009), who used a weighting function to transition between the two regimes.…”

Section: Introductionmentioning

confidence: 99%

A comparison of the performance of depth-integrated ice-dynamics solvers

Robinson

Goldberg

Lipscomb

2021

Preprint

View full text Add to dashboard Cite

Abstract. In the last decade, the number of ice-sheet models has increased substantially, in line with the growth of the glaciological community. These models use solvers based on different approximations of ice dynamics. In particular, several depth-integrated dynamics approximations have emerged as fast solvers capable of resolving the relevant physics of ice sheets at the continen- tal scale. However, the numerical stability of these schemes has not been studied systematically to evaluate their effectiveness in practice. Here we focus on three such solvers, the so-called Hybrid, L1L2-SIA and DIVA solvers, as well as the well-known SIA and SSA solvers as boundary cases. We investigate the numerical stability of these solvers as a function of grid resolution and the state of the ice sheet. Under simplified conditions with constant viscosity, the maximum stable timestep of the Hybrid solver, like the SIA solver, has a quadratic dependence on grid resolution. In contrast, the DIVA solver has a maximum timestep that is independent of resolution, like the SSA solver. Analysis indicates that the L1L2-SIA solver should behave similarly, but in practice, the complexity of its implementation can make it difficult to maintain stability. In realistic simulations of the Greenland ice sheet with a non-linear rheology, the DIVA and SSA solvers maintain superior numerical stability, while the SIA, Hybrid and L1L2-SIA solvers show markedly poorer performance. At a grid resolution of ∆x = 4 km, the DIVA solver runs approximately 15 times faster than the Hybrid and L1L2-SIA solvers. Our analysis shows that as resolution increases, the ice-dynamics solver can act as a bottleneck to model performance. The DIVA solver emerges as a clear outlier in terms of both model performance and its representation of the ice-flow physics itself.

show abstract

Numerical Modelling of Ice Sheets, Streams, and Shelves

Cited by 3 publications

References 46 publications

Towards a Digital Twin Prototype of Alpine Glaciers: Proposal for a Possible Theoretical Framework

Towards a Digital Twin Prototype of Alpine Glaciers: Proposal for a Possible Theoretical Framework

A comparison of the stability and performance of depth-integrated ice-dynamics solvers

A comparison of the performance of depth-integrated ice-dynamics solvers

Contact Info

Product

Resources

About