How caldera collapse shapes the shallow emplacement and transfer of magma in active volcanoes

Calderas are topographic depressions formed by the collapse of a partly drained magma reservoir. At volcanic edifices with calderas, eruptive fissures can circumscribe the outer caldera rim, be oriented radially and/or align with the regional tectonic stress field. Constraining the mechanisms that govern this spatial arrangement is fundamental to understand the dynamics of shallow magma storage and transport and evaluate volcanic hazard. Here we show with numerical models that the previously unappreciated unloading effect of caldera formation may contribute significantly to the stress budget of a volcano. We first test this hypothesis against the ideal case of Fernandina, Galápagos, where previous models only partly explained the peculiar pattern of circumferential and radial eruptive fissures and the geometry of the intrusions determined by inverting the deformation data. We show that by taking into account the decompression due to the caldera formation, the modeled edifice stress field is consistent with all the observations. We then develop a general model for the stress state at volcanic edifices with calderas based on the competition of caldera decompression, magma buoyancy forces and tectonic stresses. These factors control: 1) the shallow accumulation of magma in stacked sills, consistently with observations; 2) the conditions for the development of circumferential and/or radial eruptive fissures, as observed on active volcanoes. This top-down control exerted by changes in the distribution of mass at the surface allows better understanding of how shallow magma is transferred at active calderas, contributing to forecasting the location and type of opening fissures

Crustal-scale degassing due to magma system destabilization and magma-gas decoupling at Soufrière Hills Volcano, Montserrat

Activity since 1995 at Soufrière Hills Volcano (SHV), Montserrat has alternated between andesite lava extrusion and quiescence, which are well correlated with seismicity and ground deformation cycles. Large variations in SO₂ flux do not correlate with these alternations, but high and low HCl/SO₂ characterize lava dome extrusion and quiescent periods respectively. Since lava extrusion ceased (February 2010) steady SO₂ emissions have continued at an average rate of 374 tonnes/day (± 140 t/d), and incandescent fumaroles (temperatures up to 610°C) on the dome have not changed position or cooled. Occasional short bursts (over several hours) of higher (∼ 10x) SO₂ flux have been accompanied by swarms of volcano-tectonic earthquakes. Strain data from these bursts indicate activation of the magma system to depths up to 10 km. SO₂ emissions since 1995 greatly exceed the amounts that could be derived from 1.1 km³ of erupted andesite, and indicating extensive partitioning of sulfur into a vapour phase, as well as efficient decoupling and outgassing of sulfur-rich gases from the magma. These observations are consistent with a vertically extensive, crustal magmatic mush beneath SHV. Three states of the magmatic system are postulated to control degassing. During dormant periods (10³ to 10⁴ years) magmatic vapour and melts separate as layers from the mush and decouple from each other. In periods of unrest (years) without eruption, melt and fluid layers become unstable, ascend and can amalgamate. Major destabilization of the mush system leads to eruption, characterized by magma mixing and release of volatiles with different ages, compositions and sources.RSJS acknowledges an ERC advanced grant (VOLDIES). JDB acknowledges ERC advanced grant CRITMAG and a Wolfson Research Merit Award.This is the final version of the article. It first appeared from Wiley via http://dx.doi.org/10.1002/2015GC00579

Modeling the impact of melt on seismic properties during mountain building

Initiation of partial melting in the mid/lower crust causes a decrease in P-wave and S-wave velocities; recent studies imply that the relationship between these velocities and melt is not simple. We have developed a modelling approach to assess the combined impact of various melt and solid phase properties on seismic velocities and anisotropy. The modelling is based on crystallographic preferred orientation (CPO) data measured from migmatite samples, allowing quantification of the variation of seismic velocities with varying melt volumes, shapes, orientations, and matrix anisotropy. The results show non-linear behaviour of seismic properties as a result of the interaction of all of these physical properties, which in turn depend on lithology, stress regime, strain rate, pre-existing rock fabrics, and pressure-temperature conditions. This non-linear behaviour is evident when applied to a suite of samples from a traverse across a migmatitic shear zone in the Seiland Igneous Province, Northern Norway. Critically, changes in solid phase composition and CPO, and melt shape and orientation with respect to the wave propagation direction can result in huge variations in the same seismic property even if the melt fraction remains the same. A comparison with surface wave interpretations from tectonically active regions highlights the issues in current models used to predict melt percentages or partially molten regions. Interpretation of seismic data to infer melt percentages or extent of melting should, therefore, always be underpinned by robust modelling of the underlying geological parameters combined with examination of multiple seismic properties in order to reduce uncertainty of the interpretation

Post-supereruption recovery at Toba Caldera

Large calderas, or supervolcanoes, are sites of the most catastrophic and hazardous events on Earth, yet the temporal details of post-supereruption activity, or resurgence, remain largely unknown, limiting our ability to understand how supervolcanoes work and address their hazards. Toba Caldera, Indonesia, caused the greatest volcanic catastrophe of the last 100 kyr, climactically erupting ~74 ka. Since the supereruption, Toba has been in a state of resurgence but its magmatic and uplift history has remained unclear. Here we reveal that new 14 C, zircon U-Th crystallization and (U-Th)/He ages show resurgence commenced at 69.7±4.5 ka and continued until at least ~2.7 ka, progressing westward across the caldera, as reflected by post-caldera effusive lava eruptions and uplifted lake sediment. The major stratovolcano north of Toba, Sinabung, shows strong geochemical kinship with Toba, and zircons from recent eruption products suggest Toba's climactic magma reservoir extends beneath Sinabung and is being tapped during eruptions

Magma plumbing systems: a geophysical perspective

Over the last few decades, significant advances in using geophysical techniques to image the structure of magma plumbing systems have enabled the identification of zones of melt accumulation, crystal mush development, and magma migration. Combining advanced geophysical observations with petrological and geochemical data has arguably revolutionised our understanding of, and afforded exciting new insights into, the development of entire magma plumbing systems. However, divisions between the scales and physical settings over which these geophysical, petrological, and geochemical methods are applied still remain. To characterise some of these differences and promote the benefits of further integration between these methodologies, we provide a review of geophysical techniques and discuss how they can be utilised to provide a structural context for and place physical limits on the chemical evolution of magma plumbing systems. For example, we examine how Interferometric Synthetic Aperture Radar (InSAR), coupled with Global Positioning System (GPS) and Global Navigation Satellite System (GNSS) data, and seismicity may be used to track magma migration in near real-time. We also discuss how seismic imaging, gravimetry and electromagnetic data can identify contemporary melt zones, magma reservoirs and/or crystal mushes. These techniques complement seismic reflection data and rock magnetic analyses that delimit the structure and emplacement of ancient magma plumbing systems. For each of these techniques, with the addition of full-waveform inversion (FWI), the use of Unmanned Aerial Vehicles (UAVs) and the integration of geophysics with numerical modelling, we discuss potential future directions. We show that approaching problems concerning magma plumbing systems from an integrated petrological, geochemical, and geophysical perspective will undoubtedly yield important scientific advances, providing exciting future opportunities for the volcanological community

A large magmatic sill complex beneath the Toba caldera

International audienceAn understanding of the formation of large magmatic reservoirs is a key issue for the evaluation of possible strong volcanic eruptions in the future. We estimated the size and level of maturity of one of the largest volcanic reservoirs, based on radial seismic anisotropy. We used ambient-noise seismic tomography below the Toba caldera (in northern Sumatra) to observe the anisotropy that we interpret as the expression of a fine-scale layering caused by the presence of many partially molten sills in the crust below 7 kilometers. This result demonstrates that the magmatic reservoirs of present (non-eroded) supervolcanoes can be formed as large sill complexes and supports the concept of the long-term incremental evolution of magma bodies that lead to the largest volcanic eruptions

Rapid changes in magma storage beneath the Klyuchevskoy group of volcanoes inferred from time-dependent seismic tomography

We present the results of time-dependent local earthquake tomography for the Kluchevskoy group of volcanoes in Kamchatka, Russia. We consider the time period from 1999 to 2009, which covers several stages of activity of Kluchevskoy and Bezymianny volcanoes. The results are supported by synthetic tests that recover a common 3D model based on data corresponding to different time windows. Throughout the period, we observe a robust feature below 25 km depth with anomalously high Vp/Vs values (up to 2.2). We interpret this feature as a channel bringing deep mantle materials with high fluid and melt content to the bottom of the crust. This mantle channel directly or indirectly determines the activity of all volcanoes of the Kluchevskoy group. In the crust, we model complex structure that varies over time. During the pre-eruptive period, we detected two levels of potential magma storage: one in the middle crust at 10–12 km depth and one close to the surface just below Kluchevskoy volcano. In 2005, a year of powerful eruptions of Kluchevskoy and Besymiyanny volcanoes, we observe a general increase in Vp/Vs throughout the crust. In the relaxation period following the eruption, the Vp/Vs values are generally low, and no strong anomalous zones in the crust are observed. We propose that very rapid variations in Vp/Vs are most likely due to abrupt changes in the stress and deformation states, which cause fracturing and the active transport of fluids. These fluids drive more fracturing in a positive feedback system that ultimately leads to eruption. We envision the magma reservoirs beneath the Kluchevskoy group as sponge-structured volumes that may quickly change the content of the molten phases as fluids pulse rapidly through the system

Asymmetric caldera-related structures in the area of the Avacha group of volcanoes in Kamchatka as revealed by ambient noise tomography and deep seismic sounding

Avacha group includes two active and potentially dangerous volcanoes, Avachinsky and Koryaksky, located close to Petropavlovsk-Kamchatsky, the main city of Kamchatka. We present the results of two independent seismic studies of shallow crustal structures beneath the Avacha group based on passive and active source observations. The first study is based on the analysis of continuous recording by 11 seismic stations installed over the Avacha group in 2012 and 7 permanent stations in the same region. We present a series of 2D Rayleigh-wave group velocity maps based on correlation of ambient noise, that were then converted into 3D distribution of shear wave velocity. The second work was based on the reprocessing of an active source deep seismic sounding profile across the Avachinsky volcano that was shot in 1982–1984. We made the analysis of travel times of refracted waves using a 2D tomography inversion. The resulting seismic models appear to be consistent with each other and show clear low-velocity zone to the SW of the Avachinsky volcano and high velocity structures to NE. These observations also agree with the existing gravity and magnetotelluric measurements. Based on the obtained seismic models we identify two large buried calderas and large lava flows that are thought to be related to a series of large eruption episodes of Avachinsky occurred within the last 30,000 years

Asymmetric caldera-related structures in the area of the Avacha group of volcanoes in Kamchatka as revealed by ambient noise tomography and deep seismic sounding

International audienceAvacha group includes two active and potentially dangerous volcanoes, Avachinsky and Koryaksky, located close to Petropavlovsk-Kamchatsky, the main city of Kamchatka. We present the results of two independent seismic studies of shallow crustal structures beneath the Avacha group based on passive and active source observations. The first study is based on the analysis of continuous recording by 11 seismic stations installed over the Avacha group in 2012 and 7 permanent stations in the same region. We present a series of 2D Rayleigh-wave group velocity maps based on correlation of ambient noise, that were then converted into 3D distribution of shear wave velocity. The second work was based on the reprocessing of an active source deep seismic sounding profile across the Avachinsky volcano that was shot in 1982-1984. We made the analysis of travel times of refracted waves using a 2D tomography inversion. The resulting seismic models appear to be consistent with each other and show clear low-velocity zone to the SW of the Avachinsky volcano and high velocity structures to NE. These observations also agree with the existing gravity and magnetotelluric measurements. Based on the obtained seismic models we identify two large buried calderas and large lava flows that are thought to be related to a series of large eruption episodes of Avachinsky occurred within the last 30,000 years