Banding in the Margins of Basaltic Dykes Indicates Pulsatory Propagation During Emplacement

Basaltic fissure eruptions, which are the most common type of eruption on Earth, are fed by dykes which mediate magma transport through the crust. Dyke propagation processes are important because they determine the geometry of the transport pathway and the nature of any geophysical signals associated with magma ascent. Here, we investigate small‐scale (mm–cm wide) banding features at the margins of dykes in the Teno Massif (Tenerife, Spain) and the Columbia River Basalt Province (CRBP) (USA). Similar marginal bands have been reported for dykes in numerous localities around the world. Dyke margins record valuable information about propagation because they are the first material to solidify against the host rock at the propagating dyke tip. We find that the marginal bands are defined by cyclic variations in phenocryst concentration and vesicularity, and we infer that these cyclic variations in texture are a product of cyclic variations in magma flow rates and pressures within the dyke tip. This indicates that dyke emplacement occurs in pulses, with propagation repeatedly hindered by the rapid cooling and solidification of magma in the narrow dyke tip. Using a 1D conduction model, we estimate the time taken for each band to cool and solidify, which provides a timescale of several minutes to tens of minutes for the pulses. The occurrence of similar bands in various volcanic settings suggests that pulsatory propagation is a common, if not ubiquitous, process associated with dyke emplacement

The α–β phase transition in volcanic cristobalite

Cristobalite is a common mineral in volcanic ash produced from dome-forming eruptions. Assessment of the respiratory hazard posed by volcanic ash requires understanding the nature of the cristobalite it contains. Volcanic cristobalite contains coupled substitutions of Al3+ and Na+ for Si4+; similar co-substitutions in synthetic cristobalite are known to modify the crystal structure, affecting the stability of the [alpha] and [beta] forms and the observed transition between them. Here, for the first time, the dynamics and energy changes associated with the [alpha]-[beta] phase transition in volcanic cristobalite are investigated using X-ray powder diffraction with simultaneous in situ heating and differential scanning calorimetry. At ambient temperature, volcanic cristobalite exists in the [alpha] form and has a larger cell volume than synthetic [alpha]-cristobalite; as a result, its diffraction pattern sits between ICDD [alpha]- and [beta]-cristobalite library patterns, which could cause ambiguity in phase identification. On heating from ambient temperature, volcanic cristobalite exhibits a lower degree of thermal expansion than synthetic cristobalite, and it also has a lower [alpha]-[beta] transition temperature (~473 K) compared with synthetic cristobalite (upwards of 543 K); these observations are discussed in relation to the presence of Al3+ and Na+ defects. The transition shows a stable and reproducible hysteresis loop with [alpha] and [beta] phases coexisting through the transition, suggesting that discrete crystals in the sample have different transition temperatures

A validated numerical model for the growth and resorption of bubbles in magma

The rate and timing of bubble growth in magma is an important control on eruption style, determining whether or not magma fragments to produce an explosive eruption. Bubbles nucleate, grow, shrink, and de-nucleate in magma in response to changes in pressure and temperature, and these changes may be recorded in the spatial distribution and speciation of water 'frozen into' the glass in eruptive products. Accurate modelling of growth and resorption is therefore essential both for forward modelling of eruptive processes, and for inverse modelling to reconstruct pre-eruptive history. We present the first experimentally-validated numerical model for bubble growth and resorption in magma. The model includes the kinetics of speciation, allows for arbitrary temperature and pressure pathways, and accounts for the impact of spatial variations in water content on diffusivity and viscosity. We validate the model against three sets of data. (1) Continuous vesicularity-time data collected using optical dilatometry and in-situ synchrotron-source x-ray tomography of natural and synthetic magma during thermally-induced vesiculation and resorption at magmatic temperatures and ambient pressure. This represents approximately isobaric bubble growth and resorption under disequilibrium conditions. (2) Final vesicularity data from decompression experiments at magmatic temperatures and pressures. This represents isothermal, decompression-driven bubble growth from equilibrium to strongly disequilibrium conditions. (3) Speciation data from diffusion-couple experiments on synthetic haplogranites at magmatic temperatures and pressures. The numerical model closely reproduces all experimental data, providing validation against equilibrium and disequilibrium bubble growth/resorption and speciation scenarios. The validated model can be used to predict the growth and resorption of bubbles, and associated changes in magma properties, for arbitrary eruption pathways. It can also be used to reconstruct pressure-temperature-time pathways from textures and volatile contents of eruptive products. This will open up new ways of accessing the dynamics of magma ascent and eruption in unobserved volcanic eruptions

A reappraisal of explosive–effusive silicic eruption dynamics: syn-eruptive assembly of lava from the products of cryptic fragmentation

Silicic volcanic eruptions range in style from gently effusive to highly explosive, and may switch style unpredictably during a single eruption. Direct observations of subaerial rhyolitic eruptions (Chaiten 2008, Cordón Caulle 2011–2012, Chile) challenged long-standing paradigms of explosive and effusive eruptive styles and led to the formulation of new models of hybrid activity. However, the processes that govern such hybrid explosive–effusive activity remain poorly understood. Here, we bring together observations of the well-studied 2011–2012 Cordón Caulle eruption with new textural and petrologic data on erupted products, and video and still imagery of the eruption. We infer that all of the activity – explosive, effusive, and hybrid – was fed by explosive fragmentation at depth, and that effusive behaviour arose from sticking and sintering, in the shallow vent region, of the clastic products of deeper, cryptic fragmentation. We use a scaling approach to determine that there is sufficient time available, during emplacement, for diffusive pyroclast degassing and sintering to produce a degassed plug that occludes the shallow conduit, feeding clastogenic, apparently effusive, lava-like deposits. Based on evidence from Cordón Caulle, and from other similar eruptions, we further argue that hybrid explosive–effusive activity is driven by episodic gas-fracking of the occluding lava plug, fed by the underlying pressurized ash- and pyroclast-laden region. The presence of a pressurized pocket of ash-laden gas within the conduit provides a mechanism for generation of harmonic tremor, and for syn-eruptive laccolith intrusion, both of which were features of the Cordón Caulle eruption. We conclude that the cryptic fragmentation models is more consistent with available evidence than the prevailing model for effusion of silicic lava that assume coherent non-fragmental rise of magma from depth to the surface without wholesale explosive fragmentation

A model for permeability evolution during volcanic welding

Volcanic ash and pyroclasts can weld when deposited hot by pyroclastic density currents, in near-vent fall deposits, or in fractures in volcano interiors. Welding progressively decreases the permeability of the particle packs, influencing a range of magmatic and volcanic processes, including magma outgassing, which is an important control on eruption dynamics. Consequently, there is a need for a quantitative model for permeability evolution during welding of ash and pyroclasts under the range of conditions encountered in nature. Here we present in situ experiments in which hydrous, crystal-free, glassy pyroclasts are imaged via x-ray tomography during welding at high temperature. For each 3D dataset acquired, we determine the porosity, Darcian gas permeability, specific surface area, and pore connectivity. We find that all of these quantities decrease as a critical percolation threshold is approached. We develop a constitutive mathematical model for the evolution of permeability in welding volcanic systems based on percolation theory, and validate the model against our experimental data. Importantly, our model accounts for polydispersivity of the grainsize in the particle pack, the pressures acting on the pack, and changes in particle viscosity arising from degassing of dissolved H2O during welding. Our model is theoretically grounded and has no fitting parameters, hence it should be valid across all magma compositions. The model can be used to predict whether a cooling pyroclast pack will have sufficient time to weld and to degas, the scenarios under which a final deposit will retain a permeable network, the timescales over which sealing occurs, and whether a welded deposit will have disequilibrium or equilibrium H2O content. A user-friendly implementation of the model is provided

The initial phase of the 2021 Cumbre Vieja ridge eruption (Canary Islands): Products and dynamics controlling edifice growth and collapse

Tajogaite cone in the Cumbre Vieja ridge (La Palma, Canary Islands) erupted between 19 September and 13 December 2021. The tephra and lava sourced from the newly formed fissure rapidly built a pyroclastic cone. During the early days of eruption and after several small-scale landslides, the west flank of the edifice partially collapsed on 25 September, breaching the cone and emplacing a prominent raft-bearing lava flow. Our research combines direct observations, digital elevation models, thermal and visible imaging, and textural and compositional investigation of the explosive products to describe and characterize the edifice growth and collapse. The cone built over a steep slope (26°) and its failure occurred after an intense phase of lava fountaining (up to 30 m3 s−1) that produced rapid pyroclastic accumulation. We suggest that an increased magma supply, to an ascent rate of 0.30 m s−1, led to the rapid growth of the cone (at 2.4 × 106 m3 day−1). Simultaneously, the SW lava flow reactivated and formed a lava ‘seep’ that undercut the flank of the cone, triggering a lateral collapse via rotational rockslide that moved at minimum speeds of 34–70 m h−1. The lateral collapse formed a ~ 200 m wide scar, involving 5.5 × 106 m3 of material, and covered 1.17 km2 with decametric edifice portions and raft-bearing lava. The collapse produced a modest change in the vent geometry, but did not affect eruptive activity long term. A short pause in the eruption after the collapse may have been favored by rapid emptying of the shallower magma system, reducing ascent rates and increasing crystallization times. These results reveal the complex chain of events related to the growth and destruction of newly formed volcanic cones and highlight hazards when situated close to inhabited areas

Effect of glass on the frictional behavior of basalts at seismic slip rates

We performed 31 friction experiments on glassy basalts (GB) and glass-free basalts (GFB) at slip rates up to 6.5 m s−1 and normal stress up to 40 MPa (seismic conditions). Frictional weakening was associated to bulk frictional melting and lubrication. The weakening distance (Dw) was about 3 times shorter in GB than in GFB, but the steady state friction was systematically higher in GB than in GFB. The shorter Dw in GB may be explained by the thermal softening occurring at the glass transition temperature (Tg ~500°C), which is lower than the bulk melting temperature (Tm ~1250°C) of GFB. Postexperiment microanalyses suggest that the larger crystal fraction measured in GB melts results in the higher steady state friction value compared to the GFB melts. The effect of interstitial glass is to facilitate frictional instability and rupture propagation in GB with respect to GFB

Quantifying Microstructural Evolution in Moving Magma

Many of the grand challenges in volcanic and magmatic research are focused on understanding the dynamics of highly heterogeneous systems and the critical conditions that enable magmas to move or eruptions to initiate. From the formation and development of magma reservoirs, through propagation and arrest of magma, to the conditions in the conduit, gas escape, eruption dynamics, and beyond into the environmental impacts of that eruption, we are trying to define how processes occur, their rates and timings, and their causes and consequences. However, we are usually unable to observe the processes directly. Here we give a short synopsis of the new capabilities and highlight the potential insights that in situ observation can provide. We present the XRheo and Pele furnace experimental apparatus and analytical toolkit for the in situ X-ray tomography-based quantification of magmatic microstructural evolution during rheological testing. We present the first 3D data showing the evolving textural heterogeneity within a shearing magma, highlighting the dynamic changes to microstructure that occur from the initiation of shear, and the variability of the microstructural response to that shear as deformation progresses. The particular shear experiments highlighted here focus on the effect of shear on bubble coalescence with a view to shedding light on both magma transport and fragmentation processes. The XRheo system is intended to help us understand the microstructural controls on the complex and non-Newtonian evolution of magma rheology, and is therefore used to elucidate the many mobilization, transport, and eruption phenomena controlled by the rheological evolution of a multi-phase magmatic flows. The detailed, in situ characterization of sample textures presented here therefore represents the opening of a new field for the accurate parameterization of dynamic microstructural control on rheological behavior