PARTOS - Passive and Active Ray TOmography Software: description and preliminary analysis using TOMO-ETNA experiment’s dataset

In this manuscript we present the new friendly seismic tomography software based on joint inversion of active and passive seismic sources called PARTOS (Passive Active Ray TOmography Software). This code has been developed on the base of two well-known widely used tomographic algorithms (LOTOS and ATOM-3D), providing a robust set of algorithms. The dataset used to set and test the program has been provided by TOMO-ETNA experiment. TOMO-ETNA database is a large, high-quality dataset that includes active and passive seismic sources recorded during a period of 4 months in 2014. We performed a series of synthetic tests in order to estimate the resolution and robustness of the solutions. Real data inversion has been carried out using 3 different subsets: i) active data; ii) passive data; and iii) joint dataset. Active database is composed by a total of 16,950 air-gun shots during 1 month and passive database includes 452 local and regional earthquakes recorded during 4 months. This large dataset provides a high ray density within the study region. The combination of active and passive seismic data, together with the high quality of the database, permits to obtain a new tomographic approach of the region under study never done before. An additional user-guide of PARTOS software is provided in order to facilitate the implementation for new users.

[1]  J. Virieux,et al.  Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy) , 2008 .

[2]  Jean Virieux,et al.  Precise, absolute earthquake location under Somma–Vesuvius volcano using a new three-dimensional velocity model , 2001 .

[3]  E. Gordeev,et al.  Rapid changes in magma storage beneath the Klyuchevskoy group of volcanoes inferred from time-dependent seismic tomography , 2013 .

[4]  N. Shapiro,et al.  Seismic Tomography of Volcanoes , 2015 .

[5]  D. Patanè,et al.  Intrusive mechanism of the 2008–2009 Mt. Etna eruption: Constraints by tomographic images and stress tensor analysis , 2012 .

[6]  Clifford H. Thurber,et al.  Earthquake locations and three‐dimensional crustal structure in the Coyote Lake Area, central California , 1983 .

[7]  Clifford H. Thurber,et al.  Double-Difference Tomography: The Method and Its Application to the Hayward Fault, California , 2003 .

[8]  P. Gori,et al.  Time-Resolved Seismic Tomography Detects Magma Intrusions at Mount Etna , 2006, Science.

[9]  S. Sobolev,et al.  Moho depth and three-dimensional P and S structure of the crust and uppermost mantle in the Eastern Mediterranean and Middle East derived from tomographic inversion of local ISC data , 2006 .

[10]  P. Izbekov,et al.  Seismic structure changes beneath Redoubt Volcano during the 2009 eruption inferred from local earthquake tomography , 2014 .

[11]  H. Thybo,et al.  Seismic tomographic imaging of P- and S-waves velocity perturbations in the upper mantle beneath Iran , 2007 .

[12]  N. Hirata,et al.  Non‐volcanic seismic swarms triggered by circulating fluids and pressure fluctuations above a solidified diorite intrusion , 2010 .

[13]  L. Cocchi,et al.  The marine activities performed within the TOMO-ETNA experiment , 2016 .

[14]  Sonia Mota,et al.  Advances on the automatic estimation of the P-wave onset time. , 2016 .

[15]  Sonia Mota,et al.  An Automatic P-Phase Picking Algorithm Based on Adaptive Multiband Processing , 2013, IEEE Geoscience and Remote Sensing Letters.

[16]  Andreas Rietbrock,et al.  High resolution 3D P wave velocity structure beneath Tenerife Island (Canary Islands, Spain) based on tomographic inversion of active‐source data , 2012 .

[17]  Sean C. Solomon,et al.  Tomographic imaging of the shallow crustal structure of the East Pacific Rise at 9°30′N , 1994 .

[18]  J. Ibáñez,et al.  Seismic and volcanic activity during 2014 in the region involved by TOMO-ETNA seismic active experiment , 2016 .

[19]  S. Sobolev,et al.  A tomographic image of Indian lithosphere break‐off beneath the Pamir–Hindukush region , 2006 .

[20]  J. Ibáñez,et al.  Seismic tomography model reveals mantle magma sources of recent volcanic activity at El Hierro Island (Canary Islands, Spain) , 2014 .

[21]  P. Izbekov,et al.  Fluid ascent during the 2004–2005 unrest at Mt. Spurr inferred from seismic tomography , 2013 .

[22]  P. Kuznetsov,et al.  Evidence of magma activation beneath the Harrat Lunayyir basaltic field (Saudi Arabia) from attenuation tomography , 2014 .

[23]  Colin A. Zelt,et al.  Three‐dimensional seismic refraction tomography: A comparison of two methods applied to data from the Faeroe Basin , 1998 .

[24]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[25]  Mrinal K. Sen,et al.  Upper crust of the Archean Dharwar craton in southern India as revealed by seismic refraction tomography and its geotectonic implications , 2014 .

[26]  Guust Nolet,et al.  Linearized Inversion of (Teleseismic) Data , 1981 .

[27]  J. Ibáñez TOMO-ETNA experiment at Etna volcano: Activities on land , 2016 .

[28]  Apostolos Sarris,et al.  Geophysical investigations at the Istron archaeological site, eastern Crete, Greece using seismic refraction and electrical resistivity tomography , 2012 .

[29]  Nicholas Rawlinson,et al.  Simultaneous inversion of active and passive source datasets for 3‐D seismic structure with application to Tasmania , 2006 .

[30]  H. Kopp,et al.  Joint inversion of active and passive seismic data in Central Java , 2007 .

[31]  C. Vargas,et al.  Three-dimensional velocity structure of the Galeras volcano (Colombia) from passive local earthquake tomography , 2015 .

[32]  V. Chebrov,et al.  Magma sources in the mantle wedge beneath the volcanoes of the Klyuchevskoy group and Kizimen based on seismic tomography modeling , 2016 .

[33]  I. Koulakov Multiscale Seismic Tomography Imaging of Volcanic Complexes , 2012 .

[34]  M. I. Seht,et al.  Evidence for high fluid/melt content beneath Krakatau volcano (Indonesia) from local earthquake tomography , 2010 .

[35]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .

[36]  J. Ibáñez,et al.  The TOMO-ETNA experiment: an imaging active campaign at Mt. Etna volcano. Context, main objectives, working-plans and involved research projects , 2016 .

[37]  W. Rabbel,et al.  Local earthquake tomography of central Costa Rica: transition from seamount to ridge subduction , 2010 .

[38]  W. Rabbel,et al.  Shear-wave velocity structure of the Chilean subduction zone (39–40°S) based on Rayleigh wave dispersion: evidence of fluid release and melts in the mantle beneath the Villarrica volcano , 2014 .

[39]  E. Gordeev,et al.  Feeding volcanoes of the Kluchevskoy group from the results of local earthquake tomography , 2011 .

[40]  P. Kuznetsov,et al.  The three-dimensional structure beneath the Popocatépetl volcano (Mexico) based on local earthquake seismic tomography , 2014 .

[41]  Ivan Koulakov,et al.  LOTOS Code for Local Earthquake Tomographic Inversion: Benchmarks for Testing Tomographic Algorithms , 2009 .

[42]  Daria Zandomeneghi,et al.  Crustal structure of Deception Island volcano from P wave seismic tomography: Tectonic and volcanic implications , 2009 .

[43]  M. Lindenfeld,et al.  Seismic images of magmatic rifting beneath the western branch of the East African rift , 2013 .