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

Melt Inclusion Vapour Bubbles: The Hidden Reservoir for Major and Volatile Elements

Olivine-hosted melt inclusions (MIs) provide samples of magmatic liquids and their dissolved volatiles from deep within the plumbing system. Inevitable post-entrapment modifications can lead to significant compositional changes in the glass and/or any contained bubbles. Re-heating is a common technique to reverse MI crystallisation; however, its effect on volatile contents has been assumed to be minor. We test this assumption using crystallised and glassy basaltic MIs, combined with Raman spectroscopy and 3D imaging, to investigate the changes in fluid and solid phases in the bubbles before and after re-heating. Before re-heating, the bubble contains CO2 gas and anhydrite (CaSO4) crystallites. The rapid diffusion of major and volatile elements from the melt during re-heating creates new phases within the bubble: SO2, gypsum, Fe-sulphides. Vapour bubbles hosted in naturally glassy MIs similarly contain a plethora of solid phases (carbonates, sulphates, and sulphides) that account for up to 84% of the total MI sulphur, 80% of CO2, and 14% of FeO. In both re-heated and naturally glassy MIs, bubbles sequester major and volatile elements that are components of the total magmatic budget and represent a “loss” from the glass. Analyses of the glass alone significantly underestimates the original magma composition and storage parameters

Timescales of Quartz Crystallization and the Longevity of the Bishop Giant Magma Body

Supereruptions violently transfer huge amounts (100 s–1000 s km3) of magma to the surface in a matter of days and testify to the existence of giant pools of magma at depth. The longevity of these giant magma bodies is of significant scientific and societal interest. Radiometric data on whole rocks, glasses, feldspar and zircon crystals have been used to suggest that the Bishop Tuff giant magma body, which erupted ∼760,000 years ago and created the Long Valley caldera (California), was long-lived (>100,000 years) and evolved rather slowly. In this work, we present four lines of evidence to constrain the timescales of crystallization of the Bishop magma body: (1) quartz residence times based on diffusional relaxation of Ti profiles, (2) quartz residence times based on the kinetics of faceting of melt inclusions, (3) quartz and feldspar crystallization times derived using quartz+feldspar crystal size distributions, and (4) timescales of cooling and crystallization based on thermodynamic and heat flow modeling. All of our estimates suggest quartz crystallization on timescales of <10,000 years, more typically within 500–3,000 years before eruption. We conclude that large-volume, crystal-poor magma bodies are ephemeral features that, once established, evolve on millennial timescales. We also suggest that zircon crystals, rather than recording the timescales of crystallization of a large pool of crystal-poor magma, record the extended periods of time necessary for maturation of the crust and establishment of these giant magma bodies

Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia

The 2010 eruption of Merapi (VEI 4) was the volcano’s largest since 1872. In contrast to the prolonged and effusive dome-forming eruptions typical of Merapi’s recent activity, the 2010 eruption began explosively, before a new dome was rapidly emplaced. This new dome was subsequently destroyed by explosions, generating pyroclastic density currents (PDCs), predominantly consisting of dark coloured, dense blocks of basaltic andesite dome lava. A shift towards open-vent conditions in the later stages of the eruption culminated in multiple explosions and the generation of PDCs with conspicuous grey scoria and white pumice clasts resulting from sub-plinian convective column collapse. This paper presents geochemical data for melt inclusions and their clinopyroxene hosts extracted from dense dome lava, grey scoria and white pumice generated during the peak of the 2010 eruption. These are compared with clinopyroxene-hosted melt inclusions from scoriaceous dome fragments from the prolonged dome-forming 2006 eruption, to elucidate any relationship between pre-eruptive degassing and crystallisation processes and eruptive style. Secondary ion mass spectrometry analysis of volatiles (H2O, CO2) and light lithophile elements (Li, B, Be) is augmented by electron microprobe analysis of major elements and volatiles (Cl, S, F) in melt inclusions and groundmass glass. Geobarometric analysis shows that the clinopyroxene phenocrysts crystallised at depths of up to 20 km, with the greatest calculated depths associated with phenocrysts from the white pumice. Based on their volatile contents, melt inclusions have re-equilibrated during shallower storage and/or ascent, at depths of ~0.6–9.7 km, where the Merapi magma system is interpreted to be highly interconnected and not formed of discrete magma reservoirs. Melt inclusions enriched in Li show uniform “buffered” Cl concentrations, indicating the presence of an exsolved brine phase. Boron-enriched inclusions also support the presence of a brine phase, which helped to stabilise B in the melt. Calculations based on S concentrations in melt inclusions and groundmass glass require a degassing melt volume of 0.36 km3 in order to produce the mass of SO2 emitted during the 2010 eruption. This volume is approximately an order of magnitude higher than the erupted magma (DRE) volume. The transition between the contrasting eruptive styles in 2010 and 2006 is linked to changes in magmatic flux and changes in degassing style, with the explosive activity in 2010 driven by an influx of deep magma, which overwhelmed the shallower magma system and ascended rapidly, accompanied by closed-system degassing

Thoughts on the Criteria to Determine the Origin of Volcanic Unrest as Magmatic or Non-Magmatic

As our ability to detect volcanic unrest improves, we are increasingly confronted with the question of whether the unrest has a magmatic origin (magma on the move) or a non-magmatic origin from a change in the hydrothermal system (fluids that are not magma on the move) or tectonic processes. The cause of unrest has critical implications for the potential eruptive hazard (e.g. used in constructing Bayesian Event Trees), but is frequently the subject of debate, even at well-studied systems. Here, we propose a set of multi-disciplinary observations and numerical models that could be used to evaluate conceptual models about the cause of unrest. These include measurements of gas fluxes and compositions and the isotopic signature of some components (e.g. H2, He, C, SO2, H2O, CH4 and CO2), the spatial and temporal characteristics of ground deformation, thermal output, seismicity, changes in gravity, and whether there is topographic uplift or subsidence spanning hundreds to thousands of years. In several volcanic systems, both magmatic and non-magmatic unrest is occurring at the same time. While none of these observations or models is diagnostic on its own, we illustrate several examples where they have been used together to make a plausible conceptual model of one or more episodes of unrest and whether eruptions did or did not follow the unrest