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‐parallel
- Award ID(s):
- 1753125
- NSF-PAR ID:
- 10316125
- Date Published:
- Journal Name:
- Frontiers in Earth Science
- Volume:
- 9
- ISSN:
- 2296-6463
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract S wave 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
Abstract Seismic anisotropy produced by aligned olivine in oceanic lithosphere offers a window into mid‐ocean ridge (MOR) dynamics. Yet, interpreting anisotropy in the context of grain‐scale deformation processes and strain observed in laboratory experiments and natural olivine samples has proven challenging due to incomplete seismological constraints and length scale differences spanning orders of magnitude. To bridge this observational gap, we estimate an in situ elastic tensor for oceanic lithosphere using co‐located compressional‐ and shear‐wavespeed anisotropy observations at the NoMelt experiment located on ∼70 Ma seafloor. The elastic model for the upper 7 km of the mantle, NoMelt_SPani7, is characterized by a fast azimuth parallel to the fossil‐spreading direction, consistent with corner‐flow deformation fabric. We compare this model with a database of 123 petrofabrics from the literature to infer olivine crystallographic orientations and shear strain accumulated within the lithosphere. Direct comparison to olivine deformation experiments indicates strain accumulation of 250%–400% in the shallow mantle. We find evidence for D‐type olivine lattice‐preferred orientation (LPO) with fast [100] parallel to the shear direction and girdled [010] and [001] crystallographic axes perpendicular to shear. D‐type LPO implies similar amounts of slip on the (010)[100] and (001)[100] easy slip systems during MOR spreading; we hypothesize that grain‐boundary sliding during dislocation creep relaxes strain compatibility, allowing D‐type LPO to develop in the shallow lithosphere. Deformation dominated by dislocation‐accommodated grain‐boundary sliding (disGBS) has implications for in situ stress and grain size during MOR spreading and implies grain‐size dependent deformation, in contrast to pure dislocation creep.
-
more » « less
Abstract The style of convective force transmission to plates and strain‐localization within and underneath plate boundaries remain debated. To address some of the related issues, we analyze a range of deformation indicators in southern California from the surface to the asthenosphere. Present‐day surface strain rates can be inferred from geodesy. At seismogenic crustal depths, stress can be inferred from focal mechanisms and splitting of shear waves from local earthquakes via crack‐dependent seismic velocities. At greater depths, constraints on rock fabrics are obtained from receiver function anisotropy,
P n andP tomography, surface wave tomography, and splitting ofSKS and other teleseismic core phases. We construct a synthesis of deformation‐related observations focusing on quantitative comparisons of deformation style. We find consistency with roughly N‐S compression and E‐W extension near the surface and in the asthenospheric mantle. However, all lithospheric anisotropy indicators show deviations from this pattern.P n fast axes and dipping foliations from receiver functions are fault‐parallel with no localization to fault traces and match post‐Farallon block rotations in the Western Transverse Ranges. Local shear wave splitting orientations deviate from the stress orientations inferred from focal mechanisms in significant portions of the area. We interpret these observations as an indication that lithospheric fabric, developed during Farallon subduction and subsequent extension, has not been completely reset by present‐day transform motion and may influence the current deformation behavior. This provides a new perspective on the timescales of deformation memory and lithosphere‐asthenosphere interactions. -
more » « less
Abstract In approximately one fourth of worldwide subduction zones, seismic observations indicate a rotation from trench‐normal to trench‐parallel fast axis orientations in the mantle wedge. To interpret this observation we predict the evolution of crystal lattice preferred orientation in mantle wedge material as a function of the amount of water by using a model of polycrystal deformation (D‐Rex) coupled with an analytical corner flow. The resulting seismic signature is obtained from synthetic seismic wave propagation through this mantle wedge. We identify that progressive hydration produces the rotation of fast axis orientations and can generate between the two zones of trench‐parallel and trench‐normal fast axis orientations a morph zone with very small anisotropy and a related decrease in
P andS wave velocities. Such a morph zone is not produced by trench‐parallel flow, hence this signature can be used to detect water in the mantle wedge. -
more » « less
SUMMARY Determinations of seismic anisotropy, or the dependence of seismic wave velocities on the polarization or propagation direction of the wave, can allow for inferences on the style of deformation and the patterns of flow in the Earth’s interior. While it is relatively straightforward to resolve seismic anisotropy in the uppermost mantle directly beneath a seismic station, measurements of deep mantle anisotropy are more challenging. This is due in large part to the fact that measurements of anisotropy in the deep mantle are typically blurred by the potential influence of upper mantle and/or crustal anisotropy beneath a seismic station. Several shear wave splitting techniques are commonly used that attempt resolve seismic anisotropy in deep mantle by considering the presence of multiple anisotropic layers along a raypath. Examples include source-side S-wave splitting, which is used to characterize anisotropy in the deep upper mantle and mantle transition zone beneath subduction zones, and differential S-ScS and differential SKS-SKKS splitting, which are used to study anisotropy in the D″ layer at the base of the mantle. Each of these methods has a series of assumptions built into them that allow for the consideration of multiple regions of anisotropy. In this work, we systematically assess the accuracy of these assumptions. To do this, we conduct global wavefield modelling using the spectral element solver AxiSEM3D. We compute synthetic seismograms for earth models that include seismic anisotropy at the periods relevant for shear wave splitting measurements (down to 5 s). We apply shear wave splitting algorithms to our synthetic seismograms and analyse whether the assumptions that underpin common measurement techniques are adequate, and whether these techniques can correctly resolve the anisotropy incorporated in our models. Our simulations reveal some inaccuracies and limitations of reliability in various methods. Specifically, explicit corrections for upper mantle anisotropy, which are often used in source-side direct S splitting and S-ScS differential splitting, are typically reliable for the fast polarization direction ϕ but not always for the time lag δt, and their accuracy depends on the details of the upper mantle elastic tensor. We find that several of the assumptions that underpin the S-ScS differential splitting technique are inaccurate under certain conditions, and we suggest modifications to traditional S-ScS differential splitting approaches that lead to improved reliability. We investigate the reliability of differential SKS-SKKS splitting intensity measurements as an indicator for lowermost mantle anisotropy and find that the assumptions built into the splitting intensity formula can break down for strong splitting cases. We suggest some guidelines to ensure the accuracy of SKS-SKKS splitting intensity comparisons that are often used to infer lowermost mantle anisotropy. Finally, we suggest a new strategy to detect lowermost mantle anisotropy which does not rely on explicit upper mantle corrections and use this method to analyse the lowermost mantle beneath east Asia.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3389/feart.2021.644958
