Submitted as: development and technical paper
24 Oct 2023
Submitted as: development and technical paper |  | 24 Oct 2023
Status: this preprint is currently under review for the journal GMD.

Modelling chemical advection during magma ascent

Hugo Dominguez, Nicolas Riel, and Pierre Lanari

Abstract. Modelling magma transport requires robust numerical schemes for chemical advection. Current numerical schemes vary in their ability to be mass conservative, computationally efficient, and accurate. This study compares four of the most commonly used numerical schemes for advection: an upwind scheme, a weighted essentially non-oscillatory (WENO-5) scheme, a semi-Lagrangian (SL) scheme, and a marker-in-cell (MIC) method. We assess the behaviour of these schemes using the passive advection of two different magmatic compositions. This is coupled in 2D with the temporal evolution of a melt anomaly that generates porosity waves. All algorithms, except the upwind scheme, are able to predict the melt composition with reasonable accuracy. In terms of total running time, the upwind and SL schemes are the fastest, and the MIC scheme is the slowest. The WENO-5 scheme shows intermediate total running time but has the lowest amount of mass loss and therefore is best suited for this problem.

Hugo Dominguez et al.

Status: open (until 19 Dec 2023)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse

Hugo Dominguez et al.

Hugo Dominguez et al.


Total article views: 224 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
173 44 7 224 4 7
  • HTML: 173
  • PDF: 44
  • XML: 7
  • Total: 224
  • BibTeX: 4
  • EndNote: 7
Views and downloads (calculated since 24 Oct 2023)
Cumulative views and downloads (calculated since 24 Oct 2023)

Viewed (geographical distribution)

Total article views: 221 (including HTML, PDF, and XML) Thereof 221 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 06 Dec 2023
Short summary
Predicting the behaviour of magmatic systems is important for understanding Earth's matter and heat transport. Numerical modelling is a technique that can predict complex systems at different scales of space and time, by solving equations using various techniques. This study tests four algorithms to find the best way to transport the melt composition. The "Weighted Essentially Non-Oscillatory" algorithm emerges as the best choice, minimizing errors and preserving system mass.