The Northeastern Tibetan Plateau has thickened crust and is still undergoing strong active crustal shortening and deformation. Crustal anisotropy can provide clues to how the crust is currently deforming and evolving. We use an automatic method to analyse the upper-crustal anisotropy of the NE Tibetan Plateau and the adjacent region using local earthquakes recorded at 39 permanent seismic stations during the period 2009–2018. The majority of the dominant fast directions are consistent with the maximum horizontal stress orientation, suggesting that the upper-crustal anisotropy is mainly controlled by the regional or local stress field. Several fault-parallel measurements are observed for stations on or near to the main faults. These fault-parallel fast directions indicate that the main mechanism of upper-crustal anisotropy is associated with shear fabric caused by deformation. Fast directions neither fault-parallel nor stress-parallel are observed at stations lying several kilometres away from fault zones, likely reflecting the combined influence of stress-aligned microcracks and active faults. A comparison between our upper-crustal anisotropy parameters and those inferred from previous anisotropy studies that used receiver function and teleseismic shear wave splitting measurements suggests that the crust has the same deformation mechanisms as mantle anisotropy in the southern part of the Western Qinling Fault, whereas the upper-crustal anisotropic mechanism is different from those of lower crust and mantle anisotropy in the northern part of the Western Qinling Fault. These observations imply that the Western Qinling Fault may be an important boundary fault.
To discern spatial and explore possible existence of temporal variations of upper crustal anisotropy in an ∼15 km section of the San Jacinto Fault Zone (SJFZ) that is composed of the Buck Ridge and Clark faults in southern California, we conduct a systematic shear wave splitting investigation using local S‐wave data recorded by three broadband seismic stations located near the surface expression of the SJFZ. An automatic data selection and splitting measurement procedure is first applied, and the resulting splitting measurements are then manually screened to ensure reliability of the results. Strong spatial variations in crustal anisotropy are revealed by 1,694 pairs of splitting parameters (fast polarization orientation and splitting delay time), as reflected by the dependence of the resulting splitting parameters on the location and geometry of the raypaths. For raypaths traveling through the fault zones, the fast orientations are dominantly WNW‐ESE which is parallel to the faults and may be attributed to fluid‐filled fractures in the fault zones. For non‐fault‐zone crossing raypaths, the fast orientations are dominantly N–S which are consistent with the orientation of the regional maximum compressive stress. A three‐dimensional model of upper crustal anisotropy is constructed based on the observations. An increase in the raypath length normalized splitting times is observed after the 03/11/2013 M4.7 earthquake, which is probably attributable to changes in the spatial distribution of earthquakes before and after the M4.7 earthquake rather than reflecting temporal changes of upper crustal anisotropy.
more » « less- NSF-PAR ID:
- 10359821
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 4
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
SUMMARY -
more » « less
Abstract Seismic azimuthal anisotropy characterized by shear wave splitting analyses using teleseismic
X K S phases (includingS K S ,S K K S , andP K S ) is widely employed to constrain the deformation field in the Earth's crust and mantle. Due to the near‐vertical incidence of theX K S arrivals, the resulting splitting parameters (fast polarization orientations and splitting times) have an excellent horizontal but poor vertical resolution, resulting in considerable ambiguities in the geodynamic interpretation of the measurements. Here we useP ‐to‐S converted phases from the Moho and the 410‐ (d 410) and 660‐km (d 660) discontinuities to investigate anisotropy layering beneath Southern California. Similarities between the resulting splitting parameters from theX K S andP ‐to‐S converted phases from thed 660 suggest that the lower mantle beneath the study area is azimuthally isotropic. Similarly, significant azimuthal anisotropy is not present in the mantle transition zone on the basis of the consistency between the splitting parameters obtained usingP ‐to‐S converted phases from thed 410 andd 660. Crustal anisotropy measurements exhibit a mean splitting time of 0.2 ± 0.1 s and mostly NW‐SE fast orientations, which are significantly different from the dominantly E‐W fast orientations revealed usingX K S andP ‐to‐S conversions from thed 410 andd 660. Anisotropy measurements using shear waves with different depths of origin suggest that the Earth's upper mantle is the major anisotropic layer beneath Southern California. Additionally, this study demonstrates the effectiveness of applying a set of azimuthal anisotropy analysis techniques to reduce ambiguities in the depth of the source of the observed anisotropy. -
more » « less
Abstract Mature faults with large cumulative slip often separate rocks with dissimilar elastic properties and show asymmetric damage distribution. Elastic contrast across such bimaterial faults can significantly modify various aspects of earthquake rupture dynamics, including normal stress variations, rupture propagation direction, distribution of ground motions, and evolution of off‐fault damage. Thus, analyzing elastic contrasts of bimaterial faults is important for understanding earthquake physics and related hazard potential. The effect of elastic contrast between isotropic materials on rupture dynamics is relatively well studied. However, most fault rocks are elastically anisotropic, and little is known about how the anisotropy affects rupture dynamics. We examine microstructures of the Sandhill Corner shear zone, which separates quartzofeldspathic rock and micaceous schist with wider and narrower damage zones, respectively. This shear zone is part of the Norumbega fault system, a Paleozoic, large‐displacement, seismogenic, strike‐slip fault system exhumed from middle crustal depths. We calculate elastic properties and seismic wave speeds of elastically anisotropic rocks from each unit having different proportions of mica grains aligned sub‐parallel to the fault. Our findings show that the horizontally polarized shear wave propagating parallel to the bimaterial fault (with fault‐normal particle motion) is the slowest owing to the fault‐normal compliance and therefore may be important in determining the elastic contrast that affects rupture dynamics in anisotropic media. Following results from subshear rupture propagation models in isotropic media, our results are consistent with ruptures preferentially propagated in the slip direction of the schist, which has the slower horizontal shear wave and larger fault‐normal compliance.
-
more » « less
Abstract To systematically investigate seismic azimuthal anisotropy in the Sumatra subduction zone and probe mantle dynamics associated with the subduction of the Australian Plate beneath the Sunda Plate, a total of 169 pairs of teleseismic XKS (including PKS, SKKS, SKS) and 115 pairs of local
S splitting parameters are obtained using broadband seismic data recorded at ~70 stations. Additionally, crustal anisotropy in the overriding Sunda Plate is measured by analyzing the moveout ofP ‐to‐S conversions from the Moho using a sinusoidal function. Comparison between the three sets of anisotropy measurements obtained using shear waves with different depths of origin suggests that (1) the crust of the Sunda Plate is anisotropic with mostly trench‐parallel fast orientations and a mean splitting time of 0.28 ± 0.05 s; (2) the mantle wedge is azimuthally anisotropic with dominantly trench‐parallel fast orientations and splitting times ranging from 0.22 to 0.81 s, which generally increase with the focal depth; and (3) subslab anisotropy is mostly trench‐normal beneath the fore‐arc region with an averaged splitting time of 1.48 ± 0.06 s, and becomes trench‐parallel beneath the arc and back‐arc areas with a mean splitting time of 0.33 ± 0.04 s. The resulting lateral and vertical distributions of anisotropy obtained using splitting of three types of shear waves advocate the presence of an entrained subslab flow that is deflected by the mantle transition zone. The flow enters the mantle wedge through a slab window and flows horizontally parallel to the trench. -
more » « less
Abstract Seismic anisotropy beneath eastern North America, as expressed in shear wave splitting observations, has been attributed to plate motion‐parallel shear in the asthenosphere, resulting in fast axes aligned with the plate motion. However, deviations of fast axes from plate motion directions are observed near major tectonic boundaries of the Appalachians, indicating contributions from lithospheric anisotropy associated with past tectonic processes. In this study, we conduct anisotropic receiver function (RF) analysis using data from a dense seismic array traversing the New England Appalachians in Connecticut to examine anisotropic layers in the crust and upper mantle and correlate them with past tectonic processes as well as present‐day mantle flow. We use the harmonic decomposition method to separate directionally‐dependent variations of RFs and focus on features with the same harmonic signals observed across multiple stations. Within the crust, there are multiple features that may be correlated with stratification in the Hartford Basin, faults in the Taconic thrust belt, shear zones formed during Salinic/Acadian terrane accretion events, and orogen‐parallel crustal flow in the Acadian orogenic plateau. We apply a Bayesian inversion method to obtain quantitative constraints on the direction and strength of intra‐crustal anisotropy beneath the Hartford Basin. In the upper mantle, we identify a fossil shear zone possibly formed during oblique subduction of Rheic Ocean lithosphere. We also find evidence for a plate motion‐parallel flow zone in the asthenosphere that is likely disturbed by mantle upwelling near the southern margin of the Northern Appalachian Anomaly in the eastern part of the study area.
