Cnr.it INDIVIDUATO SOTTO il fondale del Mar Jonio un sistema di faglie che ha controllato l’evoluzione dell’antico oceano della Tetide e che è ancora in grado di innescare processi vulcanici e sismici. A svelarlo, uno studio coordinato dall’Istituto di scienze marine del Cnr di Bologna, in collaborazione con Università di Parma, Ingv e Geomar (Germania), pubblicato su Nature Communications.
Un sistema di spaccature profonde sta separando la Sicilia dal resto dell’Italia nella regione compresa tra lo stretto di Messina e l’Etna. Lungo queste strutture geologiche risale materiale del mantello che formava il basamento dell’oceano mesozoico, chiamato Tetide, da una profondità di circa 15-20 chilometri. Si tratta di una vera e propria finestra sotto il fondale del Mar Jonio, che consente di osservare da vicino blocchi dell’antico oceano, svelando i processi che hanno portato alla sua formazione. Lo studio Lower plate serpentinite diapirism in the Calabrian Arc subduction complex, condotto da un team di ricercatori dell’Istituto di scienze marine del Consiglio nazionale delle ricerche (Ismar-Cnr) di Bologna, dell’Università di Parma, dell’Istituto nazionale di geofisica e vulcanologia (Ingv) e del Geomar (Kiel, Germania), è stato pubblicato su Nature Communications.
“Le faglie lungo le quali risale il mantello della Tetide”, spiega Alina Polonia, ricercatrice Ismar-Cnr e coordinatrice della ricerca, “controllano anche la formazione del Monte Etna, dimostrando che si tratta di strutture in grado di innescare processi vulcanici e causare terremoti. Queste faglie, infatti, sono profonde e lunghe decine di chilometri, e separano blocchi di crosta terrestre in movimento reciproco”.
Attraverso uno studio multi-disciplinare, che integra immagini acustiche del sottosuolo, dati geofisici e campioni di sedimento, acquisiti nel corso di spedizioni scientifiche con la nave oceanografica del Cnr Urania, è stato possibile identificare le faglie, ricostruire la loro geometria e scoprire anomalie geochimiche nei sedimenti legate alla presenza di fluidi profondi. L’analisi di tutti i dati raccolti ha permesso di proporre un modello geologico che conferma l’origine profonda del materiale in risalita lungo le faglie.
“Grazie a questa scoperta”, prosegue Alina Polonia, “l’Arco Calabro, il sistema di subduzione tra Africa ed Europa nel Mar Jonio, ha un importante primato: è l’unica regione al mondo in cui sia stato descritto materiale del mantello in risalita dalla placca in subduzione. Questa scoperta avrà importanti implicazioni per capire meglio come si formano le catene montuose e come questi processi siano legati ai forti terremoti storici registrati in Sicilia e Calabria”.
Lower plate serpentinite diapirism in the Calabrian Arc subduction complex
Mantle-derived serpentinites have been detected at magma-poor rifted margins and above subduction zones, where they are usually produced by fluids released from the slab to the mantle wedge. Here we show evidence of a new class of serpentinite diapirs within the external subduction system of the Calabrian Arc, derived directly from the lower plate. Mantle serpentinites rise through lithospheric faults caused by incipient rifting and the collapse of the accretionary wedge. Mantle-derived diapirism is not linked directly to subduction processes. The serpentinites, formed probably during Mesozoic Tethyan rifting, were carried below the subduction system by plate convergence; lithospheric faults driving margin segmentation act as windows through which inherited serpentinites rise to the sub-seafloor. The discovery of deep-seated seismogenic features coupled with inherited lower plate serpentinite diapirs, provides constraints on mechanisms exposing altered products of mantle peridotite at the seafloor long time after their formation.
Tectonic activity in subduction-rollback convergent plate boundaries commonly involves backward migration of subducting slabs, arc magmatism, and plate boundary segmentation. Incipient collision, the presence of indentors, and generally lower plate topography may increase complexity and plate fragmentation. Slab tearing has a major role in developing segmented subduction zones, in particular in narrow slab segments1, triggering asthenospheric upwelling, dynamic topography2, and magmatism3,4.
The origin of the Ionian Sea is related to the opening of the Neo-Tethys basin5, which has now mostly disappeared below the Alps/Apennines. An exceptionally thick sedimentary cover (up to 6–8 km) and the lack of detailed magnetic anomalies makes it difficult to define the nature of the Ionian Sea crust, which has been interpreted as either oceanic6,7,8 or oceanic/severely stretched continental9,10, or thinned continental11,12,13.
Consumption of the Neo-Tethys lithosphere since Oligocene times has caused the emplacement of the Calabrian Arc (CA), a narrow and arcuate subduction system (Fig. 1) related to Africa/Eurasia plate convergence, and to the southeastward retreat of the Tethyan slab14. Convergence generated a 10–30 km-thick, 300 km-wide subduction system, including a submarine accretionary wedge15. The migrating Calabria trench drives the entire subduction complex outward, but different boundary conditions produce segmentation in two distinct lobes, with different structural and morphotectonic characteristics (Fig. 1). The eastern lobe (EL) collides with the Hellenic subduction system and produces basement-involved tectonics similar to a fold and thrust belt15, whereas the western lobe (WL) is free to spread into the abyssal plain of the Ionian Sea (Fig. 1) through frontal accretion and offscraping processes 16.
Two oppositely dipping fault systems, i.e., the Ionian fault (IF) and Alfeo-Etna fault (AEF), deform the WL offshore Eastern Sicily (Fig. 1) along a complex strike-slip/transtensional pattern15,17. These two systems mark the external boundary of a wide deformation zone including three additional faults (F1, F2, and F3; Fig. 1) and they might be the source of major earthquakes in the CA18. Trans-tensional deformation and deep fragmentation of the Western Ionian domain19,20,21 are in agreement with geodetic measurements22, suggesting that Calabria is moving SE and Southern Sicily is moving South in an Apulia-fixed reference frame. Moreover, mechanical/numerical models taking into account geometry and strength contrast across the Neo-Tethyan passive margin and major faults, suggest that pre-existing fracture zones and transform faults controlled the subsequent Ionian margin segmentation17. In this context, the NW portion of the IF corresponds to basement deformation associated with slab tearing, whereas the AEF might be the result of regional scale dextral shearing.
Seismic reflection profiles integrated with multibeam and potential field data reveal 13 sub-circular magnetic bodies beneath the seafloor aligned along transtensional faults segmenting the continental margin. We investigated the nature of such bodies within the accretionary wedge, their link with Mesozoic paleoceanography and with the latest Quaternary extensional phase that disrupted the subduction complex. Our multidisciplinary approach suggests that the diapirs are made of inherited Mesozoic Tethyan serpentinites, which rise along transtensive faults controlling rifting processes perpendicular to the subduction arc. They provide the first example of serpentinite diapirism from the lower plate in a subduction system. Mantle-derived serpentinites combined with geophysical data suggest that the ocean-continent transition (OCT) in the Ionian Sea has a complex geometry reflecting the 90° curvature of the continental margin. Oceanic crust, if present, should be limited to the inner accretionary wedge of the subduction complex and the AEF should be located close to the Mesozoic Tethyan OCT and/or to a major transform fault.
Rifting and diapirism along the AEF system. Multibeam, seismic reflection, magnetic, and gravity field data call for incipient deformation of four areas marked by buried sub-circular features aligned along the AEF (Fig. 2). Area-1 is offshore Mt. Etna, whereas Area-2 and -3 are located north and northeast of Alfeo Seamount; the more distal Area-4 is located on a slope terrace of the accretionary complex. Deep deformation in these areas is imaged by seismic reflection lines showing sediment tilting and up warping around acoustically blind nuclei (Fig. 3). More information can be found in Supplementary Figs. 1–7. Multibeam data show sub-circular seafloor patterns, marked by increased roughness particularly in Area-2 (Fig. 2).
Magnetic anomaly map derived from this study superimposed over a gray levels bathymetric slope map. The four major areas of mantle diapirism described in this work are indicated (black squares 1–4). The inset highlights the rough seafloor topography of Area-2, which corresponds to the magnetic anomaly. White lines correspond to MCS data shown in Fig. 3, whereas gray lines are shown in the supplementary material. Blue stars represent heat flow measurements from ref. 23 (17: 15.7 ± 1.9 mW m−2; 18: 16.2 ± 2.7 mW m−2; 29: 26.1 ± 5.2 mW m−2; 35: 30.3 ± 5.0 mW m−2; 36: 34.4 ± 4.0 mW m−2). Red squares a–drepresent gravity cores described in the manuscript (Fig. 7) and Supplementary Figs. 9–12
Line drawings of time migrated MCS lines used for structural reconstructions in Areas-1 and -2. Location of seismic profiles is shown in Fig. 2. Yellow/brown pattern: diapiric features. a ETNASEIS-6 MCS line collected orthogonal to the continental margin SE of Etna. This line crosses the entire WL of the accretionary wedge intercepting both the Alfeo-Etna and Ionian faults. Such faults control extensional processes, which are associated with a shallower Moho at about 11 s TWT. A second diapiric feature is present in the centre of the rifting zone. b CA-A MCS line collected across the WL of the Calabrian Arc subduction complex orthogonally to the AEF. A reflector at about 11 s TWT is interpreted as the Moho reflection. c MS-26 MCS line across the AEF north of the Alfeo Seamount. M, multiple. The diapir corresponds to a very rough seafloor morphology. Bold lines: major trans-tensional faults controlling rifting processes
Deep penetrating multichannel seismic lines collected in Areas-1 and -2 (Fig. 3) run across the subsiding WL of the subduction complex, from the Malta escarpment to the IF, and image crustal trans-tensional faults segmenting the accretionary wedge. The AEF forms at the boundary between the salt-bearing complex emplaced after the Messinian salinity crisis and the older wedge consisting of clastic sediments (Fig. 3b). This fault system served as main thrust fault during Messinian time, when frontal accretion resumed after the desiccation of the Mediterranean Sea15, and later accommodated transtensive deformation, as indicated by a syn-tectonic sedimentary basin in the footwall (Fig. 3b). Corresponding to the AEF, the basal detachment is disrupted, displaced and associated with large 10–20 km-wide diapirs. A reflector at about 13 s two-way travel time (TWT) in seismic line ETNASEIS 6 and 11 s (TWT) in lines CA-A is interpreted as the Moho (Fig. 3a, b). Arching and disruption of deep reflectors in correspondence to the diapiric structures suggest that the source region is deep and below the detachment.
Rising material in subduction complexes might be ascribed to different processes, such as salt/mud diapirism, magmatic intrusions, and serpentinite diapirism. We attempted to discriminate between these processes by gravity and magnetic modeling (Supplementary Fig. 8). A salt source for the diapirs can be excluded because no thick salt was deposited in the inner wedge during Messinian times15, and because the source of the rising material is deeper than the Messinian salt layer (Fig. 3). Moreover, salt and mud diapirs would not cause positive magnetic anomalies; gravity modeling also excludes a mud composition for the diapirs (Supplementary Fig. 8). As these structures are aligned offshore Mt. Etna, an origin related to magmatism needs to be explored. However, magnetic and gravity anomalies are not consistent with a volcanic or magmatic source (Supplementary Fig. 8) in agreement with heat flow data showing a lower than normal thermal regime for the Ionian basin (mean heat flow value is 31.8 ± 5.0 mW m−223) (Fig. 2). As the basal detachment of the subduction system is deformed and segmented in correspondence with these diapirs (Fig. 3), we argue that the source of the rising material is below the detachment; thus, a possible serpentinite composition for the diapirs was investigated (Fig. 4).
Forward magnetic and gravity modeling. a Observed, calculated, and misfit (error) curves for magnetic modeling along seismic line CA-A (Fig. 3b). b Observed, calculated, and misfit (error) curves for gravimetric modeling along seismic line CA-A (Fig. 3b). c Geometry of causative rocks derived from depth-migrated seismic section used as background. Average fitting errors were 5.48 nT (5.5 %) and 1.354 mGal (2.5 %), for magnetic and gravity profiles, respectively. Red line identifies the depth distribution of the local geothermal gradient
Joint interpretation of magnetic and gravity data has been achieved computing a 2.75D forward modeling crossing the diapir in Area 2. The forward model was geologically constrained in depth through seismic data. The 2.75D forward model was approached redrawing the geometry of sedimentary sequences directly in TWT sections. Time–depth conversion was achieved assuming reliable seismic velocities for each block (Fig. 4). The best fit (5.48 nT and 1.354 mGal for magnetic and gravity errors, respectively) was obtained modeling the intrusion as a serpentine diapir with low magnetization (0.015 SI) and density 2.7 g cm−324. Magnetic modeling was limited in depth considering a thermal boundary approximately at 16–17 km estimated assuming a linear increase along a normal thermal gradient.
Combined P- and S-wave velocity profiles are required to assess whether the seismic structure derived from deep penetration multichannel seismic and potential field data are compatible with a variable abundance of mafic and mantle-derived serpentinite25. P- and S- wave profiles are available in the southern region of the diapiric field26. Simultaneous inversion of P- and S-wave arrival times, collected during a 3-year ocean bottom seismometer/hydrophone (OBS/H) monitoring campaign (Fig. 5 green stars), yields one-dimensional (1D) P- and S-wave velocity models for the Ionian lithosphere. The 1D model indicates the presence in the Ionian upper mantle of two layers with high seismic P-wave velocity, low S-wave velocity, and high Vp/Vs values, interpreted as partly serpentinized peridotites with 55–65% and 15–25% serpentinization, respectively26. These values of Vp/Vs for serpentinites are in agreement with data in the literature25 and the degree of serpentinization is similar to that derived from gravity and magnetic data. Serpentinization of peridotites produces an increase in susceptibility, because magnetite forms during serpentinization by the breakdown of the iron-rich olivines27,28. Serpentinization is also responsible for a density decrease, which follows a linear trend related to the degree of alteration27. Considering that serpentinization modifies the physical parameter of the original peridotites, we use our forward model to deduce degrees of serpentinization. According to our model the density of the intruded body is of the order of 2.7 g cm−3 with susceptibility of 0.015 SI, suggesting a serpentinization of about 75%29 and 68%28, respectively. These values, obtained by two independent data sets, are similar to those obtained from Vp/Vs data26, supporting our hypothesis of serpentinite intrusions below the Ionian seafloor.
Structural map of the Western Calabrian Arc subduction complex with distribution of tectonically controlled serpentinite diapirs. Faulting associated with arc-perpendicular extensional tectonics in the accretionary wedge (modified from ref. 17) provides an opportunity for fluids to be introduced deep into the underplating plate. A-A’: MCS line CALA-02 across a diapir in the sedimentary basin formed along the Alfeo-Etna fault. A diapiric intrusion of serpentinite bearing clayey material in the Hyblean plateau suggests affinity of clays with hydrothermally modified mafic and ultramafic rocks forming the Hyblean lower crust43,78. Sections B-B’ and C-C’ are shown in Fig. 6, section D-D’ in Fig. 8. Green stars: OBS/H stations used for P- and S-wave velocity model for the Ionian lithosphere26. AEF, Alfeo-Etna fault; IF, Ionian fault; OCT, Oceanic-Continent transition. F1, F2, F3, transtensional faults belonging to the AEF-IF deformation zone 17
Comparing these results with our pre-stack depth migrated multichannel seismic data in the same region we find a good correlation between Vp/Vs-derived serpentinite layers and our interpretation of the Ionian lithosphere. In particular, the two serpentinite layers derived from the Vp/Vs data fit well with the source area of the diapirs (Fig. 6b). The geometry and rooting in the Moho of the diapirs combined with the two high Vp/Vs layers at the same source depth support a serpentinite origin for the rising material.
Structural cross-sections of the continental margin. a Sketch based on MCS line MS-107 showing the margin structure parallel to the arc where lithospheric transtensive faults trigger the rise of serpentinite material from the lower plate. b Sketch based on MCS line CROP-M2B showing the margin structure orthogonal to the arc where the geometry of the subduction thrust is imaged. The presence of serpentinite units and Vp/Vs data in the outer part of the subduction complex is taken from26. OCT, Oceanic-Continent transition. Location of seismic sections in Fig. 5
Pore water geochemistry
Pore water data from mud volcanoes and other fluid dewatering features are a common tool to constrain subsurface geochemical processes30,31,32. To check for potential geochemical anomalies that may be related to the advection of deep-source fluids, pore waters from gravity cores were analyzed from above the diapiric structure and compared with a reference site (core d) outside this range (Fig. 7 and Supplementary Figs. 9–12).
Depth profiles of selected pore water constituents in the sediment cores. Core a: SF17-04; core b: GH15-03; core c: SF17-05; core d: CQ14-02 (reference core not affected by diapirism). Core a and c from above the diapir show deviations from normal diagenetic trends (core d). In both cores, pore water Sr concentrations increase continuously with depth; in core a, this trend is accompanied by depletion in Cl (and other elements that typically behave conservative at this depth below the seafloor). These geochemical anomalies provide evidence for fluid advection (affected probably by mineral dewatering and other deep-seated geochemical processes) and suggest a direct relation to the observed diapirism. Core location in Fig. 2
Core a and core c from above the diapiric structures show deviations from normal diagenetic trends. Specifically, depletion in elements such as Cl behaving conservatively under normal diagenetic conditions, suggests freshening of pore water, e.g., caused by (clay) mineral dehydration. This freshening is typically associated with enrichment of fluid mobile elements, such as Li, B, and Sr31,32. Core c shows a significant (10%) depletion in chlorinity and a decrease in SO4 and Ca, both indicating that a major horizon of anaerobic oxidation of methane is present a few meters below the seafloor. In addition, core-a and core-c show enrichment in dissolved Sr above background values. However, as Sr-enrichment is not accompanied by a freshening trend in both cores and B and Li are not enriched, other processes such as Sr-release from recrystallized carbonates32 can be suggested. Similar geochemical anomalies have been found in fluids emanating from serpentinite mud volcanoes in the Mariana subduction zone, where slab-derived fluids drive serpentinization and mud volcanism33,34. However, serpentinization does not produce a direct typical geochemical signature in sedimentary environments, given variables as composition of the source fluid, nature of ambient sediments, reaction temperatures, etc., causing a wide range of potential variations.
In summary, the observed pore water signature indicates active fluid migration related to subsurface dewatering and associated processes. This is confirmed by sediment facies analyses showing fluid flow related features such as patchy cloudy facies35, sediment layers disruption, mud injections, and vertical fluid paths (Supplementary Fig. 10).
The CA subduction system is close to a continental collision that causes margin segmentation and disruption of the subduction system driven by pre-existing heterogeneities in the forearc basement17. Crustal velocity models obtained from deep seismic sounding (DSS) data12 account for a Moho depth decreasing from about 35 km in the EL, to some 20 km in the WL, close to the Malta escarpment. Moho depth variations occur in the region between the AEF and IF, where our seismic reflection profiles (Fig. 3) show a Moho shallower than in the surrounding regions, in agreement with other evidence, suggesting a depth between 12 and 15 km10. This geometry suggests an asymmetric rifting affecting a large corridor between the AEF and IF, and the collapse of the WL of the accretionary wedge. Thus, the AEF may represent the western master fault accommodating plate divergence along a dextral shear zone, triggering serpentinite diapirism.
Differences in plumbing systems and hydrogeological characteristics of forearc basements may trigger contrasting processes in different forearcs, including fluid flow and mud volcanism (i.e., Nankai Trough or Mediterranean Ridge36) or serpentinite diapirism (Izu-Bonin forearc37,38,39), whereas impermeable and consolidated sediments tend to inhibit fluid circulation (Japan Trench40). The CA, the narrowest subduction system on Earth, hosts both processes: mud volcanism in the EL35,41,42, serpentinite diapirism in the WL (this study), as well as on the Hyblean Plateau close to the AEF onland43.
Mantle-derived serpentinites are found along slow to ultraslow spreading ridges, at magma-poor rifted margins and at subduction zones, where the downgoing plate may trigger arc magmatism or hydrate the deep Earth44,45,46. We suggest three possible scenarios for serpentinite emplacement in the Ionian Sea: during Pleistocene shear processes and fluid migration along the AEF; as a consequence of CA subduction; or during the formation of the Tethyan Ocean.
Diapirs alignment along the AEF suggests a direct link between tectonic activity and uprising of low-density material. However, we exclude that serpentinites formed in Pleistocene times since serpentinization is not limited to the fault zone, as indicated by velocity models of the lithosphere26 (Fig. 6). Tectonic dilatation in the rifting zone of the western Ionian Sea may trigger downwards fluid flow driving mobilization of serpentinites inherited from the Mesozoic ocean (Figs. 6 and 8) and already emplaced in the underplating African plate.
Schematic block diagram and cross section of the Western Calabrian Arc. a Block diagram of the Western Calabrian Arc subduction complex showing relationships between deep slab, transtensional faults, serpentinite diapirs, and Mt. Etna volcano (modified from refs. 62,69). Faulting associated with arc-perpendicular extensional tectonics in the accretionary wedge (modified from ref. 17) is shown in the yellow corridor between the Ionian and Alfeo-Etna fault systems. Extensional processes and associated Mt. Etna volcanism could be related to vertical upwelling of the asthenosphere at the SW lateral edge of the Ionian slab62,63,79. At the edges of a retreating slab the toroidal component of asthenospheric flow may induce an upward flow and decompression melting which could explain the formation of Mt. Etna2. Transtensional processes along the rifting zone between the Alfeo-Etna and Ionian fault may drive inherited serpentinite diapirism during slab tearing processes. b Sketch of the margin structure orthogonal to the arc based on interpretation of available multichannel seismic and Vp/Vs data26. Section D-D’ (location also indicated in Fig. 5) crosses the Calabrian Arc accretionary wedge, the inherited serpentinite diapir field, and Mt. Etna implying structural control on both diapirism and Mt. Etna formation. AEF, Alfeo-Etna fault; EL, eastern lobe of the accretionary wedge; IF, Ionian fault; OCT, Oceanic-continental transition; WL, western lobe of the accretionary wedge
A different scenario for serpentinite emplacement includes the mantle wedge at the subduction zone, possibly hydrated by the subducting plate. However, this scenario is unlikely in our region, because of the external position of the diapirs not far from the African foreland, in an area where the slab is rather shallow (Fig. 5). Velocity models for the Ionian lithosphere26, suggest two upper mantle layers (3.3 and 5 km thick, respectively, at 8 and 16 km depth), with high Vp/Vs ratio (1.95–2.1), interpreted as due to partly serpentinized peridotites26. The high Vp/Vs ratio suggests low pressure-low temperature serpentinization47, different from that inferred usually for subduction settings. Serpentinites may thus characterize the basement of the southern Ionian Sea in a region as wide as hundreds km2 at the transition with the African foreland.
If serpentinites are inherited from the Tethyan ocean, they are Mesozoic in age, and could have developed along fractures located near the boundary between the oceanic crust and the adjacent thinned continental crust, a situation similar to that observed west of the Galicia Bank along the Iberian margin24,48,49. Serpentinization of peridotites is due usually to hydrothermal processes and fluid circulation associated with early stages of accretion during rifting. As the sub- Ionian Sea crust has never been sampled and its nature never unambiguously defined, we cannot exclude that this part of the Ionian basin formed after tectonic extension and exhumation of a continental mantle causing serpentinization and fluid flow. In magma-poor margins, crustal faulting and mantle serpentinization are believed to occur when the crust becomes entirely brittle after a certain amount of extension50. In both scenarios, the AEF would represent a major Mesozoic tectonic discontinuity. Formation of serpentinites, regardless of geodynamic setting, is thus a consequence of lithospheric mantle exhumation as a result of a bulk pure shear, or related to Mesozoic rifting that brought the Moho and upper mantle rocks close to the surface. In any case, serpentinization could take place if suitable fluid pathways existed, with faults cutting down from the seafloor to the mantle44; the AEF may represent one such major discontinuity.
Interesting inferences can be obtained from xenoliths sampled from basalts and diatremes of the Hyblean Plateau (SE Sicily) East of the study region (Fig. 5). The Hyblean Plateau has been considered traditionally an uplifted portion of the Africa continental plate51,52. However, among the abundant lower crustal/upper mantle Hyblean xenoliths, none are of sub-continental derivation53. Hyblean xenoliths include serpentinized harzburgites with a relict primary spinel Cr number between 0.25 and 0.30; this and other geochemical indicators are within the range of oceanic peridotites53, suggesting that the lower crust may include bodies derived from an ancient oceanic domain54. Some Hyblean xenoliths consist of metasomatized gabbroic rocks of lower crustal origin probably derived from an oceanic or OCT environment54. Hydrothermal modification of xenoliths suggests that serpentinite-hosted hydrothermal systems54, typical of slow and ultra-slow spreading ridges, were active since the Middle Triassic, as deduced from U/Pb ages of hydrothermal zircons55. These results, backed by geophysics56,57, suggest a serpentinized oceanic-type lithosphere beneath the Hyblean Plateau and may imply some connection between sub Ionian Sea and sub Hyblean lithosphere.
Based on the P-wave velocity model proposed for southeastern Sicily56, petrological evidence for the lithosphere beneath the Hyblean Plateau and the occurrence of serpentinized harzburgite xenoliths in Hyblean diatremes were used to estimate the degree of serpentinization53. The resulting values are in line with the presence of peridotites affected by different degrees of serpentinization (35–100 vol.%) at 8–19 km depth, with the Moho discontinuity as a serpentinization front at about 19 km53 in agreement with results from the Ionian Sea26, further stressing a connection between the Ionian and Hyblean lithosphere.
Serpentinite diapirs below the Ionian Sea would have important implications to reconstruct the geometry and processes that have driven Tethyan ocean rifting. Our reconstruction, pointing towards mantle serpentinization coeval with Tethyan rifting, should be considered in the frame of Mesozoic paleogeography. Widespread serpentinites suggest that the Tethyan Ocean, at least in this domain, was not a magma-rich basin; serpentinization may have occurred along the OCT or may represent the product of hydration of exhumed mantle. Accordingly, the AEF may represent one of the shear zones where extension was focused, similar to the Magnaghi and Vavilov sub-basins of the Tyrrhenian Sea, where tensional faults favor mantle exhumation and percolation of water in the deep crust58. Alternatively, if the Tethyan basin was floored by oceanic crust, the AEF might represent the OCT where the crust was thin enough to become entirely brittle during rifting, as in the Porcupine Basin and Galicia margin24,44,49. Gravity data (Supplementary Fig. 13) provide evidence of OCT geometry. We suggest that gravity highs marking the AEF trend and basement topography below the accretionary wedge16 (Figs. 2 and 6) represent the transition between continental and oceanic or continental stretched domains (Fig. 5 and Supplementary Fig. 13). The complex geometry of the OCT reflects the 90° curvature of the continental margin close to the Sirte gulf (Fig. 1). In this reconstruction, oceanic crust, if present, should be limited to the inner accretionary wedge of the subduction complex and the AEF should be located close to the Mesozoic Tethyan OCT (Fig. 5).
Our reconstruction suggests that diapirism of inherited serpentinites occurs along active transverse shear zones segmenting the CA subduction complex, which thus controls structural development and fluid flow. Mt. Etna volcano is aligned with the serpentinite diapirs; its volcanism appears to be linked to active extensional tectonics, because the volcano is located on the NW tip of the AEF and has the same age of the inception of trans-tensional deformation in this sector17. Although Multi-Channel Seismic (MCS) data suggest that both the AEF and IF develop over structural boundaries inherited by the Mesozoic Tethyan basin, their recent tectonic activity fits with the latest Quaternary extensional phase. In fact, fault inception along the AEF produced the down throw of the WL and the formation of sedimentary basins, filled by up to 700–800 m-thick sediments (Fig. 3b s.p. 11,200). This observation was used to reconstruct the age of the fault assuming a sedimentation rate derived from sediment cores collected in the region18,59, indicating that 700 m-thick sediment along the AEF might correspond to 350,000–700,000 years, in agreement with the age of Mt. Etna volcano60. This connection between large-scale offshore tectonic processes and the formation of the volcano appears thus to be supported also by the age of transtensive reactivation along the AEF, which might indicate a primary role played in Mt. Etna volcanism by a Pleistocene geodynamic re-organization.
Etna Na-alkaline geochemistry61 does not indicate a deep slab magma source, as in the Aeolian K-alkaline magmatism of the southern Tyrrhenian. It may be related to a subduction rollback driven vertical ‘slab window’ that might explain its alkaline chemistry and the anomalous external position of the volcano, not derived from the mantle wedge but in some way related to the slab2,62,63. Moreover, recent Etnean basalts composition may reflect addition of fluids to the mantle source and/or selective crustal contamination12 occurring when the AEF shifted from compression with extension during Pleistocene17. This agrees with the voluminous melting under Mount Etna having originated from below the African plate62. Although it cannot be excluded, there is no geochemical evidence for a direct role of serpentinites in the genesis of Mt. Etna volcanism. Boron isotopes of Mt. Etna lavas do not support a role of serpentinite-derived fluids in their genesis64. However, a major tectonic discontinuity has been invoked as leading to the formation of Mt. Etna17,19,64. This fault would represent a window allowing upwelling of asthenosphere from below the subduction system between the Ionian and Sicilian slabs62 (Fig. 8). On the other hand, shallow partial melting as a response to mantle decompression was suggested as a source of Mt. Etna magmas within the framework of extensional tectonics65. Both scenarios (i.e., deep and shallow source) call for a major discontinuity, i.e., the AEF. Fluid addition along the fault coupled with a pressure drop might induce generation of magmas with the peculiar geochemistry of Etnean lavas64. The same process along the wide corridor between the Ionian and AEFs might have triggered offshore inherited serpentinite diapirism.
Serpentinite diapirs in the Ionian Sea provide evidence on how altered mantle rocks may be transferred from rifting, to basin floor, to subduction66 where they can give rise to large scale diapirism, enhancing shear processes and margin disruption during the final closure of the ocean. In this reconstruction, lithospheric inheritance controls neotectonic activity and location/reactivation of tectonic structures, as also observed in Canada67 or in the Eastern Mediterranean68. The transition from weak and stretched continental crust to the relatively strong oceanic crust appears to guide present day seismicity patterns as well as plate kinematic evolution of this region. It was suggested that serpentinites control the seismic or aseismic behavior of subduction zones38. Major historical earthquakes in the CA (i.e., the Messina 1908, Catania 1693, and Sicily 1169 earthquakes) were generated in an area close to the Ionian and AEFs18,69,70. Moreover, the distribution of seismicity and the lack of thrust-type earthquakes along the CA subduction complex in the Ionian Sea may be due to sub Ionian serpentinites influencing the mechanical coupling between plates as proposed for the Japanese margin71. Low thrust-type seismic activity may be an indication that subduction and convergence is aseismic, or that plate motion is expressed by large earthquakes with long recurrence time intervals. Serpentinites, when thicker than 2–3 km, may form an efficient flow channel at the top of the subducting slab influencing seismogenesis72; this has to be taken into account for seismic hazard assessments in this densely populated region.
The discovery of serpentinite diapirs in the CA provides the first example of a new class of serpentinites in a subduction system derived directly from the lower plate through lithospheric extensional faults disrupting the accretionary wedge and driven by tectonic processes at the slab edge. These findings may lead towards a more complete understanding of the structure of the Ionian lithosphere, of the mechanisms exposing altered products of inherited mantle peridotite at the seafloor long time after their formation, and finally of the role of serpentinization in driving continental collision.
Multibeam data. The morphobathymetric map is derived from three datasets collected by ISMAR-CNR during different cruises onboard of R/V CNR Urania, integrated by the 500 m grid DTM from CIESM/IFREMER Medimap group73. Swath bathymetry during the Urania cruises (2010, 2011, and 2015) was collected using a SIMRAD EM 302 Multibeam Echosounder. During those cruises, a keel-mounted probe provided real-time estimate for sound-velocity at the surface, whereas sound velocity profiles were performed daily to correct for density, salinity, and temperature variations in the water columns. Functioning of the multibeam echosounder was controlled by the Kongsberg software SIS, which provided data quality control and information for the survey plan. Multibeam data were subsequently processed using Caris Hips and Sips suite software.
Seismic reflection data
A close-spaced grid of seismic reflection profiles, characterized by different penetrations and vertical resolutions, was used to define the structural setting of the CA subduction complex in the Ionian Sea (Supplementary Fig. 1a). The deep structure was reconstructed from interpretation of MCS reflection profiles collected during different cruises (ENI-CA, Etnaseis, CALAMARE and Mediterranean Sea—MS data sets), from the 70s to the 90s. Single-channel Sparker profiles (the J data set of ISMAR-CNR) were used for detailing shallow structure of key areas previously identified through the analyses of deep penetrating seismic data. Acquisition and processing parameters of each set of data are described in Table 1.
Final data processing and visualization, as well as structural mapping was performed using SeisPrho software74, freely distributed at http://software.bo.ismar.cnr.it/seisprho.
We used all available seismic data (CNR_ENI Deep Crust Seismic Profiles—CROP, ENI CA MCS seismic, MCS CALAMARE lines, Sparker J data set, and Chirp profiles), characterized by different vertical resolutions (Table 1), to describe geometry and kinematics of the fault systems segmenting the western Ionian Sea.
Magnetic and gravimetric data
Shipborne magnetic data were acquired during CALAMARE08 cruise16 using a SeaSpy Marine Magnetics magnetometer, towed 180 m astern of the vessel. Raw magnetic data were collected at 1 Hz sampling frequency, resulting in an average spacing between every single measurement of about 4 meters.
Positioning of the magnetic sensor was obtained by a FUGRO differential Global Positioning System (GPS) mounted onboard of the R/V Urania after performing the layback correction. Despiking and low-pass filtering were applied to the raw data to attenuate high frequency noise. Magnetic anomalies were obtained removing the International Geomagnetic Reference Field (IGRF) geomagnetic reference field. Final magnetic anomalies (Fig. 2) were gridded using a minimum curvature method at 500 m grid cell size.
Gravimetric profiles were acquired along each CALAMARE08 seismic profile track using a LaCoste and Romberg AirSea Gravity System II at 1 Hz of sampling frequency. Gravity values were smoothed at 120 s by a low-pass filter, with an Exact Blackman Window. Drifting of the gravimeter was estimated at the beginning and end of the cruise; the average drift was below 0.2 mGal per day. Accuracy in positioning and velocity measures guaranteed 1 mGal accuracy in the gravity anomaly. Gravity meter readings were tied to an absolute gravity station in Naples, by using a portable model G gravity meter. Drift, Eötvös and Latitude corrections were applied to obtain the Free Air anomaly.
Quantitative interpretation of crustal properties below the CA-A seismic line track was obtained performing a 2.75D forward modeling, which allows geometry and density/susceptibility values of crustal bodies in strike direction to be defined. Magnetic and gravity data were sampled along the CA-A track using 400 data points (about 10 sample per km). The model was constrained by geometries and crustal features of major crustal bodies determined by seismostratigraphic interpretation of depth-migrated seismic sections. Seismic velocity of each stratigraphic interval was converted into density using Gardner’s equation75, which produced a preliminary set of values for subsequent re-adjustments. Best-fit between gravity and magnetic modeling was achieved by fine-tuning susceptibility and density values. After fixing the geometry of the diapiric body and the properties of the overburden, we tested the effect of different lithologies, including fine-grained clastic sediments, gabbroic cumulates, and a typical oceanic crust of slow-spreading domains. We finally observed that all realistic lithologies, alternative to that of serpentinized peridotites, provided different gravity and magnetic patterns relative to those observed in the field (Supplementary Fig. 8).
Pore water analyses
Pore water samples were extracted using Rhizones76 from cooled cores. Concentrations of chloride and sulphate were measured by ion chromatography (METROHM 761). Acidified subsamples were analyzed by inductively coupled plasma optical emission spectrometry (VARIAN 720-ES). Details on analytical methods applied routinely at GEOMAR are available at http://www.geomar.de or in publications31.
X-ray diffraction analysis of core GH15–03
The homogenized sediment sample A1 (Supplementary Figs. 10 and 11) was split and air-dried, and glycolized subsamples were analyzed by X-ray diffraction (Cu Kα, 3–70°, 0.02° steps). Mineral components and quantitative estimates were determined from the diffraction patterns, by using the software package XPowder.
Satellite-derived vertical gravity gradient
Vertical Gravity Gradient (VGG) data are from a global marine gravity model compilation77. CryoSat-2 and Jason-1 satellite missions provided high track coverage that enhanced the gravity anomaly signal within 12–40 wavelength band and a maximum resolution of 2 mGal. VGG computation permits to remove long wavelengths of gravity field particularly related to large-scale Moho-isostatic effects77. Peculiar tectonic structures and specific seafloor fabric are visible in the VGG data, which highlight density variation at a frequency band typical of upper/middle crust. The VGG model computed starting from satellite gravity data permits to distinguish OCT’s geometry68, as well as relict ridge structures below thick crust77
Data associated with this paper are available upon request (email@example.com).
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
1. Govers, R. & Wortel, M. J. R. Lithosphere tearing at STEP faults: response to edges of subduction zones. Earth Planet. Sci. Lett. 236, 505–523 (2005).
2. Faccenna, C. et al. Topography of the Calabria subduction zone (southern Italy): clues for the origin of Mt. Etna. Tectonics 30, TC1003 (2011).
3. Rosenbaum, G., Gasparon, M., Lucente, F. P., Peccerillo, A. & Miller, M. S. Kinematics of slab tear faults during subduction segmentation and implications for Italian magmatism. Tectonics 27, TC2008 (2008).
4. Schellart, W. P. Mount Etna-Iblean volcanism caused by rollback-induced upper mantle upwelling around the Ionian slab edge: an alternative to the plume model. Geology 38, 691–694 (2010).
5. Frizon de Lamotte, D. et al. The southernmost margin of the Tethys realm during the Mesozoic and Cenozoic: Initial geometry and timing of the inversion processes. Tectonics 30, TC3002 (2011).
6. de Voogd, B. et al. Two-ship deep seismic soundings in the basins of the Eastern Mediterranean Sea (Pasiphae cruise). Geophys. J. Int. 109, 536–552 (1992).
7. Catalano, R., Doglioni, C. & Merlini, S. On the Mesozoic Ionian Basin. Geophys. J. Int. 144, 49–64 (2001).
8. Speranza, F., Minelli, L., Pignatelli, A. & Chiappini, M. The Ionian Sea: the oldest in situ ocean fragment of the world? J. Geophys. Res. Solid Earth 117, B12101 (2012).
9. Makris, J., Nicolich, R. & Weigel, W. A seismic study in the western Ionian Sea. Ann. Geophys. Ser. B 4 B, 665–678 (1986).
10. Cowie, L. & Kusznir, N. Mapping crustal thickness and oceanic lithosphere distribution in the Eastern Mediterranean using gravity inversion. Pet. Geosci. 18, 373–380 (2012).
11. Cernobori, L., Nicolich, R., Petronio, L. & Romanelli, M. Crustal image of the Ionian basin and its Calabrian margins. Tectonophysics 264, 175–189 (1996).
12. Continisio, R., Ferrucci, F., Gaudiosi, G., Lo Bascio, D. & Ventura, G. Malta escarpment and Mt. Etna: early stages of an asimmetric rifting process? Acta Vulcanol. 9, 45–53 (1997).
13. Hieke, W., Hirschleber, H. B. & Dehghani, G. A. The Ionian Abyssal Plain (central Mediterranean Sea): Morphology, subbottom structures and geodynamic history—an inventory. Mar. Geophys. Res. 24, 279–310 (2003).
14. Malinverno, A. & Ryan, W. B. F. Extension in the Tyrrhenian Sea and shortening in the Apennines as result of arc migration driven by sinking of the lithosphere. Tectonics 5, 227–245 (1986).
15. Polonia, A. et al. The Calabrian Arc subduction complex in the Ionian Sea: regional architecture, active deformation, and seismic hazard. Tectonics 30, TC5018 (2011).
16. Bortoluzzi, G. et al. Styles and rates of deformation in the frontal accretionary wedge of the Calabrian Arc (Ionian Sea): controls exerted by the structure of the lower African plate. Ital. J. Geosci. 136, 347–364 (2017).
17. Polonia, A. et al. The Ionian and Alfeo–Etna fault zones: New segments of an evolving plate boundary in the central Mediterranean Sea? Tectonophysics 675, 69–90 (2016).
18. Polonia, A. et al. Turbidite paleoseismology in the Calabrian Arc Subduction Complex (Ionian Sea). Geochem. Geophys. Geosyst. 14, 112–140 (2013).
19. Hirn, A. et al. Roots of Etna volcano in faults of great earthquakes. Earth Planet. Sci. Lett. 148, 171–191 (1997).
20. Nicolich, R., Laigle, M., Hirn, A., Cernobori, L. & Gallart, J. Crustal structure of the ionian margin of Sicily: Etna volcano in the frame of regional evolution. Tectonophysics 329, 121–139 (2000).
21. Chamot-Rooke, N., Rangin, C., Le Pichon, X. & Group, D. W. DOTMED—deep offshore tectonics of the Mediterranean: a synthesis of deep marine data in eastern Mediterranean [CD-ROM]. Mem. Soc. Geol. Fr. 177, 64 (2005).
22. Mastrolembo Ventura, B. et al. Fast geodetic strain-rates in eastern Sicily (southern Italy): new insights into block tectonics and seismic potential in the area of the great 1693 earthquake. Earth Planet. Sci. Lett. 404, 77–88 (2014).
23. Della Vedova, B. & Pellis, G. New heat flow density measurements in the Ionian Sea. in 8° Meeting GNGTS 1133–1146 (GNGTS, 1989).
24. Bronner, A., Sauter, D., Manatschal, G., Péron-Pinvidic, G. & Munschy, M. Magmatic breakup as an explanation for magnetic anomalies at magma-poor rifted margins. Nat. Geosci. 4, 549–553 (2011).
25. Carlson, R. L. & Miller, D. J. A new assessment of the abundance of serpentinite in the oceanic crust. Geophys. Res. Lett. 24, 457–460 (1997).
26. D′Alessandro, A. et al. Evidence for serpentinization of the Ionian upper mantle from simultaneous inversion of P- and S-wave arrival times. J. Geodyn. 102, 115–120 (2016).
27. Toft, P. B., Arkani-Hamed, J. & Haggerty, S. E. The effects of serpentinization on density and magnetic susceptibility: a petrophysical model. Phys. Earth Planet. Inter. 65, 137–157 (1990).
28. Maffione, M., Morris, A., Plümper, O. & van Hinsbergen, D. J. J. Magnetic properties of variably serpentinized peridotites and their implication for the evolution of oceanic core complexes. Geochem. Geophys. Geosyst. 15, 923–944 (2014).
29. Miller, D. J. & Christensen, N. I. in Proceedings of the Ocean Drilling Program, 153 Scientific Results (Ocean Drilling Program, 1997). https://doi.org/10.2973/odp.proc.sr.153.043.1997.
30. Martin, J. B., Kastner, M., Henry, P., Le Pichon, X. & Lallement, S. Chemical and isotopic evidence for sources of fluids in a mud volcano field seaward of the Barbados accretionary wedge. J. Geophys. Res. Solid Earth 101, 20325–20345 (1996).
31. Scholz, F. et al. Isotopic evidence (87Sr/86Sr, δ7Li) for alteration of the oceanic crust at deep-rooted mud volcanoes in the Gulf of Cadiz, NE Atlantic Ocean. Geochim. Cosmochim. Acta 73, 5444–5459 (2009).
32. Hensen, C. et al. Strike-slip faults mediate the rise of crustal-derived fluids and mud volcanism in the deep sea. Geology 43, 339–342 (2015).
33. Mottl, M. J. J., Wheat, C. G. G., Fryer, P., Gharib, J. & Martin, J. B. B. Chemistry of springs across the Mariana forearc shows progressive devolatilization of the subducting plate. Geochim. Cosmochim. Acta 68, 4915–4933 (2004).
34. Hulme, S. M., Wheat, C. G., Fryer, P. & Mottl, M. J. Pore water chemistry of the Mariana serpentinite mud volcanoes: a window to the seismogenic zone. Geochem. Geophys. Geosyst. 11, Q01X09 (2010).
35. Panieri, G. et al. Mud volcanoes along the inner deformation front of the Calabrian Arc accretionary wedge (Ionian Sea). Mar. Geol. 336, 84–98 (2013).
36. Kopf, A., Klaeschen, D. & Mascle, J. Extreme efficiency of mud volcanism in dewatering accretionary prisms. Earth Planet. Sci. Lett. 189, 295–313 (2001).
37. Guillot, S., Schwartz, S., Reynard, B., Agard, P. & Prigent, C. Tectonic significance of serpentinites. Tectonophysics 646, 1–19 (2015).
38. Hyndman, R. D. & Peacock, S. M. Serpentinization of the forearc mantle. Earth. Planet Sci. Lett. 212, 417–432 (2003).
39. Kamimura, A. et al. Crustal structure study at the Izu-Bonin subduction zone around 31°N: implications of serpentinized materials along the subduction plate boundary. Phys. Earth Planet. Inter. 132, 105–129 (2002).
40. Taira, A. & Pickering, K. T. Sediment deformation and fluid activity in the Nankai, Izu-Bonin and Japan forearc slopes and trenches. Philos. Trans. R. Soc. Lond. A 335, 289–313 (1991).
41. Ceramicola, S., Praeg, D., Cova, A., Accettella, D. & Zecchin, M. Seafloor distribution and last glacial to postglacial activity of mud volcanoes on the Calabrian accretionary prism, Ionian Sea. Geo-Mar. Lett. 34, 111–129 (2014).
42. Capozzi, R. et al. Neogene to Quaternary tectonics and mud diapirism in the Gulf of Squillace (Crotone-Spartivento Basin, Calabrian Arc, Italy). Mar. Pet. Geol. 35, 219–234 (2012).
43. Manuella, F. C., Carbone, S. & Barreca, G. Origin of saponite-rich clays in a fossil serpentinite-hosted hydrothermal system in the crustal basement of the Hyblean Plateau (Sicily, Italy). Clays Clay. Miner. 60, 18–31 (2012).
44. Bayrakci, G. et al. Fault-controlled hydration of the upper mantle during continental rifting. Nat. Geosci. 9, 384–388 (2016).
45. Fryer, P. Serpentinite mud volcanism: observations, processes, and implications. Ann. Rev. Mar. Sci. 4, 345–373 (2012).
46. Deschamps, F., Godard, M., Guillot, S. & Hattori, K. Geochemistry of subduction zone serpentinites: a review. Lithos 178, 96–127 (2013).
47. Ji, S. et al. Seismic velocities, anisotropy, and shear-wave splitting of antigorite serpentinites and tectonic implications for subduction zones. J. Geophys. Res. Solid Earth 118, 1015–1037 (2013).
48. Boillot, G. et al. Ocean-continent boundary off the Iberian margin: a serpentinite diapir west of the Galicia Bank. Earth. Planet. Sci. Lett. 48, 23–34 (1980).
49. Dean, S. L., Sawyer, D. S. & Morgan, J. K. Galicia Bank ocean–continent transition zone: new seismic reflection constraints. Earth Planet. Sci. Lett. 413, 197–207 (2015).
50. Pérez-Gussinyé, M., Reston, T. J. & Phipps Morgan, J. Serpentinization and magmatism during extension at non-volcanic margins: the effect of initial lithospheric structure. Geol. Soc. Lond. Spec. Publ. 187, 551–576 (2001).
51. Catalano, R., Franchino, A., Merlini, S. & Sulli, A. A crustal section from the eastern Algerian basin to the Ionian ocean (Central Mediterranean). Mem. della Soc. Geol. Ital. 55, 71–85 (2000).
52. Scarascia, S., Lozej, A. & Cassinis, R. Crustal structures of the Ligurian, Tyrrhenian and Ionian Seas and adjacent onshore areas interpreted from wide-angle seismic profiles. Boll. di Geofis. Teor. Ed. Appl. 36, 5–19 (1994).
53. Manuella, F. C., Brancato, A., Carbone, S. & Gresta, S. A crustal–upper mantle model for southeastern Sicily (Italy) from the integration of petrologic and geophysical data. J. Geodyn. 66, 92–102 (2013).
54. Scribano, V., Sapienza, G., Braga, R. & Morten, L. Gabbroic xenoliths in tuff-breccia pipes from the Hyblean Plateau: insights into the nature and composition of the lower crust underneath South-eastern Sicily, Italy. Mineral. Petrol. 86, 63–88 (2006).
55. Sapienza, G. T., Griffin, W. L., O’Reilly, S. Y. & Morten, L. Crustal zircons and mantle sulfides: Archean to Triassic events in the lithosphere beneath south-eastern Sicily. Lithos 96, 503–523 (2007).
56. Brancato, A., Hole, J. A., Gresta, S. & Beale, J. N. Determination of seismogenic structures in southeastern Sicily (Italy) by high-precision relative relocation of microearthquakes. Bull. Seismol. Soc. Am. 99, 1921–1936 (2009).
57. Giampiccolo, E., Gresta, S. & Ganci, G. Attenuation of body waves in Southeastern Sicily (Italy). Phys. Earth Planet. Inter. 135, 267–279 (2003).
58. Prada, M., Ranero, C. R., Sallarès, V., Zitellini, N. & Grevemeyer, I. Mantle exhumation and sequence of magmatic events in the Magnaghi–Vavilov Basin (Central Tyrrhenian, Italy): new constraints from geological and geophysical observations. Tectonophysics 689, 133–142 (2016).
59. Polonia, A. et al. Mediterranean megaturbidite triggered by the AD 365 Crete earthquake and tsunami. Sci. Rep. 3, 1285 (2013).
60. Branca, S., Coltelli, M. & Groppelli, G. Geological evolution of a complex basaltic stratovolcano: Mount Etna, Italy. Ital. J. Geosci. 130, 306–317 (2011).
61. Peccerillo, A. & Frezzotti, M. L. Magmatism, mantle evolution and geodynamics at the converging plate margins of Italy. J. Geol. Soc. Lond. 172, 407–427 (2015).
62. Gvirtzman, Z. & Nur, A. The formation of Mount Etna as the consequence of slab rollback. Nature 401, 782–785 (1999).
63. Doglioni, C., Innocenti, F. & Mariotti, G. Why Mt Etna? Terra Nov. 13, 25–31 (2001).
64. Tonarini, S., Armienti, P., D’Orazio, M. & Innocenti, F. Subduction-like fluids in the genesis of Mt. Etna magmas: evidence from boron isotopes and fluid mobile elements. Earth Planet. Sci. Lett. 192, 471–483 (2001).
65. Viccaro, M. & Zuccarello, F. Mantle ingredients for making the fingerprint of Etna alkaline magmas: implications for shallow partial melting within the complex geodynamic framework of Eastern Sicily. J. Geodyn. 109, 10–23 (2017).
66. Bonatti, E. & Michael, P. J. Mantle peridotites from continental rifts to ocean basins to subduction zones. Earth Planet. Sci. Lett. 91, 297–311 (1989).
67. Audet, P., Sole, C. & Schaeffer, A. J. Control of lithospheric inheritance on neotectonic activity in northwestern Canada? Geology 44, 807–810 (2016).
68. Granot, R. Palaeozoic oceanic crust preserved beneath the eastern Mediterranean. Nat. Geosci. 9, 701–705 (2016).
69. Gutscher, M.-A. et al. Tectonic expression of an active slab tear from high-resolution seismic and bathymetric data offshore Sicily (Ionian Sea). Tectonics 35, 39–54 (2016).
70. Polonia, A. et al. A depositional model for seismo-turbidites in confined basins based on Ionian Sea deposits. Mar. Geol. 384, 177–198 (2017).
71. Seno, T. Variation of downdip limit of the seismogenic zone near the Japanese islands: implications for the serpentinization mechanism of the forearc mantle wedge. Earth Planet. Sci. Lett. 231, 249–262 (2005).
72. Hilairet, N. & Reynard, B. Stability and dynamics of serpentinite layer in subduction zone. Tectonophysics 465, 24–29 (2009).
73. Loubrieu, B. et al. Morpho-bathymetry of the Mediterranean Sea. (CIESM / Ifremer Medimap Group, 2008). https://doi.org/10.12770/5924fa10-ac45-11de-8853-000086f6a603.
74. Gasperini, L. & Stanghellini, G. SeisPrho: an interactive computer program for processing and interpretation of high-resolution seismic reflection profiles. Comput. Geosci. 35, 1497–1507 (2009).
75. Gardner, G. H. F., Gardner, L. W. & Gregory, A. R. Formation velocity and density – the diagnostic basics for stratigraphic traps. Geophysics 39, 770–780 (1974).
76. Seeberg-Elverfeldt, J., Schlüter, M., Feseker, T. & Kölling, M. Rhizon sampling of porewaters near the sediment-water interface of aquatic systems. Limnol. Oceanogr. Methods 3, 361–371 (2005).
77. Sandwell, D. T., Muller, R. D., Smith, W. H. F., Garcia, E. & Francis, R. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure. Science 346, 65–67 (2014).
78. Barreca, G. Geological and geophysical evidences for mud diapirism in south-eastern Sicily (Italy) and geodynamic implications. J. Geodyn. 82, 168–177 (2014).
79. Billi, A., Presti, D., Orecchio, B., Faccenna, C. & Neri, G. Incipient extension along the active convergent margin of Nubia in Sicily, Italy: Cefalù-Etna seismic zone. Tectonics 29, TC4026 (2010).
80. Wessel, P. & Smith, W. H. F. New, improved version of generic mapping tools released. Eos Trans. Am. Geophys. Union 79, 579–579 (1998).
We are greatly indebted to A.L. Cazzola, A. Fattorini, and C. Cattaneo (ENI Spa) for providing us with MCS lines CA-A, CA-B, and CA-C, to R. Nicolich for MCS line ETNASEIS-6, to OGS (Istituto nazionale di oceanografia e di geofisica sperimentale, Trieste) for MCS lines MS-26 and MS-107; M. Ligi for fruitful discussions and acquisition of seismic reflection line CALA-02. Cores “a” and “c” were collected during the “Seismofaults 2017” cruise onboard of R/V CNR Minerva-Uno (A. Billi, M. Cuffaro, and the Seismofaults 2017 team). We thank B. Domeyer, A. Bleyer, and R. Surberg for pore water analyses, Jutta Heinze for assisting XRD measurements, and Katja Lindhorst for pore water collection in core CQ14-02. We thank C. J. Garrido and an anonymous reviewer for their detailed and constructive reviews, which have contributed to improve the final quality of the manuscript. This work was supported by the project FLOWS (EU-COST action ES 1301). Morphobathymetric and structural maps were generated using the GMT software80. We wish to dedicate this work to the memory of Giovanni Bortoluzzi who suddenly died in September 2015. The acquisition of geophysical data in the Ionian Sea was possible thanks to his valuable contribution. This is ISMAR paper n. 1944.
Institute of Marine Sciences CNR ISMAR-Bo, Via Gobetti, 101, 40129, Bologna, Italy
A. Polonia, L. Gasperini, E. Bonatti & S. Romano
Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze, 157/A, Parma, Italy
L. Torelli, A. Artoni & M. Carlini
Istituto Nazionale di Geofisica e Vulcanologia, via di Vigna Murata 605, 00143, Rome, Italy
L. Cocchi & F. Muccini
Lamont-Doherty Earth Observatory of Columbia University, Palisades, NY, 10964, USA
GEOMAR Helmholtz-Zentrum für Ozeanforschung Kiel, Wischhofstrasse 1-3, 24148, Kiel, Germany
C. Hensen & M. Schmidt
A.P. conceived and designed the research. A.P., L.T., L.G., and S.R. acquired data at sea, and completed data processing and interpretation. L.C. and F.M. acquired potential field at sea and completed potential data processing and interpretation. E.B., C.H., and M.S. performed the mineralogical observations, PW analyses, and integrated data interpretation. A.A. and M.C. performed structural interpretations on seismic data. All authors discussed the problem, methods, analyses and results, and reviewed the manuscript.