skip to main content


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
NSF-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

    Phase A (Mg7Si2O8(OH)6) is one of the important dense hydrous magnesium silicates in subducting slab, because it forms after the breakdown of antigorite serpentine and could be the dominant hydrous phase in the upper‐mantle deep slab for the water transportation into deep Earth. In this study, the compressional (P) and shear (S) wave velocities of phase A were measured at simultaneous pressure and temperature conditions up to 11 GPa and 1073 K. Combined with elastic properties of olivine and pyroxenes, we calculate the hydration effect on the velocities of harzburgite lithology in cold subduction zones throughout the depth ranges where phase A is thermodynamically stable. Our calculations suggest that the hydration increases bothP‐ andS‐wave velocities of harzburgite; at ∼5 wt% hydration, its seismic detectability is enhanced by 1%–1.5% in velocity contrasts relative to its anhydrous counterpart.

     
    more » « less
  4. 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
  5. Abstract

    Shear‐wave splitting observations of SKS and SKKS phases have been used widely to map azimuthal anisotropy, as caused by the occurrence of olivine, to constrain the dominant directions of upper mantle deformation. While SK(K)S splitting measurements at individual seismic stations are often averaged before interpretation, it is useful to consider additional information, for example, based on the variation of splitting parameters with azimuth due to the non‐vertical incidence of core‐phases. These constraints in theory enable a differentiation between various types of olivine and may allow us to infer otherwise poorly known upper mantle parameters such as stress, temperature, and water content. In this study, we predict the azimuthal variation of splitting parameters for A‐, C‐, and E‐type olivine fabrics and match them with observations from the High Lava Plains, Northwestern Basin and Range, and Western Yellowstone Snake River Plain in the Pacific Northwest US. This helps to constrain the amount of water in the upper mantle in the back‐arc of the Cascadia subduction zone, known for its consistent E‐W oriented seismic anisotropy, and particularly large splitting delay times. Our investigation renders the C‐type olivine mechanism improbable for this location; A‐ and E‐type fabrics match the observations, although differentiating between them is difficult. However, the agreement of the amplitude of backazimuthal variation of the fast orientation, plus the potential to explain large splitting delay times, suggest the occurrence of E‐type olivine, and thus the likely presence of a moderately hydrated upper mantle beneath Cascadia's back‐arc.

     
    more » « less