We use structural and microstructural observations from exhumed subduction-related rocks exposed on Syros Island (Cyclades, Greece) to provide constraints on the length scales and types of heterogeneities that occupy the deep subduction interface, with possible implications for episodic tremor and slow slip. We selected three Syros localities that represent different oceanic protoliths and deformation conditions within a subduction interface shear zone, including: (1) prograde subduction of oceanic crust to eclogite facies; (2) exhumation of oceanic crust from eclogite through blueschist-greenschist facies; and (3) exhumation of mixed mafic crust and sediments from eclogite through blueschist-greenschist facies. All three localities preserve rheological heterogeneities that reflect metamorphism of primary lithological, geochemical, and/or textural variations in the subducted protoliths and that take the form of brittle pods and lenses within a viscous matrix. Microstructural observations indicate that the matrix lithologies (blueschists and quartz-rich metasediments) deformed by distributed power-law viscous flow accommodated by dislocation creep in multiple mineral phases. We estimate bulk shear zone viscosities ranging from ∼1018 to 1020 Pa-s, depending on the relative proportion of sediments to (partially eclogitized) oceanic crust. Eclogite and coarse-grained blueschist heterogeneities within the matrix preserve multiple generations of dilational shear fractures and veins formed under high-pressure conditions. The veins commonly show coeval or overprinting viscous shear, suggesting repeated cycles of frictional and viscous strain. These geologic observations are consistent with a mechanical model of episodic tremor and slow slip (ETS), in which the deep subduction interface is a rheologically heterogeneous distributed shear zone comprising transiently brittle (potentially tremor-genic) sub-patches within a larger, viscously creeping interface patch. Based on our observations of outcrop and map areas of heterogeneous patches and the sizes, distributions, and amounts of brittle offset recorded by heterogeneities, we estimate that simultaneous brittle failure of heterogeneities could produce tremor bursts with equivalent seismic moments of 4.5 × 109–4.7 × 1014 N m, consistent with seismic moments estimated from geophysical data at active subduction zones.

We present observations from an exhumed subduction complex that resembles the environment of modern deep episodic tremor and slow slip (ETS). We focus on the Cycladic Blueschist Unit on Syros Island in Greece. Syros metabasites consist of blueschists and eclogites that record prograde deformation, with peak metamorphism of 1200-1600 MPa and 450-550 °C. Field observations reveal that coexistence of blueschist and eclogite sets up an important rheological contrast: blueschists show distributed viscous dislocation creep, whereas eclogites dispersed within the blueschist matrix show brittle shear fractures and veins. These observations are consistent with the inferred prominent role of high fluid pressures from geophysical studies, but are inconsistent with models of deep ETS that invoke changes in rate-and-state friction parameters along a narrow fault. Instead, we suggest deep ETS may be controlled by coupled brittle-viscous deformation in partially eclogitized oceanic crust embedded within high-fluid-pressure patches along the plate interface.

Exhumed high‐pressure/low‐temperature (HP/LT) metamorphic rocks provide insights into deep (∼20–70 km) subduction interface dynamics. On Syros Island (Cyclades, Greece), the Cycladic Blueschist Unit preserves blueschist‐to‐eclogite facies oceanic‐ and continental‐affinity rocks that record the structural and thermal evolution linked to Eocene subduction. Despite decades of research, the metamorphic and deformation history (P‐T‐D) and timing of subduction and exhumation are matters of ongoing discussion. We suggest that Syros comprises three coherent tectonic slices and that each slice underwent subduction, underplating, and syn‐subduction return flow along similar P‐T trajectories, but at progressively younger times. Subduction and exhumation are distinguished by lineations and ductile fold axis orientations, and are kinematically consistent with previous studies that document top‐to‐the‐S‐SW shear (prograde‐to‐peak subduction), top‐to‐the‐NE shear (blueschist facies exhumation), and then E‐W coaxial stretching (greenschist facies exhumation). Amphibole zonations record cooling during decompression, indicating return flow above a cold slab. Multi‐mineral Rb‐Sr isochrons and compiled metamorphic geochronology show that the three slices record distinct stages of peak subduction (53–52, ∼50, and 45 Ma) that young with structural depth. Retrograde blueschist and greenschist facies fabrics span ∼50–40 and ∼43–20 Ma, respectively, and also young with structural depth. Synthesized data sets support a revised tectonic framework for Syros, involving subduction of structurally distinct coherent slices and simultaneous return flow of previously accreted tectonic slices in the subduction channel shear zone. Distributed, ductile, dominantly coaxial return flow in an Eocene‐Oligocene subduction channel proceeded at rates of ∼1.5–5 mm/yr and accommodated ∼80% of the total exhumation of this HP/LT complex.