Seismic investigation of the transition from continental to oceanic subduction along the western Hellenic Subduction Zone

Abstract

The western Hellenic subduction zone (WHSZ) exhibits well-documented along-strike variations in lithosphere density (i.e., oceanic versus continental), subduction rates, and overriding plate extension. Differences in slab density are believed to drive deformation rates along the WHSZ; however, this hypothesis has been difficult to test given the limited seismic constraints on the structure of the WHSZ, particularly beneath northern Greece. Here, we present high-resolution seismic images across northern and southern Greece to constrain the slab composition and mantle wedge geometry along the WHSZ. Data from two temporary arrays deployed across Greece in a northern line (NL) and southern line (SL) are processed using a 2D teleseismic migration algorithm based on the Generalized Radon Transform. Images of P- and S-wave velocity perturbations reveal N60E dipping low-velocity layers beneath both NL and SL. The ∼8 km thick layer beneath SL is interpreted as subducted oceanic crust while the ∼20 km thick layer beneath NL is interpreted as subducted continental crust. The thickness of subducted continental crust inferred within the upper mantle suggests that ∼10 km of continental crust has accreted to the overriding plate. The relative position of the two subducted crusts implies ∼70–85 km of additional slab retreat in the south relative to the north. Overall, our seismic images are consistent with the hypothesis that faster sinking of the denser, oceanic portion of the slab relative to the continental portion can explain the different rates of slab retreat and deformation in the overriding plate along the WHSZ.