A review of a quarter century of International Workshops on Seismic Anisotropy in the crust (0IWSA–12IWSA)

In 25 years, the presence of azimuthally varying seismic anisotropy throughout the Earth’s crust has progressed from general denial to universal acceptance, so that many international geophysical meetings now have sessions on seismic anisotropy. Over this period, the proceedings of the biennial series of International Workshops in Seismic Anisotropy (IWSAs) have captured many of the notable advances in the theory, calculation, observation and interpretation of particularly shear-wave splitting (seismic birefringence) in the Earth’s crust. Shear-wave splitting is the almost-infallible indicator of seismic anisotropy along the ray path. This paper reviews 13 IWSA meetings (0IWSA–12IWSA) as a catalogue of 25 years of progress in seismic anisotropy. The evidence now suggests that shear-wave splitting monitors the low-level pre-fracturing deformation of the stress-aligned fluid-saturated microcracks pervading almost all in situ rocks in the crust. Shear-wave splitting indicates that microcracks are so closely spaced they are critical systems with all the universality, calculability, predictability, “butterfly wing’s” sensitivity, and deterministic chaos that that implies. This leads to a New Geophysics, where low-level deformation can be monitored with shear-wave splitting, future behaviour calculated–predicted with the anisotropic poro-elastic model of rock evolution, and in some circumstances even potentially controlled by feedback. We anticipate the New Geophysics will greatly invigorate IWSA.     

Seismic Velocities and Anisotropy of the Lower Continental Crust: A Review

—Seismic anisotropy is often neglected in seismic studies of the earth’s crust. Since anisotropy is a common property of many typically deep crustal rocks, its potential contribution to solving questions of the deep crust is evaluated. The anisotropic seismic velocities obtained from laboratory measurements can be verified by computations based on the elastic constants and on numerical data pertaining to the texture of rock-forming minerals. For typical lower crustal rocks the influence of layering is significantly less important than the influence of rock texture. Surprisingly, most natural lower crustal rocks show a hexagonal type of anisotropy. Maximum anisotropy is observed for rocks with a high content of aligned mica. It seems possible to distinguish between layered intrusives and metasediments on the basis of in situ measurements of anisotropy, which can thus be used to validate different scenarios of crustal evolution.     

青藏高原东北缘地壳各向异性及其动力学意义

利用中国地震科学探测台阵项目二期(ChinaArrayⅡ)81个台站的地震数据,使用时间域反褶积方法提取接收函数,挑选满足要求的高质量Ps震相,通过改进的剪切波分裂计算方法获取了53个台站共130对高质量各向异性参数对.地壳各向异性分析指出,研究区东南部地区地壳各向异性方向为NWW向,与XKS各向异性方向、GPS速度场三者近平行关系,说明该地区在青藏高原向欧亚大陆增生的过程中是一个耦合的连贯变形过程;位于研究区中、北部地区的地壳各向异性方向表现为NEE-SWW或E-W方向,与GPS速度场方向一致,而与XKS结果的偏振方向大角度相交,说明该地区受到青藏高原下地壳塑性管道流的影响,可能存在壳幔解耦作用.     

Seismic crust structure beneath the Aegean region in southwest Turkey from radial anisotropic inversion of Rayleigh and Love surface waves

The Turkish plate is covered by hundreds of accelerometer and broadband seismic stations with less than 50 km inter-station distance providing high-quality earthquake recordings within the last decade. We utilize part of these stations to extract the fundamental mode Rayleigh and Love surface wave phase and group velocity data in the period range 5–20 s to determine the crust structure beneath the Aegean region in southwest Turkey. The observed surface wave signals are interpreted using both single-station and two-station techniques. A tomographic inversion technique is employed to obtain the two-dimensional group velocity maps from the single-station group velocities. One-dimensional velocity–depth profiles under each two-dimensional mesh point, which are jointly interpreted to acquire the three-dimensional image of the shear-wave velocities underneath the study area, are attained by utilizing the least-squares inversion technique, which is repeated for both Rayleigh and Love surface waves. The isotropic crust structure cannot jointly invert the observed Rayleigh and Love surface waves where the radial anisotropic crust better describes the observed surface wave data. The intrusive magmatic activity related to the northward subducting African plate under the Turkish plate results the crust structure deformations, which we think, causing the observed radial anisotropy throughout complex pattern of dykes and sills. The magma flow resulting in the mineral alignment within dykes and sills contributes to the observed anisotropy. Due to the existence of dykes, the radial anisotropy in the upper crust is generally negative, i.e., vertically polarized S-waves (Vsv) are faster than horizontally polarized S-waves (Vsh). Due to the existence of sills, the radial anisotropy in the middle-to-lower crust is generally positive, i.e., horizontally polarized S-waves (Vsh) are faster than vertically polarized S-waves (Vsv). Similar radial anisotropic results to those of the single-station analyses are obtained by the two-station analyses utilizing the cross-correlograms. The widespread volcanic and plutonic rocks in the region are consistent with the current seismic interpretations of the crustal deformations.     

蒙古中南部地区地壳各向异性及其动力学意义

利用蒙古中南部地区布设的69套宽频带数字地震仪2011年8月—2013年7月记录的远震事件,使用时间域反褶积方法提取接收函数,并挑选高质量Pms震相,通过改进的剪切波分裂方法对研究区地壳各向异性参数进行了研究,最终获取了1473对各向异性参数.经过统计分析,有48个台站可以归纳出两个方向的各向异性,11台站得到单个方向的各向异性,而剩余10个台站各向异性方向比较发散.结果显示,各向异性在蒙古中南部地壳中呈不均匀分布,有54个台站得到了NE-SW向各向异性,快波偏振方向平均值为N58°E±16°,与最大水平主应力σHmax方向和区域内主要断层走向一致,说明这部分地壳各向异性的主要成因存在于上地壳,可能与流体填充的微裂隙有关.而NW-SE向各向异性在53个台站被观测到,各向异性方向变化范围平均N132°E±16°,与研究区大部分SKS分裂快波方向具有较好的一致性,说明下地壳成岩矿物晶体定向排列是各向异性的主要成因.研究区地壳各向异性的分层特征总体上支持岩石圈受到NE-SW向挤压的动力学模型.     

青藏高原东北缘壳幔各向异性研究:基于SKS和Pms震相分析

利用甘肃和青海两省固定宽频带地震台记录的远震波形资料,挑选高质量SKS震相,联合使用最小切向能量方法和旋转互相关方法获得230对高信噪比分裂参数;同时对接收函数中Pms震相随方位角的变化进行拟合,得到了24个台站的地壳各向异性分裂参数.整个区域SKS分裂快波方向均值为123°,Pms分裂快波方向均值为132°,且大部分区域SKS、Pms快波方向与地表构造走向相一致,说明青藏高原东北缘以岩石圈垂直连贯变形为主,地壳上地幔相互耦合.SKS、Pms分裂时差均值分别为1.0s和0.6s,显示地壳各向异性对于SKS分裂时差有较大贡献.昆仑断裂附近Pms、SKS分裂快波方向与昆仑断裂走向基本一致,说明昆仑断裂可能是岩石圈尺度深大断裂;而阿尔金断裂东缘二者快波方向显著差异意味着阿尔金断裂在东缘可能仅为地壳尺度的断裂.     

Large-scale coherent anisotropy of upper mantle beneath the Italian peninsula comparing quasi-Love waves and SKS splitting

We document quantitatively observations of quasi-Love waves obtained at permanent (Italian National Seismic Network) and temporary seismic stations deployed in Italy between 2003 and 2006 (Retreat, CAT/SCAN projects). We analyzed large earthquakes with source parameters that favor quasi-Love wave generation within this time-span, including the Sumatra–Andaman earthquake of 12/26/04. The presence or the absence of the quasi-Love phase is compared to the smoothed anisotropic pattern defined by the numerous SKS splitting measurements obtained in peninsular Italy, and to the Italian upper mantle structure as defined by seismic tomography. The large-scale anisotropic features, responsible for shear-wave splitting and documented also by Pn and surface-wave anisotropy, generally display the correct geometry to explain the scattered quasi-Love waves. Quasi-Love observations do not demand a tilted-axis anisotropic geometry. We argue instead for anisotropy with laterally-variable horizontal symmetry axis in the upper mantle below the Italian peninsula.     

Complex deep seismic anisotropy below the Scandinavian Mountains

Several seismological projects focused on the deep structure of the Scandinavian Mountains, in Norway and neighbouring Sweden. We use these recordings to study seismic anisotropy by analysing the birefringence of SKS and SKKS phases. These phases, which should be polarised radially, are split into an additional transverse component if they propagate through an anisotropic medium. Our results are directions Φ of the apparent fast shear wave polarisation and delay times δt between the split phases. For station KONO in Southern Norway, we find frequency-dependent Φ and δt values, indicating a depth-dependent anisotropy. Additionally, Φ and δt values vary with epicentre backazimuths in Norway, indicating a complex anisotropic structure in the crust and upper mantle. Stacking of the SKS/SKKS waveforms improves the signal-to-noise ratio along one station line and allows us to better determine the splitting parameters. A unique and complete model of the complex anisotropy cannot be obtained due to the limited observed backazimuth range. Near-surface tectonic structures correlate with the splitting pattern and thus the crust is one anisotropic layer in the region. Partly preferred orientations in the rock fabric at the surface can be correlated with Φ. Below one or more anisotropic layers must exist to explain the backazimuth- and frequency-dependent observations, as well as the long δt values (>2 s) which cannot be explained with crustal anisotropy alone. The spatial distribution of the splitting results indicates that different tectonics units, e.g. the Sveconorwegian, the Central and Northern Svecofennian and the Caledonian nappes, are each characterised by specific anisotropic signatures.     

青藏高原东南缘中-下地壳地震波各向异性成因——来自云南六合深源包体组构学的约束

本文通过对出露于青藏高原东南缘云南六合地区的新生代深源岩石包体(斜长角闪岩、角闪石岩和石榴石透辉岩)的显微组构和地震波各向异性的研究来约束新生代青藏高原东南缘的地壳各向异性.通过角闪石地质压力计计算得知斜长角闪岩、角闪石岩和石榴石透辉岩包体来源于地壳28～36km,为中-下地壳岩石包体.EBSD测量结果显示包体中角闪石的CPO (晶格优选定向)为Type-IV型和(100)[001]滑移,单斜辉石的CPO为SL型和(100)[001]滑移,暗示中-下地壳为高温强变形的特征.通过CPO数据计算获得斜长角闪岩、角闪石岩和石榴石透辉岩包体全岩VP各向异性为1.9%～13.3%,最大分裂的剪切波各向异性(AVS)为1.17%～8.01%.结合前人的研究结果,该地区的地壳岩石能够解释利用Pms震相测量获得的分裂延迟时间,表明云南西北地区的壳内各向异性源于中-下地壳矿物的定向排列.云南西北地区的Pms快波方向近NW-SE向分布并与SKS的快波方向相近,暗示岩石圈变形是耦合的,受控于青藏高原向东南挤出的构造背景.     

中国地震科学台阵两期观测资料近场记录揭示的南北地震带地壳剪切波分裂特征

利用中国地震科学台阵第一期(2011-01-2014-06)及部分中国地震科学台阵第二期(2013-02-2015-12)的流动地震台阵记录到的小震波形资料,运用剪切波分裂系统分析(SAM)方法,分析南北地震带的地壳各向异性,对剪切波分裂参数所反映的区域应力环境及构造特征,以及区域内主压应力方向与断裂分布的关系展开讨论.研究结果表明,南北地震带快剪切波偏振方向自北向南由NE向逐渐转变为NNW向,与南北地震带区域主压应力的方向变化具有一致性.区域内分布的大量NE及WNW或NW向断裂构造同样对快波偏振方向有比较大的影响,位于走滑断裂附近的台站,其快波方向与断裂走向大致平行,部分位于走滑断裂附近的台站其快波方向几乎垂直于断裂走向,而与构造应力场方向一致性较好.个别台站表现出复杂快波优势方向特征,反映出研究区内构造环境的复杂性.慢波时间延迟结果显示,南北地震带南段的平均时间延迟高于北段,反映了受印度板块和欧亚板块的碰撞挤压作用,南段地壳介质各向异性程度更大,构造变形更加剧烈.对比南北地震带上地幔各向异性特征,推测在川滇菱形块体内部可能存在复杂的壳幔耦合现象,地壳剪切波分裂除了反映区域应力特征,还可以揭示出区域构造信息.     

Calculated and observed seismic behaviour of exposed upper mantle

In order to determine the seismic properties of rocks representing the continental crust/upper mantle transition, a structural and petrographical study of peridotites of the Ivrea-Verbano zone (N. Italy) has been undertaken. The seismic behaviour (Vp, Vs, birefringence, orientation of the fast shear-wave polarisation plane) has been calculated from the modal composition of rock samples, together with the crystallographic preferred orientation and the single crystal elastic constants, using the methods of Mainprice (1990). It is shown that the intrinsic anisotropy of peridotites is mainly due to the lattice preferred orientation (LPO) of olivine and to the competing effects due to the LPO of orthopyroxene. If present, well-oriented hydrous phases such as phlogopite, leads to both P- and S-waves velocities becoming drastically reduced, but the anisotropy and the S-wave splitting (mainly on the foliation plane) is increased. An attempt has been made to apply this petrologic and microstructural approach to understand the seismic properties of ultramafic rocks with the aid of a mobile seismic network set up in the Ivrea-Verbano zone. Five recording units with 3-component geophones were installed on the most representative lithologies of the Ivrea-Verbano zone: kinzigites, gabbros, stronalites, and on metapelites of the adjacent outcrop of the Serie dei Laghi. During five months the seismic network recorded 49 events that were mainly regional earthquakes. The preliminary analysis of the first arrival revealed a high-velocity anomaly near the Finero ultramafic massif. Rough calculations indicated that this anomaly could be due to a shallow-depth, peridotitic body about six kilometres thick. S-wave polarisation analysis was carried out on the horizontal components of the seismograms of four regional earthquakes that showed a first arrival propagating almost vertically at the recording sites. Four earthquakes were chosen in order to reconstruct the orientation of the polarisation plane of the shear phases. The preliminary results of this analysis have shown a qualitatively good agreement between the orientation of the polarisation plane of the fast shear-wave and the strike of the foliation at the station near Finero, located on kinzigites of the Ivrea-Verbano zone, although the relationships between the local structural setting and the polarisation orientations at the other recording sites is not yet so clear.     

Preliminary analysis of crustal shear-wave splitting in the Sanjiang lateral collision zone of the southeast margin of the Tibetan Plateau and its tectonic implications

The collision of the Indian and Eurasian plates, to the east of the eastern Himalayan syntaxes, forms the Sanjiang lateral collision zone in the southeast margin of the Tibetan Plateau, where there are intense crustal deformation, active faults, earthquakes, as well as a metallogenic belt. Given the lack of adequate seismic data, shear-wave splitting in this area has not been studied. With seismic data from a temporary seismic linear array, as well as permanent seismic stations, this paper adopts the identification on microseismic event to pick more events and obtains shear-wave splitting parameters from local earthquakes. From the west to the east, the study area can be divided into three subzones. The “fast” polarization (i.e. the polarization of the fast shear wave) varies gradually from NNW to NS to NNE in these three subzones. The time delay of the slow shear wave (i.e. the time difference between the two split shear waves) also increases in the same direction, indicating the presence of seismic anisotropy above 25 km in the crust. Both shear-wave splitting parameters are closely related to stress, faults and tectonics. The scatter and the “dual” (i.e. two) dominant orientations of the fast polarizations at several stations indicate strong distortions caused by nearby faults or deep tectonics. The anisotropic parameters are found to be related to some degree to the metallogenic belt. It is worth to further analyse the link between the anisotropic pattern and the metallogenic area, which suggests that shear-wave splitting could be applied to study metallogeny. This paper demonstrates that the identification on microseismic event is a useful tool in detecting shear-wave splitting details and exploring its tectonic implications.     

按方位叠加接收函数分析青藏高原东南缘的地壳各向异性

在各向异性地壳中,来自Moho的P-to-S转换波(Pms)的到时不仅取决于入射角和地壳厚度,而且还随地震事件方位角而变化.地处青藏高原东南缘的川滇地区,地壳变形十分强烈.本文利用川滇地区的108个固定台站记录的远震三分量地震波形数据提取台站下方的P波接收函数,并把接收函数被校正到了同一参考震中距处(例如67°).然后按后方位角10°为间隔将接收函数叠加成一道信号以增强信噪比,并从叠加信号里拾取不同后方位角对应的Pms相的观测到时.在快波极化方向和分裂时间构成的解的平面上,能使观测到时与理论到时之差最小的点即为所求的分裂参数的位置.合成地震图和实际观测数据的实验表明,这个方法不但稳定性较好,而且误差估计也较小.我们从108个台中获得了96个Pms相的分裂参数,结果表明,川滇地区地壳各向异性十分强烈,Pms相分裂时间在0.05s±0.06s到1.27s±0.10s之间,平均值为0.54s±0.12s.地壳各向异性的快波极化方向与地表GPS速度场的差异性表明,印支块体的上下地壳之间是解耦的,而川滇菱形块体北部、松藩—甘孜和四川盆地的上下地壳之间是耦合的.然而,川滇菱形块体南部,地壳变形主要受控于小江断裂和金沙江—红河断裂.     

地壳剪切波分裂研究及地震定位误差对剪切波分裂分析的影响

地壳介质的各向异性主要是由大量充满液体的微裂隙的定向排列引起的, 剪切波穿过这种含有微裂隙的介质时会产生分裂。 剪切波分裂参数与地壳介质的物理性质有关, 与区域应力场有关。 在地壳剪切波分裂研究中, 地震的定位误差对剪切波分裂参数有影响。 本文简单介绍了几种影响定位精度的因素; 并对定位误差特别是深度方向的误差所引起的时间延迟变化进行了讨论, 比较了不同大小的深度误差对归一化后的时间延迟的影响, 定量分析结果显示, 深度定位误差对慢剪切波时间延迟影响很大。 剪切波分裂分析方法有多种, 本文还对不同的剪切波分裂分析方法做了简介。     

基于地壳介质各向异性分析青藏高原东北缘构造应力特征

青藏高原东北缘由于受到多个构造块体的共同约束,表现出复杂的地球物理特性和地质特性,本文利用甘肃数字地震台网(2001-2008年)的观测资料,采用系统分析方法(SAM),进行地壳剪切波分裂分析,获得研究区内18个台站共1005条记录的剪切波分裂参数.研究结果表明,青藏高原东北缘介质各向异性在空间上存在差异,慢剪切波延迟时间表明了地壳介质各向异性的强弱变化特征,快剪切波平均偏振方向则反映了本区区域构造应力的空间变化特征.分析认为,祁连山-河西走廊活动构造区直接受青藏地块与阿拉善地块间相互作用,与青藏地块构造应力一致;甘东南活动构造区的应力环境主要受到内部活动断裂的共同作用,具有局部构造应力的特征.     

青藏高原东北缘中上地壳介质各向异性及其构造意义

青藏高原东北缘(94°E—105°E,32°N-40°N)是高原北东向扩张的前沿地带,亦是研究高原生长过程的重要区域.本文利用青海省数字地震台网(2008-2014年)共7年的地震目录和波形数据,首先使用双差定位获取精定位震源位置,在此基础上,挑选位于S波窗口内(射线入射角≤45°)的地震事件,依据S波分裂分析方法(SAM),获取研究区域内共26个台站的S波分裂参数.研究结果表明:地处多个块体交汇部位的西宁及其周缘,地壳各向异性呈现两个优势偏振方向,表明该区中上地壳应力环境由区域主压应力场和活动断层共同约束;玉树地震序列的地壳各向异性优势偏振方向与区域主压应力场一致.     

CRUSTAL ANISOTROPY AND ITS TECTONIC IMPLICATIONS IN THE CHONGQING REGION

The receiver function which carries the information of crustal materials is often used to study the shear-wave velocity of the crust as well as the crustal anisotropy. However, because of the low signal-to-noise ratio in Pms(P-to-S converted phase from the Moho), the crustal anisotropy obtained by shear-wave splitting technique for a single receiver function usually has large errors in general. Recent advance in the analysis method based on Pms arrival time varying with the back-azimuth change can effectively overcome the above defects. Thus in this paper, we utilize the azimuth variations of the Pms to study the crustal anisotropy in Chongqing region for the first time. According to the earthquake catalogue provided by USGS, seismic waveform of earthquakes with magnitude larger than 5.5 and epicenter distance range of 30°~90° between January 2015 and December 2016 are collected from 14 broadband seismic stations of Chongqing seismic network. We carry out the bootstrap resampling to test the reliability of the radial maximum energy method for the observation data. In addition, we also applied the receiver function H-Kappa analysis in this paper to study the crustal thickness and Poisson's ratio. Our results show the crustal thickness ranges from 40~50km, and there is a thin and thick crust in the southern and northern Chongqing, respectively. The crustal average Poisson's ratio ranges from 0.23~0.31, the Poisson's ratio reaches the maximum value in the central part of Chongqing, while the Poisson's ratio in the northern and southern parts of Chongqing is obviously low. We obtain the crustal anisotropy from 9 stations in total. The delay time of crustal anisotropy distributes between 0.08s and 0.48s, with the average value of 0.22s. Among them, the CHS, QIJ and WAZ stations in central Chongqing have relatively large crustal delay time(>0.3s), followed by ROC station in the western Chongqing(0.25s), while the delay time in CHK station in northern Chongqing and WAS station in southern Chongqing are 0.08s, showing relatively weak crustal anisotropy. The fast polarization directions(FPDs)also change obviously from south to north. In southern Chongqing, FPDs are dominant in NNE-SSW and NEE-SWW, while the FPDs in WAZ station change to NWW-SEE, and the FPDs appear to be NW-SE in CHK in the northern Chongqing. In general, the FPDs are sub-parallel to the strikes of faults in most areas of Chongqing areas. Combined with other results from GPS observations, tectonic stress field and XKS splitting measurements, the main conclusions can be suggested as following:The cracks preferred orientation in the upper crust is not the main source of crustal anisotropy in Chongqing area. The crust and lithospheric upper mantle in the eastern Sichuan fold belt(ESFB)and Sichuan-Guizhou fault fold belt(SGFFB)are decoupled, and the deformation characteristics in the north and south parts of ESFB and SGFFB is different. The complex tectonic deformation may exist beneath the mountain-basin boundary, causing the fast directions of crustal anisotropy different from that in other areas of ESFB and SGFFB. The faults with different strikes may weaken the strength of average crustal anisotropy in some areas. The crustal deformation in southern Dabashan nappe belt(DNB)may be mainly controlled by the fault structure.     

A model for lower continental crust

This proposed model is based on geological, geophysical and geochemical data. Previous models suggested for the lower continental crust consisted of basalt, gabbro, or charnockitic rocks; however, experimental and field petrological data indicate that the bulk of crustal rocks are metamorphic. A lower crust of heterogeneous metamorphic rocks also agrees with seismic reflection results which show numerous reflections from “layering”. Geothermal conditions favor a “dry” charnockitic or gabbroic lower crust rather than an amphibolitic lower crust because heat production data imply that wet amphibolitic rocks would have a higher heat production than their dry metamorphic equivalents. Relatively high velocities from field and laboratory measurements in such low-density rocks as granite, syenite, anorthosite and granulitic rocks in general imply that the composition of the lower crust is more felsic than gabbro. Variation in seismic velocity and depths from crustal refraction studies and numerous seismic reflections all indicate a highly heterogeneous lower crust. The lower crust, which has traditionally been described as gabbroic or mafic, may consist of such diverse rocks as granite gneiss, syenite gneiss, anorthosite, pyroxene granulite, and amphibolite, interlayered on a small scale, deformed, and intruded by granite and gabbro. Interlayering of these rocks explains the presence and character of seismic reflections. Abrupt changes in dip, tight folding, disruption of layers, intrusion, and changes in layer thickness explain the characteristic discontinuity of deep reflections. Igneous intrusions may be floored by metamorphic rocks. The lower crust consists of a complex series of igneous and metamorphic rock of approximate intermediate composition.     

Shear-wave splitting beneath Yunnan area of Southwest China

Systematic analyses of seismic data recorded by the Yunnan regional seismograph network reveal significant crustal and upper mantle anisotropy. Splitting of the S phase of local earthquakes and teleseismic SKS, PKS, and SKKS phases indicates time-delays from 1.60 ms/km to 2.30 ms/km in the crust, and from 0.55 s to 1.65 s in the upper mantle which corresponds to an anisotropic layer with a thickness about between 55–165 km. The polarization orientations of fast shear waves in the crust are complicated with a predominantly north-south direction, and the mantle anisotropy has a nearly west-east direction. Our results show different deformation styles and mechanisms exist between the crust and upper mantle.     

首都圈地区远震基底PS波分裂研究

远震PS波的一个重要优势是其携带台站下方介质的各向异性信息.本文基于首都圈数字地震台网2002-2003年记录的高精度远震波形资料,从基底PS转换波入手,采用三分量PS波分裂偏振分析方法,获取各个台站下方来自基底的PS波分裂参数,包括快波偏振方向和快、慢波的时间延迟等,分析台站下方浅部地壳介质的各向异性特征.几个台站的研究结果显示,该方法保持了原始波场的信息,可以用来分析PS波分裂特征;短周期台站记录的远震来自基底的PS波具有干扰少、信噪比高的特点,存在分裂现象,各台站快慢波时差超出我们的预期值,平均在0.1~0.2 s左右.此研究可与横波分裂相对比,对于认识地壳各向异性及其内部的应力状态,分析其构造及地震活动现象有较大的意义.