Modèles numériques de coulées de lave pour une estimation et une délimitation du risque volcanique, le cas du volcan El Reventador (Equateur)

Lava flows are the most representative volcanic products of effusive eruptions and are formed whenthe magma is extruded and flows on the surface. When lava flows reach the surface they lose heat and cool.Cooling affects directly the rheology of the lava up to a point where it cannot flow anymore. Rheologicalparameters that control the dynamics of lava flows are the viscosity and the yield strength which in turndepends on the chemical composition, crystallinity and bubble content. There exist numerous models forthe rheology estimation, mostly developed for basaltic lava flows and few for andesitic ones.Lava flows can highly affect populated areas, infrastructures and environment. A way to forecastthe future damages is to developed numerical codes of the lava propagation on real volcanic topography.This challenging method combines the topography, the rheology, the heat loss, and flow dynamics tosimulate the emplacement of a particular lava flow. The numerical code VolcFlow which is based on thedepth-averaged approach is able to reproduce the main physical characteristics of the deposits likemorphology, length and thickness. Here 3 models are proposed for their implementation in VolcFlow withthe aim to simulate lava flows. One model is isothermal, the second includes cooling and the associatedrheological variations, and the third takes into account the crust formation and its effect on the flowemplacement. To check the validity of the different approaches, the models were tested with four studycases, two with basaltic compositions (molten basalt experiment of the Syracuse lava Project and the August-November, 2015 lava flow from Piton de la Fournaise, France) and two with andesitic compositions (theDecember 4th-5th lava flow from Tungurahua, Ecuador, and three lava flows from El Reventador,Ecuador). Results of the simulations shows that the isothermal model can reproduce the flows even if itdoes not consider the cooling and rheology variation. The model that includes rheological laws as functionof crystallization induced by cooling down flow can give very good results but is very sensitive to the inputdata, in particular to the fluid viscosity that is very dependent on chemical composition and temperature.Finally, the model that includes cooling and synthetic sigmoid rheological law shows good coherence for allthe cases except at Piton de la Fournaise. The model that aims to simulate the formation of a crust on thelava flow surface, lava flowing underneath and break-out mechanisms leads to the thickening of the crust.Hence, break-out mechanism is not reproduced with VolcFlow.Les coulées de laves sont les produits volcaniques les plus représentatifs des éruptions effusives. Elles sont formées quand le magma est extrudé et se répand à la surface de la Terre. Quand la lave arrive en surface, elle perd de la chaleur et refroidit. Le refroidissement affecte directement les propriétés rhéologiques de la lave, jusqu’à arrêter son écoulement. Les paramètres rhéologiques qui contrôlent la dynamique des coulées de laves sont la viscosité et le seuil de plasticité, qui dépendent eux-mêmes de la composition chimique, de la cristallinité et de la teneur en bulles. Il existe de nombreux modèles d’estimation de la rhéologie, la plupart développés pour les coulées de lave basaltiques et quelque uns pour les coulées de lave andésitiques. Les coulées de laves peuvent grandement affecter les régions peuplées, les infrastructures et l’environnement. Un moyen de prévoir les futurs dégâts est de développer des modèles numériques pour prévoir la propagation des coulées de laves sur des topographies volcaniques réelles. Cette méthode difficile combine la topographie, la rhéologie, la perte de chaleur et la dynamique de l’écoulement pour simuler l’emplacement d’une coulée de lave précise. Le code numérique VolcFlow, qui est basé sur une approche moyennée verticale, est capable de reproduire les caractéristiques principales des dépôts comme la morphologie, la longueur et l’épaisseur. Dans cette étude sont proposés trois modèles implémentés dans VolcFlow et ayant pour but de simuler des coulées de laves. Le premier est isotherme, le deuxième inclut le refroidissement et les variations rhéologiques associées, et le troisième prend en considération la déformation de la croûte à la surface de la coulée et son effet sur l’emplacement de la coulée. Afin de vérifier la validité des différentes approches, les modèles sont testés sur quatre cas d’étude : deux coulées de lave de composition basaltique (expérience de basalte fondu de Syracuse lava Project et la coulée de lave d’août novembre 2015 du Piton de la Fournaise, France) et deux de compositions andésitique (la coulée de lave du 4-5 décembre 2015 du Tungurahua et trois coulées de lave du Reventador, Equateur). Les résultats des simulations montrent que le modèle isotherme peut reproduire les coulées même s’il ne prend pas en compte les variations de rhéologie et le refroidissement. Le modèle incluant la cristallisation, induite par le refroidissement de la lave au cours de son écoulement, et les variations rhéologiques associées donne de très bons résultats mais est très sensible aux paramètres d’entrée, en particulier à la viscosité, elle-même très dépendante de la composition chimique et de la température. Enfin, le modèle prenant en compte le refroidissement et les variations de rhéologie par une loi synthétique sigmoïde montre une bonne cohérence dans tous les cas simulés, sauf pour le Piton de la Fournaise. Le modèle visant à simuler la formation d’une croûte à la surface de la lave et sa percée par l’écoulement sous-jacent amène uniquement à l’épaississement de la croûte. Le mécanisme de percée n’est pas reproduit avec VolcFlow