青藏高原東緣地應(yīng)變演化特征

青藏高原東緣地應(yīng)變演化特征

姜永濤1，2， 張永志1， 吳然1， 王帥1

(1.長安大學(xué)，陜西 西安710000； 2.南陽師范學(xué)院，河南 南陽473061)

摘要：利用青藏高原東緣1999—2013年間多期GPS水平速率觀測數(shù)據(jù)，基于多面函數(shù)擬合，計算球面坐標(biāo)系下區(qū)域不同時期的面應(yīng)變和最大剪應(yīng)變，分析地應(yīng)變的時空演化特征，結(jié)合不同時期發(fā)生的中強(qiáng)以上地震(MS>6.0)，研究期間大震分布與地應(yīng)變時空演化特征的關(guān)系，主要結(jié)論如下：(1)青藏高原東緣面應(yīng)變分布與地塊有一定的對應(yīng)關(guān)系，面應(yīng)變的差異會在塊體邊界和內(nèi)部形成不同的斷層閉鎖形式，與地震發(fā)生位置和震源機(jī)制有一定的關(guān)聯(lián)；(2)區(qū)域最大剪應(yīng)變的高值區(qū)對應(yīng)于構(gòu)造活動性較強(qiáng)的斷裂帶，這些斷裂帶鮮有地震發(fā)生；低值區(qū)對應(yīng)于活動性較弱的斷裂帶，在區(qū)域地殼運(yùn)動劇烈的背景下，在這些活動性相對較弱的斷層上易形成應(yīng)變能積累，因而會發(fā)生地震。區(qū)域絕大多數(shù)地震都發(fā)生在最大剪應(yīng)變的低值區(qū)。

關(guān)鍵詞：青藏高原東緣； 多面函數(shù)； 面應(yīng)變； 最大剪應(yīng)變； 地震

收稿日期：*2013-04-15

基金項(xiàng)目：國家自然科學(xué)基金(41374028, 41274083, 41304013)；國土資源大調(diào)查項(xiàng)目(1212010914015)；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)(CHD2012TD004)

作者簡介：姜永濤(1985-)，男，山東菏澤人，博士生，主要研究方向?yàn)榈貧ば巫兣c地球動力學(xué).E-mail：1212203jiang@sina.com

通訊作者：張永志(1961-)，男，教授，博士生導(dǎo)師，主要從事地殼形變的教學(xué)與研究.E-mail：cadxzyz@263.net

中圖分類號:P227文獻(xiàn)標(biāo)志碼：A

DOI:10.3969/j.issn.1000-0844.2015.01.0152

EvolutionCharacteristicsofCrustalStrainatthe

EasternMarginoftheTibetanPlateau

JIANGYong-tao1，2，ZHANGYong-zhi1， WU Ran1， WANG Shuai1

(1.Chang’an University，Xi’an，Shaanxi710000，China； 2.Nanyang Normal University，Nanyang，Henan473061，China)

Abstract：A study on the evolution of crustal strain in the eastern Tibetan Plateau not only provides a perspective regarding regional crustal strain characteristics but is also very useful for analyzing the relationship between earthquakes and strain rates.Here，we first calculate the crustal strain in spherical coordinates around the eastern Tibetan Plateau in different periods based on a multi-surface function using regional GPS velocities for the periods 2001 to 2004，2004 to 2007，and 1999 to 2013.Then，we study the characteristics of the evolution of surface strain and maximum shear strain in this area，analyzing its relationship with the large earthquakes (>MS6.0) that occurred during each period. It was determined that during the years 2001 to 2004，the eastern Tibetan Plateau generally suffered surface compression.The value of the surface compression of the southern part of the Longmenshan sub-block was lower than the northern part.The Sichuan-Yunnan rhombus-shaped block generally suffered surface compression，but the value of the Baoshan sub-block’s surface compression was lower than the other three sub-blocks.From the distribution and focal mechanisms of the earthquakes (>MS6.0)，we found that the three strike-slip earthquakes between the Baoshan and Dianzhong sub-blocks might be a response to the considerable difference in surface compression between the two sub-blocks.The normal-slip earthquake in the Yajiang sub-block might be a response to the difference of the N-S orientation in the surface compression of the sub-block.During the years 2004 to 2007， the emergence of two areas of strong surface compression formed a “strong-weak-strong” pattern of compression along the Longmenshan Fault，which could have triggered the May 12，2008，Wenchuan MS8.0 earthquake.The surface strain in the Sichuan-Yunnan rhombus-shaped block changed dramatically from 2004 to 2007；this may have had some impact on the occurrence of earthquakes after 2008. From the surface strain rate of 1999 to 2013，we found the evidence of strain responses to the Wenchuan earthquake，manifested as the uniform surface compression around the Xianshuihe faults.Generally，the distribution of surface strain rate agrees well with the scope of the sub-blocks，but the distribution of maximum shear strain rate agrees with the fault systems among the sub-blocks.The areas of high values of maximum shear strain rate correspond to the fault zones of strong tectonic activity, where it is difficult for the strain accumulation to grow and no large earthquakes occur.From the distribution of maximum shear strain rate during 2001 to 2004，we found that most earthquakes were located around the area of low values.This means that lower the maximum shear strain rate，greater the strain accumulation，and higher the risk of large earthquakes.During 2004 to 2007，the low-value area of maximum shear strain rate was located around the Longmenshan Fault，as was the situation during 2001 to 2004.This means that the Longmenshan Fault is not active and it can accumulate considerable strain energy over a very long period，which is an important requirement for the occurrence of the MS8.0 Wenchuan earthquake.

Keywords：easternmarginoftheTibetanPlateau；multi-surfacefunction；surfacestrain；maximumshearstrain；earthquake

0引言

受印度板塊北推碰撞歐亞大陸、青藏高原NE向擠壓和向東擠出動力環(huán)境控制，青藏高原東緣地區(qū)地質(zhì)構(gòu)造復(fù)雜、地形地貌反差顯著、深淺構(gòu)造活動強(qiáng)烈、地震發(fā)生頻度高且強(qiáng)度大，是國內(nèi)外地學(xué)專家研究的熱點(diǎn)[1]。王雙緒等[1]利用1999—2011年GPS區(qū)域站觀測資料和1970—2011年水準(zhǔn)測量資料研究了青藏高原東緣現(xiàn)今三維地殼運(yùn)動的時空分布特征及其運(yùn)動機(jī)制。ZhangPZ[2]對南北地震帶活動斷裂、GPS地殼形變特征和區(qū)域地球物理結(jié)構(gòu)進(jìn)行了綜合分析，探究了與南北地震帶構(gòu)造變形相關(guān)的構(gòu)造樣式、運(yùn)動模式以及深部物質(zhì)變動過程；JiangGY等[3]利用最小二乘配置方法，對青藏高原主要走滑斷裂的現(xiàn)今構(gòu)造應(yīng)力率進(jìn)行了計算，并通過斷層應(yīng)力率與滑動速率之間的關(guān)系探討了地震復(fù)發(fā)周期等相關(guān)問題。李杰等[4]利用2008年新疆于田地震的震前GPS速度場，研究了震前斷層滑移速率和應(yīng)變特征及其與于田地震的關(guān)系。張希等[5]利用1991—1999年、1999—2001年川滇地區(qū)GPS水平運(yùn)動資料，基于改進(jìn)的非震負(fù)位錯模型研究了川滇地區(qū)的應(yīng)變場演化特征及其與強(qiáng)震的關(guān)系。蔣鋒云等[6]基于平面均勻彈性塊體模型，利用1999—2007年和2009—2011年兩個時段GPS觀測到的地殼水平運(yùn)動場對川滇地區(qū)地殼塊體應(yīng)變狀態(tài)及其隨時間的變化特征進(jìn)行研究。吳嘯龍等[7]基于多面函數(shù)擬合，計算了橢球坐標(biāo)系下青藏高原及周邊的地應(yīng)變特征。對于青藏高原東南緣地殼形變機(jī)理，ShenF等[8]和RoydenL[9]提出了三維塑性“管道流”模式，認(rèn)為該區(qū)深部地幔存在右旋剪切，同時在黏滯性較弱的下地殼流動拖拽作用下，淺層地殼向東重力擴(kuò)展的結(jié)果。

本文利用青藏高原東緣2001—2004年、2004—2007年和1999—2013年的GPS水平速率觀測結(jié)果，基于多面函數(shù)擬合，計算了球面坐標(biāo)系下青藏高原東緣不同時期的面應(yīng)變和最大剪應(yīng)變，分析了地應(yīng)變的時空演化特征，結(jié)合不同時期發(fā)生的中強(qiáng)以上地震(MS>6.0)，研究了期間大震分布與地應(yīng)變時空演化特征的關(guān)系。

1數(shù)據(jù)與方法

1.1GPS數(shù)據(jù)資料

本文用到2001—2004年[10]、2004—2007年南北地震帶和1999—2013年[11]的GPS觀測數(shù)據(jù)。由于大震同震位移對速度場影響很大，對2001—2004年GPS數(shù)據(jù)進(jìn)行了2001年昆侖山口西MS8.1地震同震位移建模消除[10]；對1999—2013年GPS數(shù)據(jù)進(jìn)行了2001年昆侖山口西MS8.1、2008年于田MS7.3、2008年汶川MS8.0等地震同震位移建模消除[11]。消除方法是：利用地震斷層滑動模型預(yù)估每個GPS站點(diǎn)的同震位移及其不確定性，對同震水平位移預(yù)估值大于3mm的站點(diǎn)，在求解其站速度時增設(shè)了同震位移三參數(shù)[10-11]，同震位移參數(shù)以其模型預(yù)估值為初值，并施加不確定性約束。因此以上三期GPS數(shù)據(jù)均能以可靠精度反映地殼真實(shí)變形信息。

1.2多面函數(shù)理論模型

多面函數(shù)法是美國Hardy教授[12]于1977年提出的。利用多面函數(shù)法求算地應(yīng)變的基本思路，首先用多個光滑曲面(核函數(shù))以加權(quán)求和方式擬合GPS觀測位移(或速率)，然后由該擬合關(guān)系求算GPS觀測位移對經(jīng)度、緯度的偏微分，最后利用這些偏微分求算應(yīng)變率張量和應(yīng)變參數(shù)。利用多面函數(shù)擬合可采用整體擬合的方式[13]。多面函數(shù)的基本模型為

式中，Q(x，y；xj，yj)為核函數(shù)；u為所取節(jié)點(diǎn)的個數(shù)；αj(j=1，2，…，u)為待估擬合參數(shù)。

在地殼形變分析中核函數(shù)一般選取如式(2)：

其中，(x，y)為GPS點(diǎn)位坐標(biāo)；(xj，yj)為建立多面函數(shù)模型時所選擇的結(jié)點(diǎn)坐標(biāo)；δ為光滑因子；β為核函數(shù)冪次，改變δ、β的值可以改變核函數(shù)形狀，從而影響擬合效果。

利用式(1)可分別建立GPS站速度EW向分量uλ和SN向分量uφ的運(yùn)動模型：

大區(qū)域應(yīng)變計算必須使用球坐標(biāo)[14]，因此在多面函數(shù)擬合GPS站速率時選用大地坐標(biāo)。球坐標(biāo)系下的應(yīng)變率張量與站速度對經(jīng)緯度偏導(dǎo)數(shù)之間的關(guān)系[14]：

式中，λ為經(jīng)度；φ為緯度；R為地球平均曲率半徑；h為大地高。由于本文是在球面上對應(yīng)變場進(jìn)行解算，忽略徑向位移對結(jié)果的影響，故令h=0、uh=0。

利用應(yīng)變率張量，依據(jù)下式可以求算應(yīng)變參數(shù)[15]最大主應(yīng)變ε1，最小主應(yīng)變ε2，最大剪應(yīng)變γMax和面膨脹 εArea。

2青藏高原東緣構(gòu)造變形與應(yīng)變場演化特征分析

2.1青藏高原東緣地殼運(yùn)動特征

圖1是1999—2013年穩(wěn)定歐亞板塊參考框架下青藏高原東緣GPS速率場[11]。圖2為青藏高原東緣次級塊體的劃分[16]，灰色細(xì)線條表示區(qū)域斷層，其中：EKL-Fau 東昆侖斷裂，LRB-Fau 龍日壩斷裂，LMS-Fau 龍門山斷裂，PQ-Fau 平武-青川斷裂，GY-Fau 甘孜—玉樹斷裂，XSH-Fau 鮮水河斷裂，ANH-Fau 安寧河斷裂，ZMH-Fau 則木河斷裂，XJ-Fau 小江斷裂，DLS-Fau 大涼山斷裂，GL-Fau甘孜—理塘斷裂，LT-Fau 理塘斷裂，LJ-Fau 麗江斷裂，JSJ-Fau 金沙江斷裂，DZD-Fau 德欽—中甸—大具斷裂，HH-Fau 紅河斷裂，CH-Fau 程海斷裂，QJ-Fau 曲江斷裂；羅馬數(shù)字表示次級塊體：Ⅰ為阿壩，Ⅱ?yàn)辇堥T山，Ⅲ為藏東，Ⅳ為雅江，Ⅴ為香格里拉，Ⅵ為滇中，Ⅶ為保山，Ⅷ為華南，其中，阿壩次級塊體和龍門山次級塊體合稱為馬爾康塊體，雅江次級塊體、香格里拉次級塊體、滇中次級塊體和保山次級塊體東部合稱為川滇菱形塊體。

圖1顯示，以甘孜—玉樹斷裂帶和鮮水河斷裂帶為界，北側(cè)的馬爾康塊體自西向東GPS速率方向由東向轉(zhuǎn)為東南向，大小逐漸減少。在龍門山次級塊體中存在明顯的ES向地殼縮短現(xiàn)象，可能是在青藏高原物質(zhì)東向擠出作用下，脆性上地殼沿低速層(解耦層)EES滑脫[17]，并受到四川盆地堅硬地殼阻擋的結(jié)果。這種上地殼運(yùn)動特征導(dǎo)致了龍門山次級塊體的高地形和壓應(yīng)變在脆性地殼前端龍門山斷裂帶上的積聚[18]，2008年汶川8.0級地震和2013年雅安7.0級地震的發(fā)生可能是在龍門山斷裂帶上長時期積聚的應(yīng)變能的突然釋放。

白色線條為斷層，紅色箭頭為GPS站平均速度，白色圓圈是區(qū)域發(fā)生的6～7級地震 圖1　高原東緣1999-2013年GPS速率場 [11] Fig.1　The GPS velocity field of eastern Tibetan Plateau during 1999—2013.

圖2　青藏高原東緣主要斷層及塊體劃分 [16] Fig.2　The main faults and blocks of eastern Tibetan Plateau [16]

在甘孜—玉樹斷裂、鮮水河斷裂、安寧河斷裂、則木河斷裂和小江斷裂西南側(cè)，GPS速率值呈現(xiàn)繞東喜馬拉雅構(gòu)造結(jié)順時針旋轉(zhuǎn)特征，大小由南向北逐漸減少，且方向具有規(guī)律性，呈現(xiàn)“順時針地殼流”態(tài)勢。結(jié)合區(qū)域地形特征可以發(fā)現(xiàn)，在北側(cè)馬爾康塊體、東側(cè)四川盆地的阻擋作用和西側(cè)東喜馬拉雅山構(gòu)造結(jié)的鍥入作用下，青藏高原物質(zhì)流呈現(xiàn)出明顯的擠出特征；但在28°N以南地區(qū)，GPS速度方向呈現(xiàn)向東、西發(fā)散趨勢，這可能是由于兩側(cè)地殼對高原地殼流的限制作用減弱，在重力作用下坍塌的表現(xiàn)。青藏高原東南緣地殼運(yùn)動特征，Shen F等[8]和Royden L[9]的觀點(diǎn)是在弱黏滯性下地殼的流動拖拽作用下，淺層地殼向東重力擴(kuò)展的結(jié)果。

2.2地應(yīng)變演化特征分析

基于多面函數(shù)法，利用2001—2004年、2004—2007年和1999—2013年青藏高原東緣的GPS觀測數(shù)據(jù)，求算了該區(qū)的面應(yīng)變(圖3)和最大剪應(yīng)變(圖4)(圖3、4中，白色線條為區(qū)域斷層，黑色震源球表示相應(yīng)時期內(nèi)6級以上地震；紅色震源球分別表示汶川MS8.0地震和雅安MS7.0地震)。本節(jié)結(jié)合圖2對青藏高原東緣次級塊體的劃分[16]，詳細(xì)分析了該區(qū)應(yīng)變場演化特征及其與強(qiáng)震的關(guān)系。

圖3(a)為2001—2004年區(qū)域面應(yīng)變分布，可見該時段青藏高原東緣整體上呈現(xiàn)面壓縮特征，只是在阿壩次級塊體內(nèi)部存在較小的面膨脹。龍門山次級塊體表現(xiàn)為面壓縮特征，印證了該塊體存在地殼縮短，塊體南側(cè)的面壓縮大于北側(cè)，說明南側(cè)的地殼縮短速率大于北側(cè)；川滇菱形塊體整體呈現(xiàn)面壓縮特征，但保山次級塊體與滇中次級塊體交界區(qū)域的面壓縮明顯小于其他次級塊體，這種面應(yīng)變差異會引起平行于斷裂帶方向的負(fù)位錯積累，2001—2004年該斷裂帶上發(fā)生的3次走滑型地震可能是對所積累負(fù)位錯的響應(yīng)；在雅江次級塊體中出現(xiàn)的正傾滑性質(zhì)地震位于該塊體面應(yīng)變梯度帶上，說明在雅江次級塊體內(nèi)部存在南北向不同步的面壓縮特征，南部壓縮相對大于北部，因此發(fā)生了正傾滑型地震。

圖3(b)為2004—2007年區(qū)域面應(yīng)變分布，相對于2001—2004年，阿壩次級塊體中的面膨脹量值增大，范圍向西移動；龍門山次級塊體繼續(xù)呈現(xiàn)面壓縮特征，龍門山斷裂帶南北兩端出現(xiàn)了兩個面壓縮增強(qiáng)區(qū)，沿龍門山斷裂呈現(xiàn)出明顯的“強(qiáng)-弱-強(qiáng)”的面壓縮特征，這可能對汶川地震有一定的觸發(fā)作用；川滇菱形塊體內(nèi)部，相對于2001—2004年，保山次級塊體與滇中次級塊體交匯區(qū)域由面壓縮低值區(qū)變?yōu)槊媾蛎浉咧祬^(qū)，而其他次級塊體則由中等面壓縮變?yōu)閺?qiáng)面壓縮，這說明川滇菱形塊體內(nèi)部出現(xiàn)了較強(qiáng)的“震蕩”，從而可能觸發(fā)了2008年之后該區(qū)的多次地震(圖3(c)中黑色震源球)。

圖3(c)為消除了同震位移影響的1999—2013年區(qū)域面應(yīng)變，包括汶川大震震后應(yīng)變調(diào)整的平均面應(yīng)變。圖3(c)與圖3(b)的面應(yīng)變分布的最大不同點(diǎn)是鮮水河斷裂兩側(cè)的面壓縮突破界限，呈現(xiàn)出連通的特點(diǎn)，這可能是汶川地震震后龍門山斷裂帶上的應(yīng)變能大部分釋放，在斷裂帶南端弱應(yīng)變能的情況下，龍門山次級塊體東南向運(yùn)動受到較小阻力，從而加速東南向運(yùn)動，改變了龍門山次級塊體南側(cè)的震前面應(yīng)變狀態(tài)，導(dǎo)致鮮水河兩側(cè)呈現(xiàn)出一致的面壓縮特征。汶川地震后，青藏高原東緣的地震發(fā)生率明顯高于震前，說明汶川大震對積聚在青藏高原東緣多條斷裂帶上的應(yīng)變能釋放有一定的觸發(fā)作用。

圖3　青藏高原東緣面應(yīng)變演化特征 Fig.3　Evolution characteristics of surface strain in eastern Tibetan Plateau

最大剪應(yīng)變是兩個最大主應(yīng)變之差，其大小反映了地殼構(gòu)造不穩(wěn)定的程度，剪應(yīng)變越大構(gòu)造活動性越強(qiáng)。面應(yīng)變的分布特征大致與塊體范圍對應(yīng)，而從最大剪應(yīng)變的分布特征上看不出塊體范圍。最大剪應(yīng)變高值區(qū)對應(yīng)于構(gòu)造活動性較強(qiáng)的斷裂帶，斷裂帶的活動性越強(qiáng)，在該斷裂帶上就不易形成應(yīng)變能積累，因而很少會發(fā)生地震。從圖4(a)可以看出2001—2004年，青藏高原東緣最大剪應(yīng)變高值區(qū)位于甘孜—玉樹斷裂、鮮水河斷裂帶附近，表明這兩個斷裂的構(gòu)造活動性較強(qiáng)，不利于應(yīng)變能積累。2001—2004年的6級以上地震多發(fā)生在最大剪應(yīng)變的低值區(qū)，印證了構(gòu)造活動較弱的斷層易于應(yīng)變能積累，從而會有長時期積累的應(yīng)變能突然釋放，即發(fā)生地震。

圖4(b)為2004—2007年區(qū)域最大剪應(yīng)變分布。可見該時間段內(nèi)的最大剪應(yīng)變高值區(qū)分布于鮮水河斷裂帶西北段、甘孜—理塘斷裂帶、金沙江斷裂、安寧河斷裂、小江斷裂和紅河斷裂，說明這些斷裂的構(gòu)造活動性在2004—2007年較強(qiáng)。而龍門山斷裂帶從2001—2004年持續(xù)了最大剪應(yīng)變低值特征，說明龍門山斷裂帶活動性不明顯，且在2004—2007年斷裂帶中北段呈現(xiàn)最大剪應(yīng)變最低值，這可能是汶川大震前的應(yīng)變異常。

圖4(c)為1999—2013年平均最大剪應(yīng)變分布。可以看出，1999—2013年甘孜—玉樹斷裂、鮮水河斷裂帶、甘孜—理塘斷裂帶、金沙江斷裂、安寧河斷裂、小江斷裂北段呈現(xiàn)出強(qiáng)最大剪應(yīng)變特征，說明這些斷裂帶上構(gòu)造活動性較強(qiáng)，相應(yīng)的，這些斷裂帶鮮有地震發(fā)生。區(qū)域6級以上地震的發(fā)震位置大部分都處于最大剪應(yīng)變的低值區(qū)。

3結(jié)論與展望

本文計算了球面坐標(biāo)系下青藏高原東緣不同時期的面應(yīng)變和最大剪應(yīng)變，分析了地應(yīng)變的時空演化特征，結(jié)合不同時期發(fā)生的中強(qiáng)地震(MS>6.0)，研究了大震分布與地應(yīng)變時空演化特征的關(guān)系，主要結(jié)論如下：

(1) 2001—2004年保山次級塊體與滇中次級塊體邊界帶上發(fā)生的3次走滑型地震可能是對由面壓縮差異導(dǎo)致邊界斷裂上積聚負(fù)位錯的響應(yīng)，雅江次級塊體中部的正傾滑地震是對塊體南北不同步面壓縮的響應(yīng)；2004—2007年龍門山斷裂呈現(xiàn)出明顯的“強(qiáng)-弱-強(qiáng)”的面壓縮特征；從1999—2013年區(qū)域面應(yīng)變圖中可以明顯發(fā)現(xiàn)汶川大震的應(yīng)變響應(yīng)，即鮮水河斷裂兩側(cè)的面壓縮突破界限，呈現(xiàn)出連通的態(tài)勢。

圖4　青藏高原東緣最大剪應(yīng)變演化特征 Fig.4　Evolution characteristics of maximum shear strain in eastern Tibetan Plateau

(2) 面應(yīng)變的分布特征大致與塊體范圍對應(yīng)，而最大剪應(yīng)變的分布特征上看不出塊體范圍。2001—2004年區(qū)域6級以上地震多發(fā)生在最大剪應(yīng)變的低值區(qū)；2004—2007年龍門山斷裂帶延續(xù)了2001—2004年的最大剪應(yīng)變低值特征，且在斷裂帶中北段出現(xiàn)最大剪應(yīng)變最低值，可能是汶川大震的震前應(yīng)變異常。

多面函數(shù)法建立地應(yīng)變的精度一定程度上依賴于核函數(shù)和結(jié)點(diǎn)的選取，而核函數(shù)和結(jié)點(diǎn)的選取往往是依據(jù)一定的經(jīng)驗(yàn)，通過反復(fù)試算得到的。如何通過一定的算法，快速、精準(zhǔn)地反演得到與觀測值(位移)擬合效果最佳的核函數(shù)和結(jié)點(diǎn)應(yīng)是進(jìn)一步研究的問題；再者，由GPS數(shù)據(jù)獲得的地應(yīng)變是地球內(nèi)部應(yīng)力在地表的響應(yīng)，利用GPS觀測數(shù)據(jù)或其他地球物理數(shù)據(jù)獲取地球內(nèi)部，特別是活動斷裂帶內(nèi)部的應(yīng)力分布及其演化特征，對認(rèn)識區(qū)域變形構(gòu)造背景、預(yù)測強(qiáng)震危險區(qū)具有非常重要的意義。

致謝：感謝審稿專家提出的寶貴意見！

參考文獻(xiàn)(References)

[1]王雙緒, 蔣鋒云, 郝明，等. 青藏高原東緣現(xiàn)今三維地殼運(yùn)動特征研究[J].地球物理學(xué)報，2013，56(10)：3334-3345.

WANG Shuang-xu，JIANG Feng-yun，HAO Ming，et al.Investigation of Feature of Present 3D Crustal Movement in Eastern Edge of Tibet Plateau[J].Chinese J Geophys，2013，56(10)：3334-3345.(in Chinese)

[2]Zhang P Z.A Review on Active Tectonics and Deep Crustal Processes of the Western Sichuan Region，Eastern Margin of the Tibetan Plateau[J].Tectonophysics，2013，584：7-22.

[3]Jiang G Y, Xu C J, Wen Y M, et al.Contemporary Tectonic Stressing Rates of Major Strike-slip Faults in the Tibetan Plateau from GPS Observations Using Least-squares collocation[J].Tectonophysics，2014，615：85-95.

[4]李杰，王曉強(qiáng)，王琪，等.2008 年于田 7.3 級地震前西昆侖地形變的 GPS 初步研究[J].西北地震學(xué)報，2012，34(3)：256-263.

LI Jie，WANG Xiao-qiang，WANG Qi，et al.Crustal Deformation in Western Kunlun Tectonic Zone from GPS Measurements before the YutianMS7.3 Earthquake in 2008[J].Northwestern Seismological Journal，2012，34(3)：256-263.(in Chinese)

[5]張希，江在森，王琪，等.川滇地區(qū)地殼水平運(yùn)動特征與強(qiáng)震關(guān)系研究[J].大地測量與地球動力學(xué)，2003，23(3)：35-41.

ZHANG Xi，JIANG Zai-sen，WANG Qi，et al.Recent Crust Motion Features in Sichuan-Yunnan and Their Relationship to Strong Earthquakes[J].Journal of Geodesy and Geodynamics，2003，23(3)：35-41.(in Chinese)

[6]蔣鋒云，朱良玉，王雙緒，等.青藏東北緣塊體水平運(yùn)動特征研究[J].大地測量與地球動力學(xué)，2013，33(2)：13-17.

JIANG Feng-yun，ZHU Liang-yu，WANG Shuang-xu，et al.Horizontal Movement Characteristics of Northeastern Margin of Qinghai—Tibet Crustal Blocks[J].Journal of Geodesy and Geodynamics，2013，33(2)：13-17. (in Chinese)

[7]吳嘯龍，楊志強(qiáng)，武繼峰，等.基于多面函數(shù)法的青藏高原應(yīng)變特征分析[J].大地測量與地球動力學(xué)，2013，33(4)：17-21.

WU Xiao-long，YANG Zhi-qiang，WU Ji-feng，et al.Analysis on Qinghai—Tibet Plateau Strain Field Based on Multi-surface Function[J].Journal of Geodesy and Geodynamics，2013，33(4)：17-21.(in Chinese)

[8]Shen F，Royden L H，Burchfiel B C.Large-scale Crustal Deformation of the Tibetan Plateau[J].Journal of Geophysical Research：Solid Earth (1978—2012)，2001，106(B4)：6793-6816.

[9]Royden L.Coupling and Decoupling of Crust and Mantle in Convergent Orogens：Implications for Strain Partitioning in the Crust[J].Journal of Geophysical Research：Solid Earth (1978—2012)，1996，101(B8)：17679-17705.

[10]牛之俊，王敏，孫漢榮，等.中國大陸現(xiàn)今地殼運(yùn)動速度場的最新觀測結(jié)果[J].科學(xué)通報，2005，50(8)：839-840.

NIU Zhi-jun，WANG Min，SUN Han-rong，et al.Contemporary Velocity Field of Crustal Movement of Chinese Mainland from Global Positioning System Measurements[J].Chinese Science Bulletin，2005，50(8)：839-840 .(in Chinese)

[11]Liang S M，Gan W J，Shen C Z，et al.Three-dimensional Velocity Field of Present-day Crustal Motion of the Tibetan Plateau Derived from GPS Measurements[J].Journal of Geophysical Research：Solid Earth，2013，118(10)：5722-5732.

[12]Hardy R L.The Mpplication of Multi-quadric Equations and Point Mass Anomaly Models to Crustal Movement Studies.Report National Research Council[J].National Academy of Science，Washington，DC，1977.

[13]陳光保.GPS 技術(shù)監(jiān)測地殼運(yùn)動的新進(jìn)展[J].西北地震學(xué)報，2009，31(3)：302-307.

CHEN Guang-bao.Recent Progress in Research on Monitoring Crustal Movement and Deformation Using GPS[J].Northwestern Seismological Journal，2009，31(3)：302-307.(in Chinese)

[14]石耀霖，朱守彪.用 GPS 位移資料計算應(yīng)變方法的討論[J].大地測量與地球動力學(xué)，2006，26(1)：1-8.

SHI Yao-lin，ZHU Shou-biao.Discussion on Method of Calculating Strain with GPS Displacement Data[J].Journal of Geodesy and Geodynamics ，2006，26(1)：1-8.(in Chinese)

[15]武艷強(qiáng)，江在森，楊國華，等.利用最小二乘配置在球面上整體解算 GPS 應(yīng)變場的方法及應(yīng)用[J].地球物理學(xué)報，2009，52(7)：1707-1711.

WU Yan-qiang，JIANG Zai-sen，YANG Guo-hua，et al.The Application and Method of GPS Strain Calculation in Whole Mode Using Least-squares Collocation in Sphere Surface[J].Chin J Geophys，2009，52(7)：1707-1711.(in Chinese)

[16]程佳，徐錫偉，甘衛(wèi)軍，等.青藏高原東南緣地震活動與地殼運(yùn)動所反映的塊體特征及其動力來源[J].地球物理學(xué)報，2012，55(4)：1198-1212.

CHENG Jia，XU Xi-wei，GAN Wei-jun，et al.Block Model and Dynamic Implication from the Earthquake Activities and Crustal Motion in the Southeastern Margin of Tibetan Plateau[J].Chinese Journal of Geophysics，2012，55(4)：1198-1212.(in Chinese)

[17]Hubbard J，Shaw J H.Uplift of the Longmen Shan and Tibetan Plateau，and the 2008 Wenchuan (M=7.9) Earthquake[J].Nature，2009，458(7235)：194-197.

[18]滕吉文，白登海，楊輝，等.2008汶川MS8.0地震發(fā)生的深層過程和動力學(xué)響應(yīng)[J].地球物理學(xué)報，2009，52(1)：260-276.

TENG Ji-wen，BAI Deng-hai，YANG Hui，et al.Deep Processes and Dynamic Responses Associated with the WenchuanMS8.0 Earthquake of 2008[J].Chinese Journal of Geophysics，2009，52(1)：260-276.(in Chinese)