The Evolution of Paleo-Porosity in Basalts: Reversing Pore-Filling Mechanisms Using X-Ray Computed Tomography

Often carrying a high-volume fraction of vesicles, basaltic rocks can be an important reservoir horizon in petroleum systems, and are considered an excellent candidate for CO2 storage by in situ mineral trapping. The frequency of amygdaloidal basalts in many sequences highlights the prevalence of mineralisation, but when the vesicle network has been filled, the basalts can act as impermeable seals and traps. Characterising the spatial and temporal evolution of the porosity and permeability is critical to understanding the petro-physical properties and CO2 storage potential of basalts. We exploit X-ray computed tomography (XCT) to investigate the precipitation history of an amygdaloidal basalt containing a pore-connecting micro fracture network now partially filled by calcite as an analogue for CO2 mineral trapping in a vesicular basalt. The fracture network likely represents a preferential pathway for CO2-rich fluids during mineralisation. We investigate and quantify the evolution of basalt porosity and permeability during pore-filling calcite precipitation by applying novel numerical erosion techniques to “back-strip” the calcite from the amygdales and fracture networks. We provide a semi-quantitative technique for defining reservoir potential and quality through time and understanding sub-surface flow and storage. We found that permeability evolution is dependent on the precipitation mechanism and rates, as well as on the presence of micro fracture networks, and that once the precipitation is sufficient to close off all pores, permeability reaches values that are controlled by the micro fracture network. These results prompt further studies to determine CO2 mineral trapping mechanisms in amygdaloidal basalts as analogues for CO2 injections in basalt formations

Sintering of vesiculating pyroclasts

Hot volcanic pyroclasts can sinter, vesiculate, and outgas in concert – a combination of processes which remains poorly constrained. And yet this combination of processes can occur coincidently during deposition from pyroclastic density currents, in conduit-filling pyroclastic debris, and in tuffisites. In many of these settings, it is the sintering-driven evolution of permeability that is key to gas transport through the evolving deposit. Here, we experimentally and theoretically investigate the evolution of the permeable networks during sintering of hot fragmental volcanic systems, which are hydrous and oversaturated at the experimental conditions. Firstly, we find that vesiculation results in shutting of the inter-granular porous network as bubble growth drives expansion of the particles into one another, destroying interconnected pores. Secondly, we observe that degassing by diffusion out of the particle edge results in contraction of the vesicular particles, re-opening pore spaces between them. Therefore, we find that vesiculation, and diffusive outgassing compete to determine both the intra-fragment vesicularity and the permeability during sintering. The development of intra-fragment vesicularity directly impacts the inter-fragment pore space and its connectivity, which decreases during vesiculation and subsequently increases during diffusive outgassing, prompting complex, non-linear permeability evolution.The relative dominance of these processes is fragment size dependent; proportionally, fine fragments lose gas at a higher rate than coarser fragments during diffusive outgassing due to larger surface area to volume ratios. As the systems progress, larger fragments retain a higher proportion of gas and so attain greater vesicularities than finer ones – and therefore, the coarse fragmental pyroclasts experience a greater, yet transient, reduction in connected porosity and permeability. We suggest that where vesiculation is sufficient, it can lead to the complete loss of connected porosity and the sealing of permeable pathways much earlier than in a sintering-only system. Our results suggest that classical sintering models must be modified to account for these vesiculation and diffusive degassing processes, and that only a combined vesiculation, sintering, and diffusive outgassing model can resolve the evolution of permeability in hot clastic volcanic systems

Obsidian clasts as sintered remnants of agglutination processes in volcanic conduits, evidence from the Pepom tephras (Sete Cidades), São Miguel, Azores

The youngest explosive eruptions of the Sete Cidades volcano, São Miguel, Azores, are recorded by a series of relatively small-volume (<0.03 to 0.13 km3 DRE) trachytic pyroclastic deposits termed the Pepom tephra deposits. While dominated by crystal-poor to crystal-moderate (<10%) pumice clasts, these deposits also contain a suite of dense glassy clasts of broadly similar crystallinity. The obsidian clasts from a single deposit vary in texture from entirely dense to those that are moderately vesicular and typically single clasts will contain multiple textural domains. The majority (∼71%) of these dense clasts have compositions both from bulk rock and in-situ glass measurements that are identical to those of the pumice clasts within the same deposit. We interpret these dense clasts to reflect sintering of previously fragmented magma at shallow levels in the conduit prior to being re-entrained and erupted with the vesicular magma, in agreement with recent studies focussing on textural observations. Notably, across the exposed volcanic stratigraphy of São Miguel obsidian domes, flows/coulees are not preserved, arguing against the idea that the dense glass clasts within the Pepom tephras are sourced from existing surficial rocks. In contrast, the neighbouring island of Terceira exhibits domes and coulees with large obsidian bands that cut through the crystalline groundmass. Most silicic rocks of Santa Bárbara and Pico Alto volcanoes on Terceira are peralkaline, comenditic to pantelleritic in composition, and at similar conditions (e.g., temperature and water content) have lower viscosities than the trachytic Pepom obsidian clasts. However, the Santa Bárbara silicic lavas on Terceira (the less peralkaline suite) have more obsidian than the more peralkaline Pico Alto domes and coulees indicating that while peralkalinity, developed during magmatic evolution in the crust, may play a role, sintering occurring at shallow levels within the conduit likely is more important in producing obsidian

Permeability of compacting porous lavas

The highly transient nature of outgassing commonly observed at volcanoes is in part controlled by the permeability of lava domes and shallow conduits. Lava domes generally consist of a porous outer carapace surrounding a denser lava core with internal shear zones of variable porosity. Here we examine densification using uniaxial compression experiments on variably crystalline and porous rhyolitic dome lavas from the Taupo Volcanic Zone. Experiments were conducted at 900°C and an applied stress of 3MPa to 60% strain, while monitoring acoustic emissions to track cracking. The evolution of the porous network was assessed via X-ray computed tomography, He-pycnometry, and relative gas permeability. High starting connected porosities led to low apparent viscosities and high strain rates, initially accompanied by abundant acoustic emissions. As compaction ensued, the lavas evolved; apparent viscosity increased and strain rate decreased due to strain hardening of the suspensions. Permeability fluctuations resulted from the interplay between viscous flow and brittle failure. Where phenocrysts were abundant, cracks had limited spatial extent, and pore closure decreased axial and radial permeability proportionally, maintaining the initial anisotropy. In crystal-poor lavas, axial cracks had a more profound effect, and permeability anisotropy switched to favor axial flow. Irrespective of porosity, both crystalline samples compacted to a threshold minimum porosity of 17–19%, whereas the crystal-poor sample did not achieve its compaction limit. This indicates that unconfined loading of porous dome lavas does not necessarily form an impermeable plug and may be hindered, in part by the presence of crystals

Vesiculation and Quenching During Surtseyan Eruptions at Hunga Tonga-Hunga Ha'apai Volcano, Tonga

Surtseyan eruptions are shallow to emergent subaqueous explosive eruptions that owe much of their characteristic behavior to the interaction of magma with water. The difference in thermal properties between water and air affects the cooling and postfragmentation vesiculation processes in magma erupted into the water column. Here we study the vesiculation and cooling processes during the 2009 and 2014–2015 Surtseyan eruptions of Hunga Tonga‐Hunga Ha'apai volcano by combining 2‐D and 3‐D vesicle‐scale analyses of lapilli and bombs and numerical thermal modeling. Most of the lapilli and bombs show gradual textural variations from rim to core. The vesicle connectivity in the lapilli and bombs increases with vesicularity from fully isolated to completely connected and also increases from rim to core in transitional clasts. We interpret the gradual textural variations and the connectivity‐vesicularity relationships as the result of postfragmentation bubble growth and coalescence interrupted at different stages by quenching in water. The measured vesicle size distributions are bimodal with a population of small and large vesicles. We interpret this bimodality as the result of two nucleation events, one prefragmentation with the nucleation and growth of large bubbles and one postfragmentation with nucleation of small vesicles. We link the thermal model with the textural variations in the clasts—showing a dependence on particle size, Leidenfrost effect, and initial melt temperature. In particular, the cooling profiles in the bombs are consistent with the gradual textural variations from rim to core in the clasts, likely caused by variations in time available for vesiculation before quenching

Silicic conduits as supersized tuffisites:Clastogenic influences on shifting eruption styles at Cordón Caulle volcano (Chile)

Understanding the processes that drive explosive-effusive transitions during large silicic eruptions is crucial to hazard mitigation. Conduit models usually treat magma ascent and degassing as a gradual, unidirectional progression from bubble nucleation through magmatic fragmentation. However, there is growing evidence for the importance of bi-directional clastogenic processes that sinter fragmented materials into coherent clastogenic magmas. Bombs that were ejected immediately before the first emergence of lava in the 2011–2012 eruption at Cordón Caulle volcano (Chile) are texturally heterogeneous composite assemblages of welded pyroclastic material. Although diverse in density and appearance, SEM and X-ray tomographic analysis show them all to have been formed by multi-generational viscous sintering of fine ash. Sintering created discrete clasts ranging from obsidian to pumice and formed a pervasive clast-supporting matrix that assembled these clasts into a conduit-sealing plug. An evaluation of sintering timescales reveals texturally disparate bomb components to represent only minutes of difference in residence time within the conduit. Permeability modelling indicates that the plug was an effective conduit seal, with outgassing potential—even from high-porosity regions—being limited by the inability of gas to flow across tendrils of densely sintered inter-clast matrix. Contrary to traditional perspectives, declining expressions of explosivity at the surface need not be preceded or accompanied by a decline in fragmentation efficiency. Instead, they result from tips in balance between the opposing processes of fragmentation and sintering that occur in countless cycles within volcanic conduits. These processes may be particularly enhanced at silicic fissure volcanoes, which have laterally extensive subsurface plumbing systems that require complex magma ascent pathways. The textures investigated here reveal the processes occurring within silicic fissures to be phenomenologically identical to those that have been inferred to occur in tuffisite veins: silicic conduits are essentially supersized examples of edifice-penetrating tuffisites

Fracture healing and strength recovery in magmatic liquids

Cycles of fracture and healing in magma are important controls on outgassing time scales and repetitive seismicity at silicic volcanoes. Here, we experimentally drove silicate melts (at 109–1011 Pa·s) to tensile failure, measuring the strength during fracture of the otherwise liquid material. We then took the same melts with parallel contact surfaces and closed the fracture under compressive stress and recorded the evolution of tensile strength of the interface healed for different times. We provide a semi-empirical model for fracture healing time scales useful for volcanic applications. As the time available for healing is increased, strength nonlinearly recovers toward that of unfractured glass. We parameterized the healing kinetics as a three-stage process: (1) relaxation of the compressive stress, (2) fracture surface–surface wetting, and (3) diffusive removal of the interface. During welding of these surfaces in air, we observed that micropores are trapped along the fracture plane, which may inhibit complete healing and provide a textural record of relict fracture planes in volcanic glasses. We conclude that at magmatic conditions, fracture healing is efficient in crystal-poor melts, and it could rapidly seal outgassing pathways over eruptive time scales, contributing to cyclic behavior associated with recurring gas-and-ash explosions and outgassing events

Permeability Evolution in Variably Glassy Basaltic Andesites Measured Under Magmatic Conditions

Heat from inflowing magma may act to seal permeable networks that assist passive outgassing at volcanic conduit margins and in overlying domes, reducing the efficiency of overpressure dissipation. Here we present a study of the evolution of permeability—measured under magmatic conditions—with increasing temperature in glassy and glass‐poor basaltic andesites from Merapi volcano (Indonesia). Whereas the permeability of glass‐poor rocks decreases little up to a temperature of 1,010°C, glassy specimens experience a pronounced decrease in permeability above the glass transition once the viscosity of the crystal suspension is low enough to relax under external stresses. Changes in temperature alone are thus not enough to significantly modify the permeability of the glass‐poor rocks that commonly form Merapi's dome. However, the presence of glass‐rich domains in a dome may lead to local sealing of the volcanic plumbing between eruptions, exacerbating localized overpressure development that could contribute to explosivity

Explosive-effusive volcanic eruption transitions caused by sintering

Silicic volcanic activity has long been framed as either violently explosive or gently effusive. However, recent observations demonstrate that explosive and effusive behavior can occur simultaneously. Here, we propose that rhyolitic magma feeding subaerial eruptions generally fragments during ascent through the upper crust and that effusive eruptions result from conduit blockage and sintering of the pyroclastic products of deeper cryptic fragmentation. Our proposal is supported by (i) rhyolitic lavas are volatile depleted; (ii) textural evidence supports a pyroclastic origin for effusive products; (iii) numerical models show that small ash particles ≲10-5 m can diffusively degas, stick, and sinter to low porosity, in the time available between fragmentation and the surface; and (iv) inferred ascent rates from both explosive and apparently effusive eruptions can overlap. Our model reconciles previously paradoxical observations and offers a new framework in which to evaluate physical, numerical, and geochemical models of Earth's most violent volcanic eruptions

Pressure-Driven Opening and Filling of a Volcanic Hydrofracture Recorded by Tuffisite at Húsafell, Iceland:A Potential Seismic Source

The opening of magmatic hydraulic fractures is an integral part of magma ascent, the triggering of volcano seismicity, and defusing the explosivity of ongoing eruptions via outgassing magmatic volatiles. If filled with pyroclastic particles, these fractures can be recorded as tuffisites. Tuffisites are therefore thought to play a key role in both initiating eruptions and controlling their dynamics, and yet their genesis remains poorly understood. Here we characterise the processes, pressures and timescales involved in tuffisite evolution within the country rock through analysis of the sedimentary facies and structures of a large sub-horizontal tuffisite vein, 0.9 m thick and minimum 40 m in length, at the dissected Húsafell volcano, western Iceland. The vein occurs where a propagating rhyolitic sheet intrusion stalled at a depth of ∼500 m beneath a relatively strong layer of welded ignimbrite. Laminations, cross-stratification, channels, and internal injections indicate erosion and deposition in multiple fluid pulses, controlled by fluctuations in local fluid pressure and changes in fluid-particle concentration. The field evidence suggests that this tuffisite was emplaced by as many as twenty pulses, depositing sedimentary units with varying characteristics. Assuming that each sedimentary unit (∼0.1 m thick and minimum 40 m in length) is emplaced by a single fluid pulse, we estimate fluid overpressures of ∼1.9–3.3 MPa would be required to emplace each unit. The Húsafell tuffisite records the repeated injection of an ash-laden fluid within an extensive subhorizontal fracture, and may therefore represent the fossil record of a low-frequency seismic swarm associated with fracture propagation and reactivation. The particles within the tuffisite cool and compact through time, causing the rheology of the tuffisite fill to evolve and influencing the nature of the structures being formed as new material is injected during subsequent fluid pulses. As this new material is emplaced, the deformation style of the surrounding tuffisite is strongly dependent on its evolving rheology, which will also control the evolution of pressure and the system permeability. Interpreting tuffisites as the fossil record of fluid-driven hydrofracture opening and evolution can place new constraints on the cycles of pressurisation and outgassing that accompany the opening of magmatic pathways, key to improving interpretations of volcanic unrest and hazard forecasting