Preprints
https://doi.org/10.5194/se-2022-7
https://doi.org/10.5194/se-2022-7
 
28 Jan 2022
28 Jan 2022
Status: a revised version of this preprint is currently under review for the journal SE.

An efficient parallel method to compute lithostatic pressure in thermo-mechanical geodynamic models

Anthony Jourdon and Dave A. May Anthony Jourdon and Dave A. May
  • Institute of Geophysics and Planetary Physics, Scripps Institution of Oceanography, UC San Diego, La Jolla, CA

Abstract. Modelling the lithostatic pressure in the Earth's interior is a common problem in Earth sciences. In this study we propose to compute the lithostatic pressure from the conservation of momentum reduced to the hydrostatic particular case. It results a partial differential equation that can be solved using the classical numerical methods. To show the usefulness of solving a PDE to compute the lithostatic pressure we propose two 2D models, one with a deformed mesh and one with a radial gravity acceleration vector and a concentric density distribution. Moreover, we also present a 3D rift model using the lithostatic pressure as a boundary condition. This model shows a high non cylindricity resulting from the Neumann boundary condition that is accomodated by strike-slip shear zones. We compare the result of this numerical model with a simple free-slip boundary conditions model to demonstrate the first order implications of considering "open" boundary conditions in 3D thermo-mechanical models.

Anthony Jourdon and Dave A. May

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on se-2022-7', Cedric THIEULOT, 15 Feb 2022
  • RC2: 'Comment on se-2022-7', Rene Gassmoeller, 07 Mar 2022

Anthony Jourdon and Dave A. May

Anthony Jourdon and Dave A. May

Viewed

Total article views: 373 (including HTML, PDF, and XML)
HTML PDF XML Total Supplement BibTeX EndNote
277 82 14 373 27 4 3
  • HTML: 277
  • PDF: 82
  • XML: 14
  • Total: 373
  • Supplement: 27
  • BibTeX: 4
  • EndNote: 3
Views and downloads (calculated since 28 Jan 2022)
Cumulative views and downloads (calculated since 28 Jan 2022)

Viewed (geographical distribution)

Total article views: 362 (including HTML, PDF, and XML) Thereof 362 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 18 May 2022
Download
Short summary
In this study we presented a method to compute the lithostatic pressure in which we cast the problem in terms of a partial differential equation (PDE). We showed in the context of 3D models of continental rifting, using the lithostatic pressure as a boundary condition within the flow problem, resulted in non-cylindrical velocity fields which produced strain localization in the lithosphere along large scale strike-slip shear zones and the formation and evolution of triple junctions.