Multiobjective genetic algorithm inversion of ground deformation and gravity changes spanning the 1981 eruption of Etna volcano

[1] During the last few decades, joint investigations of microgravity and surface deformation measurements have played an increasingly important role in studying the internal dynamics of active volcanoes. Deformation and microgravity observations have been accomplished at Mt Etna since the eighties. Past data sets collected during important paroxysmal events can be utilized as case-studies to both (1) test the possibilities of nowadays more powerful inversion tools and improved analytical formulations to model the source-mechanisms of volcano-related deformation and gravity changes and (2) in turn obtain new insights into the functioning of the plumbing system of the volcano. Here we analyze a data set spanning the March 1981 eruption of Mt. Etna. Large horizontal displacements were evidenced on the NE and SW flanks of the volcano through electrooptical distance measurements (EDM) during two 20-month periods, both encompassing the March 1981 eruption. Elevation changes, evidenced through leveling measurements, during a 12-month period spanning the eruption, were in general smaller than horizontal displacements with important amplitudes only close to the eruptive fissure. Gravity measurements, carried out together with leveling measurements, evidenced positive changes, spatially well correlated with elevation changes, but having a larger wavelength. The joint inversion of the multimethod geophysical data is regarded as a multiobjective optimization problem and solved through a Genetic Algorithm technique of the nondominated type. We conclude that a composite intrusive mechanism with two tensile cracks, each associated to a zone where preexisting microfractures were filled with new magma, leaded to the 1981 eruption. The results of the present study highlight the advantages of multiobjective evolutionary algorithms, as a powerful tool to jointly invert multimethod geophysical data, and pose important issues on the subject of volcano-monitoring.

[1]  D. Carbone,et al.  Temporal variations in gravity at Mt. Etna (Italy) associated with the 1989 and 1991 eruptions , 1998 .

[2]  C. Negro,et al.  Elastic model for the gravity and elevation changes before the 2001 eruption of Etna volcano , 2007 .

[3]  Christoph Schwarzbach,et al.  Two‐dimensional inversion of direct current resistivity data using a parallel, multi‐objective genetic algorithm , 2005 .

[4]  L. Fortuna,et al.  Integrated inversion of ground deformation and magnetic data at Etna volcano using a genetic algorithm technique , 2007 .

[5]  Carlo Giunchi,et al.  Analytical and 3-D numerical modelling of Mt. Etna (Italy) volcano inflation , 2005 .

[6]  D. E. Goldberg,et al.  Genetic Algorithms in Search , 1989 .

[7]  S. C. Scott.  Variations in lava composition during the March 1981 eruption of mount etna and the implications of a compositional comparison with earlier historic eruptions , 1983 .

[8]  Y. Okada.  Internal deformation due to shear and tensile faults in a half-space , 1992, Bulletin of the Seismological Society of America.

[9]  G. Wadge,et al.  An accurate and efficient method for including the effects of topography in three‐dimensional elastic models of ground deformation with applications to radar interferometry , 2000 .

[10]  N. Draper,et al.  Applied Regression Analysis , 1966 .

[11]  C. Negro,et al.  Modelling of ground deformation and gravity fields using finite element method: An application to Etna volcano , 2007 .

[12]  John H. Holland,et al.  Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence , 1992 .

[13]  C. Fonseca,et al.  GENETIC ALGORITHMS FOR MULTI-OBJECTIVE OPTIMIZATION: FORMULATION, DISCUSSION, AND GENERALIZATION , 1993 .

[14]  F. Beauducel,et al.  Collection and three-dimensional modeling of GPS and tilt data at Merapi volcano, Java , 1999 .

[15]  Peter J. Fleming,et al.  Genetic Algorithms for Multiobjective Optimization: FormulationDiscussion and Generalization , 1993, ICGA.

[16]  P. Segall,et al.  Rapid deformation of Kilauea Volcano: Global Positioning System measurements between 1990 and 1996 , 2000 .

[17]  P. Segall,et al.  The Interpretation of Gravity Changes and Crustal Deformation in Active Volcanic Areas , 2004 .

[18]  M. Murray,et al.  The 1992 M = 7 Cape Mendocino, California, earthquake: Coseismic deformation at the south end of the Cascadia megathrust , 1996 .

[19]  M. Bonafede,et al.  Modelling gravity variations consistent with ground deformation in the Campi Flegrei caldera (Italy) , 1998 .

[20]  B. Behncke,et al.  Lava flow hazard at Mount Etna (Italy): New data from a GIS-based study , 2005 .

[21]  John H. Holland,et al.  Adaptation in Natural and Artificial Systems: An Introductory Analysis with Applications to Biology, Control, and Artificial Intelligence , 1992 .

[22]  Paul Segall,et al.  Estimating source parameters from deformation data, with an application to the March 1997 earthquake swarm off the Izu Peninsula, Japan , 2001 .

[23]  A. Bonaccorso.  The March 1981 Mount Etna eruption inferred through ground deformation modelling , 1999 .

[24]  P. Davis,et al.  Deformation due to a rectangular tension crack in an elastic half-space , 1986 .

[25]  Shuhei Okubo,et al.  Gravity and potential changes due to shear and tensile faults in a half-space , 1992 .

[26]  Kalyanmoy Deb,et al.  A Fast Elitist Non-dominated Sorting Genetic Algorithm for Multi-objective Optimisation: NSGA-II , 2000, PPSN.

[27]  Giuseppe Nunnari,et al.  Inverse modelling of volcanomagnetic fields using a genetic algorithm technique , 2005 .

[28]  Kalyanmoy Deb,et al.  A fast and elitist multiobjective genetic algorithm: NSGA-II , 2002, IEEE Trans. Evol. Comput..

[29]  Mike Dentith,et al.  Inversion of potential field data by genetic algorithms , 1997 .

[30]  Kristy F. Tiampo,et al.  Spherical and ellipsoidal volcanic sources at Long Valley caldera, California, using a genetic algorithm inversion technique , 2000 .

[31]  G. Berrino,et al.  Ground deformation and gravity changes accompanying the March 1981 eruption of Mount Etna , 1983 .

[32]  A. Bonaccorso,et al.  Mt. Etna Volcano Laboratory , 2004 .

[33]  N. Draper,et al.  Applied Regression Analysis: Draper/Applied Regression Analysis , 1998 .

[34]  Y. Okada.  Surface deformation due to shear and tensile faults in a half-space , 1985 .

[35]  Paola Del Carlo,et al.  Types of eruptions of Etna volcano AD 1670–2003: implications for short-term eruptive behaviour , 2005 .

[36]  Kalyanmoy Deb,et al.  Searching for Robust Pareto-Optimal Solutions in Multi-objective Optimization , 2005, EMO.

[37]  M. Moorkamp,et al.  Joint inversion of teleseismic receiver functions and magnetotelluric data using a genetic algorithm: Are seismic velocities and electrical conductivities compatible? , 2007 .

[38]  K. Tiampo,et al.  New Results at Mayon, Philippines, from a Joint Inversion of Gravity and Deformation Measurements , 2004 .

[39]  H. Rymer,et al.  Gravity fields and the interpretation of volcanic structures: Geological discrimination and temporal evolution , 1986 .

[40]  Hitoshi Mori,et al.  Elastic models for the magma intrusion associated with the 2000 eruption of Usu Volcano, Hokkaido, Japan , 2003 .