3D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes
Identifieur interne : 000005 ( Pmc/Corpus ); précédent : 000004; suivant : 0000063D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes
Auteurs : Grazia De Landro ; Vincenzo Serlenga ; Guido Russo ; Ortensia Amoroso ; Gaetano Festa ; Pier Paolo Bruno ; Marceau Gresse ; Jean Vandemeulebrouck ; Aldo ZolloSource :
- Scientific Reports [ 2045-2322 ] ; 2017.
Abstract
Seismic tomography can be used to image the spatial variation of rock properties within complex geological media such as volcanoes. Solfatara is a volcano located within the Campi Flegrei, a still active caldera, so it is of major importance to characterize its level of activity and potential danger. In this light, a 3D tomographic high-resolution P-wave velocity image of the shallow central part of Solfatara crater is obtained using first arrival times and a
Url:
DOI: 10.1038/s41598-017-03604-0
PubMed: 28611382
PubMed Central: 5469761
Links to Exploration step
PMC:5469761Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">3D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes</title>
<author><name sortKey="De Landro, Grazia" sort="De Landro, Grazia" uniqKey="De Landro G" first="Grazia" last="De Landro">Grazia De Landro</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Serlenga, Vincenzo" sort="Serlenga, Vincenzo" uniqKey="Serlenga V" first="Vincenzo" last="Serlenga">Vincenzo Serlenga</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.466609.b</institution-id>
<institution></institution>
<institution>Now at Consiglio Nazionale delle Ricerche, Istituto di Metodologie per l’Analisi Ambientale,</institution>
</institution-wrap>
Tito, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Russo, Guido" sort="Russo, Guido" uniqKey="Russo G" first="Guido" last="Russo">Guido Russo</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Amoroso, Ortensia" sort="Amoroso, Ortensia" uniqKey="Amoroso O" first="Ortensia" last="Amoroso">Ortensia Amoroso</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Festa, Gaetano" sort="Festa, Gaetano" uniqKey="Festa G" first="Gaetano" last="Festa">Gaetano Festa</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bruno, Pier Paolo" sort="Bruno, Pier Paolo" uniqKey="Bruno P" first="Pier Paolo" last="Bruno">Pier Paolo Bruno</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1796 0286</institution-id>
<institution-id institution-id-type="GRID">grid.452820.9</institution-id>
<institution></institution>
<institution>Petroleum Institute, Department of Petroleum Geosciences,</institution>
</institution-wrap>
Abu Dhabi, United Arab Emirates</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gresse, Marceau" sort="Gresse, Marceau" uniqKey="Gresse M" first="Marceau" last="Gresse">Marceau Gresse</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.5388.6</institution-id>
<institution></institution>
<institution>ISTerre, Université Savoie Mont Blanc,</institution>
</institution-wrap>
Chambéry, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Vandemeulebrouck, Jean" sort="Vandemeulebrouck, Jean" uniqKey="Vandemeulebrouck J" first="Jean" last="Vandemeulebrouck">Jean Vandemeulebrouck</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.5388.6</institution-id>
<institution></institution>
<institution>ISTerre, Université Savoie Mont Blanc,</institution>
</institution-wrap>
Chambéry, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zollo, Aldo" sort="Zollo, Aldo" uniqKey="Zollo A" first="Aldo" last="Zollo">Aldo Zollo</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">28611382</idno>
<idno type="pmc">5469761</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC5469761</idno>
<idno type="RBID">PMC:5469761</idno>
<idno type="doi">10.1038/s41598-017-03604-0</idno>
<date when="2017">2017</date>
<idno type="wicri:Area/Pmc/Corpus">000005</idno>
<idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000005</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">3D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes</title>
<author><name sortKey="De Landro, Grazia" sort="De Landro, Grazia" uniqKey="De Landro G" first="Grazia" last="De Landro">Grazia De Landro</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Serlenga, Vincenzo" sort="Serlenga, Vincenzo" uniqKey="Serlenga V" first="Vincenzo" last="Serlenga">Vincenzo Serlenga</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="Aff4"><institution-wrap><institution-id institution-id-type="GRID">grid.466609.b</institution-id>
<institution></institution>
<institution>Now at Consiglio Nazionale delle Ricerche, Istituto di Metodologie per l’Analisi Ambientale,</institution>
</institution-wrap>
Tito, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Russo, Guido" sort="Russo, Guido" uniqKey="Russo G" first="Guido" last="Russo">Guido Russo</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Amoroso, Ortensia" sort="Amoroso, Ortensia" uniqKey="Amoroso O" first="Ortensia" last="Amoroso">Ortensia Amoroso</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Festa, Gaetano" sort="Festa, Gaetano" uniqKey="Festa G" first="Gaetano" last="Festa">Gaetano Festa</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Bruno, Pier Paolo" sort="Bruno, Pier Paolo" uniqKey="Bruno P" first="Pier Paolo" last="Bruno">Pier Paolo Bruno</name>
<affiliation><nlm:aff id="Aff2"><institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1796 0286</institution-id>
<institution-id institution-id-type="GRID">grid.452820.9</institution-id>
<institution></institution>
<institution>Petroleum Institute, Department of Petroleum Geosciences,</institution>
</institution-wrap>
Abu Dhabi, United Arab Emirates</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Gresse, Marceau" sort="Gresse, Marceau" uniqKey="Gresse M" first="Marceau" last="Gresse">Marceau Gresse</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.5388.6</institution-id>
<institution></institution>
<institution>ISTerre, Université Savoie Mont Blanc,</institution>
</institution-wrap>
Chambéry, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Vandemeulebrouck, Jean" sort="Vandemeulebrouck, Jean" uniqKey="Vandemeulebrouck J" first="Jean" last="Vandemeulebrouck">Jean Vandemeulebrouck</name>
<affiliation><nlm:aff id="Aff3"><institution-wrap><institution-id institution-id-type="GRID">grid.5388.6</institution-id>
<institution></institution>
<institution>ISTerre, Université Savoie Mont Blanc,</institution>
</institution-wrap>
Chambéry, France</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Zollo, Aldo" sort="Zollo, Aldo" uniqKey="Zollo A" first="Aldo" last="Zollo">Aldo Zollo</name>
<affiliation><nlm:aff id="Aff1"><institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">Scientific Reports</title>
<idno type="eISSN">2045-2322</idno>
<imprint><date when="2017">2017</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>Seismic tomography can be used to image the spatial variation of rock properties within complex geological media such as volcanoes. Solfatara is a volcano located within the Campi Flegrei, a still active caldera, so it is of major importance to characterize its level of activity and potential danger. In this light, a 3D tomographic high-resolution P-wave velocity image of the shallow central part of Solfatara crater is obtained using first arrival times and a <italic>multiscale</italic>
approach. The retrieved images, integrated with the resistivity section and temperature and the CO<sub>2</sub>
flux measurements, define the following characteristics: 1. A depth-dependent P-wave velocity layer down to 14 m, with V<sub>p</sub>
< 700 m/s typical of poorly-consolidated tephra and affected by CO<sub>2</sub>
degassing; 2. An intermediate layer, deepening towards the mineralized liquid-saturated area (Fangaia), interpreted as permeable deposits saturated with condensed water; 3. A deep, confined high velocity anomaly associated with a CO<sub><bold>2</bold>
</sub>
reservoir. These features are expression of an area located between the Fangaia, water saturated and replenished from deep aquifers, and the main fumaroles, superficial relief of the deep rising CO<sub>2</sub>
flux. Therefore, the changes in the outgassing rate greatly affect the shallow hydrothermal system, which can be used as a “mirror” of fluid migration processes occurring at depth.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="De Vivo, B" uniqKey="De Vivo B">B De Vivo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Bianco, F" uniqKey="Bianco F">F Bianco</name>
</author>
<author><name sortKey="Del Pezzo, E" uniqKey="Del Pezzo E">E Del Pezzo</name>
</author>
<author><name sortKey="Saccorotti, G" uniqKey="Saccorotti G">G Saccorotti</name>
</author>
<author><name sortKey="Ventura, G" uniqKey="Ventura G">G Ventura</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Orsi, G" uniqKey="Orsi G">G Orsi</name>
</author>
<author><name sortKey="De Vita, S" uniqKey="De Vita S">S De Vita</name>
</author>
<author><name sortKey="Di Vito, M" uniqKey="Di Vito M">M di Vito</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zollo, A" uniqKey="Zollo A">A Zollo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aiuppa, A" uniqKey="Aiuppa A">A Aiuppa</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Natale, G" uniqKey="De Natale G">G De Natale</name>
</author>
<author><name sortKey="Pingue, F" uniqKey="Pingue F">F Pingue</name>
</author>
<author><name sortKey="Allard, P" uniqKey="Allard P">P Allard</name>
</author>
<author><name sortKey="Zollo, A" uniqKey="Zollo A">A Zollo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Del Gaudio, C" uniqKey="Del Gaudio C">C Del Gaudio</name>
</author>
<author><name sortKey="Aquino, I" uniqKey="Aquino I">I Aquino</name>
</author>
<author><name sortKey="Ricciardi, Gp" uniqKey="Ricciardi G">GP Ricciardi</name>
</author>
<author><name sortKey="Ricco, C" uniqKey="Ricco C">C Ricco</name>
</author>
<author><name sortKey="Scandone, R" uniqKey="Scandone R">R Scandone</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Saccorotti, G" uniqKey="Saccorotti G">G Saccorotti</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Gaeta, Fs" uniqKey="Gaeta F">FS Gaeta</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Battaglia, M" uniqKey="Battaglia M">M Battaglia</name>
</author>
<author><name sortKey="Troise, C" uniqKey="Troise C">C Troise</name>
</author>
<author><name sortKey="Obrizzo, F" uniqKey="Obrizzo F">F Obrizzo</name>
</author>
<author><name sortKey="De Natale, G" uniqKey="De Natale G">G De Natale</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cusano, P" uniqKey="Cusano P">P Cusano</name>
</author>
<author><name sortKey="Petrosino, P" uniqKey="Petrosino P">P Petrosino</name>
</author>
<author><name sortKey="Saccorotti, G" uniqKey="Saccorotti G">G Saccorotti</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Chiodini, G" uniqKey="Chiodini G">G Chiodini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Todesco, M" uniqKey="Todesco M">M Todesco</name>
</author>
<author><name sortKey="Rinaldi, Ap" uniqKey="Rinaldi A">AP Rinaldi</name>
</author>
<author><name sortKey="Bonafede, M" uniqKey="Bonafede M">M Bonafede</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chiodini, G" uniqKey="Chiodini G">G Chiodini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Petrillo, Z" uniqKey="Petrillo Z">Z Petrillo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Coco, A" uniqKey="Coco A">A Coco</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Afanasyev, A" uniqKey="Afanasyev A">A Afanasyev</name>
</author>
<author><name sortKey="Costa, A" uniqKey="Costa A">A Costa</name>
</author>
<author><name sortKey="Chiodini, G" uniqKey="Chiodini G">G Chiodini</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Serra, M" uniqKey="Serra M">M Serra</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Di Stefano, R" uniqKey="Di Stefano R">R Di Stefano</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Amoroso, O" uniqKey="Amoroso O">O Amoroso</name>
</author>
<author><name sortKey="Ascione, A" uniqKey="Ascione A">A Ascione</name>
</author>
<author><name sortKey="Mazzoli, S" uniqKey="Mazzoli S">S Mazzoli</name>
</author>
<author><name sortKey="Virieux, J" uniqKey="Virieux J">J Virieux</name>
</author>
<author><name sortKey="Zollo, A" uniqKey="Zollo A">A Zollo</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Zucca, Jj" uniqKey="Zucca J">JJ Zucca</name>
</author>
<author><name sortKey="Hutchings, Lj" uniqKey="Hutchings L">LJ Hutchings</name>
</author>
<author><name sortKey="Kasameyer, Pw" uniqKey="Kasameyer P">PW Kasameyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Husen, S" uniqKey="Husen S">S Husen</name>
</author>
<author><name sortKey="Smith, Rb" uniqKey="Smith R">RB Smith</name>
</author>
<author><name sortKey="Waite, Gp" uniqKey="Waite G">GP Waite</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Vanorio, T" uniqKey="Vanorio T">T Vanorio</name>
</author>
<author><name sortKey="Virieux, J" uniqKey="Virieux J">J Virieux</name>
</author>
<author><name sortKey="Capuano, P" uniqKey="Capuano P">P Capuano</name>
</author>
<author><name sortKey="Russo, G" uniqKey="Russo G">G Russo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="De Lorenzo, S" uniqKey="De Lorenzo S">S de Lorenzo</name>
</author>
<author><name sortKey="Gasparini, P" uniqKey="Gasparini P">P Gasparini</name>
</author>
<author><name sortKey="Mongelli, F" uniqKey="Mongelli F">F Mongelli</name>
</author>
<author><name sortKey="Zollo, A" uniqKey="Zollo A">A Zollo</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Bruno, Pp" uniqKey="Bruno P">PP Bruno</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Isaia, R" uniqKey="Isaia R">R Isaia</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Harris, Ajl" uniqKey="Harris A">AJL Harris</name>
</author>
<author><name sortKey="Carniel, R" uniqKey="Carniel R">R Carniel</name>
</author>
<author><name sortKey="Jones, J" uniqKey="Jones J">J Jones</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Caliro, S" uniqKey="Caliro S">S Caliro</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aiuppa, A" uniqKey="Aiuppa A">A Aiuppa</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Aster, Rc" uniqKey="Aster R">RC Aster</name>
</author>
<author><name sortKey="Meyer, Rp" uniqKey="Meyer R">RP Meyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Tramelli, A" uniqKey="Tramelli A">A Tramelli</name>
</author>
<author><name sortKey="Del Pezzo, E" uniqKey="Del Pezzo E">E Del Pezzo</name>
</author>
<author><name sortKey="Bianco, F" uniqKey="Bianco F">F Bianco</name>
</author>
<author><name sortKey="Boschi, E" uniqKey="Boschi E">E Boschi</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Caudill, M" uniqKey="Caudill M">M Caudill</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Latorre, D" uniqKey="Latorre D">D Latorre</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zollo, A" uniqKey="Zollo A">A Zollo</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lutter, Wj" uniqKey="Lutter W">WJ Lutter</name>
</author>
<author><name sortKey="Nowack, R" uniqKey="Nowack R">R Nowack</name>
</author>
<author><name sortKey="Braile, Lw" uniqKey="Braile L">LW Braile</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Akaike, H" uniqKey="Akaike H">H Akaike</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Eberhart Phillips, D" uniqKey="Eberhart Phillips D">D Eberhart-Phillips</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Michelini, A" uniqKey="Michelini A">A Michelini</name>
</author>
<author><name sortKey="Mcevilly, Tv" uniqKey="Mcevilly T">TV McEvilly</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Caliro, S" uniqKey="Caliro S">S Caliro</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Sissons, Ba" uniqKey="Sissons B">BA Sissons</name>
</author>
<author><name sortKey="Dibble, Rr" uniqKey="Dibble R">RR Dibble</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Mora, M" uniqKey="Mora M">M Mora</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Toksoz, Nm" uniqKey="Toksoz N">NM Toksöz</name>
</author>
<author><name sortKey="Cheng, Ch" uniqKey="Cheng C">CH Cheng</name>
</author>
<author><name sortKey="Timur, A" uniqKey="Timur A">A Timur</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Loke, Mh" uniqKey="Loke M">MH Loke</name>
</author>
<author><name sortKey="Dahlin, T" uniqKey="Dahlin T">T Dahlin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Loke, Mh" uniqKey="Loke M">MH Loke</name>
</author>
<author><name sortKey="Barker, Rd" uniqKey="Barker R">RD Barker</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Revil, A" uniqKey="Revil A">A Revil</name>
</author>
<author><name sortKey="Florsch, N" uniqKey="Florsch N">N Florsch</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Chiodini, G" uniqKey="Chiodini G">G Chiodini</name>
</author>
<author><name sortKey="Frondini, F" uniqKey="Frondini F">F Frondini</name>
</author>
<author><name sortKey="Raco, B" uniqKey="Raco B">B Raco</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Letort, J" uniqKey="Letort J">J Letort</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Sci Rep</journal-id>
<journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id>
<journal-title-group><journal-title>Scientific Reports</journal-title>
</journal-title-group>
<issn pub-type="epub">2045-2322</issn>
<publisher><publisher-name>Nature Publishing Group UK</publisher-name>
<publisher-loc>London</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">28611382</article-id>
<article-id pub-id-type="pmc">5469761</article-id>
<article-id pub-id-type="publisher-id">3604</article-id>
<article-id pub-id-type="doi">10.1038/s41598-017-03604-0</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Article</subject>
</subj-group>
</article-categories>
<title-group><article-title>3D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" corresp="yes"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0003-1136-1357</contrib-id>
<name><surname>De Landro</surname>
<given-names>Grazia</given-names>
</name>
<address><email>delandro@fisica.unina.it</email>
</address>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Serlenga</surname>
<given-names>Vincenzo</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
<xref ref-type="aff" rid="Aff4">4</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0001-9268-8647</contrib-id>
<name><surname>Russo</surname>
<given-names>Guido</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Amoroso</surname>
<given-names>Ortensia</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Festa</surname>
<given-names>Gaetano</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0003-2622-3037</contrib-id>
<name><surname>Bruno</surname>
<given-names>Pier Paolo</given-names>
</name>
<xref ref-type="aff" rid="Aff2">2</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Gresse</surname>
<given-names>Marceau</given-names>
</name>
<xref ref-type="aff" rid="Aff3">3</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Vandemeulebrouck</surname>
<given-names>Jean</given-names>
</name>
<xref ref-type="aff" rid="Aff3">3</xref>
</contrib>
<contrib contrib-type="author"><contrib-id contrib-id-type="orcid">http://orcid.org/0000-0002-8191-9566</contrib-id>
<name><surname>Zollo</surname>
<given-names>Aldo</given-names>
</name>
<xref ref-type="aff" rid="Aff1">1</xref>
</contrib>
<aff id="Aff1"><label>1</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0001 0790 385X</institution-id>
<institution-id institution-id-type="GRID">grid.4691.a</institution-id>
<institution>Department of physics “Ettore Pancini”,</institution>
<institution>University of Naples ”Federico II”,</institution>
</institution-wrap>
Naples, Italy</aff>
<aff id="Aff2"><label>2</label>
<institution-wrap><institution-id institution-id-type="ISNI">0000 0004 1796 0286</institution-id>
<institution-id institution-id-type="GRID">grid.452820.9</institution-id>
<institution></institution>
<institution>Petroleum Institute, Department of Petroleum Geosciences,</institution>
</institution-wrap>
Abu Dhabi, United Arab Emirates</aff>
<aff id="Aff3"><label>3</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.5388.6</institution-id>
<institution></institution>
<institution>ISTerre, Université Savoie Mont Blanc,</institution>
</institution-wrap>
Chambéry, France</aff>
<aff id="Aff4"><label>4</label>
<institution-wrap><institution-id institution-id-type="GRID">grid.466609.b</institution-id>
<institution></institution>
<institution>Now at Consiglio Nazionale delle Ricerche, Istituto di Metodologie per l’Analisi Ambientale,</institution>
</institution-wrap>
Tito, Italy</aff>
</contrib-group>
<pub-date pub-type="epub"><day>13</day>
<month>6</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="pmc-release"><day>13</day>
<month>6</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection"><year>2017</year>
</pub-date>
<volume>7</volume>
<elocation-id>3412</elocation-id>
<history><date date-type="received"><day>14</day>
<month>10</month>
<year>2016</year>
</date>
<date date-type="accepted"><day>9</day>
<month>5</month>
<year>2017</year>
</date>
</history>
<permissions><copyright-statement>© The Author(s) 2017</copyright-statement>
<license license-type="OpenAccess"><license-p><bold>Open Access</bold>
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link>
.</license-p>
</license>
</permissions>
<abstract id="Abs1"><p>Seismic tomography can be used to image the spatial variation of rock properties within complex geological media such as volcanoes. Solfatara is a volcano located within the Campi Flegrei, a still active caldera, so it is of major importance to characterize its level of activity and potential danger. In this light, a 3D tomographic high-resolution P-wave velocity image of the shallow central part of Solfatara crater is obtained using first arrival times and a <italic>multiscale</italic>
approach. The retrieved images, integrated with the resistivity section and temperature and the CO<sub>2</sub>
flux measurements, define the following characteristics: 1. A depth-dependent P-wave velocity layer down to 14 m, with V<sub>p</sub>
< 700 m/s typical of poorly-consolidated tephra and affected by CO<sub>2</sub>
degassing; 2. An intermediate layer, deepening towards the mineralized liquid-saturated area (Fangaia), interpreted as permeable deposits saturated with condensed water; 3. A deep, confined high velocity anomaly associated with a CO<sub><bold>2</bold>
</sub>
reservoir. These features are expression of an area located between the Fangaia, water saturated and replenished from deep aquifers, and the main fumaroles, superficial relief of the deep rising CO<sub>2</sub>
flux. Therefore, the changes in the outgassing rate greatly affect the shallow hydrothermal system, which can be used as a “mirror” of fluid migration processes occurring at depth.</p>
</abstract>
<custom-meta-group><custom-meta><meta-name>issue-copyright-statement</meta-name>
<meta-value>© The Author(s) 2017</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body><sec id="Sec1" sec-type="introduction"><title>Introduction</title>
<p>Solfatara is a tuff cone, which was formed between 3.8 and 4.1 ky ago<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
, located 1.5 km NE of the town of Pozzuoli and about 10 km west of the city of Naples, in Southern Italy. The sub-rectilinear NE and SW rims are cut by two normal faults that strike NW-SE, along which deep geothermal fluids can ascend. Outside the crater, two NW-SE striking faults cut the eastern part of the tuff cone<sup><xref ref-type="bibr" rid="CR2">2</xref>
</sup>
.</p>
<p>Solfatara is one of the many volcanoes located within the Campi Flegrei caldera. This is a nested, resurgent caldera, resulting from two large collapses related to the Campanian Ignimbrite (39 ka<sup><xref ref-type="bibr" rid="CR1">1</xref>
</sup>
) and the Neapolitan Yellow Tuff (14 ka<sup><xref ref-type="bibr" rid="CR3">3</xref>
</sup>
) eruptions (Orsi <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR3">3</xref>
</sup>
). The Campi Flegrei volcanic system is still active, since the last eruption occurred in 1538 A.D. at Monte Nuovo.</p>
<p>Zollo <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR4">4</xref>
</sup>
identified a deep, sill-like, mid-crustal magmatic body supplying heat to this volcanic system at a depth of 7–8 km. They have also identified another shallower interface at 2.5 km depth, associated with a discontinuity between the older caldera deposits and a fluid-saturated metamorphic rock layer. Concerning the magmatic source characteristics, Aiuppa <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR5">5</xref>
</sup>
proposed two different source models based upon the measures of Campi Flegrei fumarolic gas output. The two main degassing areas are the Solfatara and Pisciarelli, the latter being located 100 m East of the Solfatara. In particular, the authors estimated that the current Campi Flegrei fumarolic sulphur flux is low, whereas the fumarolic CO<sub>2</sub>
flux is surprisingly high for a dormant volcano in the hydrothermal stage of activity. Thus, they proposed that the current CO<sub>2</sub>
output can be supplied by either a large (0.6–4.6 km<sup>3</sup>
), deeply stored (>7 km) magmatic source with low CO<sub>2</sub>
contents (0.05–0.1 wt%) or by a small to medium-sized (0.01–0.1 km<sup>3</sup>
) but CO<sub>2</sub>
-rich (2 wt%) magma, possibly stored at pressures of 100 to 120 MPa.</p>
<p>The caldera has been characterized by periodic episodes of extended, low-rate ground subsidence and uplift, a phenomenon called bradyseism, accompanied by intense seismic and geochemical activity<sup><xref ref-type="bibr" rid="CR6">6</xref>
</sup>
. During the past century, three main episodes of bradyseism occurred in the area: 1950–1952, 1969–1972, and 1982–1984<sup><xref ref-type="bibr" rid="CR7">7</xref>
</sup>
. During the subsidence following the 1984 crisis, a series of small uplift episodes and seismic swarms occurred at the Solfatara. Saccorotti <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR8">8</xref>
</sup>
suggest that these earthquakes are likely to be associated with the vibration of a buried cavity filled with a water-vapour mixture at poor gas-volume fractions. The role of fluids in bradyseism has been recognized by many authors [e.g. refs <xref ref-type="bibr" rid="CR5">5</xref>
, <xref ref-type="bibr" rid="CR9">9</xref>
–<xref ref-type="bibr" rid="CR11">11</xref>
].</p>
<p>The Solfatara crater is characterized by intense diffuse degassing and fumarolic emissions [e.g. refs <xref ref-type="bibr" rid="CR12">12</xref>
–<xref ref-type="bibr" rid="CR14">14</xref>
]. Chiodini <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
showed that Solfatara releases about 1500 t/day of volcanic-hydrothermal CO<sub>2</sub>
, as a result of diffuse degassing through soil, during which about 3350 t/day of steam condense. This hydrothermal water mixes up with the meteoric one. The energy released by degassing at Solfatara is much higher than the energy released within the caldera during the current period, through other processes such as thermal conduction, earthquakes, and ground deformation. Moreover, the impressive magnitude of diffuse degassing process confirmed the relevance of fluid and heat transport at Solfatara<sup><xref ref-type="bibr" rid="CR15">15</xref>
</sup>
and prompted further research to improve the understanding of the hydrothermal system feeding the surface phenomenon.</p>
<p>Although the magmatic source is rather deep, the changes in its outgassing rate greatly affect the shallow hydrothermal system processes, which can therefore be used as a constraint for the fluid migration processes occurring at depth<sup><xref ref-type="bibr" rid="CR16">16</xref>
, <xref ref-type="bibr" rid="CR17">17</xref>
</sup>
.</p>
<p>The effect of deep processes on the shallow hydrothermal system was confirmed by Afanasyev <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR18">18</xref>
</sup>
through a simulation of hydrothermal activity at Campi Flegrei caldera. The authors showed that the physical quantities characterizing the deep source, i.e. gas injection rate, gas mixture composition and temperature, strongly control the shallow gas composition, temperature, CO<sub>2</sub>
flux and heat flux monitored at Solfatara.</p>
<p>This justifies the growing interest of the last years in delineating the physical properties of the shallow hydrothermal system of Solfatara, as one of the means to assess the level of potential danger characterizing this crater of the Campi Flegrei volcanic complex.</p>
<p>In the framework of the experiment RICEN, Repeated Induced Earthquake and Noise (EU Project MEDSUV), an active seismic experiment was carried out between September of 2013 and November of 2014 in order to provide time-varying high-resolution images of the structure of the Solfatara<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
. A grid of 240 receivers, placed at a distance of 5 m on 10 lines, which were 10 m distant from each other, was deployed in the Solfatara crater. A vibroseis energized the soil at the centre of almost all grid cells (Fig. <xref rid="Fig1" ref-type="fig">1</xref>
). A large and highly informative data-set was then built. The data used in this work were acquired during the first campaign of the experiment, which was carried out on September of 2013.<fig id="Fig1"><label>Figure 1</label>
<caption><p>Solfatara with station-shot configuration. (<bold>a</bold>
) A photo of the RICEN experiment, with an example of seismic sensors. (<bold>b</bold>
) Station-shot configuration. (<bold>c</bold>
) The Solfatara crater, with evidence of crater rims, fumarole, Fangaia and the station-shot configuration. The maps in Fig. 1b,c has been obtained with Google Maps 9.38.1 2016 (Map data: Google, DigitalGlobe): <italic>Solfatara, Pozzuoli, Metropolitan City of Naples, Italy</italic>
retrieved from <ext-link ext-link-type="uri" xlink:href="https://www.google.com/maps/@40.8174278,14.1393587,1292a,20y,41.32t/data=!3m1!1e3">https://www.google.com/maps/@40.8174278,14.1393587,1292a,20y,41.32t/data=!3m1!1e3</ext-link>
.</p>
</caption>
<graphic xlink:href="41598_2017_3604_Fig1_HTML" id="d29e456"></graphic>
</fig>
</p>
<p>Velocity and attenuation seismic tomography can be used to image the spatial variation of elastic/anelastic rock properties within complex geological media such as volcanoes. This approach has been widely used in complex tectonic<sup><xref ref-type="bibr" rid="CR20">20</xref>
–<xref ref-type="bibr" rid="CR23">23</xref>
</sup>
, geothermal<sup><xref ref-type="bibr" rid="CR24">24</xref>
–<xref ref-type="bibr" rid="CR26">26</xref>
</sup>
and volcanic areas<sup><xref ref-type="bibr" rid="CR27">27</xref>
–<xref ref-type="bibr" rid="CR31">31</xref>
</sup>
in order to detect the presence of fluids and track their migration within the fluid permeated crustal volume.</p>
<p>The goal of the present study is to obtain an ultra-high-resolution (metric) 3D velocity image of the shallow (down to 35 m depth) hydrothermal structure of the central part of Solfatara crater using a technique of delay-time tomography, based on the P-wave first arrival times. Due to the extremely dense acquisition lay-out, the 3D tomographic survey allows us to achieve an unprecedented spatial detail on the shallow velocity structure, which can help in understanding the complex hydrothermal degassing and condensation processes into the porous rocks media. In this way, the interpretation obtained by stratigraphic analysis and resistivity profiles<sup><xref ref-type="bibr" rid="CR32">32</xref>
, <xref ref-type="bibr" rid="CR33">33</xref>
</sup>
can be complemented. For our purposes, the results of temperature and CO<sub>2</sub>
flux measurements and of the resistivity survey carried out in the frame of the RICEN experiment will be used to further constrain our interpretation. In fact, several studies over the last few years showed how the multi-parametric analysis of combined geophysical and/or geochemical data sets has been helpful to investigate the complex dynamics of volcanic systems at different scales<sup><xref ref-type="bibr" rid="CR34">34</xref>
–<xref ref-type="bibr" rid="CR36">36</xref>
</sup>
.</p>
<p>In the following, we will briefly describe the data processing and tomographic method, then we will present and analyse the retrieved P-wave velocity model and assess uncertainty and spatial resolution. Finally, we will compare the velocity model with a high-resolution electrical resistivity tomography performed in May of 2014 to investigate and assess the mechanisms of hydrothermal fluid circulation.</p>
<sec id="Sec2"><title>Previous works</title>
<p>Bruno <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR32">32</xref>
</sup>
imaged the shallow and intermediate subsurface of Solfatara through the integration of high-resolution geophysical and hydrogeological investigations, including 2D P-wave velocity and Electrical Resistivity Tomography (ERT) profiles. Their results show that the Solfatara subsurface structure can be roughly divided into two zones: a dry, outcropping layer and an underlying saturated zone, in which faults and fractures act as preferential escape conduits for the hydrothermal fluids [refs <xref ref-type="bibr" rid="CR1">1</xref>
, <xref ref-type="bibr" rid="CR32">32</xref>
, and <xref ref-type="bibr" rid="CR33">33</xref>
].</p>
<p>Byrdina <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
used the results from the ERT survey, the mappings of diffuse CO<sub>2</sub>
flux, the ground temperature and self-potential (SP) to understand the mechanisms and paths of shallow fluid circulation. They interpreted the resistivity changes at depth, associated with surface gas flux anomalies, as a double-plume structure: a liquid-dominated conductive plume below the Fangaia mud-pool and a gas-dominated plume below the Bocca Grande fumarole (see Fig. <xref rid="Fig1" ref-type="fig">1</xref>
).</p>
<p>Isaia <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
integrated the electrical resistivity tomography investigations with the volcano-tectonic information to better constrain the subsurface structure by outlining a complex hydrothermal system. In particular, they assume that the 100 m thick upper zone of Solfatara comprises desegregated rocks and collapse breccias, post-eruptive sediments, whereas the lower sector, down to about 3–4 km, is where the gas-saturated conduit is connected to a magmatic source.</p>
<p>Referring to the Campi Flegrei, de Lorenzo <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR30">30</xref>
</sup>
and De Siena <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR31">31</xref>
</sup>
carried out two studies on the anelastic properties of the northern part of the Campi Flegrei caldera, including the Solfatara crater. Both studies retrieved a heterogeneous distribution of low-Qp and high-Qp anomalies in the investigated area, thus confirming the strong geological complexity of the subsoil. Moving closer to Solfatara, the two works reveal different features: De Siena <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR31">31</xref>
</sup>
found a high Qp and low Qs body. Their correlation with the low Vp/Vs ratio and the low Vp was interpreted as the effect of a small gas reservoir. de Lorenzo <italic>et al</italic>
,<sup><xref ref-type="bibr" rid="CR30">30</xref>
</sup>
. on the other hand, found low Qp values at shallow depths, well correlated with high Vp/Vs value<sup><xref ref-type="bibr" rid="CR38">38</xref>
</sup>
; this correlation was interpreted as produced by densely fractured, porous and fluid-filled rocks. A similar interpretation was also provided by Tramelli <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR39">39</xref>
</sup>
who found a high-scattering zone in the area of the Solfatara. A recent attenuation tomography<sup><xref ref-type="bibr" rid="CR40">40</xref>
</sup>
mainly focused on the shallowest subsurface of the Campi Flegrei offshore caldera. The retrieved anelastic images are not spatially correlated with Solfatara; however, they describe the very shallow volcanic system as an environment which is greatly affected by the heterogeneous distribution of different saturation conditions of fluids permeating sediments and rocks of caldera.</p>
<p>Using a subset of the data-bank obtained through the RICEN experiment, Serra <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
found the spatial variation of surface wave phase and group velocities. By inverting the related dispersion curves they obtained a one-dimensional S-wave model for different sub-grids. Together, the different 1-D S-wave models provided a three-dimensional description of the S-wave model in the area down to about 15 m depth. In the upper 4 m, they associated the changes of the S-wave velocity to the temperature gradient, while at greater depths the seismic images were correlated with the resistivity maps obtained from the measurements carried out during the RICEN experiment. They evidenced the presence of a water layer close to the Fangaia area (see Fig. <xref rid="Fig1" ref-type="fig">1</xref>
) and an abrupt variation in NE direction.</p>
</sec>
<sec id="Sec3"><title>Data and method</title>
<p>The RICEN experiment consisted in three successive geophysical surveys carried out at the Solfatara volcano respectively in September 2013, May and November 2014, each one lasting one week. During each experiment, the recording both of active seismic data and of the continuous ambient noise was performed. A dataset with more than 75,000 seismograms was collected during the active seismic part of the three experiments. Active seismic data were obtained using a Vibroseis Truck soil energizator, which operated in the frequency range 5–125 Hz. Seismic waveforms were recorded by 4.5 Hz vertical component geophones (GS-11D, Fig. <xref rid="Fig1" ref-type="fig">1a</xref>
).</p>
<p>In this study, we analyse the data collected during the first experiment. In particular, an area of 90 × 115 m<sup>2</sup>
was sampled by a regular grid of 240 vertical sensors, which were deployed at the crater surface (Fig. <xref rid="Fig1" ref-type="fig">1b,c</xref>
). The seismic network geometry was set up according to a two-dimensional (2D) grid with 10 lines of 24 sensors, with 5 m spacing between the stations (i.e. in-line distance). The distance between two adjacent lines (i.e. cross-line distance) was 10 m. About 100 shot-points were energized on a staggered grid relatively to the receiver grid. For the vibrational sources, both the in-line and cross-line inter-distances were 10 m. For each shot position, three consecutive energizations were performed and waveforms at each site were stacked in order to increase the signal-to-noise ratio.</p>
<p>It is well known that this volcanic area is very complex, where effects of scattering and attenuation could in general contaminate the first arrival picking. In order to check the quality requirements for an accurate picking analysis, we verified <italic>a posteriori</italic>
that the signal-to-noise ratio was sufficiently high at the first arrival onset, and we evaluated waveform/picking coherence vs. offset along the seismic section (Fig. <xref rid="Fig2" ref-type="fig">2a</xref>
). Through the analysis of the first arrival incidence angles and horizontal polarization directions measured at nearby installed, three-component stations, we confirmed the P-wave nature of the first arriving phase and validated the use of the vertical component for time picking, notwithstanding the strong medium heterogeneity and very short distance between source and receiver In particular, by considering 97 shot records, about the 90% of three-component waveforms showed an incidence angle less than 40° (60% less than 25°) in the first arrival time window, which is the evidence of dominant, near-vertical first arrivals. Moreover, both the horizontal polarization analysis, circular wave-front tracing and apparent velocity estimations, showed a first arrival phase back-azimuth strongly consistent with the expected P-wave azimuth and velocity given the relative location of the shots and the receivers (for further details see Supplementary Materials).<fig id="Fig2"><label>Figure 2</label>
<caption><p>Data and travel-time-vs-distance plot. (<bold>a</bold>
) Example of seismic section in common shot gather. The processed traces concerns the shot 194 recorded at stations 25 to 48. The green dots represents picks of first arrival P phases obtained with PROMAX and manually revised. (<bold>b</bold>
) Plot of travel-time-vs-distance (red dots). The black lines are the limit imposed for the construction of 1D initial velocity models. Blue line is the 1D initial model selected for the inversion.</p>
</caption>
<graphic xlink:href="41598_2017_3604_Fig2_HTML" id="d29e645"></graphic>
</fig>
</p>
<p>The P-wave first arrivals have been first detected through a Neural Network implemented into the ProMAX SeisSpace software<sup><xref ref-type="bibr" rid="CR41">41</xref>
–<xref ref-type="bibr" rid="CR44">44</xref>
</sup>
trained on a limited, manually picked dataset of source-receiver couples. In that way about 19,000 picked P-wave travel times have been obtained (for further details check the Supplementary Materials, “Data processing” paragraph). Thereafter, the picking dataset has been manually validated on the basis of a visual inspection of seismic signals. Data affected by a high pick of uncertainty (>0.05 s), based on a very low signal-to-noise ratio (<2), have been discarded. An example of one seismic section in the common-shot gather can be found in Fig. <xref rid="Fig2" ref-type="fig">2a</xref>
. Moreover, for each shot the travel times referring to the eight closest receivers have been excluded. In those cases, indeed, it was too difficult to correctly measure the P-wave travel time, due to the short pre-event time window available for picking. All data have been considered with the same weight, since only the best quality pickings have been considered for the tomographic analysis. The final selected dataset contained 17,418 P wave travel times relative to 94 sources and 240 receivers.</p>
<p>In order to determine the 3-D P-wave velocity model, the selected travel-time dataset was inverted by applying an iterative, linearized, tomographic approach<sup><xref ref-type="bibr" rid="CR45">45</xref>
</sup>
and by adopting a <italic>multiscale</italic>
procedure<sup><xref ref-type="bibr" rid="CR46">46</xref>
</sup>
. Several inversion runs were performed by progressively increasing the density of grid points (e.g. increasing the model complexity) describing the velocity field, and at each iteration the starting model is the one estimated in the previous run. The proposed procedure is equivalent in principle to moving from a low to a high wave number description of the velocity field. The <italic>multiscale</italic>
strategy does not depend on the scale of the application, because its basic assumption is that large wave-length anomalies in the velocity structure have a dominant amplitude when compared to the smaller ones, which we believe is reasonable at all scales of investigation within Earth. This strategy has been already used in seismic tomography<sup><xref ref-type="bibr" rid="CR46">46</xref>
, <xref ref-type="bibr" rid="CR47">47</xref>
</sup>
and in migration techniques, where it is known as “multi-scale” approach. The application of this strategy for our analysis is justified by the observed travel-time curves (T(X)) and the distribution of arrival time data vs. offset around their mean trend (Fig. <xref rid="Fig2" ref-type="fig">2b</xref>
). Indeed, it is clear that even at the small scale of the Solfatara crater, the long wavelength (low-frequency) component of the velocity model (reflected by the piece-wise linear increase in T(X)) has a dominant amplitude relative to short wavelength (high frequency) anomalies possibly producing the dispersion of arrival times around the average trend. We therefore attributed them to both the uncertainty on data picking and to small wavelength perturbations of the velocity model. Specific tests with several inversion, runs with different medium parameterization, confirmed ‘a posteriori’ that the adopted multi-scale strategy was efficient and robust to explore the multi-dimensional model parameter space and to catch the minimum norm model solution (for further details see Supplementary Materials, “Inversion procedure” paragraph).</p>
<p>Due to the source-receiver configuration we have investigated a volume of 160 × 160 × 45 m<sup>3</sup>
, the top being at 100 m a.s.l. The three-dimensional hosting medium has been discretized with a grid of regularly spaced nodes.</p>
<p>The final parameterization of the medium in the multiscale procedure is chosen by applying the corrected Akaike Criterion [AICc; ref. <xref ref-type="bibr" rid="CR48">48</xref>
], which is based on a statistical comparison between models characterized by a different number of model parameters. By introducing the minimum AICc criterion<sup><xref ref-type="bibr" rid="CR49">49</xref>
</sup>
, the problem of selecting the optimal model parameterization is solved avoiding a subjective decision. The minimum AICc value, representing the best compromise between data misfit reduction and model simplicity, is obtained with the 10 × 10 × 5 m<sup>3</sup>
grid spacing (for further details see Supplementary Materials, “Inversion procedure” paragraph), which is the final parameterization in the <italic>multiscale</italic>
procedure.</p>
<p>At first, for the <italic>multiscale</italic>
approach, a coarser parameterization was used with a node-spacing of 16×16 × 7.5 m<sup>3</sup>
.</p>
<p>In order to take into account the possible dependence of the final tomographic solution on the starting velocity model we adopted the following strategy: (1) first, we estimated a reference 1D velocity model; (2) starting from this, we generated a set of 200 1D initial models; (3) a 3D inversion is performed for each initial model.</p>
<p>The estimation of the 1D reference velocity model is obtained by a modelling procedure, which minimized the RMS of the travel time residuals (see Supplementary Materials). This procedure allows us to assess that the simplest velocity model minimizing the travel-time residuals is two-layered structure. As for the 200 1D initial models, they have been constructed as 2-layer models with the depth of the interface fixed to that of the reference 1D velocity model (10 m). The variability range of the velocities in the 2 layers was established by constraining the respective theoretical travel-time curve to lie within the limits defined by the scattering of measured data in a travel-time-vs.-distance graph (see Fig. <xref rid="Fig2" ref-type="fig">2b</xref>
and Supplementary Fig. <xref rid="MOESM1" ref-type="media">S8b</xref>
). Each 1D velocity model is used as a starting model for the data inversion. Then, by computing the average of the 200 final three-dimensional velocity models and the normalized standard deviation for each model parameter (σ/<italic>V</italic>
<sub><italic>P</italic>
</sub>
), we observe that deviations from the average model on retrieved velocity values are less than 15%, except for some grid nodes (see Supplementary Fig. <xref rid="MOESM1" ref-type="media">S9</xref>
).</p>
<p>The starting velocity model minimizing the final misfit of the residuals is chosen as the initial velocity model for the inversion with the coarser parameterization (16 × 16 × 7.5 m<sup>3</sup>
) (see Fig. <xref rid="Fig2" ref-type="fig">2b</xref>
). The three-dimensional P-wave velocity model retrieved at the sixth iteration is used as starting model for further data inversion in the final grid (spacing 10 × 10 × 5 m<sup>3</sup>
).</p>
<p>For each parameterization we selected a damping parameter using an empirical approach, by performing several inversions with different values<sup><xref ref-type="bibr" rid="CR50">50</xref>
</sup>
. The selected damping parameters are the ones providing the best compromise between the variance reduction of the residuals and the increase of the solution variance. By inspecting the trade-off curves retrieved for each parameterization, the value of 0.5 was selected for all the parameterizations graph (see Supplementary Fig. <xref rid="MOESM1" ref-type="media">S10b-c</xref>
).</p>
<p>In order to assess the reliability of the final solution, we numerically computed the resolution matrix from which we extracted the RDE (Resolution of Diagonal Elements) and the spread function<sup><xref ref-type="bibr" rid="CR51">51</xref>
</sup>
. The results of the resolution analysis allowed us to assess that the resolved area is about 100 × 120 × 35 m<sup>3</sup>
(for further details see Supplementary Figs <xref rid="MOESM1" ref-type="media">S13</xref>
and <xref rid="MOESM1" ref-type="media">S14</xref>
).</p>
</sec>
</sec>
<sec id="Sec4" sec-type="results"><title>Results</title>
<p>In Fig. <xref rid="Fig3" ref-type="fig">3</xref>
we show the 3D P-wave velocity model obtained with the finest parameterization (10 × 10 × 5 m<sup>3</sup>
). Considering all the <italic>multiscale</italic>
steps, we achieved a reduction of the root mean square of travel times residuals (rms) of about 70%, with a final rms of 4 ms (see Supplementary Fig. <xref rid="MOESM1" ref-type="media">S11</xref>
).<fig id="Fig3"><label>Figure 3</label>
<caption><p>P-wave velocity model. (<bold>a</bold>
) Horizontal slice of P-wave velocity model at different depths. The black contour delimitates the resolved area, i.e. the area for which the tree resolution parameters (RDE, Sj and <italic>DWS</italic>
) are included in a threshold value. The threshold values of <italic>S</italic>
<sub><italic>j</italic>
</sub>
and DWS are chosen in order to obtain a similar contour, binding the RDE to be higher than 0.9. The grey regions in each slice represent areas not covered by rays. (<bold>b</bold>
) P-wave velocity model projected onto the SW-NE cross-sections located in Fig. 3a. The blurred regions in each slice represent areas not resolved.</p>
</caption>
<graphic xlink:href="41598_2017_3604_Fig3_HTML" id="d29e815"></graphic>
</fig>
</p>
<p>The resolution matrix and the spread function, shown in Supplementary Fig. <xref rid="MOESM1" ref-type="media">S12</xref>
, allow us to state that the final model is well resolved down to 30–35 m depth.</p>
<p>The plane view of the model (Fig. <xref rid="Fig3" ref-type="fig">3a</xref>
) shows a depth-increasing P-wave velocity structure within the shallower 13–14 m depth, with velocity values ranging from 200 to 700 m/s. Between 13 and 15 m depth an abrupt increase of velocity is observed, from 500 m/s to 1000 m/s. Moreover, around these depths the model shows a significant lateral variation. From 15 m to 22 m depth, in the western part of the model, a low velocity anomaly trending NW-SE is imaged. Its P-wave velocity is about 800 m/s and its areal extension decreases with depth. From 22 m to 30 m depth, a NW-SE interface separates two distinct zones characterized by different compressional velocity values. The Western depth section has a velocity ranging between about 900 m/s and 1200 m/s, and the Eastern one has higher velocity values of about 1200–1500 m/s. Finally, a well confined, lens-shaped, high P-wave velocity zone is imaged in the deepest part of the model, with velocity values higher than 1500 m/s.</p>
<p>The SW-NE cross-sections shown in Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
allow us to better delineate the shape and areal extension of the features which were already identified in the plane view representation. First of all, the intermediate layer, with velocity ranging between 1000 and 1100 m/s, progressively deepens toward SW. Moreover, it quickly rises again at the interface with the mud pools (Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
). Then, the low velocity anomaly, included in the previous layer, is more extended in the central part of the model (see section B-C-D in Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
), and reduces towards the northern and southern edges of the grid (see section A-E-F in Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
). Finally, the same behaviour is observed in the deep high velocity anomaly, which becomes thinner and thinner at the edges of the grid and completely disappears in the “F” section in Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
.</p>
</sec>
<sec id="Sec5" sec-type="discussion"><title>Discussion</title>
<p>In this study we obtained a new 3D, high-resolution image of P-wave velocity of the Solfatara crater through an inversion strategy based on a <italic>multiscale</italic>
approach. We suggest that the retrieved subsoil velocity images can be explained in terms of the interaction between structural patterns and degassing dynamics within shallow hydrothermal circulation cells<sup><xref ref-type="bibr" rid="CR32">32</xref>
</sup>
. According to the geochemical analysis of the Solfatara fumarolic gases performed by Caliro <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR52">52</xref>
</sup>
this degassed flux is composed by deep magmatic CO<sub>2</sub>
-rich fluids mixed with hydrothermal liquids of meteoric origin.</p>
<p>The Solfatara shallow stratigraphy consists of eruption deposits, <italic>tephra</italic>
, which are hydrothermally altered and mainly composed of alternating fine to coarse ash deposits with limited distribution, scoria layers, and lavas<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
. In detail, the first 10–15 m are composed by recent, unconsolidated deposits; beneath this layer, the deeper 20–30 m thick deposits have a dominant composition of Astroni tephra, i.e. sandwave, ash surge and fallout deposits with variable degrees of consolidation<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
. The retrieved velocity values, V<sub>p</sub>
< 1800 m/s, are consistent with the ones found in other volcanic areas for the tephra deposits<sup><xref ref-type="bibr" rid="CR53">53</xref>
, <xref ref-type="bibr" rid="CR54">54</xref>
</sup>
. In particular, in the first 10–15 m the P-wave velocity values range from 200 m/s to 700 m/s, which correspond to aerated tephra<sup><xref ref-type="bibr" rid="CR54">54</xref>
</sup>
. On the other hand, the higher velocity values in the deeper zone (V<sub>p</sub>
up to 1800 m/s) can be related to the tephra deposit, which are more consolidated and possibly saturated<sup><xref ref-type="bibr" rid="CR53">53</xref>
, <xref ref-type="bibr" rid="CR54">54</xref>
</sup>
.</p>
<p>Taking into account that the presence of fluids and their circulation may greatly affect the rock volume, and therefore the average compressional wave velocity, we expect that the tomographic images can constrain the possible location and phase of permeating fluids<sup><xref ref-type="bibr" rid="CR55">55</xref>
</sup>
. For this purpose, we compared our seismic tomography with 2D cross sections of resistivity (ρ) and with temperature and CO<sub>2</sub>
flux measurements. Serra <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
show how the integration of velocity images with resistivity tomography provides a more complete interpretation of the complex Solfatara system. A high-resolution electrical resistivity tomography was performed with 16, 115-m-long NW–SE profiles and 24, 75-m-long NE–SW profiles. The surveys were carried out in March and May 2014 during the RICEN experiment. We used a Wenner-Schlumberger configuration with a 5 m spacing between electrodes. In addition to these short profiles, we used the 1 km-long profile performed by Byrdina <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
(labeled Pr2 in their paper) in order to increase the resolution at depth. Resistivity data were filtered by removing values with a standard deviation exceeding 5%, or when the injected current was lower than 20 mA. The 3-D resistivity inversion was performed with 3028 filtered apparent resistivity measurements using RES3DINV software<sup><xref ref-type="bibr" rid="CR56">56</xref>
</sup>
which adopts a smoothness-constrained least-square algorithm<sup><xref ref-type="bibr" rid="CR57">57</xref>
</sup>
. The resistivity tomogram shown in Fig. <xref rid="Fig4" ref-type="fig">4a</xref>
has a 7.5% rms error after 5 iterations. The same section of the sensitivity map, which shows a very good resolution up to 25–30 m depth, has been provided in Supplementary Fig. <xref rid="MOESM1" ref-type="media">S15</xref>
. The electrical conductivity can be written as the sum of the surface conductivity, which is prominent in the case of clay-rich minerals produced by hydrothermal alteration, and bulk conductivity, which in our case mainly depends on hydrothermal fluid saturation and temperature<sup><xref ref-type="bibr" rid="CR58">58</xref>
</sup>
. The CO<sub>2</sub>
flux has been measured using the Accumulation Chamber Method [for details, see ref. <xref ref-type="bibr" rid="CR59">59</xref>
], whereas the soil temperature was recorded at 30 cm depth by using a type K thermocouple. Both measurements have been carried out at each electrode location during the electrical resistivity tomography survey. The studied area includes the second largest diffuse degassing spot at Solfatara after the fumarolic area, with CO<sub>2</sub>
flux values ranging from 700 to 10000 g.m<sup>−2</sup>
.day<sup>−1</sup>
. This high CO<sub>2</sub>
flux anomaly is accompanied by a significant thermal anomaly, with soil temperature up to 80 °C at 30 cm depth, due to the latent heat transfer that occurs during steam condensation at the surface.<fig id="Fig4"><label>Figure 4</label>
<caption><p>Comparison between temperature (black line) and CO<sub>2</sub>
flux (red line) in (<bold>a</bold>
), resistivity cross section in (<bold>b</bold>
) and P-wave velocity model projected onto the cross-section D in Fig. <xref rid="Fig3" ref-type="fig">3a</xref>
in (<bold>c</bold>
). (<bold>d</bold>
) Schematic representation of geological features and fluid flux direction, discussed in the interpretation.</p>
</caption>
<graphic xlink:href="41598_2017_3604_Fig4_HTML" id="d29e993"></graphic>
</fig>
</p>
<p>The resistivity values of the SW-NE cross-section range between 1.5 and 20 Ωm, thus indicating a globally conductive subsurface for Solfatara (Fig. <xref rid="Fig4" ref-type="fig">4b</xref>
), in agreement with the previous results<sup><xref ref-type="bibr" rid="CR33">33</xref>
, <xref ref-type="bibr" rid="CR37">37</xref>
</sup>
. The model clearly highlights a sharp horizontal contrast between a resistive structure (20 Ωm) beneath the eastern part of the cross section and a conductive body (<5 Ωm) on the western part, towards the liquid-saturated Fangaia mud pool.</p>
<p>With reference to Fig. <xref rid="Fig4" ref-type="fig">4c</xref>
, the superficial layer “I”, with a P-wave velocity lower than 700 m/s, can be related to the shallow resistive zone “I” in Fig. <xref rid="Fig4" ref-type="fig">4b</xref>
, with resistivity ρ > 12 Ωm, located at the top of the resistivity section. The same resistive upper layer was obtained by Byrdina <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
and can be interpreted as unsaturated recent deposits<sup><xref ref-type="bibr" rid="CR32">32</xref>
, <xref ref-type="bibr" rid="CR33">33</xref>
</sup>
. This layer is affected by a high CO<sub>2</sub>
degassing (Fig. <xref rid="Fig4" ref-type="fig">4a</xref>
), higher than 1000 g.m<sup>−2</sup>
d<sup>−1</sup>
, and by an elevated surface temperature (>40 °C, Fig. <xref rid="Fig4" ref-type="fig">4a</xref>
). A passive seismic survey performed by Letort <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR60">60</xref>
</sup>
identified this area as the source of the ambient seismic noise generated by hydrothermal processes. The 3D velocity sections (Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
) allow us to delineate this layer, which extends through the entire central area of the crater.</p>
<p>The second layer “II” shows P-velocity values ranging from 900 m/s to 1200 m/s and is about 10 m thick. It progressively deepens toward the SW part of the crater, becoming shallower at the interface with the mud pool (Fig. <xref rid="Fig4" ref-type="fig">4c</xref>
). This structure corresponds to the eastern part of Fangaia, an area saturated with mineralized liquids<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
. Here, the resistivity image shows at 10–30 m depth a body (“II”) saturated with conductive liquid (Fig. <xref rid="Fig4" ref-type="fig">4b</xref>
),with ρ < 5 Ω m, getting deeper in the Fangaia direction like the velocity layer “II”. The slope of the liquid-saturated body is explained by a pronounced influence of local topography<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
.</p>
<p>On these grounds, we suggest that steam condensation is produced on the eastern part (“IV” and above) in a gas dominated structure that is characterized by both high diffuse degassing and high surface temperature. This steam produces hot condensate water, which is channelled within the “II” layer and finally reaches the Fangaia mud pool in the western part.</p>
<p>This fluid directionality is also inferred from the Self-Potential mapping by Byrdina <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
Indeed, in their paper Self Potential anomalies decrease from −40 mV in the eastern part to −100 mV in the Fangaia area. In a volcanic area, the streaming potential is the main source of current. It is associated with the drag of the excess of charge in the pore water of fluids. Consequently, above the isoelectric point of pH (which is the case in the Fangaia survey, with a soil pH above 5), an advective flow can be inferred from a decrease of the streaming potential.</p>
<p>In Fig. <xref rid="Fig4" ref-type="fig">4c</xref>
the low velocity anomaly “III”, with values between 700 and 900 m/s, is characterized by the same conductive properties already described for the body “II” in Fig. <xref rid="Fig4" ref-type="fig">4a</xref>
. There, we expect to find a region saturated with the liquid flow of the Fangaia area. However, inasmuch this body has velocity values lower than the “II” layer, we suggest, in agreement with the lithology<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
, that it is composed by less consolidated tephra sediments than the layer “II”. Therefore the velocity variation between the “III” and the “II” zone may identify a gentle discontinuity in the degree of consolidation of rocks.</p>
<p>The shape of this anomaly, which has the maximum extension in the central area and decreases going towards the NW and SE directions (Fig. <xref rid="Fig3" ref-type="fig">3b</xref>
), could be due to the accumulation of deposit materials in a low of the most compact layer “II”.</p>
<p>The previous results are in agreement with what was found in Serra <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
The retrieved 3D S-wave velocity model, obtained from a combination of 1D velocity models in array sub-grids and resolved up to 15 m, showed two main domains: the S-W one, closest to the Fangaia, is slower than the N-E one. This feature is common to all the depths. In order to interpret the low S-wave velocity domain at SW, the authors qualitatively hypothesized a manifestation of an unconsolidated layer at shallow depths and a water aquifer at greater depths. The P-wave velocity model obtained in this work not only has a much higher resolution, but also covers an area that is about twice the one of the S-wave velocity model obtained by Serra <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
These features allowed us to be in agreement with the previously results, but also to considerably increase the knowledge and characterization of the deeper anomalies of the Solfatara hydrothermal system. The water aquifer, mentioned by Serra <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
is imaged and better characterized through the interpretation of layer “II”, located at 15–20 m depth.</p>
<p>Finally, we will try to interpret the high velocity anomaly “IV” in Fig. <xref rid="Fig4" ref-type="fig">4c</xref>
. It is located between 25 and 30 m depths and it constitutes a novelty with respect to the S-wave velocity model by Serra <italic>et al</italic>
.<sup><xref ref-type="bibr" rid="CR19">19</xref>
</sup>
since no resolution was achieved in this cited work. This high velocity anomaly is characterized by values that range between 1500 and 1800 m/s. The anomaly “IV” does not have a perfect match in the ERT section, but could be located in resistivity tomography at depth between the resistive and conductive bodies. The high-velocity body appears at the beginning of the degassing structure (Stufe di Nerone, SN in Fig. <xref rid="Fig1" ref-type="fig">1</xref>
) on the eastern border of Solfatara crater, where major NW-SE directed fractures<sup><xref ref-type="bibr" rid="CR33">33</xref>
</sup>
enhance the up-flow of the water steam and the CO<sub>2</sub>
flux. This interpretation is consistent with the high CO<sub>2</sub>
flux and ground temperature measured in the same area<sup><xref ref-type="bibr" rid="CR37">37</xref>
</sup>
.</p>
<p>From the resistivity images (Fig. <xref rid="Fig4" ref-type="fig">4b</xref>
) it can be clearly seen that the conduit of the rising gas plume on the NE direction is located where the “II” layer saturated with condensed water is thinner (Fig. <xref rid="Fig4" ref-type="fig">4c</xref>
). Therefore, the anomaly found in our velocity images can represent an area of gas accumulation, trapped by the liquid saturated “II” layer, located immediately above and not completely resolved by the resistivity images due to its location below a conductive area.</p>
<p>In order to summarise our interpretation, we provide in Fig. <xref rid="Fig4" ref-type="fig">4d</xref>
a schematic representation of the discussed features and the direction of gas and condensed water fluxes.</p>
<p>The velocity tomographic model obtained for the Solfatara allows us to image with a high resolution the shallow area of the crater and to understand the processes taking place inside. These processes are part of a complex dynamics triggered by the interaction between structural patterns and degassing within shallow hydrothermal circulation cells. Then, at the shallow investigated depths, we can see the effect of deep processes affecting the caldera, and our investigated area can be interpreted as the surface evidence of deeper processes.</p>
</sec>
<sec id="Sec6" sec-type="conclusion"><title>Conclusion</title>
<p>The importance of this work lies in these principal aspects:<list list-type="order"><list-item><p>the Solfatara crater represents one of the main pressure release areas of the entire Campi Flegrei volcanic system, considering the impressive magnitude of the diffuse degassing process. Hence, the interest in the knowledge of this area grows, especially with the aim of assessing the level of potential danger characterizing this crater;</p>
</list-item>
<list-item><p>the 3D tomographic survey allows us to achieve an unprecedented spatial detail on the shallow velocity structure of the central part of the Solfatara crater. The 3D high-resolution tomographic images allow us to better understand, in terms of velocity anomalies and fluid type, the complex hydrothermal processes into the shallow part (30–35 m) of the volcano;</p>
</list-item>
<list-item><p>the procedure used in this work represents a new multi-parametric approach that can be used in a volcanic environment; it shows how the interpretation of velocity tomographic images can be complemented with the ones obtained by stratigraphic analysis and resistivity profiles, and, most of all, how this joint interpretation leads to a more robust and reliable interpretation of complex hydrothermal system;</p>
</list-item>
<list-item><p>the complex interactions (deep fluids, hydrothermal system, geological structures) between deep and shallow sources and structures allow us to use the shallow hydrothermal system processes as a constrain for the fluid migration processes occurring at depth.</p>
</list-item>
</list>
</p>
</sec>
<sec sec-type="supplementary-material"><title>Electronic supplementary material</title>
<sec id="Sec7"><p><supplementary-material content-type="local-data" id="MOESM1"><media xlink:href="41598_2017_3604_MOESM1_ESM.pdf"><caption><p>Supplementary Material</p>
</caption>
</media>
</supplementary-material>
</p>
</sec>
</sec>
</body>
<back><fn-group><fn><p><bold>Electronic supplementary material</bold>
</p>
<p><bold>Supplementary information</bold>
accompanies this paper at doi:10.1038/s41598-017-03604-0
</p>
</fn>
<fn><p><bold>Publisher's note:</bold>
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.</p>
</fn>
</fn-group>
<ack><title>Acknowledgements</title>
<p>This work has been a part of the Med-Suv project. MED-SUV has received funding from the European Union’s Seventh Program for research, technological development and demonstration under the grant agreement No. 308665. We also thank the CNR- IAMC for providing instrumentation and personnel for the survey. The authors wish to thank the editor Antonio Costa and the reviewers, including J.F. Lenat, for their valuable comments, which have contributed to improve the quality of the manuscript. Moreover, we thank all the members of the seismological laboratory of the Department of Physics (University of Naples Federico II) who have taken part in the RICEN experiment and, in particular, A. Caruso for fruitful discussions and comments on polarization analysis.</p>
</ack>
<notes notes-type="author-contribution"><title>Author Contributions</title>
<p>G.D.L. and V.S. worked at the inversions obtaining the 3D tomographic models. O.A. was involved in the preliminary parameterization test. G.R. obtained the 1D velocity model. G.D.L., V.S. and G.R. wrote the paper. A.Z. supervised and coordinated the different phases of the tomographic data analysis and inversion. G.D.L., V.S., G.R., G.F. and A.Z. were involved in model analysis and interpretation of results. M.G. and J.V. provide the resistivity image and were involved in the interpretation of results. P.P.B. performed the data processing and provided the first pick data-set. All co-authors were involved in the review of the manuscript.</p>
</notes>
<notes notes-type="COI-statement"><sec id="FPar1"><title>Competing Interests</title>
<p>The authors declare that they have no competing interests.</p>
</sec>
</notes>
<ref-list id="Bib1"><title>References</title>
<ref id="CR1"><label>1.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>De Vivo</surname>
<given-names>B</given-names>
</name>
<etal></etal>
</person-group>
<article-title>New constraints on the pyroclastic eruptive history of the Campanian volcanic Plain (Italy)</article-title>
<source>Mineralogy and Petrology.</source>
<year>2001</year>
<volume>73</volume>
<fpage>47</fpage>
<lpage>65</lpage>
<pub-id pub-id-type="doi">10.1007/s007100170010</pub-id>
</element-citation>
</ref>
<ref id="CR2"><label>2.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bianco</surname>
<given-names>F</given-names>
</name>
<name><surname>Del Pezzo</surname>
<given-names>E</given-names>
</name>
<name><surname>Saccorotti</surname>
<given-names>G</given-names>
</name>
<name><surname>Ventura</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>The role of hydrothermal fluids in triggering the July–August 2000 seismic swarm at Phlegraean Fields, Italy: evidence from seismological and mesostructural data</article-title>
<source>Journal of Volcanology and Geothermal Research.</source>
<year>2004</year>
<volume>133</volume>
<fpage>229</fpage>
<lpage>246</lpage>
<pub-id pub-id-type="doi">10.1016/S0377-0273(03)00400-1</pub-id>
</element-citation>
</ref>
<ref id="CR3"><label>3.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Orsi</surname>
<given-names>G</given-names>
</name>
<name><surname>De Vita</surname>
<given-names>S</given-names>
</name>
<name><surname>di Vito</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>The restless, resurgent Phlegraean Fields nested caldera (Italy): constraints on its evolution and configuration</article-title>
<source>Journal of Volcanology and Geothermal Research.</source>
<year>1996</year>
<volume>74</volume>
<issue>3–4</issue>
<fpage>179</fpage>
<lpage>214</lpage>
<pub-id pub-id-type="doi">10.1016/S0377-0273(96)00063-7</pub-id>
</element-citation>
</ref>
<ref id="CR4"><label>4.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zollo</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Seismic reflections reveal a massive melt layer feeding Phlegraean Fields caldera</article-title>
<source>Geophys. Res. Lett.</source>
<year>2008</year>
<volume>35</volume>
<fpage>L12306</fpage>
<pub-id pub-id-type="doi">10.1029/2008GL034242</pub-id>
</element-citation>
</ref>
<ref id="CR5"><label>5.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aiuppa</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>First observations of the fumarolic gas output from a restless caldera: Implications for the current period of unrest (2005–2013) at Campi Flegrei</article-title>
<source>Geochemistry, Geophysics, Geosystems.</source>
<year>2013</year>
<volume>14</volume>
<issue>10</issue>
<fpage>4153</fpage>
<lpage>4169</lpage>
<pub-id pub-id-type="doi">10.1002/ggge.20261</pub-id>
</element-citation>
</ref>
<ref id="CR6"><label>6.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>De Natale</surname>
<given-names>G</given-names>
</name>
<name><surname>Pingue</surname>
<given-names>F</given-names>
</name>
<name><surname>Allard</surname>
<given-names>P</given-names>
</name>
<name><surname>Zollo</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Geophysical and geochemical modelling of the 1982–1984 unrest phenomena at Phlegraean Fields caldera (southern Italy)</article-title>
<source>Journal of Volcanology and Geothermal Research.</source>
<year>1991</year>
<volume>48</volume>
<issue>1–2</issue>
<fpage>199</fpage>
<lpage>222</lpage>
<pub-id pub-id-type="doi">10.1016/0377-0273(91)90043-Y</pub-id>
</element-citation>
</ref>
<ref id="CR7"><label>7.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Del Gaudio</surname>
<given-names>C</given-names>
</name>
<name><surname>Aquino</surname>
<given-names>I</given-names>
</name>
<name><surname>Ricciardi</surname>
<given-names>GP</given-names>
</name>
<name><surname>Ricco</surname>
<given-names>C</given-names>
</name>
<name><surname>Scandone</surname>
<given-names>R</given-names>
</name>
</person-group>
<article-title>Unrest episodes at Phlegraean Fields: A reconstruction of vertical ground movements during 1905–2009</article-title>
<source>Journal of Volcanology and Geothermal Research.</source>
<year>2010</year>
<volume>195</volume>
<issue>1</issue>
<fpage>48</fpage>
<lpage>56</lpage>
<pub-id pub-id-type="doi">10.1016/j.jvolgeores.2010.05.014</pub-id>
</element-citation>
</ref>
<ref id="CR8"><label>8.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Saccorotti</surname>
<given-names>G</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Seismicity associated with the 2004–2006 renewed ground uplift at Phlegraean Fields Caldera, Italy</article-title>
<source>Physics of the Earth and Planetary Interiors.</source>
<year>2007</year>
<volume>165</volume>
<issue>1–2</issue>
<fpage>14</fpage>
<lpage>24</lpage>
<pub-id pub-id-type="doi">10.1016/j.pepi.2007.07.006</pub-id>
</element-citation>
</ref>
<ref id="CR9"><label>9.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gaeta</surname>
<given-names>FS</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Genesis and evolution of unrest episodes at Phlegraean Fields caldera: the role of the thermal fluid dynamical processes in the geothermal system</article-title>
<source>J. Geophys. Res.</source>
<year>1998</year>
<volume>103</volume>
<issue>20</issue>
<fpage>921</fpage>
<lpage>20,933</lpage>
</element-citation>
</ref>
<ref id="CR10"><label>10.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Battaglia</surname>
<given-names>M</given-names>
</name>
<name><surname>Troise</surname>
<given-names>C</given-names>
</name>
<name><surname>Obrizzo</surname>
<given-names>F</given-names>
</name>
<name><surname>De Natale</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Evidence for fluid migration as the source of deformation at Phlegraean Fields caldera (Italy)</article-title>
<source>Geophys. Res. Lett.</source>
<year>2006</year>
<volume>33</volume>
<fpage>L01307</fpage>
<pub-id pub-id-type="doi">10.1029/2005GL024904</pub-id>
</element-citation>
</ref>
<ref id="CR11"><label>11.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cusano</surname>
<given-names>P</given-names>
</name>
<name><surname>Petrosino</surname>
<given-names>P</given-names>
</name>
<name><surname>Saccorotti</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Hydrothermal origin for sustained Long-Period (LP) activity at Phlegraean Fields volcanic complex, Italy</article-title>
<source>J. Volcanol. Geotherm. Res.</source>
<year>2008</year>
<volume>177</volume>
<fpage>1035</fpage>
<lpage>1044</lpage>
<pub-id pub-id-type="doi">10.1016/j.jvolgeores.2008.07.019</pub-id>
</element-citation>
</ref>
<ref id="CR12"><label>12.</label>
<mixed-citation publication-type="other">Chiodini, G. <italic>et al</italic>
. Carbon dioxide diffuse degassing and estimation of heat release from volcanic and hydrothermal systems. <italic>Journal of Geophysical Research: Solid Earth</italic>
. <bold>110, </bold>
doi:10.1029/2004JB003542 (2005).</mixed-citation>
</ref>
<ref id="CR13"><label>13.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chiodini</surname>
<given-names>G</given-names>
</name>
<etal></etal>
</person-group>
<article-title>CO2/CH4 ratio in fumaroles a powerful tool to detect magma degassing episodes at quiescent volcanoes</article-title>
<source>Geophys. Res. Lett.</source>
<year>2009</year>
<volume>36</volume>
<fpage>L02302</fpage>
<pub-id pub-id-type="doi">10.1029/2008GL036347</pub-id>
</element-citation>
</ref>
<ref id="CR14"><label>14.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Todesco</surname>
<given-names>M</given-names>
</name>
<name><surname>Rinaldi</surname>
<given-names>AP</given-names>
</name>
<name><surname>Bonafede</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Modeling of unrest signals in heterogeneous hydrothermal systems</article-title>
<source>J. Geophys. Res.</source>
<year>2010</year>
<volume>115</volume>
<fpage>B09213</fpage>
<pub-id pub-id-type="doi">10.1029/2010JB007474</pub-id>
</element-citation>
</ref>
<ref id="CR15"><label>15.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chiodini</surname>
<given-names>G</given-names>
</name>
<etal></etal>
</person-group>
<article-title>CO2 degassing and energy release at Solfatara volcano, Phlegraean Fields, Italy</article-title>
<source>J. Geophys. Res.</source>
<year>2001</year>
<volume>106</volume>
<fpage>16213</fpage>
<lpage>16221</lpage>
<pub-id pub-id-type="doi">10.1029/2001JB000246</pub-id>
</element-citation>
</ref>
<ref id="CR16"><label>16.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Petrillo</surname>
<given-names>Z</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Defining a 3D physical model for the hydrothermal circulation at Phlegraean Fields caldera (Italy)</article-title>
<source>J. Volcanol. Geotherm. Res.</source>
<year>2013</year>
<volume>274</volume>
<fpage>172</fpage>
<lpage>182</lpage>
<pub-id pub-id-type="doi">10.1016/j.jvolgeores.2013.08.008</pub-id>
</element-citation>
</ref>
<ref id="CR17"><label>17.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Coco</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Numerical models for ground deformation and gravity changes during volcanic unrest: simulating the hydrothermal system dynamics of a restless caldera</article-title>
<source>Solid Earth.</source>
<year>2016</year>
<volume>7</volume>
<fpage>557</fpage>
<lpage>577</lpage>
<pub-id pub-id-type="doi">10.5194/se-7-557-2016</pub-id>
</element-citation>
</ref>
<ref id="CR18"><label>18.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Afanasyev</surname>
<given-names>A</given-names>
</name>
<name><surname>Costa</surname>
<given-names>A</given-names>
</name>
<name><surname>Chiodini</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Investigation of hydrothermal activity at Campi Flegrei caldera using 3D numerical simulations: Extension to high temperature processes</article-title>
<source>Journal of Volcanology and Geothermal Research.</source>
<year>2015</year>
<volume>299</volume>
<fpage>68</fpage>
<lpage>77</lpage>
<pub-id pub-id-type="doi">10.1016/j.jvolgeores.2015.04.004</pub-id>
</element-citation>
</ref>
<ref id="CR19"><label>19.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Serra</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A strongly heterogeneous hydrothermal area imaged by surface waves: the case of Solfatara, Phlegraean Fields, Italy</article-title>
<source>Geophysical Journal International.</source>
<year>2016</year>
<volume>205</volume>
<fpage>1813</fpage>
<lpage>1822</lpage>
<pub-id pub-id-type="doi">10.1093/gji/ggw119</pub-id>
</element-citation>
</ref>
<ref id="CR20"><label>20.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Di Stefano</surname>
<given-names>R</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Fault zone properties affecting the rupture evolution of the 2009 (Mw 6.1) L’Aquila earthquake (central Italy): Insights from seismic tomography</article-title>
<source>Geophys. Res. Lett.</source>
<year>2011</year>
<volume>38</volume>
<fpage>L10310</fpage>
<pub-id pub-id-type="doi">10.1029/2011GL047365</pub-id>
</element-citation>
</ref>
<ref id="CR21"><label>21.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Amoroso</surname>
<given-names>O</given-names>
</name>
<name><surname>Ascione</surname>
<given-names>A</given-names>
</name>
<name><surname>Mazzoli</surname>
<given-names>S</given-names>
</name>
<name><surname>Virieux</surname>
<given-names>J</given-names>
</name>
<name><surname>Zollo</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Seismic imaging of a fluid storage in the actively extending Apennine mountain belt, southern Italy</article-title>
<source>Geophys. Res. Lett.</source>
<year>2014</year>
<volume>41</volume>
<fpage>3802</fpage>
<lpage>3809</lpage>
<pub-id pub-id-type="doi">10.1002/2014GL060070</pub-id>
</element-citation>
</ref>
<ref id="CR22"><label>22.</label>
<mixed-citation publication-type="other">Amoroso <italic>et al</italic>
. From Velocity and Attenuation Tomographies to Rock Physical Modeling: Inferences on fluid-driven earthquake processes at the Irpinia fault system in Southern Italy. <italic>Geophys. Res. Lett</italic>
. <bold>44</bold>
, doi:10.1002/2016GL072346 (2017).</mixed-citation>
</ref>
<ref id="CR23"><label>23.</label>
<mixed-citation publication-type="other">De Landro <italic>et al</italic>
. High resolution 3D imaging of the Irpinia active fault zone part 2: Modeling of rock physical properties. <italic>European Seismological Commission 35th General Assembly. ESC</italic>
2016-317-1, Vienna (2016).</mixed-citation>
</ref>
<ref id="CR24"><label>24.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zucca</surname>
<given-names>JJ</given-names>
</name>
<name><surname>Hutchings</surname>
<given-names>LJ</given-names>
</name>
<name><surname>Kasameyer</surname>
<given-names>PW</given-names>
</name>
</person-group>
<article-title>Seismic velocity and attenuation structure of the Geysers geothermal field, California</article-title>
<source>Geothermics.</source>
<year>1994</year>
<volume>23</volume>
<fpage>111</fpage>
<lpage>126</lpage>
<pub-id pub-id-type="doi">10.1016/0375-6505(94)90033-7</pub-id>
</element-citation>
</ref>
<ref id="CR25"><label>25.</label>
<mixed-citation publication-type="other">Gunasekera, R. C., Foulger, G. R. & Julian, B. R. Reservoir depletion at The Geysers geothermal area, California, shown by four-dimensional seismic tomography. <italic>Journal of Geophysical Research: Solid Eart</italic>
. <bold>108, </bold>
doi:10.1029/2001JB000638 (2003).</mixed-citation>
</ref>
<ref id="CR26"><label>26.</label>
<mixed-citation publication-type="other">Amoroso, O., De Landro, G., Russo, G. & Zollo, A. 4D imaging of elastic/anelastic medium properties: application to the Geysers geothermal area. <italic>AGIS Schatzalp Workshop on Induced Seismicity</italic>
, Davos Schatzalp, Switzerland (2015).</mixed-citation>
</ref>
<ref id="CR27"><label>27.</label>
<mixed-citation publication-type="other">Pasquet, S., Holbrook, W. S., Carr, B. J. & Sims, K. W. W. Geophysical imaging of shallow degassing in a Yellowstone hydrothermal system. <italic>Geophysical Research Letters</italic>
. <bold>43</bold>
(23), doi:10.1002/2016GL071306 (2016).</mixed-citation>
</ref>
<ref id="CR28"><label>28.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Husen</surname>
<given-names>S</given-names>
</name>
<name><surname>Smith</surname>
<given-names>RB</given-names>
</name>
<name><surname>Waite</surname>
<given-names>GP</given-names>
</name>
</person-group>
<article-title>Evidence for gas and magmatic sources beneath the Yellowstone volcanic field from seismic tomographic imaging</article-title>
<source>Journal of Volcanology and Geothermal Research</source>
<year>2004</year>
<volume>131</volume>
<fpage>397</fpage>
<lpage>410</lpage>
<pub-id pub-id-type="doi">10.1016/S0377-0273(03)00416-5</pub-id>
</element-citation>
</ref>
<ref id="CR29"><label>29.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vanorio</surname>
<given-names>T</given-names>
</name>
<name><surname>Virieux</surname>
<given-names>J</given-names>
</name>
<name><surname>Capuano</surname>
<given-names>P</given-names>
</name>
<name><surname>Russo</surname>
<given-names>G</given-names>
</name>
</person-group>
<article-title>Three-dimensional seismic tomography from P wave and S wave microearthquake travel times and rock physics characterization of the Phlegraean Fields Caldera</article-title>
<source>J. Geophys. Res.</source>
<year>2005</year>
<volume>110</volume>
<fpage>B03201</fpage>
<pub-id pub-id-type="doi">10.1029/2004JB003102</pub-id>
</element-citation>
</ref>
<ref id="CR30"><label>30.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>de Lorenzo</surname>
<given-names>S</given-names>
</name>
<name><surname>Gasparini</surname>
<given-names>P</given-names>
</name>
<name><surname>Mongelli</surname>
<given-names>F</given-names>
</name>
<name><surname>Zollo</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Thermal state of the Phlegraean Fields caldera inferred from seismic attenuation tomography</article-title>
<source>Journal of geodynamics.</source>
<year>2001</year>
<volume>32</volume>
<issue>4</issue>
<fpage>467</fpage>
<lpage>486</lpage>
<pub-id pub-id-type="doi">10.1016/S0264-3707(01)00044-8</pub-id>
</element-citation>
</ref>
<ref id="CR31"><label>31.</label>
<mixed-citation publication-type="other">De Siena, L., Del Pezzo, E. & Bianco, F. Seismic attenuation imaging of Phlegraean Fields: Evidence of gas reservoirs, hydrothermal basins, and feeding systems. <italic>Journal of Geophysical Research: Solid Earth</italic>
. <bold>115</bold>
(B9), doi:10.1029/2009JB006938 (2010).</mixed-citation>
</ref>
<ref id="CR32"><label>32.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bruno</surname>
<given-names>PP</given-names>
</name>
</person-group>
<article-title>Geophysical and hydrogeological experiments from a shallow hydrothermal system at Solfatara Volcano, Phlegraean Fields, Italy: Response to caldera unrest</article-title>
<source>J. Geophys. Res.</source>
<year>2007</year>
<volume>112</volume>
<fpage>B06201</fpage>
<pub-id pub-id-type="doi">10.1029/2006JB004383</pub-id>
</element-citation>
</ref>
<ref id="CR33"><label>33.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Isaia</surname>
<given-names>R</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Stratigraphy, structure, and volcano-tectonic evolution of Solfatara maar-diatreme (Phlegraean Fields, Italy)</article-title>
<source>Geological Society of America Bulletin</source>
<year>2015</year>
<volume>127</volume>
<fpage>1485</fpage>
<lpage>1504</lpage>
<pub-id pub-id-type="doi">10.1130/B31183.1</pub-id>
</element-citation>
</ref>
<ref id="CR34"><label>34.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Harris</surname>
<given-names>AJL</given-names>
</name>
<name><surname>Carniel</surname>
<given-names>R</given-names>
</name>
<name><surname>Jones</surname>
<given-names>J</given-names>
</name>
</person-group>
<article-title>Identification of variable convectives at Erta’ Ale Lava Lake</article-title>
<source>J. Volcanol. Geotherm. Res.</source>
<year>2005</year>
<volume>142</volume>
<fpage>207</fpage>
<lpage>223</lpage>
<pub-id pub-id-type="doi">10.1016/j.jvolgeores.2004.11.011</pub-id>
</element-citation>
</ref>
<ref id="CR35"><label>35.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Caliro</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Recent activity of Nisyros volcano (Greece) inferred from structural, geochemical and seismological data</article-title>
<source>Bull. Volcanol.</source>
<year>2005</year>
<volume>67</volume>
<fpage>358</fpage>
<lpage>369</lpage>
<pub-id pub-id-type="doi">10.1007/s00445-004-0381-7</pub-id>
</element-citation>
</ref>
<ref id="CR36"><label>36.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aiuppa</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Patterns in the recent 2007–2008 activity of Mount Etna volcano investigated by integrated geophysical and geochemical observations</article-title>
<source>Geochem. Geophys. Geosyst.</source>
<year>2010</year>
<volume>11</volume>
<fpage>Q09008</fpage>
<pub-id pub-id-type="doi">10.1029/2010GC003168</pub-id>
</element-citation>
</ref>
<ref id="CR37"><label>37.</label>
<mixed-citation publication-type="other">Byrdina, S. <italic>et al</italic>
. Relations between electrical resistivity, carbon dioxide flux, and self-potential in the shallow hydrothermal system of Solfatara (Phlegrean Fields, Italy). <italic>Journal of Volcanology and Geothermal Research</italic>
, ISSN 0377–0273, <ext-link ext-link-type="uri" xlink:href="https://doi.org/10.1016/j.jvolgeores.2014.07.010">https://doi.org/10.1016/j.jvolgeores.2014.07.010</ext-link>
(2014).</mixed-citation>
</ref>
<ref id="CR38"><label>38.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aster</surname>
<given-names>RC</given-names>
</name>
<name><surname>Meyer</surname>
<given-names>RP</given-names>
</name>
</person-group>
<article-title>Three-dimensional velocity structure and hypocenter distribution in the Phlegraean Fields caldera, Italy</article-title>
<source>Tectonophysics.</source>
<year>1988</year>
<volume>149</volume>
<fpage>195</fpage>
<lpage>218</lpage>
<pub-id pub-id-type="doi">10.1016/0040-1951(88)90173-4</pub-id>
</element-citation>
</ref>
<ref id="CR39"><label>39.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tramelli</surname>
<given-names>A</given-names>
</name>
<name><surname>Del Pezzo</surname>
<given-names>E</given-names>
</name>
<name><surname>Bianco</surname>
<given-names>F</given-names>
</name>
<name><surname>Boschi</surname>
<given-names>E</given-names>
</name>
</person-group>
<article-title>3D scattering image of the Phlegraean Fields caldera (Southern Italy): new hints on the position of the old caldera rim</article-title>
<source>Physics of the Earth and Planetary Interiors.</source>
<year>2006</year>
<volume>155</volume>
<fpage>269</fpage>
<lpage>280</lpage>
<pub-id pub-id-type="doi">10.1016/j.pepi.2005.12.009</pub-id>
</element-citation>
</ref>
<ref id="CR40"><label>40.</label>
<mixed-citation publication-type="other">Serlenga, V. <italic>et al</italic>
. A three-dimensional QP imaging of the shallowest subsurface of Phlegraean Fields offshore caldera, southern Italy. <italic>Geophysical Research Letters</italic>
. <bold>43</bold>
(21), doi:10.1002/2016GL071140 (2016).</mixed-citation>
</ref>
<ref id="CR41"><label>41.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Caudill</surname>
<given-names>M</given-names>
</name>
</person-group>
<article-title>Neural networks primer, part VII</article-title>
<source>AI Expert 4.5.</source>
<year>1989</year>
<volume>5</volume>
<fpage>51</fpage>
<lpage>58</lpage>
</element-citation>
</ref>
<ref id="CR42"><label>42.</label>
<mixed-citation publication-type="other">Fahlman, S. & Lebiere, C. The Cascade Correlation Learning Architecture in <italic>Advances in Neural Information Processing</italic>
Systems <italic>2</italic>
(eds. Touretzky D. S.) 524–532 (1990).</mixed-citation>
</ref>
<ref id="CR43"><label>43.</label>
<mixed-citation publication-type="other">Rumelhart, D. E., McClelland, J. L. and PDP Research Group. In <italic>Parallel Distributed Processing: Explorations in the Microstructure of Cognition</italic>
, Vol. 1, 12–21 (MIT Press, Cambridge MA, 1986).</mixed-citation>
</ref>
<ref id="CR44"><label>44.</label>
<mixed-citation publication-type="other">Wasserman, P. D. In <italic>Neural Computing Theory and Practice</italic>
. 44–54 (Van Nostrand Reinhold, New York, 1989).</mixed-citation>
</ref>
<ref id="CR45"><label>45.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Latorre</surname>
<given-names>D</given-names>
</name>
<etal></etal>
</person-group>
<article-title>A new seismic tomography of Aigion area (Gulf of Corinth, Greece) from the 1991 data set</article-title>
<source>Geophys. J. Int.</source>
<year>2004</year>
<volume>159</volume>
<fpage>1013</fpage>
<lpage>1031</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-246X.2004.02412.x</pub-id>
</element-citation>
</ref>
<ref id="CR46"><label>46.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zollo</surname>
<given-names>A</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Bayesian estimation of 2-D P-velocity models from active seismic arrival time data: imaging of the shallow structure of Mt Vesuvius (Southern Italy)</article-title>
<source>Geophysical Journal International.</source>
<year>2002</year>
<volume>151.2</volume>
<fpage>566</fpage>
<lpage>582</lpage>
<pub-id pub-id-type="doi">10.1046/j.1365-246X.2002.01795.x</pub-id>
</element-citation>
</ref>
<ref id="CR47"><label>47.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lutter</surname>
<given-names>WJ</given-names>
</name>
<name><surname>Nowack</surname>
<given-names>R</given-names>
</name>
<name><surname>Braile</surname>
<given-names>LW</given-names>
</name>
</person-group>
<article-title>Inversion for crustal structure using reflections from the PASSCAL Ouachita experiment</article-title>
<source>J. geophys. Res.</source>
<year>1990</year>
<volume>95</volume>
<fpage>4623</fpage>
<lpage>4646</lpage>
</element-citation>
</ref>
<ref id="CR48"><label>48.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Akaike</surname>
<given-names>H</given-names>
</name>
</person-group>
<article-title>A new look at the statistical model identification</article-title>
<source>IEEE Transactions on Automatic Control.</source>
<year>1974</year>
<volume>19</volume>
<fpage>716</fpage>
<lpage>723</lpage>
<pub-id pub-id-type="doi">10.1109/TAC.1974.1100705</pub-id>
</element-citation>
</ref>
<ref id="CR49"><label>49.</label>
<mixed-citation publication-type="other">Cavanaugh, J. E. & Shumway, R. H. A bootstrap variant of AIC for state-space model selection. <italic>StatisticaSinica</italic>
, 473–496 (1997).</mixed-citation>
</ref>
<ref id="CR50"><label>50.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Eberhart-Phillips</surname>
<given-names>D</given-names>
</name>
</person-group>
<article-title>Three-dimensional velocity structure in northern California Coast Ranges from inversion of local earthquake arrival times</article-title>
<source>Bulletin of the Seismological Society of America</source>
<year>1986</year>
<volume>76</volume>
<fpage>1025</fpage>
<lpage>1052</lpage>
</element-citation>
</ref>
<ref id="CR51"><label>51.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Michelini</surname>
<given-names>A</given-names>
</name>
<name><surname>McEvilly</surname>
<given-names>TV</given-names>
</name>
</person-group>
<article-title>Seismological studies at Parkfield. I. Simultaneous inversion for velocity structure and hypocenters using cubic B-splines parameterization</article-title>
<source>Bulletin of the Seismological Society of America</source>
<year>1991</year>
<volume>81</volume>
<fpage>524</fpage>
<lpage>552</lpage>
</element-citation>
</ref>
<ref id="CR52"><label>52.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Caliro</surname>
<given-names>S</given-names>
</name>
<etal></etal>
</person-group>
<article-title>The origin of the fumaroles of La Solfatara (Phlegraean Fields, South Italy)</article-title>
<source>Geochimica et CosmochimicaActa.</source>
<year>2007</year>
<volume>71</volume>
<fpage>3040</fpage>
<lpage>3055</lpage>
<pub-id pub-id-type="doi">10.1016/j.gca.2007.04.007</pub-id>
</element-citation>
</ref>
<ref id="CR53"><label>53.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sissons</surname>
<given-names>BA</given-names>
</name>
<name><surname>Dibble</surname>
<given-names>RR</given-names>
</name>
</person-group>
<article-title>A seismic refraction experiment southeast of Ruapehu volcano</article-title>
<source>New Zealand Journal of Geology and Geophysics</source>
<year>1981</year>
<volume>24</volume>
<fpage>31</fpage>
<lpage>38</lpage>
<pub-id pub-id-type="doi">10.1080/00288306.1981.10422696</pub-id>
</element-citation>
</ref>
<ref id="CR54"><label>54.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mora</surname>
<given-names>M</given-names>
</name>
<etal></etal>
</person-group>
<article-title>Shallow velocity structure and seismic site effects at Arenal volcano, Costa Rica</article-title>
<source>Journal of Volcanology and Geothermal Research</source>
<year>2006</year>
<volume>152</volume>
<fpage>121</fpage>
<lpage>139</lpage>
<pub-id pub-id-type="doi">10.1016/j.jvolgeores.2005.09.013</pub-id>
</element-citation>
</ref>
<ref id="CR55"><label>55.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Toksöz</surname>
<given-names>NM</given-names>
</name>
<name><surname>Cheng</surname>
<given-names>CH</given-names>
</name>
<name><surname>Timur</surname>
<given-names>A</given-names>
</name>
</person-group>
<article-title>Velocities of seismic waves in porous rocks</article-title>
<source>Geophysics.</source>
<year>1976</year>
<volume>41</volume>
<fpage>621</fpage>
<lpage>645</lpage>
<pub-id pub-id-type="doi">10.1190/1.1440639</pub-id>
</element-citation>
</ref>
<ref id="CR56"><label>56.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Loke</surname>
<given-names>MH</given-names>
</name>
<name><surname>Dahlin</surname>
<given-names>T</given-names>
</name>
</person-group>
<article-title>A comparison of Gauss-Newton and quasi-Newton methods in resistivity imaging inversion</article-title>
<source>J. app. Geophys.</source>
<year>2002</year>
<volume>49</volume>
<fpage>149</fpage>
<lpage>162</lpage>
<pub-id pub-id-type="doi">10.1016/S0926-9851(01)00106-9</pub-id>
</element-citation>
</ref>
<ref id="CR57"><label>57.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Loke</surname>
<given-names>MH</given-names>
</name>
<name><surname>Barker</surname>
<given-names>RD</given-names>
</name>
</person-group>
<article-title>Practical techniques for 3D resistivity surveys and data inversion</article-title>
<source>Geophysical Prospecting.</source>
<year>1996</year>
<volume>44</volume>
<fpage>499</fpage>
<lpage>523</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-2478.1996.tb00162.x</pub-id>
</element-citation>
</ref>
<ref id="CR58"><label>58.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Revil</surname>
<given-names>A</given-names>
</name>
<name><surname>Florsch</surname>
<given-names>N</given-names>
</name>
</person-group>
<article-title>Determination of permeability from spectral induced polarization in granular media</article-title>
<source>Geophys. J. Int.</source>
<year>2010</year>
<volume>181</volume>
<issue>3</issue>
<fpage>1480</fpage>
<lpage>1498</lpage>
</element-citation>
</ref>
<ref id="CR59"><label>59.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chiodini</surname>
<given-names>G</given-names>
</name>
<name><surname>Frondini</surname>
<given-names>F</given-names>
</name>
<name><surname>Raco</surname>
<given-names>B</given-names>
</name>
</person-group>
<article-title>Diffuse emission of CO2 from the Fossa crater, Vulcano Island (Italy)</article-title>
<source>Bulletin of Volcanology</source>
<year>1996</year>
<volume>58</volume>
<fpage>41</fpage>
<lpage>50</lpage>
<pub-id pub-id-type="doi">10.1007/s004450050124</pub-id>
</element-citation>
</ref>
<ref id="CR60"><label>60.</label>
<element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Letort</surname>
<given-names>J</given-names>
</name>
<etal></etal>
</person-group>
<article-title>High-resolution shallow seismic tomography of a hydrothermal area: application to the Solfatara, Pozzuoli</article-title>
<source>Geophys. J. Int.</source>
<year>2012</year>
<volume>189</volume>
<fpage>1725</fpage>
<lpage>1733</lpage>
<pub-id pub-id-type="doi">10.1111/j.1365-246X.2012.05451.x</pub-id>
</element-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Terre/explor/NissirosV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000005 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000005 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Wicri/Terre |area= NissirosV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:5469761 |texte= 3D ultra-high resolution seismic imaging of shallow Solfatara crater in Campi Flegrei (Italy): New insights on deep hydrothermal fluid circulation processes }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:28611382" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a NissirosV1
![]() | This area was generated with Dilib version V0.6.33. | ![]() |