Search

Mathematical models of natural processes can be used as inversion tools to predict unobserved properties from measured quantities. Uncertainty in observations and model formulation impact on the efficacy of inverse modelling. We present a general methodology, history matching, that can be used to investigate the effect of observational and model uncertainty on inverse modelling studies. We demonstrate history matching on an integral model of volcanic plumes that is used to estimate source conditions from observations of the rise height of plumes during the eruptions of Eyjafjallajökull, Iceland, in 2010 and Grímsvötn, Iceland, in 2011. Sources of uncertainty are identified and quantified, and propagated through the integral plume model. A preliminary sensitivity analysis is performed to identify the uncertain model parameters that strongly influence model predictions. Model predictions are assessed against observations through an implausibility measure that rules out model inputs that are considered implausible given the quantified uncertainty. We demonstrate that the source mass flux at the volcano can be estimated from plume height observations, but the magmatic temperature, exit velocity and exsolved gas mass fraction cannot be accurately determined. Uncertainty in plume height observations and entrainment coefficients results in a large range of plausible values of the source mass flux. Our analysis shows that better constraints on entrainment coefficients for volcanic plumes and more precise observations of plume height are required to obtain tightly constrained estimates of the source mass flux.

Observations of volcanic lightning made using a lightning mapping array during the 2010 eruption of Eyjafjallajökull allow the trajectory and growth of the volcanic plume to be determined. The lightning observations are compared with predictions of an integral model of volcanic plumes that includes descriptions of the interaction with wind and the effects of moisture. We show that the trajectory predicted by the integral model closely matches the observational data and the model well describes the growth of the plume downwind of the vent. Analysis of the lightning signals reveals information on the dominant charge structure within the volcanic plume. During the Eyjafjallajökull eruption both monopole and dipole charge structures were observed in the plume. By using the integral plume model, we propose the varying charge structure is connected to the availability of condensed water and low temperatures at high altitudes in the plume, suggesting ice formation may have contributed to the generation of a dipole charge structure via thunderstorm-style ice-based charging mechanisms, though overall this charging mechanism is believed to have had only a weak influence on the production of lightning.