A Method for Magma Viscosity Assessment by Lava Dome Morphology

Lava domes form when a highly viscous magma erupts on the surface. Several types of lava dome morphology can be distinguished depending on the flow rate and the rheology of magma: obelisks, lava lobes, and endogenic structures. The viscosity of magma nonlinearly depends on the volume fraction of crystals and temperature. Here we present an approach to magma viscosity estimation based on a comparison of observed and simulated morphological forms of lava domes. We consider a two-dimensional axisymmetric model of magma extrusion on the surface and lava dome evolution, and assume that the lava viscosity depends only on the volume fraction of crystals. The crystallization is associated with a growth of the liquidus temperature due to the volatile loss from the magma, and it is determined by the characteristic time of crystal content growth (CCGT) and the discharge rate. Lava domes are modeled using a finite-volume method implemented in Ansys Fluent software for various CCGTs and volcanic vent sizes. For a selected eruption duration a set of morphological shapes of domes (shapes of the interface between lava dome and air) is obtained. Lava dome shapes modeled this way are compared with the observed shape of the lava dome (synthesized in the study by a random modification of one of the calculated shapes). To estimate magma viscosity, the deviation between the observed dome shape and the simulated dome shapes is assessed by three functionals: the symmetric difference, the peak signal-to-noise ratio, and the structural similarity index measure. These functionals are often used in the computer vision and in image processing. Although each functional allows to determine the best fit between the modeled and observed shapes of lava dome, the functional based on the structural similarity index measure performs it better. The viscosity of the observed dome can be then approximated by the viscosity of the modeled dome, which shape fits best the shape of the observed dome. This approach can be extended to three-dimensional case studies to restore the conditions of natural lava dome growth.

A numerical modelling toolbox for identifying the expression of dome-forming volcanism on exoplanets

<p>The presence of volcanism is often anecdotally used to define a “living planet”. Since dome-building volcanism on Earth occurs primarily at plate boundaries, the identification of such domes could inform on exoplanetary development. Lava domes form when extruded magma is too viscous to flow from a vent, and their morphology on Earth varies from flat, pancake lobes to steep, blocky domes. Identification of lava domes on other terrestrial planets in our Solar System indicates that they likely also exist on rocky exoplanets. Here we show, using particle-based modelling, that the diversity of lava dome morphology in our Solar System is dwarfed by the diversity expected for exoplanets. Specifically, the height-to-diameter ratio of a dome decreases as a function of increasing gravity (i.e., planetary mass and radius). For example, lava domes on high-gravity super-Earths will be extremely wide and flat and a volcanic origin may not be immediately apparent. Creating a toolbox to help identify exoplanetary volcanism will allow us to make initial estimations as to the development and habitability of these alien worlds as images become available.</p>

Characterising fracture patterns at growing lava domes

<p>Volcanic domes form when lava is too viscous to flow away from an active volcanic vent; instead, the lava accumulates into a mound consisting of a hotter, ductile core and a colder, brittle outer layer. An existing lava dome grows when new material is injected into the core of the dome, causing the  outer layer to stretch and develop tensile fractures. With continued dome growth, these weaknesses can propagate to form an extensive fracture network and the dome may fail. Collapse events often generate rock falls and debris avalanches, lahars, and high-speed pyroclastic flows, endangering populations residing at the base of a volcano. Since such fractures represent potential failure planes, in this project we aim to understand the role they have in destabilising lava domes.</p><p>This project will build on the work published by Harnett et al. (2018), which demonstrates the suitability of a discrete element modelling approach to simulate dome emplacement and evolution. Specifically, this project is designed to:</p><p>1. Use high-resolution photogrammetry to characterise the possible fracture states of a dome;</p><p>2. Establish up-scaled rock-mass properties by performing geomechanical experiments on both fractured and non-fractured samples of dome rock from prior collapses;</p><p>3. Develop a numerical model to investigate how the presence and properties of fracture networks influence dome stability.</p><p>The model, developed using PFC, will be used to identify critical fracture states that can signify a dome collapse is likely to occur. Under the current model, parallel bonds simulate the fluid magma core and flat joints simulate the solid talus material. This project will build on this original model by incorporating discrete fracture networks into the smooth-joint model to implement dome fracturing. The new model will look to investigate the effect  of a fracture network on a static dome that, when in its unfractured state, is stable under gravity. Additionally, the model will be designed such that inputs can include experimentally derived rock-mass properties. It is hoped that, by incorporating observational and experimental data into a more  complex model, the dynamic evolution of fractures in a growing lava dome can be investigated and the ongoing likelihood of a dome collapse event can be assessed.</p><p> </p><p>Harnett, C. E. et al., 2018. J. Volcanol. Geoth. Res., 359: 68-77.</p>

Effusive Monogenetic Volcanism

The study of monogenetic volcanism around Earth is rapidly growing due to the increasing recognition of monogenetic volcanic edifices in different tectonic settings. Far from the idea that this type of volcanism is both typically mafic and characteristic from intraplate environments, it occurs in a wide spectrum of composition and geological settings. This volcanism is widely known by the distinctive pyroclastic cones that represent both magmatic and phreatomagmatic explosive activity; they are known as scoria or spatter cones, tuff cones, tuff rings, maars and maar-diatremes. These cones are commonly associated with lava domes and usually accompanied by lava flows as part of their effusive eruptive phases. In spite of this, isolated effusive monogenetic emissions also appear around Earth’s surface. However, these isolated emissions are not habitually considered within the classification scheme of monogenetic volcanoes. Along with this, many of these effusive volcanoes also contrast with the belief that this volcanism is indicative of rapidly magma ascent from the asthenosphere, as many of the products are strongly evolved reflecting differentiation linked to stagnation during ascent. This has led to the understanding that the asthenosphere is not always the place that directly gives rise to the magma batches and rather, they detach from a crustal melt storage. This chapter introduces four singular effusive monogenetic volcanoes as part of the volcanic geoforms, highlights the fact that monogenetic volcanic fields can also be associated with crustal reservoirs, and outlines the processes that should occur to differentiate the magma before it is released as intermediate and acidic in composition. This chapter also provides an overview of this particular volcanism worldwide and contributes to the monogenetic comprehension for future studies.

4. Lava

‘Lava’ focuses on lava. Within a kilometre or two of the ground surface, the rise of hot magma is accelerated by the growth of gas bubbles which increase the magma’s buoyancy. If the gases can readily escape from the magma, the magma may ascend more slowly without being ripped apart, and it eventually emerges from the volcano as a lava flow. There are many different types of lava flows: pahoehoe, a’a, block-lavas, inflated lavas, and lava domes, which decrepitate to form hazardous ‘block-and-ash flows’.

The rebirth and evolution of Bezymianny volcano, Kamchatka after the 1956 sector collapse

Continued post-collapse volcanic activity can cause the rise of a new edifice. However, details of such edifice rebirth have not been documented yet. Here, we present 7-decade-long photogrammetric data for Bezymianny volcano, Kamchatka, showing its evolution after the 1956 sector collapse. Edifice rebirth started with two lava domes originating at distinct vents ~400 m apart. After 2 decades, activity became more effusive with vents migrating within ~200 m distance. After 5 decades, the activity focused on a single vent to develop a stratocone with a summit crater. We determine a long-term average growth rate of 26,400 m3/day, allowing us to estimate the regain of the pre-collapse size within the next 15 years. Numerical modeling explains the gradual vents focusing to be associated with loading changes, affecting magma pathways at depth. This work thus sheds light on the complex regrowth process following a sector collapse, with implications for regrowing volcanoes elsewhere.

Lava dome morphology inferred from numerical modelling

SUMMARY Lava domes form when highly viscous magmas erupt on the surface. Several types of lava dome morphology can be distinguished depending on the flow rate and the rheology of magma. Here, we develop a 2-D axisymmetric model of magma extrusion on the surface and lava dome evolution and analyse the dome morphology using a finite-volume method implemented in Ansys Fluent software. The magma/lava viscosity depends on the volume fraction of crystals and temperature. We show that the morphology of domes is influenced by two parameters: the characteristic time of crystal content growth (CCGT) and the discharge rate (DR). At smaller values of the CCGTs, that is, at rapid lava crystallization, obelisk-shaped structures develop at low DRs and pancake-shaped structures at high DRs; at longer CCGTs, lava domes feature lobe- to pancake-shaped structures. A thick carapace of about 70 per cent crystal content evolves at smaller CCGTs. We demonstrate that cooling does not play the essential role during a lava dome emplacement, because the thermal thickness of the evolving carapace remains small in comparison with the dome's height. A transition from the endogenic to exogenic regime of the lava dome growth occurs after a rapid increase in the DR. A strain-rate-dependent lava viscosity leads to a more confined dome, but the influence of this viscosity on the dome morphology is not well pronounced. The model results can be used in assessments of the rates of magma extrusion, the lava viscosity and the morphology of active lava domes..