Skip to content

Vesicle evolution during rhyolitic Plinian eruptions

Jong Gil Park 1, Helge M. Gonnermann1


Abstract

Plinian eruptions of rhyolitic magma are characterized by magma fragmentation, producing pyroclasts with bubble number densities (BNDs) of 1015±1 m-3 at discharge rates of 106 - 109 kg/s. What magma ascent pathways, that is conduit geometries, give rise to Plinian eruptions, and what are the resulting magma decompression rates? Our study addresses this by modeling bubble nucleation and growth under a large number of randomly chosen decompression paths between the chamber and surface, generated using the Monte Carlo method. We use a conduit model to interpret the simulated results, which allows us to derive conduit size and discharge rate and account for permeable outgassing after fragmentation. For each given decompression scenario, the model predicts whether fragmentation occurs and at what depth. Bubble number density, vesicularity, and residual water content of pyroclasts, as well as the magma discharge rate and associated conduit geometry, are all model predictions, as opposed to prior chosen parameters. We find that multiple bubble nucleation events, ranging from three to five, occur during rhyolitic Plinian eruptions, producing a wide range of size distributions. These degassing processes are explained by volcanic conduits, which narrow as they approach the surface. A bubble number density of 1015±1 m-3 results from a decompression rate of 1-2 order of magnitude in MPa/s at fragmentation. Additionally, the vesicularity and water content are 0.7-0.9 and 1.0 wt.%, respectively, aligning with observational data.