Sensitivity analysis for volcanic source modeling quality assessment and model selection

The increasing knowledge and understanding of volcanic sources has led to the development and implementation of sophisticated and complex mathematical models with the main goal of describing field and experimental data. Quantification of the model's ability in describing the data becomes fundamental for a realistic estimate of the model parameters. The analysis of sensitivity can help us in identifying the parameters that significantly affect the model's output and in assessing its quality factor. In this paper, we describe the Global Sensitivity Analysis (GSA) methods based both on Fourier Amplitude Sensitivity Test and on the Sobol' approach and discuss their implementation in a Matlab software tool (GSAT). We also introduce a new criterion for model selection based on sensitivity analysis. The proposed approach is tested and applied to quantify the fitting ability of an analytic volcanic source model on a synthetic deformation data. Results show the validity of the method, against the traditional approaches, in supporting the volcanic model selection and the flexibility of the GSAT software tool in analyzing the model sensitivity.

[1]  E. Triantaphyllou,et al.  A Sensitivity Analysis Approach for Some Deterministic Multi-Criteria Decision-Making Methods* , 1997 .

[2]  I. M. Soboĺ.  Quasi-Monte Carlo methods , 1990 .

[3]  Kristy F. Tiampo,et al.  Volcanic source inversion using a genetic algorithm and an elastic-gravitational layered earth model for magmatic intrusions , 2004, Comput. Geosci..

[4]  James H. Dieterich,et al.  Finite element modeling of surface deformation associated with volcanism , 1975 .

[5]  Saltelli Andrea,et al.  Global Sensitivity Analysis: The Primer , 2008 .

[6]  Emanuele Borgonovo,et al.  A new uncertainty importance measure , 2007, Reliab. Eng. Syst. Saf..

[7]  I. Sobol,et al.  Global sensitivity indices for nonlinear mathematical models. Review , 2005 .

[8]  I. Sobola,et al.  Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates , 2001 .

[9]  C. Fortuin,et al.  Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients. I Theory , 1973 .

[10]  K. Shuler,et al.  Study of the sensitivity of coupled reaction systems to uncertainties in rate coefficients. III. Analysis of the approximations , 1975 .

[11]  A. Saltelli,et al.  A quantitative model-independent method for global sensitivity analysis of model output , 1999 .

[12]  G. Glatting,et al.  Choosing the optimal fit function: comparison of the Akaike information criterion and the F-test. , 2007, Medical physics.

[13]  B. Matthews.  Comparison of the predicted and observed secondary structure of T4 phage lysozyme. , 1975, Biochimica et biophysica acta.

[14]  Ilya M. Sobol,et al.  Theorems and examples on high dimensional model representation , 2003, Reliab. Eng. Syst. Saf..

[15]  Stefano Tarantola,et al.  Variance‐Based Methods , 2008 .

[16]  I. Sobol.  On the distribution of points in a cube and the approximate evaluation of integrals , 1967 .

[17]  Jacob Cohen.  A Coefficient of Agreement for Nominal Scales , 1960 .

[18]  Bryan E. Shepherd,et al.  Global Sensitivity Analysis. The Primer by SALTELLI, A., RATTO, M., ANDRES, T., CAMPOLONGO, F., CARIBONI, J., GATELLI, D., SAISANA, M., and TARANTOLA, S , 2009 .

[19]  Paul M. Davis,et al.  Surface deformation due to inflation of an arbitrarily oriented triaxial ellipsoidal cavity in an elastic half-space, with reference to Kilauea volcano, Hawaii , 1986 .

[20]  Vi︠a︡cheslav Moiseevich Zobin.  Introduction to Volcanic Seismology , 2003 .

[21]  D. Parkinson,et al.  Bayesian Methods in Cosmology: Model selection and multi-model inference , 2009 .

[22]  Andrea Saltelli,et al.  Use of Sobol's Quasirandom Sequence Generator for Integration of Modified Uncertainty Importance Measure , 1995 .

[23]  Lauwerens Kuipers,et al.  Uniform distribution of sequences , 1974 .

[24]  Sergei S. Kucherenko,et al.  Derivative based global sensitivity measures and their link with global sensitivity indices , 2009, Math. Comput. Simul..

[25]  A. Saltelli,et al.  Importance measures in global sensitivity analysis of nonlinear models , 1996 .

[26]  J. B. Walsh,et al.  Surface deformation associated with volcanism , 1971 .

[27]  Robert M. Owen,et al.  Statistical Tests for Closure of Plate Motion Circuits (Paper 7L7161) , 1987 .

[28]  Harvey M. Wagner,et al.  Global Sensitivity Analysis , 1995, Oper. Res..

[29]  Thomas M. Ludden,et al.  Comparison of the Akaike Information Criterion, the Schwarz criterion and the F test as guides to model selection , 1994, Journal of Pharmacokinetics and Biopharmaceutics.

[30]  David J. Pannell,et al.  Sensitivity Analysis of Normative Economic Models: Theoretical Framework and Practical Strategies , 1997 .

[31]  Ron Kohavi,et al.  Guest Editors' Introduction: On Applied Research in Machine Learning , 1998, Machine Learning.

[32]  Moon-Hyun Chun,et al.  An uncertainty importance measure using a distance metric for the change in a cumulative distribution function , 2000, Reliab. Eng. Syst. Saf..

[33]  Salvatore Gambino,et al.  A new dyke intrusion style for the Mount Etna May 2008 eruption modelled through continuous tilt and GPS data , 2009 .

[34]  David P. Hill,et al.  Analytical modeling of gravity changes and crustal deformation at volcanoes: The Long Valley caldera, California, case study , 2009 .

[35]  Ilya M. Sobol,et al.  Sensitivity Estimates for Nonlinear Mathematical Models , 1993 .

[36]  K. Mogi.  Relations between the Eruptions of Various Volcanoes and the Deformations of the Ground Surfaces around them , 1958 .

[37]  Mimmo Palano,et al.  Imaging the multi‐level magma reservoir at Mt. Etna volcano (Italy) , 2011 .

[38]  Andrea Saltelli,et al.  Sensitivity Analysis for Importance Assessment , 2002, Risk analysis : an official publication of the Society for Risk Analysis.

[39]  M. Kramer,et al.  Sensitivity Analysis in Chemical Kinetics , 1983 .