Feedbacks between metasomatism, rheological heterogeneities and strain localization in deep subduction interface shear zones

Characterizing deep subduction dynamics is crucial for understanding processes of high-pressure-rock exhumation, fluid flow, seismicity and transient slip events. Metamorphic phase transformations at the blueschist-to-eclogite transition induce important rheological changes, commonly transitioning from a more brittle to a mixed brittle-viscous plate interface rheology. This shift may promote slip transients such as slow slip and tremor (SST) observed in modern subduction zones. Geophysical and geologic data as well as numerical models suggest that slow slip is likely accommodated along weak, fluid-rich shear zones and that accompanying tremor may represent km-scale brittle asperities embedded within localized slip zones. Here we use the geologic record exposed on Syros Island (Greece) to investigate the relationships between strain localization and fluid-rock interactions along the deep megathrust, and explore their implications for SST.

We used high-resolution drone surveys, along with microstructural, geochemical, and petrologic data, to examine a blueschist-to-eclogite facies subduction shear zone in the Kampos Belt near Grizzas locality on northern Syros. The estimated P-T conditions are comparable to the SST zone along active warm subduction margins such as Cascadia and Central Chile. Our approach involved mapping strain and lithologies, constructing a 3D geological model, and performing detailed analyses of localized shear zones and metasomatic rocks.

At the hectometre-scale, the Grizzas locality exposes a stack of progressively underplated oceanic and metasedimentary rocks. Individual slices include brittlely deformed metagabbros up to 200 m-thick, weakly-strained to undeformed igneous breccias up to 30 m-thick, and foliated quartz-mica schists. These slices are repeated along five localized shear zones composed of chlorite-tremolite and glaucophane schists that are less than 10 m-thick. Fine-scale characterization of one of these shear zones reveal several discrete intercalations of blueschists/glaucophanites, tremolite-chlorite schists and metasediments. Microstructural and petrologic analyses suggest that blueschist/glaucophanite layers formed through the transformation of a gabbro/blueschist breccia precursor, likely induced by along-dip fluid influx. This metasomatic process extensively replaced the precursor gabbro fabric with nearly pure glaucophane and also enhanced the development of high-strain zones. Geochemical analyses indicate the formation of tremolite-chlorite (+/- talc) schists through chemical exchange between metamafic and metaultramafic rocks or by the interaction with serpentinite-derived fluids. This is supported by the presence of partially-digested metagabbro pods which contain garnet and chlorite with anomalously high Cr₂O₃ contents (up to 1.2 and 2.1 Wt%, respectively) as well as omphacitites associated with glaucophane-phengite veins and glaucophane-bearing veins crosscutting the chlorite schists.

We suggest that metasomatism triggered localized deformation around gabbro blocks and permitted repeated down-slicing and underplating of subducting oceanic material on the deep subduction interface. The metasomatism likely exploited precursory features such as lithological contacts, fractures, and/or fabric heterogeneity, to transiently increase permeability and allow further fluid ingress eventually resulting in the development of major shear zones. The degree of localization in these major shear zones and the concentration of foliated phyllosilicates within them mean they may have been capable of hosting slow slip (to be explored further), and the up-to-km-scale of brittlely deformed metagabbro blocks embedded between the shear zones are compatible with tremor sources.