Numerical simulations of convection in crystal-bearing magmas: a case study of the magmatic system at Erebus, Antarctica

• Published in: Journal of Geophysical Research: Solid Earth
• Main Authors: Molina, Indira, Burgisser, Alain, Oppenheimer, Clive
• Other Authors: Institut des Sciences de la Terre d'Orléans - UMR7327 (ISTO), Bureau de Recherches Géologiques et Minières (BRGM)-Observatoire des Sciences de l'Univers en région Centre (OSUC), Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de Paris, Université Paris Sciences et Lettres (PSL)-Université Paris Sciences et Lettres (PSL)-Université d'Orléans (UO)-Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire de Paris, Université Paris Sciences et Lettres (PSL)-Université Paris Sciences et Lettres (PSL)-Université d'Orléans (UO)-Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Université d'Orléans (UO)-Centre National de la Recherche Scientifique (CNRS), University of Cambridge UK (CAM), ANT-0538414, ANT-0838817 awarded by the Office of Polar Programs (National Science Foundation)
• Format: Article in Journal/Newspaper
• Language: English
• Published: HAL CCSD 2012
• Subjects: [SDU.STU.VO]Sciences of the Universe [physics]/Earth Sciences/Volcanology; [SDE.MCG]Environmental Sciences/Global Changes; Antarc*; Antarctica
• Online Access: https://insu.hal.science/insu-00712931; https://insu.hal.science/insu-00712931/document; https://insu.hal.science/insu-00712931/file/Molina2012.pdf; https://doi.org/10.1029/2011JB008760

International audience The sustained heat and gas output from Erebus volcano reflects a regime of magma convection that we investigate here using a bi-phase (melt and crystals), fluid dynamical model. Following validity and verification tests of the model, we carried out four single-phase and three bi-phase numerical 30-year- simulations, in an idealized 2D geometry representing a lava lake cooled from above and a reservoir heated from below that are linked by a 4-to-10-m-diameter conduit. We tested the effects of crystals on convection while changing conduit size and the system boundaries from closed to open. Neglecting crystal settling yields only a limited number of features, i.e., (i) the formation of a central instability, (ii) the average temperature evolution, and (iii) the average velocity range of the surface flow motion. Biphase simulations show that while crystals are quite efficiently transported by the liquid phase a small decoupling reflecting their large size (5 cm) results in settling. This leads to more complex circulation patterns and enhances the vigor of fluid motion. A sufficiently large conduit sustains convection and retains 6 and 20% of crystals in suspension, for a closed and open system, respectively. Model outputs do not yet correspond well with field observations of Erebus lava lake (e.g., real surface velocities are much faster than those modeled), suggesting that exsolved volatiles are an important source of buoyancy.