OATAO - Open Archive Toulouse Archive Ouverte Open Access Week

3-D numerical modelling of crustal polydiapirs with volume-of-fluid methods

Louis-Napoléon, Aurélie and Gerbault, Muriel and Bonometti, Thomas and Thieulot, Cédric and Martin, Roland and Vanderhaeghe, Olivier 3-D numerical modelling of crustal polydiapirs with volume-of-fluid methods. (2020) Geophysical Journal International, 222 (1). 474-506. ISSN 0956-540X

(Document in English)

PDF (Publisher's version) - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader

Official URL: https://doi.org/10.1093/gji/ggaa141


Gravitational instabilities exert a crucial role on the Earth dynamics and in particular on its differentiation. The Earth’s crust can be considered as a multilayered fluid with different densities and viscosities, which may become unstable in particular with variations in temperature. With the specific aim to quantify crustal scale polydiapiric instabilities, we test here two codes, JADIM and OpenFOAM, which use a volume-of-fluid (VOF) method without interface reconstruction, and compare them with the geodynamics community code ASPECT, which uses a tracking algorithm based on compositional fields. The VOF method is well-known to preserve strongly deforming interfaces. Both JADIM and OpenFOAM are first tested against documented two and three-layer Rayleigh–Taylor instability configurations in 2-D and 3-D. 2-D and 3-D results show diapiric growth rates that fit the analytical theory and are found to be slightly more accurate than those obtained with ASPECT. We subsequently compare the results from VOF simulations with previously published Rayleigh–B´enard analogue and numerical experiments.We show that the VOF method is a robust method adapted to the study of diapirism and convection in the Earth’s crust, although it is not computationally as fast as ASPECT. OpenFOAM is found to run faster than, and conserve mass as well as JADIM. Finally, we provide a preliminary application to the polydiapiric dynamics of the orogenic crust of Naxos Island (Greece) at about 16 Myr, and propose a two-stages scenario of convection and diapirism. The timing and dimensions of the modelled gravitational instabilities not only corroborate previous estimates of timing and dimensions associated to the dynamics of this hot crustal domain, but also bring preliminary insight on its rheological and tectonic contexts.

Item Type:Article
Additional Information:This article has been accepted for publication in Geophysical Journal International ©: 2020 The Author(s) Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved.
HAL Id:hal-02874978
Audience (journal):International peer-reviewed journal
Uncontrolled Keywords:
Institution:French research institutions > Centre National d'Études Spatiales - CNES (FRANCE)
French research institutions > Centre National de la Recherche Scientifique - CNRS (FRANCE)
Université de Toulouse > Institut National Polytechnique de Toulouse - Toulouse INP (FRANCE)
French research institutions > Institut de Recherche pour le Développement - IRD (FRANCE)
Université de Toulouse > Université Toulouse III - Paul Sabatier - UT3 (FRANCE)
French research institutions > Bureau de Recherches Géologiques et Minières - BRGM (FRANCE)
Other partners > Utrecht University - UU (NETHERLANDS)
Laboratory name:
Deposited On:04 Jun 2020 17:55

Repository Staff Only: item control page