The fragmentation-induced fluidisation of pyroclastic density currents

Pyroclastic density currents (PDCs) are the most lethal volcanic process on Earth. Forecasting their inundation area is essential to mitigate their risk, but existing models are limited by our poor understanding of their dynamics. Here, we explore the role of evolving grain-size distribution in controlling the runout of the most common PDCs, known as block-and-ash flows (BAFs). Through a combination of theory, analysis of deposits and experiments of natural mixtures, we show that rapid changes of the grain-size distribution transported in BAFs result in the reduction of pore volume (compaction) within the first kilometres of their runout. We then use a multiphase flow model to show how the compressibility of granular mixtures leads to fragmentation-induced fluidisation (FIF) and excess pore-fluid pressure in BAFs. This process dominates the first ~2 km of their runout, where the effective friction coefficient is progressively reduced. Beyond that distance, transport is modulated by diffusion of the excess pore pressure. Fragmentation-induced fluidisation provides a physical basis to explain the decades-long use of low effective friction coefficients used in depth-averaged simulations required to match observed flow inundation.

Phase partitioning during fragmentation revealed by QEMSCAN Particle Mineralogical Analysis of volcanic ash

Volcanic ash particle properties depend upon their genetic fragmentation processes. Here, we introduce QEMSCAN Particle Mineralogical Analysis (PMA) to quantify the phase distribution in ash samples collected during activity at Santiaguito, Guatemala and assess the fragmentation mechanisms. Volcanic ash from a vulcanian explosion and from a pyroclastic density current resulting from a dome collapse were selected. The ash particles resulting from both fragmentation modes are dense and blocky, typical of open-vent dome volcanoes and have a componentry consistent with their andesitic composition. We use image analysis to compare the fraction of each phase at particle boundaries compared to the total particle fraction. Our results show that the explosion-derived ash has an even distribution of plagioclase and glass, but boundaries enriched in pyroxene and amphibole. In contrast, the ash generated during dome collapse has an increased fraction of glass and decreased fraction of plagioclase at particle boundaries, suggesting that fractures preferentially propagate through glass during abrasion and milling in pyroclastic flows. This study presents QEMSCAN PMA as a new resource to identify generation mechanisms of volcanic ash, which is pertinent to volcanology, aviation, respiratory health and environmental hazards, and highlights the need for further experimental constraints on the fragmentation mechanism fingerprint.

Attrition in the kimberlite system

The sustained transportation of particles in a suspension commonly results in particle attrition leading to grain size reduction and shape modification. Particle attrition is a well-studied phenomenon that has mainly focussed on sediments produced in aeolian or fluvial environments. Here, we present analogue experiments designed to explore processes of attrition in the kimberlite system; we focus on olivine as it is the most abundant constituent of kimberlite. The attrition experiments on olivine use separate experimental set-ups to approximate two natural environments relevant to kimberlites. Tumbling mill experiments feature a low energy system supporting near continual particle-particle contact and are relevant to re-sedimentation and dispersal processes. Experiments performed in a fluidized particle bed constitute a substantially higher energy environment pertinent to kimberlite ascent and eruption. The run-products of each experiment are analysed for grain size reduction and shape modification and these data are used to elucidate the rates and extents of olivine attrition as a function of time and energy. Lastly, we model the two experimental datasets with an empirical rate equation that describes the production of daughter products (fines) with time. Both datasets approach a fines production limit, or plateau, at long particle residence times; the fluidized system is much more efficient producing a substantially higher fines content and reaches the plateau faster. Our experimental results and models provide a way to forensically examine a wide range of processes relevant to kimberlite on the basis of olivine size and shape properties.

Ash production by attrition in volcanic conduits and plumes

Tephra deposits result from explosive volcanic eruption and serve as indirect probes into fragmentation processes operating in subsurface volcanic conduits. Primary magmatic fragmentation creates a population of pyroclasts through volatile-driven decompression during conduit ascent. In this study, we explore the role that secondary fragmentation, specifically attrition, has in transforming primary pyroclasts upon transport in volcanic conduits and plumes. We utilize total grain size distributions from a suite of natural and experimentally produced tephra to show that attrition is likely to occur in all explosive volcanic eruptions. Our experimental results indicate that fine ash production and surface area generation is fast (<15 min) thereby rapidly raising the fractal dimension of tephra deposits. Furthermore, a new metric, the Entropy of Information, is introduced to quantify the degree of attrition (secondary fragmentation) from grain size data. Attrition elevates fine ash production which, in turn, has consequences for eruption column stability, tephra dispersal, aggregation, volcanic lightening generation, and has concomitant effects on aviation safety and Earth's climate.

On the fate of pumice rafts formed during the 2012 Havre submarine eruption

Pumice rafts are floating mobile accumulations of low-density pumice clasts generated by silicic volcanic eruptions. Pumice in rafts can drift for years, become waterlogged and sink, or become stranded on shorelines. Here we show that the pumice raft formed by the impressive, deep submarine eruption of the Havre caldera volcano (Southwest Pacific) in July 2012 can be mapped by satellite imagery augmented by sailing crew observations. Far from coastal interference, the eruption produced a single >400 km² raft in 1 day, thus initiating a gigantic, high-precision, natural experiment relevant to both modern and prehistoric oceanic surface dispersal dynamics. Observed raft dispersal can be accurately reproduced by simulating drift and dispersal patterns using currents from an eddy-resolving ocean model hindcast. For future eruptions that produce potentially hazardous pumice rafts, our technique allows real-time forecasts of dispersal routes, in addition to inference of ash/pumice deposit distribution in the deep ocean.

Pumice rafts result from volcanic eruptions into and onto water, and can be extensive and potentially hazardous, but tracking their dispersal is difficult. Jutzeler et al. combine satellite imagery and an ocean model to accurately forecast pumice raft dispersal routes.

Sedimentology and geomorphology of the deposits from the August 2006 pyroclastic density currents at Tungurahua volcano, Ecuador

The deposits of the pyroclastic density currents from the August 2006 eruption of Tungurahua show three facies associations depending on the topographic setting: the massive, proximal cross-stratified, and distal cross-stratified facies. (1) The massive facies is confined to valleys on the slopes of the volcano. It contains clasts of >1 m diameter to fine ash material, is massive, and interpreted as deposited from dense pyroclastic flows. Its surface can exhibit lobes and levees covered with disk-shaped and vesicular large clasts. These fragile large clasts must have rafted at the surface of the flows all along the path in order to be preserved, and thus imply a sharp density boundary near the surface of these flows. (2) The proximal cross-stratified facies is exposed on valley overbanks on the upper part of the volcano and contains both massive coarse-grained layers and cross-stratified ash and lapilli bedsets. It is interpreted as deposited from (a) dense pyroclastic flows that overflowed the gentle ridges of valleys of the upper part of the volcano and (b) dilute pyroclastic density currents created from the dense flows by the entrainment of air on the steep upper flanks. (3) The distal cross-stratified facies outcrops as spatially limited, isolated, and wedge-shaped bodies of cross-stratified ash deposits located downstream of cliffs on valleys overbanks. It contains numerous aggrading dune bedforms, whose crest orientations reveal parental flow directions. A downstream decrease in the size of the dune bedforms, together with a downstream fining trend in the grain size distribution are observed on a 100-m scale. This facies is interpreted to have been deposited from dilute pyroclastic density currents with basal tractional boundary layers. We suggest that the parental flows were produced from the dense flows by entrainment of air at cliffs, and that these diluted currents might rapidly deposit through "pneumatic jumps". Three modes are present in the grain size distribution of all samples independently of the facies, which further supports the interpretation that all three facies derive from the same initial flows. This study emphasizes the influence of topography on small volume pyroclastic density currents, and the importance of flow transformation and flow-stripping processes.