Submitted as: development and technical paper
24 Oct 2023
Submitted as: development and technical paper |  | 24 Oct 2023
Status: a revised version of 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.

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.
Hugo Dominguez, Nicolas Riel, and Pierre Lanari

Status: final response (author comments only)

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • RC1: 'Comment on gmd-2023-189', Marcin Dabrowski, 02 Apr 2024
    • AC1: 'Reply on RC1', Hugo Dominguez, 03 May 2024
  • RC2: 'Comment on gmd-2023-189', Albert de Montserrat Navarro, 05 Apr 2024
    • AC2: 'Reply on RC2', Hugo Dominguez, 03 May 2024
Hugo Dominguez, Nicolas Riel, and Pierre Lanari
Hugo Dominguez, Nicolas Riel, and Pierre Lanari


Total article views: 488 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
363 91 34 488 23 30
  • HTML: 363
  • PDF: 91
  • XML: 34
  • Total: 488
  • BibTeX: 23
  • EndNote: 30
Views and downloads (calculated since 24 Oct 2023)
Cumulative views and downloads (calculated since 24 Oct 2023)

Viewed (geographical distribution)

Total article views: 486 (including HTML, PDF, and XML) Thereof 486 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
Latest update: 08 May 2024
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.