Three‐Dimensional Electrical Resistivity Tomography of the Solfatara Crater (Italy): Implication for the Multiphase Flow Structure of the Shallow Hydrothermal System

The Solfatara volcano is the main degassing area of the Campi Flegrei caldera, characterized by 60 years of unrest. Assessing such renewal activity is a challenging task because hydrothermal interactions with magmatic gases remain poorly understood. In this study, we decipher the complex structure of the shallow Solfatara hydrothermal system by performing the first 3-D, high-resolution, Electrical Resistivity Tomography (ERT) of the volcano. The 3-D resistivity model was obtained from the inversion of 43,432 resistance measurements performed on an area of 0.68 km2. The proposed interpretation of the multiphase hydrothermal structures is based on the resistivity model, a high-resolution infrared surface temperature image, and 1,136 soil CO2 flux measurements. In addition, we realized 27 soil Cation Exchange Capacity (CEC) and pH measurements demonstrating a negligible contribution of surface conductivity to the shallow bulk electrical conductivity. Hence, we show that the resistivity changes are mainly controlled by fluid content and temperature. The high-resolution tomograms identify for the first time the structure of the gas-dominated reservoir at 50 m depth that feeds the Bocca Grande fumarole through a ~10-m-thick channel. In addition, the resistivity model reveals a channel-like conductive structure where the liquid produced by steam condensation around the main fumaroles flows down to the Fangaia area within a buried fault. The model delineates the emplacement of the main geological structures: Mt Olibano, Solfatara crypto-dome, and tephra deposits. It also reveals the anatomy of the hydrothermal system, especially two liquid-dominated plumes, the Fangaia mud pool and the Pisciarelli fumarole, respectively.

[1]  E. Guidoboni,et al.  The Campi Flegrei caldera: historical revision and new data on seismic crises, bradyseisms, the Monte Nuovo eruption and ensuing earthquakes (twelfth century 1582 ad) , 2011 .

[2]  S. Vitale,et al.  Stratigraphy, structure, and volcano-tectonic evolution of Solfatara maar-diatreme (Campi Flegrei, Italy) , 2015 .

[3]  N. Pearce,et al.  Tephrostratigraphy and glass compositions of post-15 kyr Campi Flegrei eruptions: implications for eruption history and chronostratigraphic markers , 2011 .

[4]  A. Nur,et al.  Ultrasonic velocity measurements in volcanic rocks: correlation with microtexture , 2002 .

[5]  G. E. Archie.  The electrical resistivity log as an aid in determining some reservoir characteristics , 1942 .

[6]  A. Revil,et al.  Induced polarization of volcanic rocks – 1. Surface versus quadrature conductivity , 2017 .

[7]  G. Rolandi,et al.  New constraints on the pyroclastic eruptive history of the Campanian volcanic Plain (Italy) , 2001 .

[8]  Carlo Cardellini,et al.  Carbon dioxide Earth degassing and seismogenesis in central and southern Italy , 2004 .

[9]  G. Chiodini,et al.  Geochemical evidences of magma dynamics at Campi Flegrei (Italy) , 2014 .

[10]  M. Karaoulis,et al.  Image-guided inversion of electrical resistivity data , 2014 .

[11]  Massimo D'Antonio,et al.  Volcanism and deformation since 12,000 years at the Campi Flegrei caldera (Italy) , 1999 .

[12]  K. Aizawaa,et al.  Hydrothermal system beneath Mt . Fuji volcano inferred from magnetotellurics and electric self-potential , 2005 .

[13]  W. Evans,et al.  Deep magmatic degassing versus scrubbing: Elevated CO2emissions and C/S in the lead‐up to the 2009 eruption of Redoubt Volcano, Alaska , 2012 .

[14]  M. G. Di Giuseppe,et al.  Electromagnetic outline of the Solfatara–Pisciarelli hydrothermal system, Campi Flegrei (Southern Italy) , 2014 .

[15]  R. Fournier.  Hydrothermal systems and volcano geochemistry , 2007 .

[16]  M. Uyeshima,et al.  Hydrothermal system beneath Mt. Fuji Volcano inferred from magnetotellurics and electric self-potential , 2005 .

[17]  N. Fournier,et al.  Seismo-acoustic evidence for an avalanche driven phreatic eruption through a beheaded hydrothermal system: An example from the 2012 Tongariro eruption , 2014 .

[18]  Carlo Cardellini,et al.  First observations of the fumarolic gas output from a restless caldera: Implications for the current period of unrest (2005–2013) at Campi Flegrei , 2013 .

[19]  Douglas LaBrecque,et al.  Difference Inversion of ERT Data: a Fast Inversion Method for 3-D In Situ Monitoring , 2001 .

[20]  P. Marianelli,et al.  Caldera unrest prior to intense volcanism in Campi Flegrei (Italy) at 4.0 ka B.P.: Implications for caldera dynamics and future eruptive scenarios , 2009 .

[21]  G. Chiodini,et al.  Geophysical and hydrogeological experiments from a shallow hydrothermal system at Solfatara Volcano, Campi Flegrei, Italy: Response to caldera unrest , 2007 .

[22]  M. H. Waxman,et al.  Electrical Conductivities in Oil-Bearing Shaly Sands , 1968 .

[23]  Luigi Marini,et al.  CO2 degassing and energy release at Solfatara Volcano , 2001 .

[24]  E. Rizzo,et al.  Inner structure of La Fossa di Vulcano (Vulcano Island, southern Tyrrhenian Sea, Italy) revealed by high-resolution electric resistivity tomography coupled with self-potential, temperature, and CO2 diffuse degassing measurements , 2008 .

[25]  T. Ishido.  Electrokinetic mechanism for the “W”‐shaped self‐potential profile on volcanoes , 2004 .

[26]  A. Revil,et al.  Alteration of volcanic rocks: A new non-intrusive indicator based on induced polarization measurements , 2017 .

[27]  Viviana Piermattei,et al.  Mapping Spatial Patterns of Posidonia oceanica Meadows by Means of Daedalus ATM Airborne Sensor in the Coastal Area of Civitavecchia (Central Tyrrhenian Sea, Italy) , 2013, Remote. Sens..

[28]  Hydrothermal system beneath Aso volcano as inferred from self-potential mapping and resistivity structure , 2005 .

[29]  A. Revil,et al.  Streaming electrical potential anomaly along faults in geothermal areas , 1998 .

[30]  S. Sudarman,et al.  Geothermal resources of Sumatra , 1993 .

[31]  G. Chiodini,et al.  Geochemical and isotopic changes in the fumarolic and submerged gas discharges during the 2011–2012 unrest at Santorini caldera (Greece) , 2013, Bulletin of Volcanology.

[32]  S. Byrdina,et al.  Relations between electrical resistivity, carbon dioxide flux, and self-potential in the shallow hydrothermal system of Solfatara (Phlegrean Fields, Italy) , 2014 .

[33]  A. Revil,et al.  Importance of structural history in the summit area of Stromboli during the 2002–2003 eruptive crisis inferred from temperature, soil CO2, self-potential, and electrical resistivity tomography , 2009 .

[34]  William H. Hendershot,et al.  Soil Reaction and Exchangeable Acidity , 2007 .

[35]  Valerio Acocella,et al.  An overview of recent (1988 to 2014) caldera unrest: Knowledge and perspectives , 2015 .

[36]  F. Flandoli,et al.  Quantifying volcanic hazard at Campi Flegrei caldera (Italy) with uncertainty assessment: 2. Pyroclastic density current invasion maps , 2015 .

[37]  Luigi Marini,et al.  Soil CO2 flux measurements in volcanic and geothermal areas , 1998 .

[38]  E. Spangenberg,et al.  Effect of the water–steam phase transition on the electrical conductivity of porous rocks , 2010 .

[39]  J. D. Mcneill.  Electromagnetic Terrain Conduc-tivity Measurement at Low Induction Numbers , 1980 .

[40]  N. Geshi,et al.  Growth of phreatomagmatic explosion craters: A model inferred from Suoana crater in Miyakejima Volcano, Japan , 2011 .

[41]  L. Civetta,et al.  The present state of the magmatic system of the Campi Flegrei caldera based on a reconstruction of its behavior in the past 12 ka , 1999 .

[42]  T. Günther,et al.  Constraining 3-D electrical resistance tomography with GPR reflection data for improved aquifer characterization , 2012 .

[43]  Roelof Versteeg,et al.  Improved hydrogeophysical characterization and monitoring through parallel modeling and inversion of time-domain resistivity andinduced-polarization data , 2010 .

[44]  G. Chiodini,et al.  Magmas near the critical degassing pressure drive volcanic unrest towards a critical state , 2016, Nature Communications.

[45]  Pier Paolo Bruno,et al.  3D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes , 2017, Scientific Reports.

[46]  R. Barker,et al.  Rapid least-squares inversion of apparent resistivity pseudosections , 1994 .

[47]  D. Kulik,et al.  Differential alteration of basaltic lava flows and hyaloclastites in Icelandic hydrothermal systems , 2015, Geothermal Energy.

[48]  M. Zamora,et al.  Subsurface geophysics of the Phlegrean Fields: New insights from downhole measurements , 2003 .

[49]  D. Dingwell,et al.  Hydrothermal alteration of surficial rocks at Solfatara (Campi Flegrei): Petrophysical properties and implications for phreatic eruption processes , 2016 .

[50]  J. Gottsmann,et al.  Oscillations in hydrothermal systems as a source of periodic unrest at caldera volcanoes: Multiparameter insights from Nisyros, Greece , 2007 .

[51]  B. Raco,et al.  Diffuse emission of CO2 from the Fossa crater, Vulcano Island (Italy) , 1996 .

[52]  Carlo Cardellini,et al.  Long-term variations of the Campi Flegrei, Italy, volcanic system as revealed by the monitoring of hydrothermal activity , 2010 .

[53]  Jean Virieux,et al.  Seismic reflections reveal a massive melt layer feeding Campi Flegrei caldera , 2008 .

[54]  G. Chiodini.  CO2/CH4 ratio in fumaroles a powerful tool to detect magma degassing episodes at quiescent volcanoes , 2009 .

[55]  M. L. Sorey,et al.  Vapor‐Dominated Zones Within Hydrothermal Systems: Evolution and Natural State , 1988 .

[56]  Rosario Avino,et al.  Early signals of new volcanic unrest at Campi Flegrei caldera? Insights from geochemical data and physical simulations , 2012 .

[57]  C. Cardellini,et al.  Geochemical indicators of possible ongoing volcanic unrest at Nisyros Island (Greece) , 2002 .

[58]  R. Barker,et al.  Rapid least-squared inversion of apparent resisitivity pseudosections by a quasi-Newton method , 1996 .

[59]  A. Revil,et al.  Induced polarization of volcanic rocks. 2. Influence of pore size and permeability , 2017 .

[60]  M. Ripepe,et al.  Geophysical investigations at Stromboli volcano, Italy. Implications for ground water flow and paroxysmal activity. , 2004 .

[61]  Aldo Zollo,et al.  Geophysical Exploration of the Campi Flegrei (Southern Italy) Caldera' Interiors: Data, Methods and Results , 2006 .

[62]  D. Gibert,et al.  Volcano electrical tomography unveils edifice collapse hazard linked to hydrothermal system structure and dynamics , 2016, Scientific Reports.

[63]  Luca D'Auria,et al.  Evidence of thermal-driven processes triggering the 2005–2014 unrest at Campi Flegrei caldera , 2015 .

[64]  Donald E. White,et al.  Vapor-dominated hydrothermal systems compared with hot-water systems , 1971 .

[65]  S. Byrdina,et al.  Geophysical image of the hydrothermal system of Merapi volcano , 2017 .

[66]  Anil V. Rao,et al.  GPOPS-II , 2014, ACM Trans. Math. Softw..

[67]  P. Gasparini,et al.  History of earthquakes and vertical ground movement in Campi Flegrei caldera, Southern Italy: comparison of precursory events to the A.D. 1538 eruption of Monte Nuovo and of activity since 1968 , 1991 .

[68]  A. Zollo,et al.  A strongly heterogeneous hydrothermal area imaged by surface waves: the case of Solfatara, Campi Flegrei, Italy , 2016 .

[69]  André Revil,et al.  The Self-Potential Method: Theory and Applications in Environmental Geosciences , 2018 .

[70]  Maurizio Ripepe,et al.  The summit hydrothermal system of Stromboli. New insights from self-potential, temperature, CO2 and fumarolic fluid measurements, with structural and monitoring implications , 2003 .

[71]  G. Chiodini,et al.  Carbon dioxide diffuse degassing and estimation of heat release from volcanic and hydrothermal systems , 2005 .

[72]  R. Versteeg,et al.  Characterization of a contaminated wellfield using 3D electrical resistivity tomography implemented with geostatistical, discontinuous boundary, and known conductivity constraints , 2012 .

[73]  Hang Si,et al.  TetGen, a Delaunay-Based Quality Tetrahedral Mesh Generator , 2015, ACM Trans. Math. Softw..

[74]  S. Parolai,et al.  A 3-D algorithm based on the combined inversion of Rayleigh and Love waves for imaging and monitoring of shallow structures , 2017 .

[75]  J. Thouret,et al.  Fluid circulation and structural discontinuities inside Misti volcano (Peru) inferred from self-potential measurements , 2004 .

[76]  A. Revil,et al.  Electrical properties of zeolitized volcaniclastic materials , 2002 .

[77]  P. Agustín,et al.  Hydrothermal system of Central Tenerife Volcanic Complex, Canary Islands (Spain), inferred from self-potential measurements , 2014 .

[78]  C. Harvey,et al.  UNDERSTANDING THE RESISTIVITIES OBSERVED IN GEOTHERMAL SYSTEMS , 2000 .

[79]  M. Rosi,et al.  The phlegraean fields: Structural evolution, volcanic history and eruptive mechanisms , 1983 .

[80]  G. Orsi,et al.  Volcanic hazard assessment at the restless Campi Flegrei caldera , 2004 .

[81]  S. Byrdina,et al.  Hydrogeology of Stromboli volcano, Aeolian Islands (Italy) from the interpretation of resistivity tomograms, self-potential, soil temperature and soil CO2 concentration measurements , 2011 .

[82]  E. Del Pezzo,et al.  Subsurface structure of the Solfatara volcano (Campi Flegrei caldera, Italy) as deduced from joint seismic‐noise array, volcanological and morphostructural analysis , 2012 .

[83]  Brian P. Bonner,et al.  The effects of capillarity on electrical resistivity during boiling in metashale from scientific corehole SB-15-D, The Geysers, California, USA , 2001 .

[84]  G. P. Ricciardi,et al.  Unrest episodes at Campi Flegrei: A reconstruction of vertical ground movements during 1905-2009 , 2010 .

[85]  A. Maul,et al.  A spectrophotometric measurement of soil cation exchange capacity based on cobaltihexamine chloride absorbance , 2008 .

[86]  André Revil,et al.  Three‐dimensional resistivity tomography of Vulcan's forge, Vulcano Island, southern Italy , 2010 .