Search

Factorisation is widely used in the analysis of numerical simulations. It allows changes in properties of a system to be attributed to changes in multiple variables associated with that system. There are many possible factorisation methods; here we discuss three previously-proposed factorisations that have been applied in the field of climate modelling: the linear factorisation, the Stein and Alpert (1993) factorisation, and the Lunt et al (2012) factorisation. We show that, when more than two variables are being considered, none of these three methods possess all three properties of uniqueness, symmetry, and completeness. Here, we extend each of these factorisations so that they do possess these properties for any number of variables, resulting in three factorisations – the linear-sum factorisation, the shared-interaction factorisation, and the scaled-total factorisation. We show that the linear-sum factorisation and the shared-interaction factorisation reduce to be identical. We present the results of the factorisations in the context of studies that used the previously-proposed factorisations. This reveals that only the linear-sum/shared-interaction factorisation possesses a fourth property – boundedness, and as such we recommend the use of this factorisation in applications for which these properties are desirable.

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.