THE EVOLUTION OF BUBBLE SIZE DISTRIBUTIONS IN VOLCANIC ERUPTIONS

Abstract

We review observations of bubble size distributions (BSDs) generated during explosive volcanic eruptions and laboratory explosions, as inferred from vesicle size distributions found in the end products. Unimodal, polymodal, exponential and power law BSDs are common, even in the absence of coalescence, and both power law and exponential distributions have been generated in the same eruption. To date theoretical models have proposed incompatible mechanisms for producing the various distributions. We here present a unifying mechanism. Data from our laboratory analogue experiments suggest that power law distributions are associated with highly non-equilibrium degassing. A numerical model is developed in which bubbles nucleate repeatedly and grow in the spaces between those of previous generations, where, in a non-equilibrium degassing scenario, the volatile concentration remains high. This process causes the BSD to evolve from unimodal, through exponential, into a power law. The exponent of the power law is a measure of the number of nucleation events, or the duration of the nucleation period compared with the timescale of bubble growth. The mathematical inevitability of the evolution from unimodal (Poissonian) to power law is discussed. The findings may resolve the apparent contradiction between the equilibrium degassing conduit flow models and the non-equilibrium degassing conditions derived from bubble growth models of explosive volcanic eruptions. The process of ongoing nucleation is the mechanism whereby the volcanic system maintains near-equilibrium in the case of rapid depressurisation and slow volatile diffusion.