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: Analytical and numerical models of viscous anisotropy: a toolset to constrain the role of mechanical anisotropy for regional tectonics and fault loading
o what extent mechanical anisotropy is required to explain the dynamics of the lithosphere is an important yet unresolved question. If anisotropy affects stress and deformation, and hence processes such as fault loading, how can we quantify its role from observations? Here, we derive analytical solutions and build a theoretical framework to explore how a shear zone with linear anisotropic viscosity can lead to deviatoric stress heterogeneity, strain-rate enhancement, as well as non-coaxial principal stress and strain rate. We develop an open-source finite-element software based on FEniCS for more complicated scenarios in both 2-D and 3-D. Mechanics of shear zones with transversely isotropic and orthorhombic anisotropy subjected to misoriented shortening and simple shearing are explored. A simple regional example for potential non-coaxiality for the Leech River Schist above the Cascadia subduction zone is presented. Our findings and these tools may help to better understand, detect and evaluate mechanical anisotropy in natural settings, with potential implications including the transfer of lithospheric stress and deformation through fault loading.  more » « less
Award ID(s):
2045292 1927216 2121666
PAR ID:
10549437
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
239
Issue:
2
ISSN:
0956-540X
Page Range / eLocation ID:
950 to 963
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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,PnandPtomography, surface wave tomography, and splitting ofSKSand 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.Pnfast 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
  2. Abstract Large continental transform faults are thought to cross‐cut the crust and extend into the lithospheric mantle. However, the strain distribution associated with these faults at mantle depths is not well understood. At the heart of this question is the rheology of the lithosphere: when stressed by displacement on a trans‐crustal fault, does the uppermost mantle localize strain in a continuing narrow fault zone, or is that strain instead distributed in a shear zone that widens with increasing depth? This study uses harmonic decomposition of receiver functions to measure the spatial and depth distribution of seismic anisotropy, a proxy for viscous deformation, around the Alpine Fault in Aotearoa New Zealand. Anisotropy aligned with the fault is present in the lithosphere at least 100 km away from the fault trace, suggesting that the Alpine Fault shear zone widens at depth. 
    more » « less
  3. We perform numerical experiments of damped quasi-dynamic fault slip that include a rate-and-state behavior at steady state to simulate earthquakes and a plastic rheology to model permanent strain. The model shear zone has a finite width which represents a natural fault zone. Here we reproduce fast and slow events that follow theoretical and observational scaling relationships for earthquakes and slow slip events (SSEs). We show that the transition between fast and slow slip occurs when the friction drop in the shear zone is equal to a critical value, Δμc. With lower friction drops, SSEs use nearly all of mechanical work to accumulate inelastic strain, while with higher friction drops fast slips use some of the mechanical work to slip frictionally. Our new formulation replaces the state evolution of rate and state by the stress evolution concurrent with accumulation of permanent damage in and around a fault zone. 
    more » « less
  4. Abstract Mantle xenoliths from the Southern Alps, New Zealand, provide insight into the origin of mantle seismic anisotropy related to the Australian‐Pacific plate boundary. Most xenoliths from within 100 km lateral distance of the Alpine Fault are coarse grained, but a small number are finer grained protomylonites. The protomylonites contain connected networks of fine grains with a different crystallographic preferred orientation (CPO) to coarse porphyroclasts in the same xenolith, suggesting that protomylonites and coarse‐grained samples record different deformation kinematics. The CPOs of fine grains in protomylonites have monoclinic symmetry, with the 2‐fold rotation axis normal to a plane that contains olivine [010] and orthopyroxene [100] maxima, suggesting that the protomylonite deformation involved significant simple shear. Some coarse‐grained samples contain unconnected lenses and layers of fine grains with the same CPO as the coarse grains. Microstructures suggest that these fine grains formed by subgrain rotation recrystallization and that protomylonites may represent an up‐strain progression of this microstructure, where the connectivity of fine grains has allowed them to localize shear and develop a new Alpine Fault CPO. The samples tell us about the state of the mantle at 25 Ma, in the early history of the plate boundary. If this suite of samples is representative of the mantle beneath the Alpine Fault in the present day, then we can interpret the complex seismic anisotropy patterns in the lithospheric mantle as representative of blocks containing variably rotated older CPOs juxtaposed by narrow shear zones associated with Alpine Fault deformation. 
    more » « less
  5. Shear wave splitting of teleseismic core phases such as SKS is commonly used to constrain mantle seismic anisotropy, a proxy for convective deformation. In plate boundaries, sharp lateral variations of splitting measurements near transform faults are often linked to deformation within a lithospheric shear zone below, but potential seismic waveform effects from heterogeneous structure on small scales may influence the interpretation. Here, we explore possible finite frequency effects on shear wave splitting near fault zones in a fully three‐dimensional anisotropic setting. We find that shear zones wider than 80 km, a scale set by the Fresnel zone, can be clearly detected, but narrower zones are less distinguishable. Near the edge of the shear zone, the combined effect of anisotropy and scattering generates false splitting measurements with large delay times and fast axis orientation approaching the back‐azimuth, a bias which can only be identified when records from different back‐azimuths are analyzed together. This substantiates that back‐azimuthal variations of splitting can arise not just from vertical layering but also lateral changes of anisotropic media. We also test the effects of shear zone edge geometry, epicentral distance, filtering frequency, crustal thickness, and sediment cover. Our study delineates the ability of shear wave splitting to resolve and investigate fault zones, and emphasizes the importance of good azimuthal coverage to correctly interpret observed anisotropy. Based on revisiting previous shear wave splitting and lithospheric deformation studies, we infer that many crustal fault zones are underlain by lithospheric shear zones at least 20 km wide. 
    more » « less