趙 靜 武艷強(qiáng) 江在森 牛安福劉 杰 王麗鳳 魏文薪

1）中國(guó)北京 100045 中國(guó)地震臺(tái)網(wǎng)中心

2）中國(guó)北京 100036 中國(guó)地震局地震預(yù)測(cè)研究所

引言

2013年4月20日，在巴顏喀拉與華南地塊交界的北東向龍門(mén)山斷裂帶南段發(fā)生了四川蘆山MS7．0地震，這是繼2008年5月12日汶川MS8．0地震后，龍門(mén)山斷裂帶上發(fā)生的又一次強(qiáng)震，震源機(jī)制結(jié)果顯示蘆山地震是一次以逆沖為主，兼有少量右旋走滑分量的地震（陳運(yùn)泰等，2013；劉超等，2013）．曾祥方等（2013）利用近遠(yuǎn)震數(shù)字波形記錄，探討了蘆山地震的震源機(jī)制解．通過(guò)分析波形反演中深度和震源區(qū)地震波速度模型對(duì)斷層面傾角的影響，并結(jié)合余震機(jī)制解，確認(rèn)該地震是一個(gè)高角度逆沖型發(fā)震構(gòu)造產(chǎn)生的地震．張勇等（2013）利用全球地震臺(tái)網(wǎng)的遠(yuǎn)震地震波數(shù)據(jù)分析認(rèn)為，蘆山地震包括兩次破裂子事件，都發(fā)生在斷層面上震源附近，沒(méi)有表現(xiàn)出明顯的破裂方向性，也沒(méi)有發(fā)現(xiàn)大規(guī)模出露地表的破裂跡象．王衛(wèi)民等（2013）利用遠(yuǎn)場(chǎng)體波資料反演了蘆山MW6．7地震的震源破裂過(guò)程與地震斷層面上的滑動(dòng)量分布，結(jié)果顯示蘆山逆沖型地震的最大滑動(dòng)量為159cm．

汶川地震后，許多地震工作者對(duì)汶川地震的震前變形（杜方等，2009；江在森等，2009；李延興等，2009；武艷強(qiáng)等，2011；趙靜等，2012）、同震位移（徐錫偉等，2008；Shenetal，2009；王敏，2009；許才軍等，2009；Xuetal，2009；Wangetal，2011）、動(dòng)力機(jī)制（滕吉文等，2008；張培震等，2008，2009；朱介壽，2008）、破裂過(guò)程（陳運(yùn)泰等，2008；Nishimura，Yagi，2008；王衛(wèi)民等，2008；張勇等，2008；趙翠萍等，2009）和震后斷層應(yīng)力變化（單斌等，2009；萬(wàn)永革等，2009；邵志剛等，2010）等進(jìn)行了系統(tǒng)的研究，獲得了許多有意義的成果．多數(shù)研究結(jié)果顯示汶川地震的破裂擴(kuò)展類(lèi)型為自初始破裂點(diǎn)沿龍門(mén)山斷裂帶向北東方向單側(cè)破裂，且震后余震特別是MS≥5．0余震都發(fā)生在主震位置的北東方向斷層上，而龍門(mén)山斷裂帶南段沒(méi)有發(fā)生破裂且沒(méi)有發(fā)生較大余震，所積累的應(yīng)變能并沒(méi)有得到釋放．由于汶川地震后龍門(mén)山地震空區(qū)的剩余部分位于龍門(mén)山斷裂西南段上，因此多位學(xué)者對(duì)南段的地震危險(xiǎn)性也進(jìn)行了關(guān)注和研究（陳運(yùn)泰等，2008；聞學(xué)澤等，2009）．單斌等（2009）根據(jù)地震的靜態(tài)觸發(fā)原理，利用彈性位錯(cuò)理論和分層地殼模型計(jì)算了汶川地震引起的周邊斷層同震應(yīng)力變化，認(rèn)為彭縣—灌縣斷裂、北川—映秀斷裂南段所積累的構(gòu)造應(yīng)力未充分釋放，同震庫(kù)侖應(yīng)力的加載有可能在今后觸發(fā)該區(qū)較強(qiáng)余震；萬(wàn)永革等（2009）以汶川地震的破裂為驅(qū)動(dòng)源，計(jì)算了該地震造成周?chē)鷶鄬由系撵o態(tài)庫(kù)侖破裂應(yīng)力變化，結(jié)果表明汶川地震的發(fā)生使得龍門(mén)山斷裂最南端的庫(kù)侖破裂應(yīng)力增加；邵志剛等（2010）基于Burgers體黏滯松弛模型計(jì)算了汶川地震引起的庫(kù)侖應(yīng)力動(dòng)態(tài)演化，認(rèn)為龍門(mén)山斷裂南段的庫(kù)侖應(yīng)力顯著增加，存在較高的地震危險(xiǎn)性；聞學(xué)澤等（2011）通過(guò)對(duì)區(qū)域動(dòng)力學(xué)背景與地震發(fā)生關(guān)系的研究認(rèn)為，巴顏喀拉塊體東邊界斷裂系統(tǒng)有發(fā)生MS≥7．0地震的危險(xiǎn)性，其中包括龍門(mén)山斷裂帶南段（汶川以南—康定）．

利用GPS觀測(cè)資料反演斷層不同位置的閉鎖程度和閉鎖深度可以為判斷斷裂帶未來(lái)一段時(shí)期內(nèi)地震危險(xiǎn)性提供參考．隨著汶川地震后陸態(tài)網(wǎng)絡(luò)二期GPS連續(xù)站和流動(dòng)站逐步投入使用，龍門(mén)山斷裂帶及其附近地區(qū)GPS站點(diǎn)的分布密度得以較大提高，能夠?yàn)檎囱萦?jì)算提供更為有效的數(shù)據(jù)約束．本文利用1999—2007年和2009—2011年GPS水平速度場(chǎng)結(jié)果，采用DEFNODE負(fù)位錯(cuò)反演程序估算了蘆山地震前龍門(mén)山斷裂帶閉鎖程度和垂直斷層的滑動(dòng)虧空速率的三維變化情況，并綜合剖面結(jié)果分析了汶川地震前后斷層深淺部變形的變化特征．

1 DEFNODE負(fù)位錯(cuò)反演方法基本原理

McCaffrey（2002，2007）公布了DEFNODE負(fù)位錯(cuò)反演程序．該程序假定塊體內(nèi)部點(diǎn)的運(yùn)動(dòng)為塊體旋轉(zhuǎn)、塊體內(nèi)部整體均勻應(yīng)變及塊體邊界由于斷層閉鎖產(chǎn)生的滑動(dòng)虧空而引起的地表彈性變形之和．若塊體內(nèi)部不存在整體均勻應(yīng)變，其原理如下（McCaffrey，2002）：

式中X為GPS測(cè)站的位置，Vi（X）為測(cè)站X的速度，B為塊體數(shù)目，Δb為b塊體模型區(qū)域的子集（如果點(diǎn)X在塊體b范圍內(nèi)H＝1，否則H＝0），i為速度分量索引（x或者y），i為i方向的單位矢量，RΩb為塊體b相對(duì)于參考框架的歐拉極，hΩf＝hΩR－fΩR為斷層下盤(pán)f相對(duì)于上盤(pán)h的歐拉極，F(xiàn)為斷層數(shù)量，Nk為定義斷層k的節(jié)點(diǎn)數(shù)，Xnk為斷層k上節(jié)點(diǎn)n的位置，φnk為斷層k上節(jié)點(diǎn)n的閉鎖程度，j為斷層面上j方向的單位矢量（走向方向或垂直等深線(xiàn)方向），Gij（X，Xnk）為響應(yīng)函數(shù)，表示斷層面上節(jié)點(diǎn)Xnk處j方向上的單位滑動(dòng)速率引起的地表位置X處i方向的速度分量（McCaffrey，2002；趙靜等，2012）．

式中θ和λ分別表示余緯度和余經(jīng)度，R為區(qū)域幾何中心的平均曲率半徑，θ0為區(qū)域幾何中心余緯度，Δθ和Δλ分別為測(cè)點(diǎn)到區(qū)域幾何中心的余緯度差和余經(jīng)度差，V為速度分量，˙ε為應(yīng)變率分量．

若塊體內(nèi)部存在整體均勻應(yīng)變，應(yīng)變率的計(jì)算采用Savage等（2001）給出的公式．此模型為式（1）加上式（2）計(jì)算得到的，式（2）表示塊體內(nèi)部均勻應(yīng)變引起的速度量值的大小．

利用GPS水平速度場(chǎng)數(shù)據(jù)反演塊體旋轉(zhuǎn)、塊體內(nèi)部均勻應(yīng)變和斷層閉鎖程度等參數(shù)的過(guò)程中，參數(shù)擬合的不符值可由式（3）表征．

式中n為所有觀測(cè)數(shù)據(jù)的數(shù)量，dof為自由度（所有觀測(cè)數(shù)據(jù)數(shù)量－自由參數(shù)數(shù)量），ri為觀測(cè)數(shù)據(jù)殘差，σi為數(shù)據(jù)標(biāo)準(zhǔn)差，f為數(shù)據(jù)誤差權(quán)重因子，它可以平衡各種數(shù)據(jù)（如GPS水平速度場(chǎng)和GPS垂直速度場(chǎng)）的影響（McCaffrey，2002；趙靜等，2012）．

為得到使χ2n達(dá)到最小的一組最佳參數(shù)值，本反演程序使用網(wǎng)格搜索和模擬退火算法求解，并且對(duì)每個(gè)塊體內(nèi)部是否存在整體均勻應(yīng)變進(jìn)行檢驗(yàn)．檢驗(yàn)完成后，通過(guò)不斷改變f值大小進(jìn)行逐步擇優(yōu)，尋求最滿(mǎn)意的f值，進(jìn)而保證χ2n≈1，此時(shí)模型能夠最佳地解釋觀測(cè)數(shù)據(jù)．

2 龍門(mén)山斷裂帶計(jì)算結(jié)果

2．1 GPS數(shù)據(jù)與斷層模型

根據(jù)地質(zhì)結(jié)果與現(xiàn)今大地測(cè)量結(jié)果，本文研究區(qū)域以龍門(mén)山斷裂帶與鮮水河斷裂帶為界，包含巴顏喀拉、華南和川滇3個(gè)地塊的部分區(qū)域（圖1）．1999—2007年有98個(gè)GPS水平速度參與反演（趙靜等，2012）；2009—2011年有159個(gè)GPS測(cè)點(diǎn)參與反演（圖1），其中巴顏喀拉地塊45個(gè)、華南地塊90個(gè)、川滇地塊24個(gè)，平均誤差為2．1mm／a．采用GAMIT／GLOBK（Herringetal，2006a，b）和QOCA（Dongetal，1998）軟件解算得到速度場(chǎng)．

反演過(guò)程中，沿龍門(mén)山斷裂走向共有6條等深線(xiàn)，深度依次為0．1，6，12，16，21和24km，每條等深線(xiàn)上有10個(gè)節(jié)點(diǎn)，節(jié)點(diǎn)間距離為45—60km；根據(jù)文獻(xiàn)（王衛(wèi)民等，2008；徐錫偉等，2008；張培震等，2008；朱介壽，2008；杜方等，2009；Xuetal，2009；Wanget al，2011）的研究結(jié)果，每?jī)膳殴?jié)點(diǎn)之間的斷層傾角從地表至深部依次為55°，50°，20°，7°和7°；斷裂總長(zhǎng)度為475km，斷層面總寬度為92km（圖2）．

2．2 不同模型對(duì)比與最優(yōu)模型

由于汶川地震前龍門(mén)山斷裂帶（特別是其中北段）處于強(qiáng)閉鎖狀態(tài)，積累了足夠的彈性應(yīng)變能，因此為了使斷層閉鎖產(chǎn)生的彈性應(yīng)變能得到充分體現(xiàn)，在處理1999—2007年GPS數(shù)據(jù)時(shí)，認(rèn)為3個(gè)地塊內(nèi)只存在整體旋轉(zhuǎn)運(yùn)動(dòng)和由于斷層閉鎖而產(chǎn)生的彈性應(yīng)變，不存在整體均勻應(yīng)變，即采用式（1）的原理來(lái)進(jìn)行反演計(jì)算，然后通過(guò)反演殘余速度計(jì)算殘差應(yīng)變率．通過(guò)大量試算得到了最優(yōu)模型，其中GPS水平速度場(chǎng)數(shù)據(jù)誤差權(quán)重因子f取為1．27，χ2n＝0．996（觀測(cè)值個(gè)數(shù)為196，自由度為169）．

圖1 研究區(qū)GPS站2009—2011年水平運(yùn)動(dòng)速度場(chǎng)（相對(duì)于華南地塊）與活動(dòng)斷裂分布Fig．1 The horizontal GPS velocity field during 2009—2011（related to South China block）and active faults in the study region

圖2 斷層模型與節(jié)點(diǎn)設(shè)置（a）斷層面結(jié)構(gòu)與節(jié)點(diǎn)位置；（b）斷層剖面結(jié)構(gòu)Fig．2 Fault model and node setting（a）Fault surface geometry and nodes；（b）Vertical section of fault

汶川地震后，龍門(mén)山斷裂帶中北段附近積累的彈性應(yīng)變能得到了釋放，巴顏喀拉地塊內(nèi)靠近中北段附近的GPS測(cè)點(diǎn)運(yùn)動(dòng)速率明顯加快（圖1），在地塊內(nèi)部產(chǎn)生拉張應(yīng)變．因此在處理2009—2011年GPS數(shù)據(jù)時(shí)，對(duì)于地塊內(nèi)部是否存在整體均勻應(yīng)變進(jìn)行了8種試算，最優(yōu)模型為巴顏喀拉地塊和川滇地塊內(nèi)部存在整體均勻應(yīng)變，華南地塊內(nèi)部不存在整體均勻應(yīng)變．GPS水平速度場(chǎng)數(shù)據(jù)誤差權(quán)重因子f取為2．18，χ2n＝1．007（觀測(cè)值個(gè)數(shù)為318，自由度為285）．

3 蘆山地震前龍門(mén)山斷裂地殼變形特征分析

3．1 震前兩個(gè)時(shí)段斷層閉鎖比例與滑動(dòng)虧空速率分布

根據(jù)前面的最優(yōu)模型得到了蘆山地震前兩個(gè)時(shí)段的龍門(mén)山斷裂閉鎖程度（φ）空間分布．圖3為1999—2007年GPS速度場(chǎng)反演結(jié)果（汶川地震的起始破裂深度取19km（劉啟元等，2008））．結(jié)果表明，汶川地震前龍門(mén)山斷裂中北段處于強(qiáng)閉鎖狀態(tài)（閉鎖比例0．99），而且閉鎖深度約為21km；南段只有地表以下16km深度為強(qiáng)閉鎖狀態(tài)，在16—21km深度處閉鎖比例降低為0．62，較中北段同一深度的閉鎖程度稍弱，在21—24km深度處整條斷裂均逐步由閉鎖轉(zhuǎn)變?yōu)槿浠卣鳎ㄚw靜等，2012）．

圖3 1999—2007年龍門(mén)山斷裂閉鎖程度（引自趙靜等，2012）（a）斷層閉鎖地面投影；（b）斷層閉鎖分布Fig．3 Locking level of the Longmenshan fault during 1999—2007（a）Projection of fault locking；（b）Fault locking distribution

圖4為2009—2011年GPS速度場(chǎng)反演結(jié)果．其中蘆山地震的震中（30．286°N，102．956°E）、起始破裂深度（12．3km）和震源機(jī)制參考USGS最初給出的結(jié)果．反演結(jié)果表明，在蘆山地震前，龍門(mén)山斷裂中段和北段已經(jīng)基本解鎖（只有北段0．1—6km深度的閉鎖比例約為0．15），處于蠕滑狀態(tài)，屬于大震后早期的愈合階段；龍門(mén)山斷裂南段仍然是地表以下深度16km內(nèi)為強(qiáng)閉鎖，在深度16—21km處閉鎖比例降低為0．45，較汶川震前的0．62稍有降低，在21—24km深度處逐漸轉(zhuǎn)變?yōu)槿浠卣鳎J山地震的起始破裂點(diǎn)位于龍門(mén)山斷裂南段的強(qiáng)閉鎖位置（圖4b），分析認(rèn)為，由于汶川地震的破裂擴(kuò)展類(lèi)型為自初始破裂點(diǎn)沿龍門(mén)山斷裂向北東方向的單側(cè)破裂，且震后余震特別是MS≥5．0余震都發(fā)生在主震位置的北東方向斷層上，而在汶川地震時(shí)龍門(mén)山斷裂帶南段并沒(méi)有發(fā)生破裂，汶川地震后也沒(méi)有發(fā)生較大余震，因此其所積累的彈性應(yīng)變能并未得到釋放，一直處于較強(qiáng)的閉鎖狀態(tài)，存在發(fā)生中強(qiáng)地震背景．同時(shí)，由于龍門(mén)山斷裂帶南段的閉鎖深度為16km，較汶川地震前中北段的閉鎖深度淺，因此蘆山地震較汶川地震強(qiáng)度低、震級(jí)小、破裂范圍窄．

圖4 2009—2011年龍門(mén)山斷裂閉鎖程度（a）斷層閉鎖地面投影；（b）斷層閉鎖分布Fig．4 Locking level of the Longmenshan fault during 2009—2011（a）Projection of fault locking；（b）Fault locking distribution

圖5給出了龍門(mén)山斷裂兩個(gè)時(shí)段的垂直斷層方向的擠壓滑動(dòng)虧空速率分布．圖5a結(jié)果顯示，汶川地震前龍門(mén)山斷裂南段從地表至深度16km處滑動(dòng)虧空速率約為2．0mm／a，深度16—21km處滑動(dòng)虧空速率降低為1．2mm／a；圖5b結(jié)果顯示，汶川地震后龍門(mén)山斷裂南段從地表至深度16km處滑動(dòng)虧空速率約為10．5mm／a，深度16—21km處滑動(dòng)虧空速率降低為4．8mm／a，深度21—24km處斷裂逐漸轉(zhuǎn)變?yōu)槿浠瑺顟B(tài)，滑動(dòng)虧空值逐漸降低為0．以上結(jié)果表明，汶川地震后龍門(mén)山斷裂南段的滑動(dòng)虧空速率明顯增大，汶川地震的發(fā)生導(dǎo)致了龍門(mén)山斷裂南段應(yīng)變積累速率的提高．結(jié)合圖1速度場(chǎng)結(jié)果分析認(rèn)為，汶川地震后巴顏喀拉地塊東向運(yùn)動(dòng)顯著增強(qiáng)，導(dǎo)致仍然處于閉鎖狀態(tài)的龍門(mén)山斷裂帶南段應(yīng)變積累總量增加，因此汶川地震的發(fā)生加速了蘆山地震的孕育進(jìn)程．

3．2 剖面結(jié)果分析

根據(jù)圖1中剖面位置（剖面長(zhǎng)均為700km、中段剖面A寬180km和南段剖面B寬120km），圖6給出了龍門(mén)山斷裂中段剖面A和南段剖面B的2009—2011年GPS水平觀測(cè)值與模型值的剖面擬合結(jié)果．圖6a結(jié)果表明，位于龍門(mén)山斷裂帶附近一些測(cè)點(diǎn)的垂直斷層運(yùn)動(dòng)分量的擬合殘差較大，原因在于汶川逆沖型地震發(fā)生以后，巴顏喀拉地塊處于震后余滑與松弛狀態(tài)，因此斷層附近GPS測(cè)點(diǎn)的垂直斷層運(yùn)動(dòng)分量比遠(yuǎn)離斷層GPS測(cè)點(diǎn)的垂直斷層運(yùn)動(dòng)分量大很多；其它GPS測(cè)站的模型擬合殘差較小，基本保持在3mm／a以?xún)?nèi)．

圖5 1999—2007年（a）與2009—2011年（b）龍門(mén)山斷裂垂直斷層方向的擠壓滑動(dòng)虧空分布Fig．5 Normal compressional slip deficit rate of the Longmenshan fault during 1999—2007（a）and 2009—2011（b）

圖6 橫跨龍門(mén)山斷裂中段剖面A和南段剖面B的2009—2011年GPS站速度剖面和擬合結(jié)果（相對(duì)于華南地塊）剖面A、B位置見(jiàn)圖1所示．圖中灰色線(xiàn)為平行斷層運(yùn)動(dòng)的模型值，斜率為正表示左旋，反之為右旋；黑色線(xiàn)為垂直斷層運(yùn)動(dòng)的模型值，斜率為正表示拉張，反之為擠壓Fig．6 GPS velocity profiles A （a）and B （b）across the middle and the southern segments of the Longmenshan fault during 2009—2011，as well as model fitting results（relative to South China block）Loctions of profiles Aand Bsee Fig．1．Gray line denotes velocity parallel to fault，positive slope means sinistral slip and negative denotes dextral silp．Black line indicates velocity normal to fault，positive slope indicates extensional slip and negative indicates compressional slip

圖6a結(jié)果表明，龍門(mén)山斷裂中段的逆沖和走滑分量均存在一個(gè)明顯的階躍變化，表現(xiàn)出震后斷裂滑動(dòng)特征，表明中段已經(jīng)解鎖，處于震后蠕滑狀態(tài)；圖6b結(jié)果表明，龍門(mén)山斷裂南段的逆沖與走滑分量均表現(xiàn)為連續(xù)變形特征，斷層變形寬度較大，顯示斷層不能自由滑動(dòng)，仍然處于強(qiáng)閉鎖狀態(tài)．趙靜等（2012）的結(jié)果與本文圖6b對(duì)比可知，龍門(mén)山斷裂帶南段西北側(cè)250km范圍內(nèi)，其垂直斷層的水平縮短速率由汶川震前約2mm／a增大到震后約6mm／a，其平行斷層的右旋扭動(dòng)量由汶川震前5mm／a增大到震后約7mm／a（此處按（32．22N°，104．85E°—29．91N°，102．28E°）方位斷層線(xiàn)投影，而非蘆山地震破裂方向投影），對(duì)比表明汶川地震使龍門(mén)山斷裂帶南段西北側(cè)的巴顏喀拉地塊的逆沖與右旋走滑運(yùn)動(dòng)都有所加強(qiáng)．同時(shí)在汶川地震后，龍門(mén)山斷裂南段仍處于較強(qiáng)的閉鎖狀態(tài)．因此，汶川地震后巴顏喀拉地塊垂直斷層的水平縮短速率和平行斷層的右旋扭動(dòng)量都增加了，進(jìn)而龍門(mén)山斷裂南段的應(yīng)變積累速率也增加了．

4 討論與結(jié)論

1）DEFNODE負(fù)位錯(cuò)反演方法在國(guó)際上已得到了較為廣泛的應(yīng)用（McCaffreyetal，2000，2007；McCaffrey，2002，2005；Wallaceetal，2004a，b，2006），在國(guó)內(nèi)僅有較少文章發(fā)表．筆者在對(duì)靜態(tài)問(wèn)題進(jìn)行模擬的基礎(chǔ)上，驗(yàn)證了該方法在低傾角與高傾角條件下反演結(jié)果均具有高度可靠性（趙靜等，2013）．本文將該方法應(yīng)用于動(dòng)態(tài)資料的反演中，初步證明反演結(jié)果有效．

2）1999—2007年GPS數(shù)據(jù)反演結(jié)果表明，汶川地震前龍門(mén)山斷裂中北段處于強(qiáng)閉鎖狀態(tài)；南段只有地表以下深度16km內(nèi)為強(qiáng)閉鎖，深度16—21km處閉鎖比例降低為0．62，較中北段同一深度的閉鎖程度稍弱；深度21—24km處整條斷裂逐漸轉(zhuǎn)變?yōu)槿浠瑺顟B(tài)．2009—2011年GPS數(shù)據(jù)反演結(jié)果表明，汶川地震后龍門(mén)山斷裂中北段處于蠕滑狀態(tài)，屬于震后愈合階段；龍門(mén)山斷裂南段仍然是地表以下深度16km內(nèi)為強(qiáng)閉鎖，深度16—21km處閉鎖比例降低為0．45，深度21—24km處逐漸轉(zhuǎn)變?yōu)槿浠卣鳎氪ǖ卣鸢l(fā)生在龍門(mén)山斷裂帶中段與南段交匯的位置，震中兩側(cè)斷層閉鎖程度有所差異，地震發(fā)生時(shí)從震中向北東向單向破裂而南段沒(méi)有發(fā)生破裂；蘆山地震發(fā)生在汶川地震后依然處于較高閉鎖狀態(tài)的龍門(mén)山斷裂南段，反演破裂過(guò)程結(jié)果顯示，破裂具有各向同性——從一點(diǎn)開(kāi)始，以近似圓形向四面擴(kuò)展（王衛(wèi)民等，2013；張勇等，2013）．

3）垂直斷層方向的擠壓滑動(dòng)虧空速率結(jié)果顯示，汶川地震前，龍門(mén)山斷裂南段從地表至深度16km處滑動(dòng)虧空速率約為2．0mm／a，深度16—21km處滑動(dòng)虧空速率降低為1．2mm／a；汶川地震后，龍門(mén)山斷裂南段從地表至深度16km處滑動(dòng)虧空速率約為10．5mm／a，深度16—21km處滑動(dòng)虧空速率降低為4．8mm／a．這表明汶川地震后龍門(mén)山斷裂南段的滑動(dòng)虧空速率明顯增大，在沒(méi)有發(fā)生蠕滑的情況下，汶川地震促進(jìn)了其應(yīng)變積累速率的提高，加速了蘆山地震的孕育進(jìn)程．

4）2009—2011年GPS水平觀測(cè)值與模型值的剖面擬合結(jié)果較好．汶川地震后，龍門(mén)山斷裂南段的逆沖和走滑分量均表現(xiàn)為連續(xù)變形特征，斷層變形寬度較大，顯示斷層不能自由滑動(dòng)，仍然處于強(qiáng)閉鎖狀態(tài)；龍門(mén)山斷裂中段的逆沖和走滑分量均存在一個(gè)明顯的階躍變化，表現(xiàn)出滑動(dòng)特征，表明中段已經(jīng)解鎖，處于震后蠕滑狀態(tài)．龍門(mén)山斷裂帶南段西北側(cè)250km范圍內(nèi)，其垂直斷層的水平縮短速率由汶川震前約2mm／a增大到震后約6mm／a；其平行斷層的右旋扭動(dòng)量由汶川地震前5mm／a增大到震后約7mm／a．這表明龍門(mén)山斷裂南段來(lái)自西北側(cè)相對(duì)大尺度加載速率加大了，使得龍門(mén)山斷裂南段應(yīng)變積累速度加快．因此，汶川地震的發(fā)生促進(jìn)了龍門(mén)山斷裂南段應(yīng)變積累速率的提高．

陳運(yùn)泰，許力生，張勇，杜海林，馮萬(wàn)鵬，劉超，李春來(lái)．2008．2008年5月12日汶川特大地震震源特性分析報(bào)告［R／OL］．［2009-01-05］．http：∥www．cea-igp．a(chǎn)c．cn／汶川地震專(zhuān)題／地震情況／初步研究及考察結(jié)果（一）．

陳運(yùn)泰，楊智嫻，張勇，劉超．2013．從汶川地震到蘆山地震［J］．中國(guó)科學(xué)：地球科學(xué)，43（6）：1064-1072．

杜方，聞學(xué)澤，張培震，王慶良．2009．2008年汶川8．0級(jí)地震前橫跨龍門(mén)山斷裂帶的震間形變［J］．地球物理學(xué)報(bào)，52（11）：2729-2738．

江在森，方穎，武艷強(qiáng)，王敏，杜方，平建軍．2009．汶川8．0級(jí)地震前區(qū)域地殼運(yùn)動(dòng)與變形動(dòng)態(tài)過(guò)程［J］．地球物理學(xué)報(bào)，52（2）：505-518．

李延興，張靜華，周偉，胡新康，郭良遷，張中伏．2009．汶川MS8．0地震孕育發(fā)生的機(jī)制與動(dòng)力學(xué)問(wèn)題［J］．地球物理學(xué)報(bào)，52（2）：519-530．

劉超，許力生，陳運(yùn)泰．2013．2013年4月20日四川蘆山7．0級(jí)地震（1．0版）［EB／OL］．［2013-05-01］．http：∥www．cea-igp．a(chǎn)c．cn／tpxw／266824．shtml．

劉啟元，陳九輝，李順成，李昱，郭飚，王峻，齊少華．2008．汶川8．0級(jí)地震：川西流動(dòng)地震臺(tái)陣觀測(cè)數(shù)據(jù)的初步分析［J］．地震地質(zhì)，30（3）：584-596．

單斌，熊熊，鄭勇，刁法啟．2009．2008年5月12日MW7．9汶川地震導(dǎo)致的周?chē)鷶鄬討?yīng)力變化［J］．中國(guó)科學(xué)：D輯，39（5）：537-545．

邵志剛，周龍泉，蔣長(zhǎng)勝，馬宏生，張浪平．2010．2008年汶川MS8．0地震對(duì)周?chē)鷶鄬拥卣鸹顒?dòng)的影響［J］．地球物理學(xué)報(bào)，53（8）：1784-1795，doi：10．3969／j．issn．0001-5733．2010．08．004．

滕吉文，白登海，楊輝，閆雅芬，張洪雙，張永謙，阮小敏．2008．2008汶川MS8．0地震發(fā)生的深層過(guò)程和動(dòng)力學(xué)響應(yīng)［J］．地球物理學(xué)報(bào)，51（5）：1385-1402．

萬(wàn)永革，沈正康，盛書(shū)中，徐曉楓．2009．2008年汶川大地震對(duì)周?chē)鷶鄬拥挠绊懀跩］．地震學(xué)報(bào)，31（2）：128-139．

王敏．2009．基于GPS同震位移場(chǎng)約束反演2008年5·12汶川大地震破裂空間分布［J］．地球物理學(xué)報(bào)，52（10）：2519-2526，doi：10．3969／j．issn．0001-5733．2009．10．010．

王衛(wèi)民，郝金來(lái)，姚振興．2013．2013年4月20日四川蘆山地震震源破裂過(guò)程反演初步結(jié)果［J］．地球物理學(xué)報(bào)，56（4）：1412-1417，doi：10．6038／cjg20130436．

王衛(wèi)民，趙連鋒，李娟，姚振興．2008．四川汶川8．0級(jí)地震震源過(guò)程［J］．地球物理學(xué)報(bào)，51（5）：1403-1410．

聞學(xué)澤，張培震，杜方，龍鋒．2009．2008年汶川8．0級(jí)地震發(fā)生的歷史與現(xiàn)今地震活動(dòng)背景［J］．地球物理學(xué)報(bào)，52（2）：444-454．

聞學(xué)澤，杜方，張培震，龍鋒．2011．巴顏喀拉塊體北和東邊界大地震序列的關(guān)聯(lián)性與2008年汶川地震［J］．地球物理學(xué)報(bào)，54（3）：706-716，doi：10．3969／j．issn．0001-5733．2011．03．010．

武艷強(qiáng)，江在森，楊國(guó)華，劉曉霞，趙靜．2011．汶川地震前GPS資料反映的應(yīng)變率場(chǎng)演化特征［J］．大地測(cè)量與地球動(dòng)力學(xué)，31（5）：20-25．

徐錫偉，聞學(xué)澤，葉建青，馬保起，陳杰，周榮軍，何宏林，田勤儉，何玉林，王志才，孫昭民，馮希杰，于貴華，陳立春，陳桂華，于慎鄂，冉勇康，李細(xì)光，李陳俠，安艷芬．2008．汶川MS8．0地震地表破裂帶及其發(fā)震構(gòu)造［J］．地震地質(zhì)，30（3）：597-629．

許才軍，劉洋，溫?fù)P茂．2009．利用GPS資料反演汶川MW7．9級(jí)地震滑動(dòng)分布［J］．測(cè)繪學(xué)報(bào)，38（3）：195-201．

曾祥方，羅艷，韓立波，石耀霖．2013．2013年4月20日四川蘆山MS7．0地震：一個(gè)高角度逆沖地震［J］．地球物理學(xué)報(bào)，56（4）：1418-1424，doi：10．6038／cjg20130437．

張培震，徐錫偉，聞學(xué)澤，冉勇康．2008．2008年汶川8．0級(jí)地震發(fā)震斷裂的滑動(dòng)速率、復(fù)發(fā)周期和構(gòu)造成因［J］．地球物理學(xué)報(bào)，51（4）：1066-1073．

張培震，聞學(xué)澤，徐錫偉，甘衛(wèi)軍，王敏，沈正康，王慶良，黃媛，鄭勇，李小軍，張竹琪，馬勝利，冉勇康，劉啟元，丁志峰，吳建平．2009．2008年汶川8．0級(jí)特大地震孕育和發(fā)生的多單元組合模式［J］．科學(xué)通報(bào)，54（7）：944-953．

張勇，馮萬(wàn)鵬，許力生，周成虎，陳運(yùn)泰．2008．2008年汶川大地震的時(shí)空破裂過(guò)程［J］．中國(guó)科學(xué)：D輯，38（10）：1186-1194．

張勇，許力生，陳運(yùn)泰．2013．蘆山4·20地震破裂過(guò)程及其致災(zāi)特征初步分析［J］．地球物理學(xué)報(bào)，56（4）：1408-1411，doi：10．6038／cjg20130435．

趙翠萍，陳章立，周連慶，李志雄，康英．2009．汶川MW8．0地震震源破裂過(guò)程研究：分段特征［J］．科學(xué)通報(bào)，54（22）：3475-3482．

趙靜，江在森，武艷強(qiáng)，劉曉霞，魏文薪，李強(qiáng)．2012．汶川地震前龍門(mén)山斷裂帶閉鎖程度和滑動(dòng)虧損分布研究［J］．地球物理學(xué)報(bào)，55（9）：2963-2972，doi：10．6038／j．issn．0001-5733．2012．09．015．

趙靜，江在森，武艷強(qiáng)，劉曉霞，魏文薪，王閱兵，李強(qiáng)，徐晶．2013．Defnode負(fù)位錯(cuò)模型反演結(jié)果的可靠性和穩(wěn)定性分析［J］．大地測(cè)量與地球動(dòng)力學(xué)，33（1）：21-24．

朱介壽．2008．汶川地震的巖石圈深部結(jié)構(gòu)與動(dòng)力學(xué)背景［J］．成都理工大學(xué)學(xué)報(bào)：自然科學(xué)版，35（4）：348-356．

Dong D，Herring T A，King R W．1998．Estimating regional deformation from a combination of space and terrestrial geodetic data［J］．JGeophysRes，72：200-214．

Herring T A，King R W，McClusky S C．2006a．GAMITReferenceManual［M／OL］．Release 10．3．Massachussetts Institute Technology．［2013-04-29］．http：∥www-gpsg．mit．edu／～simon／gtgk／index．htm．

Herring T A，King R W，McClusky S C．2006b．GLOBKReferenceManual［M／OL］．Release 10．3．Massachussetts Institute Technology．［2013-04-29］．http：∥www-gpsg．mit．edu／～simon／gtgk／index．htm．

McCaffrey R．2002．Crustal block rotations and plate coupling：in plate boundary zones［J］．GeodynamicsSeries，30：101-122．

McCaffrey R．2005．Block kinematics of the Pacific-North America plate boundary in the southwestern US from inversion of GPS，seismological，and geologic data［J］．JGeophysRes，110（B7）：B07401，doi：10．1029／2004JB003307．

McCaffrey R．2007．DEFNODEUser′sManual［M／OL］．Version 2007．10．25．［2013-04-29］．http：∥web．pdx．edu／～mccaf／www／defnode／defnode＿071025．html．

McCaffrey R，Long M D，Goldfinger C，Zwick P C，Nabelek J L，Johnson C K，Smith C．2000．Rotation and plate locking at the southern Cascadia subduction zone［J］．GeophysResLett，27：3117-3120．

McCaffrey R，Qamar A I，King R W，Wells R，Khazaradze G，Williams C A，Stevens C W，Vollick J J，Zwick P C．2007．Fault locking，block rotation and crustal deformation in the Pacific Northwest［J］．GeophysJInt，169：1315-1340，doi：10．1111／j．1365-246X．2007．03371．x．

Nishimura N，Yagi Y．2008．Rupture process for May 12，2008Sichuan earthquake［EB／OL］．［2009-06-05］．http：∥www．geol．tsukuba．a(chǎn)c．jp／～nisimura／20080512／．

Savage J C，Gan W J，Svarc J L．2001．Strain accumulation and rotation in the eastern California shear zone［J］．JGeophysRes，106（B10）：21995-22007．

Shen Z K，Sun J B，Zhang P Z，Wan Y G，Wang M，Bürgmann R，Zeng Y H，Gan W J，Liao H，Wang Q L．2009．Slip maxima at fault junctions and rupturing of barriers during the 2008Wenchuan earthquake［J］．NatureGeosci，2：718-724，doi：10．1038／NGEO636．

Wallace L M，Beavan J，McCaffrey R，Berryman K，Denys P．2006．Balancing the plate motion budget in the South Island，New Zealand using GPS，geological and seismological data［J］．GeophysJInt，168（1）：332-352，doi：10．1111／j．1365-246X．2006．03183．x．

Wallace L M，Beavan J，McCaffrey R，Darby D．2004a．Subduction zone coupling and tectonic block rotations in the North Island，New Zealand［J］．JGeophysRes，109（B12）：B12406，doi：10．1029／2004JB003241．

Wallace L M，Stevens C，Silver E，McCaffrey R，Loratung W，Hasiata S，Stanaway R，Curley R，Rosa R，Taugaloidi J．2004b．GPS and seismological constraints on active tectonics and arc-continent collision in Papua New Guinea：Implications for mechanics of microplate rotations in a plate boundary zone［J］．JGeophysRes，109（B05）：B05404，doi：10．1029／2003JB002481．

Wang Q，Qiao X J，Lan Q G，Jeffrey F，Yang S M，Xu C J，Yang Y L，You X Z，Tan K，Chen G．2011．Rupture of deep faults in the 2008Wenchuan earthquake and uplift of the Longmen Shan［J］．NatureGeosci，4：634-640，doi：10．1038／NGEO1210．

Xu X W，Wen X Z，Yu G H，Chen G H，Klinger Y，Hubbard J，Shaw J．2009．Coseismic reverse and oblique-slip surface faulting generated by the 2008MW7．9Wenchuan earthquake，China［J］．Geology，37（6）：515-518，doi：10．1130／G25462A．1．