The interdisciplinary study of volcanic processes, which extend across all timescales and lengths, requires a multitude of approaches, ranging from analogue and numerical modelling to observations and fieldwork and extending to mathematics. A conference was held at the East African Institute for Fundamental Research, affiliated with the University of Rwanda, a country which, along with the Democratic Republic of Congo, presents a unique geodynamic context. Located along the East African Rift, an active seismic region, Rwanda is close to two of Africa's most active volcanoes, including Nyiragongo, which overlooks Lake Kivu, a deep volcanic lake rich in dissolved carbon dioxide and methane, the latter of which is used for electricity generation. In this context, the conference addressed many “classic” volcanological topics and their modern advances, such as multiphase lava flows, subsurface magma propagation, seismic and deformation signals from a volcano, modelling of volcanic emission dispersion, and volcanic lakes. Yet, it broadened the discussion to the volcanic particle-water interface and its impact on soils, volcanoes and climate change, and volcanoes and health. This article aims to highlight and share the richness of the integration and interconnectedness of the various questions related to a volcanic environment, as well as their impact on society. Ultimately, this conference also demonstrated the importance of promoting science in Africa and developing countries so that the next generation of African researchers is equipped to address the challenges facing their nations.

Surface deformation measured from satellites has provided useful information about the magma plumbing system at active volcanoes. Observed deformation results from complex interactions and coupling between the magma and the host rock. Fracturing of the crust during its deformation can make the pattern of surface displacement even more complex. Models taking into account both the fluid and solid phases of natural systems and linking them are a crucial next step for a better understanding of natural systems and observed deformation. We use the software MFiX (Multiphase Flow with Interphase eXchanges) which considers two phases: a fluid phase computed with a Computational Fluid Dynamics (CFD) method, and a solid phase discretized as spherical particles computed using Discrete Element Methods (DEM) method. Spherical particles are bonded together. Bonds can break at any time step, such that actual fractures can develop through the simulations. We present here the modified drag force between fluid and particles that allows us to model a bonded packing of particles impermeable to a fluid phase. Reproducing a set of analogue experiments, we simulate the injection of fluid in a spherical cavity. Rock tests implemented in MFiX allow us the precise calibration of the packing to the gelatine mechanical properties. The injected volume, the cavity dilatation, the fluid pressure evolution and the surface deformation are measured in the numerical modelling and compared to analogue experiment for benchmarking. We show that this new model has the potential to model the magmatic phase and coupling it to the elastic and brittle deformation of the surrounding rock.

Magma transport in dikes is usually modelled by means of lubrication theory, assuming that magma properties are uniform across the dike. We explore the influence of cross-dike temperature heterogeneity on the dynamics of dike propagation using a quasi-2-D model, derived from a full 2-D model with an assumption of small width to length ratio. The model couples elastic fracture mechanics with multiphase magma flow, solving the governing equations using a hybrid numerical approach that combines the Displacement Discontinuity Method for elasticity with finite volume discretization for fluid flow and heat transfer. The model includes heat exchange with wall rocks, shear heating and latent heat release. It accounts for non-equilibrium magma crystallization, implementing temperature-dependent crystallization kinetics using an Arrhenius formulation for the relaxation timescale. As a case study, we simulate the ascent of a volatile-rich dacite from a source at 30 km depth. The distribution of temperature, crystallinity and, thus, viscosity across the dike leads to a plug-like velocity profile with magma stagnation near the walls, substantially different from the parabolic Poiseuille flow assumed in classical lubrication theory. With temperature-dependent crystallization rate, rapid cooling of magma near the dike walls can generate a glassy chilled margin. The adjacent magma has higher crystallinity due to intermediate cooling rates, while the hotter core remains depleted in crystals throughout dike propagation. The dike propagates further and is thinner than predicted by (1-D) lubrication theory because the low-viscosity core continues to facilitate vertical transport while the wall zones become progressively more viscous due to cooling and crystallization. The latent heat of crystallization can have a substantial impact in slowing down cooling and prolonging propagation. Other important factors include the characteristic crystal growth time, initial magma temperature and water content. Our quasi-2-D approach bridges the gap between oversimplified 1-D models and computationally expensive 3-D simulations, providing a practical framework for investigating magma transport in silicic dikes.

Krafla volcano, located in Iceland's North Volcanic Zone, has been extensively studied using static gravity surveys since 1975. This study integrates measurements acquired between 2022 and 2023 with legacy data to produce a new gravity map of the area. We produce gravity gradient maps to delineate shallow subsurface density features and invert them to develop a 3D density model that images key subsurface structures. The combined survey coverage, extending from Fremrinámar to the Þeistareykir Volcanic Systems, reveals a range of features with contrasting densities which reflect the tectonic and volcanic processes that have shaped the area. Denser materials are associated with mafic intrusions along the principal faults related to caldera‐subsidence and ancillary faults within the caldera. In the Iceland Deep Drilling Project 1 (IDDP‐1) area, located inside the caldera, we identify a positive density anomaly at depths corresponding to magma encountered during drilling at 2.1 km depth. Negative gravity anomalies, indicative of relatively low‐density materials, are distributed within and outside the Krafla caldera. Within the caldera, these anomalies are interpreted as felsic intrusions and highly fractured geothermal zones. Beyond the caldera, a prominent negative density anomaly corresponds to the graben structure associated with the Húsavík‐Flatey Fault. Although strong near‐surface gravity anomalies and the spatial distribution of our gravity measurements limit the identification of deeper structures (>5 km), this study offers valuable insights into the distribution of the magmatic system features, as well as large tectonic characteristics in the area. These insights improve our understanding of magmatic and tectonic processes, volcanic hazards, and the future development of geothermal production in the area.

Geophysical, geochemical, and petrological studies indicate that magma storage and chemical differentiation can occur at various depths in the continental crust. Examination of exposed ancient magmatic systems reveal a continuum from a refractory lower crust to a silica-rich upper crust, prompting a debate on the mechanisms governing magma storage and differentiation in vertically extensive magmatic systems.Conceptual models for magma storage and chemical differentiation within the crust range from (i) a low crystallinity magma body, where differentiation occurs primarily by fractional crystallisation and partial melting of the crust, to (ii) a trans-crustal high crystallinity ‘mush’ reservoir consisting of a porous crystal matrix with melt in the pore space. Chemical differentiation in these mush reservoirs is driven by the compaction of the crystal matrix, buoyant upward reactive flow of melt, and partial melting of the surrounding crust.We use a numerical model to investigate magma storage and chemical differentiation in the continental crust. The model assumes a magmatic system sustained by the intrusion of mantle-derived basalt into the lower crust. We find that intrusion of basalt creates a high crystallinity reservoir in the lower crust. Reactive flow and compaction cause melt to accumulate at the top of the reservoir, creating a layer of low crystallinity, chemically differentiated magma. Buoyancy overpressure causes this magma to evacuate via dikes and we assume the magma intrudes the mid-crust to form a second high crystallinity mush reservoir. Reactive flow and compaction once again lead to the accumulation of low-crystallinity, evolved magma at the top of the reservoir that can evacuate and intrude the upper crust, driving volcanic eruptions or forming shallow plutons.The compositional evolution of magma as it ascends through the system is influenced by the flux of the parental basalt and the fertility of the crust. A higher basalt flux leads to warmer reservoirs that evacuate less evolved magma. Reservoirs in infertile crust that cannot melt are formed only of hot, parental magma, and also evacuate less evolved magma. Partial melting of fertile crust allows melt to percolate upwards and accumulate in cooler crust, leading to the evacuation of evolved (silicic) magma.In most of our simulation cases, the reservoirs remain as discrete bodies in the lower- and mid-crust, with magma transfer occurring via dikes; a trans-crustal reservoir does not form. The lower-crust reservoir evacuates mafic to intermediate magma, whilst the cooler mid-crust reservoir primarily evacuates evolved, silicic magma, consistent with vertical compositional trends in exposed, ancient magmatic systems. Both reservoirs comprise primarily low-melt fraction mush, consistent with geophysical imaging of contemporary systems. Layers of accumulated low-crystallinity magma are transient and likely too thin to resolve in geophysical data. The predicted volume, composition, and frequency of episodic magma intrusions into the upper crust are consistent with observed data from large volcanic eruptions. Our results suggest that reactive flow in multiple mush reservoirs controls magma storage and differentiation.

It is now generally accepted that large earthquakes can promote eruptions at nearby volcanoes. However, the prevalence of “triggered” eruptions, as well as the distance and timescale over which triggering occurs, remain unclear. Here, we use modern global earthquake and eruption records to compare volcanic eruption rates before and after large earthquakes with the time- averaged background eruption rate. We quantify the significance of observed deviations from the average eruption rate using Monte Carlo simulations. To integrate our findings with previous eruption triggering studies, we systematically vary the earthquake magnitudes we consider, as well as the distances and timescales used to calculate eruption rates. We also investigate the effects of earthquake depth and slip orientation. Overall, we find that post-earthquake eruption rates are around 1.25 times the average eruption rate within 750 km and one year following Mw ≥ 7 earthquakes, with above-average post-earthquake eruption rates possibly lasting for two to four years. By contrast, pre-earthquake eruption rates are around 0.9 times the average eruption rate within 750 km and182 days before Mw ≥7 earthquakes. Furthermore, deep earthquakes (≥7 0km) appear to more strongly affect eruption rates than shallow earthquakes, while earthquake slip orientation is also important. Further study of the relationships reported here represents a good opportunity to improve our understanding of tectono-magmatic relationships.

Many volcanoes show continuous but variable deformation over timescales of years to decades. Variations in uplift rate are typically interpreted as changes in magma supply rate and/or a viscoelastic response of the host rock. Here we conduct analogue experiments in the laboratory to represent the inflation of a silicic magma body at a constant volumetric flux, and measure the chamber pressure and resulting surface displacement field. We observe that dyke intrusions radiating from the magma body cause a decrease in the peak uplift rate, but do not significantly affect the spatial pattern of deformation or spatially averaged uplift rate. We identify 4 distinct phases: 1) elastic inflation of the chamber, 2) a gradual decrease in the rate of uplift and pressurisation, associated with the formation of visible cracks 3) propagation of a dyke by mode 1 failure at the crack tip and 4) a pressure decrease within the chamber. Phase 2 can be explained by either a) crack damage, which reduces the elastic moduli of the surrounding rock or b) magma filling pre-existing cracks. Thus these experiments provide alternative mechanisms to explain observed variations in uplift rate, with important implications for the interpretation of deformation patterns at volcanoes around the world.

Stratovolcanoes are common globally, with high altitude summit regions that are often glacier-clad and intersect the seasonal and perennial snow line. Explosive eruptions from stratovolcanoes can generate pyroclastic density currents (PDCs). When PDCs are emplaced onto and propagate over glacierised substrates, melt and steam are generated and incorporated into the flow, which can cause a transformation from hot, dry granular flow, to a water-saturated, sediment-laden flow, termed a lahar. Both PDCs and ice-melt lahars are highly hazardous due to their high energy during flow and long runout distances. Knowledge of the physics that underpin these interactions and the transformation to ice-melt lahar is extremely limited, preventing accurate descriptions within hazard models. To physically constrain the thermal interactions we conduct static melting experiments, where a hot granular layer was emplaced onto an ice substrate. The rate of heat transfer through the particle layer, melt and steam generation were quantified. Experiments revealed systematic increases in melt and steam with increasing particle layer thicknesses and temperatures. We also present a one-dimensional numerical model for heat transfer, calibrated against experiment data, capable of accurately predicting temperature and associated melting. Furthermore, we present similarity solutions for early-time melting which are used to benchmark our numerical scheme, and to provide rapid estimates for meltwater flux hydrographs. These data are vital for predicting melt volume and incorporation into PDCs required to facilitate the transformation to and evolution of ice-melt lahars.

SUMMARY Dyke propagation is a mechanism for more rapid ascent of felsic magmas through the crust than is possible via diapirs or percolative flow. As it ascends, the magma undergoes complex physical and chemical transformations induced by decompression and cooling. These processes dramatically change the magma density and viscosity, which in turn affect magma ascent rate and the depth at which the dyke arrests. We present a mathematical model of dyke propagation for silicic magmas taking into account the presence of multiple volatile species (H2O and CO2), bubble growth, heat advection and loss, crystallization and latent heat release. We consider conditions for dykes associated with porphyry ore deposits, which may represent an end-member in rapid ascent of felsic magmas from depth. In particular, we simulate the propagation of dykes launched from a deep (900 MPa), volatile-saturated magma source, testing the effects of the magma H2O/CO2 content, temperature and mass on its ascent rate and final emplacement depth. The model predicts short ascent times (hours to days), with a large increase in viscosity at shallow depth, leading to stagnation and solidification of the dyke. Higher initial water content, higher temperature and larger mass of the magma in the dyke promote faster propagation and shallower arrest. Volatile loss from ascending magma remains limited until the stagnation depth, providing a potential mechanism for transfer of deep volatiles to hypabyssal blind intrusions associated with porphyry ore deposits. Our findings are applicable to the problem of silicic magma ascent through the crust more generally.

Large-scale, explosive volcanic eruptions are one of the Earth’s most hazardous natural phenomena. We demonstrate that their size, frequency, and composition can be explained by processes in long-lived, high-crystallinity source reservoirs that control the episodic creation of large volumes of eruptible silicic magma and its delivery to the subvolcanic chamber where it is stored before eruption. Melt percolates upward through the reservoir and accumulates a large volume of low-crystallinity silicic magma which remains trapped until buoyancy causes magma-driven fractures to propagate into the overlying crust, allowing rapid magma transfer from the reservoir into the chamber. Ongoing melt percolation in the reservoir accumulates a new magma layer and the process repeats. Our results suggest that buoyancy, rather than crystallinity, is the key control on magma delivery from the source reservoir. They identify an optimum reservoir size for the largest silicic eruptions that is consistent with data from natural systems and explain why larger magnitude eruptions are not observed on Earth.

The Krafla Volcanic System (KVS) in the Northern Volcanic Zone (NVZ) in Iceland last erupted between 1975 and 1984, during an eruptive period called “the Krafla Fires”. The KVS is composed of a restless caldera, an array of scoria cones along a fissure swarm and is among the best-studied volcanic systems due to the exploitation of its geothermal potential. In 2009, the Icelandic Deep Drilling Project (IDDP) encountered a shallow rhyolitic magma body at 2.1 km depth beneath the caldera. To date, no geophysical method has been able to image this magma body at Krafla within the top 4 km of the crust. Here we present new micro-gravity data collected in June and July 2022 across a 14-station network of benchmarks in the KVS. Micro-gravimetry is a relative method that records changes in gravity between a reference and a series of benchmarks over both space and time to investigate subsurface mass or density changes via time-series analysis and modelling. Our 2022 survey highlights negative gravity differences of benchmarks located in the centre of the caldera with respect to a reference located to the south and outside the caldera. The most negative values are found in its eastern part. Positive gravity differences can be found south of the southern caldera wall along a set of past eruptive fissures. The next steps in data processing include data reduction for deformation effects to link the new data to previous joint deformation and micro-gravity surveys conducted at the KVS since 1965. This should enable us to quantify the long-term evolution of the KVS over more than 50 years providing unprecedented insights into its inner workings.

Knowledge of the architecture of active magmatic systems is important for both volcanic hazard assessment and evaluating potential geothermal energy production and metals recovery from magmatic fluids. Increasingly, magmatic systems are imaged using the magnetotelluric method to detect electrically conductive partial melt and/or saline magmatic fluid reservoirs. We review recent magnetotelluric studies at eight Andean volcanoes, revealing electrical conductivity anomalies with variable magnitudes and locations. Six of the studied volcanoes exhibit three main electrical conductivity anomalies, located at shallow (<3 km), intermediate (≈5 km), and deep (>10 km) depths. The shallow anomalies are often thin and laterally extensive, consistent with clay cap alteration layers, while the deep anomalies are generally interpreted as partial melt reservoirs. The intermediate depth anomalies, although also often attributed to partial melt, have less clear origins. By analysing laboratory-derived electrical conductivity relationships, we show that the intermediate depth anomalies are generally most consistent with saline magmatic fluids stored in porous rock. However, other geophysical and petrological data suggest that localised partial melt also exists at intermediate depths. Therefore, the intermediate depth anomalies likely represent mixed melt and saline magmatic fluid systems, such as those responsible for forming magmatic-hydrothermal alteration zones and copper porphyry deposits. At individual volcanoes, refining the generalised three layer model proposed here by using additional geophysical or petrological data is key to constraining the resources and/or hazard potential of the magmatic system.

Pyroclastic density current (PDC) interactions with ice are common at high altitude and latitude stratovolcanoes. When PDCs propagate over ice, melt and steam are generated. The incorporation of melt and steam into PDCs can alter the flow dynamics by reducing friction at the particle-ice interface and between individual particles. Melt incorporation can also transform a PDC into an ice-melt lahar. The hazardous and temporally unpredictable nature of these flows limits field observations. Conceptual models of PDC-ice interactions for hazard assessment and modelling exist, but quantifications of the microscale physical processes that underpin these interactions are limited. We use experiments to characterise the melting and friction reduction that occur when PDCs are emplaced onto ice. In experiment set one, a heated particle layer was rapidly emplaced onto a horizontal ice layer contained within an insulated beaker 7.3 cm in diameter. The particle types used were glass ballotini, crushed pumice, and Ruapehu PDC samples, covering a diverse range of grain characteristics. The particle layer was varied in thickness up to 45 mm and across temperatures up to 700 °C. In each experiment, the mass of melt and steam were quantified, and the time evolution of temperature through the particle layer was measured. Across all particle types, increasing particle layer mass (therefore layer thickness) and temperature increased melt and steam production. However, Ruapehu and pumice melt masses showed greater sensitivity than ballotini to particle temperature for any given layer thickness. Conversely, steam production was greater for the ballotini for any given layer thickness and was more sensitive to ballotini particle temperature. Localised steam escape, fluidisation, capillary action, and particle sinking, were observed to varying extents in the experiments. These phenomena caused melt to be incorporated into the particle layer. The rate of increase in melt generation decreases with increasing particle layer thickness. This is due to increasing steam production, the increasing temperature of incorporated meltwater, energy losses to the atmosphere, and alterations to the bulk particle diffusivity. Experiment set two characterised the mobility of particles over frozen and non-frozen substrates. Pumice and Ruapehu particles of varying temperature and layer thickness were poured into a 4.5 cm diameter alumina tube, which was rapidly lifted, allowing the particles to radially spread over the substrate. This configuration has been widely studied in experiments on granular flow mobility. The initial and final aspect ratios of the particle layer were measured, and conform to a power-law form previously interpreted as showing that frictional interactions are only important in the final stages of flow emplacement. Enhanced particle layer mobility over ice was only observed for Ruapehu particles above 400 °C, which we interpret to be due to fluidisation of the particles by rising steam. This is consistent with experiment set one, where Ruapehu particles produced more steam than pumice, and were often fluidised above 400 °C. Experimental data will be used to calibrate surface flow hazard models for PDC runout and lahar generation, enabling prediction of PDC-ice interaction hazards. These models will be tested at Mt. Ruapehu, New Zealand.

When volcanic eruptions involve interaction with external water (hydrovolcanism), the result is an ash-rich and energetic volcanic plume, as illustrated dramatically by the January 2022 Tonga eruption. The origin of the high explosive energy of these events remains an important question. We investigate this question by studying Prince Rupert’s Drops (PRDs)—tadpole-shaped glass beads formed by dripping molten glass into water—which have long fascinated materials scientists because the great strength of the head contrasts with the explosivity of the metastable interior when the tail is broken. We show that the fragment size distribution (FSD) produced by explosive fragmentation changes systematically with PRD fragmentation in air, water, and syrup. Most FSDs are fractal over much of the size range, scaling that can be explained by the repeated fracture bifurcation observed in three-dimensional images from microcomputed tomography. The shapes of constituent fragments are determined by their position within the original PRD, with platey fragments formed from the outer (compressive) shell and blocky fragments formed by fractures perpendicular to interior voids. When molten drops fail to form PRDs, the glass disintegrates by quench granulation, a process that produces fractal FSDs but with a larger median size than explosively generated fragments. Critically, adding bubbles to the molten glass prevents PRD formation and promotes quench granulation, suggesting that granulation is modulated by heterogeneous stress fields formed around the bubbles during sudden cooling and contraction. Together, these observations provide insight into glass fragmentation and potentially, processes operating during hydrovolcanism.

The scarcity of historical rhyolite eruptions means that volcanological and petrological studies of past eruptions are a key tool for assessing the potential for future hazardous activity and improving interpretations of unrest signals. For the last 18 ky, the Laguna del Maule (LdM) volcanic complex in Chile has erupted primarily rhyolites but with differing magma compositions, eruption styles and eruptive volumes. Rapid surface uplift and episodic seismic activity at LdM over the last two decades has emphasized the need to understand both the recent evolution of the magmatic system and the most likely future eruption scenarios. Using mineral composition, geothermobarometry and MELTS modeling, we assess the influence of the magma storage and ascent conditions on the magnitude and styles of three LdM rhyolites. Magmas of the first and largest Plinian-ignimbrite eruption (LdM: rdm; >17 km(3) DRE) form a distinct mineral assemblage with An(2)(8-)(60) plagioclase, amphibole and quartz derived from a magma plumbing system over a large pressure range (90-350 MPa). We suggest that the rdm eruption was triggered by magma recharge and overpressure within a similar to 90 MPa magma chamber of high-silica (>76.5 wt.% SiO2) rhyolitic melt. The rdm eruption appears to have reset the LdM storage conditions, such that subsequent rhyolite eruptions have been smaller (<4 km(3) DRE), involved less evolved melt (<75.8 wt.% SiO2) and produced (sub)Plinian-effusive eruption sequences. Of these, we have studied the earliest (Los Espejos: rle) and the most recent (Las Nieblas: rln). Both have mineral assemblages with An(20)(-3)(4) plagioclase and biotite formed in water-saturated crystal-poor rhyolite batches over limited pressure ranges (Delta P similar to 50 MPa) within a rhyodacitic mush reservoir. We suggest that rle and rln were triggered by injection of water-rich magma and volatile exsolution; the transition to effusive behavior occurred because of changes in magma temperature (rle) or magma decompression rate (rln). Similarities in the limited pressure range and high volatile content inferred for the current low-density magma body with the petrologically determined source conditions of magmas that fed past subplinian-effusive eruptions suggest that the next LdM eruption will be similar.

Detailed microstructural analysis of three basaltic sills of the Little Minch Sill Complex demonstrates that convection leaves a detectable signature in fully solidified bodies. The presence of dense clusters of equant grains of olivine and clinopyroxene in the central parts of sills can only be accounted for if they formed and were enlarged while suspended in convecting magma, with delayed settling to the sill floor. An associated stratigraphic invariance of plagioclase grain shape is consistent with growth while suspended in convecting magma. These microstructural indicators demonstrate that convection during solidification was vigorous and long-lived in the 135-m-thick picrodolerite-crinanite unit (PCU) of the composite Shiant Isles Main sill and vigorous and likely short-lived in the PCU of the composite Creagan Iar sill. In contrast, convection in the Meall Tuath sill was weak and short-lived: plagioclase grain shape in this sill varies with stratigraphic height, indicative of primarily in situ nucleation and growth at the magma-mush interface, while olivine and clinopyroxene were kept suspended in the overlying convecting magma. The magma in all three sills fractionated during solidification, permitting convection driven by the instability of an upper thermal boundary layer. The comparative vigour and longevity of convection in the Shiant Isles Main sill and the Creagan Iar sill was due to their emplacement above an earlier, still-hot, intrusion, resulting in highly asymmetric cooling.

The breaking apart of magma into fragments is intimately related to the eruptive style and thus the nature and footprint of volcanic hazards. The size and shape distributions of the fragments, in turn, affect the efficiency of heat transfer within pyroclastic plumes and currents and the settling velocity, and so the residence time, of particles in the atmosphere. Fundamental work relating the glass transition to the fragmentation of magmas remains at the heart of conceptual and numerical models of volcanic eruptions. Current fragmentation criteria, however, do not predict the sizes and shapes of the resulting fragments, or fully account for the multiphase nature of magmas or ways in which magma can break in a fluidal manner or by thermal stress. The pulsatory, non-steady state nature of some eruptions, and related interactions with these fragmentation criteria, also requires further investigation. Here, we briefly review some recent advances in the field of magma fragmentation and provide a perspective on how integrated field, experimental and numerical modelling studies can address key outstanding challenges.

Generation of silicic magmas leads to emplacement of granite plutons, huge explosive volcanic eruptions and physical and chemical zoning of continental and arc crust1–7. Whereas timescales for silicic magma generation in the deep and middle crust are prolonged8, magma transfer into the upper crust followed by eruption is episodic and can be rapid9–12. Ages of inherited zircons and sanidines from four Miocene ignimbrites in the Central Andes indicate a gap of 4.6 Myr between initiation of pluton emplacement and onset of super-eruptions, with a 1-Myr cyclicity. We show that inherited zircons and sanidine crystals were stored at temperatures <470 °C in these plutons before incorporation in ignimbrite magmas. Our observations can be explained by silicic melt segregation in a middle-crustal hot zone with episodic melt ascent from an unstable layer at the top of the zone with a timescale governed by the rheology of the upper crust. After thermal incubation of growing plutons, large upper-crustal magma chambers can form in a few thousand years or less by dike transport from the hot-zone melt layer. Instability and disruption of earlier plutonic rock occurred in a few decades or less just before or during super-eruptions. Analysis of inherited zircons and sanidines from Miocene ignimbrites in the Central Andes shows that plutons were emplaced for up to 4 million years prior to onset of volcanism and that disruption of plutonic rock occurs a few decades or less just before or during super-eruptions.

Recent studies have shown that large tectonic earthquakes are capable of triggering volcanic eruptions (i.e. increasing the number of eruptions within a defined time period) up to hundreds of kilometres away. However, the prevalence of eruption triggering is less clear, with findings ranging from little evidence for triggered eruptions, to a fourfold increase in the number of eruptions following nearby large earthquakes. Some of this variability is likely due to differences in definitions of what constitutes a triggered volcanic eruption, including a lack of consensus on the maximum distance and time lag between an earthquake and a triggered volcanic eruption, the minimum magnitude of earthquake considered, and how aftershocks are incorporated into the analysis. A further source of variability arises from the different datasets used, including regional versus global studies, and the inclusion of incomplete earthquake and eruption records from before the modern instrumental era. To help address these issues, we provide a comprehensive statistical study of how large earthquakes affect volcanic eruption rates, using complete and unbiased global datasets spanning 1960-2021. We take a systematic approach to investigating how parameters such as the maximum distance and time lag between earthquake-eruption pairs, the minimum earthquake magnitude considered, and the declustering of aftershocks affects the results. We also investigate how previously unstudied earthquake parameters such as source depth and mechanism affect the prevalence of eruption triggering. Our results are placed in statistical context through the use of Monte Carlo simulations using randomised earthquake and eruption catalogues. Preliminary results indicate that, contrary to a previous focus on large subduction megathrust earthquakes, deep normal faulting earthquakes have the greatest eruption triggering tendency. However, when compared with randomised earthquake and eruption catalogues, the overall statistical significance of observed eruption triggering is fairly low.