High-resolution seismic imaging of the western Hellenic subduction zone using teleseismic scattered waves

SUMMARY The active Hellenic subduction system has long been considered an ideal setting for studying subduction dynamics because it is easily accessible and of limited spatial extent. It has been the focus of numerous seismological studies over the last few decades but, nonetheless, the detailed structure of both the slab and the surrounding mantle remain poorly constrained in an intermediate depth range from 30 to 150 km. The objective of this paper is to fill this gap. The intermediate depth regime is of particular interest because it is pivotal for improving our understanding of the dynamic interaction between subducting lithosphere and the surrounding mantle. An interdisciplinary effort aimed at addressing this challenge is currently undertaken by the ‘Multidisciplinary Experiments for Dynamic Understanding of Subduction under the Aegean Sea’ (MEDUSA) project. As part of the MEDUSA initiative, a temporary array consisting of 40 densely spaced broad-band seismometers from the IRIS-PASSCAL pool has been deployed in southern Greece. We process the teleseismic data recorded by this array with a migration algorithm based on the generalized radon transform to obtain high-resolution images of the subduction zone in 2-D. The images reveal a sharp Mohoroviy´ discontinuity (Moho) at depths ranging from 30 km beneath the western margin of the Aegean Sea to 40 km beneath the central Peloponnesus, where it outlines the crustal root of the Hellenides. To the west of the Hellenides, the continental Moho is not identified, but we interpret a pronounced discontinuity imaged at ∼20 km depth as the contact between low-velocity sediments and high-velocity crystalline basement. The images also show the subducted oceanic crust as a low-velocity layer that plunges at a constant angle of 21 ◦ from west to east. The oceanic crust exhibits low velocities to at least 90 km depth, indicating that the bulk of fluid transfer from the subducted slab into the mantle wedge occurs below this depth. A detailed comparison of images constructed for distinct backazimuthal illuminations reveals deviations in the geometry of the subducted slab. These deviations are attributed to structural and/or compositional changes taking place directly to the north of the MEDUSA array, and are consistent with the existence of a slab tear beneath the Central Hellenic Shear Zone.

[1]  C. Papazachos,et al.  The geometry of the Wadati–Benioff zone and lithospheric kinematics in the Hellenic arc , 2000 .

[2]  H. Kárason.  Constraints on mantle convection from seismic tomography and flow modeling , 2002 .

[3]  A. Hirn,et al.  A traverse of the ionian islands front with coincident normal incidence and wide-angle seismics , 1996 .

[4]  P. Charvis,et al.  Reflection-refraction seismics in the Gulf of Corinth: hints at deep structure and control of the deep marine basin , 2004 .

[5]  Gregory Beylkin,et al.  A new slant on seismic imaging: Migration and integral geometry , 1987 .

[6]  R. Armijo,et al.  Quaternary evolution of the Corinth Rift and its implications for the Late Cenozoic evolution of the Aegean , 1996 .

[7]  Fabrizio Innocenti,et al.  Tertiary to Quaternary evolution of volcanism in the Aegean region , 1984, Geological Society, London, Special Publications.

[8]  J. Canas,et al.  Objective regionalization of Rayleigh wave dispersion data by clustering algorithms: an application to the Mediterranean basin , 2001 .

[9]  Demitris Paradissis,et al.  Global Positioning System constraints on plate kinematics and dynamics in the eastern Mediterranean and Caucasus , 2000 .

[10]  B. Zelt,et al.  Crustal velocity and Moho structure beneath the Gulf of Corinth, Greece , 2005 .

[11]  W. Spakman,et al.  Subduction and slab detachment in the Mediterranean-Carpathian region. , 2000, Science.

[12]  S. Karato.  Mapping water content in the upper mantle , 2013 .

[13]  Robert D. van der Hilst,et al.  Travel-time tomography of the European-Mediterranean mantle down to 1400 km , 1993 .

[14]  M. Pasyanos,et al.  Crustal and upper-mantle structure in the Eastern Mediterranean from the analysis of surface wave dispersion curves , 2007 .

[15]  W. Walter,et al.  Crust and upper‐mantle structure of North Africa, Europe and the Middle East from inversion of surface waves , 2002 .

[16]  Gregory Beylkin,et al.  Linearized inverse scattering problems in acoustics and elasticity , 1990 .

[17]  J. Peterson,et al.  Observations and modeling of seismic background noise , 1993 .

[18]  Matthias Hort,et al.  Serpentine and the subduction zone water cycle , 2004 .

[19]  D. Giardini,et al.  Correlated shear and bulk moduli to 1400 km beneath the Mediterranean region , 2006 .

[20]  C. Haberland,et al.  Complex patterns of fluid and melt transport in the central Andean subduction zone revealed by attenuation tomography , 2003 .

[21]  H. Kawakatsu,et al.  Seismic Evidence for Deep-Water Transportation in the Mantle , 2007, Science.

[22]  B. Endrun,et al.  Comparison of gravimetric and seismic constraints on the structure of the Aegean lithosphere in the forearc of the Hellenic subduction zone in the area of Crete , 2007 .

[23]  B. Endrun,et al.  S velocity structure and radial anisotropy in the Aegean region from surface wave dispersion , 2008 .

[24]  A. Serpetsidaki,et al.  Moho topography under central Greece and its compensation by Pn time-terms for the accurate location of hypocenters: The example of the Gulf of Corinth 1995 Aigion earthquake , 2007 .

[25]  Gregory C. Beroza,et al.  High-resolution subducting-slab structure beneath northern Honshu, Japan, revealed by double-difference tomography , 2004 .

[26]  M. Meijde,et al.  Receiver Function Study of the Hellenic Subduction Zone , 2003 .

[27]  L. Royden.  The tectonic expression slab pull at continental convergent boundaries , 1993 .

[28]  A. Kopf,et al.  The Mediterranean Ridge: A mass balance across the fastest growing accretionary complex on Earth , 2003 .

[29]  E. Chaumillon,et al.  An overview of Mediterranean Ridge collisional accretionary complex as deduced from multichannel seismic data , 1998 .

[30]  Michael G. Bostock,et al.  Seismic evidence for overpressured subducted oceanic crust and megathrust fault sealing , 2009, Nature.

[31]  Hans-Peter Harjes,et al.  Imaging crustal discontinuities and the downgoing slab beneath western Crete , 2000 .

[32]  R. O. Hansen,et al.  Study of the crustal thickness and the subducting lithosphere in Greece from gravity data , 1997 .

[33]  D. McKenzie.  Active tectonics of the Alpine—Himalayan belt: the Aegean Sea and surrounding regions , 1978 .

[34]  Cenka Christova.  Seismicity depth pattern, seismic energy and b value depth variation in the Hellenic Wadati-Benioff zone , 1992 .

[35]  H. Lyon‐Caen,et al.  Moho topography beneath the Corinth Rift area (Greece) from inversion of gravity data , 2001 .

[36]  H. Pedersen,et al.  Constraints on absolute S velocities beneath the Aegean Sea from surface wave analysis , 2005 .

[37]  P. Neumann,et al.  The Cretaceous sedimentary history of the Pindos Basin (Greece) , 2004 .

[38]  A. Hosa.  Imaging of the Hellenic subduction zone by seismic tomography , 2008 .

[39]  Simon M. Peacock,et al.  Subduction factory 1. Theoretical mineralogy, densities, seismic wave speeds, and H 2 O contents , 2003 .

[40]  G. Abers,et al.  Unusual mantle Poisson's ratio, subduction, and crustal structure in central Alaska , 2006 .

[41]  J. Dewey,et al.  Aegean and surrounding regions: Complex multiplate and continuum tectonics in a convergent zone , 1979 .

[42]  W. Spakman,et al.  The Hellenic Subduction Zone: A tomographic image and its geodynamic implications , 1988 .

[43]  A. Vafidis,et al.  Deep Seismic Imaging and Velocity Estimation in Ionian Sea , 2003 .

[44]  Xavier Le Pichon,et al.  Constraints on the evolution and vertical coherency of deformation in the Northern Aegean from a comparison of geodetic, geologic and seismologic data , 2004 .

[45]  S. Karato,et al.  Dynamic recrystallization of olivine single crystals during high‐temperature creep , 1980 .

[46]  D. McKenzie.  Active Tectonics of the Mediterranean Region , 1972 .

[47]  X. Le Pichon,et al.  The Aegean Sea , 1981, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[48]  H. Lyon‐Caen,et al.  Crustal and upper mantle structure beneath the Corinth rift (Greece) from a teleseismic tomography study , 2000 .

[49]  T. Hearn.  Uppermost mantle velocities and anisotropy beneath Europe , 1999 .

[50]  J. Makris.  The crust and upper mantle of the Aegean region from deep seismic soundings , 1978 .

[51]  G. Abers,et al.  High resolution image of the subducted Pacific (?) plate beneath central Alaska, 50–150 km depth , 2003 .

[52]  G. Panza,et al.  Shear velocity structure in the Aegean area obtained by inversion of Rayleigh waves , 2003 .

[53]  A. Sheehan,et al.  Mantle discontinuity structure from midpoint stacks of converted P to S waves across the Yellowstone hotspot track , 1997 .

[54]  J. Gerald,et al.  Relationships between dynamically recrystallized grain size and deformation conditions in experimentally deformed olivine rocks , 1993 .

[55]  Daniel E. McNamara,et al.  Seismic Noise Analysis System, Power Spectral Density Probability Density Function: Stand-Alone Software Package , 2005 .

[56]  K. Priestley,et al.  Lithospheric structure of the Aegean obtained from P and S receiver functions , 2006 .

[57]  Michael G. Bostock,et al.  Multiparameter two-dimensional inversion of scattered teleseismic body waves 3. Application to the Cascadia 1993 data set , 2001 .

[58]  G. Nolet,et al.  P and S deep velocity structure of the Hellenic area obtained by robust nonlinear inversion of travel times , 1997 .

[59]  G. Panza,et al.  The lithosphere in the central-eastern Mediterranean area , 1982 .

[60]  P. Davy,et al.  Periodic instabilities during compression or extension of the lithosphere 1. Deformation modes from an analytical perturbation method , 1992 .

[61]  X. Pichon,et al.  The hellenic arc and trench system: A key to the neotectonic evolution of the eastern mediterranean area , 1979 .

[62]  F. Marone,et al.  Joint inversion of local, regional and teleseismic data for crustal thickness in the Eurasia–Africa plate boundary region , 2003 .

[63]  D. Papanikolaou,et al.  Active faults, deformation rates and Quaternary paleogeography at Kyparissiakos Gulf (SW Greece) deduced from onshore and offshore data , 2007 .

[64]  J. Gerald,et al.  Grain-size-sensitive seismic wave attenuation in polycrystalline olivine , 2002 .

[65]  Brian J. Mitchell,et al.  Crustal thickness variations in the Aegean region and implications for the extension of continental crust , 2005 .

[66]  G. Abers,et al.  Subduction Factory 3: An Excel worksheet and macro for calculating the densities, seismic wave speeds, and H2O contents of minerals and rocks at pressure and temperature , 2004 .

[67]  B. Endrun,et al.  Modeling the influence of Moho topography on receiver functions: A case study from the central Hellenic subduction zone , 2005 .

[68]  Peter J. Clarke,et al.  Crustal strain in central Greece from repeated GPS measurements in the interval 1989–1997 , 1998 .

[69]  Mark van der Meijde,et al.  Crustal structure beneath broad-band seismic stations in the Mediterranean region , 2003 .

[70]  I. Main,et al.  3-D structure of the lithosphere in the Aegean region , 1990 .

[71]  E. Carminati,et al.  The role of slab detachment processes in the opening of the western–central Mediterranean basins: some geological and geophysical evidence , 1998 .

[72]  G. Veis,et al.  ESTABLISHMENT OF A CONTINUOUS GPS NETWORK ACROSS THE KEPHALONIA FAULT ZONE, IONIAN ISLANDS, GREECE , 1998 .

[73]  Michael G. Bostock,et al.  Multiparameter two‐dimensional inversion of scattered teleseismic body waves 2. Numerical examples , 2001 .

[74]  P. Charvis,et al.  Seismic structure and the active Hellenic subduction in the Ionian islands , 2000 .

[75]  G. Abers,et al.  Seismic imaging of subduction zone metamorphism , 2008 .

[76]  Andrea Morelli,et al.  P wave tomography of the mantle under the Alpine-Mediterranean area , 2003 .

[77]  Gregory Beylkin,et al.  Imaging of discontinuities in the inverse scattering problem by inversion of a causal generalized Radon transform , 1985 .

[78]  B. Endrun,et al.  Seismicity of the Hellenic subduction zone in the area of western and central Crete observed by temporary local seismic networks , 2004 .

[79]  Correlated shear and bulk moduli to 1400 km beneath the Mediterranean region , 2006 .

[80]  Fabrizio Innocenti,et al.  Geochronological data on recent magmatism of the Aegean Sea , 1976 .

[81]  M. Bohnhoff,et al.  Crustal investigation of the Hellenic subduction zone using wide aperture seismic data , 2001 .

[82]  D. McKenzie,et al.  Plate Tectonics of the Mediterranean Region , 1970, Nature.

[83]  A. Robertson,et al.  Sedimentary geology as a key to understanding the tectonic evolution of the Mesozoic–Early Tertiary Paikon Massif, Vardar suture zone, N Greece , 2003 .

[84]  A. Vafidis,et al.  DEEP SEISMIC REFLECTION DATA FROM OFFSHORE WESTERN GREECE: A NEW CRUSTAL MODEL FOR THE IONIAN SEA , 2005 .

[85]  E. R. Engdahl,et al.  Evidence for deep mantle circulation from global tomography , 1997, Nature.

[86]  Daniel E. McNamara,et al.  Ambient Noise Levels in the Continental United States , 2004 .

[87]  J. Brun,et al.  Aegean crustal thickness inferred from gravity inversion. Geodynamical implications , 2004 .

[88]  K. Priestley,et al.  Variations in the crustal structure beneath western Turkey , 1998 .

[89]  R. Twiss.  Theory and applicability of a recrystallized grain size paleopiezometer , 1977 .

[90]  Michael G. Bostock,et al.  Multiparameter two-dimensional inversion of scattered teleseismic body waves 1. Theory for oblique incidence , 2001 .

[91]  G. Veis,et al.  The Kephalonia Transform Fault and the rotation of the Apulian Platform: Evidence from satellite geodesy , 1993 .

[92]  D. Giardini,et al.  Teleseismic delay times and shear-wave splitting in the mediterranean region , 2004 .