No description available.

Dans un contexte de perturbations du système Terre, il est important d'améliorer notre compréhension des cycles géochimiques déterminant la composition chimique de l'atmosphère. Les aérosols de sulfate sont directement émis par les sources (naturelle ou anthropique) ou bien produit par des gaz soufrés précurseurs. Ils participent à l'augmentation de la pollution urbaine impliquant des risques sanitaires élevés et à l'augmentation des particules atmosphériques qui détermine la chimie atmosphérique, le bilan radiatif terrestre et in fine l'évolution du climat. Pour ces raisons, le cycle du soufre est un enjeu scientifique fort. Les fractionnements isotopiques dépendant et indépendant de la masse de l'oxygène et du soufre sont des traceurs géochimiques performants qui permettent d'étudier quantitativement l'origine et le devenir des sulfates. Notre étude développe en premier lieu un outil analytique de pré traitement des échantillons. Ensuite, nous avons montré que les sulfates primaires d'origine volcanique (via une étude au Stromboli) ne présentent pas d'anomalie isotopique en soufre et en oxygène. Des résultats similaires sont observés pour les sulfates adsorbés sur les cendres fraiches du Popocatépetl (Mexique) indiquant des processus de formation primaires. Les sulfates de la ville de Mexico présentent des anomalies isotopiques de l'ordre de D17O=0.9 et D33S=0.3. Les valeurs en D17O sont interprétées par des processus d'oxydation secondaires via les oxydants atmosphérique OH, H2O2, O3 et O2-TMI. Les valeurs en D33S sont encore peu contraintes. Notre étude montre finalement une influence des émissions volcaniques sur le bilan global des sulfates urbains à Mexico City. / Currently in our environmental context it is important to better constraint the geochemical cycles that are a key factor for atmospheric composition and climate change. Sulfate aerosols are emitted directly by sources (anthropogenic or natural) or produced from sulfur precursor gases. In the atmosphere, sulfate aerosols are associated with the high level of urban pollutant that impact human health and with the increase of atmospheric particle that have radiative effect. For these reasons studies of sulfate aerosols from anthropogenic and volcanic emissions are an important scientific stake. First, our study develops a method to extract and purify sulfates from natural samples which is compulsory for ambiguous isotope analyses. Then, from a field project at Stromboli, we show that primary sulfate aerosols from volcano are mass dependent (i.e. D17O=0‰). Similar results are observed for sulfates adsorbed on fresh ash from Popocatépetl volcano (Mexico) suggesting primary processes during volcanic eruptions. For urban sulfate aerosols we observe significant oxygen and sulfur mass independent fractionation (MIF) (D17O=0.9‰ and D33S=0.3‰) which are interpreted as secondary oxidation processes from atmospheric oxidant (OH, H2O2, O3, O2-TMI) for O-MIF and perhaps stratospheric input for S-MIF but the quantitative processes remained poorly constraint. Finally we discuss the influence of volcanic emissions on the urban area of Mexico City.

El presente trabajo de investigación tiene como objetivo determinar el valor fundamental de Volcan Compañía Minera S.A.A. Para calcularlo, se ha empleado el método de flujos de caja descontados, al costo promedio ponderado de capital (WACC, por sus siglas en inglés). Se efectuó una proyección de los estados financieros de la empresa por un horizonte de cuarenta años (2020 a 2059), para lo cual se utilizaron distintos supuestos respecto a producción, precios, inversiones, costos operativos y fuentes de financiamiento, entre otros aspectos. A través de esta metodología, se estimó un valor del patrimonio de $745.0 millones y un valor de S/ 0.44 por acción tipo B, lo cual representa un downside del 7.8 % respecto a la cotización al 30 de septiembre de 2019 de S/ 0.48. Asimismo, se estimó un valor de S/ 0.89 por acción de tipo A, con base en un valor del voto estimado de 2.03x. Luego de realizar el respectivo análisis de riesgos se encontró que la volatilidad del valor intrínseco de la acción es alta, principalmente debido a la alta sensibilidad del valor de la compañía a los precios del zinc y la plata.

Un des enjeux de la volcanologie réside dans la compréhension du fonctionnement du système de stockage et de transferts magmatiques associés aux éruptions. Au Piton de La Fournaise, cette étude est ici abordée par une analyse croisée des déformations, couplée à des modélisations numériques, et aux données géochimiques de la période 1998-2007. Au sein de cette activité, nous définissons, depuis 2000, un caractère cyclique ; chaque cycle est défini par une séquence d'éruptions sommitales et latérales terminées par une éruption distale à océanite. Ces éruptions successives sont caractérisées par des dykes s'initiant à des niveaux de plus en plus profonds et des laves de plus en plus primitives. Ces cycles durant lesquels l'inflation du cône est continue, pourraient témoigner d'une réalimentation continue du système de stockage superficiel. Ces résultats permettront dans le futur une meilleure prédiction de la localisation des éruptions, selon leur position dans un cycle d'activité.

Les études des déformations de surface en relation avec l'activité volcanique permettent de quantifier les phénomènes de transfert de magma qui s'opèrent dans les structures superficielles et profondes d'un édifice volcanique. Ces études s'appuient essentiellement sur l'utilisation de séries temporelles acquises par des réseaux de récepteurs GNSS installés sur les flancs de l'édifice volcanique et sur l'utilisation d'images acquises par des satellites équipés de capteurs à ouverture de synthèse. Les objectifs de ce travail ont été de mettre en œuvre sur deux des volcans les plus actifs du monde des méthodes numériques pour détecter, analyser et interpréter les déformations du sol associées à l'activité. Sur le Piton de la Fournaise, nous avons analysé l'évolution spatiale et temporelle du champ de déplacement entre l'éruption historique d'avril 2007 et octobre 2014 à partir de l'analyse de mesures continues acquises par les stations GNSS et de longues séries temporelles d'images radar Cosmos-Skymed et TerraSAR acquises en bande X. Pour le traitement des données radars, nous avons adopté une approche classique qui exploite la redondance d'information dans les interférogrammes et nous avons mis en œuvre une méthode originale de correction des effets troposphériques reposant sur la décomposition des signaux radars en valeurs singulières. La complexité spatiale et temporelle du champ de déplacement obtenu indique qu'une partie importante de l'édifice volcanique est affectée par des déformations d'origines diverses qui se superposent spatialement et temporellement. Ainsi, on observe des processus de subsidence qui ne s'accompagnent pas de déplacements horizontaux sur les coulées de lave récentes. Nous montrons qu'il existe une relation linéaire entre cette subsidence et l'épaisseur de la coulée et que son amplitude décroit avec le temps. Ces relations nous permettent de construire une loi empirique pour estimer la contribution de ce processus dans le champ de déformation. Nous observons également que le cône central subside de manière persistante durant la période étudiée. Nous interprétons cette subsidence comme l'expression d'une relaxation des contraintes provoquée par l'effondrement de plus de 350 m du Dolomieu survenu lors de l'éruption d'avril 2007. Enfin, nous montrons qu'une large partie du flanc est de l'édifice volcanique est affectée d'un mouvement lent le long de la pente entre 2007 et 2014. L'absence d'évidences sur la structure et sur la rhéologie de l'édifice nous amène à explorer différentes hypothèses pour expliquer l'origine de ce glissement qui pourrait être contrôlé par les propriétés frictionnelles des roches le long d'un ou de plusieurs plans de faille, ou bien être l'expression d'une déformation ductile dépendante de la viscosité du milieu. Le Llaima est un large strato-volcans andin dont les processus de déformations sont mal compris à cause principalement de la complexité de son mode de fonctionnement mais, également, aussi par l'absence de réseaux d'observation au sol. Dans ce contexte, les potentialités des données radar pour le suivi des déformations de surface de ces volcans constituent un intérêt scientifique majeur. / We address in this dissertation the use of Interferometric Synthetic Aperture Radar (InSAR) to measure and characterize the ground surface deformation at two volcanoes - Piton de la Fournaise (La Réunion Island, France) and Llaima (Chile). For Piton de la Fournaise, we analyzed the spatial pattern and temporal evolution of the ground displacement between the historical March-April 2007 eruption and October 2014, based on continuous measurements recorded by GNSS stations and X band COSMO-SkyMed and TerraSAR-X/TanDEM-X time series analysis. For the processing of radar data, we adopted a classical InSAR time series approach that exploits the information redundancy in the interferograms and we implemented an original method for correcting artifacts based on the principal component decomposition. The spatial and temporal complexity of the obtained deformation field indicates that an important part of the volcanic edifice is affected by deformations of various origins that overlap spatially and temporally. We observe also subsidence processes that are not accompanied by horizontal displacements in recent lava fields. We show that there exists a linear relationship between the subsidence and the thickness of lava and that the amplitude of subsidence decreases with time. These relationships allow us to construct an empirical law to estimate the contribution of post-lava emplacement process in the deformation field. We also observe that the Central Cone subsides persistently during the study period. We interpret this subsidence as the expression of a relaxation of the stresses caused by the Dolomieu collapse during the March-April 2007 eruption. Finally, we show that a widespread time-dependent moving sector on the Eastern Flank is affected by downslope motion during the 2007-2014 period. The uncertainties on both the structure and rheology parameters of the edifice leads us to explore different hypotheses to explain the origin of this flank motion which could be controlled by the frictional properties of the rocks along one or more fault planes, or be the expression of a dependent ductile deformation of the viscosity of the medium. Llaima is a large Andean stratospheric volcano whose deformation processes are poorly understood not only because of the complexity of its functioning mode but also because of the absence of observation networks on the ground. In this context, the potential of radar data for monitoring the ground deformations of these volcanoes is a main scientific interest. However, the complex environment conditions (steep slopes, snow- or ice-capped summit, dense vegetation cover, and strong tropospheric artifacts) and limited amount of available radar data make it challenging to accurately measure ground displacement with InSAR.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Le Merapi est un volcan qui a une actvité quasi-continue, L'explosion est normalement accompagnée par une série de nuées ardentes,edites merapienncs. La chronologie d'activité volcanique varie d'une activité à l'autre, il est donc difficile de trouver le niveau d'activité sismique critique pour ce volcan, En principe, il y a 2 types d'activité: avec et sans séismes volcanotectoniques de types VTA et VTB. Les données de séismes de types VTA et VTB récoltées au courant de l'année 1991 qui ont servi à cette etude, ont fait l'objet d'un traitement de routine (dépouillement, localisation de l'hypocentre, ... ), Ce traitement nous a conduit à remarquer que les séismes de type VTA et ceux de VTB sont séparés par une zone asismique à environ 1.5 km de profondeur. Ce résultat et ceIui de l'analyse pétrographique (Berthommier et al., 1992) conduisent à l'hypothèse de l'existence d'une poche magmatique à cette profondeur. Parmi les séismes de type VTB, nous avons trouvé certains de forme similaire, que l'on appelle doublets ou multiplets. En analysant la variation du délai le long du sismogramme, en utilisant la méthode inter-spectrale de la fenêtre mobile, nous avons constaté que la vitesse sismique a augmenté d'environ 1, 7 % de janvier à septembre 1991. Cette valeur est de l'ordre de dix fois celle observée sur la région tectonique. L'augmentation de vitesse est générée par la croissance de la pression du magma avant l'éruption. La surveillance de la formation du dôme de 1994 a montré que, au bout d'un moment, le processus s'arrête. Nous supposons que ceci est dû à l'équilibre qui s'établit entre la pression du magma à l'extrémité du conduit et la pression lithostatique que génère le dôme, lui-même. Si cel équilibre quasi-stable est interrompu,l'explosion se produit comme celle du 22 novembre 1994.

En 2004 et 2007, neuf structures d'échappements de fluides (SEF), constituées de volcans de boue (VdB), d'évents de gaz (GS, gas seeps en anglais) et/ou de sources, ont été échantillonnées dans la partie émergée du Bassin de la Côte Est (BCE) de l'île nord de la Nouvelle-Zélande.L'étude granulométrique indique que la boue émise par les VdB, les roches encaissantes et les niveaux de décollement voisins sont composés d'argiles et de silts en majorité. L'étude de la fraction argileuse et de la roche totale par diffraction des rayons X (DRX) a révélé de nombreuses similitudes. Les volcans de boue, les roches encaissantes et les niveaux de décollement présentent ainsi le même assemblage minéralogique, à savoir smectite, illite, chlorite, kaolinite, quartz et feldspaths. Les proportions sont variables d'un échantillon à un autre mais le couple smectite-illite est toujours majoritaire.L'étude géochimique de la fraction solide indique que les échantillons sont riches en Si02, pauvres en Fe2O3, MgO, MnO et en alcalins, à quelques exceptions près. La composition en éléments majeurs s'organise entre un pôle argileux alcalin et un pôle carbonaté. Les spectres de terres rares sont similaires et caractéristiques des argiles ; Ils présentent également un faible degré de fractionnement, lié à la formation des carbonates. Ce sont donc les mêmes minéraux qui contrôlent la chimie des échantillons.L'étude géochimique de la phase liquide montre que l'eau impliquée dans les volcans de boue est d'origine marine essentiellement, et des réactions eau/roche similaires, notamment l'altération de smectite en illite. Cette étude a permis aussi d'obtenir une estimation de la température d'équilibre, comprise entre 60 et 110°C, impliquant une profondeur d'origine de 2 à 3 km, voire plus.L'étude géophysique indique qu'à l'aplomb des VdB et des deux sources chaudes étudiés, la croûte continentale a sensiblement la même épaisseur et que la profondeur de la croûte océanique en subduction avoisine les 20 km. A cette profondeur, la fusion de la péridotite n'est pas possible et la fusion résultante de la croûte continentale, responsable du volcanisme d'arc, non plus. Le gradient géothermique mesuré à TePuia est donc influencé par un autre phénomène, mais la modélisation de la lithosphère ne nous a pas permis de trouver lequel.Ces différentes études mettent en évidence des caractéristiques géochimiques, pétrographiques et minéralogiques communes aux volcans de boue de Nouvelle-Zélande. Les fluides impliqués dans ces structures proviendraient donc d'un même niveau source, recouvert du même assemblage sédimentaire. L'étude géophysique ne nous apporte aucune information à ce sujet mais permet cependant d'établir avec certitude que le régime thermique est le même du Nord de Hawke's Bay au Sud du BCE ; la région de TePuia est un cas particulier, peut-être influencée par le complexe volcanique de Matakaoa.

Sur les volcans on enregistre fréquemment des accélérations du nombre cumulé de séismes et parfois du déplacement de la surface dans les heures, les jours ou les semaines précédant les éruptions. Expliquer une accélération du déplacement de surface par l'utilisation de modèles élastiques linéaires pour les édifices implique d'introduire une accélération de la pressurisation du réservoir et de la base du conduit magmatique avec des temps caractéristiques courts, ce qui est peu réaliste. Une autre voie de recherche consiste à considérer qu'une éruption est due à la rupture du réservoir magmatique, et que cette rupture n'est pas instantanée - elle peut être précédée par une phase d'endommagement progressif de l'édifice, pendant laquelle de petites ruptures ont lieu dans l'édifice. Ces ruptures provoquent l'affaiblissement des modules élastiques de l'édifice. Nous avons examiné l'effet que peut avoir l'affaiblissement des modules élastiques sur le déplacement de surface, ainsi que sur la pression dans le réservoir et le débit de magma entrant. Pour cela nous avons d'abord recherché une loi d'endommagement qui permette d'ajuster les déplacements de surface enregistrés en utilisant le nombre cumulé de séismes. La loi trouvée est une loi puissance d'un paramètre d'endommagement incrémental, qui est le rapport entre la longueur de rupture incrémentale caractéristique et la longueur à rompre pour obtenir une éruption. Cette loi fait intervenir, à l'exposant, le nombre cumulé de séismes.Ce modèle permet d'expliquer des vitesses et des accélérations du déplacement aussi bien faibles que fortes, selon la valeur prise par le paramètre d'endommagement, ainsi que la bifurcation entre un état stable de l'édifice, et un état instable qui mène à l'éruption. L'affaiblissement de l'édifice induit une diminution de la surpression et une accélération du déplacement de surface, la surpression diminuant plus lentement que les modules élastiques du fait de l'approvisionnement du réservoir en magma.Ce modèle peut expliquer une large gamme de comportements pré-éruptifs et éruptifs dans les volcans. Il permet d'ajuster simplement aussi bien les accélérations du déplacement de surface enregistré par les stations GPS sommitales du Piton de la Fournaise pour l'éruption du 30/03/2007, que les déplacements croissants de façon stationnaire du volcan Grimsvötn (Islande), entre 2005 et 2011. Il peut aussi être utilisé pour comprendre les comportements éruptifs plus complexes, comme dans le cas de l'éruption de 2010 du Mérapi (Indonésie). La valeur de la constante de temps du système d'approvisionnement contrôle le décalage temporel entre le maximum de l'endommagement (fin de la crise sismique volcano-tectonique) et le maximum de la déformation ou du flux de magma et donc la dynamique du processus. Une analyse préliminaire montre que le rapport entre les constantes de temps d'affaiblissement de la résistance et de la pression est un facteur de contrôle de la dynamique du processus pré-éruptif et éruptif.L'endommagement d'un édifice volcanique ne peut continuer indéfiniment, cela interdirait le stockage de magma dans l'édifice à long terme. L'édifice volcanique connait donc nécessairement des phases de consolidation. Une première approche d'étude de la consolidation reposant sur l'utilisation d'un critère basé sur les variations relatives de volume a été investiguée. Les résultats montrent cependant qu'une telle approche n'est pas suffisante pour expliquer la nécessaire augmentation des modules élastiques, et que les processus de consolidation chimique sans variation de volume doivent avoir une importance fondamentale dans l'existence de cycles éruptifs. Pour étudier ces processus, il faut disposer d'une observable géophysique indépendante du déplacement de surface - nous proposons d'utiliser l'anisotropie de propagation des ondes de cisaillement dans l'édifice volcanique. Ces variations de résistance de l'édifice peuvent contrôler la pression en profondeur. / On volcanoes accelerations of the cumulated number of earthquakes and sometimes of the surface displacement occur in the hours, days and weeks preceding the eruptions. Explaining an acceleration of surface displacement using an elastic model for the volcanic edifice implies an acceleration of the pressure increase in the magma reservoir and in the mantle, at the base of the magma conduit, at short timescales, that is not realistic. Another direction for research consists in considering that an eruption is due to the rupture of the magma reservoir, and that this rupture is not instantaneous - it can be preceded by progressive damage of the edifice, during which low-magnitude earthquakes occur. These ruptures provoke the progressive weakening of the edifice elastic moduli. In this work we have inverstigated the effect of this weakening on the surface displacements, on the reservoir pressure and on the magma flow. We have first searched for a damage law allowing fitting the measured surface displacements by using the cumulated number of recorded earthquakes. The law allowing this fit is a power-law of the incremental damage parameter, with the cumulated number of earthquakes at the exponent. This incremental damage parameter is the ratio between the characteristic incremental rupture length and the length to be ruptured for the eruption.This model allows us to explain weak and strong rates and accelerations of the surface displacement, with the value of the damage parameter. Bifurcation that occurs between a stable state of the edifice and an instable state leading to an eruption is also explained. The weakening of the edifice induces simultaneously a diminution of the reservoir overpressure and an acceleration of the surface displacement, overpressure decreasing more slowly than the elastic moduli, due to the magma feeding of the reservoir.This type of model can explain a wide range of pre-eruptive and eruptive behaviours on volcanoes. It allows the fitting of surface displacement accelerations recorded at the GPS summit stations of Piton de la Fournaise volcano for the March 30th, 2007 eruption, and explains the steady increase of the surface displacement recorded at the summit of Grimsvötn volcano (Iceland), between 2005 and 2011. It may be used for understanding the complex eruptive behaviours, comprising multiple episodes, as in the case of the October-November 2010 eruption of the Merapi volcano (Indonesia). The value of the characteristic time of the feeding system controls the time delay between the maximum of damage (end of the volcano-tectonic seismic crisis) and the maximum of the surface displacement or magma flow. A preliminary analysis shows that the ratio between the weakening characteristic times of the edifice strength and reservoir pressure participate to the control of the pre-eruptive and eruptive process dynamics.Damage of the volcanic edifice can not continue indefinitely, otherwise long-term magma storage would be impossible in the volcanic edifice. The volcanic edifice therefore necessarily experiments periods of strenghtening. A first approach for studying this strengthening relying on a criterion based on the volume relative variation of the edifice has been investigated. Results show that this approach is not sufficient to explain the necessary increase of the elastic moduli. Chemical strenghtening processes without volume changes should have a fundamental role in the existence of eruptive cycles. For improving our knowledge of these processes, we need a further geophysical observable, independent on the surface displacements - we propose to use shear-wave anisotropy in the volcanic edifice.We finally show that the surface dynamics imposed by strength changes of the edifice can influence the dynamics of the deep processes, and consider the possibility of a retroaction of the edifice strength changes on the deep source processes.

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. / 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.