The effect of nonlinear decompression history on H2O/CO2 vesiculation in rhyolitic magmas

Abstract

Magma ascent rate is one of the key parameters that control volcanic eruption style, tephra dispersion, and volcanic atmospheric impact. Many methods have been employed to investigate the magma ascent rate in volcanic eruptions, and most rely on equilibrium thermodynamics. Combining the mixed H2O-CO2 solubility model with the diffusivities of both H2O and CO2 for normal rhyolitic melt, we model the kinetics of H2O and CO2 in rhyolitic eruptions that involve nonlinear decompression rates. Our study focuses on the effects of the total magma ascent time, the nonlinearity of decompression paths, and the influence of different initial CO2/H2O content on the posteruptive H2O and CO2 concentration profiles around bubbles within the melt. Our results show that, under most circumstances, volatile diffusion profiles do not constrain a unique solution for the decompression rate of magmas during an eruption, but, instead, provide a family of decompression paths with a well-defined trade-off between ascent time and nonlinearity. An important consequence of our analysis is that the common assumption of a constant decompression rate (averaged value) tends to underestimate the actual magma ascent time.