地球物理学报  2010, Vol. 53 Issue (6): 1374-1383   PDF    
青海地区S波分裂研究
李永华 , 吴庆举 , 冯强强 , 张瑞青     
中国地震局地球物理研究所, 北京 100081
摘要: 本文利用国家地震台网及中国地震局“十五”期间在青海布设的30个宽频带地震台站记录到的远震波形数据,分别采用最小能量法和旋转相关法对SKS、SKKS和PKS波震相进行了偏振分析,计算了台站下方介质的各向异性分裂参数:快波的偏振方向(φ)和慢波延迟时间(δt).本文研究结果表明,研究区多数台站下方的地震各向异性参数都表现出随方位角变化而变化的特征,可以用双层各向异性模型来解释.其中上层各向异性的快波偏振方向位于N65°E~N95°E之间,可能与中下地壳物质的流动有关;而下层各向异性的快波偏振方向位于N105°E~N135°E之间,可能为祁连块体NEE向的推移及导致的岩石圈缩短有关.此外,我们还发现,与周边的台站下方各向异性分裂参数相比,昆仑断裂附近两个台站GOM和DAW的各向异性特征急剧变化,其快波方向都与该断层近乎平行,这很可能暗示昆仑断裂已经切穿整个岩石圈;阿尔金断裂附近两个台站(LEH和HTG)无效分裂事件的方位分布与阿尔金断裂走向缺乏相关性,我们推测研究区内阿尔金断裂可能为地壳尺度的断裂.
关键词: S波分裂      各向异性      壳幔变     
Seismic anisotropy beneath Qinghai province revealed by shear wave splitting
LI Yong-Hua, WU Qiang-Ju, FENG Qiang-Qiang, ZHANG Rui-Qing     
Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
Abstract: Shear wave splitting analyses have been carried out using teleseismic data from 30 broad-band seismograph stations deployed at permanent locations in Qinghai province, and we determined fast polarization directions φ, and delay times δt between fast and slow components using both the minimum energy method and rotation-correlation method. Our shear wave splitting measurements show that there are variations with propagation directions beneath most stations. Analysis of individual measurements from multi-station highlights variations of splitting parameters with event backazimuth and can be best explained by a two-layer anisotropic model with horizontal symmetry axes. The optimal fast direction for the upper layer is in the 65°~95° range which is related to crust flow in this region. The fast direction for the lower layer ranges from 105° to 135° and can consequently be attributed to the ongoing NEE direction lithospheric shortening. Our results also show that the fast directions beneath stations GOM and DWA are parallel with the trend of Kunlun fault. The correlation between anisotropy and the fault strike suggests that the Kunlun fault must be deep and probably extends into the mantle. At stations (HTG and LEH) lying close to Altyn-Tagh fault, we can not determine any splitting measurements and only observe many Nulls whose event backazimuths are neither parallel nor normal to the strike of Altyn-Tagh fault, suggesting that in this region the fault is a crustal scale structure..
Key words: Shear wave splitting      Anisotropy      Crust-mantle deformation     
1 引言

青海地处青藏高原东北缘,是晚新生代构造变形最强烈的地区之一(图 1)遍布全区的第四纪褶皱、断层[1~3]及地震震源机制[4, 5]等研究表明,该地区正遭受着地壳缩短,并且伴随着垂直隆升作用;地震剖面、地震活动性和钻孔资料等研究表明,柴达木盆地的褶皱和逆断层形成于白垩纪,上新世活动最为活跃,直到现今仍处于活动状态[1~5].在该区开展深部结构与变形研究是认识高原隆升机制和印度-欧亚板块碰撞远程效应的重要途径之一.

图 1 青海地区构造分区和地震台站分布 ①班公-怒江缝合带;②金沙江缝合带;③昆仑断裂带;④北祁连山逆冲带;⑤阿尔金断裂带;Ⅰ羌塘地块;Ⅱ巴颜喀拉地块;Ⅲ柴达木地块;Ⅳ祁连地块. Fig. 1 The distribution of the seismic stations in Qinghai province, and the solid circles represent the broadband seismic stations ①Bangong-Nujiang suture zone; ②Jinsha suture zone; ③Kunlun fault; ④North Qilian thrust belt; ⑤Altyn Tagh fault; Ⅰ Qaingtang block; Ⅱ Bayan Har block; Ⅲ Qaidam basin; Ⅳ Qilian block.

20世纪80年代以来,在青藏高原东北缘开展的地质调查、人工地震探测、地震台阵观测,为青藏高原东北缘深部构造研究奠定了基础[6~23]然而,关于青藏高原东北缘的壳幔变形机制仍存在争议.如Tapponnier等[3]推断印度-欧亚板块碰撞挤压而导致脆性的上地壳沿着走滑线侧向挤出,导致地表产生了与地壳缩短方向相垂直的走滑断裂系.根据该假说,边界断裂带(如,阿尔金断裂和昆仑断裂)变形异常强烈,而断裂带周边较宽的地段也会产生强的地幔变形.人工地震探测与地震台阵观测得到的地震波速度结构研究显示,研究区具有低速的中下地壳[7, 8, 10~18],似乎可以用Royden等[24]提出的“下地壳流模式”来解释.由于该模型中上地壳和岩石圈地幔之间存在一韧性层,这就意味着作用在地壳上的力就不能传递到地幔,因而上地壳和岩石圈地幔之间是解耦的.

在青藏高原东北缘这一强变形区,开展壳幔深、浅层的变形关系研究对于研究区变形机制的确立具有非常重要的意义.前人[6, 9, 19~23]曾在青海及邻区开展了S波分裂研究,部分学者[19, 20, 25, 26]通过将观测到的S波分裂快波方向与地表变形场的最大剪切方向作比较,推断青藏高原内部壳幔变形为垂直连贯变形,青藏高原外部的云南地区壳幔变形属于壳幔解耦模式,然而受台站分布限制,并没有详细讨论青藏高原东北缘的可能变形方式.中国地震局“十五”期间在青海地区布设了20余台宽频带地震台站,到目前为止已经积累了大量的观测资料[7].为此,本文利用国家地震台网及中国地震局“十五”期间在青海布设的30个宽频带地震台站记录到的远震波形数据(图 1),分别采用最小能量法[8]和旋转相关法[29]开展S波分裂研究.最后我们还结合已有的地质、地球物理证据对本文的地震各向异性分析结果及其可能的地质含义进行了讨论.

2 数据与方法

本文研究收集了青海地震台网(2007.07~2010.03)、中国国家数字地震台网(2001.01~2010.03)和中国数字地震台网(1990.01~2010.03)共30个宽频带地震台站(图 1)的远震记录.为利用远震记录中的SKS、SKKS和PKS震相(本文以后通称XKS震相)开展偏振分析,我们选取了震中距大于85°震级Mw大于5. 8级的地震事件.为提高远震记录中XKS震相的信噪比,在计算各向异性参数前对数据进行了带通滤波处理.多数情况下,我们选择0.02~1 Hz的滤波器,少数情况下( < 5%),我们也采用了0.04~0.5 Hz的滤波器.

研究中利用Splitlab软件[30],同时使用最小能量法[28]和旋转相关法[29]对单个远震记录的XKS震相开展偏振分析,求取各向异性参数对(φδt)其中φ表示快波偏振方向,δt表示快、慢波的到时差. 图 2为利用台站BAM记录的远震事件,通过对SKS震相进行偏振分析的一个例子.由图 2可以看出,同时使用多种方法开展S波分裂研究,求取的各向异性参数较任何一种单一方法得到的结果都要可靠.此外,采用多种方法用于S波偏振分析,还有助于弱各向异性或无效事件(null)的确定[31].所谓的无效测量事件指不产生分裂的地震记录,无效测量事件的出现可能与台站下方各向同性的介质相关,也可能是由于该事件的方位分布平行或垂直于快波偏振方向所致.

图 2 采用Splitlab软件[30]对BAM台远震记录的SKS震相开展S波偏振分析的实例 (a)左侧为观测到的地震波形(图中径向、切向分量分别用虚线和实线表示, 虚线为依据IASP91模型计算得到的SKS震相理论到时, 阴影部分表示用于S波分裂计算的时间窗);中间部分为波形数据信息及用三种方法的计算结果;右侧为S波分裂结果的立体投影;(b), (c)分别表示采用旋转相关法和最小切向能量法开展S波分裂测量的例子, 从左到右(1)表示经时移校正后的快波(虚线)、慢(实线)波分量;(2)各向异性校正后的径向(虚线)、切向(实线)分量;(3)各向异性校正前(虚线)、后(实线)的质点运动轨迹;(4)各向异性参数分布等值线图. Fig. 2 Example of a SKS splitting measurement using the Splitlab package[30]at BAM station (a) The left-hand panel shows the observed seismograms (dashed line, radial component; solid line, transverse component; vertical dashed lines, the phase arrival times computed using IASP91 model; gray zone, calculation window).The right-hand panel presents a stereoplot of the result.The header shows the information of teleseismic event as well as splitting parameters resulting from the three techniques.(b) and (c) display the results from the Rotation-Correlation technique and the minimum energy, respectively.From left-hand side to right-hand side:(1) Fast (dashed) and slow (solid) components, corrected for the calculated splitting delay time; (2) Corrected radial (dashed) and transverse (solid) components; (3) Particle motion before (dashed) and after (solid) correction; (4) Contour plot for the maximum value of correlation coefficient and for the energy on transverse component as function of delay time and fast polarization angle.

根据Barruol等[32]和WUstefeld等[31]关于分裂测量的判别方法,我们根据原始地震记录中XKS震相的信噪比、快慢波校正前后的波形及质点运动轨迹及不同测量方法(旋转相关法和最小能量法)得到的S波分裂参数等,将分裂结果分为高质量有效分裂结果(good)、一般(fair)、不可靠(poor)、高质量的无效分裂结果(good Null)和质量一般的无效分裂结果(fairNull)五类.一般认为,目前的S波分裂测量方法中,利用最小能量法测量得到的结果相对更为可靠、稳定[33, 34],因此在下文中关于有效分裂结果的讨论,以利用最小能量法测量得到的结果为准.

3 剪切波分裂结果

本文根据前述分析方法,对研究区30个台站的记录开展了S波分裂计算,最后得到220个(good+ fair)高质量的S波分裂结果(图 3)和216个高质量的无效分裂测量(图 4),其中用于测量的SKS、SKKS和PKS震相分别为38121和34个.遗憾的是,在YUS和TTH台没有得到可靠的分析结果(图 1图 3, 图 4).

图 3 青海地区S波分裂测量结果 线段方向代表的是快波的偏振方向,长度代表的是慢波的延迟时间. Fig. 3 Summary of splitting measurements at each station in this study, where the lines represent φ (the azimuth of the fast splt shear wave) and their lengths are proportional to the delay tmes δt
图 4 青海地震台站无效分裂地震方位分布, 说明同图 3 Fig. 4 Distribution of null measurements at seismic stations

图 3, 4可以看出,研究区多数台站下方的地震各向异性参数都表现出随方位角变化而变化的特征.这就意味着,如果我们将研究区所有台站下方的各向异性进行平均,将可能会导致对该区的变形特征分析产生严重的偏差和误导.为此,本文研究并没有试图给出台站下方的平均各向异性.

尽管研究区不同构造带地震台站下方各向异性强度和快波方向横向变化明显,但是相邻地震台站的分裂结果具有很好的相似性(图 3, 4).我们的分析结果显示:(1)青海西部的一些台站(如,BAM,GOM,DAW)表现为简单的分裂模式,其台站下方的剪切波视分裂参数不随方位变化而变化,且无效事件的方位分布与台站下方的快波偏振方向垂直或平行,其中GOM台的快波偏振方向呈近W-E向分布;而其他台站的快波偏振方向呈NW-SE或NNW-SSE向分布,即台站下方的各向异性快波偏振方向与地表断裂走向相一致,这与前人[6, 9, 19~23]的研究结果相一致;(2)相对而言,青海东部多数台站(如DUL,TOR,HUL,LJS,LYX,MAD, TIJ,XIN,XUH,LWS,QIL,QSS,DAT,MEY,QML,DCD,DLH,LED)下方的各向异性参数并不一致,这一现象在常利军等[21, 22]的研究结果中也有显示.我们的进一步研究显示这些台站下方的各向异性参数并非随机分布,而是随方位角变化而规律地变化(图 5).受观测事件方位分布的限制,部分台站(如DUL、DLH)的各向异性参数虽然没有表现出这么明显方位相关性,但是其快波偏振方向与无效事件的方位分布并无关联(图 34),这些特征都暗示台站下方具有复杂的各向异性结构.事实上,在青藏高原腹地及其边缘地带也有各向异性参数随方位变化的例子[34~38]; (3)少数台站(如ZAD,XIH)只观测到1个有效分裂,而在HTG和LEH,我们只得到了一些无效分裂事件的方位分布.对于这4个台站,我们只能参照其周边台站的各向异性特征给以判别.如,XIH台与周边台站(如LYX,TOR等)具有相似的各向异性特征;但由于ZAD、HTG和LEH这3个台站与周边其他台站相距较远,我们无法判断其台站下方的各向异性特征表现为简单的分裂方式还是复杂的方式.

4 讨论 4.1 各向异性的深度分布

在对S波分裂测量结果进行解释之前,我们首先必须回答各向异性的深度分布问题.由于分析中所使用的XKS震相的射线路径近似垂直,单独利用XK S测量无法探测台站下方各向异性的确切深度分布情况,但是我们可以根据S波分裂测量的特点及其他地球物理探测,定性地对各向异性的深度分布给出约束.(1)由S波分裂测量得到的研究区台站下方快慢波到时差一般都大于0.7 s,尽管接收函数研究显示该区存在强的壳内各向异性层[12, 39],但青藏高原地区地壳各向异性导致的快慢波到时差一般介于0. 1~0. 3 s之间[6, 40],因此本文测量得到的S波分裂延时主要来自于地幔.(2)尽管已有面波各向异性研究[41~43]关于不同周期面波方位各向异性的快波方向不很一致,但是都显示面波方位各向异性主要分布在20~100 s的周期范围内.(3)通过计算SK (K) S震相的菲涅耳带,我们能对地震各向异性的深度分布做出估计,如DAW和XIH台站间的距离约为110 km,但是其台站下方的各向异性特征明显不同;依据Alsina and Snieder[44]的方法,通过计算SK (K) S震相的菲涅耳带,我们能对地震各向异性的深度分布做出估计.为简单起见,我们取SKS震相的主频为8s,这样我们可以粗略地估计出100, 200, 300 km深处相对应的菲涅耳带直径分别为~80, 120, 140 km.因此,我们推断研究区各向异性的深度分布不超过200 km. (4)假定造成XKS分裂过程中快慢波的延迟由地幔主要造岩矿物晶格定向排列引起,那么1 s的快慢到时差大概相当于100 km厚度的各向异性层.研究区XKS分裂延迟多数介于0.4~1.5 s,由此估算出的该区各向异性层厚度范围是40~150 km,这也与前述推论相一致.

4.2 各向异性随深度的变化

前述已经提及,研究区部分台站(尤其是青海东部地区)下方的地震各向异性参数都表现出随方位角变化而变化的特征,台下方地震波各向异性分裂参数的这种变化,暗示使用单层各向异性模型来解释其地震各向异性结构是不合理的.研究人员一般将这种变化归结为各向异性随深度变化(如双层各向异性[45]、倾斜各向异性对称轴[46]或者各向异性横向不均一)所致[21, 22, 35].然而,受地震事件分布所限,我们无法逐一对每个台进行分析,以确定这种变化的原因.本文研究中,我们通过观察每个台站的S波分裂参数及其变化,发现台站DUL,TOR,HUL,LJS,LYX,MAD,TIJ,XIH,XIN,XUH,LWS,QIL,QSS,LED,DAT和MEY具有相似分裂特征(图 3~5),为此可将上述16个台站结果看作一个台站下方的观测结果.如图 5所示,视分裂参数并不是杂乱无章的分布,尤其在后方位角108°~166°范围内,快波偏振方向在N54°E~N119°E之间呈规律性变化,而分裂延迟在0.4~1.6 s之间表现出逐渐增大的趋势,其周期为90°为此我们将测试是否可以用双层水平对称轴介质各向异性模型解释本文的观测结果.

图 5 快波偏振方向(a)、快慢波延迟(b)与方位角的关系 其中的视分裂参数为对西宁附近16个台站进行偏振分析得到的结果,虚线为平均各向异性;实线为根据最优双层各向异性模型合成得到的各向异性分裂参数(Radjusted2=0.55;φlow=N115 E,δtlow=0.9 s,φupp=N72° E,δtupp=0.5 s).垂直阴影部分为无效事件的分布范围. Fig. 5 Plots of the apparent fast polarization direction (a) and delay time (b) as a function of BAZ The measurements are from sixteen stations (detaH in text). Apparent splitting parameters predicted for the optimum two-layer model (Radjusted2=0.55;φlow=N115 E, δtlow=0.9 s, φupp=N72° E, δtupp=0.5 s) are plotted for frequencies of 0. 125 Hz (solid line). Dotted horizontal lines show the mean splitting parameters. Light grey vertical bars illustrate the BAZ distribution of nulls.

以下研究中,我们仅将有效分裂观测用于双层各向异性模型参数的确定.假定XKS震相的主频为8s,我们让上、下两层各向异性的快波方向别以1°的步长在0°~180°之间变化,分裂延迟以0. 1 s的步长在0~1.6 s之间变化(1. 6 s为所有视分裂延迟的最大值).首先采用Silver和Savage[45]提出的方法合成不同双层各向异性模型导致的S波视分裂参数.然后,按照Walker等[47]和Fontaine等[48]描述的方法,对比计算合成模型与实际分裂观测值的误差及拟合程度.其中残差最小的模型为φupp=N72° E,δtupp=0.5 s; φlow=N115 E,δtlow=0.9 s.由图 5可以看出,该模型可以很好地解释本文的S波视分裂参数(Radjusted2=0.55, 即拟合度约为55%),由该模型预测的无效分裂事件分布方位为N100°E,这也与上述台站的无效分裂事件分布范围相一致.

已有研究[34, 37, 45~48]表明,仅仅依赖于视分裂参数无法确定唯一的双层各向异性模型,当视分裂参数方位分布不均匀时更是如此,为此给出模型的误差范围显得尤为必要.根据Walker等[47]的方法,我们从前述8294400个模型中选出了2500个分布于95%置信区间的最优模型(Radjusted2>0.47)(图 6),其上、下两层各向异性的快波方向分别位于N22°E-N98°E和N103°E-N177°E; 而上、下两层的分裂延时都分别位于0.4~1.6 s之间,即上、下两层的各向异性快波方向误差为±35°,分裂延时的误差为±0. 6 s,这意味着我们无法很好地约束上、下两层的各向异性参数.但我们注意到,其中~70%模型的上、下两层各向异性快波方向分别集中分布于N65°E~N95°E和N105°E~N135°E之间,而上、下两层的分裂延时都分别位于0. 4~0. 9 s之间.为此,我们认为上述台站下方各向异性分裂参数可以用双层各向异性模型来解释,其中上、下层各向异性的快波方向分别为NEE向和SEE向,如图 5图 6所示.尽管我们没法得到唯一的双层各向异性模型,但与单层各向异性模型相比较,我们提出的双层各向异性模型能更好地解释该区的各向异性参数分布特征.

图 6 双层各向异性模型中上、下层各向异性分裂参数统计分布 其中2500个最优模型(Radjusted2>0.47)均分布于95%置信区间, 直方图(a)、(b)、(c)、(d)的横坐标分别为上层、下层各向异性分裂参数; 纵坐标为模型个数. Fig. 6 Distribution of the four splitting parameters for the two-layer models (with misfit reduction Radjusted2 > 0.47) included in the 95 per cent confidence region from multi-station two-layer modeling approach. (a) Fast direction in the upper layer, (b) Delay time in the upper layer, (c) Fast direction in the lower layer; (d) Delay time in the lower layer.

我们发现QML、DCD、DUL台也具有复杂的各向异性特征,但由于分裂事件有限,我们采用上述方法没能很好地对这两个台站下方的各向异性特征给出约束.我们的结果表明,基于上述台站观测得到的双层各向异性模型并不能很好地解释QML、DCD、DUL台站的S波分裂观测结果,这也暗示该区的各向异性特征远比我们想象的要复杂的多.

4.3 各向异性的地球动力学意义

通常认为,地幔各向异性是由于地幔变形导致地幔主要造岩矿物(如橄榄岩、斜方辉石)晶格定向排列引起的,但关于导致地幔主要造岩矿物晶格产生定向排列的动力学机制还存在争议.Silver[49] (1996)提出可以采用简单软流圈流动模型和垂直连贯变形模型两种端元模型来解释大陆地区的地幔各向异性.已有研究[19, 20, 25, 26]通过对3波分裂得到的快波方向与GPS、第四纪断裂活动速率推断的地表变形场作比较,推断青藏高原内部壳幔变形为垂直连贯变形,青藏高原外部的云南地区壳幔变形属于壳幔解耦模式,但是受台站分布限制,他们并没有详细讨论青藏高原东北缘的变形方式,从Wang等[20]的结果来看,研究区的壳幔变形方式似乎应属于壳幔解耦一壳幔耦合的过渡地带,但壳幔解耦一壳幔耦合过渡地带的变形方式究竟表现如何尚不清楚.

前述提及青海东部西宁(XIN)附近16个台站的S波视分裂参数可以用双层模型来解释,其中上、下层各向异性的快波偏振方向分别位于N65°E~N95°E和N105°E~N135°E之间.其中上层各向异性的快波方向与该区的最大主压应力方向[45]、GPS运动方向[50]一致.由于研究区的缝合带和断裂带多呈SEE-SE向展布,因此,我们不能将上层各向异性归结为地壳的裂隙或线理构造所致.大量的天然地震和人工地震测深研究显示,该区存在显著的中下地壳低速层[7, 8, 10~18],壳内低速物质的定向流动也能产生强的地壳各向异性,因此,我们推断该层各向异性可能与中下地壳物质的流动有关.而下层各向异性的快波偏振方向与祁连山的走向近乎一致,据此我们推断下层各向异性的形成可能与祁连块体NEE向的推移及导致的岩石圈缩短有关?

我们知道,仅仅依赖于XKS测量和模拟无法解决各向异性的深度分布问题,因此,关于上、下两层各向异性的深度分布和形成原因的判别还需要借助其他地球物理探测结果的约束.Shaprio等[41]的Rayleigh波群速度方位各向异性研究结果显示,青藏高原东北缘周期20 s与50 s的快波方向分别为NEE向和SEE向,这也说明青藏高原东北缘上地壳、下地壳及上地幔的运动存在一定程度的解耦,这一结果也再次印证了本文前述的推论.事实上,这一推论绝不仅仅出现在青藏高原的边缘地带,如Gao和Liu[34]的结果显示,拉萨块体也具有类似的复杂各向异性特征,并不能用垂直连贯变形模型来解释.

研究区的昆仑走滑断裂和阿尔金走滑断裂是研究区内两条非常重要的断裂带,这两条断裂带究竟属于地壳断裂还是岩石圈断裂,与研究区壳幔解耦一耦合变形关系密切[1, 3, 24].本文研究发现,与周边的台站下方各向异性分裂参数相比,昆仑断裂附近两个台站GOM,DAW的各向异性特征急剧变化,其快波方向都与该断层近乎平行.事实上,已有研究[6, 9, 19~22, 51]也显示昆仑断裂与其附近台站下方各向异性快波方向具有很强的相关性,这很可能暗示昆仑断裂已经切穿整个岩石圈[51].

很遗憾,阿尔金断裂附近两个台站(LEH和HTG)没有S波分裂结果,但是两个台站无效分裂事件的方位分布与阿尔金断裂走向既不平行也不垂直,因此,我们无法支持Witlmger等[7]关于阿尔金断裂中断是岩石圈断裂这一推论.与之形成鲜明对比的是,姜枚等[9]的研究结果揭示,90°E以西阿尔金断裂附近地震台站的各向异性快波方向与该断裂近乎平行.阿尔金断裂不同地段S波分裂快波方向的差异,可能暗示阿尔金东、西段穿切深度不同.地质调查表明,在阿尔金断裂带西部的阿什库勒盆地与普鲁附近有新生代玄武质火山岩的出露[2],但在且末以东的阿尔金断裂东段却没有新生代玄武质火山岩喷发的痕迹,这也暗示阿尔金断裂西段可能为岩石圈断裂,而阿尔金断裂东段则为地壳断裂[2],但我们推测其分界线可能并不在且末(~85°E)附近,而在~90°E附近,这较Yn等2所推测位置更靠东一些.

5 结论

本文利用国家地震台网及中国地震局“十五”期间在青海布设的30个宽频带地震台站记录到的远震波形数据,分别采用最小能量法和旋转相关法对XKS波震相进行了偏振分析.本文研究结果显示,研究区多数台站下方的地震各向异性参数都表现出随方位角变化而变化的特征,不能用单层各向异性模型来解释.本文研究测试表明,尽管我们不能给出惟一的双层各向异性模型,但与单层各向异性模型相比较,我们提出的双层各向异性模型能更好地解释该区的各向异性参数分布特征.其中上层各向异性的快波偏振方向位于N65°E~N95°E之间,可能与中下地壳物质的流动有关;而下层各向异性的快波偏振方向位于N105°E~N135°E之间,可能为祁连块体NEE向的推移及导致的岩石圈缩短有关,这也暗示研究区的上地壳、下地壳及上地幔的运动存在一定程度的解耦.

本文结果还显示,与周边的台站下方各向异性分裂参数相比,昆仑断裂附近两个台站GOM,DAW的各向异性特征急剧变化,其台站下方的剪切波视分裂参数不随方位变化而变化,且快波方向都与昆仑断裂带近乎平行,暗示昆仑断裂已经切穿整个岩石圈.这很可能意味着昆仑断裂周边地带的壳幔变形很可能属于强耦合模式.但阿尔金断裂附近两个台站(LEH和HTG)无效分裂事件的方位分布与阿尔金断裂走向并不相关.结合已有的S波分裂结果和地质调查研究,我们初步推断阿尔金断裂西段(~85°E以西)属于超壳断裂,而阿尔金断裂东段(~90°E以东)则为地壳尺度的断裂.

综上所述,研究区可能属于壳幔解耦一壳幔耦合过渡地带,仅用壳幔解耦或者壳幔耦合都不能很好地解释该区的S波分裂观测,其变形方式的确立有待于更多资料和更多地球物理观测的约束和证实.

致谢

感谢中国地震局地球物理研究所“国家数字测震台网数据备份中心”和中国数字地震台网为本研究提供地震波形数据.在资料的搜集过程中,郑秀芬博士提供了无私的帮助.本文使用了Wiistefeld的软件splitlab, 在该软件的使用过程中Wiistefeld给予了耐心的指导和帮助;本文成文过程中与王椿镛研究员的请教与探讨获益匪浅;两位审稿人提出的修改意见和建议对本文的完善和提高有很大帮助,在此向他们表示感谢.

参考文献
[1] Meyer B, Tapponnier P, Bourjot L, et al. Crustal thickening in Gansu-Qinghai, lithospheric mantle subduction, and oblique, strike-slip controlled growth of the Tibet plal. eau. Geophys. J. Int. , 1998, 135: 1-47. DOI:10.1046/j.1365-246X.1998.00567.x
[2] Yin A, Harrison T M. Geologic evolution of the Himalayan-Tibetan orogen. Annu. Rev. Earth Planet. Sci. , 2000, 28: 211-280. DOI:10.1146/annurev.earth.28.1.211
[3] Tapponnier P, Xu Z Q, Roger F, et al. Oblique stepwise rise and growth of the Tibet plal. eau. Science , 2001, 294: 1671-1677.
[4] 许忠淮. 东亚地区现今构造应力图的编制. 地震学报 , 2001, 23(5): 492–501. Xu Z H. A present-day tectonic stress map for eastern Asia region. AcLa Seismologica Sinica (in Chinese) , 2001, 23(5): 492-501.
[5] 徐纪人, 赵志新, 石川有三. 中国大陆地壳应力场与构造运动区域特征研究. 地球物理学报 , 2008, 51(3): 770–781. Xu J R, Zhao Z X, Ishikawa Yuzo. Regional characteristics of crustal stress field and tectonic motions in and around Chinese mainland. Chinese J. Geophys. (in Chinese) , 2008, 51(3): 770-781.
[6] McNamara D E, Owens T J, Silver P G, et al. Shear wave anisotropy beneal. h the Tibetan Plal.eau. J. Geophys. Res. , 1994, 99: 13655-13665. DOI:10.1029/93JB03406
[7] Wittlinger G, Tapponnier P, Poupinet G, et al. Tomographic evidence for localized lithospheric shear along the Altyn Tagh fault. Science , 1998, 282: 74-76. DOI:10.1126/science.282.5386.74
[8] 髙锐, 李朋武, 李秋生, 等. 青藏髙原北缘碰撞变形的深部过程-深地震探测成果之启示. 中国科学(D辑) , 2001, 44(增刊): 66–71. Gao R, Li P W, Li Q S, et al. Deep process of the collision and deformation on the northern margin of the Tibetan Plal.eau: Revelal.ion from investigal.ion of the deep seismic profiles. Sci. China Ser. D-Earth Sci. (in Chinese) , 2001, 44(增刊): 66-71.
[9] 姜枚, 许志琴, HirnA, 等. 青藏髙原及部分邻区地震各向异性和上地幔特征. 地球学报 , 2001, 22(2): 111–116. Jiang M, Xu Z Q, Hirn A, et al. Teleseismic anisotropy and corresponding features of the upper mantle in Tibet plateau and its neighboring areas. Acta GeoscienLia Sinica (in Chinese) , 2001, 22(2): 111-116.
[10] 李松林, 张先康, 张成科, 等. 玛沁-兰州-靖边地震测深剖面速度结构的初步研究. 地球物理学报 , 2002, 45(2): 210–217. Li S L, Zhang X K, Zhang C K, et al. A preliminary study on the crustal velocity structure of Maqin-Lanzhou-Jingbian by means of deep seismic sounding profile. Chinese. J. Geophys. (in Chinese) , 2002, 45(2): 210-217.
[11] Vergne J, Wittlinger G, Qian H, et al. Seismic evidence for stepwise thickening of the crust across the NE Tibetan plal. eau. Earth Planet. Sci. Lett. , 2002, 203: 25-53. DOI:10.1016/S0012-821X(02)00853-1
[12] Vergne J, Wittlinger G, Farra V, et al. Evidence for upper crustal anisotropy in the Songpan-Ganze (northeastern Tibet) terrane. Geophys. Res. Lett. , 2003, 30(11): 1552. DOI:10.1029/2002GL016847
[13] 周民都.青藏髙原东北缘地壳上地幔速度结构的地震层析成像研究[博士论文].北京:中国地震局地球物理研究所, 2005 Zhou M D. The seismic tomography on the velocity structure of the crust and upper mantle in the northeastern margin of the Qinghai-Tibet plal.eau [Ph. D. thesis]. Beijing: Institute of Geophysics, China Earthquake Administral.ion, 2005
[14] 陈九辉, 刘启元, 李顺成, 等. 青藏髙原东北缘一鄂尔多斯地块地壳上地幔S波速度结构. 地球物理学报 , 2005, 48(2): 333–342. Chen J H, Liu Q Y, Li S C, et al. Crust and upper mantle S wave velocity structure across Northeastern Tibetan plal. eau and Ordos block. Chinese J. Geophys. (in Chinese) , 2005, 48(2): 333-342.
[15] 李永华, 吴庆举, 安张辉, 等. 青藏髙原东北缘地壳S波速度结构与泊松比及其意义. 地球物理学报 , 2006, 49(5): 1359–1368. Li Y H, Wu Q J, An Z H, et al. The Poisson ral. io and crustal structure across the NE Tibetan Plal.eau determined from receiver functions. Chinese J. Geophys. (in Chinese) , 2006, 49(5): 1359-1368.
[16] Liu M J, Mooney W D, Li S L, et al. Crustal structure of the northeastern margin of the Tibetan plal. eau from the Songpan-Ganzi terrane to the Ordos basin. Tectonophysics , 2006, 420: 253-266.
[17] 段永红, 张先康, 刘志, 等. 阿尼玛卿缝合带东段地壳结构的接收函数研究. 地震学报 , 2007, 29(5): 513–522. Duan Y H, Zhang X K, Liu Z, et al. Crustal structure using receiver function in the east part of Anymaqn suture belt. AcLa Seismologica Sinica (in Chinese) , 2007, 29(5): 513-522.
[18] 张先康, 嘉世旭, 赵金仁, 等. 西秦岭-东昆仑及邻近地区地壳结构-深地震宽角反射/折射剖面结果. 地球物理学报 , 2008, 51(2): 439–450. Zhang X K, Jia S X, Zhao J R., et al. Crustal structures beneath west Qinling-east Kunlun orogen and its adjacent area-result of wide angle seismic reflection and refraction experiment. Chinese. J. Geophys (in Chinese) , 2008, 51(2): 439-450.
[19] 王椿镛, 常利军, 吕智勇, 等. 青藏髙原东部上地幔各向异性及相关的壳幔耦合型式. 中国科学(D辑) , 2007, 50(4): 1150–1160. Wang C Y, Chang L J, Lu Z Y, et al. Seismic anisotropy of upper mantle in eastern Tibetan Plateau and related crust-mantle coupling pal. tern. Sci. China Ser. DDEarih Sci. (in Chinese) , 2007, 50(4): 1150-1160.
[20] Wang C Y, Flesch L M, Silver P G, et al. Evidence for mechanically coupled lithosphere in central Asia and resulting implications. Geology , 2008, 36(5): 363-366. DOI:10.1130/G24450A.1
[21] 常利军, 王椿镛, 丁志峰, 等. 青藏髙原东北缘上地幔各向异性研究. 地球物理学报 , 2008, 51(2): 431–438. Chang L J, Wang C Y, Ding Z F, et al. Seismic anisotropy of upper mantle in the northeastern margin of the Tibet Plateau. ChfnssJ. Geophys. (in Chinese) , 2008, 51(2): 431-438.
[22] 常利军, 王椿镛, 丁志峰. 四川及邻区上地幔各向异性研究. 中国科学D辑:地球科学 , 2008, 51(12): 1683–1693. Chang L J, Wang C Y, Ding Z F. Seismic anisotropy of upper mantle in Sichuan and adjacent regions. Sci. China Sr. D-Earth Sci (in Chinese) , 2008, 51(12): 1683-1693. DOI:10.1007/s11430-008-0147-8
[23] 髙原, 滕吉文. 中国大陆地壳与上地幔地震各向异性研究. 地球物理学进展 , 2005, 20(1): 180–185. ao Y, Teng J W. Studies on seismic anisotropy in the crust and mantle on Chinese mainland. Progress in Geophysics (in Chinese) , 2005, 20(1): 180-185.
[24] Royden L H. Coupling and decoupling of crust and mantle in convergent orogens: implications for strain partitioning in the crust. J.Geophys. Res. , 1996, 101: 17679-17705. DOI:10.1029/96JB00951
[25] Holt W E. Correlated crust and mantle strain field in Tibet. Geology , 2000, 28: 67-70. DOI:10.1130/0091-7613(2000)28<67:CCAMSF>2.0.CO;2
[26] Flesch L M, Holt W E, Silver P G, et al. Constraining the extent of crust-mantle coupling in Central Asia using GPS, geologic, and shear-wave splitting data. Earth Planet. Sci. Lett. , 2005, 238: 248-268. DOI:10.1016/j.epsl.2005.06.023
[27] 郑秀芬, 欧阳飚, 张东宁, 等. 国家数字测震台网数据备份中心"技术系统建设及其对汶川大地震研究的数据支撑. 地球物理学报 , 2009, 52(5): 1412–1417. Zheng X F, Ouang B, Zhang D N, et al. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake. Chinese J. Geophys. (in Chinese) , 2009, 52(5): 1412-1417.
[28] Silver P G, Chan W W. Shear wave splitting and subcontinental mantle deformation. J. Geophys. Res. , 1991, 96: 16429-16454. DOI:10.1029/91JB00899
[29] Bowman J R, Ando M. Shear-wave splitting in the upper mantle wedge above the Tonga subduction zone. Geophys. J. R. Astron. Soc. , 1987, 88: 25-41. DOI:10.1111/j.1365-246X.1987.tb01367.x
[30] Wüstefeld A, Bokelmann G H R, Zaroli C, et al. SplitLab: a shear-wave splitting environment in Matlab. Comput. Geosci. , 2008, 34(5): 515-528. DOI:10.1016/j.cageo.2007.08.002
[31] Wüstefeld A, Bokelmann G H R. Null detection in shearwave splitting measurements. Bull. Seismol. Soc. Am. , 2007, 97(4): 1204-1211. DOI:10.1785/0120060190
[32] Barruol G, Silver P G, Vauchez A. Seismic anisotropy in the eastern US: deep structure of a complex continental plate. J. Geophys. Res. , 1997, 102: 8329-8348. DOI:10.1029/96JB03800
[33] Vecsey L, Plomerová J, Babuska V. Shear-wave splitting measurements-Problems and solutions. Tectonophysics , 2008, 462: 178-196. DOI:10.1016/j.tecto.2008.01.021
[34] Gao S S, Liu K H. Significant seismic anisotropy beneath the southern Lhasa Terrane, Tibetan Plateau. Geochem. Geophys. Geosyst. , 2009, 10: Q02008. DOI:10.1029/2008GC002227
[35] 廖武林, 丁志峰, 曾融生, 等. 喜马拉雅地区S波分裂研究. 地球物理学报 , 2007, 50(5): 1437–1447. Liao W L, Ding Z F, Zeng R S, et al. Shear wave splitting in Himalaya. Chinese J. Geophys. (in Chinese) , 2007, 50(5): 1437-1447.
[36] Liu K H, Gao S S, Gao Y, et al. Shear wave splitting and mantle flow associated with the deflected Pacific slab beneath northeast Asia. J. Geophys. Res. , 2008, 113: B01305. DOI:10.1029/2007JB005178
[37] Li Y H, Wu Q J, Jiang L J, et al. Complex seismic anisotropic structure beneath the central Tien Shan revealed by shear wave splitting. Geophys. J. Int. , 2010, 181: 1678-1686.
[38] 江丽君, 李永华, 吴庆举. 中天山及邻区S波分裂研究及其动力学意义. 地球物理学报 , 2010, 53(6): 1399–1408. Jiang L J, Li Y H, Wu Q J. The shear wave splitting of Central Tien Shan and its implications. Chinese J. Geophys. (in Chinese) , 2010, 53(6): 1399-1408.
[39] Sherrington H F, Zandt G, Frederiksen A. Crustal fabric in the Tibetan Plateau based on waveform inversions for seismic anisotropy parameters. J. Geophys. Res. , 2004, 109: B02312. DOI:10.1029/2002JB002345
[40] Herquel G, Wittlinger G, Guilbert J. Anisotropy and crustal thickness of northern-Tibet. New constraints for tectonic modeling. Geophys. Res. Lett. , 1995, 22(14): 1925-1928.
[41] Shapiro N M, Ritzwoller M H, Smith D, et al. Seismic anisotropy within the Tibetan crust. AGU Fall Meeting, San Francisco, California, USA, 2003
[42] 苏伟, 王椿镛, 黄忠贤. 青藏髙原及邻区的Rayleigh面波的方位各向异性. 中国科学D辑:地球科学 , 2008, 51(12): 1717–1725. Su W, Wang C Y, IIuang Z X. Azimuthal anisotropy of Rayleigh waves beneath the Tibetan plateau and adjacent areas. Sci. China Ser. (in Chinese) , 2008, 51(12): 1717-1725.
[43] 易桂喜, 姚华建, 朱介寿, 等. 用Rayleigh面波方位各向异性研究中国大陆岩石圈形变特征. 地球物理学报 , 2010, 53(2): 256–268. Yi G X, Yao H J, Zhu J S, et al. Lithospheric deformation of continental China from Rayleigh wave azimuthal anisotropy. ChfnssJ. Geophys. (in Chinese) , 2010, 53(2): 256-268.
[44] Alsina D, Snieder R. Small-scale sublithospheric continental mantle deformation: constraints from SKS splitting observations. Geophys. J. Int. , 1995, 123: 431-448. DOI:10.1111/gji.1995.123.issue-2
[45] Silver P G, Savage M. The interpretation of shear-wave splitting parameters in the presence of two anisotropic layers. Geophys. J. Int. , 1994, 119: 949-963. DOI:10.1111/gji.1994.119.issue-3
[46] Hartog R, Schwartz Y S. Subduction-induced strain in the upper mantle east of the Mendocino triple junction, California. J. Geophys. Res. , 2000, 105: 7909-7930. DOI:10.1029/1999JB900422
[47] Walker K T, Bokelmann G H R, Klemperer S L, et al. Shear wave splitting around the Eifel hotspot: evidence for a mantle upwelling. Geophys. J. nt. , 2005, 163: 962-980. DOI:10.1111/gji.2005.163.issue-3
[48] Fontaine F R, Barruol G, Tommasi A, et al. Uppermantle flow beneath French Polynesia from shear-wave splitting. Geophys. J. Int. , 2007, 170: 1262-1288. DOI:10.1111/gji.2007.170.issue-3
[49] Silver P G. Seismic anisotropy beneath the continents: Probing the depth of geology. Annu. Rev. Earth Planet. Sci. , 1996, 24: 385-432. DOI:10.1146/annurev.earth.24.1.385
[50] Zhang P Z, Shen Z K, Wang M, et al. Continuous deformation of the Tibetan Plateau from global positioning system data. Geology , 2004, 32(9): 809-812. DOI:10.1130/G20554.1
[51] Lave T, Avouac J P, Lacassin R, et al. Seismic anisotropy beneath Tibet: evidence for eastward extrusion of the Tibetan lithosphere. Earth Planet. Sci. Lett. , 1996, 140: 83-96. DOI:10.1016/0012-821X(96)00045-3
[52] 邓万明. 青藏髙原北部新生代板内火山岩. 北京: 地质出版社, 1998 . Deng W M. Cenozoic Intraplate Volcanic Rocks in the Northern Qinghai-Xizang Plateau (in Chinese). Beijing: Geological Pablishing House, 1998 .