Istituto Nazionale di Geofisica e Vulcanologia Osservatorio Etneo, Catania, Italy 2

Transcription

1 GNGTS 2016 Sessione 1.2 The Tethyan roots of southeastern Sicily (Italy): An interdisciplinary insight from new 3D seismic tomographies and petrophysical inferences E. Giampiccolo 1, A. Brancato 1, F.C. Manuella 2, S. Carbone 2, S. Gresta 2, V. Scribano 2 1 Istituto Nazionale di Geofisica e Vulcanologia Osservatorio Etneo, Catania, Italy 2 Dip. Scienze Biologiche, Geologiche e Ambientali, Sezione di Scienze della Terra, Università di Catania, Italy A number of geological and geophysical investigations were carried out to determine the structure and composition of lithosphere of the Hyblean Plateau (south-eastern Sicily) (e.g., Finetti et al., 2005 and references therein). Detailed datasets were gathered for the Hyblean upper crust, consisting of a Meso-Cenozoic sedimentary-volcanic sequence (Bianchi et al., 1987). However, no direct information is known for the buried Permo-Triassic basement of this region (Bianchi et al., 1987; Finetti et al., 2005), which is believed the uplifted emerged portion of the Pelagian block, as the northern offshoot of the Africa continental plate (Burollet et al., 1978, Ben-Avraham et al., 1990; Lentini et al., 1994), extending from Tunisia to Sicily. On the other hand, a 30-year study of deep-seated xenoliths embedded in tuff-breccia deposits of Hyblean diatremes (e.g., Scribano 1987; Scribano et al., 2006a,b), Late Tortonian-Lower Messinian in age (Suiting and Schmincke, 2010), led Manuella et al. (2013, 2014, 2015) to propose a new lithospheric model for southeastern Sicily, and neighboring areas. This model, which is based on an interdisciplinary study of petrologic and geochemical data from Hyblean xenoliths and volcanic rocks, with recent geophysical datasets (Giampiccolo et al., 2003; Brancato, 2005), is consistent with the hypothesis proposed by Scribano et al. (2006a,b) and the palaeogeographic model developed by Vai (1994, 2003) and Catalano et al. (1996), attesting that Hyblean Permo- Triassic basement, is a remnant of the slow-spreading Permian Tethys Ocean (Manuella et al., 2015a) as a part of the Oman-Iraq-Levantine-Sicily-Texas oceanic seaway. 210

2 GNGTS 2016 Sessione 1.2 Geophysical investigations of the crustal structure of the Hyblean Plateau include several velocity tomography studies, in which different techniques and datasets were considered (e.g. Di Stefano et al., 1999; Barberi et al., 2004; Scarfì et al., 2007; Brancato et al., 2009; Musumeci et al., 2014). Fig. 1 P-waves velocity model resulting from the 3D inversion. S-N and W-E vertical sections along the traces reported at 0 km level are also shown. On the 0 km layer the main structural features are reported with grey lines. The zones with RDE 0.4 SF 3.5; DWS 100 are shown with thick black, dotted white and dashed white contour lines, respectively. 211

3 GNGTS 2016 Sessione 1.2 Fig. 2 VP/VS model resulting from the 3D inversion. S-N and W-E vertical sections along the traces reported at 0 km level are also shown. The zones with RDE 0.4 SF 3.5; DWS 100 are shown with thick black, dotted white and dashed white contour lines, respectively. In the present work we propose new 3D VP and VP/VS models and, for the first time, a 3D P-wave attenuation image (QP) for southeastern Sicily, obtained by using a large dataset of local seismic events that occurred in the time span Measurements of attenuation of seismic waves are important indicators of Earth heterogeneities, not easily amenable to study using only seismic velocities (i.e., Fehler et al., 1992; Mitchell, 1995). We used the SIMULPS12 (Thurber, 1983). algorithm to jointly invert a total of 4448 P and 2876 S absolute arrival times. Based on the P and S ray paths of the selected data obtained by 212

4 GNGTS 2016 Sessione 1.2 Fig. 3 QP model resulting from the 3D inversion. S-N and W-E vertical sections along the traces reported at 0 km level are also shown. The zones with RDE 0.35 SF 4.5 and DWS > 100 shown with thick black, dotted white and dashed white contour lines, respectively. the pseudo-bending method (Um and Thurber, 1987), we used a regular horizontal grid with 10 km 10 km node spacing, covering an area of 110 km 110 km. In the vertical directions we used a grid spacing of 2 km, covering a depth range from the surface to 24 km, as in Brancato et al. (2009). A minimum 1D velocity model was used as input for the 3D inversion, with the initial VP/VS fixed at a constant value of 1.78 at all nodes, based on previous studies in the investigated area (Musumeci et al., 2003; Piana Agostinetti and Amato, 2009). Figgs. 1 and 2 show the horizontal slices at different depths and the N-S and W-E cross 213

5 GNGTS 2016 Sessione 1.2 sections of the V P, and V P models obtained, respectively. In the obtained V P image (Fig. 1) well resolved areas (RDE 0.4 SF 3.5; DWS 100) are mostly located in the center of the model. From the surface down to 4 km, we observe a central low V P volume that well correlates with the known geology consisting of a package of Meso-Cenozoic carbonate rocks and thinner marl levels (Bianchi et al., 1987). From 6 to about 16 km, the most significant feature is a central eastern high V P anomaly that occurs with a remarkable variation from 4.8 to 6.9 km/s, showing a dome-shaped basement (see cross sections in Fig. 1). As for V P, the V P model (Fig. 2), is well resolved (RDE 0.4 SF 3.5 DWS 100) only in the central eastern part of the investigated region. From the surface down to 4 km depth, we notice clear V P contrasts well matching with V P anomalies. From about 6 km to 19 km, the central eastern area shows high V P overlapping the high V P region. The attenuation along a ray path is quantified by the t* operator which is evaluated, following Eberhart-Phillips and Chadwick (2002), by fitting the amplitude spectrum decay of P-waves (e.g. De Gori et al., 2005). A total of 2596 t* observations were inverted for the 3D Q P structure, by using the SIMULPS12 algorithm modified for attenuation by Rietbrock (2001). The medium was parameterized with the same 3D grid nodes and velocity values obtained by the 3D V P tomography, assuming an initial Q P = 300, which is reported in literature as the average value in the crust of south-eastern Sicily (e.g., Giampiccolo et al., 2003; de Lorenzo et al., 2004). Well resolved areas, delimited by RDE 0.35 SF 4.5 and DWS > 100, are located in the central eastern sector of the study area (Fig. 3). The most prominent feature is a high Q P (> 300) dome-shaped volume that extends from 6 km to about 20 km. Above 6 km, a very small central portion is characterized by low Q P that correlates with low V P and high V P. At greater depths, in the eastern sector, high Q P values are mainly found in coincidence with high V P and low V P /Vs volumes. Conversely, in the central area high Q P values are associated with high V P and high V P /Vs regions (Figs. 1 and 2). The joint interpretation of both tomographies and the petrophysical interpretations of velocity and attenuation anomalies provide some constraints for the nature and composition of the lithosphere beneath the region. In particular, the present tomographies allow us to extend the punctual outlook deriving from Hyblean xenoliths found in diatremes, representing natural inverse drillings (from the bottom to the surface), to the whole Hyblean lithosphere. Moreover, although the oceanic nature of the Hyblean lithosphere is currently strongly debated, new density and degree of serpentinization models retrieved here from V P values, suggest that the structure and composition of the Hyblean lithosphere may differ from the traditional Africa continental plate. References Barberi G., Cosentino M.T., Gervasi A., Guerra I., Neri G. and Orecchio B.; 2004: Crustal seismic tomography in the Calabrian Arc region, south Italy, Phys. Earth Planet. Int., 147, Ben-Avraham Z., Boccaletti M., Cello G., Grasso M., Lentini F., Torelli L. and Tortorici L.; 1990: Principali domini strutturali originatisi dalla collisione nogenico-quaternaria nel Mediterraneo centrale, Mem. Soc. Geol. It., 45, Bianchi F., Carbone S., Grasso M., Invernizzi G., Lentini F., Longaretti G., Merlini S. and Mostardini F.; 1987: Sicilia orientale: Profilo geologico Nebrodi-Iblei, Mem. Soc. Geol. It., 38, Brancato A.; 2005: Reconstruction of seismogenic structures by high-precise relative location of microearthquake hypocenters occurred in southeastern Sicily (Italy). Ph.D. thesis, Department of Earth Science, University of Naples Federico II, Naples, Italy. Brancato A., Hole J.A., Gresta S. and Beale J.N.; 2009: Determination of seismogenic structures in southeastern Sicily (Italy) by high-precision relative relocation of microearthquakes, Bull. Seismol. Soc. Amer., 99, Burollet F.P., Mugniot J.M. and Sweeney P.; 1978: The geology of the Pelagian Block: the Margins and Basins of Southern Tunisia and Tripolitania, in The Ocean Basins and Margins, edited by A.E.M. Nairm, W.H. Kanes, and F.G. Stehli, pp , The Western Mediterranean. Plenum Press, New York. Catalano R., Di Stefano P., Sulli A. and Vitale F.P.; 1996: Paleogeography and structure of the central Mediterranean: Sicily and its offshore area, Tectonophysics, 260,

6 GNGTS 2016 Sessione 1.2 De Gori P., Chiarabba C. and Patanè D.; 2005: Qp structure of Mount Etna: Constraints for the physics of the plumbing system, J. Geophys. Res.: Solid Earth, 110, B5, doi: /2003JB de Lorenzo S., Di Grazia G., Giampiccolo E., Gresta S., Langer H., Tusa G. and Ursino A.; 2004: Source and Q P parameters from pulse width inversion of microearthquake data in southeastern Sicily, Italy. J. Geophys. Res.: Solid Earth, 109, doi: /2003JB Di Stefano R., Chiarabba C., Lucente F. and Amato A.; 1999: Crustal and uppermost mantle structure in Italy from the inversion of P-wave arrival times: geodynamic implications, Geophys. J. Int., 139, Eberhart-Phillips D. and Chadwick M.; 2002: Three-dimensional attenuation model of the shallow Hikurangi subduction zone in the Raukumara Peninsula, New Zealand, J. Geophys. Res., 107, B2, doi: /2000JB Fehler M., Hoshiba M., Sato H. and Obara K.; 1992: Separation of scattering and intrinsic attenuation for the Kanto- Tokai region, using measurements of S-wave energy versus hypocentral distance, Geophys. J. Int., 108, Finetti I., Lentini F., Carbone S., Del Ben A., Di Stefano A., Forlin E., Guarnieri P., Pipan M. and Prizzon A.; 2005: Geological outline of Sicily and lithospheric tectono-dynamics of its Tyrrhenian margin from new CROP seismic data, in Seismic Exploration of the Mediterranean region, CROP Deep, edited by I.R., Chapter 15, pp , Elsevier. Giampiccolo E., Gresta S. and Ganci G.; 2003: Attenuation of body waves in Southeastern Sicily (Italy), Phys. Earth Planet. Int., 135, Lentini F., Carbone S. and Catalano S.; 1994: Main structural domains of the Central Mediterranean region and their Neogene tectonic evolution, Boll, Geofis, Teor, Appl,, 36, Manuella F.C., Brancato A., Carbone C. and Gresta S.; 2013: A crustal-upper mantle model for southeastern Sicily (Italy) from the integration of petrologic and geophysical data, J. Geodyn., 66, Manuella F.C., Scribano V., Carbone S. and Brancato A.; 2015: The Hyblean xenolith suite (Sicily): an unexpected legacy of the Ionian Tethys realm, Int. J. Earth Sci., 104, Manuella F.C., Brancato A., Carbone S. and Gresta S.; 2014: Reply to Comments on the paper A crustal-upper mantle model for southeastern Sicily (Italy) from the integration of Mitchell B.J.; 1995: Anelastic structure and evolution of the continental crust and upper mantle from seismic surface wave attenuation, Rev. Geophys., 33, Musumeci C., Di Grazia G. and Gresta S.; 2003: Minimum 1-D velocity model in Southeastern Sicily (Italy) from local earthquake data: an improvement in location accuracy, J. Seismol., 7, Musumeci C., Scarfì L., Palano M. and Patanè D.; 2014: Foreland segmentation along an active convergent margin: New constraints in southeastern Sicily (Italy) from seismic and geodetic observations, Tectonophysics, 630, Piana Agostinetti N. and Amato A.; 2009: Moho depth and Vp/Vs ratio in peninsular Italy from teleseismic receiver functions, J. Geophys. Res., 114, B06303, doi: /2008jb Rietbrock A.; 2001: P wave attenuation structure in the fault area of the 1995 Kobe earthquake, J. Geophys. Res.: Solid Earth, 106, doi: /2000JB Scarfì L., Giampiccolo E., Musumeci C., Patanè D. and Zhang H.; 2007: New insights on 3D crustal structure in southeastern Sicily (Italy) and tectonic implications from an adaptive mesh seismic tomography, Phys. Earth Planet. Int., 161, Scribano V.; l987: Origin of websterite nodules from some alkaline volcanic rocks of Hyblean Plateau (South-Eastern Sicily), Per. Mineral., 56, Scribano V., Ioppolo S. and Censi P.; 2006: Chlorite/smectite-alkali feldspar metasomatica xenoliths from Hyblean Miocenic diatremes (Sicily, Italy): evidence for early interaction between hydrothermal brines and ultramafic/ mafic rocks at crustal levels, Ofioliti, 31, Scribano V., Sapienza G.T., Braga R. and Morten L.; 2006b: Gabbroic xenoliths in tuff-breccia pipes from the Hyblean Plateau: insights into the nature and composition of the lower crust underneath Southeastern Sicily, Italy, Mineral. Petrol., 86, Suiting I., Schmincke H.U.; 2010: Iblean diatremes 2: shallow marine volcanism in the Central Mediterranean at the onset of the Messinian Salinity Crisis (Iblean Mountains, SE-Sicily) a multidisciplinary approach, Int. J. Earth Sci., 99, Thurber C.H.; 1983: Earthquake locations and three-dimensional crustal structure in the Coyote Lake area, central California, J. Geophys. Res.: Solid Earth, 88, Um J., and Thurber C.; 1987: A fast algorithm for two-point seismic ray tracing, Bull. Seismol. Soc. Am., 77, Vai G.B.; 1994: Crustal evolution and basement elements in the Italian area: palaeogeography and characterization, Boll. Geofis. Teor. Appl., 36, Vai G.B.; 2003: Development of the palaeogeography of Pangaea from Late Carboniferous to Early Permian, Palaeogeogr. Palaeocl. Palaeoec., 196,