The analog model, the numerical model and the Piton de la Fournaise : tale of a propagating dike
Authors: Furst Séverine , Pinel Virginie , Maccaferri Francesco 
Published: March 2022

&lt;p&gt;The transport of magma through the crust may sometimes result in volcanic eruptions at the surface, feeding a central conduit, opening fissures on the volcano flank, or at new locations in a volcanic field. Magma travels in the brittle crust by opening its way through the surrounding rock. In addition to the fracturation of the medium, the process of diking is also controlled by the magma flow, fluid-gas phase transitions, and the heat exchange. Representing the propagation of magmatic intrusions using analog and numerical model is essential to understand the physical processes occurring in nature and to mitigate the volcanic hazard linked to the emplacement of magmatic intrusions.&lt;/p&gt;&lt;p&gt;&lt;span&gt;In this context, w&lt;/span&gt;&lt;span&gt;e performed analog experiments of &lt;/span&gt;&lt;span&gt;air and &lt;/span&gt;&lt;span&gt;silicon oil injections in a solidified gelatin block. &lt;/span&gt;&lt;span&gt;Using three cameras, we monitored the propagation of the &lt;/span&gt;&lt;span&gt;oil-filled cracks &lt;/span&gt;&lt;span&gt;from the front, side and top views of the tank. The processing of time lapsed pictures enables to access the crack shape (dimension and orientation), trajectory and velocity. This analog modeling technique is routinely used to simulate magmatic dike propagation in the crust. &lt;/span&gt;&lt;span&gt;Then, taking advantage of these &lt;/span&gt;&lt;span&gt;well constrained experiments, we could validate &lt;/span&gt;&lt;span&gt;a novel 2D boundary element model for crack propagation &lt;/span&gt;&lt;span&gt;coupling brittle-elastic and fluid-dynamic equations&lt;/span&gt;&lt;span&gt;. &lt;/span&gt;&lt;span&gt;To do s&lt;/span&gt;&lt;span&gt;o, we initiate the input and boundary conditions of our numerical simulations, using &lt;/span&gt;&lt;span&gt;gelatin &lt;/span&gt;&lt;span&gt;and &lt;/span&gt;&lt;span&gt;oil &lt;/span&gt;&lt;span&gt;parameters &lt;/span&gt;&lt;span&gt;from the analog experiments&lt;/span&gt;&lt;span&gt;. &lt;/span&gt;&lt;span&gt;The ou&lt;/span&gt;&lt;span&gt;tputs &lt;/span&gt;&lt;span&gt;of the model include &lt;/span&gt;&lt;span&gt;the &lt;/span&gt;&lt;span&gt;crack shape, &lt;/span&gt;&lt;span&gt;trajectory, and&lt;/span&gt;&lt;span&gt; velocity, &lt;/span&gt;&lt;span&gt;that&lt;/span&gt;&lt;span&gt; is computed according to an energy conservation equation, under the assumption that fluid viscous forces are limiting the crack propagation velocity.&lt;/span&gt; &lt;span&gt;Numerical simulations are&lt;/span&gt;&lt;span&gt; faced with&lt;/span&gt;&lt;span&gt; the observations from &lt;/span&gt;&lt;span&gt;our&lt;/span&gt; &lt;span&gt;air and oil-filled crack propagation&lt;/span&gt;&lt;span&gt; experiments. &lt;/span&gt;&lt;span&gt;E&lt;/span&gt;&lt;span&gt;ventually, we applied the numerical model to the 1998 magmatic intrusion at Piton de la Fournaise volcano (La R&amp;#233;union Island), confronting the timing of the of the propagation with the migration of volcano-tectonic events.&lt;/span&gt;&lt;/p&gt;