We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

An experimentally validated numerical model for bubble growth in magma.

Coumans, J.P. and Llewellin, E.W. and Wadsworth, F.B. and Humphreys, M.C.S. and Mathias, S.A. and Yelverton, B.M. and Gardner, J.E. (2020) 'An experimentally validated numerical model for bubble growth in magma.', Journal of volcanology and geothermal research., 402 . p. 107002.


Volcanic eruptions are driven by the growth of gas bubbles in magma. The timing and rate of bubble growth are important because they determine whether enough gas pressure can develop to fragment the melt. Bubbles grow in response to decompression and diffusive transport of dissolved volatiles (predominantly H2O) that exsolve into the bubbles. Growth is resisted by the viscosity of the melt. Both melt viscosity and H2O diffusivity have non-linear dependence on the concentration of H2O dissolved in the melt, which necessitates a numerical approach to modelling bubble growth. Several bubble growth models have previously been published and applied, but none of them has been validated against continuous, in situ experimental data or provided as a user-friendly tool. Here we present a numerical bubble growth model, implemented in MATLAB, which allows for arbitrary temperature and pressure pathways, and accounts for the impact of spatial variations in dissolved H2O concentration on viscosity and diffusivity. We validate the model against two sets of experimental data: (1) New continuous data for gas-volume fraction as a function of time, collected using optical dilatometry of vesiculating hydrous obsidian samples which were heated from 930 °C to 1000 °C at atmospheric pressure. This dataset captures isobaric, isothermal bubble growth under strongly disequilibrium conditions. (2) Discrete data from published decompression experiments at 825 °C and pressures from 200 MPa to ~5 MPa with decompression rates from 0.1 MPa s−1 to 10 MPa s−1. These experiments represent isothermal, decompression-driven bubble growth spanning equilibrium to strongly disequilibrium conditions. The numerical model closely reproduces the experimental data across all conditions, providing validation against contrasting bubble growth scenarios. The validated model has a wide range of potential volcanological applications, including forward modelling of bubble growth phenomena, and inverse modelling to reconstruct pressure–temperature–time pathways from textures and volatile contents of eruptive products.

Item Type:Article
Full text:(VoR) Version of Record
Available under License - Creative Commons Attribution.
Download PDF
Publisher Web site:
Publisher statement:© 2020 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license (
Date accepted:12 July 2020
Date deposited:04 August 2020
Date of first online publication:19 July 2020
Date first made open access:04 August 2020

Save or Share this output

Look up in GoogleScholar