Explosive expansion of a slowly decompressed magma analogue: Evidence for delayed bubble nucleation

While ascending in the plumbing system of volcanoes, magma undergoes decompression at rates spanning several orders of magnitude and set by a number of factors internal and external to the volcano. Slow decompression generally results in an effusive or mildly explosive expansion of the magma, but counterexamples of sudden switches from effusive to explosive eruptive behavior have been documented at various volcanoes worldwide. The mechanisms involved in this behavior are currently debated, in particular for basaltic magmas. Here, we explore the interplay between decompression rate and vesiculation vigor by decompressing a magma analogue obtained by dissolving pine resin into acetone in varying proportions. Analogue experiments allow direct observations of the processes of bubble nucleation and growth, flow dynamics, and fragmentation that is not currently possible with magmatic systems. Our mixtures contain solid particles, and upon decompression, nucleation of acetone bubbles is observed. We find that mixtures with a high acetone content, containing smaller and fewer solid particles, experience strong supersaturation and fragment under very slow decompressions, despite having low viscosity, while mixtures with lower acetone content, with more and larger solid particles, degas efficiently without fragmentation. We interpret our results in terms of delayed bubble nucleation due to a lack of efficient nucleation sites. We discuss how a similar mechanism might induce violent, explosive expansion in volatile‐rich and poorly crystalline low‐silica magmas, by analogy with the behavior of rhyolitic magmas

Sloshing of a bubbly magma reservoir as a mechanism of triggered eruptions

open4siLarge earthquakes sometimes activate volcanoes both in the near field as well as in the far field. One possible explanation is that shaking may increase the mobility of the volcanic gases stored in magma reservoirs and conduits. Here experimentally and theoretically we investigate how sloshing, the oscillatory motion of fluids contained in a shaking tank, may affect the presence and stability of bubbles and foams, with important implications for magma conduits and reservoirs. We adopt this concept from engineering: severe earthquakes are known to induce sloshing and damage petroleum tanks. Sloshing occurs in a partially filled tank or a fully filled tank with density-stratified fluids. These conditions are met at open summit conduits or at sealed magma reservoirs where a bubbly magma layer overlays a newly injected denser magma layer. We conducted sloshing experiments by shaking a rectangular tank partially filled with liquids, bubbly fluids (foams) and fully filled with density-stratified fluids; i.e., a foam layer overlying a liquid layer. In experiments with foams, we find that foam collapse occurs for oscillations near the resonance frequency of the fluid layer. Low viscosity and large bubble size favor foam collapse during sloshing. In the layered case, the collapsed foam mixes with the underlying liquid layer. Based on scaling considerations, we constrain the conditions for the occurrence of foam collapse in natural magma reservoirs. We find that seismic waves with lower frequencies < 1 Hz, usually excited by large earthquakes, can resonate with magma reservoirs whose width is > 0.5 m. Strong ground motion > 0.1 m s− 1 can excite sloshing with sufficient amplitude to collapse a magma foam in an open conduit or a foam overlying basaltic magma in a closed magma reservoir. The gas released from the collapsed foam may infiltrate the rock or diffuse through pores, enhancing heat transfer, or may generate a gas slug to cause a magmatic eruption. The overturn in the magma reservoir provides new nucleation sites which may help to prepare a following/delayed eruption. Mt. Fuji erupted 49 days after the large Hoei earthquake (1707) both dacitic and basaltic magmas. The eruption might have been triggered by magma mixing through sloshing.Acknowledgments We thank Ralf Bauz, Stefan Mikulla, and Peter Neuendorf for their help in the GFZ bubble-lab, and Matthias Rosenau for providing the glucose syrup. Helpful and careful reviews by David P. Hill and an anonymous reviewer, and the journal editor, Lionel Wilson, and Maggie Mangan are deeply appreciated. This work was supported by KAKENHI 24681035, the European Union through the ERC StG project CCMP-POMPEI, grant number 230583, and the European Supersite project MED-SUV, grant number 308665.openNamiki, Atsuko; Rivalta, Eleonora; Woith, Heiko; Walter, Thomas R.Namiki, Atsuko; Rivalta, Eleonora; Woith, Heiko; Walter, Thomas R

Graben formation and dike arrest during the 2009 Harrat Lunayyir dike intrusion in Saudi Arabia: Insights from InSAR, stress calculations and analog experiments

Detailed spatial and temporal accounts of propagating dikes from crustal deformation data, including their interplay with faulting, are rare, leaving many questions about how this interplay affects graben formation and the arrest of dikes unanswered. Here we use interferometric synthetic aperture radar (InSAR) observations, stress calculations, and analog experiments to investigate the interaction between an intruding dike and normal faulting during the 2009 Harrat Lunayyir dike intrusion in western Saudi Arabia. We generated five displacement maps from InSAR data to unravel the temporal evolution of deformation covering the majority of the intrusion. We find that the observed surface displacements can be modeled by a ~2 m thick dike with an upper tip ~2\u2009km below the surface on 16 May 2009, 4\u2009weeks after the onset of seismic unrest. In the following three days, the dike propagated to within ~1\u2009km of the surface with graben\u2010bounding normal faulting dominating the near\u2010field deformation. The volume of the dike doubled between mid\u2010May and mid\u2010June. We carried out analog experiments that indicate that the wedge\u2010shaped graben grew outward with the faulting style changing progressively from normal faulting to oblique. Coulomb failure stress change calculations show that the intruding dike caused two zones of shallow horizontal tension on both sides of the dike, producing two zones of fissuring and normal faulting at the surface. In return, the faulting provoked compression around the upper tip of the dike, holding back its vertical propagation

Detection of Spatial and Temporal Stress Changes During the 2016 Central Italy Seismic Sequence by Monitoring the Evolution of the Energy Index

We consider approximately 23,000 microearthquakes that occurred between 2005 and 2016 in central Italy to investigate the crustal strength before and after the three largest earthquakes of the 2016 seismic sequence (i.e., the Mw 6.2, 24 August 2016 Amatrice, the Mw 6.1, 26 October 2016 Visso, and the Mw 6.5, 30 October 2016 Norcia earthquakes). We monitor the spatiotemporal deviations of observed radiated energy, ES, with respect to theoretical values, ESt, derived from a scaling model between ES and M0 calibrated for background seismicity in central Italy. These deviations, defined here as Energy Index (EI), allow us to identify in the years following the Mw 6.1, 2009 L’Aquila earthquake a progressive evolution of the dynamic properties of microearthquakes and the existence of high EI patches close to the Amatrice earthquake hypocenter. We show the existence of a crustal volume with high EI even before the Mw 6.5 Norcia earthquake. Our results agree with the previously suggested hypothesis that the Norcia earthquake nucleated at the boundary of a large patch, highly stressed by the two previous mainshocks of the sequence. We highlight the mainshocks interaction both in terms of EI and of the mean loading shear stress associated to microearthquakes occurring within the crustal volumes comprising the mainshock hypocenters. Our study shows that the dynamic characteristics of microearthquakes can be exploited as beacons of stress change in the crust and thus be exploited to monitor the seismic hazard of a region and help to intercept the preparation phase of large earthquakes

The impact of unloading stresses on post-caldera magma intrusions

Calderas represent morphological depressions several kilometers in diameter, and the unloaded crustal stresses they produce can form rapidly (e.g. Pinatubo, 1990) or slowly (e.g. Hawaii, 2018). Active calderas are known as sites of persistent magma intrusions, and yet the dynamics of their shallow plumbing system is not well constrained. We use scaled laboratory experiments to study how experimental intrusions are created beneath a caldera by injecting dyed water (magma analogue) into the base of an elastic gelatin solid (crust analogue) with a cylindrical cavity in its surface to mimic a caldera-like topography. The evolving dike geometry and stress field were qualitatively determined using polarized light, and digital image correlation allowed the incremental and total strain to be quantified by tracking passive-tracer particles in the gelatin that fluoresced in a thin 2D vertical laser sheet. Our results show that the unloaded stress field from a caldera can cause a divergence of vertical dikes, and leads to circumferential dikes and cone sheets. When the caldera was large the initially vertical dike became arrested, then grew laterally via circumferentially-propagating en echelon segments; these eventually joined to complete a cone sheet that was parallel to, but extended outside and beneath, the large caldera. When the caldera was small, a circumferential dike erupted, producing a short fissure which was outside, but parallel to, the caldera. We suggest that the distinct curved geometry, velocity, strain and stress characteristics of circumferential dikes and cone sheets can be used to interpret the origin and growth of post-caldera magmatism and the likelihood of eruption in caldera systems

Seismic Amplitude Ratio Analysis of the 2014-2015 Bárðarbunga-Holuhraun Dike Propagation and Eruption

Magma is transported in brittle rock through dikes and sills. This movement may be accompanied by the release of seismic energy that can be tracked from the Earth's surface. Locating dikes and deciphering their dynamics is therefore of prime importance in understanding and potentially forecasting volcanic eruptions. The Seismic Amplitude Ratio Analysis (SARA) method aims to track melt propagation using the amplitudes recorded across a seismic network without picking the arrival times of individual earthquake phases. This study validates this methodology by comparing SARA locations (filtered between 2 and 16 Hz) with the earthquake locations (same frequency band) recorded during the 2014–2015 Bár urn:x-wiley:jgrb:media:jgrb52508:jgrb52508-math-0003arbunga‐Holuhraun dike intrusion and eruption in Iceland. Integrating both approaches also provides the opportunity to investigate the spatiotemporal characteristics of magma migration during the dike intrusion and ensuing eruption. During the intrusion SARA locations correspond remarkably well to the locations of earthquakes. Several exceptions are, however, observed. (1) A low‐frequency signal was possibly associated with a subglacial eruption on 23 August. (2) A systematic retreat of the seismicity was also observed to the back of each active segment during stalled phases and was associated with a larger spatial extent of the seismic energy source. This behavior may be controlled by the dike's shape and/or by dike inflation. (3) During the eruption SARA locations consistently focused at the eruptive site. (4) Tremor‐rich signal close to ice cauldrons occurred on 3 September. This study demonstrates the power of the SARA methodology, provided robust site amplification; Quality Factors and seismic velocities are available.The authors thank both reviewers, the Associate Editor, and the Editor for their insightful comments and suggestions that greatly improved this study. Seismometers were borrowed from the Natural Environment Research Council (NERC) SEIS-UK facility (loans 968 and 1022), with funding by research grants from the NERC and the European Community’s Seventh Framework Programme grant 308377 (Project FUTUREVOLC), and graduate studentships from the NERC. C. Caudron benefited from a Fondation Wiener Anspach postdoctoral fellowship, then from a FNRS Chargé de Recherche postdoctoral grantPeer Reviewe

Cyclical geothermal unrest as a precursor to Iceland’s 2021 Fagradalsfjall eruption

Understanding and constraining the source of geodetic deformation in volcanic areas is an important component of hazard assessment. Here, we analyse deformation and seismicity for one year before the March 2021 Fagradalsfjall eruption in Iceland. We generate a high-resolution catalogue of 39,500 earthquakes using optical cable recordings and develop a poroelastic model to describe three pre-eruptional uplift and subsidence cycles at the Svartsengi geothermal field, 8 km west of the eruption site. We find the observed deformation is best explained by cyclic intrusions into a permeable aquifer by a fluid injected at 4 km depth below the geothermal field, with a total volume of 0.11 ± 0.05 km3 and a density of 850 ± 350 kg m–3. We therefore suggest that ingression of magmatic CO2 can explain the geodetic, gravity and seismic data, although some contribution of magma cannot be excluded

Massive earthquake swarm driven by magmatic intrusion at the Bransfield Strait, Antarctica

An earthquake swarm affected the Bransfield Strait, Antarctica, a unique rift basin in transition from intra-arc rifting to ocean spreading. The swarm, counting ~85,000 volcanotectonic earthquakes since August 2020, is located close to the Orca submarine volcano, previously considered inactive. Simultaneously, geodetic data reported up to ~11 cm northwestward displacement over King George Island. We use a broad variety of geophysical data and methods to reveal the complex migration of seismicity, accompanying the intrusion of 0.26-0.56 km³ of magma. Strike-slip earthquakes mark the intrusion at depth, while shallower normal faulting the ~20 km long lateral growth of a dike. Seismicity abruptly decreased after a Mw 6.0 earthquake, suggesting the magmatic dike lost pressure with the slipping of a large fault. A seafloor eruption is likely, but not confirmed by sea surface temperature anomalies. The unrest documents episodic magmatic intrusion in the Bransfield Strait, providing unique insights into active continental rifting

Aseismic transient driving the swarm-like seismic sequence in the Pollino range, Southern Italy

Tectonic earthquake swarms challenge our understanding of earthquake processes since it is difficult to link observations to the underlying physical mechanisms and to assess the hazard they pose. Transient forcing is thought to initiate and drive the spatio-temporal release of energy during swarms. The nature of the transient forcing may vary across sequences and range from aseismic creeping or transient slip to diffusion of pore pressure pulses to fluid redistribution and migration within the seismogenic crust. Distinguishing between such forcing mechanisms may be critical to reduce epistemic uncertainties in the assessment of hazard due to seismic swarms, because it can provide information on the frequency–magnitude distribution of the earthquakes (often deviating from the assumed Gutenberg–Richter relation) and on the expected source parameters influencing the ground motion (for example the stress drop). Here we study the ongoing Pollino range (Southern Italy) seismic swarm, a long-lasting seismic sequence with more than five thousand events recorded and located since October 2010. The two largest shocks (magnitude Mw = 4.2 and Mw = 5.1) are among the largest earthquakes ever recorded in an area which represents a seismic gap in the Italian historical earthquake catalogue. We investigate the geometrical, mechanical and statistical characteristics of the largest earthquakes and of the entire swarm. We calculate the focal mechanisms of the Ml > 3 events in the sequence and the transfer of Coulomb stress on nearby known faults and analyse the statistics of the earthquake catalogue. We find that only 25 per cent of the earthquakes in the sequence can be explained as aftershocks, and the remaining 75 per cent may be attributed to a transient forcing. The b-values change in time throughout the sequence, with low b-values correlated with the period of highest rate of activity and with the occurrence of the largest shock. In the light of recent studies on the palaeoseismic and historical activity in the Pollino area, we identify two scenarios consistent with the observations and our analysis: This and past seismic swarms may have been ‘passive’ features, with small fault patches failing on largely locked faults, or may have been accompanied by an ‘active’, largely aseismic, release of a large portion of the accumulated tectonic strain. Those scenarios have very different implications for the seismic hazard of the area