Lava Dome Morphology and Viscosity Inferred From Data-Driven Numerical Modeling of Dome Growth at Volcán de Colima, Mexico During 2007-2009

Magma extrusion, lava dome growth, collapse of domes, and associated pyroclastic flow hazards are among important volcanological studies. In this paper, we analyze the influence of the magma viscosity and discharge rates on the lava dome morphology at Volcán de Colima in Mexico during a long dome-building episode lasting from early 2007 to fall 2009 without explosive dome destruction. Camera images of the lava dome growth together with recorded volumes of the erupted lava have been used to constrain numerical modeling and hence to match the history of the dome growth by nudging model forecasts to observations. Our viscosity model incorporates crystal growth kinetics and depends on the characteristic time of crystal content growth (or CCGT) and the crystal-free magma viscosity. Initially, we analyze how this viscosity, CCGT, and the rate of lava extrusion influence the morphology of the growing dome. Several model scenarios of lava dome growth are then considered depending on the crater geometry, the conduit location, the effective viscosity of dome carapace, and the extrusion rates. These rates are determined either empirically by optimizing the fit between the morphological shape of modeled domes and that of the observed dome or from the recorded lava dome volumes. The maximum height of the modeled lava dome and its horizontal extent are in a good agreement with observations in the case of the empirically-derived extrusion rates. It is shown that the topography of the crater at Volcán de Colima is likely to be inclined toward the west. The viscosity of the modeled lava dome (∼1012 Pa s) is in a good agreement with the effective viscosity estimated experimentally from lavas of Volcán de Colima. Due to the interplay between the lava extrusion and the gravity forces, the dome reaches a height threshold, and after that a horizontal gravity spreading starts to play an essential role in the lava dome evolution. The model forecasts that the dome carapace of higher viscosity (∼1014 Pa s) influences the dome growth and its morphology during long dome-building episodes by retarding horizontal advancement and developing steep-sided eastern edge of the dome at the volcano. The developed model can be used in assessments of future effusive eruptions and lava dome growth at Volcán de Colima or elsewhere. History matching modeling of lava dome growth sheds a light on dynamic processes inside the dome and may assist in assessing stress state in the dome carapace and in forecasting the dome failures.

Syn-Eruptive Conditions of the AD 1530 Sub-Plinian Eruption of La Soufrière of Guadeloupe (Lesser Antilles)

A prerequisite in refining volcanic hazard at explosive volcanoes is a better quantification of the timescales of the syn-eruptive processes, such as magma degassing and crystallization prior to eruption. To this aim, new data on the matrices (microlites, residual glasses, and bubbles) of pumice, scoria, and dense clasts erupted during the AD 1530 andesitic eruption of La Soufrière of Guadeloupe are combined with published data from phase-equilibrium and kinetic experiments, in order to estimate pressures of microlite crystallization and magma ascent rates. From the timescale data, we infer that the AD 1530 eruption started with phreatomagmatic explosions tapping magmas that ascended during about 1 month (decompression rate of ∼50 Pa/s) from the coldest parts of the reservoir (∼825°C and a 74 wt% SiO2 melt). These magmas continuously crystallized microlites (∼25 vol% plagioclase, orthopyroxene, clinopyroxene, magnetite, quartz, and apatite), but did not outgas likely due to sealed conduit margins. The flank collapse (preexisting “cold” edifice) that followed the phreatomagmatic phase triggered a sub-Plinian eruption that progressively tapped the hotter main reservoir (∼875°C and 71 wt% SiO2 interstitial melt), emitting banded and homogeneous pumice. The banded pumice did not significantly outgas and mostly lack decompression-induced microlites, suggesting short ascent durations of the order of 0.5–1 day (decompression rates of 1,400–4,000 Pa/s). The following Strombolian phase emitted dark scoria that did not significantly outgas and only crystallized rare microlites, suggesting ascent duration of the order of 2 days (decompression rates of ∼550 Pa/s). The terminal lava dome growth involved fully outgassed magmas ascended during more than 1 month, giving time for microlite crystallization (∼40 vol% plagioclase, orthopyroxene, clinopyroxene, magnetite, and cristobalite). The detection of any shallow new magmatic intrusion is crucial, since it can trigger a sequence of conduit processes leading to an eruption marked by a succession of different and disastrous eruptive styles, following a scenario similar to the AD 1530 eruption. Overall, we provide a combined approach of petrological, geochemical, and experimental data that may be used to infer ascent conditions and rates at other volcanic systems.

Assessing the Use of Optical Satellite Images to Detect Volcanic Impacts on Glacier Surface Morphology

Globally, about 250 Holocene volcanoes are either glacier-clad or have glaciers in close proximity. Interactions between volcanoes and glaciers are therefore common, and some of the most deadly (e.g., Nevado del Ruiz, 1985) and most costly (e.g., Eyjafjallajökull, 2010) eruptions of recent years were associated with glaciovolcanism. An improved understanding of volcano-glacier interactions is therefore of both global scientific and societal importance. This study investigates the potential of using optical satellite images to detect volcanic impacts on glaciers, with a view to utilise detected changes in glacier surface morphology to improve glacier-clad volcano monitoring and eruption forecasting. Roughly 1400 optical satellite images are investigated from key, well-documented eruptions around the globe during the satellite remote sensing era (i.e., 1972 to present). The most common observable volcanic impact on glacier morphology (for both thick and thin ice-masses) is the formation of ice cauldrons and openings, often associated with concentric crevassing. Other observable volcanic impacts include ice bulging and fracturing due to subglacial dome growth; localized crevassing adjacent to supraglacial lava flows; widespread glacier crevassing, presumably, due to meltwater-triggered glacier acceleration and advance. The main limitation of using optical satellite images to investigate changes in glacier morphology is the availability of cloud- and eruption-plume-free scenes of sufficient spatial- and temporal resolution. Therefore, for optimal monitoring and eruption prediction at glacier-clad volcanoes, optical satellite images are best used in combination with other sources, including SAR satellite data, aerial images, ground-based observations and satellite-derived products (e.g., DEMs).

Landscape and hazard evolution during the Montserrat volcanic crisis 1995–2010: an integrated simulation with r.avaflow

<p>An active phase of Soufrière Hills Volcano (Montserrat, Lesser Antilles) has started in 1995 and had its most intense period between 1995 and 2010, when phases of lava dome growth were interrupted by dome collapses triggering ash clouds and different types of pyroclastic flows. These flows were released in various directions, so that two thirds of the island were left in an inhabitable state. The material deposited was later remobilized through lahar flows, burying the centre of the former capital town of Plymouth. In the present work, we attempt to back-calculate the sequences of dome growth – pyroclastic flows, and the subsequent lahar flows, in an integrated way, using the mass flow simulation tool r.avaflow. Thereby, we build on the reconstruction of the pre-event topography as well as on various reference data obtained from the large amount of available literature – mainly, the peak elevation and volumes of the lava domes, the impact areas of the flow processes, and ash fall characteristics. Most observations are successfully reproduced with physically plausible, though calibrated, parameter sets and temporal scaling of lava dome growth. Due to the complexity and multi-stage nature of the volcanic crisis, a number of simplifications had to be introduced, such as considering only the twelve largest sequences of dome growth and pyroclastic flows, and evaluating some of the results on the basis of aggregated impact areas for more than one event. Consequently, the results reflect a strong conceptual component, but can - at least in part - be considered useful for predictive modelling of similar events. Another scope of the simulation results, however, is its educational use. Appropriately presented, they greatly facilitate the generation of a better understanding for complex chains of volcanic processes and their consequences to learners at various levels in different educational contexts, from school and university all the way to targeted awareness-building campaigns.</p>

Long-term evolution of the gas composition of Popocatepetl's plume

<p>Changes in the eruptive dynamics are mainly controlled by the magma gas content, and the degassing processes impacting the magma viscosity and ascending speed. The progressive exsolution of the gas species, their release at different depths, their mutual interaction and the eventual assimilation of crustal rocks are reflected in the volcanic plume composition changes. Combining long-term ground-based FTIR and UV remote measurements of the Popocatepetl's plume, seismic data and visual monitoring, we explore the relationship between the gas composition changes in the volcanic plume and the transition between extrusive and passive degassing regimes.</p><p>SO2, HCl, HF, BrO, SiF4 and CO2 are simultaneously measured in the volcanic plume since 2013 from the Altzomoni observatory, located 12 km north of the crater. We capture several phases of lava dome growth, different types of explosions and passive degassing periods. The evolution of the gas species ratios through these events allows deciphering the degassing processes.</p>

Using glaciers to identify, monitor, and predict volcanic activity

<p>Many (about 250) volcanoes worldwide are occupied by glaciers. Often glaciers are regarded as problematic for volcano monitoring, since glacier ice potentially masks evidence of volcanic activity. The most devastating volcanic eruptions of the last 100 years involved volcano-glacier interactions. The 1985 eruption of Nevado del Ruiz killed 23000 people, and the 2010 eruption of Eyjafjallajökull led to the closure of many European airports. Therefore, it is imperative to minimize these impacts on society by improving methods for monitoring of glacier-clad volcanoes. Amongst several methods, optical satellite remote sensing techniques are perhaps most auspicious, since they frequently have a relatively high temporal and spatial resolution, and are mostly freely available. They often clearly show the effects of volcanic activity on glaciers, including ice cauldron formation, ice fracturing and glacier terminus changes potentially due to subglacial melt or subglacial dome growth. This study has the objective to link pre-, syn- and post-eruption glacier behaviour to the type and timing of volcanic activity, and to develop a satellite based predictive tool for monitoring future eruptions. Despite several studies that link volcanic activity and changing glacier behaviour, the potential of using the latter to predict the former has yet to be systematically tested. Our approach is to observe how glaciers responded to past volcanic events using mostly, but not exclusively optical satellite imagery, and to build a database of examples for potential automated detection and forecasting on a global scale.</p>

Insights into lava dome and spine extrusion using analogue sandbox experiments

<p>Lava dome formation is a common process at stratovolcanoes involving the shallow intrusion or extrusion of viscous lava and may lead to the rise of spines. Spines are protrusions observed to extrude episodically during lava dome growth, yet the structural and mechanical factors controlling their formation are only partially understood. Here, we provide new, detailed insight into lava dome growth and the production of spines using a novel set of analogue experiments extruding sand-plaster mixtures from a fixed-diameter conduit under isothermal conditions. We trace displacement and strain with photogrammetric methods for precise and detailed monitoring of the extrusion process. Results show initial dome growth forming a steep-sided and flat-topped shape through extrusion of new material, leading to slumping of oversteepening slopes, forming a talus. Spines are found to protrude at a later stage through the dome surface along discrete circular faults that originate from the conduit walls, starting a cycle of spine growth and collapse. As our spines only appear after prolonged extrusion, we relate their appearance to the compaction and strengthening of material within the conduit. We find that spine diameter, height and volume are positively correlated with increasing cohesion and therefore material strength. The spine diameter was also observed to be smaller or equal to the diameter of the underlying conduit, as shear extrusion occurs along vertical to outward-dipping fault planes. For natural domes, our findings imply that spine growth may be the consequence of compaction and densification via porosity loss, shearing and/or outgassing of conduit magma during ascent. More efficient compaction will yield wider and taller spines as a result of increasing rock strength. Our study further highlights the relevance of analogue experiments in the study of lava domes and spines, which remain one of the most hazardous and unpredictable features at dome-forming volcanoes worldwide.</p>

Anatomy of the Bezymianny volcano merely before an explosive eruption on 20.12.2017

Strong explosive eruptions of volcanoes throw out mixtures of gases and ash from high-pressure underground reservoirs. Investigating these subsurface reservoirs may help to forecast and characterize an eruption. In this study, we compare seismic tomography results with remote sensing and petrology data to identify deep and subaerial manifestations of pre-eruptive processes at Bezymianny volcano in Kamchatka shortly before its violent explosion on December 20, 2017. Based on camera networks we identify precursory rockfalls, and based on satellite radar data we find pre-eruptive summit inflation. Our seismic network recorded the P and S wave data from over 500 local earthquakes used to invert for a 3D seismic velocity distribution beneath Bezymianny illuminating its eruptive state days before the eruption. The derived tomography model, in conjunction with the presence of the high-temperature-stable SiO2 polymorph Tridymite in juvenile rock samples , allowed us to infer the coexistence of magma and gas reservoirs revealed as anomalies of low (1.5) and high (2.0) Vp/Vs ratios, respectively, located at depths of 2–3 km and only 2 km apart. The reservoirs both control the current eruptive activity: while the magma reservoir is responsible for episodic dome growth and lava flow emplacements, the spatially separated gas reservoir may control short but powerful explosive eruptions of Bezymianny.