Salamat sa pagbisita sa Nature.com. Limitado ang suporta para sa CSS sa bersyon ng browser na iyong ginagamit. Para sa pinakamahusay na karanasan, inirerekomenda namin na gumamit ka ng na-update na browser (o i-off ang compatibility mode sa Internet Explorer). Samantala, upang matiyak ang patuloy na suporta, ipapakita namin ang site nang walang mga estilo at JavaScript.
Iniuulat namin ang ebidensya ng aktibong pagtaas ng sahig ng dagat at mga emisyon ng gas ilang kilometro mula sa pampang mula sa daungan ng Naples (Italya). Ang mga pockmark, tambak, at bunganga ay mga katangian ng sahig ng dagat. Ang mga pormasyong ito ay kumakatawan sa mga tuktok ng mababaw na istruktura ng crustal, kabilang ang mga pagoda, fault, at fold na nakakaapekto sa sahig ng dagat ngayon. Naitala nila ang pagtaas, pressurization, at paglabas ng helium at carbon dioxide sa mga reaksyon ng decarbonization ng mga natutunaw na mantle at mga batong crustal. Ang mga gas na ito ay malamang na katulad ng mga nagpapakain sa mga hydrothermal system ng Ischia, Campi Flegre, at Soma-Vesuvius, na nagmumungkahi ng isang pinagmumulan ng mantle na may halong crustal fluid sa ilalim ng Golpo ng Naples. Ang paglawak at pagkapunit sa ilalim ng dagat na dulot ng proseso ng pagtaas at pressurization ng gas ay nangangailangan ng overpressure na 2-3 MPa. Ang mga pagtaas ng sahig ng dagat, mga fault, at mga emisyon ng gas ay mga manipestasyon ng mga pag-aalsa na hindi bulkan na maaaring magpahiwatig ng mga pagsabog ng sahig ng dagat at/o mga pagsabog ng hydrothermal.
Ang mga paglabas ng hydrothermal (mainit na tubig at gas) sa malalim na dagat ay isang karaniwang katangian ng mga tagaytay sa kalagitnaan ng karagatan at mga gilid ng convergent plate (kabilang ang mga nakalubog na bahagi ng mga arko ng isla), samantalang ang mga malamig na paglabas ng mga gas hydrate (chlatrate) ay kadalasang katangian ng mga continental shelf at passive margin1, 2,3,4,5. Ang paglitaw ng mga paglabas ng hydrothermal sa sahig ng dagat sa mga lugar sa baybayin ay nagpapahiwatig ng mga pinagmumulan ng init (mga imbakan ng magma) sa loob ng continental crust at/o mantle. Ang mga paglabas na ito ay maaaring mauna sa pag-akyat ng magma sa pinakamataas na mga patong ng crust ng Daigdig at magtatapos sa pagsabog at paglalagay ng mga bulkan sa dagat6. Samakatuwid, ang pagkilala sa (a) mga morpolohiya na nauugnay sa aktibong deformasyon ng seabed at (b) mga emisyon ng gas na malapit sa mga mataong lugar sa baybayin tulad ng rehiyon ng bulkan ng Naples sa Italya (~1 milyong naninirahan) ay mahalaga para sa pagtatasa ng mga posibleng bulkan. Mababaw na pagsabog. Bukod pa rito, habang ang mga katangiang morpolohikal na nauugnay sa mga emisyon ng hydrothermal o hydrate sa malalim na dagat ay medyo kilala dahil sa kanilang mga geological at biological na katangian, ang mga eksepsiyon ay mga katangiang morpolohikal na nauugnay sa mas mababaw na tubig, maliban sa mga nangyayari sa In Lake 12, medyo kakaunti ang mga tala. Dito, ipinapakita namin ang mga bagong datos mula sa bathymetric, seismic, water column, at geochemical para sa isang rehiyon sa ilalim ng tubig, morpolohikal at estruktural na masalimuot na apektado ng mga emisyon ng gas sa Golpo ng Naples (Timog Italya), humigit-kumulang 5 km mula sa daungan ng Naples. Ang mga datos na ito ay nakolekta noong panahon ng SAFE_2014 (Agosto 2014) na paglalakbay sakay ng R/V Urania. Inilalarawan at binibigyang-kahulugan namin ang mga istruktura sa ilalim ng dagat at ilalim ng lupa kung saan nagaganap ang mga emisyon ng gas, sinisiyasat ang mga pinagmumulan ng mga likidong naglalabas ng hangin, tinutukoy at kinikilala ang mga mekanismo na kumokontrol sa pagtaas ng gas at kaugnay na deformasyon, at tinatalakay ang mga epekto ng bulkanolohiya.
Ang Golpo ng Naples ay bumubuo sa kanlurang gilid na Plio-Quaternary, ang NW-SE ay pahabang Campania tectonic depression13,14,15.EW ng Ischia (mga 150-1302 AD), bunganga ng Campi Flegre (mga 300-1538) at Soma-Vesuvius (mula <360-1944). Ang pagkakaayos ay naglilimita sa look sa hilaga AD)15, habang ang timog ay nasa hangganan ng Sorrento Peninsula (Fig. 1a). Ang Golpo ng Naples ay apektado ng umiiral na NE-SW at pangalawang NW-SE na mahahalagang fault (Fig. 1)14,15. Ang Ischia, Campi Flegre at Somma-Vesuvius ay nailalarawan sa pamamagitan ng mga hydrothermal manifestation, ground deformation, at shallow seismicity16,17,18 (hal., ang magulong kaganapan sa Campi Flegrei noong 1982-1984, na may pagtaas ng 1.8 m at libu-libong lindol). Ang mga kamakailang pag-aaral19,20 ay nagmumungkahi na Maaaring may kaugnayan sa dinamika ng Soma-Vesuvius at ng Campi Flegre, na posibleng nauugnay sa 'malalim' na iisang reservoir ng magma. Kinontrol ng aktibidad ng bulkan at mga osilasyon sa antas ng dagat sa huling 36 ka ng Campi Flegrei at 18 ka ng Somma Vesuvius ang sedimentary system ng Golpo ng Naples. Ang mababang antas ng dagat sa huling glacial maximum (18 ka) ay humantong sa regresyon ng offshore-shallow sedimentary system, na kalaunan ay napunan ng mga transgressive na kaganapan noong Late Pleistocene-Holocene. Natukoy ang mga emisyon ng submarine gas sa paligid ng isla ng Ischia at sa baybayin ng Campi Flegre at malapit sa Mount Soma-Vesuvius (Fig. 1b).
(a) Morpolohikal at istruktural na kaayusan ng continental shelf at ng Golpo ng Naples 15, 23, 24, 48. Ang mga tuldok ay mga pangunahing sentro ng pagsabog ng submarino; ang mga pulang linya ay kumakatawan sa mga pangunahing fault.(b) Bathymetry ng Look ng Naples na may mga nakitang fluid vent (mga tuldok) at mga bakas ng mga seismic lines (mga itim na linya). Ang mga dilaw na linya ay ang mga trajectory ng mga seismic lines na L1 at L2 na iniulat sa Figure 6. Ang mga hangganan ng mga istrukturang parang simboryo ng Banco della Montagna (BdM) ay minarkahan ng mga asul na putol-putol na linya sa (a,b). Ang mga dilaw na parisukat ay minarkahan ang mga lokasyon ng mga profile ng acoustic water column, at ang mga CTD-EMBlank, CTD-EM50 at ROV frame ay iniulat sa Fig. 5. Ang dilaw na bilog ay minarkahan ang lokasyon ng sampling gas discharge, at ang komposisyon nito ay ipinapakita sa Table S1. Ang Golden Software (http://www.goldensoftware.com/products/surfer) ay gumagamit ng mga graphics na nabuo ng Surfer® 13.
Batay sa datos na nakuha noong SAFE_2014 (Agosto 2014) cruise (tingnan ang Mga Paraan), isang bagong Digital Terrain Model (DTM) ng Golpo ng Naples na may 1 m na resolusyon ang naitayo. Ipinapakita ng DTM na ang seafloor sa timog ng Port of Naples ay nailalarawan sa pamamagitan ng isang bahagyang pahilig na nakaharap sa timog (slope ≤3°) na ibabaw na naputol ng isang 5.0 × 5.3 km na parang simboryo na istraktura, na lokal na kilala bilang Banco della Montagna (BdM). Fig. 1a,b). Ang BdM ay nabubuo sa lalim na humigit-kumulang 100 hanggang 170 metro, 15 hanggang 20 metro sa ibabaw ng nakapalibot na sahig ng dagat. Ang simboryo ng BdM ay nagpakita ng mala-bundok na morpolohiya dahil sa 280 subcircular hanggang oval na mga burol (Larawan 2a), 665 cone, at 30 hukay (Mga Larawan 3 at 4). Ang burol ay may pinakamataas na taas at circumference na 22 m at 1,800 m, ayon sa pagkakabanggit. Ang circularity [C = 4π(area/perimeter2)] ng mga burol ay bumaba kasabay ng pagtaas ng perimeter (Larawan 2b). Ang mga axial ratio para sa mga burol ay nasa pagitan ng 1 at 6.5, kung saan ang mga burol na may axial ratio >2 ay nagpapakita ng mas gustong N45°E + 15° strike at mas nakakalat na sekondarya, mas nakakalat na N105°E hanggang N145°E strike (Larawan 2c). May mga single o nakahanay na cone sa BdM plane at sa ibabaw ng mound (Fig. 3a,b). Ang mga conical na kaayusan ay sumusunod sa kaayusan ng mga mound kung saan sila matatagpuan. Ang mga pockmark ay karaniwang matatagpuan sa patag na seabed (Fig. 3c) at paminsan-minsan sa mga mound. Ang spatial densities ng mga cone at pockmark ay nagpapakita na ang nangingibabaw na NE-SW alignment ay naglilimita sa mga hangganan ng hilagang-silangan at timog-kanluran ng BdM dome (Fig. 4a,b); ang hindi gaanong pinahabang NW-SE na ruta ay matatagpuan sa gitnang rehiyon ng BdM.
(a) Digital terrain model (1 m cell size) ng simboryo ng Banco della Montagna (BdM).(b) Perimeter at pagiging bilog ng mga tambak ng BdM.(c) Axial ratio at anggulo (orientasyon) ng pangunahing axis ng pinakamahusay na akmang ellipse na nakapalibot sa tambak. Ang standard error ng Digital Terrain model ay 0.004 m; ang standard errors ng perimeter at pagiging bilog ay 4.83 m at 0.01, ayon sa pagkakabanggit, at ang standard errors ng axial ratio at anggulo ay 0.04 at 3.34°, ayon sa pagkakabanggit.
Ang mga detalye ng mga natukoy na kono, bunganga, tambak, at hukay sa rehiyon ng BdM na kinuha mula sa DTM ay nasa Figure 2.
(a) Mga kono ng pagkakahanay sa patag na seabed; (b) mga kono at bunganga sa mga payat na tambak mula NW-SE; (c) mga bakas ng butas sa bahagyang nakalubog na ibabaw.
(a) Distribusyon sa espasyo ng mga natukoy na bunganga, hukay, at aktibong paglabas ng gas. (b) Densidad sa espasyo ng mga bunganga at hukay na iniulat sa (a) (bilang/0.2 km2).
Natukoy namin ang 37 gaseous emissions sa rehiyon ng BdM mula sa mga imahe ng ROV water column echo sounder at direktang obserbasyon sa seafloor na nakuha noong SAFE_2014 cruise noong Agosto 2014 (Mga Larawan 4 at 5). Ang mga acoustic anomalies ng mga emissions na ito ay nagpapakita ng mga patayong pahabang hugis na umaangat mula sa seafloor, na may taas na patayo sa pagitan ng 12 at humigit-kumulang 70 m (Larawan 5a). Sa ilang mga lugar, ang mga acoustic anomalies ay bumuo ng halos tuluy-tuloy na "tren." Ang naobserbahang mga bubble plume ay lubhang nag-iiba: mula sa tuluy-tuloy, siksik na daloy ng bubble hanggang sa panandaliang phenomena (Karagdagang Pelikula 1). Ang inspeksyon ng ROV ay nagbibigay-daan para sa visual na pag-verify ng paglitaw ng mga seafloor fluid vent at nagtatampok ng maliliit na pockmark sa seabed, minsan ay napapalibutan ng pula hanggang orange na sediments (Larawan 5b). Sa ilang mga kaso, ang mga ROV channel ay muling nagpapagana ng mga emisyon. Ang morphology ng vent ay nagpapakita ng isang pabilog na butas sa itaas na walang flare sa water column. Ang pH sa water column sa itaas lamang ng discharge point ay nagpakita ng isang makabuluhang pagbaba, na nagpapahiwatig ng mas acidic na mga kondisyon sa lokal (Larawan 5c, d). Sa partikular, ang pH sa itaas ng BdM gas Ang paglabas sa lalim na 75 m ay bumaba mula 8.4 (sa lalim na 70 m) patungong 7.8 (sa lalim na 75 m) (Larawan 5c), samantalang ang ibang mga lugar sa Golpo ng Naples ay may mga halaga ng pH sa pagitan ng 0 at 160 m sa pagitan ng lalim na 8.3 at 8.5 (Larawan 5d). Kulang ang mga makabuluhang pagbabago sa temperatura at kaasinan ng tubig-dagat sa dalawang lugar sa loob at labas ng lugar ng BdM ng Golpo ng Naples. Sa lalim na 70 m, ang temperatura ay 15 °C at ang kaasinan ay humigit-kumulang 38 PSU (Larawan 5c,d). Ang mga sukat ng pH, temperatura, at kaasinan ay nagpapahiwatig ng: a) ang pakikilahok ng mga acidic fluid na nauugnay sa proseso ng pag-alis ng gas sa BdM at b) ang kawalan o napakabagal na paglabas ng mga thermal fluid at brine.
(a) Acquisition window ng acoustic water column profile (echometer Simrad EK60). Patayong berdeng banda na tumutugma sa gas flare na natukoy sa EM50 fluid discharge (mga 75 m sa ibaba ng antas ng dagat) na matatagpuan sa rehiyon ng BdM; ang mga signal sa ilalim at seafloor multiplex ay ipinapakita rin (b) na kinolekta gamit ang isang remote-controlled na sasakyan sa rehiyon ng BdM. Ang nag-iisang larawan ay nagpapakita ng isang maliit na bunganga (itim na bilog) na napapalibutan ng pula hanggang orange na sediment. (c,d) Multiparameter probe CTD data na pinoproseso gamit ang SBED-Win32 software (Seasave, bersyon 7.23.2). Mga pattern ng mga piling parameter (salinity, temperatura, pH at oxygen) ng water column sa itaas ng fluid discharge EM50 (panel c) at sa labas ng Bdm discharge area panel (d).
Nakakolekta kami ng tatlong sample ng gas mula sa lugar ng pag-aaral sa pagitan ng Agosto 22 at 28, 2014. Ang mga sample na ito ay nagpakita ng magkatulad na komposisyon, na pinangungunahan ng CO2 (934-945 mmol/mol), na sinundan ng mga kaugnay na konsentrasyon ng N2 (37-43 mmol/mol), CH4 (16-24 mmol/mol) at H2S (0.10 mmol/mol) -0.44 mmol/mol), habang ang H2 at He ay hindi gaanong sagana (<0.052 at <0.016 mmol/mol, ayon sa pagkakabanggit) (Fig. 1b; Table S1, Supplementary Movie 2). Sinukat din ang medyo mataas na konsentrasyon ng O2 at Ar (hanggang 3.2 at 0.18 mmol/mol, ayon sa pagkakabanggit). Ang kabuuan ng mga light hydrocarbon ay mula 0.24 hanggang 0.30 mmol/mol at binubuo ng mga C2-C4 alkane, aromatics (pangunahing benzene), propene at mga compound na naglalaman ng sulfur (thiophene). Ang halaga ng 40Ar/36Ar ay naaayon sa hangin. (295.5), bagama't ang sample na EM35 (BdM dome) ay may halagang 304, na nagpapakita ng bahagyang labis na 40Ar. Ang δ15N ratio ay mas mataas kaysa sa hangin (hanggang +1.98% vs. Hangin), habang ang mga halaga ng δ13C-CO2 ay mula -0.93 hanggang 0.44% vs. V-PDB. Ang mga halaga ng R/Ra (pagkatapos itama ang polusyon sa hangin gamit ang 4He/20Ne ratio) ay nasa pagitan ng 1.66 at 1.94, na nagpapahiwatig ng pagkakaroon ng isang malaking bahagi ng mantle He. Sa pamamagitan ng pagsasama-sama ng helium isotope sa CO2 at sa matatag nitong isotope 22, ang pinagmulan ng mga emisyon sa BdM ay maaaring higit pang linawin. Sa mapa ng CO2 para sa CO2/3He kumpara sa δ13C (Larawan 6), ang komposisyon ng gas na BdM ay inihambing sa mga fumarole ng Ischia, Campi Flegrei at Somma-Vesuvius. Iniuulat din ng Larawan 6 ang teoretikal na paghahalo. mga linya sa pagitan ng tatlong magkakaibang pinagmumulan ng carbon na maaaring kasangkot sa produksyon ng BdM gas: mga dissolved mantle-derived melts, mga organic-rich sediments, at mga carbonates. Ang mga sample ng BdM ay nasa mixing line na inilalarawan ng tatlong bulkan ng Campania, ibig sabihin, ang paghahalo sa pagitan ng mga mantle gas (na ipinapalagay na bahagyang mayaman sa carbon dioxide kumpara sa mga classical MORB para sa layunin ng pag-aakma sa datos) at mga reaksyon na dulot ng crustal decarbonization Ang nagresultang gas rock.
Ang mga hybrid na linya sa pagitan ng komposisyon ng mantle at mga dulong bahagi ng limestone at mga organikong sediment ay iniulat para sa paghahambing. Ang mga kahon ay kumakatawan sa mga fumarole area ng Ischia, Campi Flegrei at Somma-Vesvius 59, 60, 61. Ang BdM sample ay nasa mixed trend ng bulkang Campania. Ang endmember gas ng mixed line ay nagmula sa mantle, na siyang gas na nalilikha ng decarburization reaction ng mga carbonate mineral.
Ang mga seismic section na L1 at L2 (Figs. 1b at 7) ay nagpapakita ng transisyon sa pagitan ng BdM at ng mga distal stratigraphic sequence ng mga bulkanong rehiyon ng Somma-Vesuvius (L1, Fig. 7a) at Campi Flegrei (L2, Fig. 7b). Ang BdM ay nailalarawan sa pamamagitan ng pagkakaroon ng dalawang pangunahing seismic formations (MS at PS sa Fig. 7). Ang nasa itaas (MS) ay nagpapakita ng mga subparallel reflector na may mataas hanggang katamtamang amplitude at lateral continuity (Fig. 7b,c). Kasama sa layer na ito ang mga marine sediment na hinihila ng Last Glacial Maximum (LGM) system at binubuo ng buhangin at clay23. Ang pinagbabatayang PS layer (Fig. 7b–d) ay nailalarawan sa pamamagitan ng isang magulong hanggang transparent na phase sa hugis ng mga haligi o hourglass. Ang itaas ng mga PS sediment ay bumuo ng mga seafloor mound (Fig. 7d). Ang mga diapir-like geometries na ito ay nagpapakita ng pagpasok ng PS transparent na materyal sa pinakamataas na MS deposits. Ang Uplift ay responsable para sa pagbuo ng mga fold at fault na nakakaapekto sa MS layer at nakapatong sa kasalukuyang mga sediment ng BdM seafloor (Fig. 7b–d). Ang MS stratigraphic interval ay malinaw na nakahiwalay sa ENE na bahagi ng seksyong L1, habang ito ay pumuputi patungo sa BdM dahil sa presensya ng isang gas-saturated layer (GSL) na natatakpan ng ilang panloob na antas ng MS sequence (Fig. 7a). Ang mga gravity core na nakolekta sa tuktok ng BdM na katumbas ng transparent seismic layer ay nagpapahiwatig na ang pinakamataas na 40 cm ay binubuo ng buhangin na idineposito kamakailan hanggang sa kasalukuyan; )24,25 at mga piraso ng pumice mula sa pagsabog ng Campi Flegrei ng “Naples Yellow Tuff” (14.8 ka)26. Ang transparent na yugto ng PS layer ay hindi maipaliwanag sa pamamagitan lamang ng magulong proseso ng paghahalo, dahil ang magulong mga layer na nauugnay sa mga landslide, daloy ng putik at mga pyroclastic flow na matatagpuan sa labas ng BdM sa Golpo ng Naples ay acoustically opaque21,23,24. Napagpasyahan namin na ang naobserbahang BdM PS seismic facies pati na rin ang hitsura ng subsea outcrop PS layer (Fig. 7d) ay sumasalamin sa pagtaas ng natural gas.
(a) Single-track seismic profile L1 (navigation trace sa Fig. 1b) na nagpapakita ng isang columnar (pagoda) spatial arrangement. Ang pagoda ay binubuo ng magulong deposito ng pumice at buhangin. Ang gas-saturated layer na nasa ilalim ng pagoda ay nag-aalis ng continuity ng mas malalalim na formations. (b) Single-channel seismic profile L2 (navigation trace sa Fig. 1b), na nagtatampok ng incision at deformation ng mga seafloor mound, marine (MS), at pumice sand deposits (PS). (c) Ang mga detalye ng deformation sa MS at PS ay iniulat sa (c,d). Sa pag-aakalang ang bilis ay 1580 m/s sa pinakamataas na sediment, ang 100 ms ay kumakatawan sa humigit-kumulang 80 m sa vertical scale.
Ang mga katangiang morpolohikal at estruktural ng BdM ay katulad ng iba pang mga subsea hydrothermal at gas hydrate field sa buong mundo2,12,27,28,29,30,31,32,33,34 at kadalasang iniuugnay sa mga uplift (vault at mound) at gas Discharge (cones, pits). Ang mga cones at pits na nakahanay sa BdM at mga pahabang mound ay nagpapahiwatig ng structurally controlled permeability (Mga Larawan 2 at 3). Ang spatial arrangement ng mga mound, pits at aktibong vent ay nagmumungkahi na ang kanilang distribusyon ay bahagyang kinokontrol ng NW-SE at NE-SW impact fractures (Fig. 4b). Ito ang mga ginustong strike ng mga fault system na nakakaapekto sa mga bulkanic area ng Campi Flegrei at Somma-Vesuvius at sa Golpo ng Naples. Sa partikular, ang istruktura ng nauna ay kumokontrol sa lokasyon ng hydrothermal discharge mula sa bunganga ng Campi Flegrei35. Samakatuwid, napagpasyahan namin na ang mga fault at fracture sa Golpo ng Naples ay kumakatawan sa ginustong ruta para sa paglipat ng gas patungo sa ibabaw, isang katangiang ibinahagi ng iba pang structurally controlled hydrothermal. mga sistema36,37. Kapansin-pansin, ang mga BdM cone at hukay ay hindi palaging nauugnay sa mga tambak (Fig. 3a,c). Ipinahihiwatig nito na ang mga tambak na ito ay hindi kinakailangang kumakatawan sa mga precursor sa pagbuo ng hukay, tulad ng iminungkahi ng ibang mga may-akda para sa mga gas hydrate zone32,33. Sinusuportahan ng aming mga konklusyon ang hypothesis na ang pagkagambala ng mga dome seafloor sediment ay hindi palaging humahantong sa pagbuo ng mga hukay.
Ang tatlong nakolektang emisyon ng gas ay nagpapakita ng mga kemikal na katangian na tipikal ng mga hydrothermal fluid, lalo na ang CO2 na may makabuluhang konsentrasyon ng mga reducing gas (H2S, CH4 at H2) at mga light hydrocarbon (lalo na ang benzene at propylene)38,39, 40, 41, 42, 43, 44, 45 (Talahanayan S1). Ang pagkakaroon ng mga atmospheric gas (tulad ng O2), na hindi inaasahang makikita sa mga submarine emission, ay maaaring dahil sa kontaminasyon mula sa hangin na natunaw sa tubig-dagat na nakikipag-ugnayan sa mga gas na nakaimbak sa mga plastik na kahon na ginagamit para sa sampling, dahil ang mga ROV ay kinukuha mula sa sahig ng karagatan patungo sa dagat upang mag-alsa. Sa kabaligtaran, ang mga positibong halaga ng δ15N at isang mataas na N2/Ar (hanggang 480) na mas mataas kaysa sa ASW (air-saturated water) ay nagmumungkahi na ang karamihan sa N2 ay ginawa mula sa mga extra-atmospheric na mapagkukunan, na naaayon sa pangunahing hydrothermal na pinagmulan ng mga gas na ito. Ang hydrothermal-volcanic na pinagmulan ng BdM gas ay kinukumpirma ng mga nilalaman ng CO2 at He at ng kanilang isotopic. mga lagda. Ang mga carbon isotope (δ13C-CO2 mula -0.93% hanggang +0.4%) at mga halaga ng CO2/3He (mula 1.7 × 1010 hanggang 4.1 × 1010) ay nagmumungkahi na ang mga sample ng BdM ay kabilang sa isang halo-halong trend ng mga fumarole sa paligid ng mga dulo ng mantle ng Golpo ng Naples at ang ugnayan sa pagitan ng mga gas na nalikha ng reaksyon (Larawan 6). Mas partikular, ang mga sample ng gas na BdM ay matatagpuan sa kahabaan ng mixing trend sa halos parehong lokasyon tulad ng mga likido mula sa katabing mga bulkan ng Campi Flegrei at Somma-Veusivus. Mas crustal ang mga ito kaysa sa mga fumarole ng Ischia, na mas malapit sa dulo ng mantle. Ang Somma-Vesuvius at Campi Flegrei ay may mas mataas na mga halaga ng 3He/4He (R/Ra sa pagitan ng 2.6 at 2.9) kaysa sa BdM (R/Ra sa pagitan ng 1.66 at 1.96; Talahanayan S1). Ipinahihiwatig nito na ang pagdaragdag at akumulasyon ng radyogenic. Nagmula siya sa parehong pinagmumulan ng magma na nagpakain sa mga bulkan ng Somma-Vesuvius at Campi Flegrei. Ang kawalan ng mga nakikitang organic carbon fractions sa mga emisyon ng BdM ay nagmumungkahi na ang mga organic sediment ay hindi kasangkot sa proseso ng pag-aalis ng gas sa BdM.
Batay sa datos na iniulat sa itaas at mga resulta mula sa mga eksperimental na modelo ng mga istrukturang parang simboryo na nauugnay sa mga rehiyong mayaman sa subsea gas, ang deep gas pressurization ay maaaring responsable para sa pagbuo ng mga kilometrong-scale na BdM dome. Upang matantya ang overpressure na Pdef na humahantong sa BdM vault, naglapat kami ng isang thin-plate mechanics model33,34 na ipinapalagay, mula sa nakalap na morphological at seismic data, na ang BdM vault ay isang subcircular sheet na may radius na mas malaki kaysa sa isang deformed soft viscous deposit. Ang vertical maximum displacement w at kapal h ng (Supplementary Fig. S1). Ang Pdef ay ang pagkakaiba sa pagitan ng kabuuang presyon at rock static pressure kasama ang water column pressure. Sa BdM, ang radius ay humigit-kumulang 2,500 m, ang w ay 20 m, at ang h maximum na tinantyang mula sa seismic profile ay humigit-kumulang 100 m. Kinakalkula namin ang Pdef 46Pdef = w64D/a4 mula sa relasyon, kung saan ang D ay ang flexural stiffness; Ang D ay ibinibigay ng (E h3)/[12(1 – ν2)], kung saan ang E ay ang Young's modulus ng deposito, ang ν ay ang Poisson's ratio (~0.5)33. Dahil hindi masukat ang mga mekanikal na katangian ng mga sediment ng BdM, itinakda namin ang E = 140 kPa, na isang makatwirang halaga para sa mga sediment na mabuhangin sa baybayin 47 na katulad ng BdM14,24. Hindi namin isinasaalang-alang ang mas mataas na mga halaga ng E na iniulat sa literatura para sa mga deposito ng silty clay (300 < E < 350,000 kPa)33,34 dahil ang mga deposito ng BDM ay pangunahing binubuo ng buhangin, hindi silt o silty clay24. Nakukuha namin ang Pdef = 0.3 Pa, na naaayon sa mga pagtatantya ng mga proseso ng pagtaas ng seafloor sa mga kapaligiran ng gas hydrate basin, kung saan ang Pdef ay nag-iiba mula 10-2 hanggang 103 Pa, na may mas mababang mga halaga na kumakatawan sa mababang w/a at/o ano. Sa BdM, ang pagbawas ng stiffness dahil sa lokal na gas saturation ng sediment at/o ang hitsura Ang mga dati nang bali ay maaari ring mag-ambag sa pagkabigo at kasunod na paglabas ng gas, na nagpapahintulot sa pagbuo ng mga naobserbahang istruktura ng bentilasyon. Ang mga nakolektang reflected seismic profile (Fig. 7) ay nagpapahiwatig na ang mga PS sediment ay umangat mula sa GSL, na nagtutulak pataas sa nakapatong na mga MS marine sediment, na nagresulta sa mga mound, fold, fault, at sedimentary cut (Fig. 7b,c). Ipinahihiwatig nito na ang 14.8 hanggang 12 ka taong gulang na pumice ay nakapasok sa mas batang MS layer sa pamamagitan ng pataas na proseso ng transportasyon ng gas. Ang mga katangiang morpolohikal ng istrukturang BdM ay makikita bilang resulta ng overpressure na nilikha ng fluid discharge na ginawa ng GSL. Dahil ang aktibong discharge ay makikita mula sa seafloor hanggang sa mahigit 170 m bsl48, ipinapalagay namin na ang fluid overpressure sa loob ng GSL ay lumampas sa 1,700 kPa. Ang pataas na paglipat ng mga gas sa mga sediment ay mayroon ding epekto ng pagkuskos ng materyal na nakapaloob sa MS, na nagpapaliwanag sa pagkakaroon ng mga magulong sediment sa mga gravity core na na-sample sa BdM25. Bukod dito, ang overpressure ng GSL ay lumilikha ng kumplikadong sistema ng bali (polygonal fault sa Fig. 7b). Sama-sama, ang morpolohiya, istruktura, at stratigraphic settlement na ito, na tinutukoy bilang "mga pagoda"49,50, ay orihinal na iniuugnay sa mga pangalawang epekto ng mga lumang pormasyon ng glacial, at kasalukuyang binibigyang-kahulugan bilang mga epekto ng tumataas na gas31,33 o mga evaporite50. Sa continental margin ng Campania, ang mga evaporative sediment ay kakaunti, kahit man lang sa loob ng pinakamataas na 3 km ng crust. Samakatuwid, ang mekanismo ng paglaki ng mga BdM pagoda ay malamang na kontrolado ng pagtaas ng gas sa mga sediment. Ang konklusyong ito ay sinusuportahan ng mga transparent na seismic facies ng pagoda (Fig. 7), pati na rin ang datos ng gravity core gaya ng naunang naiulat24, kung saan ang kasalukuyang buhangin ay sumasabog kasama ang 'Pomici Principali'25 at 'Naples Yellow Tuff'26 Campi Flegrei. Bukod dito, ang mga deposito ng PS ay sumalakay at nag-deform sa pinakamataas na layer ng MS (Fig. 7d). Ang istrukturang kaayusan na ito ay nagmumungkahi na ang pagoda ay kumakatawan sa isang uprising structure at hindi lamang isang gas pipeline. Kaya, Dalawang pangunahing proseso ang namamahala sa pagbuo ng pagoda: a) ang densidad ng malambot na sediment ay bumababa habang pumapasok ang gas mula sa ibaba; b) ang pinaghalong gas-sediment ay tumataas, na siyang naobserbahang pagtiklop, faulting at fracture na nagiging sanhi ng mga deposito ng MS (Larawan 7). Isang katulad na mekanismo ng pagbuo ang iminungkahi para sa mga pagoda na nauugnay sa mga gas hydrate sa South Scotia Sea (Antarctica). Ang mga BdM pagoda ay lumitaw nang pangkat-pangkat sa mga maburol na lugar, at ang kanilang patayong lawak ay may average na 70-100 m sa two-way travel time (TWTT) (Larawan 7a). Dahil sa pagkakaroon ng mga MS undulations at isinasaalang-alang ang stratigraphy ng BdM gravity core, hinuhulaan namin na ang edad ng pagbuo ng mga istruktura ng pagoda ay mas mababa sa humigit-kumulang 14-12 ka. Bukod dito, ang paglaki ng mga istrukturang ito ay aktibo pa rin (Larawan 7d) dahil ang ilang mga pagoda ay sumalakay at nag-deform sa nakapatong na kasalukuyang buhangin na BdM (Larawan 7d).
Ang pagkabigo ng pagoda na makatawid sa kasalukuyang seabed ay nagpapahiwatig na (a) ang pagtaas ng gas at/o lokal na paghinto ng paghahalo ng gas-sediment, at/o (b) ang posibleng lateral flow ng pinaghalong gas-sediment ay hindi nagpapahintulot ng isang lokalisadong proseso ng overpressure. Ayon sa diapir theory model52, ang lateral flow ay nagpapakita ng negatibong balanse sa pagitan ng rate ng supply ng pinaghalong putik-gas mula sa ibaba at ang rate kung saan ang pagoda ay gumagalaw pataas. Ang pagbawas sa rate ng supply ay maaaring nauugnay sa pagtaas ng density ng pinaghalong dahil sa pagkawala ng supply ng gas. Ang mga resultang nakabuod sa itaas at ang buoyancy-controlled rise ng pagoda ay nagbibigay-daan sa amin na tantyahin ang taas ng air column na hg. Ang buoyancy ay ibinibigay ng ΔP = hgg (ρw – ρg), kung saan ang g ay gravity (9.8 m/s2) at ang ρw at ρg ay ang mga densidad ng tubig at gas, ayon sa pagkakabanggit. Ang ΔP ay ang kabuuan ng dating kinalkulang Pdef at ang lithostatic pressure na Plith ng sediment plate, i.e. ρsg h, kung saan ang ρs ay ang densidad ng sediment. Sa kasong ito, ang halaga ng hg na kinakailangan para sa ninanais na buoyancy ay ibinibigay ng hg = (Pdef + Plith)/[g (ρw – ρg)]. Sa BdM, itinakda natin ang Pdef = 0.3 Pa at h = 100 m (tingnan sa itaas), ρw = 1,030 kg/m3, ρs = 2,500 kg/m3, ang ρg ay bale-wala dahil ρw ≫ρg. Nakukuha natin ang hg = 245 m, isang halagang kumakatawan sa lalim ng ilalim ng GSL. Ang ΔP ay 2.4 MPa, na siyang overpressure na kinakailangan upang mabasag ang BdM seafloor at bumuo ng mga vent.
Ang komposisyon ng BdM gas ay naaayon sa mga pinagmumulan ng mantle na binago ng pagdaragdag ng mga pluido na nauugnay sa mga reaksyon ng decarbonization ng mga crustal rock (Larawan 6). Ang magaspang na pagkakahanay ng EW ng mga BdM dome at mga aktibong bulkan tulad ng Ischia, Campi Flegre, at Soma-Vesuvius, kasama ang komposisyon ng mga gas na inilalabas, ay nagmumungkahi na ang mga gas na inilalabas mula sa mantle sa ilalim ng buong rehiyon ng bulkan ng Naples ay magkakahalo. Parami nang parami ang mga crustal fluid na lumilipat mula kanluran (Ischia) patungong silangan (Somma-Vesuivus) (Mga Larawan 1b at 6).
Napagpasyahan namin na sa Look ng Naples, ilang kilometro mula sa daungan ng Naples, mayroong isang 25 km2 na lapad na parang simboryo na istraktura na apektado ng isang aktibong proseso ng degassing at sanhi ng paglalagay ng mga pagoda at mga tambak. Sa kasalukuyan, ang mga lagda ng BdM ay nagmumungkahi na ang non-magmatic turbulence53 ay maaaring nauna pa sa embryonic volcanism, ibig sabihin, ang maagang paglabas ng magma at/o mga thermal fluid. Dapat ipatupad ang mga aktibidad sa pagsubaybay upang suriin ang ebolusyon ng mga phenomena at upang matukoy ang mga geochemical at geophysical signal na nagpapahiwatig ng mga potensyal na kaguluhan sa magmatic.
Ang mga acoustic water column profile (2D) ay nakuha noong SAFE_2014 (Agosto 2014) cruise sa R/V Urania (CNR) ng National Research Council Institute of Coastal Marine Environment (IAMC). Ang acoustic sampling ay isinagawa gamit ang isang siyentipikong beam-splitting echo sounder na Simrad EK60 na tumatakbo sa 38 kHz. Ang acoustic data ay naitala sa average na bilis na humigit-kumulang 4 km. Ang mga nakolektang imahe ng echosounder ay ginamit upang matukoy ang mga fluid discharge at tumpak na tukuyin ang kanilang lokasyon sa collection area (sa pagitan ng 74 at 180 m bsl). Sukatin ang mga pisikal at kemikal na parameter sa water column gamit ang mga multiparameter probe (conductivity, temperature and depth, CTD). Ang data ay kinolekta gamit ang isang CTD 911 probe (SeaBird, Electronics Inc.) at pinoproseso gamit ang SBED-Win32 software (Seasave, bersyon 7.23.2). Ang visual inspection ng seabed ay isinagawa gamit ang isang "Pollux III" (GEItaliana) ROV device (remotely operated vehicle) na may dalawang (low at high definition) na kamera.
Isinagawa ang pagkuha ng datos mula sa multibeam gamit ang isang 100 KHz Simrad EM710 multibeam sonar system (Kongsberg). Ang sistema ay nakakonekta sa isang differential global positioning system upang matiyak ang mga sub-metric error sa pagpoposisyon ng beam. Ang acoustic pulse ay may frequency na 100 KHz, firing pulse na 150° degrees at isang buong bukana ng 400 beam. Sukatin at ilapat ang mga profile ng bilis ng tunog sa real time habang kinukuha. Ang datos ay pinoproseso gamit ang PDS2000 software (Reson-Thales) ayon sa pamantayan ng International Hydrographic Organization (https://www.iho.int/iho_pubs/standard/S-44_5E.pdf) para sa nabigasyon at pagwawasto ng pagtaas at pagbaba ng tubig. Ang pagbabawas ng ingay dahil sa mga aksidenteng pagtalsik ng instrumento at mababang kalidad na pagbubukod ng beam ay isinagawa gamit ang mga tool sa pag-edit ng banda at pag-de-spiking. Ang patuloy na pagtukoy ng bilis ng tunog ay isinasagawa ng isang istasyon ng keel na matatagpuan malapit sa multi-beam transducer at kumukuha at naglalapat ng mga real-time na profile ng bilis ng tunog sa hanay ng tubig tuwing 6-8 oras upang magbigay ng real-time na bilis ng tunog para sa wastong... beam steering. Ang buong dataset ay binubuo ng humigit-kumulang 440 km2 (0-1200 m ang lalim). Ginamit ang datos upang magbigay ng isang high-resolution digital terrain model (DTM) na nailalarawan sa pamamagitan ng 1 m na laki ng grid cell. Ang pangwakas na DTM (Fig. 1a) ay ginawa gamit ang datos ng terrain (>0 m sa ibabaw ng antas ng dagat) na nakuha sa 20 m na laki ng grid cell ng Italian Geo-Military Institute.
Isang 55-kilometrong high-resolution single-channel seismic data profile, na nakolekta noong mga ligtas na cruise sa karagatan noong 2007 at 2014, ang sumasaklaw sa humigit-kumulang 113 kilometro kuwadrado, kapwa sa R/V Urania. Ang mga Marisk profile (hal., L1 seismic profile, Fig. 1b) ay nakuha gamit ang IKB-Seistec boomer system. Ang acquisition unit ay binubuo ng isang 2.5 m catamaran kung saan nakalagay ang source at receiver. Ang source signature ay binubuo ng isang positibong peak na nailalarawan sa frequency range na 1-10 kHz at nagbibigay-daan upang ma-resolve ang mga reflector na pinaghihiwalay ng 25 cm. Ang mga ligtas na seismic profile ay nakuha gamit ang isang 1.4 Kj multi-tip Geospark seismic source na naka-interface sa Geotrace software (Geo Marine Survey System). Ang sistema ay binubuo ng isang catamaran na naglalaman ng 1–6.02 KHz source na tumatagos hanggang 400 milliseconds sa malambot na sediment sa ilalim ng seabed, na may theoretical vertical resolution na 30 cm. Ang parehong Safe at Marsik device ay nakuha sa isang rate na 0.33 shots/seg na may vessel velocity na <3 Kn. Ang datos ay pinoproseso at ipinakita gamit ang Geosuite Allworks software na may sumusunod na daloy ng trabaho: dilation correction, water column muting, 2-6 KHz bandpass IIR filtering, at AGC.
Ang gas mula sa fumarole sa ilalim ng tubig ay kinolekta sa sahig ng dagat gamit ang isang plastik na kahon na may goma na diaphragm sa itaas na bahagi nito, na inilagay nang patiwarik ng ROV sa ibabaw ng vent. Kapag ang mga bula ng hangin na pumapasok sa kahon ay ganap nang napalitan ang tubig-dagat, ang ROV ay babalik sa lalim na 1 m, at inililipat ng maninisid ang nakolektang gas sa pamamagitan ng isang goma na septum papunta sa dalawang pre-evacuated na 60 mL na glass flasks na may Teflon stopcocks kung saan ang isa ay pinuno ng 20 mL ng 5N NaOH solution (Gegenbach-type flask). Ang mga pangunahing acid gas species (CO2 at H2S) ay natutunaw sa alkaline solution, habang ang mga low solubility gas species (N2, Ar+O2, CO, H2, He, Ar, CH4 at light hydrocarbons) ay nakaimbak sa sampling bottle headspace. Ang mga inorganic low solubility gas ay sinuri sa pamamagitan ng gas chromatography (GC) gamit ang isang Shimadzu 15A na may 10 m na haba na 5A molecular sieve column at isang thermal conductivity detector (TCD) 54. Ang Argon at O2 ay Sinuri gamit ang isang Thermo Focus gas chromatograph na may 30 m na haba na capillary molecular sieve column at TCD. Ang methane at light hydrocarbons ay sinuri gamit ang isang Shimadzu 14A gas chromatograph na may 10 m na haba na stainless steel column na may Chromosorb PAW 80/100 mesh, pinahiran ng 23% SP 1700 at isang flame ionization detector (FID). Ang liquid phase ay ginamit para sa pagsusuri ng 1) CO2, as, na tinitrate gamit ang 0.5 N HCl solution (Metrohm Basic Titrino) at 2) H2S, as, pagkatapos ng oksihenasyon gamit ang 5 mL H2O2 (33%), sa pamamagitan ng ion chromatography (IC) (IC) (Wantong 761). Ang analytical error ng titration, GC at IC analysis ay mas mababa sa 5%. Pagkatapos ng mga karaniwang pamamaraan ng pagkuha at paglilinis para sa mga pinaghalong gas, ang 13C/12C CO2 (ipinahayag bilang δ13C-CO2% at V-PDB) ay sinuri gamit ang isang Finningan Delta S mass spectrometer55,56. Ang mga pamantayang ginamit upang tantyahin ang panlabas na katumpakan ay ang Carrara at San Vincenzo marble (panloob), NBS18 at NBS19 (internasyonal), habang ang analytical error at reproducibility ay ±0.05% at ±0.1%, ayon sa pagkakabanggit.
Ang mga halaga ng δ15N (ipinahayag bilang % vs. Air) at 40Ar/36Ar ay natukoy gamit ang isang Agilent 6890 N gas chromatograph (GC) na nakakabit sa isang Finnigan Delta plusXP continuous flow mass spectrometer. Ang error sa pagsusuri ay: δ15N±0.1%, 36Ar<1%, 40Ar<3%. Ang ratio ng He isotope (ipinahayag bilang R/Ra, kung saan ang R ay 3He/4He na nasukat sa sample at ang Ra ay parehong ratio sa atmospera: 1.39 × 10−6)57 ay natukoy sa laboratoryo ng INGV-Palermo (Italya). Ang 3He, 4He at 20Ne ay natukoy gamit ang isang dual collector mass spectrometer (Helix SFT-GVI)58 pagkatapos ng paghihiwalay ng He at Ne. Ang error sa pagsusuri ay ≤ 0.3%. Ang karaniwang mga blangko para sa He at Ne ay <10-14 at <10-16 mol, ayon sa pagkakabanggit.
Paano banggitin ang artikulong ito: Passaro, S. et al. Ang pagtaas ng sahig ng dagat na dulot ng proseso ng degassing ay nagpapakita ng namumuong aktibidad ng bulkan sa baybayin.science.Rep. 6, 22448; doi: 10.1038/srep22448 (2016).
Aharon, P. Ang heolohiya at biyolohiya ng moderno at sinaunang mga pagtagas at butas ng hydrocarbon sa sahig ng dagat: isang panimula. Geographic Ocean Wright.14, 69–73 (1994).
Paull, CK at Dillon, WP Ang pandaigdigang paglitaw ng mga gas hydrate. Sa Kvenvolden, KA at Lorenson, TD (mga patnugot) 3–18 (Mga natural na gas hydrate: Paglitaw, distribusyon at pagtuklas. American Geophysical Union Geophysical Monograph 124, 2001).
Fisher, AT Mga limitasyong heopisikal sa sirkulasyong hidrotermal. Sa: Halbach, PE, Tunnicliffe, V. & Hein, JR (mga patnugot) 29–52 (Ulat ng Durham Workshop, Paglilipat ng Enerhiya at Masa sa mga Sistemang Hydrothermal ng Dagat, Durham University Press, Berlin (2003)).
Coumou, D., Driesner, T. & Heinrich, C. Kayarian at dinamika ng mga sistemang hydrothermal ng mid-ocean ridge. Science 321, 1825–1828 (2008).
Boswell, R. & Collett, TS Mga Kasalukuyang Pananaw sa mga Yaman ng Gas Hydrate.Enerhiya.At Agham Pangkapaligiran.4, 1206–1215 (2011).
Evans, RJ, Davies, RJ & Stewart, SA Panloob na istruktura at kasaysayan ng pagsabog ng isang sistema ng bulkang putik na may sukat na kilometro sa South Caspian Sea. Basin Reservoir 19, 153–163 (2007).
Leon, R. et al. Mga katangian ng sahig ng dagat na nauugnay sa pagtagas ng mga hydrocarbon mula sa mga tambak ng putik na carbonate sa malalim na tubig sa Golpo ng Cadiz: mula sa daloy ng putik hanggang sa mga sediment na carbonate. Heograpiya March. Wright.27, 237–247 (2007).
Moss, JL at Cartwright, J. 3D na representasyon ng seismic ng mga pipeline ng pagtakas ng pluido sa iskala kilometro sa laot ng Namibia. Basin Reservoir 22, 481–501 (2010).
Andresen, KJ Mga katangian ng daloy ng pluwido sa mga sistema ng pipeline ng langis at gas: Ano ang sinasabi ng mga ito sa atin tungkol sa ebolusyon ng basin? March Geology.332, 89–108 (2012).
Ho, S., Cartwright, JA & Imbert, P. Bertikal na ebolusyon ng istruktura ng paglabas ng Neogene Quaternary fluid kaugnay ng mga gas flux sa Lower Congo Basin, sa laot ng Angola.March Geology.332–334, 40–55 (2012).
Johnson, SY et al. Aktibidad na hidrotermal at tektoniko sa hilagang Yellowstone Lake, Wyoming.geology.Socialist Party.Yes.bull.115, 954–971 (2003).
Patacca, E., Sartori, R. & Scandone, P. Ang Tyrrhenian Basin at ang Apennine Arc: Mga Ugnayang Kinematiko Simula noong Huling Totonian. Mem Soc Geol Ital 45, 425–451 (1990).
Milia et al. Tektoniko at istrukturang crustal sa continental margin ng Campania: kaugnayan sa aktibidad ng bulkan.mineral.gasoline.79, 33–47 (2003)
Piochi, M., Bruno PP at De Astis G. Ang relatibong papel ng rift tectonics at mga proseso ng magmatic uplift: hinuha mula sa geophysical, structural, at geochemical data sa rehiyon ng bulkan ng Naples (timog Italya). Gcubed, 6(7), 1-25 (2005).
Dvorak, JJ & Mastrolorenzo, G. Mga mekanismo ng kamakailang patayong paggalaw ng crustal sa bunganga ng Campi Flegrei sa timog Italya.geology.Socialist Party.Yes.Specification.263, pp. 1-47 (1991).
Orsi, G. et al. Panandaliang deformasyon ng lupa at seismicity sa nested Campi Flegrei crater (Italy): isang halimbawa ng aktibong mass recovery sa isang densely populated na lugar. J. Volcano.geothermal.reservoir.91, 415–451 (1999)
Cusano, P., Petrosino, S., at Saccorotti, G. Mga pinagmulang hidrotermal ng patuloy na pangmatagalang aktibidad na 4D sa bulkaniko ng Campi Flegrei sa Italya. J. Volcano.geothermal.reservoir.177, 1035–1044 (2008).
Pappalardo, L. at Mastrolorenzo, G. Mabilis na pagkakaiba-iba sa mga mala-sill na imbakan ng magmatic: isang pag-aaral ng kaso mula sa bunganga ng Campi Flegrei.science.Rep. 2, 10.1038/srep00712 (2012).
Walter, TR et al. Ang InSAR time series, correlation analysis, at time-correlation modeling ay nagpapakita ng posibleng pagkakaugnay ng Campi Flegrei at Vesuvius.J. Volcano.geothermal.reservoir.280, 104–110 (2014).
Milia, A. & Torrente, M. Estruktural at estratigrapikong istruktura ng unang kalahati ng Tyrrhenian graben (Golpo ng Naples, Italya). Constructive Physics 315, 297–314.
Sano, Y. & Marty, B. Mga Pinagmumulan ng Carbon sa Abo ng Bulkan mula sa Island Arcs.Chemical Geology.119, 265–274 (1995).
Milia, A. Dohrn Canyon stratigraphy: Mga tugon sa pagbaba ng lebel ng dagat at tectonic uplift sa panlabas na continental shelf (Silangang Tyrrhenian margin, Italy). Geo-Marine Letters 20/2, 101–108 (2000).