2008年新疆烏恰6.8級地震序列震源特征及帕米爾東北緣應力場研究

唐蘭蘭，趙翠萍，王海濤

1 新疆維吾爾自治區地震局，烏魯木齊 830011

2 中國地震局地震預測研究所，北京 100036

2008年新疆烏恰6.8級地震序列震源特征及帕米爾東北緣應力場研究

唐蘭蘭1，趙翠萍2＊，王海濤1

1 新疆維吾爾自治區地震局，烏魯木齊 830011

2 中國地震局地震預測研究所，北京 100036

本文使用新疆區域數字地震臺站記錄的寬頻帶長周期數字波形資料，在時間域反演了2008年10月5日新疆烏恰6.8級地震的強余震及其周圍先后發生的52次中等強度地震的矩張量解，結合Harvard大學在該區域的地震矩張量結果，研究了帕米爾東北緣的應力場分區特征.研究結果顯示，位于印度板塊向歐亞板塊推擠的前緣及向北凸出的弧型構造的最北緣的卡茲克阿爾特弧形活動褶皺－逆斷裂帶，以逆沖推覆活動為主，并有部分走滑類型的地震，基本不存在正斷層類型的地震；該弧型構造近東西走向的頂部（文中的西區）與其北西走向的東側（文中的東區）的局部應力場最大主壓應力方向不同，分別為NW、NNE方向，顯示出在承受印度板塊向歐亞板塊俯沖作用的同時，東區也更多的受到了塔里木塊體順時針旋轉作用的影響.位于帕米爾陸內俯沖和變形作用強烈、碰撞造成深源地震帶東段的南區，地震以走滑錯動為主，逆斷、正斷層都有，顯示出相對復雜的應力狀態.位于帕米爾高原內部的西區和南區的應力場最大主壓應力方向一致，由北向南，由最大主壓應力軸接近水平，過渡為最大主張應力軸接近水平，一定程度揭示了板塊俯沖的狀態.結合南區和西區的地震深度差異及機制解中斷層面的傾角，推測在中帕米爾的東部，由北向南的板塊俯沖至150～170km深度，俯沖角度為60°左右.

帕米爾地區，地震矩張量反演，應力場特征

1 引 言

帕米爾東北側是天山褶皺帶、西昆侖造山帶和塔里木塊體3個地質構造單元的交接帶，在印度板塊和歐亞板塊的碰撞和持續匯聚作用下，作為中國大陸受印度板塊動力作用最強烈的地區之一，帕米爾的陸內俯沖和變形作用非常強烈［1－2］.地質學及地球物理觀測數據表明歐亞大陸巖石圈深部向帕米爾俯沖至少有200km，使帕米爾的地殼比正常的地殼增厚了大約2倍［3－4］.帕米爾構造帶則沿其前緣的帕米爾北緣推覆構造帶向北仰沖到歐亞大陸的穩定地塊之上，南天山也因受到來自印度板塊和歐亞板塊碰撞的影響，向南仰沖到穩定的塔里木板塊北緣之上，該地區現今的縮短速率達到19～20mm／a［5－6］，幾乎占去印度板塊向北推移速率的一半［1］.印度板塊由南向北強烈的推擠，在位于印度板塊向歐亞板塊推擠前緣的帕米爾東北緣，形成了一系列壯觀的向北凸出的弧型推覆構造，這些弧型推覆構造在東西兩側由正向逆沖漸變為斜沖－走滑大斷裂，顯示出東側右旋、西側左旋的斷層活動特征.弧型推覆構造的東側、NW走向的逆斷層與斜切天山NW向的塔拉斯－費爾干那深大右旋走滑逆沖斷裂帶、NW向的西昆侖北緣走滑斷裂構造在新疆烏恰地區交匯，這里發生的1985年烏恰7.3級強震破裂機制復雜，顯示出上述幾條斷裂的影響［7］.

上述向北凸出的帕米爾北緣弧形推覆構造帶，由艾卡爾特弧形活動褶皺－逆斷裂帶（簡稱艾卡爾特弧）與卡茲克阿爾特弧形活動褶皺－逆斷裂帶（簡稱卡茲克阿爾特弧）南、北兩條巨型邊緣弧形構造帶及其間的推覆構造構成.頻繁發生的中強地震沿這兩個活動褶皺－逆斷裂帶形成了南北兩條極為顯著的深、淺源地震帶（圖1）.其中卡茲克阿爾特弧形活動褶皺－逆斷裂帶自西向南東由近東西走向的卡巴加特弧和NW走向的烏帕爾兩個次級弧型構造構成，是6級以上地震頻繁發生的場所.2008年10月5日新疆烏恰MS6.8地震（圖1中5號地震）位于卡茲克阿爾特斷裂的中西段，在1974年8月11日烏恰西南MS7.4地震（圖1中1號地震）以東20km左右，其主要余震沿該斷裂呈北北東向分布.而1985年烏恰MS7.3強震（圖1中3號地震）則發生在斷裂東端的NW走向的烏帕爾弧，這里還發生了1983年2月13日MS6.7地震（圖1中2號地震）、1990年4月17日疏附MS6.4地震（圖1中4號地震）.

獨特的構造環境使得帕米爾東北緣成為地球動力學研究的典型地區［2，7，8－16］.然而，前人對研究區域或其附近開展的應力場研究，都是基于Harvard大學的矩張量解［17］或者利用初動解［18］，前者只有1970年以來4.8級以上地震，后者的質量受觀測條件的局限.本文利用“十五”以來新建的新疆區域臺網中寬頻帶臺站記錄的三分量地震波形數據，在時間域反演了烏恰MS6.8地震震源區及其附近區域2006年以來的52次MS3.7以上地震的矩張量解.在結合Harvard大學的矩張量解基礎上，研究了帕米爾東北緣這個典型構造區域近年來的震源機制類型及局部應力場特征.

2 研究方法

獲取地震的震源機制解是認識地震發震斷層及破裂特征的主要手段.地震的震源應力場攜帶著區域構造應力場及構造運動的信息，因此震源機制解還為構造應力場的研究提供了基礎數據.由于中等地震的能量相對較小，通常難以利用遠場波形反演地震矩張量解.利用P波初動求震源機制解，當近場沒有足夠的臺站分布時，觀測所能提供的初動資料數量以及這些初動點在震源球上分布的均勻程度對結果有相當大的影響，因此獲得的機制解常常可靠性較差.隨著新建的“十五”國家和區域數字地震臺網的運行，利用區域震源寬頻帶記錄進行中小地震的震源參數的研究，成為近年來一個活躍的研究領域［19－21］.許多研究表明，由于區域長周期地震波對速度結構的橫向變化及密度的非均勻性相對不敏感，如果震源位置比較準確，地震滿足震源的同步假設，波形信噪比高，則用區域范圍稀疏臺站記錄的三分量長周期波形就足以反演得到穩定的矩張量解.另外，由P波初動得到的機制解僅僅是初始破裂面，由波形反演結果得到的則是整個破裂過程的信息，由于初始破裂方向并不一定與地震的斷層錯動完全一致，因此矩張量解得到的震源機制結果相對更加完整和可靠［22－26］.

本文使用了Dreger等［22］提出的利用區域范圍長周期體波三分量波形在時間域反演地震矩張量解的程序，方法原理詳見文獻［26］.該方法目前在南加州臺網及全球多個國家臺網中運行以提供近實時地震矩張量解，在我國伽師及汶川地震序列研究中也得到了應用［26－27］.理論地震圖的正演計算使用經Saikia［28］改進的離散波數積分方法.當使用長周期波形資料時，地震滿足震源為點源的假設，震源時間函數直接取δ函數.每個參加反演的臺站可以分別使用不同的地殼速度、密度等.由于部分寬頻帶臺站儀器的自身低頻噪音較為顯著，長周期濾波后觀測波形受到影響，本文使用三分量速度波形進行反演.

具體反演過程是：首先對觀測波形進行去均值、去傾校正，反褶積儀器傳遞函數，將記錄變為實際速度值（m／s），再旋轉到Z、R、T分量.然后對觀測資料使用Butterworth帶通濾波器濾波到需要的長周期頻段，以剔除較低頻和較高頻噪聲而不破壞用于矩張量反演的信息.對計算的Green函數也濾波到與觀測數據同樣的頻率范圍.最后采用迭代擬合相關系數的方法對觀測波形和理論波形進行擬合.

對于烏恰MS6.8主震，由于該方法對于6級地震要求區域范圍的臺站在0.01～0.05Hz幅頻響應平坦，能滿足條件的僅有喀什臺，所以本研究采用傳統的P波初動方法求解該地震的震源機制解.

本文進一步采用了Gephart和Forsyth［29－30］提出的一種被廣泛使用的確定應力場的反演程序FMSI，對研究區域的應力場主軸方向和傾角進行了反演.該方法假設斷層面上的剪切應力方向與斷層的滑動方向一致，然后利用網格搜索方法尋找一組地震的最佳擬合應力張量.每個地震的殘差定義為使兩個節面上的滑動角與應力模型預測的滑動角一致的最小旋轉角.平均殘差表示研究區域內的應力張量的非均勻程度，平均殘差小于6°代表研究區域較均勻的應力條件，而大于9°代表應力場有較高的非均勻性［31－32］.

3 數據及處理

2003年以來，“十五”新疆區域數字地震臺網在帕米爾東北緣附近的寬頻帶數字地震臺站有喀什（KSH）、烏恰（WUQ）、阿圖什（ATS）、八盤水磨（BPM）、英吉沙（YJS）、葉城（YCH）、塔什庫爾干（TAG）.本文利用上述7個臺站記錄的寬頻帶三分量地震波形數據的長周期波形，反演2006－2010年期間，除2008年10月5日MS6.8地震外的MS3.7以上地震的矩張量解.

近年來，隨著我國區域臺網儀器布設數量劇增，數字地震觀測數據大量產出，應用這些數字地震資料的研究得到了快速發展，也對地震數據的高質量產出提出了要求.Niu等［33］通過對寬頻帶儀器遠場直達P波三分量的極性分析，檢測了全國寬頻帶地震臺在運行期間可能出現的儀器方位角以及初至波極性倒轉等問題，該研究結果顯示，本研究使用的7個臺站未出現極性反轉現象，其中的阿圖什、喀什、塔什庫爾干3個臺地震計“BHN”方向與地理北極方位角偏差分別為10°、11°、－10°，略大于認為可以接受的8°偏差.因此，在進行波形反演前，我們進行了方位角偏差校正實驗，發現經過校正后的波形得到的矩張量反演結果與校正前得到的結果差異很小，可以忽略其方位角偏差造成的影響.

反演地震矩張量的重要步驟之一是計算格林函數，其計算結果在很大程度上影響反演結果的可靠性，因此速度模型的建立是一件非常基礎而關鍵的工作.本文計算理論地震圖使用了Crust2.0的2°×2°分層的速度和密度結構模型，Crust2.0模型是一個被廣泛使用的模型，分辨率2°×2°.作者將該模型與研究區已有的速度結構進行了對比，差別不大.同時由于反演方法使用區域范圍長周期的波形數據，反演結果對結構的依賴性不強.研究中根據各個臺站的位置，最終使用了3種速度模型：烏恰（WUQ）、喀什（KSH）、阿圖什（ATS）使用一種速度模型，英吉沙（YJS）、八盤水磨（BPM）使用一種速度模型，葉城（YCH）、塔什庫爾干（TAG）使用一種速度模型.

本研究中，受區域臺網儀器觀測頻帶的限制，未能反演2008年10月5日MS6.8地震的矩張量解.但我們自IRIS臺網的臺站及新疆區域臺網中，仔細讀取了具有清晰初動的133個臺站的P波初動，彌補了新疆臺網均分布于震源東側的缺陷，采用初動方法確定了其機制解.臺站分布情況如圖2所示.

圖1 研究區1964—2008年期間MS＞5中強地震震中分布圖（據ISC地震目錄）不同顏色的圓圈代表不同的地震深度.Fig.1 Distribution of MS＞5.0earthquake in the studied region from 1964to 2008（ISC bulletin）Different color circles marks earthquake with different depth.

圖2 計算烏恰6.8級地震震源機制解所用的133個地震臺站分布圖Fig.2 Distribution of 133stations for calculating the focal mechanism of Wuqia MS6.8earthquake

4 烏恰MS6.8地震及其余震震源機制解

4.1 烏恰MS6.8地震震源機制解

圖3a為利用133個初動數據得到的烏恰6.8級地震的震源機制解.Harvard大學在震后即給出了2008年10月5日MS6.8地震的矩心矩張量解（圖3b）.由圖3可見，哈佛的解存在一定的走滑分量，兩個結果斷層走向和震源應力場的P、T軸一致，均表明此次地震是逆斷層性質的地震.在南傾的節面附近出現了矛盾符號的臺站，其中的一個節面與余震分布及發震位置所處的斷層幾何參數一致，表明是卡茲克阿爾特斷裂帶中西段發生破裂的結果.根據6.8級地震后的余震分布情況，判斷走向75°、傾角70°、滑動角90°的節面為此次地震的發震斷層面.

圖3 2008年10月5日烏恰MS6.8級地震雙力偶震源機制解（a）本文利用133個P波初動解；（b）Harvard矩張量解；圖中陰影區為壓縮區，小實心圈表示P波初動向上，小空心圈表示P波初動向下，大實心圈為P軸，大空心圈為T軸.Fig.3 The double－couple focal mechanism for Oct，5，2008 MS6.8event（a）The obtained solution using 133P－wave first motions；（b）The Harvard CMT solution.Compressional quadrants are shaded，the smaller solid circle marks positive P－wave first motions，whereas the smaller open circle marks negative ones，the larger solid circle represents Paxis，the larger open circle represents Taxis.

表1 2008年10月5日烏恰6.8級地震雙力偶震源機制解Table 1 The double－couple focal mechanism for Oct，5，2008 MS6.8event

4.2 烏恰MS6.8地震余震序列及其周緣地震的矩張量解

研究中得到2006年以來烏恰6.8級地震余震序列及其周緣共52次地震的矩張量解，在52次地震中，多數解是由2個以上臺站記錄的三分量速度波形數據反演得到，理論地震圖與觀測地震圖的相關系數大于0.7的有46個，占86.8%；附表1給出了最佳雙力偶分解得到的雙力偶機制解，圖4給出了三次地震的矩張量解及使用臺站的三分量波形擬合情況.全球很多研究機構對地震矩張量進行反演研究，其中哈佛大學在此方面被公認為最權威的研究機構，他們對全球絕大部分Mw＞4.8地震矩張量作了反演.為了進一步分析本文矩張量反演結果的可靠性，本文選取哈佛大學在該地區和時間段內也做出結果的8個地震進行了對比（圖5），圖5中僅繪制了矩張量反演得到的最佳雙力偶分解得到的雙力偶機制解的兩個截面和P、T軸，結果顯示本文的計算結果和哈佛大學的結果較為一致，尤其是斷層走向和震源應力場的P、T軸.

本研究計算的53次地震中，有31次地震為烏恰6.8級地震的余震，圖6中給出了包括烏恰6.8級地震在內的53次地震震中分布及雙力偶震源機制解的下半球投影，由圖6可見，這些余震基本沿卡茲克阿爾特斷裂的中西段向西展布.余震序列震源機制解結果顯示出以與主震一致的逆沖型為主，伴有走滑破裂類型.

圖4 理論波形與實際波形擬合及反演結果（實線為記錄波形，虛線為理論波形）由上至下分別為：（a）2008年10月13日MS5.2；（b）2008年10月14日MS5.4；（c）2009年7月28日MS4.3.圖中VR代表理論與觀測波形的擬合方差減小.Fig.4 Examples of inversion result，with comparison of the synthetic data with the observed data（Solid line represents observed data，dashed line represents the synthetic waveform）（a）Obc.13，2008，MS5.2；（b）Obc.14，2008，MS5.4；（c）Jul.28，2009，MS4.3.VR is the variance reduction of the observed waveforms with the synthetic ones.

5 帕米爾弧東北緣震源機制及應力場特征

由圖6可見，本文研究的地震在空間上分為三叢，北面的2叢沿卡茲克阿爾特弧分布，分別位于卡茲克阿爾特斷裂的中西段（卡巴加特逆斷裂帶）及東段前緣的烏帕爾弧段，前者是此次6.8級地震余震區，后者是1985年烏恰7.3級地震的震源區；南邊的一叢位于中帕米爾中深源地震帶的東段.我們將Harvard大學給出的1976—2008年MS≥5.0以上49個地震的矩張量震源機制解一起進行統計分析，發現2008年10月5日烏恰6.8級地震震源附近的地震大都是逆沖和走滑性質，45個機制解中，逆斷層占60%，走滑占36%；1985年烏恰7.3級地震震源區的地震逆斷類型比走滑類型稍多，30個地震中逆斷層占57%，走滑占40%；而南部的地震以走滑為主，22個深震中走滑地震占54.5%，逆斷、正斷層各占約23%.由附表1可見，西區斷層節面的傾角優勢分布于30°～60°之間，而東區的則更加直立，優勢分布于60°～90°之間，顯示出沿卡茲克阿爾特弧形斷裂中西部到東端部，斷面傾角逐漸增大且趨于直立的特征，與地質考察結果基本一致［7］.

圖5 本文矩張量反演結果（b）與哈佛CMT結果（a）對比圖，其中紅色實心圈代表P軸，白色實心圈代表T軸Fig.5 The moment tensor solution obtained by the author（b）and the Harvard CMT solution（a）.Red solid circle represent Paxis，white solid circle represent Taxis.

圖6 53次地震的震源機制解及臺站位置藍色三角形為臺站位置，機制解上方的數字代表發震年月日，參見附表1.Fig.6 Distribution of the obtained 53focal mechanisms and the location of the stations（blue triangle）The number above the focal mechanisms means time when the event happened，referring to the index in the attached Table 1.

我們把每個地震的P、T、N軸用線段繪于圖7，線段的方向代表方位角，線段的長短表示傾角，越短表示越直立，越長表示越水平.顯然，研究區三叢地震的P、T、N軸分布也顯示出顯著的空間分區特征.西區地震的P軸走向NNW，傾角接近水平，N軸較緩，T軸方向NNE，多數地震的傾角較陡；東區的地震以走滑為主，P軸走向NNE，傾角較緩，N軸傾角較陡，走向NNW，T軸走向近EW，傾角較緩；南區地震的P軸走向NNW，傾角較緩，T軸走向NNE，傾角接近水平.

根據上述分區特征，我們利用FMSI方法計算了各區的應力場主軸方向（表2、圖8）.從三個區的應力場反演結果（表2）可以看出，三個研究區的應力場存在差異，西區與南區的最大主壓應力方向基本一致（N30°W左右），但南區的最大主壓應力軸存在31°的傾角，而西區最大主壓應力軸基本水平，東區的最大主壓應力方向為N25°E，最大主壓應力軸與最大主張應力軸都具有22°的傾角.可見西區水平擠壓作用最為明顯，東區水平擠壓和拉張作用基本相當，而南區的水平拉張作用強于擠壓作用.

圖7 P、T、N軸空間分布圖中紅色為本次研究得到的結果，藍色為Harvard大學CMT解的結果，線段長短代表傾角的大小（90°－傾角），越短表示越直立.Fig.7 Distribution of the P，T，NaxisRed lines are resolutions from this study，blue lines are 42 resolutions from Harvard CMT.The line length is proportional to plunge angle（90°－plunge angle），the shorter means more vertical.

圖8 震源區西區（a）、東區（b）和南區（c）的應力場反演結果圖中較大的圓形、三角形、方形分別代表最大、中等及最小壓應力軸，其中較小的圓形、三角形、方形分別代表P、N、T軸投影位置.Fig.8 The result of stress field inversion（a）、（b）and（c）give the stress inversion result for the western，eastern and southern region respectively.The larger circle，triangle and rectangle representσ1，σ2，andσ3stress axis respectively；smaller circle，triangle and rectangle represent the projection location of the P，Nand Taxis respectively.

表2 應力場反演結果Table 2 The inversion results of the stress field

6 討論與結論

本文研究的西區位于卡茲克阿爾特弧形活動褶皺－逆斷裂帶的中東部，是印度板塊向歐亞板塊推擠的前緣及向北凸出的弧型構造的中間部位，這里斷層走向近東西向；而東區位于卡茲克阿爾特弧形構造向東的轉彎處，這里斷層方向由其西部的近東西向轉為NW向.東、西區地震的震源深度大多在30 km以內.根據矩張量反演結果，東、西區均有近60%的地震為逆斷層，近40%的地震為走滑性質，基本沒有正斷層地震.西區最大主壓應力軸N36°W，基本水平；東區N25°E，傾角22°.可見在帕米爾北緣卡茲克阿爾特弧形活動褶皺－逆斷裂帶的中、東部，以逆沖推覆活動為主，并有部分走滑類型的地震，基本不存在正斷層類型的地震；在帕米爾東北緣的卡茲克阿爾特弧型構造東、西兩側，局部應力場最大主壓應力方向分別為NNE、NW方向.研究區的東區處于西昆侖造山帶、天山褶皺帶以及塔里木塊體的交匯部位，構造位置特殊，是帕米爾東北緣弧形斷裂帶與塔拉斯－費爾干那深大斷裂、西昆侖斷裂帶的交界處，我們推斷，該區域在承受印度板塊向歐亞板塊的俯沖作用的同時，也更多地受到了塔里木塊體的順時針旋轉作用的影響.

樓小挺等［4］對整個興都庫什－帕米爾中深源地震帶開展了研究，本文的南區位于該研究中定義為P帶PNE段的東側，地震深度達150km左右，這里的地震分布為NE向，且再往東，出現了弧形轉向.根據本文結果，地震以走滑錯動為主，20個深震中走滑地震占一半，逆斷、正斷層各占約1／4.應力場反演結果顯示該區的最大主壓應力方向為N28°W，最大主張應力軸方向NE，基本水平.與樓小挺等［4］的研究結果中，PNE段地震的T軸近水平并沿地震帶走向的結論一致.

樓小挺等［4］的研究認為，帕米爾俯沖帶的中段俯沖深度淺，東段沒有類似于興都庫什俯沖帶以及帕米爾西南段的斷裂或錯動.根據本文結果，南區的最大主壓應力方向與其北側的卡茲克阿爾特弧形推覆構造頂部（西區）一致，由北向南，由最大主壓應力軸接近水平，過渡為最大主張應力軸接近水平，在一定程度揭示了板塊俯沖的動力學過程.張家聲等［15］認為帕米爾中源地震群與帕米爾深俯沖巖板向下由緩傾變陡的深度大體相當，結合南區和西區的地震深度差異及機制解中斷層面的傾角，推測在中帕米爾的東部，由北向南的板塊俯沖至150～170km深度，俯沖角度為60°左右.

本文得到的應力場方向與中亞地區重復GPS測量結果得到的位移方向一致［14，34－35］.研究結果顯示出，在大尺度的構造運動及動力作用下，局部應力場受所處局部構造影響，尤其是在本文所研究的典型的特殊構造部位.

致 謝 研究過程中使用了美國伯克利地震實驗室的矩張量反演軟件、IRIS的波形數據，大部分圖件使用GMT繪制，特此申明.論文形成過程中，田勤儉研究員、沈軍研究員、高國英研究員在地質構造方面提供了幫助，在此一并致謝！

（References）

［1］ Molnar P.Continental tectonics in the aftermath of plate tectonics.Nature，1988，335（6186）：131－137.

［2］ Burtman V S，Molnar P.Geological and geophysical evidence for deep subduction of continental crust beneath the Pamir.Geol.Soc.Am.Spec.Pap.，1993，281：1－76.

［3］ 汪素云，時振梁，環文林.1980年2月14日新疆葉城地震—一次破壞性的中源地震.地震學報，1992，14（2）：137－143.Wang S Y，Shi Z L，Huan W L.A destructive intermediatedepth earthquake of Feb 14 1980in Yecheng of Xinjiang.Acta Seismologica Sinica（in Chinese），1992，14（2）：137－143.

［4］ 樓小挺，刁桂苓，葉國揚等.帕米爾—興都庫什地區中源地震的空間分布和震源機制解特征.地球物理學報，2007，50（5）：1448－1455.Lou X T，Diao G L，Ye G Y，et al.Spatial distribution of intermediate earthquakes and focal mechanism solutions in the Pamir－Hindu Kush region.Chinese J.Geophys.（in Chinese），2007，50（5）：1448－1455.

［5］ Abdrakhmatov K Y，Aldazhanov S A，Hager B H，et al.Relatively recent construction of the Tien Shan inferred from GPS measurements of present－day crustal deformation rates.Nature，1996，384（6608）：450－453.

［6］ 王琪，丁國瑜，喬學軍等.用GPS研究南天山（伽師）地區現今地殼變形.地震學報，2000，22（3）：263－270.Wang Q，Ding G Y，Qiao X J，et al.Research on present crustal deformation in the southern Tianshan（Jiashi），Chinaby GPS geodesy.Acta Seismologica Sinica（in Chinese），2000，22（3）：263－270.

續附表1

［7］ 陳杰，曲國勝，胡軍等.帕米爾北緣弧形推覆構造帶東段的基本特征與現代地震活動.地震地質，1997，19（4）：301－312.Chen J，Qu G S，Hu J，et al.Arcuate thrust tectonics and its contemporary seismicity in the eastern section of the external zone of the Pamir.Seismology and Geology（in Chinese），1997，19（4）：301－311.

［8］ 趙俊猛，劉閏棟，趙國澤等.天山造山帶與準噶爾盆地殼幔過渡帶及其動力學含義.中國科學，2001，31（4）：272－282.Zhao J M，Liu G D，Zhao G Z，et a1.Crust－mantle transitional zone of Tianshan orogenic belt and Junggar Basin and its geodynamic implication.Science in China（Series D）（in Chinese），2001，31（4）：272－282.

［9］ 趙俊猛，李純植，馬宗晉.天山分段性的地球物理學分析.地學前緣，2003，10（特刊）：125－131.Zhao J M，Li C Z，Ma Z J.Geophysical evidence for segmentation of the Tianshan.Earth Science Frontiers（in Chinese），2003，10（Special issue）：125－131.

［10］ 胥頤，劉建華，劉福田等.天山—帕米爾結合帶的地殼速度結構及地震活動研究.地球物理學報，2006，49（2）：469－472.Xu Y，Liu J H，Liu F T，et al.Crustal velocity structure and seismic activity in the Tianshan－Pamir conjunctive zone.Chinese J.Geophys.（in Chinese），2006，49（2）：469－472.

［11］ 沈軍，汪一鵬，趙瑞斌等.帕米爾東北緣及塔里木盆地西北部弧形構造的擴展特征.地震地質，2001，23（3）：381－389.Shen J，Wang Y P，Zhao R B，el a1.Propagation of Cenozoic arcuate structures in northeast Pamir and northwest Tarim basin.Seismology and Geology（in Chinese），2001，23（3）：381－389.

［12］ 張先康，張成科，嘉世旭等.天山地殼結構及其動力學意義.∥周光召.西部大開發科教先行與可持續發展.北京：中國科學技術出版社，2000：556.Zhang X K，Zhang C K，Jia S X，et a1.Crust structure of Tianshan Mountains and its geodynamic implication.∥Zhou G Z.On Top Priority of Science and Education in the General Development of Western China and Sustainable Development of the Project（in Chinese）.Beijing：Science and Technology Press，2000：556.

［13］ 張先康，趙金仁，張成科等.帕米爾東北側地殼結構研究.地球物理學報，2002，45（5）：665－671.Zhang X K，Zhao J R，Zhang C K，et al.Crust structure at the northeast side of the Pamirs.Chinese J.Geophys.（in Chinese），2002，45（5）：665－671.

［14］ 劉志，張先康，周雪松等.帕米爾東北側地殼物性結構及其發震環境探討.地震學報，2003，25（3）：242－249.Liu Z，Zhang X K，Zhou X S，et al.Study of crustal structure in terms of physical property and discussion on seismogenic environment in the northeastern Pamirs.Acta Seismologica Sinica（in Chinese），2003，25（3）：242－249.

［15］ 張家聲，單新建，李建華等.帕米爾地區現今大陸深俯沖：地震構造和動力學解釋.巖石學報，2005，21（4）：1215－1227.Zhang J S，Shan X J，Li J H.Recent deep subducting of continental crust in Pamier：An interpretation on seismotectonics and geodynamics.Acta Petrologica Sinica（in Chinese），2005，21（4）：1215－1227.

［16］ 刁桂苓，王海濤，高國英等.伽師強震系列應力場的轉向過程.地球物理學報，2005，48（5）：1062－1068.Diao G L，Wang H T，Gao G Y，et al.A deflection process for the stress field of the Jiashi strong earthquake sequence.Chinese J.Geophys.（in Chinese），2005，48（5）：1062－1068.

［17］ 孫文斌，和躍時，常征等.帕米爾—興都庫什地區板塊俯沖及其應力狀態.地震地質，2009，31（2）：207－217.Sun W B，He Y S，Chang Z，et al.The plate subduction and stress state in the Pamir－Hindu Kush region.Seismology and Geology（in Chinese），2009，31（2）：207－217.

［18］ 高國英，溫和平，聶曉紅.1991—2002年新疆中強震震源機制解分析.地震，2005，25（1）：81－87.Gao G Y，Wen H P，Nie X H.Analysis of focal－mechanism solution of moderately strong earthquakes in Xinjiang from 1991to 2002.Earthquake（in Chinese），2005，25（1）：81－87.

［19］ 周榮茂，陳運泰，吳忠良.由矩張量反演得到的北部灣地震的震源機制.地震學報，1999，21（6）：561－569.Zhou R M，Chen Y T，Wu Z L.Moment tensor inversion for focal mechanism of the Beibuwan earthquakes.Acta Seismologica Sinica（in Chinese），1999，21（6）：561－569.

［20］ 倪江川，陳運泰，王鳴等.云南祿勸地震部分余震的矩張量反演.地震學報，1991，13（4）：4l2－4l9.Ni J C，Chen Y T，Wang M，et al.Moment tensor inversion of some aftershocks of the April 18，1985，Luquan，Yunnan，China，Earthquake.Acta Seismologica Sinica（in Chinese），1991，13（4）：412－419.

［21］ 吳忠良，牟其鐸.關于中國數字地震學發展戰略的討論.國際地震動態，1994，181（1）：5－8.Wu Z L，Mou Q D.Discussion on the strategy of developments of digital seismology in china.Recent Developments in World Seismology（in Chinese），1994，181（1）：5－8.

［22］ Dreger D S，Helmberger D V.Determination of source parameters at regional distances with three－components sparse network data.J.Geophys.Res.，1993，98（B5）：8107－8125.

［23］ Romanowicz B，Dreger D，Pasyanos M，et al.Monitoring of strain release in central and northern California using broadband data.Geophys.Res.Let.，1993，20（15）：1643－1646.

［24］ Randall G E，Ammom C J，Owens T J.Moment tensor estimation using regional seismograms from a Tibetan Plateau protable network deployment.J.Geophys.Res.，1995，22（13）：1665－1668.

［25］ Tan Y，Zhu L P，Helmberger D V，et al.Locating and modeling regional earthquakes with two stations.J.Geophys.Res.，2006，111：B01306，doi：10.1029／2005JB003775.

［26］ 趙翠萍，陳章立，鄭斯華等.伽師震源區中等強度地震矩張量反演及其應力場特征.地球物理學報，2008，51（3）：782－792.Zhao C P，Chen Z L，Zheng S H，et al.Moment inversion of moderate earthquakes and the locally perturbed stress field in the Jiashi source region.Chinese J.Geophys.（in Chinese），2008，51（3）：782－792.

［27］ 王勤彩，陳章立，鄭斯華.汶川大地震余震序列震源機制的空間分段特征.科學通報，2009，54（16）：2348－2354.Wang Q C，Chen Z L，Zheng S H.Spatial segmentation characteristic of focal mechanism of aftershock sequence of Wenchuan Earthquake.Chinese Science Bulletin（in Chinese），2009，54（13）：2263－2270，doi：10.1007／s11434－009－0367－0.

［28］ Saikia C K.Modified frequency－wavenumber algorithm for regional seismograms using Filon′s quadrature：modelling of Lg waves in eastern North America.Geophys.J.Int.，1994，118（1）：142－158.

［29］ Gephart J W，Forsyth D W.An improved method for determining the regional stress tensor using earthquake focal mechanism data：Application to the San Fernando earthquake sequence.J.Geophys.Res.，1984，89（B11）：9305－9320，doi：10.1029／JB089iB11p09305.

［30］ Gephart J W.Stress and the direction of slip on fault planes.Tectonics，1990，9（4）：845－858，doi：10.1029／TC009i004p00845.

［31］ Wyss M，Lu Z.Plate boundary segmentation by stress directions：Southern San Andreas fault，California.Geophys.Res.Lett.，1995，22（5）：547－550，doi：10.1029／95GL00074.

［32］ Wyss M，Liang B Y，Tanigawa W R，et al.Comparison of orientations of stress and strain tensors based on fault plane solutions in Kaoiki，Hawaii.J.Geophys.Res.，1992，97（B4）：4769－4790，doi：10.1029／91JB02968.

［33］ Niu F L，Li J.Component azimuths of the CEArray stations estimated from P－wave particle motion.Earthquake Science，2011，24（1）：3－13.

［34］ 王琪，丁國瑜，喬學軍等.天山現今地殼快速縮短與南北地塊的相對運動.科學通報，2000，45（14）：1543－1547.Wang Q，Ding G Y，Qiao X J，et al.The contemporary crustal shortening of Tienshan and its movement relative to the South－North block.Chinese Science Bulletin（in Chinese），2000，45（14）：1543－1547.

［35］ 張培震，王琪，馬宗晉.中國大陸現今構造運動的GPS速度場與活動地塊.地學前緣，2002，9（2）：430－441.Zhang P Z，Wang Q，Ma Z J.GPS velocity field and active crustal blocks of contemporary tectonic deformation in continental China.Earth Science Frontiers（in Chinese），2002，9（2）：430－441.

Study on the source characteristics of the 2008 MS6.8Wuqia，Xinjiang earthquake sequence and the stress field on the northeastern boundary of Pamir

TANG Lan－Lan1，ZHAO Cui－Ping2＊，WANG Hai－Tao1
1 Seismological Bureau of Xinjiang Uygur Autonomous Region，Urumqi 830011，China
2 Institute of Earthquake Science，China Earthquake Administration，Beijing100036，China

Moment tensor solutions of 52earthquakes occurred in the 2008 MS6.8Wuqia source region，Xinjiang，were estimated in time domain by using the long period waveforms recorded by the regional broadband stations.Then combined with the Harvard solutions，the characteristic of stress field in northeastern margin of Pamir are studied.The result reveals that overthrusting isthe main pattern on the curved Kazkeaerte active reverse fault fold zone，which is located in front of the collision zone between Indian plate and Eurasian plate，as well as the northern edge of the northward bulge arc structure.And there are also a few strike－slip faults，normal fault almost does not exist.The trend of the local maximum principal stress is different on the top of the arc（the western region in the paper）and the eastern boundary（the eastern region in the paper），which is NNE and NNW respectively，suggesting the influence of clockwise rotation of Tarim block，while the whole Pamir undergoing the force of the India plate crushing to the Eurasian plate.The types of focal mechanism solutions in the southern region，which is located on the eastern section of Pamir deep seismic zone formed by intensive collision，subduction，and deformation process，are mainly strike－slip，and with some normal and reverse faults as well.The trending of the principal compressive stress in western and southern Pamir is the same，but the dip angle in western region is horizontal and the principal tensional stress in the southern region is horizontal，which may reveal that the subduction depth of the Eurasian continental lithosphere in the eastern region of central Pamir is about 150to 170km deep and the subducting angle is about 60°.

Pamir region，Seismic moment tensor inversion，Stress field features

10.6038／j.issn.0001－5733.2012.04.018

P315

2011－05－05，2012－03－13收修定稿

大震應急流動觀測與科學產品加工關鍵技術研究與應用（201108014）、中亞大陸強震構造格局及其動力環境（2008CB425703）、亞歐大陸中、新生代構造變動及其與地震災害和油氣資源關系的研究（2008DFA20700）、2010年度震情跟蹤合同制定向工作任務（2010020107）聯合資助.

唐蘭蘭，女，1981年生，碩士研究生，主要從事數字地震學研究.E－mail：tllfeg1981＠163.com

＊通訊作者趙翠萍，女，研究員，主要從事數字地震學應用研究.E－mail：zhaocp＠seis.ac.cn

唐蘭蘭，趙翠萍，王海濤.2008年新疆烏恰6.8級地震序列震源特征及帕米爾東北緣應力場研究.地球物理學報，2012，55（4）：1228－1239，

10.6038／j.issn.0001－5733.2012.04.018.

Tang L L，Zhao C P，Wang H T.Study on the source characteristics of the 2008 MS6.8Wuqia，Xinjiang earthquake sequence and the stress field on the northeastern boundary of Pamir.Chinese J.Geophys.（in Chinese），2012，55（4）：1228－1239，doi：10.6038／j.issn.0001－5733.2012.04.018.

（本文編輯 胡素芳）