29 May 2024
 | 29 May 2024
Status: this preprint is open for discussion.

Isogeometric analysis of the lithosphere under topographic loading: Igalith v1.0.0

Rozan Rosandi, Yudi Rosandi, and Bernd Simeon

Abstract. This paper presents methods from isogeometric finite element analysis for numerically solving problems in geoscience involving partial differential equations. In particular, we consider the numerical simulation of shells and plates in the context of isostasy. Earth's lithosphere is modeled as a thin elastic shell or plate floating on the asthenosphere and subject to topographic loads. We demonstrate the computational methods on the isostatic boundary value problem posed on selected geographic locations. For Europe, the computed lithospheric depression is compared with available Mohorovičić depth data. We also perform parameter identification for the effective elastic thickness of the lithosphere, the rock density, and the topographic load that are most plausible to explain the measured depths. An example of simulating the entire lithosphere of the Earth as a spherical shell using multi-patch isogeometric analysis is presented, which provides an alternative to spherical harmonics for solving partial differential equations on a spherical domain.

Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this preprint. The responsibility to include appropriate place names lies with the authors.
Rozan Rosandi, Yudi Rosandi, and Bernd Simeon

Status: open (until 24 Jul 2024)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
Rozan Rosandi, Yudi Rosandi, and Bernd Simeon
Rozan Rosandi, Yudi Rosandi, and Bernd Simeon


Total article views: 155 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
134 14 7 155 9 5
  • HTML: 134
  • PDF: 14
  • XML: 7
  • Total: 155
  • BibTeX: 9
  • EndNote: 5
Views and downloads (calculated since 29 May 2024)
Cumulative views and downloads (calculated since 29 May 2024)

Viewed (geographical distribution)

Total article views: 187 (including HTML, PDF, and XML) Thereof 187 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 12 Jun 2024
Short summary
We model Earth's lithosphere as a thin elastic shell and present numerical methods of isogeometric finite element analysis to simulate its deformation in isostatic equilibrium using technologies from computer-aided design. The simulations also serve as a basis for identifying parameters of the model that are most plausible to explain observed data. This research has been done to showcase the capabilities of isogeometric analysis in solving higher-order problems in geoscientific applications.