skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Serpentinization, Deformation, and Seismic Anisotropy in the Subduction Mantle Wedge
Abstract Antigorite is a hydrous sheet silicate with strongly anisotropic seismic and rheological properties. Hydrous minerals such as antigorite have been invoked to explain numerous geologic observations within subduction zones including intermediate‐depth earthquakes, arc volcanism, the persistent weakness of the subduction interface, trench‐parallelSwave splitting, and episodic tremor and slip. To understand how the presence of antigorite‐bearing rocks affects observations of seismic anisotropy, three mylonites from the Kohistan palaeo‐island arc in Pakistan were analysed using electron backscatter diffraction. A fourth sample, which displayed optical evidence for crystallographically controlled replacements of olivine, was also investigated using electron backscatter diffraction to identify potential topotactic relationships. The resulting data were used to model the bulk seismic properties of antigorite‐rich rocks. The mylonitic samples exhibit incredibly strong bulk anisotropy (10–20% for the antigorite + olivine). Within the nominally undeformed protomylonite, two topotactic relationships were observed: (1) (010)ant//(100)ol with [100]ant//[001]ol and (2) (010)ant//(100)ol with [100]ant//[010]ol. However, the strength of a texture formed by topotactic replacement is markedly weaker than the strength of the textures observed in mylonitic samples. Since antigorite is thought to be rheologically weak, we hypothesise that microstructures formed from topotactic reactions will be progressively overprinted as deformation is localised in regions with greater percentages of serpentine. Regions of highly sheared serpentine, therefore, have the potential to strongly influence seismic wave speeds in subduction settings. The presence of deformed antigorite in a dipping structure is one explanation for observations of both the magnitude and splitting pattern of seismic waves in subduction zones.  more » « less
Award ID(s):
1848824
PAR ID:
10457043
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geochemistry, Geophysics, Geosystems
Volume:
21
Issue:
4
ISSN:
1525-2027
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Iron‐rich olivine is mechanically weaker than olivine of mantle composition, ca. Fo90, and thus is more amenable to study under a wide range of laboratory conditions. To investigate the effects of iron content on deformation‐produced crystallographic preferred orientation (CPO) and grain size, we analyzed the microstructures of olivine samples with compositions of Fo70, Fo50, and Fo0that were deformed in torsion under either anhydrous or hydrous conditions at 300 MPa. Electron backscatter diffraction (EBSD) observations reveal a transition in CPO from D‐type fabric, induced by dislocation glide on both the (010)[100] and the (001)[100] slip systems, at low strains, to A‐type fabric, caused by dislocation glide on the (010)[100] slip system, at high strains for all of our samples, independent of iron content and hydrous/anhydrous conditions. A similar evolution of fabric with increasing strain is also reported to occur for Fo90. Radial seismic anisotropy increases with increasing strain, reaching a maximum value of ∼1.15 at a shear strain of ∼3.5 for each sample, demonstrating that the seismic anisotropy of naturally deformed olivine‐rich rocks can be well approximated by that of iron‐rich olivine. Based on EBSD observations, we derived a piezometer for which recrystallized grain size decreases inversely with stress to the ∼1.2 power. Also, recrystallized grain size increases with increasing iron content. Our experimental results contribute to understanding the microstructural evolution in the mantle of not only Earth but also Mars, where the iron content in olivine is higher. 
    more » « less
  2. Subduction zones are often characterized by the presence of strong trench-parallel seismic anisotropy and large delay times. Hydrous minerals, owing to their large elastic anisotropy and strong lattice preferred orientations (LPOs), are often invoked to explain these observations. However, the elasticity and the LPO of chloritoid, which is one of such hydrous phases relevant in subduction zone settings, are poorly understood. In this study, we measured the LPO of polycrystalline chloritoid in natural rock samples, obtained the LPO-induced seismic anisotropy, and evaluated the thermodynamic stability field of chloritoid in subduction zones. The LPO of chloritoid aggregates displayed a strong alignment of the [001] axes subnormal to the rock foliation, with a girdle distribution of the [100] axes and the (010) poles subparallel to the foliation. New elasticity data of single-crystal chloritoid showed a strong elastic anisotropy of chloritoid with 47% for S-waves (V S ) and 22% for P-waves (V P ), respectively. The combination of the LPO and the elastic anisotropy of the chloritoid aggregates produced a strong S-wave anisotropy with a maximum AV S of 18% and a P-wave anisotropy with an AV P of 10%. The role of chloritoid LPO in seismic anisotropy was evaluated in natural rock samples and a hypothetical blueschist. Our results indicate that the strong LPO of chloritoid along the subduction interface and in subducting slabs can influence the trench-parallel seismic anisotropy in subduction zones with “cold” geotherms. 
    more » « less
  3. Abstract Shear‐wave splitting observations can provide insight into mantle flow, due to the link between the deformation of mantle rocks and their direction‐dependent seismic wave velocities. We identify anisotropy in the Cook Inlet segment of the Alaska subduction zone by analyzing splitting parameters of S waves from local intraslab earthquakes between 50 and 200 km depths, recorded from 2015–2017 and emphasizing stations from the Southern Alaska Lithosphere and Mantle Observation Network experiment. We classify 678 high‐quality local shear‐wave splitting observations into four regions, from northwest to southeast: (L1b) splitting measurements parallel to Pacific plate motion, (L1a) arc‐perpendicular splitting pattern, (L2) sharp transition to arc‐parallel splitting, and (L3) splitting parallel to Pacific plate motion. Forward modeling of splitting from various mantle fabrics shows that no one simple model fully explains the observed splitting patterns. An A‐type olivine fabric with fast direction dipping 45° to the northwest (300°)—aligned with the dipping slab—predicts fast directions that fit L1a observations well, but not L2. The inability of the forward model fabrics to fit all the observed splitting patterns suggests that the anisotropy variations are not due to variable ray angles, but require distinct differences in the anisotropy regime below the arc, forearc, and subducting plate. 
    more » « less
  4. Abstract Using numerical models, we compute the evolution of the mantle flow field and the crystal preferred orientation (CPO) of mineral aggregates in the mantle wedge of generic subduction systems from their nascent to mature stage and investigate shear wave splitting (SWS) through the forearc mantle wedge corner and overriding crust. Upon subduction initiation, the maximum depth of slab‐mantle decoupling (MDD) is relatively shallow (∼20 km depth), resulting in mantle flow and CPO development in the wedge corner. As subduction continues, the MDD deepens, the wedge corner cools and stagnates, and the olivine CPO becomes frozen‐in. In the cool wedge corner, antigorite can form if water is available. In non‐deforming mantle, antigorite CPO develops relative to the host olivine CPO through topotactic growth. We calculate splitting parameters of synthetic local S waves based on the model‐predicted A‐ and B‐type olivine CPOs and topotactically grown antigorite CPO that replaces A‐type olivine CPO in the wedge corner. The fast direction is trench‐normal for A‐type olivine and antigorite CPOs and trench‐parallel for B‐type. When the delay times are long enough (>0.1 s), we find them positively correlated with the thickness of the mantle wedge corner. In NE Japan, where the results of detailed analyses on the spatial variation of the SWS parameters are available, such correlation is not observationally reported. However, the addition of an anisotropic overriding crust provides delay times (∼0.1 s) and trench‐normal fast directions that are consistent with the local SWS observations. 
    more » « less
  5. Abstract Seismic anisotropy constitutes a useful tool for imaging the structure along the plate interface in subduction zones, but the seismic properties of mafic blueschists, a common rock type in subduction zones, remain poorly constrained. We applied the technique of electron backscatter diffraction (EBSD) based petrofabric analysis to calculate the seismic anisotropies of 14 naturally deformed mafic blueschists at dry, ambient conditions. The ductilely deformed blueschists were collected from terranes with inferred peak P‐T conditions applicable to subducting slabs at or near the plate interface in active subduction zones. Epidote blueschists display the greatestPwave anisotropy range (AVp ∼7%–20%), while lawsonite blueschist AVp ranges from ∼2% to 10%.Swave anisotropies generate shear wave splitting delay times up to ∼0.1 s over a thickness of 5 km. AVp magnitude increases with glaucophane abundance (from areal EBSD measurements), decreases with increasing epidote or lawsonite abundance, and is enhanced by glaucophane crystallographic preferred orientation (CPO) strength. Two‐phase rock recipe models provide further evidence of the primary role of glaucophane, epidote, and lawsonite in generating blueschist seismic anisotropy. The symmetry ofPwave velocity patterns reflects the deformation‐induced CPO type in glaucophane—an effect previously observed for hornblende on amphibolitePwave anisotropy. The distinctive seismic properties that distinguish blueschist from other subduction zone rock types and the strong correlation between anisotropy magnitude/symmetry and glaucophane CPO suggest that seismic anisotropy may be a useful tool in mapping the extent and deformation of blueschists along the interface, and the blueschist‐eclogite transition in active subduction zones. 
    more » « less